資源描述:
《用MATLAB解插值和曲線擬合問題.doc》由會員上傳分享,免費在線閱讀,更多相關(guān)內(nèi)容在教育資源-天天文庫。
1、注:上機作業(yè)文件夾以自己的班級姓名學(xué)號命名,文件夾包括如下上機報告和Matlab程序。上機報告模板如下:佛山科學(xué)技術(shù)學(xué)院上機報告課程名稱數(shù)學(xué)應(yīng)用軟件上機項目用MATLAB解插值和曲線擬合問題專業(yè)班級一.上機目的熟練掌握多種插值方法:線性插值,三次樣條插值,三次插值和最近鄰點插值(linear、spline、cubic、nearest)Interp1(x,y,cx,’method’)在MATLAB的線性最小二乘擬合中,用得較多的是多項式擬合,其命令為:A=polyfit(x,y,m)二.上機內(nèi)容1、在飛機的機翼加工時,由于機翼尺寸很大,通常在圖紙上只能標(biāo)出部分關(guān)鍵點的數(shù)
2、據(jù).某型號飛機的機翼上緣輪廓線的部分數(shù)據(jù)如下。用插值法求x每改變0.1時的y值,畫出圖形表示。x04.749.051938577695114133152171190y05.238.111.9716.1517.116.3414.6312.166.697.033.9902、已知觀測數(shù)據(jù)點如表所示:xy0-0.4470.11.9780.23.280.36.160.47.080.57.340.67.660.79.560.89.480.99.3111.2用3次多項式函數(shù)擬合這些數(shù)據(jù)點,畫出圖形。3、教材習(xí)題1,用griddata插值函數(shù),三次插值。4、教材習(xí)題4,初始值為V0=
3、0,=0.05三.上機方法與步驟第1題:要用插值法求x每改變0.1時的y值,可以選用程序y=interp1(xheng,yzhong,x,'spline');xheng=[04.749.051938577695114133152171190];yzhong=[05.238.111.9716.1517.116.3414.6312.166.697.033.990];x=0:0.1:190;y=interp1(xheng,yzhong,x,'spline');plot(xheng,yzhong,'+',x,y,xheng,yzhong,'r:')xlabel('xhen')
4、,ylabel('yzhon')截圖:第2題:要用3次多項式函數(shù)擬合這些數(shù)據(jù)點,并畫出圖形,選用程序A=polyfit(x,y,3),z=polyval(A,x);plot(x,y,'k+',x,z,'r')編程求解即可。x=0:0.1:1;y=[-0.4471.9783.286.167.087.347.669.569.489.311.2];A=polyfit(x,y,3)z=polyval(A,x);plot(x,y,'k+',x,z,'r')截圖:第3題:用griddata插值函數(shù),三次插值,選用cK=griddata(T,P,K,cT,cP,'cubic')編程
5、。T=[68688787106106140140];P=[9.798113.3249.007813.3559.791814.2779.656312.463];K=[0.08480.08970.07620.08070.06960.07530.06110.0651];cT=99;cP=10.3;cK=griddata(T,P,K,cT,cP,'cubic')截圖:第4題:要確定t,v,選用x=lsqcurvefit('curvefun4',x0,t,v),f=curvefun4(x,t)functionf=curvefun4(x,t)f=10-(10-x(1))*exp(
6、-t/x(2));t=[0.51234579];v=[6.366.487.268.228.668.999.439.63];x0=[0,0.05];x=lsqcurvefit('curvefun4',x0,t,v)f=curvefun4(x,t)截圖:四.上機結(jié)果第1題:截圖:第2題:A=Columns1through316.0758-33.924529.3246Column4-0.6104截圖:第3題:cK=0.0724截圖:第4題:x=5.55773.5002f=Columns1through56.14906.66167.49138.11478.5832Column
7、s6through88.93539.39879.6604截圖: