《矩阵与数值分析》由会员分享，可在线阅读，更多相关《矩阵与数值分析（4页珍藏版）》请在金锄头文库上搜索。
1、本文格式为Word版，下载可任意编辑矩阵与数值分析 矩阵与数值分析 学 院 专 业 班 级 学 号 姓 名 电子信息与电气工程学部 生物医学工程 刘江涛 1：考虑计算给定向量的范数；输入向量x?(x1,x2,?,xn)T，输出x1,x2,x?，请编制一个通用程序，并用你编制的程序计算如下向量的范数： 1?11Tx?1,?,?,y?1,2,?,n? n?23对n?10,100,1000甚至更大的n计算其范数，你会察觉什么结果？你能否修改你的程序使得计算结果相对精确呢？ 通用求范数程序： function NORM(x) y1=sum(abs(x); y2=(sum(x.2)(1/2); y3=m
2、ax(abs(x); fprintf(1-范数=%g； 2-范数= %g； inf-范数=%gn,y1,y2,y3); 例题的运行程序： function xianglaing(n) x=; y=; for i=1:n x(i)=1/i; y(i)=i; end disp(x 的范数：); NORM(x); disp( ) disp(y 的范数：); NORM(y); 运行结果如下表： T n 范数 10 100 1000 x1 2.92897 5.18738 7.48547 x2 x? 1 1 1 1.2449 1.27866 1.28216 n 范数 10 100 1000 y1 55 5
3、050 500500 y2 y? 10 100 1000 19.6214 581.679 18271.1 根据上述的两个表的运行结果，我们可以得知无论n的值如何变化，对于x?1恒成立；y?n恒成立，其1-范数与2-范数随着n的增大而增大，但是其变化越来越小，这是由于计算在举行数值计算时有误差存在，对于表达式(1)当n很大时 1却很n小，会展现“大数吃小数的现象”；修改方案：当n很大时我们制止用n做除数，由于当n分外大时 1?0成立；所以在求解其范数时我们从小数开头相加，无穷个分外n小的数值相加也可能是个很大的数，从而可以制止两个数相加时展现“大数吃小数”的现象； 2：考虑y?f(x)?ln(1
4、?x)，其中定义f(0)?1，此时f(x)是连续函数，用此公x式计算当x?10?15,10?15时的函数值，画出图像。另一方面，考虑下面算法： d?1?x;ifd?1theny?1elsey?lnd/(d?1)endif 用此算法计算x?10?15,10?15时的函数值，画出图像，对比一下发生了什么？ 程序： x=-10(-15):10(-20):10(-15); if (x=0) f=1; else f=log(1+x)/x; end figure(1) plot(x,f); d=1+x; if d=1 y=1; else y=log(d)/(d-1); end figure(2) Plot
5、(x,y); 有图可知，直接用公式f(x)?ln(1?x)计算x?10?15,10?15的函数值时，除了在xx?0出的值为1，其他的值都是无限趋近于1；而利用算法二算出的结果全为1；出 现这这处境的理由是x的取值分外接近于0，在用公式d?1?x求d得过程中展现了大数吃小数的处境，所以在用计算机计算时d?1恒成立，从而使y?1恒成立； 3： 首先编写一个利用秦九韶算法计算一个多项式在定点的函数值的通用程序，你的程序包括输入多项式的系数以及定点，输出函数值，利用你编写的程序计算 f(x)?(x?2)9?x9?18x8?144x7?672x6?2022x5?4032x4?5376x3?4608x2?2304x?512在x?2邻域邻近的值，画出p(x)在x?1.95,20.5上地图像。 秦九韶算法的通用程序： %A为多项式的以升幂排列的系数，x为初始值 function p=qinjiushao(A,x) a=A; ,n=size(a); n=n-1; S=; S(n+1)=a(n+1); for k=n:-1:1 S(k)=x.*S(k+1)+a(k); 4
