埃尔米特(Hermite)插值.docVIP

实验二 埃尔米特（Hermite）插值 掌握埃尔米特插值算法原理； 使用C语言编程实现埃尔米特插值算法。 二、实验准备： 阅读《数值分析》2.4节 二、实验要求： 某人从甲地开车去乙地，每隔一段时间对行车距离和速率进行一次采样，得到在 n+1 个采样时刻点 ti 的里程 si 和速率 vi（i=0, 1, ..., n）。要求编程构造埃尔米特插值多项式 H2n+1(t)，满足 H2n+1(ti)=si，H2n+1(ti)=vi，对所有i=0, 1, ..., n成立，并据此计算m个给定时刻的里程和速率。 函数接口定义： void Hermite_Interpolation( int N, double t[], double s[], double v[], int m, double ht[], double hs[], double hv[] ); 其中N为采样点个数（注意这个N不是公式中的最大下标n，而是等于n+1），采样时刻点ti、里程si、速率vi分别通过t、s、v传入；m是需要估算的给定时刻的个数，ht传入给定的时刻点，相应计算出的里程和速率应分别存储在hs和hv中。 裁判程序如下： #includestdio.h#define MAXN 5 /* 最大采样点个数 */#define MAXM 10 /* 最大估算点个数 */void Hermite_Interpolation( int N, double t[], double s[], double v[], int m, double ht[], double hs[], double hv[] );int main(){ int N, m; double t[MAXN], s[MAXN], v[MAXN]; /* 用于构造 的数据 */ double ht[MAXM], hs[MAXM], hv[MAXM]; /* 用 估算的数据 */ int i;  while ( scanf(%d, N) != EOF ) { for ( i=0; iN; i++ ) scanf(%lf, t[i]); for ( i=0; iN; i++ ) scanf(%lf, s[i]); for ( i=0; iN; i++ ) scanf(%lf, v[i]); scanf(%d, m); for ( i=0; im; i++ ) scanf(%lf, ht[i]); Hermite_Interpolation( N, t, s, v, m, ht, hs, hv ); for ( i=0; im; i++ ) printf(%.4lf , hs[i]); printf(\n); for ( i=0; im; i++ ) printf(%.4lf , hv[i]); printf(\n\n); } return 0;} 裁判输入数据： 0.0 1.0 0.0 1.0 0.0 0.0 5 0.0 0.2 0.5 0.8 1.0 3 0.0 0.5 1.0 100.0 170.0 200.0 30.0 150.0 0.0 5 0.0 0.25 0.5 0.75 1.0 5 0.0 1.0 2.0 3.0 4.0 0.0 60.0 160.0 260.0 300.0 5.0 70.0 100.0 120.0 20.0 10 0.5 1.0 1.5 2.0 2.5 3.0 3.5 3.8 3.95 4.0 标准输出数据： 0.0000 0.9600 1.5000 0.9600 0.0000 100.0000 127.9297 170.0000 195.9766 200.0000 30.0000 165.4688 150.0000 52.9688 0.0000 30.2222 60.0000 105.9303 160.0000 206.3438 260.0000 307.9764 305.7687 299.9796 300.0000 62.6024 70.0000 109.0488 100.0000 92.9745 120.0000 41.2374 -44.8421 -16.2783 20.0000 #includestdio.h #define MAXN 5 /* 最大采样点个数 */ #define MAXM 10 /* 最大估算点个数 */ void Hermite_Interpolation( int N, double t[], doubl