多项式插值法

问题描述

有一函数$y=f(x)$(被插值函数)，已知一系列的点$x_0,x_1,\cdots,x_n$(插值节点)，则

$x$	$x_0$	$x_1$	$\cdots$	$x_n$
$y$	$y_0$	$y_1$	$\cdots$	$y_n$

求简单函数$P(x)$(插值函数)，满足$P(x_i)=y_i$（插值条件）。
即用函数$P(x)$近似表示$f(x)$。

多项式插值

即$P(x)=a_nx^n+\cdots+a_1x+a_0$。

$P(x)的存在性$

$P(x)=a_nx^n+\cdots+a_1x+a_0$
则有，

$$\begin{cases} a_nx_0^n+\cdots+a_1x_0+a_0=y_0\\ a_nx_1^n+\cdots+a_1x_1+a_0=y_1\\ \vdots\\ a_nx_n^n+\cdots+a_1x_n+a_0=y_n\\ \end{cases}$$
上述线性方程组（这里 $a$为变量）的系数矩阵A为：
$$A= \left[ \begin{matrix} x_0^n & x_0^{n-1} & \cdots & x_0 & 1 \\ x_1^n & x_1^{n-1} & \cdots & x_1 & 1 \\ \vdots\\ x_n^n & x_n^{n-1} & \cdots & x_n & 1 \\ \end{matrix} \right]\tag{1}$$
A为范德蒙德行列式，所以|A|≠0，该方程组的解存在并且唯一。

拉格朗日插值多项式

• 当$n=1$时，即$(x_0,y_0),(x_1,y_1)$
  直线点斜式，$y-y_0$=$\dfrac{y_0-y_1}{x_0-x_1}(x-x_0)$
  即，$y=\dfrac{x-x_1}{x_0-x_1}y_0+\dfrac{x-x_0}{x_1-x_0}y_1$
  令$l_0(x)=\frac{x-x_1}{x_0-x_1},l_1(x)=\frac{x-x_0}{x_1-x_0}$
  则$P_1(x)=y_0l_0(x)+y_1l_1(x)$
  且满足
  $$\begin{cases} l_0(x_0)=1,l_1(x_0)=0\\ l_0(x_1)=0,l_1(x_1)=1\\ \end{cases}$$
• 当$n=2$时，$P_2(x)=y_0l_0(x)+y_1l_1(x)+y_2l_2(x)$
  且
  $$\begin{cases} l_0(x_0)=1,l_1(x_0)=0,l_2(x_0)=0\\ l_0(x_1)=0,l_1(x_1)=1,l_2(x_1)=0\\ l_0(x_2)=0,l_1(x_2)=0,l_2(x_2)=1\\ \end{cases}$$
  则 $$\begin{cases} l_0(x)=\dfrac{(x-x_1)(x-x_2)}{(x_0-x_1)(x_0-x_2)}\\ l_1(x)=\dfrac{(x-x_0)(x-x_2)}{(x_1-x_0)(x_1-x_2)}\\ l_2(x)=\dfrac{(x-x_0)(x-x_1)}{(x_2-x_0)(x_2-x_1)}\\ \end{cases}$$
• 当有$n+1$个节点时，即$n=n$时，
  $P_n(x)=y_0l_0(x)+\cdots+y_nl_n(x)$
  则有
  $l_0(x)=\dfrac{(x-x_1)\cdots(x-x_n)}{(x_0-x_1)\cdots(x_0-x_n)}$
  $l_k(x)=\dfrac{(x-x_0)\cdots(x-x_{k-1})(x-x_{k+1})\cdots(x-x_n)}{(x_k-x_0)\cdots(x_k-x_{k-1})(x_k-x_{k+1})\cdots(x_k-x_n)}$（插值基函数）
  记为
  $l_k(x)=\dfrac{w_{n+1}(x)}{(x-x_k)w'_{n+1}(x_k)}$（这里仅仅只是记作），其中$w_{n+1}(x)=(x-x_0)(x-x_1)\cdots(x-x_n)$

插值余项（即误差）

$R_n(x)=f(x)-L_n(x)=\dfrac{f^{(n+1)}(ξ)}{(n+1)!}w_{n+1}(x)$,$a<ξ<b$
证:$R_n(x)$在点$x_0,x_1,\cdots,x_n$处均为零，有
$R_n(x)=K(x)·(x-x_0)(x-x_1)\cdots(x-x_n)=K(x)·w_{n+1}(x)$
构造函数$Φ(t)=f(t)-L_n(t)-K(x)·(t-x_0)(t-x_1)\cdots(t-x_n)$
- 在点$x_0,x_1,\cdots,x_n$处，$f(t)-L_n(t)=0,K(x)·(t-x_0)(t-x_1)\cdots(t-x_n)=0$,则$Φ(t)=0$。
- 在点$x$处，根据$R_n(x)$的定义，$Φ(x)=f(x)-L_n(x)-R_n(x)=0$
则$Φ(t)$在区间$(a,b)$内共有$n+2$个零点，根据罗尔定理，$Φ'(t)$在区间$(a,b)$内至少有$n+1$个零点，依次类推，$Φ^{(n+1)}(t)$在区间$(a,b)$内至少有1个零点，记$a<ξ<b$,使得
$Φ^{(n+1)}(ξ)=f^{(n+1)(ξ)}-0-(n+1)!K(x)=0$
从而$K(x)=\dfrac{f^{(n+1)}(ξ)}{(n+1)!}$
即$R_n(x)=K(x)w_{n+1}(x)=\dfrac{f^{(n+1)}(ξ)}{(n+1)!}w_{n+1}(x)$,$a<ξ<b$
由于ξ的位置不确定，一般而言$R_n(x)$求不出具体的值，若在区间$(a,b)$内，$max|f^{(n+1)}(x)|=M_{n+1}$,则可知其截断误差为：
$|R_n(x)|=\dfrac{|f^{(n+1)}(ξ)|}{(n+1)!}|w_{n+1}(x)|≤\dfrac{M_{n+1}}{(n+1)!}|w_{n+1}(x)|$
- 若$f(x)=x^k且k≤n,则L_n(x)=$$\sum_{i=0}^n x_i^kl_i(x)$$,f^{(n+1)(x)=0}$,插值余项为$R_n(x)=0$，从而，$f(x)=L_n(x)$,即

$$\sum_{i=0}^n x_i^kl_i(x)=x^k,k≤n$$
- 当 $k=0$时，即为，
$$\sum_{i=0}^n l_i(x)=1,k≤n$$

拉格朗日插值多项式优缺点

优点：使用基函数表示，简洁直接；
缺点：当节点增加之后，修改起来不方便。