1. 程式人生 > >拉格朗日插值

拉格朗日插值

表達 有一個 rec med 做到 pan n+1 IT body

一,介紹

  學過FFT的人都應該知道什麽叫做插值,插值的意思就是說將一個多項式從點值表達轉變成系數表達。

  在FFT的插值中為什麽可以做到n log n,是因為單位復數根的關系。

  那對於普通的插值應該怎麽辦呢?解方程是一種方法,但是這個在計算機中十分不現實。

  所以有許多種插值的方法,其中比較普及的就是拉格朗日插值。

二,定義

   對某個多項式函數,已知有給定的k + 1個取值點:

技術分享圖片

   其中技術分享圖片對應著自變量的位置,而技術分享圖片對應著函數在這個位置的取值。

   假設任意兩個不同的xj都互不相同,那麽應用拉格朗日插值公式所得到的拉格朗日插值多項式為:

技術分享圖片

   其中每個技術分享圖片為拉格朗日基本多項式(或稱插值基函數),其表達式為:

技術分享圖片[3]

   拉格朗日基本多項式技術分享圖片的特點是在技術分享圖片上取值為1,在其它的點技術分享圖片上取值為0。

三,例子

   假設有某個二次多項式函數技術分享圖片,已知它在三個點上的取值為:

   技術分享圖片 技術分享圖片 技術分享圖片

   要求技術分享圖片的值。

   首先寫出每個拉格朗日基本多項式:

技術分享圖片
技術分享圖片
技術分享圖片

   然後應用拉格朗日插值法,就可以得到技術分享圖片的表達式(技術分享圖片為函數技術分享圖片的插值函數):

技術分享圖片
技術分享圖片
技術分享圖片

   此時代入數值技術分享圖片就可以求出所需之值:技術分享圖片

四,證明唯一性

  就是說n+1個點對應的n次多項式只有一個,這個需要證明,但是一般都會直接當成結論,所以不需要去記憶。

五,優點和缺點

  拉格朗日插值法的公式結構整齊緊湊,在理論分析中十分方便,然而在計算中,當插值點增加或減少一個時,所對應的基本多項式就需要全部重新計算,

  於是整個公式都會變化,非常繁瑣[5]。這時可以用重心拉格朗日插值法或牛頓插值法來代替。此外,當插值點比較多的時候,拉格朗日插值多項式的次數

  可能會很高,因此具有數值不穩定的特點,也就是說盡管在已知的幾個點取到給定的數值,但在附近卻會和“實際上”的值之間有很大的偏差(如右下圖)[6]。

  這類現象也被稱為龍格現象,解決的辦法是分段用較低次數的插值多項式。

拉格朗日插值