MATLAB軟件中軟件擬合與插值運算的方法.doc

《MATLAB軟件中軟件擬合與插值運算的方法.doc》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在應用文檔-天天文庫。

1、注：報告內(nèi)容根據(jù)具體實驗課程或?qū)嶒烅椖康囊蟠_定，一般包括實驗目的、實驗儀器、原理摘要、數(shù)據(jù)記錄及結(jié)果分析等。如紙張不夠請自行加紙。實驗目的學會MATLAB軟件中軟件擬合與插值運算的方法。實驗內(nèi)容>>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：

2、將上題改為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;>>se

3、rvice=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

4、,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:

5、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.42861.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.41431.3271-0.0652>>x1=0.1:0.01:3.0;>>y1=polyval(p,x1)

6、;>>plot(x,y,'*r',x1,y1,'b')>>五階>>p=polyfit(x,y,5)p=0.00020.0366-0.23400.05790.97870.0024>>x1=0.1:0.01:3.0;>>y1=polyval(p,x1);>>plot(x,y,'*r',x1,y1,'b')可以看出，階數(shù)越高擬合效果越好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];>>y

7、1=interp1(x,y,x1,'cubic')plot(x,y,'x',x1,y1)實驗內(nèi)容成績&評語評：指導老師：年月日成績