二、NURBS基础回顾：NURBS曲线定义、控制点、节点向量、权因子等核心概念

在深入NURBS插补参数优化之前，必须对NURBS曲线的数学本质与几何构成建立清晰认知。NURBS（Non-Uniform Rational B-Spline，非均匀有理B样条）之所以成为CAD/CAM与数控加工中描述自由曲面的行业标准，是因为它能够用统一的数学形式精确表达解析曲线（如圆弧、椭圆）与自由曲线，同时具备局部修改性和凸包性等优良几何特性。本节将系统梳理NURBS曲线的定义、控制点、节点向量、权因子四大核心要素，并阐明它们在插补参数优化中的实际意义。

上图展示了NURBS曲线各核心要素之间的逻辑关系：控制点提供几何骨架，权因子调节局部“引力”，节点向量定义参数域划分，而B样条基函数则负责将上述要素组合为最终曲线。在插补参数优化中，这四个要素直接影响进给速度规划、弓高误差控制和加减速策略的制定。下面逐一展开分析。

2.1 NURBS曲线数学定义与有理表示

核心原理：NURBS曲线是B样条曲线的有理推广，其数学定义为：

C(u) = Σ(i=0→n) [N_i,p(u) · w_i · P_i] / Σ(i=0→n) [N_i,p(u) · w_i]

其中，P_i为控制点（二维或三维坐标），w_i为对应权因子（正实数），N_i,p(u)为定义在节点向量U上的p次B样条基函数。当所有权因子w_i=1时，NURBS退化为非有理B样条曲线；当节点向量采用特定分布时，可精确表示圆弧、椭圆等二次曲线。有理形式带来的关键优势在于：投影变换不变性——对曲线进行透视投影变换后，其数学形式保持不变，这为CAD/CAM中的视图变换提供了极大便利。

判断方法：在实际工程中，判断一条曲线是否为NURBS，可检查其数据文件是否包含控制点坐标、权因子数组和节点向量三个独立数组。STEP文件（ISO 10303-42）中，B_SPLINE_CURVE_WITH_KNOTS实体即对应NURBS定义；IGES文件类型126则直接存储有理B样条曲线参数。

2.2 控制点与凸包性：几何直觉与工程约束

核心原理：控制点构成曲线的“控制多边形”，曲线本身并不通过控制点（除首末点外），但始终位于控制点形成的凸包内。这一凸包性具有重要工程意义：它保证了插补轨迹不会出现超出控制点范围的剧烈波动，从而限制了加工路径的几何边界。控制点的局部修改性意味着移动一个控制点，仅影响曲线的一部分（由节点向量中对应支撑区间决定），这为CAM编程中的局部修形提供了理论依据。

操作流程：在数控编程中，读取NURBS刀轨时，应首先检查控制点序列的连续性。若相邻控制点间距突变超过阈值（如平均间距的5倍），则可能预示曲线存在尖角或自交风险，需在插补前进行几何验证。

// 控制点间距检查伪代码
for i = 1 to n:
    dist = |P_i - P_{i-1}|
    if dist > 5 * avg_dist:
        warning("控制点间距异常，检查曲线几何")
    if dist < 0.01 * avg_dist:
        warning("控制点重合，可能导致基函数奇异")

案例拆解：某航空叶片加工项目中，原始CAM刀轨包含120个控制点，其中第45和第46个控制点间距仅为平均值的1/50。插补时，该区域出现进给率剧烈波动，表面粗糙度超标。经分析，原因是控制点过密导致基函数数值计算出现病态。解决方案：在CAM软件中合并这两个控制点，并重新计算权因子，使曲线几何偏差控制在0.005mm以内，插补稳定性显著改善。

2.3 节点向量：参数域划分与局部支撑性

核心原理：节点向量U是一个非递减的实数序列，将参数区间[u_0, u_m]划分为若干子区间。节点向量的分布决定了B样条基函数的支撑区间，进而控制曲线在各参数段的形状。节点向量分为三种类型：均匀（节点等距分布）、准均匀（两端节点重复p+1次，内部均匀）、非均匀（节点间距不等）。在NURBS插补中，节点向量直接决定了参数u与弧长s之间的映射关系，这是进给速度规划的基础。

判断方法：检查节点向量中重复节点的数量。若内部节点重复度达到p，则曲线在该节点处可能只具有C⁰连续（出现尖角）；若重复度为p-1，则C¹连续。对于高速高精加工，应确保关键路径段至少C²连续，即内部节点重复度不超过p-2。

实战场景：在五轴加工中，刀尖点轨迹与刀轴矢量轨迹通常采用不同的节点向量。若两者节点向量不一致，插补时需进行同步处理。某模具加工案例中，刀尖轨迹节点向量包含38个节点，而刀轴轨迹包含52个节点，导致插补周期内刀轴角度突变。解决方案：采用节点插入算法（Boehm算法）将两组节点向量统一到并集，保证每个插补周期内刀尖与刀轴同步更新。

2.4 权因子：有理调节与几何变形

核心原理：权因子w_i是NURBS区别于普通B样条的关键参数。增大w_i会使曲线更靠近控制点P_i，减小则远离。当w_i趋于无穷时，曲线通过该控制点；当w_i=0时，该控制点完全失效。权因子为正值是保证曲线不出现奇异点的必要条件。在圆弧的NURBS表示中，权因子取cos(θ/2)（θ为圆弧半角），这是有理表示精确描述二次曲线的经典应用。

操作流程：在插补参数优化中，权因子影响曲线曲率分布，进而影响弓高误差。对于曲率较大的区域，应适当增大权因子以“拉紧”曲线，减少弓高误差；但过大的权因子会导致速度规划中的加速度突变。建议权因子取值范围为0.5~2.0，超出此范围需进行几何验证。

// 权因子影响评估
function evaluateWeightEffect(P_i, w_i, u_range):
    for u in u_range:
        C_old = computeNURBS(u, w_i)
        C_new = computeNURBS(u, w_i * 1.1)
        deviation = |C_new - C_old|
        if deviation > tolerance:
            warning("权因子变化导致几何偏差超限")

案例拆解：某汽车覆盖件模具的NURBS刀轨中，一个关键控制点的权因子被CAM软件设置为3.5，导致该区域曲率半径仅为周边平均值的1/8。插补时，弓高误差达到0.03mm，超出0.01mm的工艺要求。分析发现，权因子过大使得曲线在该处过度“内凹”。优化方案：将权因子从3.5降至2.0，同时微调相邻控制点位置，使最大弓高误差降至0.008mm，且进给速度波动减少40%。

2.5 基函数递推与数值稳定性

核心原理：B样条基函数通过Cox-de Boor递推公式定义：

N_i,0(u) = 1 若 u_i ≤ u < u_{i+1}，否则 0
N_i,p(u) = (u-u_i)/(u_{i+p}-u_i) · N_i,p-1(u) + (u_{i+p+1}-u)/(u_{i+p+1}-u_{i+1}) · N_{i+1,p-1}(u)

该递推公式在数值计算中需注意分母为零的情况（当节点重复时）。工程实现中，通常采用改进的算法，在分母小于阈值时直接置零，避免除零错误。基函数的局部支撑性（N_i,p(u)仅在[u_i, u_{i+p+1}]区间非零）是NURBS局部修改性的数学基础。

判断方法：评估基函数计算效率时，可统计每个参数点u处非零基函数的数量。对于p次曲线，最多有p+1个非零基函数。若实际计算中非零基函数数量超过p+1，说明节点向量存在异常重复，需检查数据。

实战场景：在实时插补器中，基函数计算占用了约30%的CPU时间。为优化性能，可采用预计算策略：将参数域离散为N个区间，预先计算每个区间端点的基函数值及导数，插补时通过线性插值近似。某商用数控系统采用此方法，将基函数计算时间从12μs降至3μs，同时保持弓高误差精度在0.1%以内。