资源描述：
《2016矩阵与数值分析报告上机作业满分.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、第一题考虑计算给定向量的数：输入向量,输出。请编制一个通用程序，并用你编制的程序计算如下向量的数：对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(i);y1=y1+abs(Y(i)
2、);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.93
3、x的2数是1.24x的无穷数是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数是500500.00y的2数是18271.11y的无穷数是1000.00第二题考虑,其中定义f(0)=1，此时f(x)是连续函数。用此公式计算当时的函数值，画出图像。另一方面，考虑下面算法
4、：d=1+xifd=1theny=1elsey=lnd/(d-1)endif用此算法计算时的函数值，画出图像。比较一下发生了什么？代码：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+
5、i;if(d~=1)a(t)=log(d)/(d-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结果：第三题首先编制一个利用秦九韶算法计算一个多项式在给定点的函数值的通用程序，你的程序包括输入多项式的系数以及给定点，输出函数值。利用你编制的程序计算在x=2邻域附近的值。画出上
6、的图像。代码：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);结果
7、：第四题编制计算给定矩阵A的LU分解和PLU分解的通用程序，然后用你编制的程序完成下面两个计算任务：(1)考虑自己取定，并计算b=Ax。然后用你编制的不选主元和列主元的Gauss消去法求解该方程组，记你计算出的解为。对n从5到30估计计算解的精度。(2)对n从5到30计算其逆矩阵。代码：%LU分解functionproblem4()clc%生成矩阵Aforn=5:30A=-ones(n);A(:,n)=ones(n,1);foria=1:nA(ia,ia)=1;forja=(ia+1):(n-1)A(ia,ja)=
8、0;endendx=rand(n,1);b=A*x;%LU分解[n,n]=size(A);L=eye(n);M=eye(n);fori=1:n-1ifA(i,i)==0fprintf('Error:A(%d,%d)=0!',i,i);return;endforj=i+1:nL(j,i)=A(j,i)/A(i,i);endfork=i+1:nforl=i:
