分段三次Hermite插值

摘要用函数来表示变量间的数量关系广泛应用于各学科领域，但是在实际问题中，往往是通过实验、观测以及计算等方法，得到的是函数在一些点上的函数值。如何通过这些离散数据找到函数的一个满足精度要求且便于使用的近似表达式，是经常遇到的问题。对于这类问题我们解决的方法为插值法，而最常用也最简单的插值方法就是多项式插值。当然用插值法得到的近似表达式必须满足插值条件即假设给定了n+1个点的自变量的值以及函数值，近似函数必须要过这n+1个点。多项式插值，从几何角度看，就是寻求n次代数曲线y=Pn（x）通过n+1个点作为f（x）的近似。但是随着插值节点个数的增加，高次插值多项式的近似效果并不理想。根据大量实验得出，在进行高次多项式插值时，会出现龙格现象。因此，为了解决这样的一个问题，我们可以通过缩小插值区间的办法达到减小误差的目的。但是当在每个小区间上用一次函数进行插值时，有很好的收敛性但是光滑度不够，因此本实验将用三次Hermite进行插值，做具体的讨论和学习。关键词：龙格现象 分段差值 三次Hermite进行插值1、实验目的1) 通过对分段三次Hermite插值算法程序的编写，提高自己编写程序的能力2) 体会分段三次Hermite插值比分段线性插值优越在哪里3) 用实验报告的形式展现，提高自己在写论文方面的能力2、算法流程分段线性插值多项式S(x)在插值区间a,b上只能保证连续性，而不光滑。要想得到在插值区间上光滑的分段线性插值多项式，可采用分段埃尔米特（Hermite）插值，这里我们考虑在整个a,b上用分段三次埃尔米特插值多项式来逼近f(x)。一般的将带有导数的插值多项式称为Hermite插值多项式。如果已知函数y=f(x)在节点a=x0x1xn=b处的函数的值和导数值：yi=fxi,yi=fxi,i=0,1,2,n则在小区间xi-1,xi上有四个插值条件：yi-1=fxi-1,yi=fxiyi-1=fxi-1,yi=fxi故能构造一个三次多项式Hi(x)，并称为三次Hermite插值多项式。这时在整个a,b上可以用分段三次Hermite插值多项式来逼近f(x)。Hx=H1x,xx0,x1 H2x,xx1,x2 Hnx,xxn-1,xn其中Hix,xxi-1,xi满足条件：Hixi-1=fxi-1=yi-1,Hixi=fxi=yiHixi-1=fxi-1=yi-1,Hixi=fxi=yi关于Hix的构造，我们可以通过基函数来进行，这时令Hix=yi-1i-1x+yiix+yi-1i-1x+yiix其中i-1x、ix、i-1x和ix均为三次多项式，并称为三次Hermite插值多项式的基函数。对上式两边关于x求导，得到Hix=yi-1i-1x+yiix+yi-1i-1x+yiix则由插值条件可以分别给出基函数满足的条件：i-1xi-1=1,ixi-1=0,i-1xi-1=0,ixi-1=0i-1xi=0, ixi=1,i-1xi=0,ixi=0i-1xi-1=0,ixi-1=0,i-1xi-1=1,ixi-1=0i-1xi=0,ixi=0,i-1xi=1,ixi=1下面具体求解基函数i-1x、ix、i-1x和ix。由上面的条件的第一列可以得到i-1x满足条件： i-1xi-1=1, i-1xi=0,i-1xi-1=0,i-1xi=0 （1）由上式中的第二、第四个条件可知i-1x应该具有形式 i-1x=x-xi2(ax+b) （2）这时 i-1x=2x-xiax+b+a(x-xi)2 （3）再由（1）式中的第一、第三个条件分别带入（2）式（3）式得到hi2axi-1+b=1 -2hiaxi-1+b+ahi2=0解此线性方程组得到a=2hi3,b=1hi2-2xi-1hi3将a、b代入（2）式得到i-1x=x-xi22hi3x+1hi2-2xi-1hi3=(1+2x-xi-1hi)x-xi2hi2类似地有ix=(1-2x-xi-1hi)x-xi2hi2i-1x=1hi2(x-xi-1) x-xi2ix=1hi2(x-xi) x-xi-12因此将得到Hix=hi+2(x-xi-1)x-xi2hi3yi-1+hi-2(x-xi)x-xi-12hi3yi+(x-xi-1)x-xi2hi2yi-1+(x-xi)x-xi-12hi2yi这样，便求出了分段三次Hermite插值多项式：Hx=H1x,xx0,x1 H2x,xx1,x2 Hnx,xxn-1,xn3、数值算例已知下列的条件xi1 2yi2 3yi1 -1通过分段三次Hermite插值法，求解当x=1.5时的y值。解：具体的程序如下所示：#include stdafx.hfloat Hermite(float x,float y,float z,float x1,int len)int i=0;float s=0;float h=0;float L1=0;float L2=0;float L3=0;float L4=0;for(i=0;i=xi & x1xi+1)break;i=i+1;h=xi-xi-1;L1=(h+2*(x1-xi-1)*(x1-xi)*(x1-xi)/(h*h*h);L2=(h+2*(x1-xi)*(x1-xi-1)*(x1-xi-1)/(h*h*h);L3=(x1-xi-1)*(x1-xi)*(x1-xi)/(h*h);L4=(x1-xi-1)*(x1-xi-1)*(x1-xi)/(h*h);s=L1*yi-1+L2*yi+L3*zi-1+L4*zi;return s;float Hermite(float x,float y,float x1,int len);void main()float x=1,2;float y=2,3;float z=1,-1;int len=sizeof(x)/sizeof(x0);float x1=0;float s=0;printf(请输入要求解的x1的值:n);scanf(%f,&x1);s=Hermite(x,y,z,x1,len);printf(经过分段三次Hermite插值的结果为:n);pr