插值余项
若在[ a , b a,b a,b]上用 L n ( x ) L_n(x) Ln(x)近似 f ( x ) f(x) f(x),则其截断误差为 R n ( x ) = f ( x ) − L n ( x ) R_n(x)=f(x)-L_n(x) Rn(x)=f(x)−Ln(x), R n ( x ) R_n(x) Rn(x)也称为插值多项式的余项或插值余项.
定理
关于插值余项估计有以下定理:
设 f ( n ) ( x ) f^{(n)}(x) f(n)(x)在 [ a , b ] [a,b] [a,b]上连续 , f ( n + 1 ) ( x ) ,f^{(n+1)}(x) ,f(n+1)(x)在 ( a , b ) (a,b) (a,b)内存在,节点 a ⩽ x 0 < x 0 a\leqslant x_0<x_0 a⩽x0<x0, < ⋅ ⋅ ⋅ < x n ⩽ b , L n ( x ) <\cdotp\cdotp\cdotp<x_n\leqslant b,L_n(x) <⋅⋅⋅<xn⩽b,Ln(x)是满足节点条件的插值多项式,则对于任何 x ∈ [ a , b ] x\in[a,b] x∈[a,b],插值余项
R n ( x ) = f ( x ) − L n ( x ) = f ( n + 1 ) ( ξ ) ( n + 1 ) ! ω n + 1 ( x ) . R_n(x)=f(x)-L_n(x)=\frac{f^{(n+1)}(\xi)}{(n+1)!}\omega_{n+1}(x). Rn(x)=f(x)−Ln(x)=(n+1)!f(n+1)(ξ)ωn+1(x).
这里 , ξ ∈ ( a , b ) ,\xi\in(a,b) ,ξ∈(a,b)且依赖于 x x x。
证明:
由给定条件知 , R n ( x ) ,R_n(x) ,Rn(x)在节点 x k ( k = 0 , 1 , ⋅ ⋅ ⋅ , n ) x_k(k=0,1,\cdotp\cdotp\cdotp,n) xk(k=0,1,⋅⋅⋅,n)上为零,即
R n ( x k ) = 0 ( k = 0 , 1 , ⋅ ⋅ ⋅ , n ) R_n\left ( x_k\right ) = 0 \left ( k= 0, 1, \cdotp \cdotp \cdotp , n\right ) Rn(xk)=0(k=0,1,⋅⋅⋅,n)
于是, R n ( x ) = K ( x ) ( x − x 0 ) ( x − x 1 ) ⋯ ( x − x n ) = K ( x ) ω n + 1 ( x ) R_n( x) = K( x) ( x- x_0) ( x- x_1) \cdots ( x- x_n) = K( x) \omega _{n+ 1}( x) Rn(x)=K(x)(x−x0)(x−x1)⋯(x−xn)=K(x)ωn+1(x)
其中 K ( x ) K(x) K(x)是与 x x x 有关的待定函数。
现把 x x x 看成[ a , b ] a,b] a,b]上一个固定点,作函数
φ ( t ) = f ( t ) − L n ( t ) − K ( x ) ( t − x 0 ) ( t − x 1 ) ⋯ ( t − x n ) \varphi(t)=f(t)-L_n(t)-K(x)(t-x_0)(t-x_1)\cdots(t-x_n) φ(t)=f(t)−Ln(t)−K(x)(t−x0)(t−x1)⋯(t−xn)
根据插值条件及余项定义,可知 φ ( t ) \varphi(t) φ(t)在点 x 0 , x 1 , ⋅ ⋅ ⋅ , x n x_0,x_1,\cdotp\cdotp\cdotp,x_n x0,x1,⋅⋅⋅,xn 及 x x x 处均为零,故 φ ( t ) \varphi(t) φ(t)在 [ a , b ] [a,b] [a,b]上有 n + 2 n+2 n+2 个零点。
根据 Rolle 定理 , φ ′ ( t ) ,\varphi^{\prime}(t) ,φ′(t)在 φ ( t ) \varphi(t) φ(t)的两个零点间至少有一个零点,故 φ ′ ( t ) \varphi^\prime(t) φ′(t)在 [ a , b ] [a,b] [a,b]上至少有 n + 1 n+1 n+1 个零点。
对 φ ′ ( t ) \varphi^\prime(t) φ′(t)再应用 Rolle 定理,可知 φ ′ ′ ( t ) \varphi^\prime\prime(t) φ′′(t)在 [ a , b ] [a,b] [a,b]上至少有 n n n 个零点。
依此类推 , φ ( n + 1 ) ( t ) ,\varphi^{(n+1)}(t) ,φ(n+1)(t)在 ( a , b ) (a,b) (a,b)内至少有一个零点,记作 ξ ∈ \xi\in ξ∈ ( a , b ) (a,b) (a,b),使
φ ( n + 1 ) ( ξ ) = f ( n + 1 ) ( ξ ) − ( n + 1 ) ! K ( x ) = 0 \varphi^{(n+1)}\left(\xi\right)^=f^{(n+1)}\left(\xi\right)-\left(n+1\right)!K\left(x\right)=0 φ(n+1)(ξ)=f(n+1)(ξ)−(n+1)!K(x)=0
于是, K ( x ) = f ( n + 1 ) ( ξ ) ( n + 1 ) ! K( x) = \frac {f^{( n+ 1) }( \xi ) }{( n+ 1) ! } K(x)=(n+1)!f(n+1)(ξ), ξ ∈ ( a , b ) \xi \in ( a, b) ξ∈(a,b)且依赖于 x x x
将上式代人
R n ( x ) = K ( x ) ( x − x 0 ) ( x − x 1 ) ⋯ ( x − x n ) = K ( x ) ω n + 1 ( x ) R_n( x) = K( x) ( x- x_0) ( x- x_1) \cdots ( x- x_n) = K( x) \omega _{n+ 1}( x) Rn(x)=K(x)(x−x0)(x−x1)⋯(x−xn)=K(x)ωn+1(x)
就得到余项表达式。