资源描述：
《为三次样条插值函数.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、2.5埃尔米特插值有些实际的插值问题不但要求在节点上函数值相等，下面只讨论函数值与导数值个数相等的情况.满足这种要求的插值多项式就是埃尔米特插值多项式.而且还要求对应的导数值也相等，甚至要求高阶导数也相等.1(5.1)这里共有个插值条件，可唯一确定一个次数不超过的多项式，问题是求插值多项式，设在节点上，现在仍采用求拉格朗日插值多项式的基函数方法.满足条件其形式为2将满足条件（5.1）的插值多项式写成用插值基函数表示的形式(5.3)先求出个插值基函数及，每一个基函数都是次多项式，且满足条件(5.2)3令由条件（5.2），有由插值基函数所满足的条件
2、（5.2），有下面的问题就是如何求出这些基函数及利用拉格朗日插值基函数4解出由于整理得5于是(5.4)两端取对数再求导，得同理，可得(5.5)6可以证明满足条件（5.1）的插值多项式是惟一的.用反证法，假设及均满足条件（5.1），这样，有重根，但是不高于次的多项式，于是在每个节点上的值及导数值均为零，即为二重根.故惟一性成立.7其中且与有关.若在内的阶导数存在，则其插值余项(5.6)仿照拉格朗日插值余项的证明方法，可以证明：8插值多项式（5.3）的重要特例是的情形.这时可取节点为及，插值多项式为，(5.7)相应的插值基函数为它们满足条件满足91
3、0(5.8)(5.9)根据及的一般表达式（5.4）及（5.5），可得到11(5.10)其余项，于是满足条件（5.7）的插值多项式是由（5.6）得12求满足及由给定的4个条件，可确定次数不超过3的插值多项式.由于此多项式通过点的插值多项式及其余项表达式.例4故其形式为13待定常数，可由条件确定，其中为待定函数.为了求出余项的表达式，通过计算可得可设14显然故在内有5个零点（二重根算两个）.反复应用罗尔定理，得在内至少有一个零点ξ，构造且故有15(5.11)式中位于和所界定的范围内.余项表达式为于是162.6分段低次插值2.6.1高次插值的病态性质
4、这是因为对任意的插值节点，当时，不一定收敛到.在次数增加时逼近的精度不一定也增加.根据区间上给出的节点做出的插值多项式17所构造的拉格朗日插值多项式为以上的个等距节点考虑函数，它在上的各阶导数均存在.令则18表2-5列出了时的的计算结果及在上的误差19可见，随的增加，的绝对值几乎成倍增加.这说明当时在上是不收敛的.Runge证明了，存在一个常数，使得当时，而当时发散.20取根据计算画出及在上的图形，见图2-5.图2-521从图上看到，在附近,与偏离很远，这说明用高次插值多项式近似效果并不好.通常不用高次插值，而用分段低次插值.22下图是用Mat
5、lab完成的Lagrange插值（附程序）：23附：Lagrange插值程序n=11;m=61;x=-5:10/(m-1):5;y=1./(1+x.^2);z=0*x;x0=-5:10/(n-1):5;y0=1./(1+x0.^2);y1=lagr1(x0,y0,x);plot(x,z,’r’,x,y,’k:’,x,y1,’r’)gtext(‘Lagr.’),gtext(‘y=1/(1+x^2)’)title(‘Lagrange’)24附：Lagrange插值子程序lagr1:functiony=lagr1(x0,y0,x)n=length(x
6、0);m=length(x);fori=1:mz=x(i);s=0.0;fork=1:np=1.0;forj=1:nifj~=kp=p*(z-x0(j))/(x0(k)-x0(j));endends=p*y0(k)+s;endy(i)=s;end252.6.2分段线性插值所谓分段线性插值就是通过插值点用折线段连接起来逼近由于升高插值多项式的阶数有时并不能达到提高精度的效果,所以实际中往往采用分段插值的思想.分段插值的基本思想是将插值区间划分为若干个小区间,然后在每个小区间上做满足一定条件的低阶插值.26设已知节点上的函数值记求一折线函数,满足：
7、在每个小区间上是线性函数.则称为分段线性插值函数.27由定义可知在每个小区间上可表示为(6.1)若用插值基函数表示，则在整个区间上为(6.2)其中基函数满足条件其形式是28(6.3)利用插值余项（2.17）得到分段线性插值的误差估计29或写成(6.4)其中30当时，故另一方面，这时这种性质称为局部非零性质.分段线性插值基函数只在附近不为零，在其他地方均为零，利用的局部非零性质及知，31现在证明,这里是函数在区间上的连续模，即对任意两点，只要就有考虑32称为在上的连续模，当时，就有由前式可知，当时有因此，只要，就有在上一致成立，故在上一致收敛到.
8、33下图是用Matlab完成的分段线性插值（附程序）：34附：分段线性插值程序n=11;m=61;x=-5:10/(m-1):5;y=1./(1+x.
