理论教育 插值函数概述及应用-偏微分方程数值解法:土建类

插值函数概述及应用-偏微分方程数值解法:土建类

时间:2023-10-26 理论教育 版权反馈
【摘要】:,xn称为插值基点或简称为基点。可以证明,在所有次数不超过n的多项式中,满足插值函数要求的Pn存在且唯一。,ln,称为插值基函数。,n)于是,插值多项式H=H2n+1写成用插值基函数表示的形式根据限定条件,显然有H2n+1=yk,H2′n+1=mk(k=0,1,…表3.4.1 插值基函数由表3.4.1可知,α0=α0′=0,故α0应含有2因子,又α0是次数不过3的多项式,因而可将它写成α0=[a+b]2式中,a,b为待定常数。

插值函数概述及应用-偏微分方程数值解法:土建类

我们知道,在许多实际问题中,函数y=fx)在某区间[ab]上存在,但函数关系往往很复杂,甚至没有明显的解析表达式,例如,通过实验或观测得到一系列数据,这些数据是与自变量的某些点xi相应的函数值yi,而我们需计算未观测到的点的函数值。根据已观测数据来构造一个既能反映函数的特征又便于计算的而且较为简单的函数px)来替代fx),此即插值问题。

设函数y=fx)定义在区间[ab]上,x0x1,…,xn是[ab]上的n+1个互异点,且仅仅在这些点处函数值yi=fxi)为已知,要构造一个函数px),使得

pxi=yii=0,1,…,n

并要求误差

Rx=fx)-px

的绝对值|Rx)|在区间[ab]上任意一点或整个区间[ab]上比较小,即说明px)较好地逼近了fx)。点x0x1,…,xn称为插值基点或简称为基点。px)称为插值函数,fx)为被插函数,求fx)的插值函数的方法称为插值法Rx)称为插值公式的余项。插值函数px)在n+1个插值基点xii=0,1,…,n)处的值与fx)相等,在其他点就用px)的值作为fx)的近似值。

我们用px)作为fx)的插值函数,除要求px)在某些意义上更好地逼近fx)外,还希望fx)是较简单的函数,或者它便于计算机计算,因此,我们常用多项式、有理分式和三角多项式作为插值函数。选择不同的函数类作为插值函数逼近fx),其效果是不同的,需根据实际问题中被插函数fx)的特性,来选择合适的插值函数,例如,可选取多项式作为插值函数,此即所谓的插值多项式。若px)是次数不超过n的多项式,用Pnx)表示,即

Pnx01x+α2x2+nxn

则称Pnx)为n次插值多项式。可以证明,在所有次数不超过n的多项式中,满足插值函数要求的Pnx)存在且唯一。

我们考虑最简单的插值问题,设函数在区间[ab]上n+1个互异基点x0x1,…,xn上的函数值为

978-7-111-44528-9-Chapter03-95.jpg

求插值多项式lix),满足条件

lixjijij=0,1,…,n

由上式知,x0x1,…,xi-1xi+1,…,xnlix)的根,于是可令

lix=Aix-x0)(x-x1)…(x-xi-1)(x-xi+1)…(x-xn

再由lixi)=1,得

978-7-111-44528-9-Chapter03-96.jpg

于是

978-7-111-44528-9-Chapter03-97.jpg

这样就得到n+1个n次多项式l0x),l1x),…,lnx),称为插值基函数。

n=1时的一次基函数为

978-7-111-44528-9-Chapter03-98.jpg

n=2时的二次基函数为

978-7-111-44528-9-Chapter03-99.jpg

一般地,设函数在区间[ab]上n+1个互异基点x0x1,…,xn上的函数值分别为y0y1,…,yn,若n次插值多项式Inx),满足条件

Lnxj=yjj=0,1,…,n

则有拉格朗日插值多项式

978-7-111-44528-9-Chapter03-100.jpg

其中,l0x),l1x),…,lnx)为n次插值基函数。

拉格朗日插值多项式的插值条件只要求在插值基点上插值函数与被插值函数的函数值相等,即Lnxi)=yi,有时不仅要求插值多项式在插值基点上与被插值函数的函数值相等,还要求插值多项式的导数在这些点上与被插函数的导数值相等,即要求满足插值条件

978-7-111-44528-9-Chapter03-101.jpg

的次数不超过2n+1的插值多项式H2n+1x),这就是埃尔米特(Hermite)插值。下面对埃尔米特插值作简要讨论。

设在结点ax0x1<…<xnb上,yj=fxj),mj=f′xj)(j=0,1,…,n

要求插值多项式Hx)满足条件

978-7-111-44528-9-Chapter03-102.jpg

这里给出了2n+2个条件,可唯一确定一个次数不超过2n+1的多项式

H2n+1x=Hx

其形式为

H2n+1x=a0+a1x++a2n+1x2n+1

显然,直接根据Hx)满足的条件确定2n+2个系数a0a1,…,a2n+1将非常复杂,为此,我们仍采用拉格朗日插值多项式的基函数方法,先求插值基函数αjx)和βjx)(j=0,1,…,n),共有2n+2个,每一个基函数都是2n+1次多项式,且满足条件

978-7-111-44528-9-Chapter03-103.jpg(www.daowen.com)

βjxk=0,βjxkjkjk=0,1,…,n

于是,插值多项式Hx)=H2n+1x)写成用插值基函数表示的形式

978-7-111-44528-9-Chapter03-104.jpg

根据限定条件,显然有

H2n+1xk=ykH2n+1xk=mkk=0,1,…,n)现在的问题就是求满足条件的基函数αjx)和βjx),为此,可利用拉格朗日插值函数ljx)。

下面我们考虑只有两个结点的三次埃尔米特插值。设插值点为(x0y0),(x1y1),要求次数不超过3的多项式H3x),满足下列条件

978-7-111-44528-9-Chapter03-105.jpg

式中,mi=f′xi)(i=0,1)。

类似于拉格朗日插值多项式的构造过程,仍采用基函数的方法来构造H3x),令

H3x=y0α0x+y1α1x+m0β0x+m1β1x

式中,α0x),α1x),β0x),β1x)是插值基函数,为了满足插值条件,它们应满足下列条件(见表3.4.1),且α0x),α1x),β0x),β1x)均为次数不超过3的多项式。

3.4.1 插值基函数

978-7-111-44528-9-Chapter03-106.jpg

由表3.4.1可知,α0x1)=α0x1)=0,故α0x)应含有(x-x1)2因子,又α0x)是次数不过3的多项式,因而可将它写成

α0x=[a+bx-x0)](x-x1)2式中,ab为待定常数。由α0x0)=0,得

978-7-111-44528-9-Chapter03-107.jpg

α′0x0)=0,得

978-7-111-44528-9-Chapter03-108.jpg

于是

978-7-111-44528-9-Chapter03-109.jpg

类似地,可得到α1x)为

978-7-111-44528-9-Chapter03-110.jpg

对于β0x),注意到β0x0)=β0x1)=β0x1)=0,于是β0x)应含有(x-x0)(x-x1)2因子,又β0x)是次数不过3的多项式,因而可将它写成

β0x=cx-x0)(x-x1)2

式中,c为待定常数。

β0x0)=1,得

978-7-111-44528-9-Chapter03-111.jpg

于是

978-7-111-44528-9-Chapter03-112.jpg

类似地,得到

978-7-111-44528-9-Chapter03-113.jpg

显然α0x),α1x),β0x),β1x)可简单地表示为

α0x=[1+2l1x)]l20x

α1x=[1+2l0x)]l21x

β0x=x-x0l02x

β1x=x-x1l21x

其中,l0x),l1x)为以(x0y0),(x0y0)为插值点的拉格朗日一次基函数。

978-7-111-44528-9-Chapter03-114.jpg

与拉格朗日插值类似,可以证明,三次埃尔米特插值多项式存在且唯一。

免责声明:以上内容源自网络,版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。

我要反馈