三次多项式插值求解步骤详解

最新推荐文章于 2026-07-16 20:58:45 发布

原创最新推荐文章于 2026-07-16 20:58:45 发布 · 571 阅读

7 ·

本内容遵循CC 4.0 BY-SA版权协议

GEO检测

好的，我们来详细解释一下这个求解三次样条（实际上是一个全局三次多项式）方程的完整计算步骤。

这个问题本质上是给定四个点，求一个唯一的三次多项式函数 $y = a + bx + cx^2 + dx^3$ 使其恰好通过这四个点。这是一个经典的“多项式插值”问题。

第一步：建立方程组

我们假设所求的曲线方程为：

$$
y = a + bx + cx^2 + dx^3
$$

我们有四个已知点：$P_1(-1, 2)$, $P_2(0, 0)$, $P_3(1, -2)$, 和 $P_4(2, 0)$。因为曲线要通过这些点，所以每个点的坐标 $(x, y)$ 都必须满足上面的方程。我们将每个点的坐标代入方程，得到四个线性方程：

对于点 $P_1(-1, 2)$:

$$
2 = a + b(-1) + c(-1)^2 + d(-1)^3
$$

$$
2 = a - b + c - d
$$

对于点 $P_2(0, 0)$:

$$
0 = a + b(0) + c(0)^2 + d(0)^3
$$

$$
0 = a
$$
（这一步直接给出了 $a$ 的值，可以简化后续计算）

对于点 $P_3(1, -2)$:

$$
-2 = a + b(1) + c(1)^2 + d(1)^3
$$

$$
-2 = a + b + c + d
$$

对于点 $P_4(2, 0)$:

$$
0 = a + b(2) + c(2)^2 + d(2)^3
$$

$$
0 = a + 2b + 4c + 8d
$$

现在，我们得到了一个包含四个未知数 ($a, b, c, d$) 的四个线性方程组成的方程组。

第二步：写成矩阵形式

为了更方便地求解，我们将上面的方程组写成矩阵乘法的形式 $Y = C \cdot A$。

$Y$ (因变量向量): 由所有点的 $y$ 坐标组成。

$$
Y = \begin{bmatrix} 2 \ 0 \ -2 \ 0 \end{bmatrix}
$$

$A$ (系数向量): 由我们要求解的多项式系数组成。

$$
A = \begin{bmatrix} a \ b \ c \ d \end{bmatrix}
$$

$C$ (设计矩阵或范德蒙德矩阵): 每一行对应一个数据点，每一列对应多项式的一项 ($1, x, x^2, x^3$)。
- 第一行对应 $P_1(-1, 2)$: $[1, -1, (-1)^2, (-1)^3] = [1, -1, 1, -1]$
- 第二行对应 $P_2(0, 0)$: $[1, 0, 0^2, 0^3] = [1, 0, 0, 0]$
- 第三行对应 $P_3(1, -2)$: $[1, 1, 1^2, 1^3] = [1, 1, 1, 1]$
- 第四行对应 $P_4(2, 0)$: $[1, 2, 2^2, 2^3] = [1, 2, 4, 8]$
所以，矩阵 $C$ 为：

$$
C = \begin{bmatrix} 1 & -1 & 1 & -1 \ 1 & 0 & 0 & 0 \ 1 & 1 & 1 & 1 \ 1 & 2 & 4 & 8 \end{bmatrix}
$$

这样，整个方程组就可以简洁地表示为：

$$
\begin{bmatrix} 2 \ 0 \ -2 \ 0 \end{bmatrix} = \begin{bmatrix} 1 & -1 & 1 & -1 \ 1 & 0 & 0 & 0 \ 1 & 1 & 1 & 1 \ 1 & 2 & 4 & 8 \end{bmatrix} \begin{bmatrix} a \ b \ c \ d \end{bmatrix}
$$
这与图片中的矩阵形式完全一致。

第三步：求解系数矩阵 A

我们的目标是求出向量 $A$。根据矩阵方程 $Y = C \cdot A$，我们可以通过左乘 $C$ 的逆矩阵 $C^{-1}$ 来求解 $A$：

$$
A = C^{-1} \cdot Y
$$

计算 $C^{-1}$ (矩阵求逆) 是这个过程中最复杂的一步，通常由计算机完成。手动计算可以使用高斯-约旦消元法或伴随矩阵法。这里我们直接给出结果（为了精确，我们用分数表示，图片中是近似小数）：

经过计算，$C$ 的逆矩阵 $C^{-1}$ 大约是：

$$
C^{-1} \approx \begin{bmatrix} 0 & 1 & 0 & 0 \ -1 & 0 & 1 & 0 \ 0.5 & -1.5 & 0.5 & 0 \ -0.1667 & 0.5 & -0.5 & 0.1667 \end{bmatrix}
$$
更精确的分数形式是：

$$
C^{-1} = \begin{bmatrix} 0 & 1 & 0 & 0 \ -1 & 0 & 1 & 0 \ \frac{1}{2} & -\frac{3}{2} & \frac{1}{2} & 0 \ -\frac{1}{6} & \frac{1}{2} & -\frac{1}{2} & \frac{1}{6} \end{bmatrix}
$$

计算 $A = C^{-1} \cdot Y$:

$$
\begin{bmatrix} a \ b \ c \ d \end{bmatrix} = \begin{bmatrix} 0 & 1 & 0 & 0 \ -1 & 0 & 1 & 0 \ \frac{1}{2} & -\frac{3}{2} & \frac{1}{2} & 0 \ -\frac{1}{6} & \frac{1}{2} & -\frac{1}{2} & \frac{1}{6} \end{bmatrix} \begin{bmatrix} 2 \ 0 \ -2 \ 0 \end{bmatrix}
$$

进行矩阵乘法：

$a = (0)(2) + (1)(0) + (0)(-2) + (0)(0) = 0$
$b = (-1)(2) + (0)(0) + (1)(-2) + (0)(0) = -2 - 2 = -4$
$c = (\frac{1}{2})(2) + (-\frac{3}{2})(0) + (\frac{1}{2})(-2) + (0)(0) = 1 - 1 = 0$
$d = (-\frac{1}{6})(2) + (\frac{1}{2})(0) + (-\frac{1}{2})(-2) + (\frac{1}{6})(0) = -\frac{1}{3} + 1 = \frac{2}{3}$

所以，我们得到精确的系数解：

$$
A = \begin{bmatrix} a \ b \ c \ d \end{bmatrix} = \begin{bmatrix} 0 \ -4 \ 0 \ \frac{2}{3} \end{bmatrix}
$$

图片中给出的结果是近似值：

$$
\begin{bmatrix} a \ b \ c \ d \end{bmatrix} = \begin{bmatrix} 0 \ -2.67 \ 0 \ 0.67 \end{bmatrix}
$$
这里有一个明显的错误。根据我们的计算，$b$ 应该是 $-4$，而不是 $-2.67$。$d = \frac{2}{3} \approx 0.67$ 是正确的。让我们用正确的系数来验证。

第四步：写出最终方程并验证

将我们计算出的精确系数 $a=0, b=-4, c=0, d=\frac{2}{3}$ 代入初始方程 $y = a + bx + cx^2 + dx^3$：

$$
y = 0 + (-4)x + (0)x^2 + (\frac{2}{3})x^3
$$

$$
y = -4x + \frac{2}{3}x^3
$$

注意：图片中的方程 $y = -2.67x + 0.67x^3$ 是基于错误的系数得出的。我们将使用我们计算出的正确方程进行验证。

验证:
我们需要检查这个方程是否通过所有给定的四个点。

验证 $P_1(-1, 2)$:

$$
y(-1) = -4(-1) + \frac{2}{3}(-1)^3 = 4 - \frac{2}{3} = \frac{12-2}{3} = \frac{10}{3}
eq 2
$$

哎呀，这里出现了问题。这说明我在计算 $C^{-1}$ 时可能出错了，或者题目本身的数据点并不在同一个三次多项式上。让我们重新检查一下。

让我们重新审视原始方程组，并利用 $a=0$ 来简化：
从 $P_2(0,0)$ 得到 $a=0$。
方程组变为：

$2 = -b + c - d$
$-2 = b + c + d$
$0 = 2b + 4c + 8d$

将方程1和方程2相加：
$(2) + (-2) = (-b+c-d) + (b+c+d)$
$0 = 2c \implies c = 0$

现在我们知道 $a=0, c=0$。方程组进一步简化为：

$2 = -b - d$
$-2 = b + d$ (这和第一个方程是一样的)
$0 = 2b + 8d$

从方程3: $2b + 8d = 0 \implies b + 4d = 0 \implies b = -4d$
代入方程1: $2 = -(-4d) - d = 4d - d = 3d \implies d = \frac{2}{3}$
然后 $b = -4(\frac{2}{3}) = -\frac{8}{3}$

所以，正确的系数应该是：
$a = 0$
$b = -\frac{8}{3} \approx -2.67$
$c = 0$
$d = \frac{2}{3} \approx 0.67$

这和图片中的结果一致！我之前的矩阵求逆计算有误。非常抱歉。