北科大matlab数学实验报告第五次精华版.doc

《北科大matlab数学实验报告第五次精华版.doc》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、《数学实验》报告实验名称 数学实验第五次报告学院专业班级 姓名    学  号     2013年12月一、【实验目的】学会用MATLAB拟合曲线,学习二维拟合与peaks函数拟合.二、【实验任务】完成课后第6章第8,10,11题。三、【实验程序】8，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('one'),f1=poly2str(p1,'x')disp('three'),f3=poly2str(p3

2、,'x')disp('five'),f5=poly2str(p5,'x')x1=0:0.01:0.95;y1=polyval(p1,x1);y3=polyval(p3,x1);y5=polyval(p5,x1);plot(x,y,'rp',x1,y1,'--',x1,y3,'k-',x1,y5,'g-')legend('nihe','one','three','five')10,x=10:5:30;y=[25.229.831.231.729.4];x1=10:1:30;y1=interp1(x,y,x1,'nearest');y2=interp1(x,y,x1,'linear');y

3、3=interp1(x,y,x1,'spline');y4=interp1(x,y,x1,'cubic');subplot(2,2,1),plot(x,y,'rp',x1,y1),title('nearest')subplot(2,2,2),plot(x,y,'rp',x1,y1),title('linear')subplot(2,2,3),plot(x,y,'rp',x1,y1),title('spline')subplot(2,2,4),plot(x,y,'rp',x1,y1),title('cubic')disp('浓度X=18的抗压强度值')a=interp1(x,y,18,

4、'spline')disp('浓度X=26的抗压强度值')b=interp1(x,y,26,'cubic')11，[x,y]=meshgrid(-3:.5:3);z=peaks(x,y);[x1,y1]=meshgrid(-2:.1:2);z1=x1.*exp(-x1.^2-y1.^2);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,'*neare

5、st');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(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('立方插值')一、【实验结果】

6、8，10，11，一、【实验总结】一个小问题就可以引发很大的问题，只有细心做，才能做好。加油好好学习，才能有所收获。