数值计算(二)——插值计算(1)拉格朗日插值法

插值问题

设函数$f$定义在$[a,b]$上，设$x_1,x_2,...,x_n\in [a,b]$是$n+1$个相异节点，令$y_i=f(x_i),i=0,1,...,n$插值法就是构造一个便于计算的简单函数$\phi$来近似地替代$f$，并满足：

$\phi (x_i)=f(x_i),i=0,1,...,n$

对任意的$x\in [a,b],\phi (x)$作为$f(x)$的近似值，通常将$f$记作被插值函数，$x_1,x_2,...,x_n$称为插值节点，$\phi$称为插值函数，$\phi (x_i)=f(x_i),i=0,1,...,n$为插值条件，用多项式作为插值函数的插值法叫做多项式插值法，相应的多项式称为插值多项式。

拉格朗日插值法

设函数$f$的$n+1$个相异节点$x_1,x_2,...,x_n$上的值$f(x_0),f(x_1),f(x_2),...,f(x_n)$是已知的，令${\wp }_n$使得
$P(x_i)=f(x_i),i=0,1,...,n$，
对于给定的$n+1$个节点$x_i,i=0,1,...,n$记
$l_i(x)=\prod _{i=0,j\ne i}^n\frac{x-x_j}{x_i-x_j},i=0,1,2,...,n,l_i,i=0,1,2,...,n$，是n次多项式并满足：
$l_i(x_i)=\{\begin{array}{ll}1& j=i\\ 0& j\ne i\end{array}$，令$L_n(x)=\sum _{i=0}^{n}f(x_i)l_i(x),$显然，$L_n={\wp }_n$并满足插值条件：$L_n(x_j)=f(x_i),j=0,1,2,...n.$,写成上面连乘形式的n次多项式$L_n$被称为n次$Lagrange$插值多项式，$l_i\in {\wp }_n(i=0,1,2,...n)$称为n次$Lagrange$插值基函数。

>定理1：设$p\in {\wp }_n,p$的零点大于n$(p的l重零点x^{\ast }计为l个零点)$，那么p恒为零。

>定理2：给定$n+1$个相异的节点$x_0,x_1,...x_n\in [a,b]和n+1个$函数值$f(x_0),f(x_1),...,f(x_n)$那么存在唯一$n$次多项式$p\in {\wp }_n$满足插值条件$P(x_i)=f(x_i),i=0,1,...,n$，
证明过程：

插值余项：

设$x_0,x_1,...x_n\in [a,b]$是相异节点，$f\in C^{n+1}[a,b],L_n$为满足插值条件的n次插值多项式，对于任何$x\in [a,b]$存在$\xi =\xi (x)\in (a,b)$使得$R_n(x)=\frac{f^{n+1}(\xi )}{(n+1)!}{\omega }_{n+1}(x)$,其中${\omega }_{n+1}(x)=(x-x_0)(x-x_1)...(x-x_n)$.

推论

推论一：在上述条件下，若$\underset{a<x<b}{\max }\mid f^{n+1}(x)\mid \le M$,那么有$\mid R_{(}x)\mid \le \frac{M}{(n+1)!}\mid {\omega }_{n+1}(x)\mid$
推论二：设$a=x_0<x_1<x_2...<x_n=b,h=\underset{1\le j\le n}{\max }(x_j-x_{j-1}),f\in C^{n+1}[a,b],L_n$为$f$的$n$次插值多项式，那么有：$\Vert f-L_n{\Vert }_{\infty }\le \frac{h^{n+1}}{4(n+1)}\Vert f^{n+1}{\Vert }_{\infty }$
有兴趣的自己去证明吧，哈哈哈哈！！！！！

整道基础题巩固基础知识点

欢乐的时光总是短暂的，让我们下一次再见！！！

good good study,day day up! (study hard, improve every day)

预知后事，请听下回分解！！！！