资源描述：
《大连理工大学矩阵与数值分析上机作业》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、矩阵与数值分析上机作业学校：大连理工大学学院：班级：姓名：学号：授课老师：50注：编程语言Matlab程序：Norm.m函数functions=Norm(x,m)%求向量x的范数%m取1,2,inf分别表示1,2，无穷范数n=length(x);s=0;switchmcase1%1-范数fori=1:ns=s+abs(x(i));end50case2%2-范数fori=1:ns=s+x(i)^2;ends=sqrt(s);caseinf%无穷-范数s=max(abs(x));end计算向量x，y的范数Test1.mclear
2、all;clc;n1=10;n2=100;n3=1000;x1=1./[1:n1]';x2=1./[1:n2]';x3=1./[1:n3]';y1=[1:n1]';y2=[1:n2]';y3=[1:n3]';disp('n=10时');disp('x的1-范数:');disp(Norm(x1,1));disp('x的2-范数:');disp(Norm(x1,2));disp('x的无穷-范数:');disp(Norm(x1,inf));disp('y的1-范数:');disp(Norm(y1,1));disp('y的2-范数
3、:');disp(Norm(y1,2));disp('y的无穷-范数:');disp(Norm(y1,inf));disp('n=100时');disp('x的1-范数:');disp(Norm(x2,1));disp('x的2-范数:');disp(Norm(x2,2));disp('x的无穷-范数:');disp(Norm(x2,inf));disp('y的1-范数:');disp(Norm(y2,1));disp('y的2-范数:');disp(Norm(y2,2));50disp('y的无穷-范数:');disp(N
4、orm(y2,inf));disp('n=1000时');disp('x的1-范数:');disp(Norm(x3,1));disp('x的2-范数:');disp(Norm(x3,2));disp('x的无穷-范数:');disp(Norm(x3,inf));disp('y的1-范数:');disp(Norm(y3,1));disp('y的2-范数:');disp(Norm(y3,2));disp('y的无穷-范数:');disp(Norm(y3,inf));运行结果：n=10时x的1-范数:2.9290；x的2-范数:1
5、.2449；x的无穷-范数:1y的1-范数:55；y的2-范数:19.6214；y的无穷-范数:10n=100时x的1-范数:5.1874；x的2-范数:1.2787；x的无穷-范数:1y的1-范数:5050；y的2-范数:581.6786；y的无穷-范数:100n=1000时x的1-范数:7.4855；x的2-范数:1.2822；x的无穷-范数:1y的1-范数:500500；y的2-范数:1.8271e+004；y的无穷-范数:100050程序Test2.mclearall;clc;n=100;%区间h=2*10^(-15
6、)/n;%步长x=-10^(-15):h:10^(-15);%第一种原函数f1=zeros(1,n+1);fork=1:n+1ifx(k)~=0f1(k)=log(1+x(k))/x(k);elsef1(k)=1;endendsubplot(2,1,1);plot(x,f1,'-r');axis([-10^(-15),10^(-15),-1,2]);legend('原图');%第二种算法f2=zeros(1,n+1);fork=1:n+1d=1+x(k);if(d~=1)f2(k)=log(d)/(d-1);elsef2(k
7、)=1;50endendsubplot(2,1,2);plot(x,f2,'-r');axis([-10^(-15),10^(-15),-1,2]);legend('第二种算法');运行结果：显然第二种算法结果不准确，是因为计算机中的舍入误差造成的，当时，，计算机进行舍入造成恒等于1，结果函数值恒为1。程序：秦九韶算法:QinJS.mfunctiony=QinJS(a,x)50%y输出函数值%a多项式系数，由高次到零次%x给定点n=length(a);s=a(1);fori=2:ns=s*x+a(i);endy=s;计算p(
8、x):test3.mclearall;clc;x=1.6:0.2:2.4;%x=2的邻域disp('x=2的邻域：');xa=[1-18144-6722016-40325376-46082304-512];p=zeros(1,5);fori=1:5p(i)=QinJS(a,x(i));en
