第3講牛頓插值公式

2023-01-07 06:48:06 字數 1911 閱讀 6451

§1.4 差商與差分及其性質

1 差商的概念:

稱為函式f(x)的一階差商;

稱為函式f(x)的二階差商;

一般地,稱為函式f(x)的n階差商;

特別地,定義為函式f(x)關於xo的零階差商。

由此可知,高階差商總是由比它低一階的的兩個差商組合而成。

2 差商性質

(a) 性質1:n階差商可以表示成n+1個函式值的線性組合,即

該性質說明:k階差商計算是由函式值f(x0),f(x1),…f(xk)線性組合而。

如:;(b) 性質2(對稱性): 差商與節點的順序無關。即

,這一點可以從性質1看出。

3 利用差商表計算差商

利用差商的遞推定義,可以用遞推來計算差商。

差商表:

如要計算四階差商,應再增加乙個節點,表中還要增加一行。

4 差分的概念

定義設函式y=f(x)在等距節點上的函式值f(xi)=fi,其中,h為常數稱作步長。稱

△fi=fi+1-fi

▽fi=fi-fi-1

δfi=f(xi+h/2)-f(xi-h/2)=

分別為f(x)在處以h為步長的一階向前差分,一階向後差分和一階中心差分。稱符號△、▽、δ分別為向前差分運算元,向後差分運算元和中心差分運算元。

在節點等距情況下,差商可用差分表示,設步長,有

一般形式(數學歸納法可證)

§1.5 牛頓插值公式

1. 牛頓插值公式的構造

lagrange 插值雖然易算,但若要增加乙個節點時,全部基函式 li(x) 都需重新算過。本節介紹另外一種方法-牛頓插值法,並用它解決上面所述問題。

由線性插值

,令二次插值能否寫成

由條件得

推廣得,

其中,為待定係數。如何求?

解: 因為,

所以 (0)

又,有 (1)又2)

一般地,

(n)將式(n)代入式(n-1), ...,式(2)代入式 (1),式(1)代入式 (0),

如此可得:

尤為注意的是:最後一項中,差商部分含有,乃是餘項部分,記作;而前面n+1項中,差商部分都不含有,因而前面n+1項是關於的n次多項式,記作,這就是牛頓插值公式。

2 算例

例1:當n=1時,

,其中,

這就是牛頓一次插值多項式,也就是點斜式直線方程。

當n=2時,

這就是牛頓二次插值多項式。顯然,

,。即滿足二次插值條件。

例2: 已知

求滿足以上插值條件的牛頓型插值多項式。

解: 由於:

,,,;

則牛頓三次插值多項式為

。1 拉格朗日插值與牛頓插值的比較

(1)和均是n次多項式,

且均滿足插值條件:

。由多項式的唯一性,,因而,兩個公式的餘項是相等的,即

(2)當插值多項式從n-1次增加到n次時,拉格朗日型插值必須重新計算所有的基本插值多項式;而對於牛頓型插值,只需用**再計算乙個n階差商,然後加上一項即可。

4 等距牛頓插值公式

插值節點為等距節點:

,,如下圖:

h h hh

牛頓插值公式

設等距節點,

記.當,令,. 例如(下圖)

x在x2,x3的中點時,。

將牛頓插值公式中的差商用差分代替,而

從而,牛頓插值公式在等距插值節點下的形式為:

餘項為這是等距牛頓向前插值公式。

例4: 設插值節點為,相應的函式值如下表,求f(2.2)。

解:精確值f(2.2)=e2.2=9.025011。

此時[xk, xk+1],x=2.2=1+2.4h 故t=2.4,於是

求時,在後加一項:

,所以求時,在後再加一項:,所以

實驗1拉格朗日插值與牛頓插值

西華數學與計算機學院上機實踐報告 一 目的 1 通過本實驗加深對拉格朗日插值和牛頓插值法構造過程的理解 2 能對上述兩種插值法提出正確的演算法描述程式設計實現。二 內容與設計思想 編制乙個程式,分別用拉格朗日插值法和牛頓插值法求解某點的函式近似值。已知y f x 的資料表如下,求t 0.63處的函式...

第3講遞推

遞推是一種應用非常廣泛的常用演算法之一,與下一章的遞迴有著密切的聯絡。本章 遞推在求解數列 數陣以及計數等應用案例方面的應用。在紛繁變幻的世界,所有事物都隨時間的流逝發生著微妙的變化。許多現象的變化是有規律可循的,這種規律往往呈現出前因後果的關係。某種現象的變化結果與緊靠它前面變化的乙個或一些結果緊...

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

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