1.一种多轴联动轮廓误差实时检测方法，其特征在于，操作步骤如下：(1)假设期望轨迹的参数方程为P(u)，当前位控周期的参考点为P(u0)，读取当前位控周期数控机床上各轴的实际点Pa；计算实际点Pa与期望轨迹P(u)上任一点之间的距离f(u)，计算所述距离f(u)关于参数u的导数f’(u)；

所述数控机床为内齿珩轮强力珩齿机床或五轴数控加工机床或三轴数控加机机床；

具体操作如下：

假设期望轨迹的参数方程为P(u)，轨迹对应的参数取值范围是[us,ue]，当前位控周期的参考点为P(u0)，读取当前位控周期数控机床上各轴的实际点Pa；计算实际点Pa与期望轨迹上任一点之间的距离f(u)，所述距离f(u)＝|P(u)‑Pa|；并且计算所述距离f(u)关于参数u的导数f’(u)；

(2)为了搜寻导数f’(u)的零点，首先以u＝u0为中心，往两边同时搜索f’(u)符号发生变化的参数值u，搜索步距为Δu，最终搜索到导数f’(u1)·f’(u2)<0，其中|u1‑u2|＝Δu；

具体操作如下：

为了搜寻导数f’(u)的零点，首先以u＝u0为中心，计算导数f’(u)在当前时刻参考点参数u0处的值f’(u0)，然后以Δu为搜索步距，朝两边搜索导数f’(u)的值f’(u0+kΔu)和f’(u0‑kΔu)，并在每步搜索中判断f’(u0+kΔu)·f’(u0+(k‑1)Δu)和f’(u0‑kΔu)·f’(u0‑(k‑1)Δu)是否小于零，搜索次数为M，搜索次数M大于3；

(a)若首次出现f’(u0+kΔu)·f’(u0+(k‑1)Δu)<0，则说明导数f’(u)的零点在[u0+(k‑

1)Δu,u0+kΔu]之间，此时记u1＝u0+(k‑1)Δu，u2＝u0+kΔu，进入步骤(3)；

(b)若首次出现f’(u0‑kΔu)·f’(u0‑(k‑1)Δu)<0，则说明导数f’(u)的零点在[u0‑kΔu,u0‑(k‑1)Δu]之间，此时记u1＝u0‑kΔu，u2＝u0‑(k‑1)Δu，进入步骤(3)；

(c)若在M次搜索之后没有找到导数f’(u)符号变化的参数u，则取参数uc＝{u|min(f’(u0+MΔu),f’(u0‑MΔu))}作为期望轨迹上离当前时刻实际点最近的点，计算当前时刻的轮廓误差为ε＝|P(uc)‑Pa|，结束计算；

(3)接着使用二分法在u1和u2之间继续搜索导数f’(u)的零点，直到精度满足要求或达到指定的迭代次数；计算当前时刻的轮廓误差为ε＝|P(uc)‑Pa|，从而得到实际点Pa与期望轨迹P(u)之间的最短距离，即轮廓误差；

具体操作如下：

设置二分法的迭代次数为N，迭代次数为N大于10，或者使导数f’(u)＝0的u的精度范围为δ；令umid＝(u1+u2)/2，计算f’(umid)，若f’(umid)·f’(u1)<0，则令u2＝umid；若f’(umid)·f’(u2)<0，则令u1＝umid；如果迭代次数没达到或者u的精度没达到，继续重复步骤(3)，否则取参数uc＝(u1+u2)/2作为期望轨迹上离当前时刻实际点最近的点；计算当前时刻的轮廓误差为ε＝|P(uc)‑Pa|；在每个控制周期内执行上述步骤(1)‑(3)，即可实现参数曲线轮廓误差实时估计。