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

矩阵与数值分析上机作业 学校： 大连理工大学 学院： 班级： 姓名： 学号： 授课老师： 注：编程语言 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.m clear 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 的无穷-范数:1 y 的 1-范数:55； y 的 2-范数:19.6214； y 的无穷-范数:10 n=100 时 x 的 1-范数:5.1874；x 的 2-范数: 1.2787； x 的无穷-范数:1 y 的 1-范数:5050； y 的 2-范数:581.6786； y 的无穷-范数:100 n=1000 时 x 的 1-范数:7.4855； x 的 2-范数:1.2822； x 的无穷-范数:1 y 的 1-范数: 500500； y 的 2-范数:1.8271e+004；y 的无穷-范数:1000 程序 Test2.m clear 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; end end subplot(2,1,1); plot(x,f1,-r); axis([-10^(-15),10^(-15),-1,2]); legend(