基于轴不变量的通用7R机械臂逆解建模与解算方法
阅读：843发布：2021-01-09
专利汇可以提供基于轴不变量的通用7R机械臂逆解建模与解算方法专利检索，专利查询，专利分析的服务。并且本 发明 公开了一种基于轴不变量的通用7R 机械臂 逆解建模与解算方法，设定有7个转动轴，拾取点位于第7轴轴线上，且第5轴与第6轴不共轴的机械臂为通用7R机械臂；通用7R机械臂通过前6轴控制第7轴与期望 位置 及 姿态 对齐，使第7轴无限转动或控制第7轴满足径向对齐；将7R机械臂运动学方程采用居‑吉布斯四元数表达式进行表达，第6轴取距拾取点一定距离的点为名义拾取点，先计算通用6R机械臂的逆解，再应用数值 迭代 法，完成通用7R机械臂的运动规划与逆解计算，解决了 现有技术 中无法计算7R机械臂逆解的问题。，下面是基于轴不变量的通用7R机械臂逆解建模与解算方法专利的具体信息内容。
1.一种基于轴不变量的通用7R机械臂逆解建模与解算方法，其特征是，
设定有7个转动轴，拾取点位于第7轴轴线上，且第5轴与第6轴不共轴的机械臂为通用
7R机械臂；通用7R机械臂通过前6轴控制第7轴与期望位置及姿态对齐，使第7轴无限转动或控制第7轴满足径向对齐；
将给定第7轴期望位置 及第7轴期望姿态 的逆解问题与给定第6轴期望位置 及第
6轴期望姿态 的逆解问题等价，将7R机械臂运动学方程采用居-吉布斯四元数表达式进行表达；
基于轴不变量的通用机械臂运动规划，建立通用7R机械臂增量位姿方程，第6轴取距拾取点一定距离的点为名义拾取点，先计算通用6R机械臂的逆解，再应用数值迭代法，完成通用7R机械臂的运动规划与逆解计算。
2.根据权利要求1所述的基于轴不变量的通用7R机械臂逆解建模与解算方法，其特征是，
若给定7R轴链il7＝(i,1:7]]，il1＝03，期望位置矢量及Ju-Gibbs四元数分别记为 及则由轴不变量表征的7R机械臂运动学多项式方程为：
其中：
系统结构参数及期望Ju-Gibbs姿态四元数构成的矩阵表示为
式中，\为续行符； 分别表示轴5至轴6、轴6至轴7的零位矢量、径向矢
量； 是轴不变量 的叉乘矩阵；03＝[0 0 0]T；
表示系统结构参数的4
×4矩阵； 表示取 的第一行元素，依次类推， 表示取 的第k+1行元素；右上角标表达形式[]表示取行或列，表达形式[·]表示取所有列； 表示取3E5的第3行、第所有列；
4n5为杆件4到杆件5的坐标矢量，其是轴不变量； 为轴不变量4n5的叉乘矩阵，其余杆件同理；
消去 及 后的位姿方程(200)，是4个“4元2阶”多项式方程。
3.根据权利要求1所述的基于轴不变量的通用7R机械臂逆解建模与解算方法，其特征是，
基于轴不变量的通用机械臂运动规划，具体包括以下步骤：
【1】建立通用7R轴链的增量位姿方程；
【2】对基于偏速度迭代的通用7R机械臂进行运动规划。
4.根据权利要求3所述的基于轴不变量的通用7R机械臂逆解建模与解算方法，其特征是，
定义Ju-Gibbs增量四元数：
Ju-Gibbs规范四元数为
其中：轴不变量 为关节变量；若用 表示属性占位，则式中的表达形式 表示成员访问符；式中的表达形式幂符 表示 的x次幂；右上角角标∧或表示分隔符；
定义居-吉布斯增量四元数
其中：
5.根据权利要求4所述的基于轴不变量的通用7R机械臂逆解建模与解算方法，其特征是，
步骤【1】中，若给定6R轴链il7＝(i,1:7]，il1＝03，位置矢量及Ju-Gibbs增量四元数分别记为 及 则由Ju-Gibbs增量四元数表征的通用6R机械臂增量位姿方程表示为：
3
其中： n4为杆件3到杆件4的坐标矢量，其
是轴不变量； 为轴不变量3n4的叉乘矩阵，其余杆件同理；，\为续行符； 分别表示轴5至轴6、轴6至轴7的零位矢量、径向矢量； 是轴不变量 的叉乘矩阵；03＝[0 0
0]T；
由行四元数构成的系统结构参数矩阵表示为
其中，表示系统结构参数的4×4矩阵；右上角标表达形式[]表示取行或列，表达形式[·]表示取所有列； 表示取3E5的第i行、第所有列。
6.根据权利要求5所述的基于轴不变量的通用7R机械臂逆解建模与解算方法，其特征是，
分析通用7R机械臂的增量位姿逆解，式(239)是关于{εl|l∈[1:4}的线性方程；将式(239)重新表示为
A·[ε1 ε2 ε3 ε4]T＝b； (247)
若A-1存在，解式(245)得
[ε1 ε2 ε3 ε4]T＝A-1·b； (248)
至此，得到全部逆解。
7.根据权利要求6所述的基于轴不变量的通用7R机械臂逆解建模与解算方法，其特征是，
通过增量位置矢量 及Gibbs增量四元数 应用迭代逼近算法使通用7R机械臂位姿对齐期望的位姿。
8.根据权利要求7所述的基于轴不变量的通用7R机械臂逆解建模与解算方法，其特征是，
步骤【2】中，基于偏速度迭代优化的步骤如下：
记运动链为 l∈(i,1:6]，有
记 记期望位姿分别为 及 且有
将式(251)合写为
由式(252)应用梯度下降法得
其中：步长Step>0，Step→0；有
选择Step步长，由初态 开始迭代，直至终态
(1)确定目标函数
Goal表示 及 的方差；
(2)选择步长；
(3)进行迭代计算，得到稳态解，即为通用机械臂的位姿逆解。
9.根据权利要求8所述的基于轴不变量的通用7R机械臂逆解建模与解算方法，其特征是，
应用构造法确定步长
仅当 时，Step→0，Goal→0；
其中：ε(i,6]＝[ε1 ε2 … ε6]；由式(257)及式(253)确定步长
其中：
10.根据权利要求8所述的基于轴不变量的通用7R机械臂逆解建模与解算方法，其特征是，
迭代过程
取式(258)之步长Step，由式(259)完成迭代计算：
其中：当 时，迭代过程结束；
对于式(259)的迭代过程，则必有
δGoal≤0， (260)
即式(259)的迭代过程一定收敛；
当 时，得到稳态解φ(i,6]，即为通用机械臂的位姿逆解。
基于轴不变量的通用7R机械臂逆解建模与解算方法
技术领域
背景技术
[0002] 自主机器人研究的一个重要方面是需要解决变拓扑结构机器人的运动学建模问题。在MAS中，具有动态的图结构(Dynamic Graph Structure)，可以动态地建立基于运动轴的有向Span树，为研究可变拓扑结构(Variable Topology Structure)的机器人建模与控制奠定了基础。为此，需要提出基于轴不变量的通用机械臂逆解原理，既要建立包含坐标系、极性、结构参数、关节变量的完全参数化的正运动学模型，又要实时地计算位姿方程；一方面，可以提高机器人的自主性，另一方面，可以提高机器人位姿控制的绝对精度。
[0003] 由于通用6R机械臂不存在共点约束，现有技术中其逆解计算已经十分困难，在工程上不得不屈从于解耦约束，该约束既增加了机械臂加工及装配难度，又降低了机械臂绝对定位精度。通用7R机械臂逆解由于计算复杂度极高，在现有技术条件下根本无法实现。
发明内容
[0004] 本发明所要解决的技术问题是提供一种基于轴不变量的通用7R机械臂逆解建模与解算方法，解决了现有技术中无法计算7R机械臂逆解的问题。
[0005] 为解决上述技术问题，本发明采用以下技术方案：
[0006] 一种基于轴不变量的通用7R机械臂逆解建模与解算方法，其特征是，
[0007] 设定有7个转动轴，拾取点位于第7轴轴线上，且第5轴与第6轴不共轴的机械臂为通用7R机械臂；通用7R机械臂通过前6轴控制第7轴与期望位置及姿态对齐，使第7轴无限转动或控制第7轴满足径向对齐；
[0008] 将给定第7轴期望位置 及第7轴期望姿态 的逆解问题与给定第6轴期望位置及第6轴期望姿态 的逆解问题等价，将7R机械臂运动学方程采用居-吉布斯四元数表达式进行表达；
[0009] 基于轴不变量的通用机械臂运动规划，建立通用7R机械臂增量位姿方程，第6轴取距拾取点一定距离的点为名义拾取点，先计算通用6R机械臂的逆解，再应用数值迭代法，完成通用7R机械臂的运动规划与逆解计算。
[0010] 若给定7R轴链il7＝(i,1:7]，il1＝03，期望位置矢量及Ju-Gibbs四元数分别记为及 则由轴不变量表征的7R机械臂运动学多项式方程为：
[0011]
[0012] 其中：
[0013]
[0014] 系统结构参数及期望Ju-Gibbs姿态四元数构成的矩阵表示为
[0015]
[0016]
[0017]
[0018]
[0019] 式中，\为续行符； 分别表示轴5至轴6、轴6至轴7的零位矢量、径向矢量； 是轴不变量 的叉乘矩阵；03＝[0 0 0]T；
表示系统结构参数的4
×4矩阵； 表示取 的第一行元素，依次类推， 表示取 的第k+1行元素；右上角标表达形式[ ]表示取行或列，表达形式[·]表示取所有列； 表示取3E5的第3行、第所有列；4n5为杆件4到杆件5的坐标矢量，其是轴不变量； 为轴不变量4n5的叉乘矩阵，其余杆件同理。
表示系统结构参数的4
×4矩阵； 表示取 的第一行元素，依次类推， 表示取 的第k+1行元素；右上角标表达形式[ ]表示取行或列，表达形式[·]表示取所有列； 表示取3E5的第3行、第所有列；4n5为杆件4到杆件5的坐标矢量，其是轴不变量； 为轴不变量4n5的叉乘矩阵，其余杆件同理。
[0020] 消去τ5及τ6后的位姿方程(200)，是4个“4元2阶”多项式方程。
[0021] 基于轴不变量的通用机械臂运动规划，具体包括以下步骤：
[0022] 【1】建立通用7R轴链的增量位姿方程；
[0023] 【2】对基于偏速度迭代的通用7R机械臂进行运动规划。
[0024] 定义Ju-Gibbs增量四元数：
[0025] Ju-Gibbs规范四元数为
[0026]
[0027] 其中：轴不变量 为关节变量；若用 表示属性占位，则式中的表达形式表示成员访问符；式中的表达形式幂符 表示 的x次幂；右上角角标∧或 表示分隔符。
[0028] 定义居-吉布斯增量四元数
[0029]
[0030] 其中：
[0031]
[0032] 步骤【1】中，若给定6R轴链il7＝(i,1:7]，il1＝03，位置矢量及Ju-Gibbs增量四元数分别记为 及 则由Ju-Gibbs增量四元数表征的通用6R机械臂增量位姿方程表示为：
[0033]
[0034]
[0035] 其中： 3n4为杆件3到杆件4的坐标矢量，其是轴不变量； 为轴不变量3n4的叉乘矩阵，其余杆件同理；，\为续行符； 分
别表示轴5至轴6、轴6至轴7的零位矢量、径向矢量； 是轴不变量 的叉乘矩阵；03＝[0 0
0]T；
别表示轴5至轴6、轴6至轴7的零位矢量、径向矢量； 是轴不变量 的叉乘矩阵；03＝[0 0
0]T；
[0036] 由行四元数构成的系统结构参数矩阵表示为
[0037]
[0038] 其中，表示系统结构参数的4×4矩阵；右上角标表达形式[ ]表示取行或列，表达形式[·]表示取所有列； 表示取3E5的第i行、第所有列。
[0039] 分析通用7R机械臂的增量位姿逆解，式(239)是关于{εl|l∈[1:4}的线性方程。将式(239)重新表示为
[0040] A·[ε1 ε2 ε3 ε4]T＝b； (247)
[0041] 若A-1存在，解式(245)得
[0042] [ε1 ε2 ε3 ε4]T＝A-1·b； (248)
[0043] 至此，得到全部逆解。
[0045] 步骤【2】中，基于偏速度迭代优化的步骤如下：
[0046] 记运动链为 l∈(i,1:6]，有
[0047]
[0048] 记 记期望位姿分别为 及 且有将式(251)合写为
[0049]
[0050] 由式(252)应用梯度下降法得
[0051]
[0052] 其中：步长Step>0，Step→0；有
[0053]
[0054] 选择Step步长，由初态 开始迭代，直至终态
[0055] (1)确定目标函数
[0056]
[0057] Goal表示 及 的方差；
[0058] (2)选择步长；
[0059] (3)进行迭代计算，得到稳态解，即为通用机械臂的位姿逆解。
[0060] 应用构造法确定步长
[0061]
[0062] 仅当 时，Step→0，Goal→0；
[0063]
[0064] 其中：ε(i,6]＝[ε1 ε2 … ε6]；由式(257)及式(253)确定步长
[0065]
[0066] 其中：
[0067] 迭代过程
[0068] 取式(258)之步长Step，由式(259)完成迭代计算：
[0069]
[0070] 其中：当 时，迭代过程结束；
[0071] 对于式(259)的迭代过程，则必有
[0072] δGoal≤0， (260)
[0073] 即式(259)的迭代过程一定收敛；
[0074] 当 时，得到稳态解φ(i,6]，即为通用机械臂的位姿逆解。
[0075] 本发明所达到的有益效果：
[0076] 本发明建立了通用7R机械臂增量位姿方程，并进行实时逆解计算。特征在于：
[0077] 具有简洁、优雅的运动链符号系统，具有伪代码的功能，具有迭代式结构，保证系统实现的可靠性及机械化演算。
[0078] 具有基于轴不变量的迭代式，保证计算的实时性；实现坐标系、极性及系统结构参量的完全参数化，基于轴不变量的可逆解运动学具有统一的表达及简洁的结构化层次模型，保证位姿分析逆解的通用性。
[0080] 图1自然坐标系与轴链；
[0081] 图2固定轴不变量；
[0082] 图3为定轴转动示意图；
[0083] 图4为轴不变量的导出不变量。
具体实施方式
[0084] 下面对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。
[0085] 定义1自然坐标轴：称与运动轴或测量轴共轴的，具有固定原点的单位参考轴为自然坐标轴，亦称为自然参考轴。
[0086] 定义2自然坐标系：如图1所示，若多轴系统D处于零位，所有笛卡尔体坐标系方向一致，且体坐标系原点位于运动轴的轴线上，则该坐标系统为自然坐标系统，简称自然坐标系。
[0087] 自然坐标系优点在于：(1)坐标系统易确定；(2)零位时的关节变量为零；(3)零位时的系统姿态一致；(4)不易引入测量累积误差。
[0088] 由定义2可知，在系统处于零位时，所有杆件的自然坐标系与底座或世界系的方向一致。系统处于零位即 时，自然坐标系 绕轴矢量 转动角度 将 转至F[l] [l]
； 在 下的坐标矢量与 在F 下的坐标矢量 恒等，即有
； 在 下的坐标矢量与 在F 下的坐标矢量 恒等，即有
[0089]
[0090] 由上式知， 或 不依赖于相邻的坐标系 及F[l]；故称 或 为轴不变量。在不强调不变性时，可以称之为坐标轴矢量(简称轴矢量)。 或 表征的是体 与体l共有的参考单位坐标矢量，与参考点 及Ol无关。体 与体l即为杆件或轴。
[0091] 轴不变量与坐标轴具有本质区别：
[0092] (1)坐标轴是具有零位及单位刻度的参考方向，可以描述沿该方向平动的位置，但不能完整描述绕该方向的转动角度，因为坐标轴自身不具有径向参考方向，即不存在表征转动的零位。在实际应用时，需要补充该轴的径向参考。例如：在笛卡尔系F[l]中，绕lx转动，需以ly或lz为参考零位。坐标轴自身是1D的，3个正交的1D参考轴构成3D的笛卡尔标架。
[0093] (2)轴不变量是3D的空间单位参考轴，其自身就是一个标架。其自身具有径向参考轴，即参考零位。空间坐标轴及其自身的径向参考轴可以确定笛卡尔标架。空间坐标轴可以反映运动轴及测量轴的三个基本参考属性。
[0094] 已有文献将无链指标的轴矢量记为 并称之为欧拉轴(Euler Axis)，相应的关节角称为欧拉角(Euler Angle)。本申请之所以不再沿用欧拉轴，而称之为轴不变量，是因为轴不变量具有以下属性：
[0095] 【1】给定旋转变换阵 因其是实矩阵，其模是单位的，故其有一个实特征值λ1及两个互为共轭的复特征值λ2＝eiφ及λ3＝e-iφ；其中：i为纯虚数。因此，|λ1|·||λ2||·||λ3||＝1，得λ1＝1。轴矢量 是实特征值λ1＝1对应的特征矢量，是不变量；
[0096] 【2】是3D参考轴，不仅具有轴向参考方向，而且具有径向参考零位，将在3.3.1节予以阐述。
[0097] 【3】在自然坐标系下： 即轴不变量 是非常特殊的矢量，它对时间的导数也具有不变性，且有非常优良的数学操作性能；
[0098] 对轴不变量而言，其绝对导数就是其相对导数。因轴不变量是具有不变性的自然参考轴，故其绝对导数恒为零矢量。因此，轴不变量具有对时间微分的不变性。有：
[0099]
[0100] 【4】在自然坐标系统中，通过轴矢量 及关节变量 可以直接描述旋转坐标阵没有必要为除根之外的杆件建立各自的体系。同时，以唯一需要定义的根坐标系为参考，可以提高系统结构参数的测量精度；
[0101] 【5】应用轴矢量 的优良操作，将建立包含拓扑结构、坐标系、极性、结构参量及力学参量的完全参数化的统一的多轴系统运动学及动力学模型。
[0102] 因基矢量el是与F[l]固结的任一矢量，基矢量 是与 固结的任一矢量，又 是F[l]及 共有的单位矢量，故 是F[l]及 共有的基矢量。因此，轴不变量 是F[l]及 共有的参考基。轴不变量是参数化的自然坐标基，是多轴系统的基元。固定轴不变量的平动与转动与其固结的坐标系的平动与转动等价。
[0103] 在系统处于零位时，以自然坐标系为参考，测量得到坐标轴矢量 在运动副运动时，轴矢量 是不变量；轴矢量 及关节变量 唯一确定运动副 的转动关系。
[0104] 因此，应用自然坐标系统，当系统处于零位时，只需确定一个公共的参考系，而不必为系统中每一杆件确定各自的体坐标系，因为它们由轴不变量及自然坐标唯一确定。当进行系统分析时，除底座系外，与杆件固结的其它自然坐标系只发生在概念上，而与实际的测量无关。自然坐标系统对于多轴系统(MAS)理论分析及工程作用在于：
[0105] (1)系统的结构参数测量需要以统一的参考系测量；否则，不仅工程测量过程烦琐，而且引入不同的体系会引入更大的测量误差。
[0106] (2)应用自然坐标系统，除根杆件外，其它杆件的自然坐标系统由结构参量及关节变量自然确定，有助于MAS系统的运动学与动力学分析。
[0107] (3)在工程上，可以应用激光跟踪仪等光学测量设备，实现对固定轴不变量的精确测量。
[0108] (4)由于运动副R及P、螺旋副H、接触副O是圆柱副C的特例，可以应用圆柱副简化MAS运动学及动力学分析。
[0109] 定义3不变量：称不依赖于一组坐标系进行度量的量为不变量。
[0110] 定义4转动坐标矢量：绕坐标轴矢量 转动到角位置 的坐标矢量 为
[0111]
[0112] 定义5平动坐标矢量：沿坐标轴矢量 平动到线位置 的坐标矢量 为
[0113]
[0114] 定义6自然坐标：以自然坐标轴矢量为参考方向，相对系统零位的角位置或线位置，记为ql，称为自然坐标；称与自然坐标一一映射的量为关节变量；其中：
[0115]
[0116] 定义7机械零位：对于运动副 在初始时刻t0时，关节绝对编码器的零位 不一定为零，该零位称为机械零位；
[0117] 故关节 的控制量 为
[0118]
[0119] 定义8自然运动矢量：将由自然坐标轴矢量 及自然坐标ql确定的矢量 称为自然运动矢量。其中：
[0120]
[0121] 自然运动矢量实现了轴平动与转动的统一表达。将由自然坐标轴矢量及关节确定的矢量，例如 称为自由运动矢量，亦称为自由螺旋。显然，轴矢量 是特定的自由螺旋。
[0122] 定义9关节空间：以关节自然坐标ql表示的空间称为关节空间。
[0123] 定义10位形空间：称表达位置及姿态(简称位姿)的笛卡尔空间为位形空间，是双矢量空间或6D空间。
[0124] 定义11自然关节空间：以自然坐标系为参考，通过关节变量 表示，在系统零位时必有 的关节空间，称为自然关节空间。
[0125] 如图2所示，给定链节 原点Ol受位置矢量 约束的轴矢量 为固定轴矢量，记为 其中：
[0126]
[0127] 轴矢量 是关节自然坐标的自然参考轴。因 是轴不变量，故称 为固定轴不变量，它表征了运动副 的结构关系，即确定了自然坐标轴。固定轴不变量 是链节 结构参数的自然描述。
[0128] 定义12自然坐标轴空间：以固定轴不变量作为自然参考轴，以对应的自然坐标表示的空间称为自然坐标轴空间，简称自然轴空间。它是具有1个自由度的3D空间。
[0129] 如图2所示， 及 不因杆件Ωl的运动而改变，是不变的结构参考量。 确定了轴l相对于轴 的五个结构参数；与关节变量ql一起，完整地表达了杆件Ωl的6D位形。给定时，杆件固结的自然坐标系可由结构参数 及关节变量 唯一确定。称轴不变量 固定轴不变量 关节变量 及 为自然不变量。显然，由固定轴不变量及关节变量 构成的关节自然不变量 与由坐标系 至F[l]确定的空间位形
具有一一映射关系，即
具有一一映射关系，即
[0130]
[0131] 给定多轴系统D＝{T,A,B,K,F,NT}，在系统零位时，只要建立底座系或惯性系，以及各轴上的参考点Ol，其它杆件坐标系也自然确定。本质上，只需要确定底座系或惯性系。
[0132] 给定一个由运动副连接的具有闭链的结构简图，可以选定回路中任一个运动副，将组成该运动副的定子与动子分割开来；从而，获得一个无回路的树型结构，称之为Span树。T表示带方向的span树，以描述树链运动的拓扑关系。
[0133] I为结构参数；A为轴序列，F为杆件参考系序列，B为杆件体序列，K为运动副类型序列，NT为约束轴的序列即非树。 为取轴序列 的成员。转动副R，棱柱副P，螺旋副H，接触副O是圆柱副C的特例。
[0134] 描述运动链的基本拓扑符号及操作是构成运动链拓扑符号系统的基础，定义如下：
[0135] 【1】运动链由偏序集合(]标识。
[0136] 【2】A[l]为取轴序列A的成员；因轴名l具有唯一的编号对应于A[l]的序号，故A[l]计算复杂度为O(1)。
[0137] 【3】 为取轴l的父轴；轴 的计算复杂度为O(1)。计算复杂度O()表示计算过程的操作次数，通常指浮点乘与加的次数。以浮点乘与加的次数表达计算复杂度非常烦琐，故常采用算法循环过程中的主要操作次数；比如：关节位姿、速度、加速度等操作的次数。
[0138] 【4】 为取轴序列 的成员； 计算复杂度为O(1)。
[0139] 【5】llk为取由轴l至轴k的运动链，输出表示为 且 基数记为|llk|。llk执行过程：执行 若 则执行 否则，结束。llk计算复杂度为O(|llk|)。
[0140] 【6】ll为取轴l的子。该操作表示在 中找到成员l的地址k；从而，获得轴l的子A[k]。因 不具有偏序结构，故ll的计算复杂度为
[0141] 【7】lL表示获得由轴l及其子树构成的闭子树，lL为不含l的子树；递归执行ll，计算复杂度为
[0142] 【8】支路、子树及非树弧的增加与删除操作也是必要的组成部分；从而，通过动态Span树及动态图描述可变拓扑结构。在支路llk中，若 则记即 表示在支路中取成员m的子。
[0143] 定义以下表达式或表达形式：
[0144] 轴与杆件具有一一对应性；轴间的属性量 及杆件间的属性量 具有偏序性。
[0145] 约定： 表示属性占位；若属性p或P是关于位置的，则 应理解为坐标系 的原点至F[l]的原点；若属性p或P是关于方向的，则 应理解为坐标系 至F[l]。
[0146] 及 应分别理解为关于时间t的函数 及 且 及 是t0时刻的常数或常数阵列。但是正体的 及 应视为常数或常数阵列。
[0147] 本申请中约定：在运动链符号演算系统中，具有偏序的属性变量或常量，在名称上包含表示偏序的指标；要么包含左上角及右下角指标，要么包含右上角及右下角指标；它们的方向总是由左上角指标至右下角指标，或由右上角指标至右下角指标，本申请中为叙述简便，有时省略方向的描述，即使省略，本领域技术人员通过符号表达式也可以知道，本申请中采用的各参数，对于某种属性符，它们的方向总是由偏序指标的左上角指标至右下角指标，或由右上角指标至右下角指标。例如： 可简述为(表示由k至l)平动矢量；rlk表示(由k至l的)线位置；krl表示(由k至l的)平动矢量；其中：r表示“平动”属性符，其余属性符对应为：属性符φ表示“转动”；属性符Q表示“旋转变换矩阵”；属性符l表示“运动链”；属性符u表示“单位矢量”；属性符ω表示“角速度”；角标为i表示惯性坐标系或大地坐标系；其他角标可以为其他字母，也可以为数字。
[0148] 本申请的符号规范与约定是根据运动链的偏序性、链节是运动链的基本单位这两个原则确定的，反映了运动链的本质特征。链指标表示的是连接关系，右上指标表征参考系。采用这种符号表达简洁、准确，便于交流与书面表达。同时，它们是结构化的符号系统，包含了组成各属性量的要素及关系，便于计算机处理，为计算机自动建模奠定基础。指标的含义需要通过属性符的背景即上下文进行理解；比如：若属性符是平动类型的，则左上角指标表示坐标系的原点及方向；若属性符是转动类型的，则左上角指标表示坐标系的方向。
[0149] (1)lS－杆件l中的点S；而S表示空间中的一点S。
[0150] (2) －杆件k的原点Ok至杆件l的原点Ol的平动矢量；
[0151] 在自然坐标系F[k]下的坐标矢量，即由k至l的坐标矢量；
[0152] (3) －原点Ok至点lS的平动矢量；
[0153] 在F[k]下的坐标矢量；
[0154] (4) －原点Ok至点S的平动矢量；
[0155] 在F[k]下的坐标矢量；
[0157] －运动副 的轴矢量；
[0158] 及 分别在 及F[l]下的坐标矢量； 是轴不变量，为一结构常数；
[0159] 为转动矢量，转动矢量/角矢量 是自由矢量，即该矢量可自由平移；
[0160] (6) －沿轴 的线位置(平动位置)，
[0161] －绕轴 的角位置，即关节角、关节变量，为标量；
[0162] (7)左下角指标为0时，表示机械零位；如：
[0163] －平动轴 的机械零位，
[0164] －转动轴 的机械零位；
[0165] (8)0－三维零矩阵；1－三维单位矩阵；
[0166] (9)约定：“\”表示续行符； 表示属性占位；则
[0167] 幂符 表示 的x次幂；右上角角标∧或 表示分隔符；如： 或 为 的x次幂。
[0168] 表示 的转置，表示对集合转置，不对成员执行转置；如：
[0169] 为投影符，表示矢量或二阶张量对参考基的投影矢量或投影序列，即坐标矢量或坐标阵列，投影即是点积运算“·”；如：位置矢量 在坐标系F[k]中的投影矢量记为[0170] 为叉乘符；如： 是轴不变量 的叉乘矩阵；给定任一矢量 的叉乘矩阵为叉乘矩阵是二阶张量。
[0171] 叉乘符运算的优先级高于投影符 的优先级。投影符 的优先级高于成员访问符或 成员访问符 优先级高于幂符
[0172] (10)单位矢量在大地坐标系的投影矢量 单位零位矢量
[0173] (11) －零位时由原点 至原点Ol的平动矢量，且记 表示位置结构参数。
[0174] (12)iQl，相对绝对空间的旋转变换阵；
[0175] (13)以自然坐标轴矢量为参考方向，相对系统零位的角位置或线位置，记为ql，称为自然坐标；关节变量 自然关节坐标为φl；
[0176] (14)对于一给定有序的集合r＝[1,4,3,2]T，记r[x]表示取集合r的第x行元素。常记[x]、[y]、[z]及[ω]表示取第1、2、3及4列元素。
[0177] (15)ilj表示由i到j的运动链；llk为取由轴l至轴k的运动链；
[0178] 给定运动链 若n表示笛卡尔直角系，则称 为笛卡尔轴链；若n表示自然参考轴，则称 为自然轴链。
[0179] (16)Rodrigues四元数表达形式：
[0180] 欧拉四元数表达形式：
[0181] 不变量的四元数(也称为轴四元数)表达形式
[0182] 分块矩阵的高维行列式计算：
[0183] 记<1:n>表示自然数[1:n]的全排列，共有n！个实例。给定属于数域的大小为n×n的矩阵M，其j行i列元素记为 根据行列式定义得
[0184]
[0185] 其中：I[i1,…in]表示排列的逆序个数。式(1)计算复杂度为：n！次n个数积及n！次加法，具有指数计算复杂度，只能适用于维度较小的行列式。对于维度较大的行列式，通常应用Laplace公式进行递规运算，记 为 的伴随矩阵(Adjugate Matrix)，则有
[0186]
[0187] 更简单的算法通常应用高斯消去法或LU分解法，先通过初等变换将矩阵变为三角阵或三角阵的乘积，后计算行列式。上述针对数域的行列式计算方法不适用于高维度的多项式矩阵，需要引入分块矩阵的行列式计算方法。计算矢量多项式(Vector Polynomial)的行列式是一个特定的分块矩阵行列式的计算问题，它在矢量层次上表达了矢量与行列式的内在联系。而分块矩阵行列式计算则从矩阵层次上表达分块矩阵与行列式的内在规律。
[0188] 若给定矢量多项式 其中： 及 为3D坐标矢量， 为多项式变量序列；若约定
[0189]
[0190] 则有
[0191]
[0192] 上式的推导步骤为：因
[0193]
[0194] 故式(4)成立。
[0195] 式(3)及式(4)可以推广至n维空间。式(3)有助于从矢量层次上分析行列式的内在规律；比如，当任两矢量平行或三矢量共面时，对应的行列式为零。式(4)表明：矢量多项式的行列式易导致“组合爆炸”。
[0196] 实施例1
[0197] 给定2个2维行矢量多项式一方面，由式(4)得
[0198]
[0199] 另一方面，
[0200]
[0201] 上面的结果验证了式(4)的正确性。
[0202] 给出分块矩阵的行列式计算定理：
[0203] 若记大小为(n+m)·(n+m)的方阵为M，大小为n·n的矩阵 是方阵M的前n行及任意n列元素构成的子矩阵，大小为m·m的矩阵 是方阵M后m行及剩余m列元素构成的子矩阵；由升序排列的矩阵列序号构成的序列cn及cm是序列[1:m+n]的子集，[cn,cm]∈<1:n+m>，且有cm∪cn＝[1:m+n]；则方阵M行列式与分块矩阵 及 的行列式关系为
[0204]
[0205] 对行列式进行行阶梯化计算原理：
[0206] 对于S×S矩阵，其每一项是关于τ1的n阶多项式。计算该矩阵的行列式时，可通过初等行变换将原行列式变为上三角行列式，再将非零的对角线元素相乘，得到行列式的多项式表达式。因该式为0，求得到τ1的所有解。
[0207] 阶梯化的具体方法为，先对行列式第一列的最高阶次由高到低进行排序，再进行最多(S-1)×n次初等行变换消元，得到第一列只有第一个元素不为0的行列式。再对该行列式第1行及1列的余子式进行初等行变换消元，依次迭代求解。
[0208] 实施例2
[0209] 通过矩阵的初等行变换，得到 的行阶梯矩阵。
[0210] 步骤为：rk代表第k行。得
[0211]
[0212] 则得
[0213] 基于“N进位字”的N阶多项式系统：
[0214] 若n个“n元1阶”多项式幂积 中独立变量重复出现N次，则得到n个“n元N阶”多项式系统 “n元N阶多项式系统”与“n位N进位字” 同构。
[0215]
[0216]
[0217]
[0218]
[0219] n个“n元N阶”多项式系统的Dixon多项式：
[0220] 引入辅助变量[y2,y3,…,yn]，且有
[0221]
[0222] 在多元多重多项式(8)中，用辅助变量Ym的前m个依次替换原变量(Original Variables)Xn中的m个变量，记“|”为替换操作符，得到增广的(Extended)多项式[0223]
[0224]
[0225] 式中右上角标α、α表示幂；
[0226] 由式(6)及式(12)得
[0227]
[0228] 其中：
[0229]
[0230] 定义可分离组合变量 及 如下：
[0231]
[0232] 由式(14)及式(15)可知：替换式 是关于 及 的双重线性型。相应地，用辅助变量替换的多项式系统记为
[0233]
[0234] 给定n个“n元N阶”多项式系统 定义其Dixon多项式为
[0235]
[0236] 由式(17)得
[0237]
[0238] 式(15)中分离变量与文献不同：原变量Xn-1被辅助变量Yn-1替换的次序不同，Dixon多项式也不同。考虑式(13)及式(18)得该多项式的Dixon行列式
[0239]
[0240] 在笛卡尔空间下，由位置矢量或转动矢量构成的行列式表示矢量张成空间的容积(Volume)；在不同笛卡尔空间下具有容积的不变性。其中：
[0241]
[0242] n个“n元N阶”多项式的Dixon行列式的阶次及替换变量项数分别为：
[0243]
[0244] n个“n元N阶”Dixon矩阵：
[0245] 给定n个“n元N阶”多项式系统Fn(Yn-1|Xn-1)，n≥2；存在与消去变量x2,…,xn无关的Dixon矩阵SΘS(x1)，其Dixon多项式 表示为分离变量 及 的双重线性型：
[0246]
[0247] α[l]∈[0,N·(n-l+1)-1]],l∈[2:n] (23)
[0248] 为大小为S×S的Dixon矩阵，其第[i][j]成员为单变量x1的N阶多项式：
[0249]
[0250] 其中：
[0251]
[0252] 若
[0253]
[0254] 则有
[0255]
[0256] n个“n元N阶”多项式的Dixon消元与求解
[0257] 考虑式(22)，若 故得
[0258] Det(SΘS(x1))＝0。 (28)
[0259] 称式(28)中“n个n元”为Dixon消元的必要条件，从而获得可行解。若SΘS存在零行或零列向量，则无法建立x1的多项式方程；此时，通过除标量积之外的初等变换，将SΘS变为行阶梯(Row Echelon)矩阵Ech(SΘS)；在计算该矩阵的杻轴(Pivot)的积之后得方阵即在SΘS中选取S′个独立的列向量。
[0260] 任一个n个“n元N阶”多项式系统 的实例(简称多项式)记为其中： 且有 根据
的多项式确定Dixon矩阵、分离变量 及 选取 及 满足
的多项式确定Dixon矩阵、分离变量 及 选取 及 满足
[0261]
[0262] 确定双线性型
[0263]
[0264]
[0265] 其中： 中与 对应的各列线性独立。因 由式(22)及式(25)得
[0266]
[0267] 称其为结式或消去式。式(32)是单变量x1的多项式方程；消去了n-1个未知量；从而，可以获得单变量x1的可行解。若x1同时满足
[0268]
[0269] 则x1为正确解。将已解的x1代入式(34)，因式(32)成立且 任意，故得
[0270]
[0271] 即有
[0272]
[0273] 若有必要条件
[0274]
[0275] 成立，解式(35)，得被消去变量 的解；否则，需要结合式(16)得到全部解。考虑式(25)，因式(22)两边的x1阶次相等，故必有
[0276]
[0277] 若同时满足
[0278]
[0279] 则由式(35)能解得 中n-1个互不相同的组合变量；从而，得到所有独立变量的解。
[0280] 给定n个“n元N阶”多项式 Dixon矩阵计算步骤如下：
[0281] ①确定系统结构。方程数及独立变量数记为n；独立变量记为Xn；多项式复合变量记为 替换变量记为 替换变量数为n-1；大小为S·S的Dixon矩阵记为 其成员系数如式(24)所示，其中：S由式(32)确定；待消去变量为x1。
[0282] ②由式(8)得xα与 对应关系，表达式(11)中 至多有S项。
[0283] ③根据式(19)及Sarrus规则，计算Dixon(Fn(Yn-1|Xn-1))；根据 对应的N进位字运算结果，完成多项式合并。
[0284] ④Dixon矩阵成员如式(32)所示，由式(32)计算Dixon矩阵SΘS的(n+1)·S2个系数。
[0285] ⑤当满足式(37)及式(38)直接解判别准则时，由式(34)及式(35)得全部数值解。
[0286] 实施例3
[0287] 对多项式系统(39)进行Dixon消元。
[0288]
[0289] 步骤为：该式是4个“4元1阶”多项式系统，满足Dixon消元条件。由式(19)及式(22)，得
[0290]
[0291] 其中：
[0292]
[0293] 由式(34)及式(40)得5个解：
[0294]
[0295]
[0296]
[0297] 其中： 不是该方程组的解。将其它解分别代入式(35)。当 时，由式(35)得
[0298]
[0299] 解得：τ3＝1，τ4＝-2。将 τ3及τ4代入式(39)得τ2＝1。同样，可得其他三组解。显然，因变量不满足式(26)，式(40)所示的Dixon矩阵不对称。该例表明Dixon行列式为零对于多重线性多项式系统是充分的。
[0300] 基于轴不变量的定轴转动
[0301] 如图3所示，给定轴矢量 及与其固结的单位矢量 在转动前，对于单位矢量对系统零位轴 的投影矢量为 对系统径向轴 的矩矢量为 径向矢量为
[0302] 轴矢量 相对于杆件 及Ωl或自然坐标系 及F[l]是固定不变的，故称该转动为定轴转动。单位矢量 绕轴 转动 后，转动后的零位矢量 对系统零位轴 的投影矢量为 转动后的零位矢量 对系统径向轴 的矩矢量为轴向分量为 故得具有链指标的Rodrigues矢量方程
[0303]
[0304] 因单位矢量 是任意的且 得具有链指标的Rodrigues转动方程
[0305]
[0306] 若 由式(42)，得 若 即坐标系 与F[l]的方向一致，由式[l]
(42)可知：反对称部分 必有 因此，系统零位是自然坐标系 与F
重合的充分必要条件，即初始时刻的自然坐标系方向一致是系统零位定义的前提条件。利用自然坐标系可以很方便地分析多轴系统运动学和动力学。
(42)可知：反对称部分 必有 因此，系统零位是自然坐标系 与F
重合的充分必要条件，即初始时刻的自然坐标系方向一致是系统零位定义的前提条件。利用自然坐标系可以很方便地分析多轴系统运动学和动力学。
[0307]
[0308] 式(43)是关于 和 的多重线性方程，是轴不变量 的二阶多项式。给定自然零位矢量 作为 的零位参考，则 及 分别表示零位矢量及径向矢量。式
(43)即为 对称部分 表示零位轴张量，反对称部分
表示径向轴张量，分别与轴向外积张量 正交，从而确定三维自然轴空间；
式(43)仅含一个正弦及余弦运算、6个积运算及6个和运算，计算复杂度低；同时，通过轴不变量 及关节变量 实现了坐标系及极性的参数化。
(43)即为 对称部分 表示零位轴张量，反对称部分
表示径向轴张量，分别与轴向外积张量 正交，从而确定三维自然轴空间；
式(43)仅含一个正弦及余弦运算、6个积运算及6个和运算，计算复杂度低；同时，通过轴不变量 及关节变量 实现了坐标系及极性的参数化。
[0309] 对于轴链 有
[0310]
[0311] 由式(44)及式(43)得 则 是 及 的多重线性型，其中：l∈ilk。式(43)可表示为
[0312]
[0313] 称(45)为改进的Cayley变换。即有
[0314]
[0315] 规范的位置方程为
[0316]
[0317] “居-吉布斯”四元数的确定：
[0318] 对于任意杆件 定义与欧拉四元数同构的“居-吉布斯”(Ju-Gibbs)规范四元数：
[0319]
[0320] 其中： 为Gibbs矢量。Gibbs共轭四元数为：
[0321]
[0322] 其中：
[0323]
[0324] 显然， 为 模的平方。因居-吉布斯四元数是四元数，故满足四元数乘法运算[0325]
[0326] 其中：
[0327]
[0328] 由式(52)得
[0329]
[0330] 习惯上，单关节及运动链的期望姿态以规范的Ju-Gibbs四元数(简称规范Ju-Gibbs四元数，即“标部”为1的四元数)表示；但是它们积运算通常是不规范的，即其标部不为1。由式(53)可知：只有给定轴l及 的规范Ju-Gibbs四元数，且两轴正交， 才为规范四元数。
[0331] 由式(53)得
[0332]
[0333] 由四维复数性质得
[0334]
[0335] 记 由式(52)得
[0336]
[0337] 故 为单位Ju-Gibbs四元数。
[0338] 由式(48)至式(50)及式(55)得
[0339]
[0340] 由式(50)、式(54)及式(57)得
[0341]
[0342] 类DCM及性质：
[0343] 对于轴链 规范的姿态方程为：
[0344]
[0345] 由式(59)得
[0346]
[0347] 式中， 为旋转变换矩阵； 表示用辅助变量yl的前l个依次替换原变量τl中的l个变量，记“|”为替换操作符；
[0348] 其中：
[0349]
[0350] 由式(61)可知：iQn及 是关于τk的n重2阶多项式。由式(60)可知：因 与 类似，故称之为类DCM(DCM，方向余弦矩阵)。由式(62)得
[0351]
[0352] 显然，类DCM可以通过Ju-Gibbs四元数表达。因此，式(59)姿态方程及式(47)位置方程是关于Ju-Gibbs四元数的表达式。
[0353] 分块方阵的逆：
[0354] 若给定可逆方阵K、B及C，其中B及C分别为l×l、c×c的方阵；A、D分别为l×c、c×l的矩阵，且
[0355]
[0356] 则有
[0357]
[0358] 基于Ju-Gibbs四元数的指向对齐原理
[0359] 考虑轴链ill，其中 若使轴矢量 与期望轴矢量 对齐，则至少存在一个Ju-Gibbs方向四元数
[0360]
[0361] 其中：
[0362]
[0363]
[0364] 且有
[0365]
[0366] 上式的具体建立步骤为：
[0367] 定轴转动的Cayley正变换为
[0368]
[0369] 由式(69)得
[0370]
[0371] 由式(70)得
[0372]
[0373] 及
[0374]
[0375] 因 及 为单位矢量，若 且 则有
[0376]
[0377] 式(73)表明 及 相互正交。由式(72)及(73)得最优的轴矢量
[0378]
[0379] 及
[0380]
[0381] 由式(74)及(75)得式(66)。若 或 则由式(70)得
[0382]
[0383] 由式(76)得
[0384]
[0385] 因 故由式(77)得
[0386]
[0387] 及
[0388]
[0389] 由式(78)及式(79)分别得式(67)及式(68)。证毕。
[0390] 基于Ju-Gibbs四元数的指向对齐原理表明：至少存在一个期望的Ju-Gibbs四元数使单位矢量 与期望单位矢量 对齐。
[0391] 实施例4
[0392] 考虑轴链il6，由基于Ju-Gibbs四元数的指向对齐原理得
[0393]
[0394] 基于Ju-Gibbs四元数的2R方向逆解
[0395] 以Ju-Gibbs四元数指向对齐为基础，阐述2R方向逆解。
[0396] 若给定6R转动链il6＝(i,1:6]]，记第5轴关节Ju-Gibbs方向四元数期望为 及第3轴关节Ju-Gibbs规范四元数为 则有指向对齐时的逆解
[0397]
[0398]
[0399] 其中：
[0400]
[0401]
[0402] Ju-Gibbs方向四元数 满足
[0403]
[0404] 上式的具体建立步骤为：
[0405] 首先，考虑基于欧拉四元数的姿态对齐。由
[0406]
[0407] 得
[0408]
[0409] 由式(86)得
[0410]
[0411] 其中：
[0412]
[0413] 由式(64)及式(88)得
[0414]
[0415] 由式(87)得
[0416]
[0417] 其中：
[0418]
[0419] 由式(89)及式(90)得
[0420]
[0421] 由式(92)得
[0422]
[0423] 若 式(92)与式(93)两边对应相除得
[0424]
[0425] 其次，考虑Ju-Gibbs四元数的指向对齐。因 故得(82)。由式(58)得
[0426]
[0427] 以规范Ju-Gibbs四元数表征关节变量，由式(53)得
[0428]
[0429] 由式(96)得
[0430]
[0431] 因3n4及4n5独立，由式(83)可知，3E5必存在。显然， 由3n4及4n5唯一确定。将式(96)、式(97)及式(91)代入式(95)得
[0432]
[0433] 若 由式(98)第1行得
[0434]
[0435] 将式(99)代入式(98)得
[0436]
[0437] 由式(94)及式(100)可知两种原理等价。由式(100)第2、3行得
[0438]
[0439] 由式(101)可知式(80)成立。因式存(100)在4个等式，2个独立变量，由式(101)及式(100)中第4行得约束方程
[0440]
[0441] 若 由式(93)得C4 C5＝0；由式(87)得
[0442]
[0443] 由式(103)得
[0444]
[0445] 显然， 当 时，若 由式(102)得
[0446]
[0447] 若 由式(102)得
[0448]
[0449] 由式(106)可知：要么 要么 由式(106)及式(101)可知，式(81)亦成立。当 时，若 则式(80)亦成立。证毕。
[0450] Ju-Gibbs方向四元数
[0451] 由式(80),(81)和(99)得
[0452]
[0453] 且
[0454]
[0455] 把 代入以上两个方程中得到式(84)，它显示的是特定的Ju-Gibbs四元数，称为Ju-Gibbs方向四元数。
[0456] 实施例5
[0457] 继实施例4，考虑轴链il6，且有3n4＝1[x]，4n5＝1[y]，由式(83)得3E5＝1。由式(80)及式(81)得
[0458]
[0459]
[0460]
[0461] 由式(81)得：
[0462]
[0463] 基于Ju-Gibbs四元数的2R方向逆解原理表明Ju-Gibbs四元数与欧拉四元数同构；同时，式(63)表明以Ju-Gibbs四元数表示的类DCM与DCM同构。因此，应用Ju-Gibbs四元数可以完整表达位姿关系。
[0464] 若给定运动链iln，k∈iln，期望规范Ju-Gibbs四元数 及期望位置矢量 考虑式(47)及式(95)；则位置及指向对齐关系表示为
[0465]
[0466] 且具有模不变性
[0467]
[0468] 与欧拉四元数及对偶四元数相比，Ju-Gibbs四元数表征的位姿对齐不存在冗余方程；通过指向对齐，可以求解第4轴及第5轴的关节变量，为6R及7R机械臂逆解奠定了基础。
[0469] 基于轴不变量的Dixon行列式计算原理：
[0470] 下面基于轴不变量，提出径向不变量及运动链的Dixon行列式基本性质，为基于轴不变量的机器人逆运动学分析奠定基础。
[0471] 【1】轴不变量
[0472] 首先，轴不变量与坐标轴具有本质区别：坐标轴是具有零位及单位刻度的参考方向，可以描述沿轴向平动的线位置，但不能完整描述绕轴向的角位置，因为坐标轴自身不具有径向参考方向，即不存在表征转动的零位。在实际应用时，需要补充坐标轴的径向参考。坐标轴自身是1D的，3个正交的坐标轴构成3D的笛卡尔标架；轴不变量是3D空间单位参考轴(简称3D参考轴)，具有径向参考零位。“3D参考轴”及其径向参考零位可以确定对应的笛卡尔系。以自然坐标系为基础的轴不变量可以准确地反映运动轴及测量轴的“共轴性”、“极性”与“零位”三个基本属性。
[0473] 其次，轴不变量与欧拉轴具有本质的区别：方向余弦矩阵(DCM)是实矩阵，轴矢量是DCM的特征值1对应的特征矢量，是不变量；固定轴不变量是“3D参考轴”，不仅具有原点及轴向，也有径向参考零位；在自然坐标系下，轴不变量不依赖于相邻固结的自然坐标系，即在相邻固结的自然坐标系下具有不变的自然坐标；轴不变量具有幂零特性等优良的数学操作功能；在自然坐标系统中，通过轴不变量及关节坐标，可以唯一确定DCM及参考极性；没有必要为每一个杆件建立各自的体系，可以极大地简化建模的工作量。
[0474] 同时，以唯一需要定义的笛卡尔直角坐标系为参考，测量轴不变量，可以提高结构参数的测量精度。基于轴不变量的优良操作及属性，可以建立包含拓扑结构、坐标系、极性、结构参量及动力学参量的迭代式的运动学及动力学方程。
[0475] 由式(59)及式(47)可知：多轴系统的姿态及位置方程本质上是多元二阶多项式方程，其逆解本质上归结于多元二阶多项式的消元问题，包含Dixon矩阵及Dixon行列式计算的两个子问题。用式(47)的表达3R机械臂位置方程，是3个“3元2阶”多项式，应用Dixon消元方法计算逆解，有两个替换变量，在计算8×8的Dixon行列式时，最大可能的阶次为16。由式(4)可知：行列式计算是一个排列过程，面临着“组合爆炸”的难题。
[0476] 所有的不在确定的多项式时间内可解的问题称为NP问题。非确定性算法将问题分解为“猜测”与“验证”两个阶段：算法的“猜测”阶段具有非确定性，算法的“验证”阶段具有确定性，通过验证来确定猜测的解是否正确。假如可以在多项式时间内计算出来，就称为多项式非确定性问题。多元多项式的消元通常被认为是NP问题。通常应用 基进行多元多项式的消元，不得不求助于启发式的“猜测”与“验证”来解决问题。
[0477] 【2】径向不变量
[0478] 结构参数 及 是链节l的结构参量，在系统零位时，它们可以通过外部测量得到。如图4所示，零位矢量、径向矢量及轴向矢量是与转动角无关的不变量。其中，零位矢量是特定的径向矢量。
[0479] 任一个矢量可以分解为零位矢量及轴向矢量，故有
[0480]
[0481] 其中：
[0482]
[0483]
[0484] 考虑链节 其D-H参数有
[0485]
[0486] 显然，是轴l及 的公垂线或公共径向矢量，是轴l的轴向矢量。由式(112)可知：任一个结构参数矢量 可分解为与坐标系为无关的零位不变量 及轴向不变量 它们的径向矢量记为 结构参数矢量 及轴不变量 唯一确定径向坐标系，具有2个独立维度。
若两个轴向不变量 及 共线，则记为
若两个轴向不变量 及 共线，则记为
[0487]
[0488] 若两个零位不变量 及 与任两个径向不变量 及 共面，则记为
[0489]
[0490] 因此，称式(110)所示的轴向不变量及零位不变量是结构参数矢量对自然轴的分解。
[0491] 由式(113)及式(114)可知：同一个轴的三个径向矢量的行列式为零；同一个轴的任意两个轴向矢量的行列式为零。可以用轴不变量及其导出的不变量来简化Dixon行列式计算。
[0492] 由轴不变量导出的零位矢量、径向矢量及轴向矢量具有以下关系：
[0493]
[0494]
[0495]
[0496] 称式(115)为零位矢量的反转公式；称式(116)为零位矢量与径向矢量的互换公式；称式(117)为径向矢量不变性公式。由式(109)、式(115)至式(117)得
[0497]
[0498]
[0499] 由式(118)得
[0500]
[0501] 因 是 的对称部分的结构常数，故称式(118)为矢量 的对称分解式。因 是的反对称部分的结构常数，故称式(119)为矢量 的反对称分解式。称式(120)为归零等式。
[0502] 【3】运动链Dixon行列式性质
[0503] 定义
[0504]
[0505] 由式(52)得
[0506]
[0507] 其中：
[0508]
[0509] 考虑式(122)，若M为4·4的矩阵，则有
[0510]
[0511]
[0512]
[0513] 且有
[0514]
[0515] 由式(62)及式(110)得
[0516]
[0517] 由式(126)证得
[0518]
[0519] 由式(129)得
[0520]
[0521] 式(127)可以将 及 可以转化为关于 的多重线性型。同时，对yl及τl具有对称(轮换)性。由式(111)、式(118)及式(119)得
[0522]
[0523] 式(129)由三个导出的独立结构参量 及一个运动变量τl构成。由式(129)得
[0524]
[0525]
[0526] 由式(127)及式(131)得
[0527]
[0528] 由式(127)及式(132)得
[0529]
[0530] 基于类DCM的2R方向逆解
[0531] 给定6R轴链il6＝(i,1:6]，轴矢量3n4及4n5，期望第5轴的DCM为 期望第3轴的DCM为 方向矢量5l6与期望方向 对齐的逆解需要满足以下方程：
[0532]
[0533] 上述方程的建立步骤为：
[0534] 方向矢量5l6与期望方向 对齐需满足
[0535]
[0536] 由式(60)得
[0537]
[0538] 即
[0539]
[0540] 式(137)重新表示为式(134)。
[0541] 基于轴不变量的Cayley变换
[0542] 当给定角度 后，其正、余弦及其半角的正、余弦均是常数；为方便表达，记[0543]
[0544] 由式(138)得
[0545]
[0546] 定义
[0547]
[0548] 故有
[0549]
[0551] 一、基于轴不变量的3R机械臂位置逆解方法
[0552] 给定3R转动链 及期望姿态 轴不变量序列 求关节变量序列 这是3R姿态逆解问题。
[0553] 【1】根据机械臂n元3D矢量位姿方程，获得n个“n元2阶”多项式方程。
[0554] 由式(47)得3R运动学方程
[0555]
[0556] 由式(142)得
[0557]
[0558] 由式(143)得
[0559]
[0560] 若记
[0561]
[0562] 则由式(61)及式得式(145)
[0563]
[0564] 由式(144)及式(145)得
[0565]
[0566] 下面，阐述3R机械臂运动学方程的Dixon行列式的结构模型及特点。由式(147)得3R运动学多项式方程
[0567]
[0568] 多项式系统F3(Y2|T2)，根据双线性型行列式通式
[0569]
[0570] 则有
[0571]
[0572] 其中：
[0573]
[0574]
[0575] 由式(18)、式(147)及式(148)得
[0576]
[0577] 由式(22)及式(153)可知式(151)成立。由式(127)及式(145)得
[0578]
[0579]
[0580] 由式(145)、式(154)及式(155)得
[0581]
[0582] 其中：应用式(133)计算
[0583] 显然，式(156)中的y2阶次β2∈[0:3]及y3阶次β3∈[0:1]。考虑式(153)后三项：中的y2阶次β2∈[0:3]及y3阶次β3∈[0:1]；
中的y2阶次β2∈[0:2]及y3阶次β3∈[0:1]；
中的y2的阶次β2∈[0:3]及y3的阶次β3∈[0:
1]。由上可知：式(153)中的y2阶次β2∈[0:3]及y3的阶次β3∈[0:1]。故有S＝8。
中的y2阶次β2∈[0:2]及y3阶次β3∈[0:1]；
中的y2的阶次β2∈[0:3]及y3的阶次β3∈[0:
1]。由上可知：式(153)中的y2阶次β2∈[0:3]及y3的阶次β3∈[0:1]。故有S＝8。
[0584] 由式(145)、式(153)至式(156)可知： 中组合变量系数为独立的列向量，故选取的系数来构成方阵 剩余列向量一定与 的各列相关。故式(147)成立。
[0585] 【2】应用“基于轴不变量的Dixon行列式计算”方法及“分块矩阵的高维行列式计算”方法简化行列式计算。
[0586] 根据运动链Dixon行列式性质，由式(127)及式(145)得
[0587]
[0588]
[0589]
[0590] 分别表示轴2至轴3、轴3至轴3S的零位矢量、径向矢量及轴向矢量。
[0591] 由式(157)得
[0592]
[0593] 由式(158)得
[0594]
[0595] 由式(159)得
[0596]
[0597] 由式(153)得
[0598]
[0599] 将式(160)至式(162)代入式(163)得
[0600]
[0601] 【3】应用n个“n元N阶”多项式的Dixon消元与求解原理完成位姿逆解计算，其中：根据Dixon矩阵的行列式为0，得到一元高阶多项式方程，应用基于友阵的一元高阶多项式方程求解一元高阶多项式方程的解。
[0602] 一元n阶多项式p(x)＝a0+a1x+…an-1xn-1+xn具有n个解。若能找到一个矩阵A，满足|A-λl·1n|·vl＝0，其中：l∈[1:n]，λl为该矩阵的特征值，vl为对应的特征矢量。若矩阵A的特征方程为 则称该矩阵为多项式p(x)的友矩阵(Companion Matrix，简称友阵)，因此，多项式方程p(λl)＝0的解为其友阵A的特征方程|A-λl·1n|＝0的解。
[0603] 若多项式p(x)的友阵为
[0604]
[0605] 则由矩阵A的特征向量构成的矩阵为范德蒙德(Vandermonde)矩阵为
[0606]
[0607] 且有
[0608] p(λl)＝|A-λl·1n|＝0 (167)。
[0609] 由式(28)、式(151)及式(152)得
[0610]
[0611] 因S＝8，应用式(1)计算 的复杂度为8·8！＝322560；而应用式(5)进行二次分块的行列式计算，其中：2·2分块矩阵计算复杂度为4！(2·2！+2·2！+1)/(2！2！)＝30，4·4分块分矩阵计算复杂度为8！(30+30+1)/(4！4！)＝4270。一般情况下，式(168)是关于τ1的16阶单项式方程。
[0612] 二、基于轴不变量的通用6R机械臂位姿逆解方法
[0613] 设定有6个转动轴，拾取点位于第6轴轴线上，且第4轴与第5轴不共轴的机械臂为通用6R机械臂。通过前5轴控制第6轴与期望的位置及指向对齐，第6轴可以无限转动或控制第6轴满足径向对齐。
[0614] 基于轴不变量的通用6R机械臂位姿逆解：
[0615] 给定第6轴期望位置矢量 及期望姿态 的逆解问题与给定第6轴期望位置矢量及第5轴期望姿态 的逆解问题等价；通用6R机械臂本质上是5R轴链系统。
[0616] 若给定6R轴链il6＝(i,1:6]，il1＝03，第6轴期望位置矢量为 及第5轴期望姿态第3轴关节Ju-Gibbs规范四元数为 其他轴表达方式同理；则由轴不变量表征的6R机械臂运动学多项式方程为：
[0617]
[0618] 其中
[0619]
[0620] 式中，\为续行符； 分别表示轴4至轴5、轴5至轴6的零位矢量、径向矢量； 是轴不变量 的叉乘矩阵；03＝[0 0 0]T；
[0621] 系统结构参数及期望Ju-Gibbs姿态四元数构成的矩阵表示为
[0622]
[0623]
[0624]
[0625] 其中，表示系统结构参数的4×4矩阵； 表示取 的第一行元素，依次类推，表示取 的第k+1行元素；右上角标表达形式[ ]表示取行或列，表达形式[·]表示取所有列； 表示取3E5的第3行、第所有列。
[0626] 并进行如下定义：
[0627]
[0628] 下面给出上式的建立过程：
[0629] 期望 与第5轴期望姿态 对齐，若给定运动链iln，k∈iln，期望规范Ju-Gibbs四元数 及期望位置矢量 则位置及指向对齐关系表示为
[0630]
[0631] 由式(175)得 进而，由式(58)得
[0632]
[0633] 由式(52)得式(170)，由式(176)得
[0634]
[0635] 由式(53)得
[0636]
[0637] 其中：
[0638]
[0639] 由式(177)及式(178)得
[0640]
[0641] 由式(64)得
[0642]
[0643] 将式(181)代入式(177)得
[0644]
[0645] 其中：
[0646]
[0647] 由式(183)得式(172)。由式(172)及式(182)得
[0648]
[0649]
[0650] 式(184)是关于 期望姿态 及4轴、5轴结构参数的约束方程。由式(62)得
[0651]
[0652] 一方面，由式(184)、式(185)及式(186)得
[0653]
[0654] 另一方面，由式(176)、(184)及式(187)得
[0655]
[0656] 故得
[0657]
[0658] 由式(185)得(174)及
[0659]
[0660] 由式(184)及式(174)得
[0661]
[0662] 式(184)至式(190)用于后续方程简化，带有上下标的C是结构常数矩阵。考虑式(182)两边2范数得
[0663]
[0664] 考虑当il1＝03时的位置矢量对齐关系得
[0665]
[0666] 由式(47)及式(193)得
[0667]
[0668] 进而，得
[0669]
[0670] 即有
[0671]
[0672] 显然，有
[0673]
[0674] 由式(110)、式(176)及式(195)得式(194)之左式
[0675]
[0676] 结构参数 及 是链节l的结构参量，在系统零位时，它们可以通过外部测量得到。零位矢量、径向矢量及轴向矢量是与转动角无关的不变量。其中，零位矢量是特定的径向矢量。
[0677] 由式(110)、式(185)、式(190)及式(188)得
[0678]
[0679] 由式(190)及式(197)得
[0680]
[0681] 由式(184)、式(190)及式(198)得
[0682]
[0683] 将式(196)及式(199)代入(194)，且消去两边的 得式(169)。
[0684] 消去τ4及τ5后的位置方程(169)是3个“3元2阶”多项式方程，可以等同于3R机械臂问题，采用上述基于轴不变量的3R机械臂位置逆解方法进行求解，为实时计算通用6R轴机械臂的逆解奠定了基础。一方面，将有利于提高6R机械臂的绝对定位精度；另一方面，在结构上可以使传统解耦机械臂的第4轴及第5轴向根方向移动，不仅可以优化机械臂的结构，而且有利于提高6R机械臂避让障碍的灵活性。
[0685] 三、基于轴不变量的通用7R机械臂位姿逆解
[0686] 称有7个转动轴，拾取点位于第7轴轴线上，且第5轴与第6轴不共轴的机械臂为通用7R机械臂。通用7R机械臂通过前6轴控制第7轴与期望位置及姿态对齐，第7轴无限转动或控制第7轴满足径向对齐。因此，7R通用机械臂比6R机械臂具有更强的空间操作的灵活性。
[0687] 给定期望位置 及期望姿态 的逆解问题与给定期望位置 及期望姿态 的逆解问题等价，即7R机械臂本质上是6R轴链系统，可根据上述基于轴不变量的6R机械臂位置逆解方法进行求解。下面，先阐述通用7R机械臂逆解过程。
[0688] 若给定7R轴链il7＝(i,1:7]，il1＝03，期望位置矢量及Ju-Gibbs四元数分别记为及 则由轴不变量表征的7R机械臂运动学多项式方程为：
[0689]
[0690] 其中：
[0691]
[0692] 系统结构参数及期望Ju-Gibbs姿态四元数构成的矩阵表示为
[0693]
[0694]
[0695]
[0696]
[0697] 上式建立步骤为：
[0698] 期望 与姿态 对齐，由式(175)得 进而，由式(58)得
[0699]
[0700] 由式(52)得式(201)，由式(206)得
[0701]
[0702] 由式(53)得
[0703]
[0704] 其中：
[0705]
[0706] 由式(207)及式(208)得
[0707]
[0708] 由式(64)得
[0709]
[0710] 将式(211)代入式(210)得
[0711]
[0712] 其中：
[0713]
[0714] 由式(213)得式(203)。由式(203)及式(212)得
[0715]
[0716]
[0717] 式(214)是关于 姿态 及5轴、6轴结构参数的约束方程。由式(203)得(205)及[0718]
[0719] 由式(214)及式(205)得
[0720]
[0721] 式(214)至式(217)用于后续方程简化，C是结构常数矩阵。考虑式(216)两边2范数得
[0722]
[0723] 由式(62)得
[0724]
[0725] 一方面，由式(216)、式(215)及式(219)得
[0726]
[0727] 另一方面，由式(206)、(214)及式(220)得
[0728]
[0729] 故得
[0730]
[0731] 由式(222)得式(200)中的姿态方程。
[0732] 考虑当il1＝03时的位置矢量对齐关系得
[0733]
[0734] 由式(47)及式(223)得
[0735]
[0736] 进而，得
[0737]
[0738] 即有
[0739]
[0740] 显然，有
[0741]
[0742] 由式(206)及式(225)得式(224)之左式
[0743]
[0744] 由式(110)、式(215)、式(216)及式(221)得
[0745]
[0746] 由式(216)及式(227)得
[0747]
[0748] 由式(214)、式(217)及式(228)得
[0749]
[0750] 将式(226)及式(229)代入(224)，且消去两边的 得(200)中位置方程。
[0751] 消去τ5及τ6后的位姿方程(200)，是4个“4元2阶”多项式方程。
[0752] 1.Ju-Gibbs增量四元数及性质
[0753] 通用7R机械臂位姿逆解的计算复杂度较高，需要解决实时计算的技术问题。在工程上，计算精度是一个相对的概念，只需要保证数值计算精度高于系统结构参数精度的4至6倍就能满足工程要求。式(59)所示的姿态方程及式(47)所示的位置方程本质上是Ju-Gibbs四元数的表达式。只要保证式(59)及式(47)具有足够的计算精度就能满足工程精度要求。下面，首先提出“居-吉布斯”增量四元数(Delta-quaternion)，再建立通用7R机械臂增量位姿方程，并进行实时逆解计算。
[0754] 【1】“居-吉布斯”增量四元数定义
[0755] 定义“居-吉布斯”增量四元数
[0756]
[0757] 其中：
[0758]
[0759] 显然，“居-吉布斯”增量四元数是四维复数，且有
[0760]
[0761] 【2】“居-吉布斯”增量四元数性质
[0762] 由式(48)及式(230)得
[0763]
[0764] 由式(51)得
[0765]
[0766] 由式(234)得
[0767]
[0768] 由式(50)得
[0769]
[0770] 由式(236)及式(54)得
[0771]
[0772] 由式(62)得
[0773]
[0774] 2.基于轴不变量的通用7R机械臂运动规划
[0775] 通用7R机械臂逆解由于计算复杂度极高，在现有技术条件下无法实现。但是，通常位于第7轴上的拾取点与第6轴的距离很小。因此，在第6轴取距拾取点较近的点为名义拾取点，先计算通用6R机械臂的逆解；以之为基础，再应用数值迭代法，完成通用7R机械臂的运动规划与逆解计算。下面，探讨通用7R机械臂的增量(Delta)位姿方程建立及逆解的问题。
[0776] 【1】通用7R轴链的增量位姿方程
[0777] 首先陈述由Gibbs增量四元数表征的通用7R轴链的增量位姿方程，最后求其逆解。
[0778] 若给定6R轴链il7＝(i,1:7]，il1＝03，位置矢量及Ju-Gibbs增量四元数分别记为及 则由Ju-Gibbs增量四元数表征的通用6R机械臂增量位姿方程表示为：
[0779]
[0780]
[0781] 其中： 3n4为杆件3到杆件4的坐标矢量，其是轴不变量； 为轴不变量3n4的叉乘矩阵，其余杆件同理；，\为续行符；
分别表示轴5至轴6、轴6至轴7的零位矢量、径向矢量； 是轴不变量 的叉乘
矩阵；03＝[0 0 0]T；
分别表示轴5至轴6、轴6至轴7的零位矢量、径向矢量； 是轴不变量 的叉乘
矩阵；03＝[0 0 0]T；
[0782] 由行四元数构成的系统结构参数矩阵表示为
[0783]
[0784]
[0785]
[0786]
[0787]
[0788]
[0789] 下面，分析通用7R机械臂的增量位姿逆解。显然，式(239)是关于{εl|l∈[1:4}的线性方程。将式(239)重新表示为
[0790] A·[ε1 ε2 ε3 ε4]T＝b (247)
[0791] 若A-1存在，解式(245)得
[0792] [ε1 ε2 ε3 ε4]T＝A-1·b。 (248)
[0793] 由式(235)及式(248)得 由式(215)得ε5及ε6。至此，得到全部逆解。
[0794] 通用7R机械臂位姿逆解的作用在于：可以通过增量位置矢量 及Gibbs增量四元数 应用迭代逼近(Iterative Approximation)的算法使通用7R机械臂位姿对齐期望的位姿。
[0795] 【2】基于偏速度迭代的通用7R机械臂运动规划
[0796] 当k∈iln时，由式(3.304)至式(3.306)可知
[0797]
[0798]
[0799] 记运动链为 由式(249)及式(250)得
[0800]
[0801] 记 记期望位姿分别为 及 且有将式(251)合写为
[0802]
[0803] 由式(252)应用梯度(Gradient Descent Method/GDM)下降法得
[0804]
[0805] 其中：步长Step>0，Step→0。显然有
[0806]
[0807] 选择Step步长，由初态 开始迭代，直至终态 基于偏速度的迭代优化步骤如下：
[0808] (1)确定目标函数
[0809]
[0810] 显然，Goal表示 及 的方差。
[0811] (2)选择步长
[0812] 一方面，应用构造法得确定步长
[0813]
[0814] 则由式(253)及式(256)可知：仅当 时，Step→0，Goal→0。
[0815] 另一方面，由式(231)及式(252)得
[0816]
[0817] 其中：ε(i,6]＝[ε1 ε2 … ε6]T。由式(257)及式(253)确定步长
[0818]
[0819] (3)迭代过程
[0820] 一方面，若取式(258)之步长Step，由式(259)完成迭代计算：
[0821]
[0822] 其中：当 时，迭代过程结束。对于式(259)的迭代过程，则必有
[0823] δGoal≤0， (260)
[0824] 即式(259)的迭代过程一定收敛。
[0825] 迭代计算式的推导过程为：由式(252)及式(255)得
[0826]
[0827] 由式(253)、式(256)及式(261)得
[0828]
[0829] 因式(256)与式(258)在理论上等价，故可用式(258)替代式(256)。但是式(256)与式(258)的计算过程不同：因计算机字长有限，当 时，前者精度越来越差，而后者越来越高；同时，后者计算量相对较小。因此，在工程应用时，应用式(258)中的步长更好。
[0830] 当 时，得到稳态解φ(i,6]，即为通用机械臂的位姿逆解。
[0831] 基于偏速度的通用7R机械臂运动规划特点在于：通过迭代，逐步逼近期望的位姿，可以得到一条由初始位姿至期望位姿的路径。因为该方法是趋向目标的优化过程，所以实时性较差。若在迭代过程中，控制关节增量，则可以满足关节速度的约束；因此，在完成运动规划的同时，也获得对应期望位姿的一组逆解。
[0832] 以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 内啮合强制运动链立式滚磨机传动伸缩机构 | 2020-05-15 | 541 |
| 内啮合强制运动链立式滚磨机 | 2020-05-11 | 698 |
| 含带传动的变驱动运动链及含该运动链的并联机器人 | 2020-05-11 | 792 |
| 一种负反馈自适应蚁群算法的机构运动链同构识别方法 | 2020-05-18 | 726 |
| 一种基于免疫遗传混合算法的机构运动链同构识别方法 | 2020-05-19 | 908 |
| 基于编号矩阵的运动链同构判定方法 | 2020-05-13 | 885 |
| 一种基于拓扑特征环路码的含复铰运动链的同构识别方法 | 2020-05-18 | 986 |
| 用于电力推动车辆的运动链 | 2020-05-12 | 263 |
| 强制运动链伞齿轮变交角可偏置滚磨机滚磨轮夹具 | 2020-05-20 | 570 |
| 内啮合强制运动链滚磨机减速挂轮机构 | 2020-05-14 | 1031 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
分析报告
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
QQ群二维码
意见反馈
意见反馈