2016矩阵与数值分析上机作业满分.docx

《2016矩阵与数值分析上机作业满分.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第一题考虑计算给定向量的范数：输入向量x=x1,x2,…,xnT,输出x1,x2,x∞。请编制一个通用程序，并用你编制的程序计算如下向量的范数：x=1,12,13,…,1nT,y=1,2,…,nT.对n=10，100，1000甚至更大的n计算其范数，你会发现什么结果?你能否修改你的程序使得计算结果相对精确呢？代码：functionProblem1()clcN=input('inputN=');X=zeros(1,N);Y=1:N;fori=1:NX(i)=1/i;endx1=0;x2=0;y1=0;y2=0;fori=1:Nx1=x1+abs(X(i));x2=x2+X(i)*X(

2、i);y1=y1+abs(Y(i));y2=y2+Y(i)*Y(i);endx1;x2=sqrt(x2);x3=max(abs(X));y1;y2=sqrt(y2);y3=max(abs(Y));fprintf('x的1范数是%.2f',x1);fprintf('x的2范数是%.2f',x2);fprintf('x的无穷范数是%.2f',x3);fprintf('y的1范数是%.2f',y1);fprintf('y的2范数是%.2f',y2);fprintf('y的无穷范数是%.2f',y3);结果：inputN=10x的1范数是2.93x的2范数是1.24

3、x的无穷范数是1.00y的1范数是55.00y的2范数是19.62y的无穷范数是10.00inputN=100x的1范数是5.19x的2范数是1.28x的无穷范数是1.00y的1范数是5050.00y的2范数是581.68y的无穷范数是100.00inputN=1000x的1范数是7.49x的2范数是1.28x的无穷范数是1.00y的1范数是.00y的2范数是18271.11y的无穷范数是1000.00第二题考虑y=fx=ln⁡(1+x)x,其中定义f(0)=1，此时f(x)是连续函数。用此公式计算当x∈-10-15,10-15时的函数值，画出图像。另一方面，考虑下面算法：d=1+

4、xifd=1theny=1elsey=lnd/(d-1)endif用此算法计算x∈-10-15,10-15时的函数值，画出图像。比较一下发生了什么？代码：functionproblem2()clcN=1000;n=2*10^(-15)/N;%公式计算f=zeros(1,N+1);t=1;fori=-10^(-15):n:10^(-15)if(i~=0)f(t)=log(1+i)/i;elsef(t)=1;endt=t+1;end%算法计算a=zeros(1,N+1);t=1;fori=-10^(-15):n:10^(-15)d=1+i;if(d~=1)a(t)=log(d)/(d-

5、1);elsea(t)=1;endt=t+1;end%画图，左边是公式算的，右边是算法算的subplot(1,2,1);plot(-10^(-15):n:10^(-15),f);holdonsubplot(1,2,2);plot(-10^(-15):n:10^(-15),a);holdon结果：第三题首先编制一个利用秦九韶算法计算一个多项式在给定点的函数值的通用程序，你的程序包括输入多项式的系数以及给定点，输出函数值。利用你编制的程序计算px=x-29=x9-18x8+144x7-672x6+2016x5-4032x4+5376x3-4608x2-512在x=2邻域附近的值。画出p

6、x在x∈1.95,20.5上的图像。代码：functionproblem3()%秦九韶算法clcn=9;c=[1,-18,144,-672,2016,-4032,5376,-4608,2304,-512];N=100;m=(2.05-1.95)/N;y=zeros(1,N+1);ydex=0;forx=1.95:m:2.05ydex=ydex+1;dex=1;temp=c(dex);whiledex<=ndex=dex+1;temp=temp*x+c(dex);endy(ydex)=temp;endplot(1.95:m:2.05,y);结果：第四题编制计算给定矩阵A的LU分解和P

7、LU分解的通用程序，然后用你编制的程序完成下面两个计算任务：(1)考虑A=10⋯-1⋱⋱⋮⋱⋱0⋮01⋮1-1⋯-111-1⋯-1-11∈Rn×n,自己取定x∈Rn，并计算b=Ax。然后用你编制的不选主元和列主元的Gauss消去法求解该方程组，记你计算出的解为x。对n从5到30估计计算解的精度。(2)对n从5到30计算其逆矩阵。代码：%LU分解functionproblem4()clc%生成矩阵Aforn=5:30A=-ones(n);A(:,n)=ones(n,1)