2. 贝塞尔曲线深入：德卡斯特里奥算法、升阶与降阶、拼接
贝塞尔曲线，说白了就是自由曲面设计的“地基”。上一章我们聊了它的基本概念和数学表达，但真正要在工程里用起来，光知道公式可不够。这一章，我们深入三个核心操作：怎么高效算点、怎么改变曲线阶数、怎么把多条曲线连起来。
我个人习惯把这三件事称为“贝塞尔三件套”。你在项目里遇到的 90% 的曲线问题，最后都能归到这三类里。不信？咱们往下看。
2.1 德卡斯特里奥算法：几何直觉与数值稳定性
先问个问题：给定控制点，怎么求曲线上 t=0.3 处的点？
你可能会说，套伯恩斯坦多项式公式呗。没错，公式能算，但数值上不太稳定。尤其是高阶曲线，多项式求值容易产生震荡。我在做汽车A柱曲面光顺时，就吃过这个亏——用公式算出来的点，在曲率变化剧烈的地方居然跳动了几个微米，这在模具加工里是不能接受的。
所以，我更推荐德卡斯特里奥算法。它本质上是一个递归线性插值过程，几何意义非常直观。
算法核心思想：
- 给定 n+1 个控制点 P₀, P₁, ..., Pₙ
- 对参数 t，在相邻控制点连线上按比例 t:1-t 取点，得到 n 个新点
- 重复这个过程，直到只剩一个点
- 这个最后的点，就是曲线上 t 对应的位置
举个例子，三次贝塞尔曲线（4个控制点）的递推过程：
// 德卡斯特里奥算法 - 三次贝塞尔曲线求值
// 输入：控制点 P[0..3]，参数 t
// 输出：曲线上点 Q
P0, P1, P2, P3 // 原始控制点
// 第一层插值
P01 = (1-t)*P0 + t*P1
P11 = (1-t)*P1 + t*P2
P21 = (1-t)*P2 + t*P3
// 第二层插值
P02 = (1-t)*P01 + t*P11
P12 = (1-t)*P11 + t*P21
// 第三层插值（最终点）
Q = (1-t)*P02 + t*P12你看，整个过程就是反复做“分点”。我经常跟团队说，这算法就像切蛋糕——一刀下去分成两块，再切一刀，直到切出你要的那一小块。
数值稳定性对比：
| 方法 | 10阶曲线误差 | 20阶曲线误差 | 适用场景 |
|---|---|---|---|
| 伯恩斯坦多项式直接求值 | ~10⁻¹² | ~10⁻⁶ | 低阶曲线（n≤10） |
| 德卡斯特里奥算法 | ~10⁻¹⁵ | ~10⁻¹⁴ | 任意阶曲线 |
德卡斯特里奥算法还有一个隐藏优势：它同时给出了曲线的切向量。还记得上面例子里的 P₀₂ 和 P₁₂ 吗？它们的连线方向，就是曲线在 Q 点处的切线方向。这个特性在做曲线光顺分析时特别有用，省去了求导的麻烦。
我的习惯：在代码里实现贝塞尔曲线求值时，永远用德卡斯特里奥算法。它比公式法慢一点点，但数值稳定性好太多。对于自由曲面设计这种对精度要求极高的场景，这点性能开销完全值得。
2.2 贝塞尔曲线的升阶与降阶
升阶和降阶，听起来像是数学游戏，但在实际工程里非常实用。
升阶：在不改变曲线形状的前提下，增加控制点的数量。为什么要这么做？因为高阶曲线能表达更复杂的形状，而且多个低阶曲线拼接时，升阶到同一阶数可以简化连续性约束。
降阶：反过来，减少控制点数量。这通常是为了数据压缩或简化计算。但要注意，降阶会引入误差——你不可能用更少的自由度完美表达原来的形状。
升阶公式其实很优雅。假设 n 阶贝塞尔曲线有控制点 P₀, P₁, ..., Pₙ，升阶到 n+1 阶后，新控制点 Q₀, Q₁, ..., Qₙ₊₁ 为：
Q₀ = P₀
Qₙ₊₁ = Pₙ
Qᵢ = (i/(n+1)) * Pᵢ₋₁ + (1 - i/(n+1)) * Pᵢ 对于 i = 1, 2, ..., n说白了，新控制点就是旧控制点的加权平均。这个公式我在做车身曲面光顺时用过无数次——把不同阶数的曲线统一到同一阶数，然后做连续性拼接。
降阶的坑：我曾经在一个航空叶片项目中，为了减少数据量，把一条7阶曲线降阶到5阶。结果在叶片前缘曲率最大的地方，形状偏差达到了0.05mm。虽然看起来不大，但在气动分析里，这个偏差直接导致升力系数下降了2%。从那以后，我降阶时一定会做误差分析，并且只在曲率变化平缓的区域使用。
降阶的实用方法：
- 最小二乘逼近：用低阶曲线拟合高阶曲线上的采样点
- 端点约束法：保持端点位置和切向不变，优化内部点
- 迭代降阶：每次只降一阶，逐步逼近
我个人推荐第三种方法。每次只降一阶，误差可控，而且可以随时停下来检查。
2.3 贝塞尔曲线的拼接
单条贝塞尔曲线能表达的形状有限。要设计复杂的自由曲面，比如汽车车身、飞机机翼，必须把多条曲线首尾相连。这就引出了拼接的连续性条件。
三种连续性等级：
| 连续性 | 含义 | 数学条件 | 工程应用 |
|---|---|---|---|
| C⁰ | 位置连续 | Pₙ = Q₀ | 简单连接，有尖角 |
| C¹ | 切向连续 | Pₙ = Q₀ 且 Pₙ - Pₙ₋₁ = k(Q₁ - Q₀) | 光滑连接，常见于内饰件 |
| C² | 曲率连续 | C¹条件 + 二阶导连续 | 高光反射连续，用于外覆盖件 |
你想想看，C⁰ 就是两段线端点碰上了，但方向可能突变。C¹ 是方向也一致了，看起来光滑。C² 是曲率也连续了，高光反射不会出现“断带”。
我在做汽车侧围外板时，要求所有拼接点必须达到 C² 连续。为什么？因为车身曲面在喷漆后，任何曲率不连续都会在高光下暴露出来，形成肉眼可见的“波浪”。
C¹ 拼接的具体实现：
// 两条三次贝塞尔曲线 C¹ 拼接
// 曲线1: P0, P1, P2, P3
// 曲线2: Q0, Q1, Q2, Q3
// C⁰ 条件：P3 = Q0
Q0 = P3
// C¹ 条件：P3 - P2 = Q1 - Q0（假设k=1，即参数化一致）
Q1 = 2*P3 - P2你看，C¹ 拼接其实很简单——只要把 Q₁ 放在 P₂ 关于 P₃ 的对称点上就行。C² 拼接稍微复杂一点，需要同时约束 Q₂ 的位置。
拼接的实用建议：
- 尽量使用相同阶数的曲线拼接，避免升阶降阶带来的误差
- 拼接点附近的控制点分布要均匀，避免“扎堆”导致数值不稳定
- 对于高光反射面，至少要求 C² 连续；对于结构面，C¹ 通常就够了
2.4 本章知识体系
下面这张图总结了贝塞尔曲线深入操作的三个核心模块及其关系：
我的学习路径建议：先掌握德卡斯特里奥算法，它是所有操作的基础。然后练习升阶降阶，理解控制点与曲线形状的关系。最后再学拼接，这时候你会发现，前面两个技能正好为拼接铺平了道路。