第三节：三次样条的定义——分段多项式、连续性与光滑性
好，咱们今天聊聊三次样条的定义。说实话，我刚入行那会儿，看到“样条”这俩字，第一反应是“这玩意儿是不是跟弹簧有关？”后来做多了才知道，它确实有点像——一根有弹性的细木条，穿过几个点，自然弯曲成一条光滑的曲线。嗯，这就是样条名字的由来。
三次样条，说白了就是每两个数据点之间用一个三次多项式来连接。你想想看，如果只用一条直线去连，那肯定又硬又丑；如果用高次多项式去拟合所有点，又容易“龙格现象”——两头翘得飞起。三次样条刚好是个折中，既光滑又稳定。
核心定义：给定 n+1 个数据点 (x₀, y₀), (x₁, y₁), ..., (xₙ, yₙ)，三次样条 S(x) 满足：
- 在每个子区间 [xᵢ, xᵢ₊₁] 上，S(x) 是一个三次多项式
- S(x) 在整个区间 [x₀, xₙ] 上具有二阶连续导数
- S(xᵢ) = yᵢ，即通过所有给定点
分段多项式表示
我个人习惯把三次样条想象成“拼接的艺术”。每个小区间上都是一个三次多项式，但相邻区间之间要“无缝衔接”。
假设我们把区间 [xᵢ, xᵢ₊₁] 上的三次多项式写成：
Sᵢ(x) = aᵢ + bᵢ(x - xᵢ) + cᵢ(x - xᵢ)² + dᵢ(x - xᵢ)³这里 i = 0, 1, ..., n-1。你看，每个区间有 4 个未知系数（aᵢ, bᵢ, cᵢ, dᵢ），总共 n 个区间，那就是 4n 个未知数。要解出这些系数，我们需要 4n 个方程。
我在项目中遇到过一个问题：有人直接用高次多项式去拟合几十个点，结果曲线在端点处剧烈震荡，加工出来的零件表面全是波纹。后来换成三次样条，问题立马解决。说白了，分段的好处就是“局部调整，不影响全局”。
连续性条件
连续性条件，就是保证曲线在连接点处“不断开”。具体来说：
- C⁰ 连续：函数值连续。即 Sᵢ(xᵢ₊₁) = Sᵢ₊₁(xᵢ₊₁) = yᵢ₊₁
- C¹ 连续：一阶导数连续。即 Sᵢ'(xᵢ₊₁) = Sᵢ₊₁'(xᵢ₊₁)
- C² 连续：二阶导数连续。即 Sᵢ''(xᵢ₊₁) = Sᵢ₊₁''(xᵢ₊₁)
为什么会要求到二阶导数连续？你想想看，如果只有一阶连续，曲线虽然看起来没尖角，但曲率变化可能很突兀。对于数控加工来说，刀具路径的曲率突变会导致加速度突变，轻则震动，重则撞刀。我曾经就因为这个吃过亏——当时只保证了 C¹ 连续，结果加工时机床抖得跟筛子似的。
避坑指南：我曾经在调试五轴机床时发现，如果只保证 C¹ 连续，刀具在连接点处的加速度会跳变，导致表面有刀痕。后来强制加上 C² 连续条件，问题才解决。所以，三次样条的 C² 连续不是数学上的“炫技”，而是工程上的刚需。
光滑性条件
光滑性条件，其实就是连续性条件的延伸。三次样条之所以叫“三次”，是因为它最高能保证 C² 连续。如果要求 C³ 连续，那就得用四次样条了，但工程上很少用——计算量太大，而且容易过拟合。
我们来看一下具体的方程数量：
| 条件类型 | 方程数量 | 说明 |
|---|---|---|
| 插值条件 | n+1 | S(xᵢ) = yᵢ，每个数据点一个方程 |
| C⁰ 连续 | n-1 | 内部节点处函数值相等 |
| C¹ 连续 | n-1 | 内部节点处一阶导数相等 |
| C² 连续 | n-1 | 内部节点处二阶导数相等 |
| 边界条件 | 2 | 端点处额外约束 |
| 总计 | 4n | 正好等于未知数个数 |
你看，4n 个方程，4n 个未知数，刚好可解。边界条件通常有两种：自然边界（二阶导数为零）和固定边界（指定一阶导数）。我个人习惯用自然边界，因为它在大多数情况下表现最稳定。
注意：边界条件的选择会影响整条曲线的形状。自然边界适合“两端自由”的场景，比如自由曲线设计；固定边界适合“两端有约束”的场景，比如刀具路径的切入切出角度控制。选错了边界条件，曲线可能会在端点处“翘尾巴”。
知识体系结构图
下面这张图，是我自己总结的三次样条知识体系。你看，从定义出发，分支出分段表示、连续性、光滑性三个核心概念，最后汇聚到工程应用。嗯，这样学起来就清晰多了。
嗯，到这里，三次样条的定义就讲完了。说白了，它就是一套“分段三次多项式 + 连续光滑约束”的组合拳。你只要记住：每个区间一个三次式，相邻区间在连接点处函数值、一阶导、二阶导都相等，再加上两个边界条件，齐活。
下次你写样条插补代码的时候，不妨先画一画这个结构图，思路会清晰很多。我当年就是这么干的——先理清数学关系，再动手写代码，事半功倍。