1、数学实验报告实验名称 拟合和插值 学 院 机械工程学院 专业班级 姓 名 学 号 2010 年 10 月一、 【实验目的】 二、 【实验任务】8、10、12三、 【实验程序】8、x=0.1,0.3,0.4,0.55,0.7,0.8,0.95;y=2.1,3.2,2.1,2.5,3.2,3.5,3.4,4.1,4.7,5.0,4.8p1=polyfit(x,y,1) p3=polyfit(x,y,3);p5=polyfit(x,y,5);disp(一次拟合函数),f2=poly2str(p1,x)disp(三次拟合函数),f3=poly2str(p3,x);disp(五次拟合函数),f5=pol
2、y2str(p5,x);x1= 0.1,0.3,0.4,0.55,0.7,0.8,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;p7=polyfit(x,y,7);disp(七次拟合函数),f7=poly2str(p7,x);x1=10,15,20,25,30;y7=polyval(p7,x1);p
3、lot(x,y,rp,x1,y7);y18=polyval(p7,18)y26=polyval(p7,26)legend(拟合点,七次拟合点)y18 =30.1838y26=30.306212、x,y=meshgrid(-3:0.5:3);z=x.2/16-y.2/9;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,yi,*spline);zi4=interp2(x,y,z,xi,yi,*cubic);figure(2)subplot(2,2,1),mesh(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(立方插值)四、 【实验结果】五、 【实验总结】
