在某些插值問題中,不僅給出了插值點的值,也給出了插值點的導數值,這便可以運用厄公尺特插值公式進行插值。
厄公尺特插值公式如下:
下面給出**:
//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...