数值分析计算实习题答案.docxVIP

数值分析计算实习题答案 【篇一：数值分析 (第五版 )计算实习题第三章】 t＞第二次作业： 题一： x=-1:0.2:1;y=1./(1+25.*x42); f1=polyfit(x,y,3) f=poly2sym(f1) y1=polyval(f1,x) x2=linspace(-1,1,10) y2=interp1(x,y,x2) plot(x,y,r*-,x,y1,b-) holdon plot(x2,y2,k) legend( 数据点 ,3 次拟合曲线 ,3 次多项式插值 ) xlabel(x),ylabel(y) 输出： f1 = 0.0000-0.5752 0.0000 0.4841 (4591875547102675*xT)/81129638414606681695789005144064 (1469057404776431*x)/20282409603651670423947251286016 4360609662300613/9007199254740992 y1 = -0.0911 0.1160 0.2771 0.3921 0.4611 0.4841 0.2771 0.1160-0.0911 (3305*xA2)/5746 + +0.4611 0.3921 x2 = -1.0000-0.7778-0.5556-0.3333-0.1111 0.1111 0.3333 0.5556 0.7778 1.0000 y2 = 0.0385 0.0634 0.1222 0.3000 0.7222 0.7222 0.3000 0.1222 0.0634 0.0385 题二： x=[0.0 0.1 0.2 0.3 0.5 0.8 1.0]; y=[1.0 0.41 0.50 0.61 0.91 2.02 2.46]; p1=polyfit(x,y,3) p2=polyfit(x,y,4) y1=polyval(p1,x) y2=polyval(p2,x) plot(x,y,r*,x,y1,b-.,x,y2,g--) p3=polyfit(x,y,2) y3=polyval(p3,x) f1=poly2sym(p1) f2=poly2sym(p2) f3=poly2sym(p3) plot(x,y,r*,x,y1,b-.,x,y2,g--,x,y3,m--) legend( 数据点 ,3 次多项式拟合 ,4 次多项式拟合 ,2 次多项式拟合 ) xlabel(x 轴),ylabel(y 轴) 输出： p1 = -6.622112.8147-4.6591 0.9266 p2 = 2.8853 -12.334816.2747-5.2987 0.9427 y1 = 0.9266 0.5822 0.4544 0.5034 0.9730 2.0103 2.4602 y2 = 0.9427 0.5635 0.4399 0.5082 1.0005 1.9860 2.4692 p3 = 3.1316-1.2400 0.7356 y3 = 0.7356 0.6429 0.6128 0.6454 0.8984 1.7477 2.6271 f1 = -(7455778416425075*xT)/1125899906842624 + (1803512222945435*x9)/140737488355328 - (40981580032809*x)/8796093022208 + 8345953784399011/9007199254740992 f2 = (1624271450198125*xA4)/562949953421312 - (3471944732519173*xA3)/281474976710656 + (4580931990070659*xA2)/281474976710656 - (1491459232922115*x)/281474976710656 + 1061409433081293/1125899906842624 f3 = (18733*xA2)/5982 - (74179*x)/59820 + 73337/99700 题三： 建立三角插值函数的 m 文件 function [a,b,y1,rm]=sanjiaobijin(x,y,x1,m)%a b 分别是 m 阶三角 多项式 tm (x )的系数 aj，bj (j=1,2,...,m )的系数矩阵， y1 是 tm (x )在x1处的值，x y数据点,rm为均方误差 n=length(x)-1;max1=fix((n-1)/2); if mmax1 m=max1; end a=zeros(1,m+1);b=zeros(1,m+1); ym=(y(1)+y(n+1))/