11级数值方法计算实习题.docVIP

信计、数应专业11级数值方法计算实习题 要求：1、用Matlab语言或你熟悉的其他语言编写程序，使之尽可能具有通用性；2、根据上机计算实践，对所使用的数值方法的特点、性质、有效性、误差和收敛性等方面进行必要的讨论和分析；3、完成计算后写出实验报告，内容包括：课题名称、解决的问题、采用的数值方法、算法程序、数值结果、对实验结果的讨论和分析等；4、特别说明：严禁抄袭，否则一经发现，所有雷同实验报告最多评为及格。 一、在区间[-1,1]上分别取用两组等距节点对龙格函数作多项式插值及三次样条插值，对每个值，分别画出插值函数即的图形。 解：作m文件，文件名：lagrange.m function lagrange(n) x0=linspace(-1,1,n+1); y0=1./(1+25*x0.^2); x=linspace(-1,1,100); y=1./(1+25*x.^2); y3=interp1(x0,y0,x,spline); for k=1:length(x) im=ones(1,n+1); for i=1:n+1 for j=1:n+1 if x0(i)~=x0(j) im(i)=im(i)*(x(k)-x0(j))/(x0(i)-x0(j)); end end end yp(k)=im*y0; end plot(x0,y0,*,x,y,b,x,y3,-,x,yp,--); legend(原插值点,原函数,三次样条插值,多项式插值) 输入lagrange(10)调用函数，得到n=10时的插值函数与f（x）的函数图形 输入lagrange(20)调用函数，得到n=20时的插值函数与f（x）的函数图形 n=10时的插值函数与f（x）的函数图形 n=20时的插值函数与f（x）的函数图形 二、已知Wilson矩阵，且向量，则方程组有准确解。 ⑴用Matlab 内部函数求，的所有特征值和； 解： M=det(A) D=eig(A) N=cond(A) M = D = N = 1 0.0102 2.9841e+003 0.8431 3.8581 30.2887 ⑵令，解方程组，并求出向量和，从理论结果和实际计算结果两方面分析方程组解的相对误差与的相对误差的关系； 解：为了方便MATLAB的编辑，令B=A+δA;M=δA;X=x+δx;t=δx; B=[10 7 8.1 7.2;7.08 5.04 6 5;8 5.98 9.89 9;6.99 4.99 9 9.98] ; b=[32;23;33;31]; X=B\b norm(t) X = ans = -81.0000 163.9695 % ‖δx‖的结果 137.0000 norm(t)/norm(x) -34.0000 ans = 22.0000 81.9848 %的结果 (3.7347e+003)/(-11.7943) m=B-A;norm(m) ans = ans = -316.6530 0.2308 %||δA ||的结果 X=B\b norm(m)/norm(A) X = ans = -81.0000 0.00