计算方法第5章_插值法.pptVIP

* 某些插值问题不仅给出了插值节点的坐标位置，还给出了某些插值节点处的若干阶导数值。它不仅要求插值曲线经过插值节点，还要求插值曲线在插值节点处的各阶导数值与给出的各阶导数值相等。满足这种要求的多项式插值称为埃尔米特（Hermite）插值。 5.5 埃尔米特插值 一、埃尔米特插值简介 定义：若存在不超过2n＋1次的多项式函数H(x)，使H(x)与被插函数f(x)在插值节点x=x0,x1,……,xn（设x0＜x1＜……＜xn）处的函数值和1阶导数值相等，即H(xi)=f(xi)且H(xi)=f(xi) ，其中i=0,1,2,……,n，则称H(x)为带1阶导数的埃尔米特插值函数，称这种插值为带1阶导数的埃尔米特插值。 带1阶导数的埃尔米特插值函数是唯一的。 由埃尔米特插值的定义，可以由待定系数法求解埃尔米特插值函数。 * 2点3次埃尔米特插值是一种使用广泛的插值方法。仅有2个插值节点的带1阶导数的埃尔米特插值函数是1个3次代数多项式函数，因此又称为2点3次埃尔米特插值。 5.5 埃尔米特插值 二、2点3次埃尔米特插值 设插值节点在x=x0,x1处，对被插函数f(x)的2点3次埃尔米特插值函数的一般形式为： H(x)=h0(x)f(x0)+h1(x)f(x1)+?0(x)f(x0)+?1(x)f(x1) = f(x0)+ f(x1) ＋(x-x0) f(x0)＋(x-x0) f(x1) * 上述2点3次埃尔米特插值可以推广到n＋1个插值节点的带1阶导数的埃尔米特插值。 5.5 埃尔米特插值 三、带1阶导数的埃尔米特插值 设插值节点在x=x0,x1,……,xn,处，对被插函数f(x)的带1阶导数埃尔米特插值函数的一般形式为： H(x)= + 其中辅助函数hi(x)=(1-2li(xi)(x-xi))li2(x)，?i(x)=(x-xi)li2(x) li(x)为拉格郎日插值基函数，即li(x)= 式中i=0,1,2,……,n。 * 5.5 埃尔米特插值 n个插值节点带1阶导数的埃尔米特插值的算法 输入插值节点的个数n。 输入插值节点的坐标x[n],f[n]及1阶导数df[n]，插值点_x。 _y=0; for(i=0;i=n-1;i++) l=1; ① 计算li(x) for(j=0;j=n-1;j++) j≠i Y N l*=(_x-x[j])/(x[i]-x[j]); dl=0; ② 计算li’(xi) for(j=0;j=n-1;j++) j≠i Y N dl+=1/(x[i]-x[j]); _y+=(1-2*dl*(_x-x[i]))*l*l*f[i]+(_x-x[i])*l*l*df[i]; 输出插值点(_x,_y)。 * 5.5 埃尔米特插值 n个插值节点带1阶导数的埃尔米特插值的程序（1/2） #include stdio.h #define MAXSIZE 50 void input(double x[],double f[],double df[],long n); void main(void) { double x[MAXSIZE],f[MAXSIZE],df[MAXSIZE],_x,_y,l,dl;