資源描述：
《實(shí)驗(yàn)報(bào)告 數(shù)據(jù)插值與擬合實(shí)驗(yàn)》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、實(shí)驗(yàn)報(bào)告§3.8數(shù)據(jù)插值與擬合實(shí)驗(yàn)?一、實(shí)驗(yàn)?zāi)康募耙饬x[1]了解插值、最小二乘擬合的基本原理[2]掌握用MATLAB計(jì)算一維插值和兩種二維插值的方法；[3]掌握用MATLAB作最小二乘多項(xiàng)式擬合和曲線擬合的方法。二、實(shí)驗(yàn)內(nèi)容1．針對實(shí)際問題，試建立數(shù)學(xué)模型。用MATLAB計(jì)算一維插值和兩種二維插值的方法求解；1．用MATLAB中的函數(shù)作一元函數(shù)的多項(xiàng)式擬合與曲線擬合，作出誤差圖；2．用MATLAB中的函數(shù)作二元函數(shù)的最小二乘擬合，作出誤差圖；3．針對預(yù)測和確定參數(shù)的實(shí)際問題，建立數(shù)學(xué)模型，并求解。三、實(shí)驗(yàn)步驟1．開啟軟件平臺(tái)——MAT
2、LAB，開啟MATLAB編輯窗口；2．根據(jù)各種數(shù)值解法步驟編寫M文件3．保存文件并運(yùn)行；4．觀察運(yùn)行結(jié)果(數(shù)值或圖形)；5．根據(jù)觀察到的結(jié)果寫出實(shí)驗(yàn)報(bào)告，并淺談學(xué)習(xí)心得體會(huì)。四、實(shí)驗(yàn)要求與任務(wù)根據(jù)實(shí)驗(yàn)內(nèi)容和步驟，完成以下具體實(shí)驗(yàn)，要求寫出實(shí)驗(yàn)報(bào)告（實(shí)驗(yàn)?zāi)康摹鷨栴}→數(shù)學(xué)模型→算法與編程→計(jì)算結(jié)果→分析、檢驗(yàn)和結(jié)論→心得體會(huì)）1.山區(qū)地貌：在某山區(qū)測得一些地點(diǎn)的高程如下表3.8。平面區(qū)域?yàn)椋?200<=x<=4000,1200<=y<=3600)試作出該山區(qū)的地貌圖和等高線圖，并對幾種插值方法進(jìn)行比較。表3.8某山區(qū)高程表yx120016
3、00200024002800320036004000120011301250128012301040900500700160013201450142014001300700900850200013901500150014009001100106095024001500120011001350145012001150101028001500120011001550160015501380107032001500155016001550160016001600155036001480150015501510143013001200980解
4、：x0=1200:400:4000;y0=1200:400:3600;z0=[1130,1250,1280,1230,1040,900,500,700;1320,1450,1420,1400,1300,700,900,850;1390,1500,1500,1400,900,1100,1060,950;1500,1200,1100,1350,1450,1200,1150,1010;1500,1200,1100,1550,1600,1550,1380,1070;1500,1550,1600,1550,1600,1600,1600,1550
5、;1480,1500,1550,1510,1430,1300,1200,980];meshc(x0,y0,z0)xlabel('x');ylabel('y');zlabel('z')title('原始圖')分段線性插值輸入命令：[xi2,yi2]=meshgrid(1200:10:4000,1200:10:3600);zi2=interp2(x0,y0,z0,xi2,yi2,'linear');meshc(xi2,yi2,zi2)xlabel('x');ylabel('y');zlabel('z')title('插值后')三次樣條插值
6、輸入命令：[xi,yi]=meshgrid(1200:10:4000,1200:10:3600);zi=interp2(x0,y0,z0,xi,yi,'spline');meshc(xi,yi,zi),title('三次樣條插值表面圖')xlabel('x');ylabel('y');zlabel('z')雙三次插值輸入命令：[xi,yi]=meshgrid(1200:10:4000,1200:10:3600);zi=interp2(x0,y0,z0,xi,yi,'cubic');meshc(xi,yi,zi),title('雙三次插
7、值后的表面圖')xlabel('x軸');ylabel('y軸');zlabel('z軸')畫等高線圖：n=9;pcolor(x0,y0,z0)shadinginterpzmax=max(max(z0));zmin=min(min(z0));caxis([zmin,zmax])colorbarholdonC=contour(x0,y0,z0,n,'k:');clabel(C)title('未擬合圖形')holdoffpcolor(xi,yi,zi)shadinginterpzmax1=max(max(zi));zmin1=min(mi
8、n(zi));caxis([zmin1,zmax1])colorbarholdonC=contour(xi,yi,zi,n,'k:');clabel(C)title('擬合圖形')holdoff相應(yīng)的圖形中都畫出了地貌圖
