拉格朗日插值程式設計實現

2021-09-07 04:16:16 字數 962 閱讀 8190

拉格朗日插值原理:

拉格朗日插值的具體介紹**:

翻譯**話就是,該曲線是由多個n次多項式的和構成的,n是參與插值的點的個數。每個n次多項式的計算方法如上圖所示。轉化成程式的話,就是要儲存每個多項式中(x-xi)中的每一項xi。然後就是係數yi/(xj-x0)(xj-x1)...(xj-xk)

程式實現:

#include #include 

using

namespace

cv;template

struct

point

point(t a,t b)

point

operator+(point p)

};/*

拉格朗日曲線,包含多個拉格朗日多項式和曲線的x座標最大最小值 */

struct

lagrangecurve

;void callagrangepolynomial(lagrangecurve&curve,point*curvepoints,int

pointnum)}}

}double getinterpolation(int x,lagrangecurve&lagrangecurve)

}y+=temp;

}returny;}

int main(int argc, const

char *argv)

imshow(

"img

", img);

waitkey(-1

);

return0;

}

上述**使用到了opencv畫圖,執行結果圖如下:

拉格朗日插值

拉格朗日插值基函式 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 個點,利用...