北科大matlab实验报告第五次作业matlab拟合插值

《北科大matlab实验报告第五次作业matlab拟合插值》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、《数学实验》报告实验名称  第五堂课练习学院东凌经济管理学院专业班级    管信112     姓  名    杨忠璟    学  号    41171061     2013年5月6一、【实验目的】1.掌握多项式及其相关运算2.掌握多项式的基本运算3.掌握多项式的拟合4.掌握多项式的插值二、【实验任务】练习：课本6.7,6.8,6.10,6.12.三、【实验程序】1)6.7程序：分别用2,3,4,6阶多项式拟合函数y=cos(x)，并将拟合曲线与函数曲线y=cos(x)进行比较clf;x=0:pi/50:2*pi;y=cos(x)p2=polyfit(x,y,2);p3=polyfit

2、(x,y,3);p4=polyfit(x,y,4);p6=polyfit(x,y,6);disp('二阶拟合函数'),f2=poly2str(p2,'x')disp('三阶拟合函数'),f3=poly2str(p3,'x')disp('四阶拟合函数'),f4=poly2str(p4,'x')disp('六阶拟合函数'),f6=poly2str(p6,'x')x1=0:pi/50:2*pi;y2=polyval(p2,x1);y3=polyval(p3,x1);y4=polyval(p4,x1);y6=polyval(p6,x1);plot(x,y,'r-',x1,y2,'--',x1,y

3、3,'k.',x1,y4,'m*',x1,y6,'gs');legend('原函数','二阶拟合','三阶拟合','四阶拟合','六阶拟合')62)6.8程序：在钢线碳含量对于电阻的效应的研究中，得到以下数据，分别用一次、三次、五次多项式拟合曲线来拟合这组数据并画出图形clf;x=[0.100.300.400.550.700.800.95];y=[1518192122.623.826];p1=polyfit(x,y,1);p3=polyfit(x,y,3);p5=polyfit(x,y,5);disp('一阶拟合函数'),f1=poly2str(p1,'x')disp('三阶拟合函数')

4、,f3=poly2str(p3,'x')disp('五阶拟合函数'),f5=poly2str(p5,'x')y1=polyval(p1,x);y3=polyval(p3,x);y5=polyval(p5,x);6plot(x,y,'r-',x,y1,'--',x,y3,'k.',x,y5,'m*');legend('原函数','一阶拟合','三阶拟合','五阶拟合')1)6.10程序：在某种添加剂的不同浓度之下对铝合金进行抗拉强度实验，得到数据如下，现分别使用不同的插值方法，对其中间没有测度的浓度进行推测，并估算出浓度X=18及X=26时的抗压强度Y的值.x=[10,15,20,25,3

5、0];y=[25.2,29.8,31.2,31.7,29.4];xi=10:.5:30;yi1=interp1(x,y,xi,'*nearest');yi2=interp1(x,y,xi,'*linear');yi3=interp1(x,y,xi,'*spline');yi4=interp1(x,y,xi,'*cubic');plot(x,y,'ro',xi,yi1,'--',xi,yi2,'-',xi,yi3,'k.-',xi,yi4,'m:'),gridonlegend('原始数据','最近点插值','线性插值','样条插值','立方插值')2)6.12程序：求下列线性方程组通解：[

6、x,y]=meshgrid(-3:.5:3);z=x.^2/16-y.^2/9;[x1,y1]=meshgrid(-3:.1:3);z1=x1.^2/16-y1.^2/9;figure(1)subplot(1,2,1),mesh(x,y,z),title('数据点')subplot(1,2,2),mesh(x1,y1,z1),title('函数图象')[xi,yi]=meshgrid(-3:.125:3);zi1=interp2(x,y,z,xi,yi,'*nearest');zi2=interp2(x,y,z,xi,yi,'*linear');zi3=interp2(x,y,z,xi,

7、yi,'*spline');zi4=interp2(x,y,z,xi,yi,'*cubic');figure(2)subplot(221),mesh(xi,yi,zi1),title('最近点插值')subplot(222),mesh(xi,yi,zi2),title('线性插值')subplot(223),mesh(xi,yi,zi3),title('样条插值')subplot(224),mesh(xi,yi,zi4),title