分段三次Hermite 插值

《分段三次Hermite 插值》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、摘要用函数来表示变量间的数量关系广泛应用于各学科领域，但是在实际问题中，往往是通过实验、观测以及计算等方法，得到的是函数在一些点上的函数值。如何通过这些离散数据找到函数的一个满足精度要求且便于使用的近似表达式，是经常遇到的问题。对于这类问题我们解决的方法为插值法，而最常用也最简单的插值方法就是多项式插值。当然用插值法得到的近似表达式必须满足插值条件即假设给定了n+1个点的自变量的值以及函数值，近似函数必须要过这n+1个点。多项式插值，从几何角度看，就是寻求n次代数曲线y=Pn（x）通过n+1个点作为f（x）的近似。但是随着插值节

2、点个数的增加，高次插值多项式的近似效果并不理想。根据大量实验得出，在进行高次多项式插值时，会出现龙格现象。因此，为了解决这样的一个问题，我们可以通过缩小插值区间的办法达到减小误差的目的。但是当在每个小区间上用一次函数进行插值时，有很好的收敛性但是光滑度不够，因此本实验将用三次Hermite进行插值，做具体的讨论和学习。关键词：龙格现象分段差值三次Hermite进行插值51、实验目的1)通过对分段三次Hermite插值算法程序的编写，提高自己编写程序的能力2)体会分段三次Hermite插值比分段线性插值优越在哪里3)用实验报告的形

3、式展现，提高自己在写论文方面的能力2、算法流程分段线性插值多项式S(x)在插值区间a,b上只能保证连续性，而不光滑。要想得到在插值区间上光滑的分段线性插值多项式，可采用分段埃尔米特（Hermite）插值，这里我们考虑在整个a,b上用分段三次埃尔米特插值多项式来逼近f(x)。一般的将带有导数的插值多项式称为Hermite插值多项式。如果已知函数y=f(x)在节点a=x0

4、=fxiy'i-1=f'xi-1,yi'=f'xi故能构造一个三次多项式Hi(x)，并称为三次Hermite插值多项式。这时在整个a,b上可以用分段三次Hermite插值多项式来逼近f(x)。Hx=H1x,x∈x0,x1H2x,x∈x1,x2⋯⋯⋯⋯⋯⋯⋯Hnx,x∈xn-1,xn其中Hix,x∈xi-1,xi满足条件：Hixi-1=fxi-1=yi-1,Hixi=fxi=yiH'ixi-1=f'xi-1=y'i-1,H'ixi=f'xi=y'i关于Hix的构造，我们可以通过基函数来进行，这时令Hix=yi-1φi-1x+yiφ

5、ix+y'i-1ψi-1x+y'iψix其中φi-1x、φix、ψi-1x和ψix均为三次多项式，并称为三次Hermite插值多项式的基函数。对上式两边关于x求导，得到5H'ix=yi-1φ'i-1x+yiφ'ix+y'i-1ψ'i-1x+y'iψ'ix则由插值条件可以分别给出基函数满足的条件：φi-1xi-1=1,φixi-1=0,ψi-1xi-1=0,ψixi-1=0φi-1xi=0,φixi=1,ψi-1xi=0,ψixi=0φ'i-1xi-1=0,φ'ixi-1=0,ψ'i-1xi-1=1,ψ'ixi-1=0φ'i-1x

6、i=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解此线性方

7、程组得到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-xi2hi2ψi-1x=1hi2(x-xi-1)x-xi2ψix=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-xi2hi2y'i-1+(x-xi)x-xi-12hi2y'i这样

8、，便求出了分段三次Hermite插值多项式：5Hx=H1x,x∈x0,x1H2x,x∈x1,x2⋯⋯⋯⋯⋯⋯Hnx,x∈xn-1,xn3、数值算例已知下列的条件xi12yi23y'i1-1通过分段三次Hermite插值法，求解当x=1.5时的y值。解：具体的程序