购买
下载掌阅APP,畅读海量书库
立即打开
畅读海量书库
扫码下载掌阅APP

2.2 拉格朗日插值

2.2.1 线性插值

对给定的插值节点,为求得插值多项式p n (x)可以有各种不同的方法,下面先讨论n=1的简单情形。

已知

求作一次多项式L 1 (x)=a 0 +a 1 x,使它满足条件:L 1 (x 0 )=y 0 ,L 1 (x 1 )=y 1 ,这种插值称为线性插值。

将已知条件代入,得线性方程组

解方程组得

所以

为了便于推广,记

则函数l 0 (x)和l 1 (x)具有性质

函数l 0 (x)和l 1 (x)称为线性插值基函数。于是可将求得的一次插值多项式改写为

注: 用线性插值进行近似计算,当插值区间较小时,近似程度较高。

例2.1 已知f(x)的函数表

的近似值。

利用线性插值公式(2.9)得

于是

2.2.2 抛物线插值

下面讨论n=2的情况。

已知

设二次多项式L 2 (x)=a 0 +a 1 x+a 2 x 2 满足L(x i )=y i (i=0,1,2),则根据定理2.1可求出a 0 ,a 1, a 2 ,得到相应的唯一的二次插值多项式,但计算量较大。下面仿照线性插值,用插值基函数构造抛物线插值多项式。

设L 2 (x)=l 0 (x)y 0 +l 1 (x)y 1 +l 2 (x)y 2 ,其中

下面以l 0 (x)为例说明基函数的求法。由于l 0 (x)为二次多项式,由式(2.10)可设

l 0 (x)=A(x-x 1 )(x-x 2

并由l 0 (x 0 )=1求出 所以

同理

于是

2.2.3 拉格朗日插值

(1)拉格朗日插值基函数

由线性插值式(2.9),抛物线插值式(2.11)可知,插值多项式可看作某些函数的线性组合,组合系数为插值节点的函数值y i

定义2.1 称n次多项式l 0 (x),l 1 (x),l 2 (x),…,l n (x)为在n+1个节点x i (i=0,1,…,n)上的n次插值基函数,其中

根据定义,插值基函数的性质:

②l k (x)(k=0,1,…,n)为由插值节点x 0 ,x 1 ,x 2 ,…,x n 唯一确定的n次多项式。

③基函数和每一个节点有关,与被插函数无关。

(2)拉格朗日插值多项式

定理2.3 对于给定的n+1个插值节点x 0 ,x 1 ,…,x n

为插值问题的n次多项式插值函数,其中

称式(2.13)为多项式插值的拉格朗日形式,简称为拉格朗日插值多项式。

因为l k (x)为n次多项式,所以L n (x)为次数不超过n的代数多项式。

取x=x j ,由于 ,从而 ,满足插值条件。于是 为满足已知n+1个插值节点函数值插值问题的n次多项式插值函数,即n次拉格朗日插值多项式。

注: 线性插值和抛物线插值分别是1次和2次拉格朗日插值多项式。

例2.2 设x 0 ,x 1 ,x 2 ,…,x n 为插值节点,l k (x)(k=0,1,…,n)(n≥1)为拉格朗日插值基函数,则

取f(x)≡1,y i ≡1,(i=0,1,2,…,n)

则f(x)的n次拉格朗日插值为

插值余项R n (x)=f(x)-L n (x),而余项

例2.3 已知y=lnx的函数表

分别用拉格朗日线性和抛物线插值求ln11.5的近似值,并估计误差。

线性插值,其两节点x 0 =11,x 1 =12,插值基函数为

l 0 (x)=-(x-12),l 1 (x)=x-11

所以

L 1 (x)=-2.3979(x-12)+2.4849(x-11)

ln11.5≈L 1 (11.5)=2.4414

余项

由于

抛物线插值,节点x 0 =11,x 1 =12,x 2 =13,则可得

同理可得

例2.4 设f(x)=x 4 ,试用拉格朗日余项定理写出以-1,0,1,2为插值节点的三次插值多项式。

设所求多项式为L 3 (x),由余项定理得

于是

L 3 (x)=f(x)-(x+1)(x-0)(x-1)(x-2)=2x 3 +x 2 -2x lz67vnJWV7fTZofLN74i+b8K/bE4foC4p0XPBd6UqKBOpFYBAMTY3EWAB1P7ElMc

点击中间区域
呼出菜单
上一章
目录
下一章
×