北科大MATLAB第五次作业

《数学实验》报告

实验名称 MATLAB绘图

学院东凌经济管理学院

专业班级工管111 姓名李伊童

学号 41171034

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、

五、【实验总结】

此次作业比较综合，与前面二维绘图和三维绘图相关。