北科大MATLAB第五次作业.doc

《数学实验》报告 实验名称?? MATLAB绘图 学 院 东凌经济管理学院 专业班级?????????工管111 姓??名???????? 李伊童 学??号??????? ?? 2013年 5月 【实验目的】 插值与拟合 【实验任务】 第六章8、10、12 【实验程序】 8、x=[0.10,0.20,0.30,0.55,0.70,0.80,0.95]; y=[15,18,19,21,22.6,23.8,26]; p1=polyfit(x,y,1); p3=polyfit(x,y,3); p5=polyfit(x,y,5); disp(‘一阶拟合函数),f1=poly2str(p1,x); disp(‘三阶拟合函数),f3=poly2str(p3,x); disp(‘五阶拟合函数),f5=poly2str(p5,x); x1=0.10: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); legend(‘拟合点’,一次拟合,三次拟合,五次拟合); 10、x=[10,15,20,25,30]; y=[25.2,29.8,31.2,31.7,29.4]; xi=10:0.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:), grid on; legend(原始数据,最近值插点,线性插点,样条插点,立方插点); 12、[x,y]=meshgrid(-3:0.3:3); z=x.^2/16-y.^2/9; [x1,y1]=meshgrid(-3:0.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:0.05:3); zi1=interp2(x,y,z,xi,yi,*nearest); 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(立方插值) 【实验结果】 8、 10、 12、 【实验总结】 此次作业比较综合，与前面二维绘图和三维绘图相关。 5