源之:
第2讲 拉格朗日插值多项式
一、拉格朗日插值多项式 为了便于学习,我们首先考虑两种简单情况。 1. 线性插值多项式 设函数f(x)在区间[a,b]上连续,在给定节点x0,x1上的值分别为y0,y1构造一个不超过一次的多项式: L1(x)=a0+a1x 使之满足条件 :L1(x0)=y0, L1(x1)=y1 它的几何意义是,求通过给定点M0(x0,y0), M1(x1,y1)的一直线y=L1(x)近似曲线y=f(x),如图2.1所示 ![]() 图2-1 由解析几何知,过给定两点M0(x0,y0), M1(x1,y1)的直线方程为: ![]() 若记 ![]() 则(2.1)式可写成 ![]() 把 ![]() ![]() ![]() ![]() 线性插值多项式则为节点 ![]() 2. 抛物线插值多项式 设函数 ![]() ![]() ![]() ![]() ![]() 它的几何意义是,求通过给定点 ![]() ![]() ![]() ![]() 图2-2 像前面线性插值那样,设 ![]() | |
共08页 |
下面我们求二次插值基函数 ![]() ![]() 可得 ![]() 表2-1
先求x0点的二次插值基函数 ![]() 零点,于是可把 ![]() ![]() 再利用条件 ![]() ![]() 这样就求得了x0点的二次插值基函数 ![]() 用同样的方法,可以求出x1,x2点的二次插值基函数,它们分别为: ![]() ![]() 于是得求的不超过二次插值多项式为: ![]() | |||||||||
<< 上一页 02 下一页 >> | 共08页 |
(3)n次插值多项式 设函数f(x)在区间[a,b]上连续,在给定n+1个不同节点x0,x1,…,xn上的值分别为y0,y1,…,yn ,要求构造一个不超过n次的多项式: Ln(x)=a0+a1x+…+ anxn , 去近似函数y=f(x),使之满足插值条件 : ![]() 它的几何意义是,用过n+1个不同点M0(x0,y0), M1(x1,y1), M2(x2,y2),…,Mn(xn,yn)的n次代数曲线Ln(x)近似地代替曲线y=f(x),见图2-3 ![]() 图2-3 现在采取上述构造一、二次插值多项式的方法,来构造一般插值多项式Ln(x),设在给定n+1个不同节点x0,x1,…,xn上的n次插值基函数分别为l0(x),l1(x),…,ln(x),则所求不超过n次的插值多项式表示为: ![]() 下面求插值基函数l0(x),l1(x),…,ln(x).由条件:Ln(xi)=yj,i=0,1,…,n. 可得插值基函数l0(x),l1(x),…,ln(x)在各节点上函数值表为: 表2-2
| |||||
<< 上一页 03 下一页 >> | 共08页 |
作为任一xi点的n次插值基函数分别为li(x),由于在除 ![]() 再由条件li(xi)=1得 ![]() ![]() ![]() ![]() (2.3) 其中 ![]() 将(2.3)式称为函数f(x)的拉格朗日插值多项式。 | |
<< 上一页 04 下一页 >> | 共08页 |
二、插值多项式的余项及误差估计式 很显然(2.3)式给出的多项式满足插值条件:Ln(xi)=yj,i=0,1,…,n. 假设另外存在一满足插值条件的不超过n次的多项式Qn(x),设R n(x)= Ln(x)- Qn(x)则R n(x)是不超过n次的多项式,且满足R n(xi)= Ln(xi)- Qn(xi)=yi- yi =0 ( i=0,1,…,n)这说明一个不超过n次的多项式R n(x)具有n+1个不同零点,这只有R n(x)≡0,即 : L(x)= Qn(x) 这就证明了满足插值条件,不超过n次的插值多项式是唯一的。 设Ln(x)是函数y= f(x),不超过n次的拉格朗日插值多项式, 则当x=xi时,有Ln(x)=yj= f(xi),(i=0,1,…,n). 当 x≠xi时,Ln(x) 与f(x)一般是不相等的,其差函数 R n(x)= f(x)- Ln(x)称为拉格朗日插值余项(或拉格朗日插值误差)。 定理2.1 函数f(x) 在区间[a,b]上存在n+1阶导数,则有拉格朗日插值余项为 ![]() 证明 由 ![]() 即x0,x1,…,xn皆为Rn(x)的零点知Rn(x)含有因子ω(x),设Rn(x)=K(x) ω(x) 其中K(x)为待定函数.为求出K(x),我们引进辅助函数g(t)=f(t)- Ln(t)- K(x) ω(t) 则g(t)在[a,b]上至少有n+2个互异零点x0,x1,…,xn,x,且存在n+1阶导数.在n+2个零点构成的n+1个子区间上,应用微分学罗尔定理,则每个子区间中至少有g/(t)的一个零点 ,从而在[a,b]内至少有g/(t)的n+1个零点.在这n+1个零点构成的n个子区间上,再应用罗尔定理得在[a,b]内至少有g//(t)的n个零点,…在[a,b]内至少有g(n+1)(t)的一个零点,设为ξ,则有 ![]() 于是可得 ![]() ![]() | |
<< 上一页 05 下一页 >> | 共08页 |
若记 ![]() 则有 ![]() (2.5)式称为拉格朗日插值误差估计式。 对于线性插值,其余项为 ![]() 对于二次插值,其余项为 ![]() 例1 已知 ![]() ![]() 解 取最接近x=5的两点x0=4,x1=9为插值节点,两个插值基函数分别为 ![]() ![]() 故有 ![]() ![]() 下面估计误差 ![]() 因为 ![]() 所以 ![]() | |
<< 上一页 06 下一页 >> | 共08页 |
例2 已知的函数表
求在[0,2]内的零点近似值。 解 因为yi关于x严格单调减少,用反插值法求f(x) 零点的近似值比较简单, 具体作法如下: 先作反函数表
将节点x0=8,x1=-7.5,x2=-18及对应函数值y0=0,y1=1,y2=2代入二次拉格朗日插值多项式(2.2),再令x=0,得 ![]() 于是得f(x)在[0,2]内零点 ![]() 值得注意的是,只有所给函数(或函数表)在[a,b]上严格单调情况下,才能使用反插值方法,否则可能得出错误结果。 | |||||||||||||||||
<< 上一页 07 下一页 >> | 共08页 |
练习 2.1 1.已知函数表
求抛物插值多项式 ![]() 2.已知插值基函数 ![]() 证明:当 ![]() ![]() 习题2 3. 已知连续函数 ![]()
求 ![]() ![]() | |||||||||||||
<< 上一页 08 首 页 >> | 共08页 |