《数值分析(第五版)计算实习题第三章(共7页).docx》由会员分享，可在线阅读，更多相关《数值分析(第五版)计算实习题第三章(共7页).docx（7页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精选优质文档-倾情为你奉上数值分析计算实习题第三章第二次作业：题一： x=-1:0.2:1;y=1./(1+25.*x.2);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-)hold onplot(x2,y2,k)legend(数据点,3次拟合曲线,3次多项式插值)xlabel(X),ylabel(Y)输出：f1 = 0.0000 -0.5752 0.0000 0.4841 f = (02675*x3)/ - (3305*x2)/5
2、746 + (76431*x)/ + 00613/40992 y1 = -0.0911 0.1160 0.2771 0.3921 0.4611 0.4841 0.4611 0.3921 0.2771 0.1160 -0.0911x2 = -1.0000 -0.7778 -0.5556 -0.3333 -0.1111 0.1111 0.3333 0.5556 0.7778 1.0000y2 =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
3、=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轴
4、)输出：p1 = -6.6221 12.8147 -4.6591 0.9266p2 = 2.8853 -12.3348 16.2747 -5.2987 0.9427Y1 = 0.9266 0.5822 0.4544 0.5034 0.9730 2.0103 2.4602Y2 = 0.9427 0.5635 0.4399 0.5082 1.0005 1.9860 2.4692p3 = 3.1316 -1.2400 0.7356Y3 = 0.7356 0.6429 0.6128 0.6454 0.8984 1.7477 2.6271 f1 = - (25075*x3)/42624 + (45435
5、*x2)/5328 - (809*x)/08 + 99011/40992 f2 = (98125*x4)/1312 - (19173*x3)/0656 + (70659*x2)/0656 - (22115*x)/0656 + 81293/42624 f3 = (18733*x2)/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数据点 ,
6、Rm为均方误差n=length(X)-1;max1=fix(n-1)/2);if mmax1 m=max1;endA=zeros(1,m+1);B=zeros(1,m+1);Ym=(Y(1)+Y(n+1)/2;Y(1)=Ym;Y(n+1)=Ym;A(1)=2*sum(Y)/n;for i=1:m B(i+1)=sin(i*X)*Y; A(i+1)=cos(i*X)*Y;endA=2*A/n;B=2*B/n;A(1)=A(1)/2;Y1=A(1);for k=1:m Y1=Y1+A(k+1)*cos(k*X1)+B(k+1)*sin(k*X1); Tm=A(1)+A(k+1).*cos(k*X)
7、+B(k+1).*sin(k*X);k=k+1;endY,Tm,Rm=(sum(Y-Tm).2)/n输出： X=-pi:2*pi/33:pi; Y=X.2.*cos(X);A,B,Y1,Rm=sanjiaobijin(X,Y,X1,16)输出：A = 1 至 12 列 -0.1397 4.4002 -2.8326 1.2355 -0.9128 0.7914 -0.7319 0.6982 -0.6773 0.6635 -0.6541 0.6474 13 至 17 列 -0.6426 0.6393 -0.6370 0.6355 -0.6348B = 1.0e-15 * 1 至 12 列 0 -0.0194 -0.0150 -0.0044 -0.0300 0.0105 0.0627 -0.0821 -0.0599 -0.0133 -0.0211 0.0297 13 至 17 列 0.0178 0.0962 -0.1049 0.0328 -0.0122即可得16插值多项式的值X1=-pi:0.001:pi;A,B,Y1,Rm=sanjiaobijin(X,Y,X1,16)plot(X,Y,r*,X1,Y1,b-.)legend(数据点,16次三角插值多项式)xlabel(X轴),ylabel(Y轴)专心-专注-专业