插值法和擬合實驗報告(數(shù)值計算)

1、=1=1=1=1插值法和擬合實驗報告一、實驗目的通過進行不同類型的插值，比較各種插值的效果，明確各種插值的優(yōu)越性通過比較不同次數(shù)的多項式擬合效果，了解多項式擬合的原理；利用matlab編程，學會matlab命令；掌握拉格朗日插值法；掌握多項式擬合的特點和方法。二、實驗題目、插值法實驗x將區(qū)間-5,510等分，對下列函數(shù)分別計算插值節(jié)點Xk的值，進行不同類型的插值，作出插值函數(shù)的圖形并與yf(X)的圖形進行比較：f(f(x)arctanx;f(x)做拉格朗日插值做分段線性插值做三次樣條插值.2、擬合實驗給定數(shù)據(jù)點如下表所示：X-1.5-1.0-0.50.00.51.01.5yi1-4.4-0.4

2、0.550.05-0.40.544.55分別對上述數(shù)據(jù)作三次多項式和五次多項式擬合，并求平方誤差，作出離散函數(shù)(X，yi)和擬合函數(shù)的圖形。三、實驗原理、插值法實驗L(x)=工l(x)y=l(x)yl(x)yniio0nni=0i-1設L(x)=工l(x)y=c(x一x)(x一x)(x一x)(x一x)=c(x一i-1ci+1nj再由,(x-x)=1,ji=0,1,n=1=1=1=1=1=1=1=11c=(x-x)j故,l(xl(x)=iix一xix一xij=1=1=1=12、擬合實驗取基函數(shù)X）=1；卩二乂*孤（x）=工枕由于財仇心W廣1（于y3=1程組把上式代人下面的法方程組0s得到0s得到

3、陽仏、心?。②郌卩J%（產卩L）夠0=1=1=1=1=1=1=1=1E兒11!=111Z羽兒i=inVxnVxm.i=inV:r2=1n/:ri=1)1)四、實驗內容1.、插值法實驗1.1實驗步驟：打開matlab軟件，新建一個名為chazhi.m的M文件，編寫程序(見1.2實驗程序)，運行程序，記錄結果。1.2實驗程序：x=-5:1:5;xx=-5:0.05:5;y1=1./(1+x.2);L=malagr(x,y1,xx);L1=interp1(x,y1,x,linear);S=maspline(x,y1,0.0148,-0.0148,xx);holdon;plot(x,y1,b*);p

4、lot(xx,L,r);plot(x,L1,g);plot(xx,S,k);figurex=-5:1:5;xx=-5:0.05:5;y2=atan(x);L=malagr(x,y2,xx);L1=interp1(x,y2,x,linear);S=maspline(x,y2,0.0385,0.0385,xx);holdon;plot(x,y2,b*);plot(xx,L,r);plot(x,L1,g);plot(xx,S,k);figurex=-5:1:5;xx=-5:0.05:5;y3=x.2./(1+x.4);L=malagr(x,y3,xx);L1=interp1(x,y3,x,linea

5、r);S=maspline(x,y3,0.0159,-0.0159,xx);holdon;1)1)plot(x,y3,b*);plot(xx,L,r);plot(x,L1,g);plot(xx,S,k);1.3實驗設備：matlab軟件。2、擬合實驗2.1實驗步驟：新建一個名為nihe.m的M文件，編寫程序(見2.2實驗源程序)，運行程序，記錄結果。2.2實驗程序：x=-1.5-1.0-0.500.51.01.5;y=-4.45-0.450.550.05-0.440.544.55;a1=mafit(x,y,3)x1=-1.5:0.05:1.5;y1二a1(4)+a1(3)*x1+a1(2)*x

6、1.2+a1(1)*x1.3;holdonplot(x,y,b*);plot(x1,y1,r);p1=polyval(a1,x);s1=norm(y-p1)figurea2=mafit(x,y,5)x2=-1.5:0.05:1.5;y2二a2(6)+a2(5)*x2+a2(4)*x2.2+a2(3)*x2.3+a2(2)*x2.4+a2(1)*x2.5;holdonplot(x,y,b*);plot(x2,y2,r);p2=polyval(a2,x);s2=norm(y-p2)2.3實驗設備：matlab軟件。五、實驗結果1.、插值法實驗(2)(3)(1)(1)(1)(1)2、擬合實驗(2)(

7、2)(2)(2)平方誤差：si=0.0136輸入程序得到：al=2.0000-0.0014-1.50070.0514si=0.0136平方誤差：s2=0.0069輸入程序得到：a2=0.01200.00481.9650-0.0130-1.48200.0545s2=0.0069六、實驗結果分析1、插值法實驗結果分析：（1）由插值結果曲線圖可見，拉格朗日插值在節(jié)點附近誤差很小，但在兩端有振蕩現(xiàn)象；分段線性插值具有良好的收斂性，但在節(jié)點處不光滑；而三次樣條插值在直觀上與原函數(shù)曲線吻合得最好；（2）分析可知，均勻插值時（拉格朗日插值），會出現(xiàn)多項式插值的Runge現(xiàn)象，當進行非等距節(jié)點插值時（分段線性插值、三次樣條插值），其近似效果明顯要比均勻插值要好，原因是非均勻插值時，在遠離原點處的插值節(jié)點比較密集，所以其插值近似