2. 插值问题的数学描述

在上一节中，我们已经建立了“用简单函数穿过已知离散点”这一直观目标。要把这个目标变成可以推导、可以计算、可以评价的数学问题，必须把“数据点”“插值条件”“解的存在唯一性”“误差”和“光滑性”这五个要素严格地写清楚。本节的任务就是完成这一形式化过程，并说明每个要素在后续样条构造中扮演什么角色。

一、给定数据点：插值问题的输入

插值问题的原始输入是一组离散数据。设区间 [a, b] 上给定 n+1 个互异节点

a = x₀ < x₁ < … < x_n = b

以及对应的函数值 y₀, y₁, …, y_n。这组有序对 (x_i, y_i) 就是插值数据。需要强调三点：第一，节点必须互异，否则同一个 x 对应两个 y，插值函数无法成为单值函数；第二，节点顺序通常按升序排列，这为分段构造和误差估计提供便利；第三，数据可能来自实测、实验采样或某个已知函数的离散取值，来源不同会影响我们对误差的理解，但不改变插值问题本身的数学结构。

在实际工程中，数据点往往还带有精度信息。例如某温度传感器在 0、10、20、30、40 秒采样得到 20.1、22.4、25.0、27.9、30.5 ℃，若传感器分辨率为 0.1 ℃，则插值曲线在节点处的偏差不应被解释为模型误差，而应视为数据本身的噪声。这一点在后续讨论误差时非常重要。

二、插值条件：函数类与约束

仅有数据点还不足以定义问题，还必须指定“用什么函数去穿过这些点”。一般地，选取一个函数空间 S（例如次数不超过 n 的多项式空间 P_n，或三次样条空间），要求找到 s ∈ S 满足

s(x_i) = y_i, i = 0, 1, …, n

这组等式称为插值条件。它是 n+1 个等式约束，恰好对应 n+1 个数据点。插值条件的形式看似简单，但它隐含了两个关键选择：其一是函数空间 S 的维数，其二是约束的个数是否与自由度匹配。若 S 的维数大于 n+1，则满足条件的函数不唯一，需要追加额外条件（如样条中的端点条件）；若维数小于 n+1，则一般无解。因此，“插值条件”与“函数空间”必须成对出现，缺一不可。

三、存在唯一性：解是否可靠

多项式插值有一个经典结论：对 n+1 个互异节点，次数不超过 n 的多项式插值函数存在且唯一。证明思路是构造范德蒙德矩阵 V，其元素为 x_i^j，插值条件化为线性方程组 Vc = y。由于节点互异，V 的行列式非零，故系数向量 c 唯一确定。这个结论给出了插值问题“良定”的数学保证。

但存在唯一性并不等于数值可靠。当 n 较大时，范德蒙德矩阵条件数急剧增大，直接求解会带来严重舍入误差。这正是实际计算中不直接解范德蒙德系统、而采用拉格朗日形式、牛顿差商形式或分段低次插值的原因。对于样条插值，存在唯一性依赖于边界条件的选取：三次样条在给定两端一阶或二阶导数条件后，线性系统是三对角严格对角占优的，解存在且唯一，且可以用追赶法稳定求解。

四、误差与光滑性：评价插值质量的两把尺子

插值函数 s 与被插值函数 f 之间的偏差用误差刻画。多项式插值的经典误差公式为

f(x) − p_n(x) = f⁽ⁿ⁺¹⁾(ξ) / (n+1)! · ∏_i=0ⁿ(x − x_i)

其中 ξ 位于包含 x 与所有节点的区间内。该式说明误差由两部分决定：一是 f 的高阶导数，反映被插函数本身的“弯曲程度”；二是节点多项式 ω(x) = ∏(x − x_i) 的绝对值，反映节点分布。这解释了两个现象：高次多项式在区间边缘容易出现龙格振荡，因为 ω(x) 在端点附近偏大；而等距节点并非最优，切比雪夫节点能压低 ω(x) 的最大值。

光滑性要求则是另一把尺子。它约束插值函数在节点处的连续性阶数。分段线性插值只有 C⁰ 连续，节点处出现折角；三次样条要求 C² 连续，即函数、一阶导数、二阶导数都连续，因而视觉上“顺滑”，且二阶导数连续意味着曲率连续，适合描述机械运动轨迹和曲面轮廓。光滑性越高，所需附加条件越多，线性系统的规模与耦合也越强，这是光滑性与计算代价之间的基本权衡。

五、把五要素落到一个可计算的例子上

设节点为 x = 0, 1, 2, 3，函数值为 y = 1, 2, 0, 3。若取 S = P₃，则插值条件给出四个方程，三次多项式 p(x) = a₀ + a₁x + a₂x² + a₃x³ 的四个系数由线性系统唯一确定。用牛顿差商法可依次算出零阶差商 1、2、0、3，一阶差商 1、−2、3，二阶差商 −1.5、2.5，三阶差商 4/3，于是

p(x) = 1 + 1·x − 1.5·x(x−1) + (4/3)·x(x−1)(x−2)

展开后即可在任意 x 处求值。若改用三次样条，则除四个插值条件外，还需两个端点条件，例如自然边界 s''(0) = s''(3) = 0，此时六个条件确定四个区间上的三次多项式共 16 个系数中的独立部分，通过 C² 连续性方程联立求解。两条路线都满足“过点”，但多项式插值在区间外可能剧烈发散，而样条在节点间更平稳，这正是误差与光滑性要求共同作用的结果。