大连理工大学矩阵与数值分析上机作业

矩阵与数值分析上机作业 学校： 大连理工大学 学院： 班级： 姓名： 学号： 授课老师： 注：编程语言 Matlab 程序：Norm.m 函数function s=Norm(x,m)%求向量 x 的范数%m 取 1,2,inf 分别 表示 1,2，无穷范数n=length(x);s=0;switch m case 1 %1-范数 for i=1:n s=s+abs(x(i)); end case 2 %2-范数 for i=1:n s=s+x(i)^2; end s=sqrt(s); case inf %无穷-范数 s=max(abs(x));end 计算向量 x，y 的范数Test1.mclear 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-范数:');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.mclear all;clc;n=100;%区间h=2*10^(-15)/n;%步长x=-10^(-15):h:10^(-15);%第一种原函数f1=zeros(1,n+1);for k=1:n+1 if x(k)~=0 f1(k)=log(1+x(k))/x(k); else f1(k)=1; endendsubplot(2,1,1);plot(x,f1,'-r');axis([-10^(-15),10^(-15),-1,2]);l...