数值分析复习：Lagrange插值

本篇文章适合个人复习翻阅，不建议新手入门使用

拉格朗日（Lagrange）插值

拉格朗日（Lagrange）插值，又称L插值

插值条件

$n+1$ 个插值节点$x_0,x_1,\dots ,x_n$处函数值相同（插值节点均在区间$[a,b]$上）

基函数

{ l i ( x ) } i = 0 n \{l_i(x)\}_{i=0}^n {li​(x)}i=0n​，其中$$l_i(x)=\frac{\prod _{j\ne i}(x-x_j)}{\prod _{j\ne i}(x_i-x_j)}$$其性质：$l_i(x_j)={\delta }_{ij}$

插值多项式

$$\begin{array}{rl}{L}_n(x)& =\sum _{i=0}^n{y}_i{l}_i(x)\\ & =\sum _{i=0}^n{y}_i\frac{{\omega }_{n+1}(x)}{(x-x_i){\omega }_{n+1}^{\prime }(x_i)}\end{array}$$其中 ${\omega }_{n+1}(x)=(x-x_0)\cdots (x-x_n)$ 为节点多项式

余项估计

$$|R_n(x)|\le \frac{C_{n+1}}{(n+1)!}|{\omega }_{n+1}(x)|$$其中 $R_n(x)=f(x)-L_n(x)$为插值余项，$C_{n+1}$为$f(x)$的$n+1$阶导函数的界

证明：（构造辅助函数+罗尔（Rolle）中值定理）
设$$E(t)=R_n(t)-\frac{R_n(x)}{{\omega }_{n+1}(x)}{\omega }_{n+1}(t)$$则 $E(t)$ 在$[a,b]$上$n+1$次连续可微，且有$n+2$个零点$x_0,x_1,\dots ,x_n,x$

由Rolle中值定理，存在 $\xi \in [a,b],s.t.E^{(n+1)}(\xi )=0$，又
$$\begin{array}{rl}{E}_{n+1}(\xi )& =R_n^{(n+1)}(\xi )-\frac{R_n(x)}{{\omega }_{n+1}(x)}(n+1)!\\ & =f^{(n+1)}(\xi )-\frac{R_n(x)}{{\omega }_{n+1}(x)}(n+1)!\\ & =0\end{array}$$
故$$|R_n(x)|=|\frac{f^{(n+1)}(\xi )}{(n+1)!}\cdot {\omega }_{n+1}(x)|\le \frac{C_{n+1}}{(n+1)!}|{\omega }_{n+1}(x)|$$

参考书籍：《数值分析》李庆扬 王能超 易大义 编