我们知道,在许多实际问题中,函数y=f(x)在某区间[a,b]上存在,但函数关系往往很复杂,甚至没有明显的解析表达式,例如,通过实验或观测得到一系列数据,这些数据是与自变量的某些点xi相应的函数值yi,而我们需计算未观测到的点的函数值。根据已观测数据来构造一个既能反映函数的特征又便于计算的而且较为简单的函数p(x)来替代f(x),此即插值问题。
设函数y=f(x)定义在区间[a,b]上,x0,x1,…,xn是[a,b]上的n+1个互异点,且仅仅在这些点处函数值yi=f(xi)为已知,要构造一个函数p(x),使得
p(xi)=yi(i=0,1,…,n)
并要求误差
R(x)=f(x)-p(x)
的绝对值|R(x)|在区间[a,b]上任意一点或整个区间[a,b]上比较小,即说明p(x)较好地逼近了f(x)。点x0,x1,…,xn称为插值基点或简称为基点。p(x)称为插值函数,f(x)为被插函数,求f(x)的插值函数的方法称为插值法,R(x)称为插值公式的余项。插值函数p(x)在n+1个插值基点xi(i=0,1,…,n)处的值与f(x)相等,在其他点就用p(x)的值作为f(x)的近似值。
我们用p(x)作为f(x)的插值函数,除要求p(x)在某些意义上更好地逼近f(x)外,还希望f(x)是较简单的函数,或者它便于计算机计算,因此,我们常用多项式、有理分式和三角多项式作为插值函数。选择不同的函数类作为插值函数逼近f(x),其效果是不同的,需根据实际问题中被插函数f(x)的特性,来选择合适的插值函数,例如,可选取多项式作为插值函数,此即所谓的插值多项式。若p(x)是次数不超过n的多项式,用Pn(x)表示,即
Pn(x)=α0+α1x+α2x2+…+αnxn
则称Pn(x)为n次插值多项式。可以证明,在所有次数不超过n的多项式中,满足插值函数要求的Pn(x)存在且唯一。
我们考虑最简单的插值问题,设函数在区间[a,b]上n+1个互异基点x0,x1,…,xn上的函数值为
求插值多项式li(x),满足条件
li(xj)=δij(i,j=0,1,…,n)
由上式知,x0,x1,…,xi-1,xi+1,…,xn是li(x)的根,于是可令
li(x)=Ai(x-x0)(x-x1)…(x-xi-1)(x-xi+1)…(x-xn)
再由li(xi)=1,得
于是
这样就得到n+1个n次多项式l0(x),l1(x),…,ln(x),称为插值基函数。
当n=1时的一次基函数为
当n=2时的二次基函数为
一般地,设函数在区间[a,b]上n+1个互异基点x0,x1,…,xn上的函数值分别为y0,y1,…,yn,若n次插值多项式In(x),满足条件
Ln(xj)=yj(j=0,1,…,n)
则有拉格朗日插值多项式
其中,l0(x),l1(x),…,ln(x)为n次插值基函数。
拉格朗日插值多项式的插值条件只要求在插值基点上插值函数与被插值函数的函数值相等,即Ln(xi)=yi,有时不仅要求插值多项式在插值基点上与被插值函数的函数值相等,还要求插值多项式的导数在这些点上与被插函数的导数值相等,即要求满足插值条件
的次数不超过2n+1的插值多项式H2n+1(x),这就是埃尔米特(Hermite)插值。下面对埃尔米特插值作简要讨论。
设在结点a≤x0<x1<…<xn≤b上,yj=f(xj),mj=f′(xj)(j=0,1,…,n)
要求插值多项式H(x)满足条件
这里给出了2n+2个条件,可唯一确定一个次数不超过2n+1的多项式
H2n+1(x)=H(x)
其形式为
H2n+1(x)=a0+a1x+…+a2n+1x2n+1
显然,直接根据H(x)满足的条件确定2n+2个系数a0,a1,…,a2n+1将非常复杂,为此,我们仍采用拉格朗日插值多项式的基函数方法,先求插值基函数αj(x)和βj(x)(j=0,1,…,n),共有2n+2个,每一个基函数都是2n+1次多项式,且满足条件
(www.daowen.com)
βj(xk)=0,βj′(xk)=δjk(j,k=0,1,…,n)
于是,插值多项式H(x)=H2n+1(x)写成用插值基函数表示的形式
根据限定条件,显然有
H2n+1(xk)=yk,H2′n+1(xk)=mk(k=0,1,…,n)现在的问题就是求满足条件的基函数αj(x)和βj(x),为此,可利用拉格朗日插值函数lj(x)。
下面我们考虑只有两个结点的三次埃尔米特插值。设插值点为(x0,y0),(x1,y1),要求次数不超过3的多项式H3(x),满足下列条件
式中,mi=f′(xi)(i=0,1)。
类似于拉格朗日插值多项式的构造过程,仍采用基函数的方法来构造H3(x),令
H3(x)=y0α0(x)+y1α1(x)+m0β0(x)+m1β1(x)
式中,α0(x),α1(x),β0(x),β1(x)是插值基函数,为了满足插值条件,它们应满足下列条件(见表3.4.1),且α0(x),α1(x),β0(x),β1(x)均为次数不超过3的多项式。
表3.4.1 插值基函数
由表3.4.1可知,α0(x1)=α0′(x1)=0,故α0(x)应含有(x-x1)2因子,又α0(x)是次数不过3的多项式,因而可将它写成
α0(x)=[a+b(x-x0)](x-x1)2式中,a,b为待定常数。由α0(x0)=0,得
由α′0(x0)=0,得
于是
类似地,可得到α1(x)为
对于β0(x),注意到β0(x0)=β0(x1)=β0′(x1)=0,于是β0(x)应含有(x-x0)(x-x1)2因子,又β0(x)是次数不过3的多项式,因而可将它写成
β0(x)=c(x-x0)(x-x1)2
式中,c为待定常数。
由β0′(x0)=1,得
于是
类似地,得到
显然α0(x),α1(x),β0(x),β1(x)可简单地表示为
α0(x)=[1+2l1(x)]l20(x)
α1(x)=[1+2l0(x)]l21(x)
β0(x)=(x-x0)l02(x)
β1(x)=(x-x1)l21(x)
其中,l0(x),l1(x)为以(x0,y0),(x0,y0)为插值点的拉格朗日一次基函数。
与拉格朗日插值类似,可以证明,三次埃尔米特插值多项式存在且唯一。
免责声明:以上内容源自网络,版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。