在離散資料的基礎上補插出連續函數。是計算數學中最基本和常用的手段,是函數逼近的重要方法。利用它可通過函數在有限個點處的取值狀況,估算該函數在別處的值。早在西元6世紀,中國劉焯已將等距二次插值法用於天文計算。17世紀,I.牛頓和J.格雷果裏建立瞭等距結點上的一般插值公式。18世紀,J.-L.拉格朗日給出瞭更一般的非等距結點上的插值公式。在近代,插值法是觀測資料處理和函數製表所常用的工具,又是導出其他許多數值方法(例如數值積分、非線性方程求根、微分方程數值解等等)的依據。
插值問題的提法是:假定已知區間[α,b]上的實值函數 f(x)在該區間中n+1個互不相同的點x0,x1,…,xn處的值是f(x0),f(x1),…,f(xn),要求估算f(x)在[α,b]中某點x=x處的值。插值的作法是:在事先選定的一個由簡單函數所構成的含n+1個參數C0,C1,…,Cn的函數類φ(C0,C1,…,Cn)中求出滿足條件
![](/img3/2098.gif)
的函數p(x),並以p(x)作為 f(x)的估值。此處,函數f(x)稱為被插函數;x0,x1,…,xn稱為插值結點;φ(C0,C1,…,Cn) 稱為插值函數類;式(1)稱為插值條件。φ(C0,C1,…,Cn)中滿足插值條件(1)的函數p(x)稱為插值函數。誤差函數
![](/img3/2099.gif)
多項式插值 插值函數類取成代數多項式類的情形,是最常用的一種插值。此時對[α,b]上的任何實值函數f(x)都相應地有惟一的次數不超過n多項式p(x) 滿足插值條件(1)。p(x)稱為f(x)的插值多項式。當f(x)在[α,b]上n+1次可微時,插值餘項為
![](/img3/2100.gif)
下面兩種插值公式是p(x)的具體表達式:
① 拉格朗日插值公式
![](/img3/2101.gif)
![](/img3/2102.gif)
![](/img3/2103.gif)
![](/img3/2104.gif)
② 牛頓插值公式
![](/img3/2105.gif)
式中f(x0,x1,…,xk)為函數f(x)在點x0,x1,…,xk上的k階差商(或均差)。各階差商由下列遞推方式定義:
![](/img3/2106.gif)
![](/img3/2107.gif)
拉格朗日插值公式和牛頓插值公式是同一插值多項式p(x)的不同表現形式。前者結構緊湊、意義清晰和便於理論分析;後者在實際計算時較為方便:若要增加新的插值結點,隻需相應地添加新的項即可。
對於等距的插值結點,即當
![](/img3/2108.gif)
![](/img3/2109.gif)
此處Δk表示步長為h的k階差分算子,其定義是:
![](/img3/2110.gif)
埃爾米特插值 插值條件帶微商的插值,其插值條件為
![](/img3/2111.gif)
![](/img3/2112.gif)
![](/img3/2113.gif)
古典的有限項泰勒展開式也可看作是埃爾米特插值多項式,其所有的插值條件都加在一個點上。反之,也可將埃爾米特插值多項式看作是多中心泰勒展開式。
分段插值 在實用中很少采用高次多項式插值(例如7、8次以上的多項式插值),因為在被插函數不夠光滑或插值結點選擇不當時,高次插值多項式常常在被插函數附近激烈地擺動,不能逼近被插函數;再者,高次插值多項式常常將插值條件的數據中含有的誤差過分地放大和擴散。因此,在實用中,往往是先將全區間分成許多小區間,然後在每個小區間上,采用低次插值(例如一次、二次或三次插值)。通常稱這樣的方法為分段插值法。如將小區間的端點取為插值結點,則相鄰區間的兩插值多項式在公共結點處將取相同的值,即兩段多項式曲線在公共結點處銜接。實踐表明,用分段的低次插值多項式逼近被插函數往往比在全區間上用高次插值多項式逼近效果好。
樣條插值 一種非局部性的分段插值。在每個分段點處相鄰插值曲線段的銜接具有一定的光滑度。最常用的一種是三次多項式樣條插值:給定結點x0<x1<…<xn和想應的函數值f(x0),f(x1),…,f(xn),取內結點x1,x2,…,xn-1作為分段點,尋求一個插值函數S(x),它在每個小區間[xi-1,xi](i=1,2,…,n)上分別都是三次多項式,在結點處滿足插值條件
![](/img3/2114.gif)
![](/img3/2115.gif)
三角插值 插值函數類為三角多項式類的情形。當被插函數f(x)是以2π為周期的函數時,通常用n階三角多項式
![](/img3/2116.gif)
![](/img3/2117.gif)
![](/img3/2118.gif)
有理插值 插值函數類取成有理函數類的情形。當被插函數具有極點或其他奇點時,用有理插值往往很有效。假定已知被插函數f(x)在m+n+1個互不相同的點xn,x1,…,xm+n處的值為f(x0),f(x1),…,f(xm+n)。有理插值就是求一個形如
![](/img3/2119.gif)
![](/img3/2120.gif)
![](/img3/2121.gif)
![](/img3/2122.gif)
多元插值 一元插值法的多元推廣。其插值函數類可取為多元多項式,多元三角函數,多元樣條,多元有理分式等。
參考書目
胡祖熾編:《計算方法》,高等教育出版社,北京,1959。
J.F.Steffensen,Interpolation,2nd ed.,Chelsea,New York,1950.
P.J.Davis,Interpolation and ApproxiMation,Blaisdell,New York,1963.