MATLAB软件中软件拟合与插值运算的方法.doc

注：报告内容根据具体实验课程或实验项目的要求确定，一般包括实验目的、实验仪器、原理摘要、数据记录及结果分析等。如纸张不够请自行加纸。 实验目的 学会MATLAB软件中软件拟合与插值运算的方法。 实验内容 x=[0.5,1.0,1.5,2.0,2.5,3.0]; y=[1.75,2.45,3.81,4.80,7.00,8.60]; p=polyfit (x,y,2); x1=0.5:0.05:3.0; y1=polyval(p,x1); plot(x,y,*r’, ’-b’) plot(x,y,r,x1,y1,b) 问题1.32：将上题改为6阶拟合，拟合效果又如何？ 例： product=[75.995,91.972,105.711,123.203,131.669,249.633,256.344,267.893,300,310,320,330]; year=1900:10:2010; p2005=interp1(year,product,2005); x=1900:10:2010; y=interp1(year,product,x,cubic); plot(year,product,o,x,y) 例： years=1950:10:1990; service=10:10:30; wage=[150.697,199.592,187.625;179.323,195.072,250.287;203.212,179.092,322.767;226.505,153.706,426.730;249.636,120.281,598.243]; w=interp2(service,years,wage,15,1975) w = 6672/35 x=1:6;y=1:4; t=[12,10,11,11,13,15;16,22,28,35,27,20;18,21,26,32,28,25;20,25,30,33,32,30]; subplot(121) x1=1:0.1:6; y1=1:0.1:4; [x2,y2]=meshgrid(x1,y1); t1=interp2(x,y,t,x2,y2,cubic); subplot(1,2,2) mesh(x1,y1,t1) subplot(121) mesh(x,y,t) 解答：1） x=[1.2,1.8,2.1,2.4,2.6,3.0,3.3];y=[4.85,5.2,5.6,6.2,6.5,7.0,7.5]; p=polyfit (x,y,6); x1=1.2:0.01:3.3; y1=polyval(p,x1); plot(x,y,*r,x1,y1,-b) 2） 二阶： x=[0.1:0.01:3.0]; y=[sin(x)]; p=polyfit(x,y,2) p = -0.4286 1.3454 -0.0706 x1=0.1:0.01:3.0; y1=polyval(p,x1)； plot(x,y,*r,x1,y1,b) 三阶 p=polyfit(x,y,3) p = -0.0031 -0.4143 1.3271 -0.0652 x1=0.1:0.01:3.0; y1=polyval(p,x1); plot(x,y,*r,x1,y1,b) 五阶 p=polyfit(x,y,5) p = 0.0002 0.0366 -0.2340 0.0579 0.9787 0.0024 x1=0.1:0.01:3.0; y1=polyval(p,x1); plot(x,y,*r,x1,y1,b) 可以看出，阶数越高拟合效果越好 3 x=[0.1,0.8,1.3,1.9,2.5,3.1];y=[1.2,1.6,2.7,2.0,1.3,0.5]; p2=interp1(x,y,2) p2 = 1.8833 x1=[0.1:0.1:3.1]; y1=interp1(x,y,x1,cubic) plot(x,y,x,x1,y1) 实验内容 成绩评语 评： 指导老师： 年 月 日 成绩