目录

在数值分析中,拉格朗日插值法是以法国十八世纪数学家约瑟夫 · 拉格朗日命名的一种多项式插值方法。许多实际问题中都用函数来表示某种内在联系或规律,而不少函数都只能通过实验和观测来了解。如对实践中的某个物理量进行观测,在若干个不同的地方得到相应的观测值,拉格朗日插值法可以找到一个多项式,其恰好在各个观测的点取到观测到的值。这样的多项式称为拉格朗日(插值)多项式。数学上来说,拉格朗日插值法可以给出一个恰好穿过二维平面上若干个已知点的多项式函数。拉格朗日插值法最早被英国数学家爱德华 · 华林于 1779 年发现 [1],不久后(1783 年)由莱昂哈德 · 欧拉再次发现。1795 年,拉格朗日在其著作《师范学校数学基础教程》中发表了这个插值方法,从此他的名字就和这个方法联系在一起 [2]。

对于给定的若 $n+1$ 个点 $( x _ { 0 } , y _ { 0 } ) , ( x _ { 1 } , y _ { 1 } ) , \cdots , ( x _ { n } , y _ { n } )$,对应于它们的次数不超过 $n$ 的拉格朗日多项式 $L$ 只有一个。如果计入次数更高的多项式,则有无穷个,因为所有与 $L$ 相差 $\lambda ( x - x _ { 0 } ) ( x - x _ { 1 } ) \cdots ( x - x _ { n } )$ 的多项式都满足条件。例子:

已知平面上四个点:$(-9, 5),(-4, 2), (-1, -2), (7, 9)$,拉格朗日多项式:$L(x)$(黑色)穿过所有点。而每个基本多项式:$\mathrm { y } _ { 0 } \ell _ { 0 } ( x ) , \quad \mathrm { y } _ { 1 } \ell _ { 1 } ( x ) , \quad \mathrm { y } _ { 2 } \ell _ { 2 } ( x )$以及 $y_a\ell_3(x)$各穿过对应的一点,并在其它的三个点的 $x$ 值上取零。

定义

对某个多项式函数,已知有给定的 $k + 1$ 个取值点:$\left( x _ { 0 } , y _ { 0 } \right) , \ldots , \left( x _ { k } , y _ { k } \right)$ 其中 $x_j$对应着自变量的位置,而 $y_j$ 对应着函数在这个位置的取值。

假设任意两个不同的 $x_j$ 都互不相同,那么应用拉格朗日插值公式所得到的拉格朗日插值多项式为:

$$ L ( x ) : = \sum _ { j = 0 } ^ { k } y _ { j } \ell _ { j } ( x ) $$

其中每个 $\ell_j(x)$ 为拉格朗日基本多项式(或称插值基函数),其表达式为:

$$ \ell _ { j } ( x ) : = \prod _ { i = 0 , i \neq j } ^ { k } \frac { x - x _ { i } } { x _ { j } - x _ { i } } = \frac { \left( x - x _ { 0 } \right) } { \left( x _ { j } - x _ { 0 } \right) } \cdots \frac { \left( x - x _ { j - 1 } \right) } { \left( x _ { j } - x _ { j - 1 } \right) } \frac { \left( x - x _ { j + 1 } \right) } { \left( x _ { j } - x _ { j + 1 } \right) } \cdot \cdots \frac { \left( x - x _ { k } \right) } { \left( x _ { j } - x _ { k } \right) } $$

[3]拉格朗日基本多项式 $\ell_j(x)$ 的特点是在 $x_j$上取值为 $1$,在其它的点 $x_{i}, i\neq j$上取值为 $0$。

范例