插值法和擬合實驗報告 數值計算

2022-06-13 16:12:02 字數 2850 閱讀 2089

一、 實驗目的

1.通過進行不同型別的插值,比較各種插值的效果,明確各種插值的優越性;

2.通過比較不同次數的多項式擬合效果,了解多項式擬合的原理;

3.利用matlab程式設計,學會matlab命令;

4.掌握拉格朗日插值法;

5.掌握多項式擬合的特點和方法。

二、 實驗題目

1.、插值法實驗

將區間[-5,5]10等分,對下列函式分別計算插值節點的值,進行不同型別的插值,作出插值函式的圖形並與的圖形進行比較:

(1) 做拉格朗日插值;

(2) 做分段線性插值;

(3) 做三次樣條插值.

2、擬合實驗

給定資料點如下表所示:

分別對上述資料作三次多項式和五次多項式擬合,並求平方誤差,作出離散函式和擬合函式的圖形。

三、 實驗原理

1.、插值法實驗

2、擬合實驗

四、 實驗內容

1.、插值法實驗

1.1實驗步驟:

開啟matlab軟體,新建乙個名為的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);

hold on;

plot(x,y1,'b*');

plot(xx,l,'r');

plot(x,l1,'g');

plot(xx,s,'k');

figure

x=-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);

hold on;

plot(x,y2,'b*');

plot(xx,l,'r');

plot(x,l1,'g');

plot(xx,s,'k');

figure

x=-5:1:5;

xx=-5:0.05:5;

y3=x.^2./(1+x.^4);

l=malagr(x,y3,xx);

l1=interp1(x,y3,x,'linear');

s=maspline(x,y3,0.0159,-0.0159,xx);

hold on;

plot(x,y3,'b*');

plot(xx,l,'r');

plot(x,l1,'g');

plot(xx,s,'k');

1.3實驗裝置: matlab軟體。

2、擬合實驗

2.1.實驗步驟:

新建乙個名為的m檔案,編寫程式(見2.2實驗源程式),執行程式,記錄結果。

2.2實驗程式:

x=[-1.5 -1.0 -0.5 0 0.5 1.0 1.5];

y=[-4.45 -0.45 0.55 0.05 -0.44 0.54 4.55];

a1=mafit(x,y,3)

x1=[-1.5:0.05:1.5];

y1=a1(4)+a1(3)*x1+a1(2)*x1.^2+a1(1)*x1.^3;

hold on

plot(x,y,'b*');

plot(x1,y1,'r');

p1=polyval(a1,x);

s1=norm(y-p1)

figure

a2=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;

hold on

plot(x,y,'b*');

plot(x2,y2,'r');

p2=polyval(a2,x);

s2=norm(y-p2)

2.3實驗裝置: matlab軟體。

五、 實驗結果

1.、插值法實驗

(1)(2)

(3)2、擬合實驗

(1)平方誤差:s1 =0.0136

輸入程式得到:

a1 =

2.0000 -0.0014 -1.5007 0.0514

s1 =

0.0136

(2)平方誤差:s2 = 0.0069

輸入程式得到:

a2 =

0.0120 0.0048 1.9650 -0.0130 -1.4820 0.0545

s2 =

0.0069

>>六、 實驗結果分析

1.、插值法實驗

結果分析:

(1) 由插值結果曲線圖可見,拉格朗日插值在節點附近誤差很小,但在兩端有振盪現象;分段線性插值具有良好的收斂性,但在節點處不光滑;而三次樣條插值在直觀上與原函式曲線吻合得最好;

(2)分析可知,均勻插值時(拉格朗日插值),會出現多項式插值的runge現象,當進行非等距節點插值時(分段線性插值、三次樣條插值),其近似效果明顯要比均勻插值要好,原因是非均勻插值時,在遠離原點處的插值節點比較密集,所以其插值近似效果要比均勻插值時的效果要好。

2、擬合實驗

結果分析:

可能是原始資料太少的問題,在擬合結果曲線圖看不出三次和五次有什麼差別,但由於三次多項式擬合的平方誤差大於五次多項式擬合的平方誤差,因此五次多項式擬合比三次多項式擬合效果要好。

牛頓插值法的實驗報告,還有實驗流程圖

資訊工程與自動化學院學生實驗報告 2005 2006 學年第1 學期 課程名稱 計算方法實驗室 2301 2005年10 月日 一,實驗目的及內容 1 實驗目的 1 體會並了解拉格朗日插值法。2 體會並了解牛頓插值法。2 實驗內容 1 拉格朗日插值法 2 牛頓插值法 二,程式流程圖 1 拉格朗日插值...

數值計算 分析 實驗報告

南昌航空大學 數學與資訊科學學院 實驗報告 課程名稱 數值計算方法 實驗名稱 曲線擬合 實驗型別 驗證性 綜合性 設計性 實驗室名稱 數學實驗室 班級學號 09072113 學生姓名邢憲平 任課教師 教師簽名 成績實驗日期 2012 3 25 一 實驗目的 實驗目的 實驗目的 了解函式逼近與曲線擬合...

拉格朗日插值實驗報告

實驗名稱 實驗一拉格朗日插值 我們在生產生活中常常會遇到這樣的問題 某個實際問題中,函式f x 在區間 a,b 上存在且連續,但卻很難找到其表示式,只能通過實驗和觀測得到有限點上的函式表。顯然,根據這些點的函式值來求其它點的函式值是非常困難的。有些情況雖然可以寫出表示式,但結構複雜,使用不方便。所以...