拉格朗日插值法

1. 拉格朗日插值法简介

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

2.定义

对某个多项式函数，已知有给定的$k+1$个取值点：$(x_0,y_0),\dots ,(x_k,y_k)$，其中$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\ne j}^k\frac{x-x_i}{x_j-x_i}=\frac{(x-x_0)}{(x_j-x_0)}\cdots \frac{(x-x_{j-1})}{(x_j-x_{j-1})}\frac{(x-x_{j+1})}{(x_j-x_{j+1})}\cdots \frac{(x-x_k)}{(x_j-x_k)}$$

拉格朗日基本多项式${\ell }_j(x)$的特点是在$x_j$上取值为1，在其它的点$x_i,i\ne j$上取值为0。

即：对于给定的$n+1$个点$(x_0,y_0),(x_1,y_1),\dots ,(x_n,y_n)$，对应于它们的次数不超过$n$的拉格朗日多项式$L$只有一个。

3. 证明

3.1 存在性证明

对于给定的$k+1$个点：$(x_0,y_0),\dots ,(x_k,y_k)$，拉格朗日插值法的思路是找到一个在一点$x_j$取值为1，而在其他点取值都是0的多项式${\ell }_j(x)$。这样，多项式$y_j{\ell }_j(x)$在点$x_j$取值为$y_j$，而在其他点取值都是0。而多项式$L(x):={\sum }_{j=0}^k{y}_j{\ell }_j(x)$就可以满足：
$$L(x_j)=\sum _{i=0}^k{y}_i{\ell }_i(x_j)=0+0+\cdots +y_j+\cdots +0=y_j$$

在其它点取值为0的多项式容易找到，例如：
$$(x-x_0)\cdots (x-x_{j-1})(x-x_{j+1})\cdots (x-x_k)$$
它在点$x_j$取值为：$(x_j-x_0)\cdots (x_j-x_{j-1})(x_j-x_{j+1})\cdots (x_j-x_k)$。由于已经假定$x_i$两两互不相同，因此上面的取值不等于0。于是，将多项式除以这个取值，就得到一个满足“在$x_j$取值为1，而在其他点取值都是0的多项式”：

$${\ell }_j(x):=\prod _{i=0,i\ne j}^k\frac{x-x_i}{x_j-x_i}=\frac{(x-x_0)}{(x_j-x_0)}\cdots \frac{(x-x_{j-1})}{(x_j-x_{j-1})}\frac{(x-x_{j+1})}{(x_j-x_{j+1})}\cdots \frac{(x-x_k)}{(x_j-x_k)}$$
这就是拉格朗日基本多项式。

3.2 唯一性证明

次数不超过$k$的拉格朗日多项式至多只有一个，因为对任意两个次数不超过$k$的拉格朗日多项式：$P_1$和$P_2$，它们的差$P_1-P_2$在所有$k+1$个点上取值都是0，因此必然是多项式$(x-x_0)(x-x_1)\cdots (x-x_k)$的倍数。因此，如果这个差$P_1-P_2$不等于0，次数就一定不小于$k+1$。但是$P_1-P_2$是两个次数不超过k的多项式之差，它的次数也不超过k。所以$P_1-P_2=0$，也就是说$P_1=P_2$。这样就证明了唯一性。

4.拉格朗日插值法

 思考的问题：对于某一个未知函数的一组观测或者实验数据，寻找一个多项式函数，使这个多项式函数能够过这些点。

例如，已知下面这几个点，我想找到一根穿过它们的曲线：

假设这根曲线是一个二次多项式：$$y=a_0+a_{1}x+a_2{x}^2$$
因为已知有三个点，可以通过下列方程组解出这个二次多项式：

$$\{\begin{array}{l}{y}_1=a_0+a_1{x}_1+a_2{x}_1^2\\ y_3=a_0+a_1{x}_2+a_2{x}_2^2\\ y_3=a_0+a_1{x}_3+a_2{x}_3^2\end{array}$$

不过这里不打算通过解方程来得到这根二次曲线，我们来看看拉格朗日是怎么解出这根曲线的？

4.1 拉格朗日的思考

首先，这是二次曲线。
其次，拉格朗日认为可以通过三根二次曲线相加来达到目标。那么，这三根二次曲线长什么样呢？

第一根二次曲线$f_1(x)$，在点$x_1$处取值为1，其余两点取值为0：

第二根二次曲线$f_2(x)$，在点$x_2$处取值为1，其余两点取值为0：

第三根二次曲线$f_3(x)$，在点$x_3$处取值为1，其余两点取值为0：

这就是朗格朗日需要的三根曲线。

• $y_1{f}_1(x)$可以保证，在$x_1$点处，取值为$y_1$，其余两点取值为0。
• $y_2{f}_2(x)$可以保证，在$x_2$点处，取值为$y_2$，其余两点取值为0。
• $y_3{f}_3(x)$可以保证，在$x_3$点处，取值为$y_3$，其余两点取值为0。

则有：$$f(x)=y_1{f}_1(x)+y_2{f}_2(x)+y_3{f}_3(x)$$

一个动画来展示一一穿过这三个点：

拉格朗日伯爵说，看，这三根曲线就可以组成我在寻找的曲线：

4.2 插值法推导

用符号来表示这三根曲线：$f_i(x_j),i=1,2,3,j=1,2,3$。

首先，$f_i$必须是二次函数。其次，需要满足的条件是
$$f_i(x_j)=\{\begin{array}{ll}1,& \text{i=j}\\ 0,& \text{i=j}\end{array}$$

那么，如下构造$f_1(x)$很显然可以满足上述条件：
$$f_1(x)=\frac{(x-x_2)(x-x_3)}{(x_1-x_2)(x_1-x_3)}$$

不失一般性，写成：
$$f_i(x)=\prod _{i\ne j}^{1\le j\le 3}\frac{x-x_j}{x_i-x_j}$$

因此，最终得到：
$$f(x)=\sum _{i=1}^3{y}_i{f}_i(x)$$

这就是拉格朗日插值法。上面的思路要推广到更多点的插值也非常容易。

4.2 解答开始提出的方程组

如何解下面这个三元一次方程组：
$$\{\begin{array}{l}{y}_1=a_0+a_1{x}_1+a_2{x}_1^2\\ y_3=a_0+a_1{x}_2+a_2{x}_2^2\\ y_3=a_0+a_1{x}_3+a_2{x}_3^2\end{array}$$

$(x_1,y_1),(x_2,y_2),(x_3,y_3)$三个点是已知的，因此上面实际是一个线性方程组：
$$\underset{\vec{y}}{\underbrace{\left(\begin{array}{c}{y}_1\\ y_2\\ y_3\end{array}\right)}}=\underset{A}{\underbrace{\left(\begin{array}{ccc}1& x_1& x_1^2\\ 1& x_2& x_2^2\\ 1& x_3& x_3^2\end{array}\right)}}\underset{\vec{a}}{\underbrace{\left(\begin{array}{c}{a}_0\\ a_1\\ a_2\end{array}\right)}}$$

简写即：
$$\vec{y}=A\vec{a}$$

$A$就是范德蒙矩阵，如果$|A|\ne 0$，那么方程就只有唯一的解。

下面就求$|A|$，先把$r_2-r_1,r_3-r_1$，（$r_i$表示第$i$行）得到：
$$|A|=\left|\begin{array}{ccc}1& x_1& x_1^2\\ 0& x_2-x_1& x_2^2-x_1^2\\ 0& x_3-x_1& x_3^2-x_1^2\end{array}\right|=\left|\begin{array}{cc}{x}_2-x_1& x_2^2-x_1^2\\ x_3-x_1& x_3^2-x_1^2\end{array}\right|=(x_3-x_1)(x_3-x_2)(x_2-x_1)$$

由于$x_1,x_2,x_3$都不相等，所以$|A|\ne 0$。

到此为止，也说明“对于给定的$n+1$个点$(x_0,y_0),(x_1,y_1),\dots ,(x_n,y_n)$，对应于它们的次数不超过$n$的拉格朗日多项式$L$只有一个。”

5. 参考

[1] https://www.zhihu.com/question/58333118
[2] https://www.cnblogs.com/Wolfycz/p/10624678.html

欢迎关注公众号：区块链之美 。
致力于区块链技术研究，传播区块链技术和解决方案、区块链应用落地、区块链行业动态等。