拉格朗日插值法與牛頓插值法的比較_第1頁
拉格朗日插值法與牛頓插值法的比較_第2頁
拉格朗日插值法與牛頓插值法的比較_第3頁
拉格朗日插值法與牛頓插值法的比較_第4頁
拉格朗日插值法與牛頓插值法的比較_第5頁
已閱讀5頁,還剩2頁未讀 繼續免費閱讀

下載本文檔

版權說明:本文檔由用戶提供并上傳,收益歸屬內容提供方,若內容存在侵權,請進行舉報或認領

文檔簡介

1、拉格朗日插值法與牛頓插值法的比較摘 要在生產和科研中出現的函數是多樣的。對于一些函數很難找出其解析表達式。即使在某些情況下,可以寫出函數的解析表達式,但由于解析表達式的結構相當復雜,使用起來很不方便。插值法即是解決此類問題的一種古老的、然而卻是目前常用的方法,它不僅直接廣泛地應用于生產實際和科學研究中,而且也是進一步學習數值計算方法的基礎。拉格朗日插值法和牛頓插值法則是二種常用的簡便的插值法。本文即是討論拉格朗日插值法和牛頓插值法的理論及二者的比較。關鍵詞 拉格朗日插值 牛頓插值 插值多項式 比較一、 背景在工程和科學研究中出現的函數是多種多樣的。常常會遇到這樣的情況:在某個實際問題中,雖然可

2、以斷定所考慮的函數在區間上存在且連續,但卻難以找到它的解析表達式,只能通過實驗和觀測得到在有限個點上的函數值(即一張函數表)。顯然,要利用這張函數表來分析函數的性態,甚至直接求出其他一些點上的函數值可能是非常困難的。面對這些情況,總希望根據所得函數表(或結構復雜的解析表達式),構造某個簡單函數作為的近似。這樣就有了插值法,插值法是解決此類問題目前常用的方法。如設函數在區間上連續,且在個不同的點上分別取值。插值的目的就是要在一個性質優良、便于計算的函數類中,求一簡單函數,使 而在其他點上,作為的近似。通常,稱區間為插值區間,稱點為插值節點,稱式為插值條件,稱函數類為插值函數類,稱為函數在節點處的

3、插值函數。求插值函數的方法稱為插值法。插值函數類的取法不同,所求得的插值函數逼近的效果就不同。它的選擇取決于使用上的需要,常用的有代數多項式、三角多項式和有理函數等。當選用代數多項式作為插值函數時,相應的插值問題就稱為多項式插值。本文討論的拉格朗日插值法與牛頓插值法就是這類插值問題。在多項式插值中,最常見、最基本的問題是:求一次數不超過的代數多項式 使,其中,為實數。拉格朗日插值法即是尋求函數(拉格朗日插值多項式)近似的代替函數。相似的,牛頓插值法則是通過(牛頓插值多項式)近似的求得函數的值。二、 理論基礎(一)拉格朗日插值法在求滿足插值條件次插值多項式之前,先考慮一個簡單的插值問題:對節點中

4、任一點,作一n次多項式,使它在該點上取值為1,而在其余點上取值為零,即上式表明個點都是次多項式的零點,故可設其中,為待定系數。由條件立即可得故 由上式可以寫出個次插值多項式。我們稱它們為在個節點上的次基本插值多項式或次插值基函數。利用插值基函數立即可以寫出滿足插值條件的次插值多項式 根據條件,容易驗證上面多項式在節點處的值為,因此,它就是待求的次插值多項式。形如的插值多項式就是拉格朗日插值多項式,記為,即作為常用的特例,令,由上式即得兩點插值公式 ,這是一個線性函數,故又名線性插值。若令,則又可得到常用的三點插值公式這是一個二次函數,故又名二次插值或拋物插值。(二)牛頓插值法由線性代數知,任何

5、一個不高于次多項式,都可以表示成函數的線性組合。既可以吧滿足插值條件的次插值多項式寫成如下形式其中,為待定系數。這種形式的插值多項式稱為牛頓插值多項式,記為,即 因此,牛頓插值多項式是插值多項式的另一種表示形式。設函數在等距節點處的函數值為已知,其中是正常數,稱步長。我們稱兩個相鄰點和處函數之差為函數在點處以為步長的一階向前差分,記作,即于是,函數在各節點處的一階差分依次為又稱一階差分的差分為二階差分。一般的,定義函數在點處的階差分為。在等距節點情況下,可以利用差分表示牛頓插值多項式的系數。事實上,由插值條件可得;再由插值條件可得;一般的,由插值條件可得。于是,滿足插值條件的插值多項式為三、

6、二者的比較拉格朗日插值法與牛頓插值法都是二種常用的簡便的插值法。但牛頓法插值法則更為簡便,與拉格朗日插值多項式相比較,它不僅克服了“增加一個節點時整個計算工作必須重新開始”(見下面例題)的缺點,而且可以節省乘、除法運算次數。同時,在牛頓插值多項式中用到的差分與差商等概念,又與數值計算的其他方面有著密切的關系。現用一實例比較拉格朗日插值法與牛頓插值法例 已知函數表如下:x0.10.20.30.40.50.6sinx0.099830.198670.295520.389420.479430.56464計算sin(0.12)的值。利用拉格朗日插值法計算過程如下:(計算程序代碼見附件) 因為0.12位于

7、0.1與0.2之間,故取節點利用線性插值所求的近似值為計算結果如下圖利用拋物插值所求的近似值為計算結果如下圖利用牛頓插值法計算過程如下:構造差分表如下:xsinx0.10.20.30.40.099830.198670.295520.389420.098840.096850.09390-0.00199-0.00295-0.00096利用線性插值所求的近似值為利用拋物插值所求的近似值為從上面的計算過程可以看出,拉格朗日插值法的線性插值與拋物插值的計算過程沒有繼承性,即增加一個節點時整個計算工作必須重新開始。而牛頓插值則避免了這一問題,這樣大量的節省了乘、除法運算次數,減少了計算的時間。因此,對于一

8、些結構相當復雜的函數,牛頓插值法比拉格朗日插值法要占優勢。參考文獻1易大義,沈云寶,李有法編.計算方法.杭州:浙江大學出版社,20022馮康等編.數值計算方法.北京:國防工業出版社,19873李慶陽,王能超,易大義編.數值分析(第四版).北京:清華大學出版社,施普林格出版社,20014Burden R L,Faires J D,Reynolds A C. Numerical Analysis. Alpine Press,19815易大義,陳道琦編.數值分析引論.杭州:浙江大學出版社,1998 Comparison between Lagrange interpolation method an

9、d Newton interpolation methodAbstract In the production and scientific researches, there appears a variety of functions. For some function, it is difficult to find out its analytical expression. Though in some cases, the analytical expressions of the structure can be worked out, it is inconvenient t

10、o use them because of the complexity of structure. Interpolation method is a kind of old way to solve such problems, which is now commonly used. It is not only applied in the actual production or scientific researches directly and widely, but also become the foundation of further study of numerical

11、calculation method. Lagrange interpolation method and Newton interpolation law are two commonly used simple interpolation methods. This paper is a discussion of theory and the comparison between Lagrange interpolation method and Newton interpolation method. Key Words Lagrange interpolation ,Newton interpolation ,Interpolation polynomials,comparison 附件:#include <stdio.h>void main()float x6=0.1,0.2,0.3,0.4,0.5,0.6;int n,k,j;float f6=0.09983,0.19867,0.29552,0.38942,0.47943,0.56464;float p,a,sum=0;printf("輸入

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯系上傳者。文件的所有權益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網頁內容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
  • 4. 未經權益所有人同意不得將文件中的內容挪作商業或盈利用途。
  • 5. 人人文庫網僅提供信息存儲空間,僅對用戶上傳內容的表現方式做保護處理,對用戶上傳分享的文檔內容本身不做任何修改或編輯,并不能對任何下載內容負責。
  • 6. 下載文件中如有侵權或不適當內容,請與我們聯系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論