資源描述：
《matlab軟件擬合與插值運算實驗報告》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫。
1、實驗6數(shù)據(jù)擬合&插值一.實驗?zāi)康膶W會MATLAB軟件中軟件擬合與插值運算的方法。二.實驗內(nèi)容與要求在生產(chǎn)和科學實驗中，自變量X與因變量y=f(x)的關(guān)系式有吋不能直接寫出表達式，而只能得到函數(shù)在若干個點的函數(shù)值或?qū)?shù)值。當要求知道觀測點之外的函數(shù)值時，需要估計函數(shù)值在該點的值。要根據(jù)觀測點的值，構(gòu)造一個比較簡單的函數(shù)yn(x),使函數(shù)在觀測點的值等于已知的數(shù)值或?qū)?shù)值，尋找這樣的函數(shù)l(x),辦法是很多的。根據(jù)測量數(shù)據(jù)的類型有如下兩種處理觀測數(shù)據(jù)的方法。(1)測量值是準確的，沒有誤差，一般用插值。
2、(2)測量值與真實值有誤差，一般用曲線擬合。MATLAB屮提供了眾多的數(shù)據(jù)處理命令，有插值命令，擬合命令。1?曲線擬合?x二[0.5,1.0,1.5,2.0,2.5,3.0];?y=[l.75,2.45,3.81,4.80,7.00,&60];?p=polyfit(x,y,2);?xl=0.5:0.05:3.0;?yl=polyval(p,xl);?plot(x,y；*r,,xl,yl；-b,)2.一維插值?year=[1900,1910,1920,1930,1940,1990,2000,2010
3、];?product=[75.995,91.972,105.711,123.203,131.669,249.633,256.344,267.893J;?p2005=interpl(year,product,2005)p2005=262.1185?y=interpl(year,products,'cubic');?plot(year,product,'o*,x,y)3?二維插值?years=1950:10:1990;?service=10:10:30;?wage=[150.697,199.592,1
4、87.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=190.6288[例1.98]x=l:6;y=l:4;t二[12,10,ll,ll,13,15;16,22,28,35,27,20;18,21,26,32,28,25;20,25,30,33,32,30];subp
5、lot(1,2,1)mesh(x,y,t)xl=l:0.1:6;yl=l:0.1:4;[x2,y2]=meshgrid(xl,yl);11=interp2(x,y,t,x2,y2,'cubic1);subplot(1,2,2)mesh(xhyl.tl)三，練習與思考1)已知x二[1.2,1.&2」,2.4,2.6,3.0,3.3],y=[4.85,5.2,5.6,6.2,6.5,7.0,7.5],求對x和y進行6階多項式擬合的系數(shù).x=[l.2,1.8,2.1,2.4,2.6,3.0,331；y=
6、[4?85,5?2,5?6,6?2,6?5,7?0,7?5];?p=polyfit(x,y,6)P=-2.010729.0005-170.6763523.2180-878.3092763.9307-263.4667xl=0.5:0.05:3.0;?yl=polyval(p,xl);?plot(x,y嚴門xl,yl1)分別用2,3,4,5階多項式來逼近［0,3］上的正弦函數(shù)sinx,并做出擬合曲線及sinx函數(shù)曲線圖，了解多項式的逼近程度和有效擬合區(qū)間隨多項式的階數(shù)有何變化.(2)2階：?x=0:0
7、.01:3;?y=sin(x);?p=polyfit(x,y,2);?xl=0:0.01:3;?yl=polyval(p,xl);?plot(x,y；*rxl,yl；-b,)3階：?p=polyfit(x,y,3)；?xl二0:0.01:3;?yl=polyval(p,xl);?plot(x,y；*rxl,yl；-b,)4階：?p=polyfit(x,y,4);?xl=0:0.01:3;?yl=polyval(p,xl);?plot(x,y；*rxhyl；-br)5階：?p=polyfit(
8、x,y,5);?xl二0:0.01:3;?yl=polyval(p,xl);?plot(x,y；*rxl,yl；-b,)1)已知x二[0丄0?&1?3丄925,3?1]$二[12162720丄3,0?5],用不同的方法求x=2點的插值,并分析所得結(jié)果有何不同.?x=[0.1,0.&l?3,l?9,2?5,3?l];y=[l?2,l?6,2?7,2?0,l?3,0?5];?p=interpl(x,y,2)L8833?x=[0?l,0?8,l?3,l?9,2.5,3?l];y二[l?
