2. 数学基础：贝塞尔曲线与曲面的原理与参数化表示
做自由曲面设计，绕不开贝塞尔。说实话，我刚入行那会儿，觉得这东西就是一堆数学公式，离实际工程很远。直到有一次，我负责一个汽车A级曲面的项目，用传统方法怎么调都顺不了，最后回头啃透了贝塞尔的基础，才把问题搞定。嗯，从那以后，我再也不敢小看这些“基础数学”了。
2.1 从一条曲线说起：贝塞尔曲线的直觉
贝塞尔曲线，说白了就是用几个控制点来“拽”出一条光滑曲线。你想想看，你给四个点，曲线并不穿过中间两个，而是被它们“吸引”过去。这种设计方式，特别适合我们做造型——你不需要直接操作曲线上的每个点，而是通过调整控制点来改变整体形状。
核心思想： 贝塞尔曲线是控制点的加权平均。权重由 Bernstein 基函数决定。
数学上，一条 n 次贝塞尔曲线定义为：
C(t) = Σᵢ₌₀ⁿ Pᵢ · Bᵢ,ₙ(t), t ∈ [0, 1]其中 Pᵢ 是控制点，Bᵢ,ₙ(t) 是 Bernstein 基函数：
Bᵢ,ₙ(t) = C(n, i) · tⁱ · (1-t)ⁿ⁻ⁱ这里 C(n, i) 是组合数。我个人习惯把 t 想象成“时间”——从 0 到 1，曲线从起点走到终点。
2.2 三次贝塞尔：工程中最常用的形式
为什么是三次？因为二次曲线只能画弧线，自由度不够；高次曲线计算量大，而且容易抖动。我在项目中用过五次贝塞尔，结果控制点稍微一动，曲线就“抽风”，调试起来特别痛苦。三次贝塞尔刚好平衡了灵活性和稳定性。
三次贝塞尔有 4 个控制点（P₀, P₁, P₂, P₃），展开公式如下：
C(t) = (1-t)³·P₀ + 3(1-t)²t·P₁ + 3(1-t)t²·P₂ + t³·P₃你看，当 t=0 时，只有 P₀ 起作用；t=1 时，只有 P₃ 起作用。中间两个点 P₁ 和 P₂ 决定了曲线的“走势”。
| 参数 | 含义 | 工程影响 |
|---|---|---|
| P₀, P₃ | 曲线端点 | 决定曲线起始和结束位置 |
| P₁ | 起点切向控制点 | P₀→P₁ 方向决定起点切线方向 |
| P₂ | 终点切向控制点 | P₂→P₃ 方向决定终点切线方向 |
| |P₀P₁|, |P₂P₃| | 切向长度 | 长度越大，曲线越“冲”向该方向 |
我的经验： 做曲面拼接时，先保证边界曲线的 G¹ 连续（切线方向一致），再调 G²（曲率连续）。G² 连续对视觉质量影响很大，但计算量也大。我一般先用 G¹ 快速搭骨架，最后再精修 G²。
2.3 从曲线到曲面：张量积贝塞尔曲面
曲面怎么来？很简单，把曲线“拉伸”成面。贝塞尔曲面是两条贝塞尔曲线的“乘积”——一个方向用 u 参数，另一个方向用 v 参数。
数学定义：
S(u, v) = Σᵢ₌₀ᵐ Σⱼ₌₀ⁿ Pᵢ,ⱼ · Bᵢ,ₘ(u) · Bⱼ,ₙ(v), u,v ∈ [0, 1]其中 Pᵢ,ⱼ 是一个 (m+1)×(n+1) 的控制点网格。我刚开始学的时候，觉得这个公式很抽象。后来我把它理解成“先沿 u 方向算曲线，再沿 v 方向算曲线”——说白了就是两次一维贝塞尔。
注意： 控制点网格必须是拓扑矩形的。也就是说，行数和列数都是固定的。你不能在中间突然多一个点或少一个点。这是贝塞尔曲面的一个硬约束，也是它的局限性。
2.4 参数化表示：为什么是 [0, 1]？
你可能注意到了，t、u、v 的取值范围都是 [0, 1]。为什么不是 [0, 100] 或者 [-1, 1]？
我个人觉得，[0, 1] 的好处是归一化。不管你的曲面有多大，参数空间都是固定的。这样在做曲面拼接、求交、渲染时，计算起来特别方便。你想想看，如果每个曲面都用不同的参数范围，那拼接时还得做参数变换，多麻烦。
另外，[0, 1] 区间也让 Bernstein 基函数有很好的数值稳定性。我曾经试过用 [0, 100] 的参数范围，结果高次项的计算误差大得离谱，曲面都变形了。
2.5 贝塞尔曲面的核心性质
做自由曲面设计，必须掌握这几个性质。它们直接决定了你的曲面好不好调、能不能拼：
- 端点插值： 曲面的四个角点就是控制点网格的四个角。你移动角点，曲面角点跟着走。
- 凸包性： 整个曲面都在控制点网格的凸包内。这意味着你只要控制好控制点的范围，曲面就不会“飞出去”。
- 变差缩减性： 曲面不会比控制点网格“更弯”。这个性质在造型时很有用——你想让曲面更平缓，就把控制点排得更平。
- 仿射不变性： 对曲面做旋转、平移、缩放，等价于对控制点做同样的变换。这个性质让我们可以在世界坐标系里直接操作控制点。
避坑指南： 我曾经在做一个大曲率过渡面时，为了追求曲面“更饱满”，把中间控制点拉得很远。结果曲面虽然穿过了目标点，但出现了明显的“鼓包”。后来发现是违反了凸包性的约束——控制点太远，曲面虽然还在凸包内，但曲率分布已经不均匀了。
2.6 知识结构总览
下面这张图，是我梳理的本章知识体系。你可以把它当作一个“地图”，随时回来对照：
2.7 工程中的参数化技巧
在实际项目中，参数化不只是数学问题。我总结了几条实用经验：
- 参数均匀化： 如果控制点分布不均匀，参数 t 的等步长采样会导致曲面上点分布也不均匀。我一般会在预处理阶段做一次弦长参数化，让采样点更均匀。
- 避免退化： 不要让相邻控制点重合或共线。否则曲面会出现“尖点”或“褶皱”。我曾经吃过这个亏——一个控制点网格有两行几乎重合，结果曲面在那一带出现了明显的“折痕”。
- 边界处理： 曲面边界就是控制点网格的边界。如果你需要和另一个曲面拼接，先对齐边界上的控制点，再调内部点。
- 数值稳定性： 高次贝塞尔（n > 10）的 Bernstein 基函数计算容易溢出。我建议用 de Casteljau 算法递推计算，既稳定又直观。
小技巧： de Casteljau 算法不仅可以求曲线上点的位置，还能同时求出切线方向。你只需要在递推过程中保留中间点的差值即可。这个技巧在做曲面法向量计算时特别有用。
2.8 贝塞尔曲面的局限与应对
贝塞尔曲面不是万能的。它的最大问题是：全局控制。你移动一个控制点，整个曲面都会受影响。这在做局部微调时非常痛苦。
怎么解决？两个方向：
- 分段贝塞尔： 把大曲面切成多个小曲面片，每个片用低次贝塞尔。片与片之间保证连续性。这就是 B 样条曲面的雏形。
- 升阶与降阶： 如果当前次数不够灵活，可以升阶（增加控制点）而不改变曲线形状。反过来，如果次数太高想简化，可以降阶——但会有误差，需要控制误差范围。
嗯，说到 B 样条，那是下一阶段的内容了。但你要记住，贝塞尔是基础。你把它吃透了，后面学 NURBS、细分曲面都会轻松很多。