Matlab實(shí)驗(yàn)報(bào)告六(三次樣條與分段線性插值)

1、數(shù)學(xué)與信息科學(xué)系實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)名稱插值與擬合所屬課程數(shù)學(xué)軟件與實(shí)驗(yàn)實(shí)驗(yàn)類型綜合型實(shí)驗(yàn)專業(yè)信息與計(jì)算科學(xué)班級(jí)學(xué)號(hào)姓名1扌旨導(dǎo)教師 ODDDOO【實(shí)驗(yàn)?zāi)康摹繉W(xué)會(huì)在環(huán)境下使用幾種不同的插值法和擬合兩種方法構(gòu)造函數(shù)依據(jù)已經(jīng)知道的某些特殊點(diǎn)來推測(cè)實(shí)際問題中需要知道但又不便于測(cè)量出來的量?！緦?shí)驗(yàn)原理】z=interp2(x0,y0,z0,x,y,method)yx,yn;n3.linearspline【實(shí)驗(yàn)環(huán)境】MatlabR2010bODDDOO問題對(duì)函數(shù)x0,y0 xyxx0,y0,X,5,5,分別用分段線性插值和三次樣條插值作插值(其中插值節(jié)點(diǎn)不少于20),并分別作出每種插值方法的誤差曲線1.2.x=

2、-5:0.5:5;y=1./(1+x.入2);x1=-5:0.01:5;yl=interp1(x,y,x1,linear);ys=interp1(x,y,x1,spline);plot(x,y,*,x1,yl,r,x1,y2,b)y0=1./(1+x1.入2);y3=yl_y0;y4=ys_y0;holdonplot(x1,y3,y,x1,y4,g) 3.-4-3-2-1-54.問題山區(qū)地貌圖在某山區(qū)（平面區(qū)域口0,2800）（0,2400）內(nèi)，單位：米）測(cè)得一 #些地點(diǎn)的高程（單位：米）如表1,試作出該山區(qū)的地貌圖表24001430145014701320128012001080940200

3、0145014801500155015101430130012001600146015001550160015501600160016001200137015001200110015501600155013808001270150012001100135014501200115040012301390150015001400900110010600118013201450142014001300700900Y/X040080012001600200024002800 # 1.2. x=0:400:2800;y=0:400:2400;z=1430145014701320128012001080940;14501480150015501510143013001200;14601500155016001550160016001600;13701500120011001550160015501380;12701500120011001350145012001150;1230139015001500140090011001060;118013201450142014001300700900;mesh(x,y,z)xi=0:20:2800;yi=0:20:2400;zi=interp2(x,y,z,xi,yi,cubic);mesh(xi,yi,zi)3.1800-16001400-12