矩阵与数值分析上机作业学校：大连理工大学学院：班级：姓名：学号：授课老师：注：编程语言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));endcase2%2-范数fori=1:ns=s+x(i)^2;ends=sqrt(s);caseinf%无穷-范数s=max(abs(x));end计算向量x，y的范数Test1.mclearall;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-范数:');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));disp('y的无穷-范数:');disp(Norm(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.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的无穷-范数:1000程序Test2.mclearall;clc;n=100;%区间h=2*10^(-15)/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)=1;endendsubplot(2,1,2);plot(x,f2,'-r');axis([-10^(-15),10^(-15),-1,2]);legend('第二种算法');运行结果：显然第二种算法结果不准确，是因为计算机中的舍入误差造成的，当]10,10[1515x时，xd1，计算机进行舍入造成d恒等于1，结果函数值恒为1。程序：秦九韶算法:QinJS.mfunctiony=QinJS(a,x)%y输出函数值%a多项式系数，由高次到零次%x给定点n=length(a);s=a(1);fori=2:ns=s*x+a(i);endy=s;计算p(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));enddisp('相应多项式p值：');pxk=1.95:0.01:20.5;nk=length(xk);pk=zeros(1,nk);k=1;fork=1:nkpk(k)=QinJS(a,xk(k));endplot(xk,pk,'-r');xlabel('x');ylabel('p(x)');运行结果：x=2的邻域：x=1.60001.80002.00002.20002.4000相应多项式p值：p=1.0e-003*-0.2621-0.000500.00050.2621p(x)在x[1.95,20.5]上的图像程序：LU分解，LUDe..mfunction[L,U]=LUDe.(A)%不带列主元的LU分解N=size(A);n=N(1);L=eye(n);U=zeros(n);fori=1:nU(1,i)=A(1,i);L(i,1)=A(i,1)/U(1,1);endfori=2:nforj=i:nz=0;fork=1:i-1z=z+L(i,k)*U(k,j);endU(i,j)=A(i,j)-z;endforj=i+1:nz=0;fork=1:i-1z=z+L(j,k)*U(k,i);endL(j,i)=(A(j,i)-z)/U(i,i);endendPLU分解，PLUDe..mfunction[P,L,U]=PLUDe.(A)%带列主元的LU分解[m,m]=size(A);U=A;P=eye(m);L=eye(m);fori=1:mforj=i:mt(j)=U(j,i);fork=1:i-1t(j)=t(j)-U(j,k)*U(k,i);endenda=i;b=abs(t(i));forj=i+1:mifb
