浅谈拉格朗日插值

jcx106

2026-06-06 16:59:58

发布于：浙江

17阅读

0回复

0点赞

一、数学基础

理解拉格朗日插值，首先需要明确多项式函数的基本性质。对于一个 $m$ 次多项式 $f(x)$ ，其一般形式可表示为：

$f(x)=a_0+a_1x+a_2x^2+\cdots+a_mx^m$

其中 $a_i$ 为多项式的系数， $x$ 为自变量。该函数共有 $m+1$ 个待定系数，根据多项式插值的存在唯一性定理：给定平面上 $m+1$ 个横坐标互不相同的点 $(x_0, y_0),(x_1, y_1),\cdots,(x_m, y_m)$ ，我们可以唯一确定一个次数不超过 $m$ 的多项式函数 $f(x)$ ，使其图像恰好经过这 $m+1$ 个点。这正是拉格朗日插值法的核心理论基础。

二、基本原理

拉格朗日插值的核心思想是构造基函数并进行线性组合：我们可以将待求的插值函数 $P(x)$ 拆解为多个简单的基函数的加权和。我们为每个已知点 $(x_i, y_i)$ 构造一个对应的基函数 $l_i(x)$ ，要求这个基函数满足：在 $x=x_i$ 处取值为 $1$ ，而在其余所有已知点 $x_j\ (j\neq i)$ 处取值为 $0$ 。

这样一来，我们就可以把插值函数 $P(x)$ 表示为这些基函数的线性组合：

$P(x)=y_0l_0(x)+y_1l_1(x)+y_2l_2(x)+\cdots+y_ml_m(x)$

此时，求插值函数的问题就转化为了求解这一系列拉格朗日基函数的问题。

三、拉格朗日基函数

接下来我们来具体构造满足要求的基函数 $l_i(x)$ 。根据多项式因式分解的性质：若一个 $m$ 次多项式在 $m$ 个互不相同的点处取值为0，那么它一定可以写成 $f(x)=a\prod_{j=0}^m (x-b_j)$ 的形式，其中 $a$ 为非零常数， $b_0,b_1,...,b_m$ 就是该多项式的根。

基于这个性质，我们可以先构造一个初步的多项式：

$l_i'(x)=(x-x_0)(x-x_1)\cdots(x-x_{i-1})(x-x_{i+1})\cdots(x-x_m)$

这个多项式的特点是：对于所有 $j\neq i$ ，都有 $l_i'(x_j)=0$ ，正好满足基函数在其余点处取0的要求。但我们很快会发现，这个多项式在 $x_i$ 处的取值并不一定是1，因此我们需要对它进行归一化处理。

我们可以通过代入计算确定这个归一化常数：将 $x=x_i$ 代入上述初步多项式，可得：

$l_i'(x_i)=(x_i-x_0)(x_i-x_1)\cdots(x_i-x_{i-1})(x_i-x_{i+1})\cdots(x_i-x_m)$

为了让基函数在 $x_i$ 处的取值为1，我们只需要将初步多项式除以这个值即可。因此，最终的拉格朗日基函数可以表示为：

$l_i(x) = \frac{l_i'(x)}{l_i'(x_i)} = \prod_{\substack{0 \leq j \leq m \\ j \neq i}} \frac{x - x_j}{x_i - x_j}$

得到基函数后，我们将其代入插值函数的表达式，就可以得到完整的拉格朗日插值公式：

$P(x)=\sum_{i=0}^m y_i \cdot l_i(x) = \sum_{i=0}^m y_i \prod_{\substack{0 \leq j \leq m \\ j \neq i}} \frac{x - x_j}{x_i - x_j}$

四、示例验证

我们可以通过一个具体示例来验证上述公式的正确性：假设我们需要构造一个多项式，使其穿过点 $(0,1),(1,2),(2,1)$ 。根据插值的基本理论，3个横坐标互不相同的点可以唯一确定一个次数不超过2的多项式，也就是二次函数。

我们首先分别计算每个点对应的基函数：

对于点 $(0,1)$ ，基函数为：
$l_0(x)=\frac{(x-1)(x-2)}{(0-1)(0-2)}=\frac{1}{2}x^2-\frac{3}{2}x+1$
对于点 $(1,2)$ ，基函数为：
$l_1(x)=\frac{(x-0)(x-2)}{(1-0)(1-2)}=-x^2+2x$
对于点 $(2,1)$ ，基函数为：
$l_2(x)=\frac{(x-0)(x-1)}{(2-0)(2-1)}=\frac{1}{2}x^2-\frac{1}{2}x$

接下来我们将基函数与对应的函数值加权求和，得到插值多项式：

$\begin{align*} P(x) &= 1\cdot l_0(x) + 2\cdot l_1(x) + 1\cdot l_2(x) \\ &= \left(\frac{1}{2}x^2-\frac{3}{2}x+1\right) + 2\left(-x^2+2x\right) + \left(\frac{1}{2}x^2-\frac{1}{2}x\right) \\ &= -x^2+2x+1 \end{align*}$