4. NURBS曲线求导：一阶导数、二阶导数、曲率计算
好，咱们进入正题。NURBS曲线求导，说白了就是搞清楚曲线在某个点到底是怎么走的。你想想看，刀具沿着曲线跑，它不光要知道往哪走，还得知道怎么拐弯、拐多急。这就是一阶导、二阶导和曲率要干的事。
我在项目里调试五轴机床时，就遇到过因为曲率计算不准，导致刀具在拐弯处剧烈抖动的情况。那次折腾了我整整三天，最后发现是求导公式里一个分母项写错了。嗯，从那以后我对求导这块就格外上心。
4.1 一阶导数：速度与方向
一阶导数，其实就是曲线在参数u处的速度向量。它告诉刀具：你现在该往哪个方向走，走多快。
对于NURBS曲线，一阶导数的公式长这样：
C'(u) = [A'(u) - w'(u) * C(u)] / w(u)其中A(u)是加权控制点的和，w(u)是权因子函数。说白了，就是分子分母分别求导，再做个除法修正。
我个人习惯把一阶导数拆成两部分来理解：
- 方向分量：由控制点位置决定，告诉刀具往哪切
- 速度分量：由权因子变化率决定，影响走刀的快慢
重要提醒：一阶导数不能为零。如果C'(u)=0，说明曲线在该点出现了"停顿"，这在加工中是要避免的。我见过有人把控制点排得太密，结果导数归零，刀具直接卡在那不动了。
4.2 二阶导数：加速度与弯曲趋势
二阶导数描述的是曲线弯曲的趋势。你可以把它想象成刀具的加速度向量——它告诉你刀具下一步是加速拐弯还是减速直行。
公式稍微复杂一点：
C''(u) = [A''(u) - 2*w'(u)*C'(u) - w''(u)*C(u)] / w(u)为什么会这样？因为NURBS是有理形式，分母不能直接消掉。你想想看，如果分母是常数（比如B样条），那公式就简单多了。但NURBS多了权因子这个"调节器"，求导就得小心处理。
我的经验：实际编程时，我建议先算w(u)、w'(u)、w''(u)这三个标量，再算A(u)、A'(u)、A''(u)这三个向量。最后组合起来，不容易出错。我曾经图省事直接套公式，结果调试了整整一个下午才发现符号写反了。
4.3 曲率计算：拐弯有多急
曲率，这才是咱们做高速高精插补最关心的东西。它直接决定了刀具在某个点能跑多快而不产生过切。
曲率公式：
κ(u) = |C'(u) × C''(u)| / |C'(u)|³注意，这里是叉积的模除以一阶导数模的三次方。为什么是三次方？嗯，你可以这么理解：一阶导决定速度，二阶导决定加速度，曲率则是这两者的"配合结果"。
| 曲率值 | 物理含义 | 加工建议 |
|---|---|---|
| κ = 0 | 直线段 | 可以全速前进 |
| κ 较小 | 缓弯 | 适当减速，注意轮廓误差 |
| κ 较大 | 急弯 | 必须大幅减速，否则过切 |
| κ → ∞ | 尖点 | 需要特殊处理，通常要打断曲线 |
避坑指南：我曾经在加工一个模具的尖角时，曲率计算值突然变得极大。后来发现是因为控制点权因子设置不合理，导致曲线在数学上出现了"几乎尖点"。从那以后，我每次都会检查曲率的连续性——如果曲率突变超过50%，我就会重新调整控制点布局。
4.4 实用计算流程
在实际的数控系统中，我们不会每次都从头算起。我一般按这个流程走：
- 预处理阶段：把NURBS的节点向量、控制点、权因子全部读入内存
- 基函数求导：用de Boor-Cox算法递推计算基函数及其导数
- 组合计算：按公式组装一阶导、二阶导
- 曲率输出：用叉积公式算出曲率值
这里有个小技巧：基函数求导可以用递推公式，不用每次都从头算。我习惯把基函数值存成一个二维数组，这样后面算导数时直接查表，能省不少时间。
4.5 代码实现要点
给你看一段核心代码的思路：
// 计算NURBS曲线在参数u处的一阶导、二阶导和曲率
void NURBS_CurveDerivatives(double u,
Vector3d& C, // 曲线点
Vector3d& dC, // 一阶导
Vector3d& ddC, // 二阶导
double& kappa) // 曲率
{
// 1. 计算基函数及其导数
double N[degree+1], dN[degree+1], ddN[degree+1];
BasisFunctionDerivatives(u, knot_vector, degree, N, dN, ddN);
// 2. 计算加权控制点组合
Vector3d A(0,0,0), dA(0,0,0), ddA(0,0,0);
double w = 0, dw = 0, ddw = 0;
for(int i = 0; i <= degree; i++)
{
int idx = span - degree + i;
A += N[i] * w_control_points[idx];
dA += dN[i] * w_control_points[idx];
ddA += ddN[i] * w_control_points[idx];
w += N[i] * weights[idx];
dw += dN[i] * weights[idx];
ddw += ddN[i] * weights[idx];
}
// 3. 计算曲线点
C = A / w;
// 4. 计算一阶导
dC = (dA - dw * C) / w;
// 5. 计算二阶导
ddC = (ddA - 2*dw*dC - ddw*C) / w;
// 6. 计算曲率
double cross_norm = dC.cross(ddC).norm();
double dC_norm = dC.norm();
kappa = cross_norm / (dC_norm * dC_norm * dC_norm);
}个人建议：实际工程中，我建议把曲率计算单独封装成一个函数。因为插补过程中，每个插补周期都要算一次曲率，如果每次都重新算一阶导、二阶导，计算量太大了。我一般会缓存最近几个点的导数结果，用差分法快速估算曲率变化趋势。
4.6 知识体系总览
下面这张图是我自己整理的NURBS求导知识结构，你看一眼就能明白各个概念之间的关系：
从这张图你能看出来，一阶导和二阶导是曲率的基础，而曲率又直接决定了速度规划策略。我在实际项目中，就是按照这个结构来组织代码的——先保证导数算得准，再谈曲率控制。
好了，关于NURBS求导的内容就这些。记住，求导不是目的，目的是让刀具跑得又快又准。下次你调试机床时，如果发现某个拐弯处有异常抖动，不妨先检查一下曲率计算——八成是这里出了问题。