4. NURBS曲线核心：有理B样条的定义、权因子的影响、NURBS曲线的插值与拟合

好，咱们进入第四章。这一章要啃的，是NURBS——非均匀有理B样条。说实话，我刚入行那会儿，看到“有理”两个字就头大。后来在项目里被逼着用，才慢慢摸透了它的脾气。说白了，NURBS就是给B样条曲线加了个“权重”调节器，让咱们能精确描述圆锥曲线这类形状。

4.1 有理B样条的定义

先看定义。有理B样条曲线，其实就是把B样条的基函数做了一次“有理化”处理。公式长这样：

C(u) = Σᵢ₌₀ⁿ N_{i,p}(u) · wᵢ · Pᵢ  /  Σⱼ₌₀ⁿ N_{j,p}(u) · wⱼ

其中：

• N_{i,p}(u) —— p次B样条基函数
• wᵢ —— 控制点Pᵢ的权因子（权重）
• Pᵢ —— 控制点坐标
• 分母是所有权重基函数的和，保证曲线是“有理”的

你可能会问：为什么非要加个分母？我当初也这么想。后来在做一个汽车A柱截面时，用普通B样条死活拟合不出标准的圆弧，误差总是差那么一丁点。换成NURBS，把权重调一调，圆弧就精确了。原因就在于分母引入了有理分式，让曲线能表达圆锥曲线。

4.2 权因子的影响

权因子wᵢ，是NURBS的灵魂。我习惯把它想象成控制点的“质量”——质量越大，引力越强，曲线就被拉得越近。

具体影响规律：

• wᵢ → ∞：曲线无限逼近控制点Pᵢ，甚至穿过它
• wᵢ = 0：该控制点完全失效，曲线忽略它
• wᵢ = 1：退化为普通B样条（非有理）
• wᵢ < 0：曲线被“推开”，产生凹陷效果（慎用，容易出奇点）

我曾经在船体曲面设计里踩过一个坑：为了把曲线拉得更贴近某个型值点，我把权重调到了10以上。结果曲线在局部出现了“鼓包”，曲率突变，加工时刀具直接撞上了。后来学乖了——权重一般控制在0.5~5之间，超过这个范围就要检查曲线质量。

4.3 NURBS曲线的插值

插值，就是让曲线精确通过一组给定的型值点。这在逆向工程里特别常见——你扫描了一堆点云，要生成一条穿过这些点的曲线。

步骤我总结为三步：

1. 参数化：给每个型值点分配一个参数值u（常用弦长参数化或向心参数化）
2. 节点向量生成：根据参数值生成节点向量（一般用平均值法）
3. 求解控制点：建立线性方程组，解出控制点坐标和权重

方程组长这样：

对于每个型值点Qₖ：
Qₖ = C(uₖ) = Σᵢ N_{i,p}(uₖ) · wᵢ · Pᵢ  /  Σⱼ N_{j,p}(uₖ) · wⱼ

整理成矩阵形式：
[ A ] · [ P ] = [ Q ]

这里A矩阵是基函数矩阵，P是未知控制点，Q是已知型值点。解这个方程组，就能得到控制点坐标。

4.4 NURBS曲线的拟合

拟合和插值不同。拟合不要求曲线精确通过每个点，而是让曲线在整体上“最接近”这些点。适合处理带噪声的测量数据。

拟合的核心是最小二乘法：

目标函数：min Σₖ || Qₖ - C(uₖ) ||²

加上正则项（可选）：min Σₖ || Qₖ - C(uₖ) ||² + λ · ∫ ||C''(u)||² du

λ是光滑因子。λ越大，曲线越光滑，但偏离数据点越多。λ越小，曲线越贴近数据，但可能振荡。

我记得有一次做风机叶片截面拟合，原始点云有10万个点。直接用插值，曲线振荡得像心电图。后来改用拟合，λ取0.01，曲线既光滑又贴合，误差控制在0.02mm以内。嗯，这就是拟合的魅力。

方法	适用场景	优点	缺点
插值	精确通过型值点（如模具型面）	位置精度高	对噪声敏感，易振荡
拟合	带噪声的测量数据（如点云）	抗噪性强，曲线光滑	不精确通过每个点

4.5 知识结构图

下面这张图，把NURBS的核心脉络理清楚了。你一看就明白：权重怎么影响曲线，插值和拟合分别走哪条路。

4.6 实战建议

最后，给几个我自己的经验：

• 权重初始化：所有控制点先设w=1，跑一遍B样条，再根据偏差局部调整权重
• 插值前先检查数据：如果型值点间距差异过大，先做一次均匀化重采样，否则曲线会扭曲
• 拟合时别贪心：误差目标设得太小（比如0.001mm），曲线会过拟合，加工时反而出问题。留0.01~0.05mm的余量，更实用
• 曲率梳是照妖镜：每次调完权重或控制点，立刻看曲率梳。出现尖刺或突变，说明参数不合理