北科大数学实验第五次报告

《北科大数学实验第五次报告》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、《数学实验》报告实验名称  拟合与插值学院专业班级    姓  名       学  号     2013年5月5一、【实验目的】学习和掌握MATLAB中实现多项式拟合和多项式插值。二、【实验任务】P130第9题P130第10题P130第12题三、【实验程序】P130第9题t=0:.5:3.0;v=[0,0.4794,0.8415,0.9975,0.9093,0.5985,0.1411];s=[1,1.5,2,2.5,3,3.5,4];p1=polyfit(t,v,2);p2=polyfit(t,s,1);p3=polyfit(v,

2、s,1);y1=polyval(p1,t);y2=polyval(p2,t);y3=polyval(p3,t);plot(t,v,'rp',t,y1,'--',t,y2,'k-.',t,y3);legend('拟合点','v、t关系','s、t关系','s、v关系')P130第10题X=10:5:30;Y=[25.2,29.8,31.2,31.7,29.4];Xi=10:1:30;Yi1=interp1(X,Y,Xi,'*nearest');Yi2=interp1(X,Y,Xi,'*linear');Yi3=interp1(X,Y,

3、Xi,'*spline');Yi4=interp1(X,Y,Xi,'*cubic');plot(X,Y,'ro',Xi,Yi1,'--',Xi,Yi2,'-',Xi,Yi3,'k.-',Xi,Yi4,'m:')legend('原始数据','最近点插值','线性插值','样条插值','立方插值')P130第12题[x,y]=meshgrid(-3:.5:3);z=x.^2/16-y.^2/9;[x1,y1]=meshgrid(-3:0.1:3);z1=x1.^2/16-y1.^2/9;5figure(1)subplot(1,2,1),

4、mesh(x,y,z),title('Êý¾Ýµã')subplot(1,2,2),mesh(x1,y1,z1),title('º¯ÊýͼÏó')[xi,yi]=meshgrid(-3:0.125:3);zi1=interp2(x,y,z,xi,yi,'*nexarest');zi2=interp2(x,y,z,xi,yi,'*linear');zi3=interp2(x,y,z,xi,yi,'*spline');zi4=interp2(x,y,z,xi,yi,'*cubic');figure(2)subplot(2,2,1),m

5、esh(xi,yi,zi1),title('×î½ü²åÖµµã')subplot(2,2,2),mesh(xi,yi,zi2),title('ÏßÐÔ²åÖµ')subplot(2,2,3),mesh(xi,yi,zi3),title('ÑùÌõ²åÖµ')subplot(2,2,4),mesh(xi,yi,zi4),title('Á¢·½²åÖµ')三、【实验结果】P130第9题P130第10题5由图像可以估算出浓度X=18时抗压强度Y=30.8；浓度X=26时抗压强度Y=31.6。P130第12题5三、【实验总结】学习到并掌握

6、了MATLAB中求多项式拟合和多项式插值的方法。5