n次多项式空间是有限维线性空间
◉ { x 0 , ⋯ , x n } \{x_0,\cdots,x_n\} { x0,⋯,xn}互异, k ∈ { 0 , ⋯ , n } k \in \{0,\cdots,n\} k∈{ 0,⋯,n}, 求证
- ∑ i = 0 n x i k l i ( x ) ≡ x k , k ∈ { 0 , ⋯ , n } \sum\limits_{i=0}^{n} x_i^k l_i(x) \equiv x^k, \quad k \in \{0,\cdots,n\} i=0∑nxikli(x)≡xk,k∈{ 0,⋯,n}
- ∑ i = 0 n ( x i − x ) k l i ( x ) ≡ 0 , k ∈ { 1 , ⋯ , n } \sum\limits_{i=0}^{n} (x_i-x)^k l_i(x) \equiv 0, \quad k \in \{1,\cdots,n\} i=0∑n(xi−x)kli(x)≡0,k∈{ 1,⋯,n}
1.
x k x^k xk是一个不超过 n n n次的多项式.
大前提: 给定 x k x^k xk上的 n + 1 n+1 n+1个相异的点进行插值, 由插值的唯一性可知插值多项式恒等于 x k x^k xk.
小前提: 给定 x k x^k xk上的 n + 1 n+1 n+1个相异的点( x i ( i ∈ { 0 , ⋯ , n } ) x_i(i \in \{0,\cdots,n\}) xi(i∈{ 0,⋯,n}))进行(Lagrange)插值, ∑ i = 0 n x i k l i ( x ) \sum\limits_{i=0}^{n} x_i^k l_i(x) i=0∑nxikli(x)是插值多项式.
结论: ∑ i = 0 n x i k l i ( x ) \sum\limits_{i=0}^{n} x_i^k l_i(x) i=0∑nxikli(x)恒等于 x k x^k xk.
2.
k ∈ { 0 , ⋯ , n } k \in \{0,\cdots,n\} k∈{
0,⋯,n}, ∑ i = 0 n ( x i − t ) k l i ( x ) ≡ ( x − t ) k \sum\limits_{i=0}^{n} (x_i-t)^k l_i(x) \equiv (x-t)^k i=0∑n(xi−t)kli(x)≡(x−t)k.
令 t = x t = x t=x, ∑ i = 0 n ( x i − x ) k l i ( x ) ≡ ( x − x ) k ≡ 0 \sum\limits_{i=0}^{n} (x_i-x)^k l_i(x) \equiv (x-x)^k \equiv 0 i=0∑n(xi−x)kli(x)≡(x−x)k≡0.
插值余项
◉ 给出分段两点三次Hermite插值误差限的理论分析.
2 n + 1 2n+1 2n+1次插值余项
a = x 0 < ⋯ < x n = b a=x_0<\cdots<x_n=b a=x0<⋯<xn=b
π ^ ( x ) = ∏ i = 0 n ( x − x i ) 2 \hat\pi(x) = \prod\limits_{i=0}^{n}(x-x_{i})^2 π^(x)=i=0∏n(x−xi)2
R 2 n + 1 ( x ) = f ( x ) − H 2 n + 1 ( x ) R_{2n+1}(x) = f(x) - H_{2n+1}(x) R2n+1(x)=f(x)−H2n+1(x)
K ( x ) = R 2 n + 1 ( x ) π ^ ( x ) K(x) = \frac{R_{2n+1}(x)}{\hat\pi(x)} K(x)=π^(x)R2n+1(x)
ϕ ( t ) = R 2 n + 1 ( t ) − K ( x ) π ^ ( t ) \phi(t) = R_{2n+1}(t) - K(x)\hat\pi(t) ϕ(t)=R2n+1(t)−K(x)π^(t)
注意到 ϕ ( t ) \phi(t) ϕ(t)有 n + 2 n+2 n+2个零点 { x 0 , ⋯ , x n , x } \{x_0,\cdots,x_n,x\} { x0,⋯,xn,x}, ϕ ′ ( t ) \phi'(t) ϕ′(t)也有 n + 1 n+1 n+1个零点 { x 0 , ⋯ , x n } \{x_0,\cdots,x_n\} { x0,⋯,xn}, 所以 ϕ ( t ) \phi(t) ϕ(t)有 2 n + 3 2n+3 2n+3个根( n + 1 n+1 n+1个二重根 { x 0 , ⋯ , x n } \{x_0,\cdots,x_n\} { x0,⋯,xn}, 1 1 1个一重根 { x } \{x\} { x}). 迭用Rolle定理得到 ϕ ( 2 n + 2 ) ( ξ ) = 0 \phi^{(2n+2)}(\xi)=0 ϕ(2n+2)(ξ)=0(其中 ξ ∈ ( a , b ) \xi\in(a,b) ξ∈(a,b)).
R 2 n + 1 ( 2 n + 2 ) ( ξ ) = f ( 2 n + 2 ) ( ξ ) R_{2n+1}^{(2n+2)}(\xi) = f^{(2n+2)}(\xi) R2n+1(2n+2)(ξ)=f(2n+2)(ξ)
π ^ ( 2 n + 2 ) ( ξ ) = ( 2 n + 2 ) ! \hat\pi^{(2n+2)}(\xi) = (2n+2)! π^(2n+2)(ξ)=(2n+2)!
ϕ ( 2 n + 2 ) ( ξ ) = f ( 2 n + 2 ) ( ξ ) − K ( x ) ( 2 n + 2 ) ! = 0 \phi^{(2n+2)}(\xi) = f^{(2n+2)}(\xi) - K(x)(2n+2)! = 0 ϕ(2n+2)(ξ)=f(2n+2)(ξ)−K(x)(2n+2)!=0
K ( x ) = R 2 n + 1 ( x ) π ^ ( x ) = f ( 2 n + 2 ) ( ξ ) ( 2 n + 2 ) ! K(x) = \frac{R_{2n+1}(x)}{\hat\pi(x)} = \frac{f^{(2n+2)}(\xi)}{(2n+2)!} K(x)=π^(x)R2n+1(x)=(2n+2