拉格朗日插值演算法分析

2021-07-23 23:40:57 字數 1272 閱讀 5905

這幾天一直研究拉格朗日多項式,今天將自己對拉格朗日多項式的理解寫在這裡,方便大家交流。

在數值分析中,拉格朗日常用於多項式插值。假定提供一組資料點[xi,yi],拉格朗日插值多項式就是由這些

資料的線性運算得到的。

其中基本的多項式有以下公式計算得到

注意1.第一提供的xi應該是沒有相同的,否則不能應用此演算法

2.對於每乙個xi,yi都只對應乙個值

3.對於i~=j,lj(x)包括x-xi,因此,整個表示式的乘積將是0。

4.相反

5.因此yi*li(x)=yi,因此在每乙個點xi,l(xi) = yi+0+0+........+0=yi,表明這個拉格朗日插值演算法的正確性。

舉例說明

x = [1 2 3 ]  y = f(x)=[1 4 9]

則l(x) = 1*(x-2)*(x-3)/+4*(x-1)(x-3)/+9(x-1)(x-2)/=x^2

下面看一下matlab的具體**實現過程。

雖然樣本資料點只有21個,但足以看出插值曲線足以和滿足標準了,插值效果還是很好的

拉格朗日插值

拉格朗日插值基函式 li x x x 0 x xi 1 x xi 1 x x n x i x0 xi xi 1 xi xi 1 xi xn 拉格朗日差值函式 ln x i 0 nyil i x 其中,x為缺失值對應的下表序號,ln x 為缺失值的插值結果,xi 為缺失值yi 的下表序號。對全部缺失值...

拉格朗日插值

function s larg1 x,y,xi m length x 求出插值節點向量長度 n length y if m n error 向量x與y的長度必須一致 這裡肯定一致,只是為了消除直接選取x,y資料的時候出錯而設定的 end s 0 for i 1 n z ones 1,length x...

拉格朗日插值

存在性和唯一性的證明以後再補。拉格朗日插值,emmmm,名字挺高階的 joy 它有什麼應用呢?我們在fft中講到過 設 n 1 次多項式為 y sum a i x i 有乙個顯然的結論 如果給定 n 個互不相同的點 x,y 則該 n 1 次多項式被唯一確定 那麼如果給定了這互不相同的 n 個點,利用...