pre define

typedef struct Operators {std::vector<float> Operatos;}Operators;// 定義一個點 class PointF { public:PointF(){X = Y = 0.0f;}PointF(IN const PointF &point){X = point.X;Y = point.Y;}PointF(IN float x,IN float y){X = x;Y = y;}PointF operator+(IN const PointF& point) const{return PointF(X + point.X,Y + point.Y);}PointF operator-(IN const PointF& point) const{return PointF(X - point.X,Y - point.Y);}BOOL Equals(IN const PointF& point){return (X == point.X) && (Y == point.Y);}public:float X;float Y; };

fit lagrange

void CWinaTVWaveformDlg::fitLagrange(std::vector<PointF> pointList) {std::vector<Operators> OpList;// result std::vector<float> Xs;std::vector<float> Ys;try{if (pointList.size() > 0){//compute lagrange operator for each X coordinatefor (int x = 1; x < 2000; x++){//list of float to hold the Lagrange operators, Init the list with 1'sstd::vector<float> L(pointList.size(), 1); for (int i = 0; i < L.size(); i++){for (int k = 0; k < pointList.size(); k++){if (i != k)L[i] *= (float)(x - pointList[k].X) / (pointList[i].X - pointList[k].X);}}Operators o;o.Operatos = L;OpList.push_back(o);Xs.push_back(x);}//Computing the Polynomial P(x) which is y in our curvestd::vector<Operators>::iterator iter;for ( iter = OpList.begin(); iter != OpList.end(); iter++){float y = 0;for (int i = 0; i < pointList.size(); i++){y += iter->Operatos[i] * pointList[i].Y;}Ys.push_back(y);}//Drawing the curve in the simplest way}else{AfxMessageBox(_T("Lagrange curve fitting, add some points"));}}catch (CException* e){return;} }

總結(jié)

以上是生活随笔為你收集整理的拉格朗日插值拟合的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。