全區間厄公尺特插值演算法

2021-07-25 04:37:53 字數 833 閱讀 5287

在某些插值問題中,不僅給出了插值點的值,也給出了插值點的導數值,這便可以運用厄公尺特插值公式進行插值。

厄公尺特插值公式如下:

下面給出**:

//hermite插值

/* *函式宣告double hermit(double* x,double* y,double* dy,int n,double t)

*x為存放給定n個結點的指標

*y為存放n個節點上的值的指標

*dy為n個結點的導數的指標

*n為指定插值點的個數

*t為求值的插值點

*函式的返回值是double型別

*/#include

using

namespace

std;

double hermite(double*x, double* y, double* dy, int n, double t)

return z;

}int main()

; double y[10] = ;

double dy[10];

for (int i = 0; i < 10; i++)

dy[i] = -y[i];

double t = 0.356;

double z = hermite(x, y, dy, 10, t);

cout

<< "插值點的值為 "

<< z << endl;

return

0;}

MATLAB之埃公尺特插值

一 演算法原理 1 問題引入 在實際問題中,對所構造的插值多項式,不僅要求函式值重合,而且要求若干階導數也重合,即要求插值函式 p x 滿足 此類問題被稱為埃公尺特插值,相應的插值多項式成為埃公尺特插值多項式。2 考慮函式值與導數值個數相等的情況 設在節點a x0要求插值多項式h x h2n 1 x...