无锡开发网站建设搜外
函数插值
已知函数 f ( x ) f(x) f(x)在区间[a,b]上n+1个互异节点 { x i } i = 0 n \{{x_i}\}_{i=0}^{n} {xi}i=0n处的函数值 { y i } i = 0 n \{{y_i}\}_{i=0}^{n} {yi}i=0n，若函数集合 Φ \Phi Φ中函数 ϕ ( x ) \phi(x) ϕ(x)满足条件
ϕ ( x i ) = y i ( i = 0 , 1 , 2 , ⋯ , n ) \phi\left(x_{i}\right)=y_{i} \quad(i=0,1,2, \cdots, n) ϕ(xi)=yi(i=0,1,2,⋯,n)
称 ϕ ( x ) \phi(x) ϕ(x)为 f ( x ) f(x) f(x)在 Φ \Phi Φ中关于节点 { x i } i = 0 n \{{x_i}\}_{i=0}^{n} {xi}i=0n的插值函数。
R n ( x ) = f ( x ) − ϕ ( x ) R_{n}(x)=f(x)-\phi(x) Rn(x)=f(x)−ϕ(x)为插值余项。
选取 1 , x , x 2 , . . . , x n 1,x,x^2,...,x^n 1,x,x2,...,xn作为n次多项式空间的一组基函数，对于多项式插值，有
R n ( x ) = f ( x ) − ϕ ( x ) = f ( n + 1 ) ( ξ ) ( n + 1 ) ! ω n + 1 ( x ) R_{n}(x)=f(x)-\phi(x)=\frac{f^{(n+1)}(\xi)}{(n+1)!} \omega_{n+1}(x) Rn(x)=f(x)−ϕ(x)=(n+1)!f(n+1)(ξ)ωn+1(x)
其中 ξ = ξ ( x ) ∈ ( a , b ) , ω n + 1 ( x ) = ( x − x 0 ) ( x − x 1 ) ⋯ ( x − x n ) \xi=\xi(x) \in(a, b), \omega_{n+1}(x)=\left(x-x_{0}\right)\left(x-x_{1}\right) \cdots\left(x-x_{n}\right) ξ=ξ(x)∈(a,b),ωn+1(x)=(x−x0)(x−x1)⋯(x−xn)
Lagrange插值
线性插值
过两点 ( x 0 , y 0 ) (x_0,y_0) (x0,y0)和 ( x 1 , y 1 ) (x_1,y_1) (x1,y1)的直线的方程为
y − y 0 = y 1 − y 0 x 1 − x 0 ( x − x 0 ) y - y_0 =\frac{y_1 - y_0}{x_1 - x_0}( x - x_0 ) y−y0=x1−x0y1−y0(x−x0)
整理可得，
y = x − x 1 x 0 − x 1 y 0 + x − x 0 x 1 − x 0 y 1 y =\frac{x - x_1}{x_0 - x_1}y_0 + \frac{x - x_0}{x_1 - x_0}y_1 y=x0−x1x−x1y0+x1−x0x−x0y1
引入记号，
l 0 ( x ) = x − x 1 x 0 − x 1 , l 1 ( x ) = x − x 0 x 1 − x 0 l_0( x )=\frac{x - x_1}{x_0 - x_1},\quad l_1( x )=\frac{x - x_0}{x_1 - x_0} l0(x)=x0−x1x−x1,l1(x)=x1−x0x−x0
易知 l 0 l_0 l0和 l 1 l_1 l1线性无关，
则有 L 1 ( x ) = l 0 ( x ) y 0 + l 1 ( x ) y 1 L_1( x ) = l_0( x ) y_0 + l_1( x ) y_1 L1(x)=l0(x)y0+l1(x)y1为Lagrange插值多项式。
对于二次多项式，可写出
L 2 ( x ) = l 0 ( x ) y 0 + l 1 ( x ) y 1 + l 2 ( x ) y 2 L_2( x ) = l_0( x ) y_0 + l_1( x ) y_1 + l_2( x ) y_2 L2(x)=l0(x)y0+l1(x)y1+l2(x)y2
需满足条件，
( l 0 ( x 0 ) = 1 l 0 ( x 1 ) = 0 l 0 ( x 2 ) = 0 ) \left( \begin{matrix}l_0( x_0 ) = 1 \\l_0( x_1 ) = 0 \\l_0( x_2 ) = 0 \end{matrix}\right) l0(x0)=1l0(x1)=0l0(x2)=0
( l 1 ( x 0 ) = 0 l 1 ( x 1 ) = 1 l 1 ( x 2 ) = 0 ) \left(\begin{matrix}l_1( x_0 ) = 0 \\l_1( x_1 ) = 1 \\l_1( x_2 ) = 0 \end{matrix}\right) l1(x0)=0l1(x1)=1l1(x2)=0
( l 2 ( x 0 ) = 0 l 2 ( x 1 ) = 0 l 2 ( x 2 ) = 1 ) \left(\begin{matrix}l_2( x_0 ) = 0 \\l_2( x_1 ) = 0 \\l_2( x_2 ) = 1 \end{matrix}\right) l2(x0)=0l2(x1)=0l2(x2)=1
对于n+1个节点的情形，
L n ( x ) = l 0 ( x ) y 0 + l 1 ( x ) y 1 + l 2 ( x ) y 2 + ⋯ + l n ( x ) y n = ∑ i = 0 n l i ( x ) y i L_n( x ) = l_0( x ) y_0 + l_1( x ) y_1 + l_2( x ) y_2 + \cdots + l_n( x ) y_n = \sum_{i = 0}^n l_i( x ) y_i Ln(x)=l0(x)y0+l1(x)y1+l2(x)y2+⋯+ln(x)yn=i=0∑nli(x)yi
满足条件，
l i ( x j ) = { 1 , j = i , 0 , 0 ⩽ j ⩽ n , j ≠ i , l_i( x_j ) = \left\{\begin{matrix}1 ,& j = i ,\\0 ,& 0 \leqslant j \leqslant n ,j \neq i ,\end{matrix}\right. li(xj)={1,0,j=i,0⩽j⩽n,j=i,
l i ( x ) l_i(x) li(x)为
l i ( x ) = ( x − x 0 ) ( x − x 1 ) ⋯ ( x − x i − 1 ) ( x − x i + 1 ) ⋯ ( x − x n ) ( x i − x 0 ) ( x i − x 1 ) ⋯ ( x i − x i − 1 ) ( x i − x i + 1 ) ⋯ ( x i − x n ) . l_i( x )=\frac{( x - x_0 ) ( x - x_1 )\cdots( x - x_{i-1} ) ( x - x_{i+1} )\cdots( x - x_n )}{( x_i - x_0 ) ( x_i - x_1 )\cdots( x_i - x_{i-1} ) ( x_i - x_{i+1} )\cdots( x_i - x_n )}. li(x)=(xi−x0)(xi−x1)⋯(xi−xi−1)(xi−xi+1)⋯(xi−xn)(x−x0)(x−x1)⋯(x−xi−1)(x−xi+1)⋯(x−xn).
Newton插值法
过两点 ( x 0 , y 0 ) (x_0,y_0) (x0,y0)和 ( x 1 , y 1 ) (x_1,y_1) (x1,y1)的直线方程记为
p 1 ( x ) = y 0 + y 1 − y 0 x 1 − x 0 ( x − x 0 ) = y 0 + c 1 ( x − x 0 ) . p_1( x ) = y_0 + \frac{y_1 - y_0}{x_1 - x_0}( x - x_0 ) = y_0 + c_1( x - x_0 ) . p1(x)=y0+x1−x0y1−y0(x−x0)=y0+c1(x−x0).
增加一个节点，
p 2 ( x ) = p 1 ( x ) + c 2 ( x − x 0 ) ( x − x 1 ) . p_2( x )=p_1( x ) + c_2( x - x_0 ) ( x - x_1 ) . p2(x)=p1(x)+c2(x−x0)(x−x1).
节点由k个增加到k+1个，满足
p k ( x ) = p k − 1 ( x ) + c k ( x − x 0 ) ( x − x 1 ) ⋯ ( x − x k − 1 ) . p_k( x )=p_{k-1}( x ) + c_k( x - x_0 ) ( x - x_1 )\cdots( x - x_{k-1} ) . pk(x)=pk−1(x)+ck(x−x0)(x−x1)⋯(x−xk−1).
由此可得，
p k ( x ) = y 0 + c 1 ( x − x 0 ) + ⋯ + c k ( x − x 0 ) ( x − x 1 ) ⋯ ( x − x k − 1 ) . p_{k}( x )=y_{0}+c_{1}( x-x_{0} )+\cdots+c_{k}( x-x_{0} ) ( x-x_{1} )\cdots( x-x_{k-1} ) . pk(x)=y0+c1(x−x0)+⋯+ck(x−x0)(x−x1)⋯(x−xk−1).
是以 1 , x − x 0 , ( x − x 0 ) ( x − x 1 ) , ⋯ , ( x − x 0 ) ( x − x 1 ) ⋯ ( x − x k − 1 ) 1,x-x_0,(x-x_0)(x-x_1),\cdots,(x-x_0)(x-x_1)\cdots(x-x_{k-1}) 1,x−x0,(x−x0)(x−x1),⋯,(x−x0)(x−x1)⋯(x−xk−1)为基函数的插值多项式，这种方法被称为Newton插值法。
分段线性插值
高次多项式插值是不稳定的。
三次样条插值
