基于轴不变量的通用3R机械臂逆解建模与解算方法
阅读：838发布：2021-01-17
专利汇可以提供基于轴不变量的通用3R机械臂逆解建模与解算方法专利检索，专利查询，专利分析的服务。并且本 发明 公开了一种基于轴不变量的通用3R 机械臂 逆解建模与解算方法，应用n个n元N阶多项式的Dixon消元与求解原理，进行 位姿 逆解计算，根据机械臂n元3D矢量位姿方程，获得n个n元2阶多项式方程；应用基于轴不变量的Dixon行列式计算式及分 块 矩阵的行列式计算式简化行列式计算；应用n个n元N阶多项式的Dixon消元与求解原理完成位姿逆解计算，根据Dixon矩阵的行列式为0，得到一元高阶多项式方程，应用基于友阵的一元高阶多项式方程求解一元高阶多项式方程的解。本发明的方法可提高机械臂的绝对 定位 精度 ；与D‑H参数相比，求解过程具有通用性，可以获得系统全部逆解。，下面是基于轴不变量的通用3R机械臂逆解建模与解算方法专利的具体信息内容。
1.一种基于轴不变量的通用3R机械臂逆解建模与解算方法，其特征是，
应用n个“n元N阶”多项式的Dixon消元与求解原理，进行位姿逆解计算，主要包括以下步骤：
【1】根据机械臂n元3D矢量位姿方程，获得n个“n元2阶”多项式方程；
【2】应用基于轴不变量的Dixon行列式计算式、分块矩阵的行列式计算式或对行列式进行行阶梯化计算式简化行列式计算；
【3】应用n个“n元N阶”多项式的Dixon消元与求解原理完成位姿逆解计算，其中：根据Dixon矩阵的行列式为0，得到一元高阶多项式方程，应用基于友阵的一元高阶多项式方程求解一元高阶多项式方程的解；
对于任意杆件l, 定义与欧拉四元数同构的居-吉布斯即Ju-Gibbs规范四元数：
其中： 为Gibbs矢量；Gibbs共轭四元数为：
其中：
式中， 为居-吉布斯规范四元数 模的平方；表达形式幂符 表示 的x次幂；右上角角标∧或 表示分隔符；轴不变量 是轴不变量 的叉乘矩阵； 是Gibbs矢量 的叉乘矩阵；若用 表示属性占位，则式中的表达形式 表示成员访问符；
步骤【1】中，
对于轴链 有
建立规范的姿态方程为：
建立规范的定位方程：
式中，k, 为任意杆件，表达形式 表示 的x次幂；右上角角标∧或表示分隔符； 是轴不变量 的叉乘矩阵，杆件l, 为杆件k, 时同理替换；1为三维单位矩阵；iQn表示姿态；
为沿矢量轴 的线位置； 为零位时由原点 至原点Ol的平动矢量；为投影符， 为在大地坐标系的投影矢量。
2.根据权利要求1所述的基于轴不变量的通用3R机械臂逆解建模与解算方法，其特征是，
步骤【2】中，基于轴不变量的Dixon行列式计算式为：
根据运动链Dixon行列式性质有：
并记：
式中， 为旋转变换矩阵； 表示用辅助变量yl的前l个依次替换原变量τl中的l个变量，记“|”为替换操作符；
式(80)将 及 转化为关于 的多重线性型；同时， 对yl及τl具有
对称性；
由式(47)得3R运动学方程
由式(90)得
由(91)式得
记
则由式(51)及式得(93)
由式(92)及式(93)得
由式(95)得3R运动学多项式方程
多项式系统F3(Y2|T2)，根据双线性型行列式通式
则有
其中：
中组合变量系数为独立的列向量，故选取 的系数来构成方阵 剩余列向量一定与 的各列相关；
由式(80)及式(93)得
分别表示轴2至轴3、轴3至轴3S的零位矢量、径向矢量及轴
向矢量；
得简化的3元N阶Dixon行列式为
式中， 为大小为S
×S的Dixon矩阵，其第[i][j]成员为单变量τ1的N阶多项式。
3.根据权利要求2所述的基于轴不变量的通用3R机械臂逆解建模与解算方法，其特征是，
步骤【2】中，分块矩阵的行列式计算式为：
若记大小为(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行列式与分块矩阵 及 的行列式关系为
4.根据权利要求2所述的基于轴不变量的通用3R机械臂逆解建模与解算方法，其特征是，
步骤【2】中，对行列式进行行阶梯化计算原理：
对于S×S矩阵，其每一项是关于τ1的n阶多项式；计算该矩阵的行列式时，可通过初等行变换将原行列式变为上三角行列式，再将非零的对角线元素相乘，得到行列式的多项式表达式；该式为0，得到τ1的所有解；
行阶梯化的具体方法为，先对行列式第一列的最高阶次由高到低进行排序，再进行最多(S-1)×n次初等行变换消元，得到第一列只有第一个元素不为0的行列式；再对该行列式第1行及1列的余子式进行初等行变换消元，依次迭代求解。
5.根据权利要求1所述的基于轴不变量的通用3R机械臂逆解建模与解算方法，其特征是，
步骤【3】中，n个“n元N阶”多项式系统的Dixon多项式构建步骤为：
引入辅助变量[y2,y3,…,yn]，且有
对于多元多重多项式 用辅助变量Ym的前m个依次替换原变量Xn中的m个变量，记“|”为替换操作符，得到增广的多项式
得
其中：
定义可分离组合变量 及 如下：
由式(14)及式(15)知：替换式 是关于 及 的双重线性型；相应地，用辅助变量替换的多项式系统记为
给定n个“n元N阶”多项式系统 定义其Dixon多项式为
由式(17)得
考虑式(13)及式(18)得该多项式的Dixon行列式
在笛卡尔空间下，由位置矢量或转动矢量构成的行列式表示矢量张成空间的容积(Volume)；在不同笛卡尔空间下具有容积的不变性；其中：
给定n个“n元N阶”多项式系统Fn(Yn-1|Xn-1)，n≥2；存在与消去变量x2,…,xn无关的Dixon矩阵SΘS(x1)，其Dixon多项式 表示为分离变量 及 的双重线
性型：
α[l]∈[0,N·(n-l+1)-1],l∈[2:n]； (23)
为大小为S×S的Dixon矩阵，其第[i][j]成员为单变量x1的N阶多项式：
其中：
考虑式(22)，若 故得
Det(SΘS(x1))＝0； (28)
称式(28)中“n个n元”为Dixon消元的必要条件，从而获得可行解。
6.根据权利要求2所述的基于轴不变量的通用3R机械臂逆解建模与解算方法，其特征是，
由式(28)、式(99)及式(100)得
式(116)是关于τ1的16阶单项式方程，应用式(5)进行二次分块的行列式计算。
基于轴不变量的通用3R机械臂逆解建模与解算方法
技术领域
背景技术
[0002] 自主机器人研究的一个重要方面是需要解决变拓扑结构机器人的运动学建模问题。在MAS中，具有动态的图结构(Dynamic Graph Structure)，可以动态地建立基于运动轴的有向Span树，为研究可变拓扑结构(Variable Topology Structure)的机器人建模与控制奠定了基础。为此，需要提出基于轴不变量的通用机械臂逆解原理，既要建立包含坐标系、极性、结构参数、关节变量的完全参数化的正运动学模型，又要实时地计算位姿方程；一方面，可以提高机器人的自主性，另一方面，可以提高机器人位姿控制的绝对精度。
[0003] 3R机械臂位置逆解是指：给定3R机械臂结构参数及期望位置，计算3个关节变量，使腕心位置与期望位置对齐。现有的基于D-H参数的3R机械臂位置逆解方法存在如下缺点：建立D-H系及D-H参数的过程不自然，应用繁琐；需要处理由计算方法导致的奇异性问题；在应用时，易引入系统测量误差。基于D-H参数的3R机械臂逆解原理不具有普适性，难以推广来解决通用6R机械臂的逆解问题。
发明内容
[0004] 本发明所要解决的技术问题是提供一种基于轴不变量的通用3R机械臂逆解建模与解算方法，可提高机械臂的绝对定位精度；与D-H参数相比，求解过程具有通用性，可以获得系统全部逆解。
[0005] 为解决上述技术问题，本发明采用以下技术方案：
[0006] 一种基于轴不变量的通用3R机械臂逆解建模与解算方法，其特征是，
[0007] 应用n个“n元N阶”多项式的Dixon消元与求解原理，进行位姿逆解计算，主要包括以下步骤：
[0008] 【1】根据机械臂n元3D矢量位姿方程，获得n个“n元2阶”多项式方程；
[0009] 【2】应用基于轴不变量的Dixon行列式计算式、分块矩阵的行列式计算式或对行列式进行行阶梯化计算式简化行列式计算；
[0010] 【3】应用n个“n元N阶”多项式的Dixon消元与求解原理完成位姿逆解计算，其中：根据Dixon矩阵的行列式为0，得到一元高阶多项式方程，应用基于友阵的一元高阶多项式方程求解一元高阶多项式方程的解。
[0011] 对于任意杆件 定义与欧拉四元数同构的居-吉布斯即Ju-Gibbs规范四元数：
[0012]
[0013] 其中： 为Gibbs矢量；Gibbs共轭四元数为：
[0014]
[0015] 其中：
[0016]
[0017] 式中， 为居-吉布斯规范四元数 模的平方；表达形式幂符 表示 的x次幂；右上角角标∧或表示分隔符；轴不变量 是轴不变量 的叉乘矩阵； 是Gibbs矢量的叉乘矩阵；若用 表示属性占位，则式中的表达形式 表示成员访问符。
[0018] 步骤【1】中，
[0019] 对于轴链 有
[0020]
[0021] 建立规范的姿态方程为：
[0022]
[0023] 建立规范的定位方程：
[0024]
[0025] 式中， 为任意杆件，表达形式 表示 的x次幂；右上角角标∧或 表示分隔符；是轴不变量 的叉乘矩阵，杆件l, 为杆件k, 时同理替换；1为三维单位矩阵；iQn表示姿态； 为沿矢量轴 的线位置； 为零位时由原点 至原点Ol的平动矢量；为投影符，为 在大地坐标系的投影矢量。
[0026] 步骤【2】中，基于轴不变量的Dixon行列式计算式为：
[0027] 根据运动链Dixon行列式性质有：
[0028]
[0029] 并记：
[0030]
[0031]
[0032] 式中， 为旋转变换矩阵； 表示用辅助变量yl的前l个依次替换原变量τl中的l个变量，记“|”为替换操作符；
[0033] 式(80)将 及 转化为关于 的多重线性型；同时， 对yl及τl具有对称性；
[0034] 由式(47)得3R运动学方程
[0035]
[0036] 由式(90)得
[0037]
[0038] 由(91)式得
[0039]
[0040] 记
[0041]
[0042] 则由式(51)及式得(93)
[0043]
[0044] 由式(92)及式(93)得
[0045]
[0046] 由式(95)得3R运动学多项式方程
[0047]
[0048] 多项式系统F3(Y2|T2)，根据双线性型行列式通式
[0049]
[0050] 则有
[0051]
[0052] 其中：
[0053]
[0054]
[0055] 中组合变量系数为独立的列向量，故选取 的系数来构成方阵 剩余列向量一定与 的各列相关；
[0056] 由式(80)及式(93)得
[0057]
[0058]
[0059]
[0060] 分别表示轴2至轴3、轴3至轴3S的零位矢量、径向矢量及轴向矢量；
[0061] 得简化的3元N阶Dixon行列式为
[0062]
[0063] 式中， 为大小为S×S的Dixon矩阵，其第[i][j]成员为单变量τ1的N阶多项式。
[0064] 步骤【2】中，分块矩阵的行列式计算式为：
[0065] 若记大小为(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行列式与分块矩阵 及 的行列式关系为
[0066]
[0067] 步骤【2】中，对行列式进行行阶梯化计算原理：
[0068] 对于S×S矩阵，其每一项是关于τ1的n阶多项式；计算该矩阵的行列式时，可通过初等行变换将原行列式变为上三角行列式，再将非零的对角线元素相乘，得到行列式的多项式表达式；该式为0，得到t1的所有解；
[0069] 行阶梯化的具体方法为，先对行列式第一列的最高阶次由高到低进行排序，再进行最多(S-1)×n次初等行变换消元，得到第一列只有第一个元素不为0的行列式；再对该行列式第1行及1列的余子式进行初等行变换消元，依次迭代求解。
[0070] 步骤【3】中，n个“n元N阶”多项式系统的Dixon多项式构建步骤为：
[0071] 引入辅助变量[y2,y3,…,yn]，且有
[0072]
[0073] 对于多元多重多项式 用辅助变量Ym的前m个依次替换原变量Xn中的m个变量，记“|”为替换操作符，得到增广的多项式
[0074]
[0075]
[0076] 得
[0077]
[0078] 其中：
[0079]
[0080] 定义可分离组合变量 及 如下：
[0081]
[0082] 由式(14)及式(15)知：替换式 是关于 及 的双重线性型；相应地，用辅助变量替换的多项式系统记为
[0083]
[0084] 给定n个“n元N阶”多项式系统 定义其Dixon多项式为
[0085]
[0086] 由式(17)得
[0087]
[0088] 考虑式(13)及式(18)得该多项式的Dixon行列式
[0089]
[0090] 在笛卡尔空间下，由位置矢量或转动矢量构成的行列式表示矢量张成空间的容积(Volume)；在不同笛卡尔空间下具有容积的不变性。其中：
[0091]
[0092] 给定n个“n元N阶”多项式系统Fn(Yn-1|Xn-1)，n≥2；存在与消去变量x2,…,xn无关的Dixon矩阵SΘS(x1)，其Dixon多项式 表示为分离变量 及 的双重线性型：
[0093]
[0094] α[l]∈[0,N·(n-l+1)-1],l∈[2:n]； (23)
[0095] 为大小为S×S的Dixon矩阵，其第[i][j]成员为单变量x1的N阶多项式：
[0096]
[0097] 其中：
[0098]
[0099] 考虑式(22)，若 故得
[0100] Det(SΘS(x1))＝0； (28)
[0101] 称式(28)中“n个n元”为Dixon消元的必要条件，从而获得可行解。
[0102] 由式(28)、式(99)及式(100)得
[0103]
[0104] 式(116)是关于τ1的16阶单项式方程，应用式(5)进行二次分块的行列式计算。
[0105] 本发明所达到的有益效果：
[0106] 本发明的方法提出了基于轴不变量的通用3R姿态逆解方法。特征在于：
[0107] 具有简洁、优雅的运动链符号系统，具有伪代码的功能，具有迭代式结构，保证系统实现的可靠性及机械化演算。
[0108] 具有基于轴不变量的迭代式，保证计算的实时性；实现坐标系、极性及系统结构参量的完全参数化，基于轴不变量的可逆解运动学具有统一的表达及简洁的结构化层次模型，保证位姿分析逆解(Analytical Inverse Solution to Position and Attitude)的通用性。
[0110] 由于轴不变量可以精确测量，有助于提高机械臂的绝对定位精度；由于关节变量的范围覆盖完整的一周，消除了D-H计算原理导致的奇异性；与D-H参数相比，求解过程具有通用性，可以获得系统全部逆解。附图说明
[0111] 图1自然坐标系与轴链；
[0112] 图2固定轴不变量；
[0113] 图3为定轴转动示意图；
[0114] 图4为轴不变量的导出不变量。
具体实施方式
[0115] 下面对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。
[0116] 定义1自然坐标轴：称与运动轴或测量轴共轴的，具有固定原点的单位参考轴为自然坐标轴，亦称为自然参考轴。
[0117] 定义2自然坐标系：如图1所示，若多轴系统D处于零位，所有笛卡尔体坐标系方向一致，且体坐标系原点位于运动轴的轴线上，则该坐标系统为自然坐标系统，简称自然坐标系。
[0118] 自然坐标系优点在于：(1)坐标系统易确定；(2)零位时的关节变量为零；(3)零位时的系统姿态一致；(4)不易引入测量累积误差。
[0119] 由定义2可知，在系统处于零位时，所有杆件的自然坐标系与底座或世界系的方向一致。系统处于零位即 时，自然坐标系 绕轴矢量 转动角度 将 转至F[l]； 在 下的坐标矢量与 在F[l]下的坐标矢量 恒等，即有
[0120]
[0121] 由上式知， 或 不依赖于相邻的坐标系 及F[l]；故称 或 为轴不变量。在不强调不变性时，可以称之为坐标轴矢量(简称轴矢量)。 或 表征的是体 与体l共有的参考单位坐标矢量，与参考点 及Ol无关。体 与体l即为杆件或轴。
[0122] 轴不变量与坐标轴具有本质区别：
[0123] (1)坐标轴是具有零位及单位刻度的参考方向，可以描述沿该方向平动的位置，但不能完整描述绕该方向的转动角度，因为坐标轴自身不具有径向参考方向，即不存在表征转动的零位。在实际应用时，需要补充该轴的径向参考。例如：在笛卡尔系F[l]中，绕lx转动，需以ly或lz为参考零位。坐标轴自身是1D的，3个正交的1D参考轴构成3D的笛卡尔标架。
[0124] (2)轴不变量是3D的空间单位参考轴，其自身就是一个标架。其自身具有径向参考轴，即参考零位。空间坐标轴及其自身的径向参考轴可以确定笛卡尔标架。空间坐标轴可以反映运动轴及测量轴的三个基本参考属性。
[0125] 已有文献将无链指标的轴矢量记为 并称之为欧拉轴(Euler Axis)，相应的关节角称为欧拉角(Euler Angle)。本申请之所以不再沿用欧拉轴，而称之为轴不变量，是因为轴不变量具有以下属性：
[0126] 【1】给定旋转变换阵 因其是实矩阵，其模是单位的，故其有一个实特征值λ1及两个互为共轭的复特征值λ2＝eiφ及λ3＝e-iφ；其中：i为纯虚数。因此，|λ1|·||λ2||·||λ3||＝1，得λ1＝1。轴矢量 是实特征值λ1＝1对应的特征矢量，是不变量；
[0127] 【2】是3D参考轴，不仅具有轴向参考方向，而且具有径向参考零位，将在3.3.1节予以阐述。
[0128] 【3】在自然坐标系下： 即轴不变量 是非常特殊的矢量，它对时间的导数也具有不变性，且有非常优良的数学操作性能；
[0129] 对轴不变量而言，其绝对导数就是其相对导数。因轴不变量是具有不变性的自然参考轴，故其绝对导数恒为零矢量。因此，轴不变量具有对时间微分的不变性。有：
[0130]
[0131] 【4】在自然坐标系统中，通过轴矢量 及关节变量 可以直接描述旋转坐标阵没有必要为除根之外的杆件建立各自的体系。同时，以唯一需要定义的根坐标系为参考，可以提高系统结构参数的测量精度；
[0132] 【5】应用轴矢量 的优良操作，将建立包含拓扑结构、坐标系、极性、结构参量及力学参量的完全参数化的统一的多轴系统运动学及动力学模型。
[0133] 因基矢量el是与F[l]固结的任一矢量，基矢量 是与 固结的任一矢量，又 是F[l]及 共有的单位矢量，故 是F[l]及 共有的基矢量。因此，轴不变量 是F[l]及 共有的参考基。轴不变量是参数化的自然坐标基，是多轴系统的基元。固定轴不变量的平动与转动与其固结的坐标系的平动与转动等价。
[0134] 在系统处于零位时，以自然坐标系为参考，测量得到坐标轴矢量 在运动副运动时，轴矢量 是不变量；轴矢量 及关节变量 唯一确定运动副 的转动关系。
[0135] 因此，应用自然坐标系统，当系统处于零位时，只需确定一个公共的参考系，而不必为系统中每一杆件确定各自的体坐标系，因为它们由轴不变量及自然坐标唯一确定。当进行系统分析时，除底座系外，与杆件固结的其它自然坐标系只发生在概念上，而与实际的测量无关。自然坐标系统对于多轴系统(MAS)理论分析及工程作用在于：
[0136] (1)系统的结构参数测量需要以统一的参考系测量；否则，不仅工程测量过程烦琐，而且引入不同的体系会引入更大的测量误差。
[0137] (2)应用自然坐标系统，除根杆件外，其它杆件的自然坐标系统由结构参量及关节变量自然确定，有助于MAS系统的运动学与动力学分析。
[0138] (3)在工程上，可以应用激光跟踪仪等光学测量设备，实现对固定轴不变量的精确测量。
[0139] (4)由于运动副R及P、螺旋副H、接触副O是圆柱副C的特例，可以应用圆柱副简化MAS运动学及动力学分析。
[0140] 定义3不变量：称不依赖于一组坐标系进行度量的量为不变量。
[0141] 定义4转动坐标矢量：绕坐标轴矢量 转动到角位置 的坐标矢量 为
[0142]
[0143] 定义5平动坐标矢量：沿坐标轴矢量 平动到线位置 的坐标矢量 为
[0144]
[0145] 定义6自然坐标：以自然坐标轴矢量为参考方向，相对系统零位的角位置或线位置，记为ql，称为自然坐标；称与自然坐标一一映射的量为关节变量；其中：
[0146]
[0147] 定义7机械零位：对于运动副 在初始时刻t0时，关节绝对编码器的零位 不一定为零，该零位称为机械零位；
[0148] 故关节 的控制量 为
[0149]
[0150] 定义8自然运动矢量：将由自然坐标轴矢量 及自然坐标ql确定的矢量 称为自然运动矢量。其中：
[0151]
[0152] 自然运动矢量实现了轴平动与转动的统一表达。将由自然坐标轴矢量及关节确定的矢量，例如 称为自由运动矢量，亦称为自由螺旋。显然，轴矢量 是特定的自由螺旋。
[0153] 定义9关节空间：以关节自然坐标ql表示的空间称为关节空间。
[0154] 定义10位形空间：称表达位置及姿态(简称位姿)的笛卡尔空间为位形空间，是双矢量空间或6D空间。
[0155] 定义11自然关节空间：以自然坐标系为参考，通过关节变量 表示，在系统零位时必有 的关节空间，称为自然关节空间。
[0156] 如图2所示，给定链节 原点Ol受位置矢量 约束的轴矢量 为固定轴矢量，记为 其中：
[0157]
[0158] 轴矢量 是关节自然坐标的自然参考轴。因 是轴不变量，故称 为固定轴不变量，它表征了运动副 的结构关系，即确定了自然坐标轴。固定轴不变量 是链节 结构参数的自然描述。
[0159] 定义12自然坐标轴空间：以固定轴不变量作为自然参考轴，以对应的自然坐标表示的空间称为自然坐标轴空间，简称自然轴空间。它是具有1个自由度的3D空间。
[0160] 如图2所示， 及 不因杆件Ωl的运动而改变，是不变的结构参考量。 确定了轴l相对于轴 的五个结构参数；与关节变量ql一起，完整地表达了杆件Ωl的6D位形。给定时，杆件固结的自然坐标系可由结构参数 及关节变量 唯一确定。称轴不变量 固定轴不变量 关节变量 及 为自然不变量。显然，由固定轴不变量及关节变量 构成的关节自然不变量 与由坐标系 至F[l]确定的空间位形 具
有一一映射关系，即
有一一映射关系，即
[0161]
[0162] 给定多轴系统D＝{T,A,B,K,F,NT}，在系统零位时，只要建立底座系或惯性系，以及各轴上的参考点Ol，其它杆件坐标系也自然确定。本质上，只需要确定底座系或惯性系。
[0163] 给定一个由运动副连接的具有闭链的结构简图，可以选定回路中任一个运动副，将组成该运动副的定子与动子分割开来；从而，获得一个无回路的树型结构，称之为Span树。T表示带方向的span树，以描述树链运动的拓扑关系。
[0164] I为结构参数；A为轴序列，F为杆件参考系序列，B为杆件体序列，K为运动副类型序列，NT为约束轴的序列即非树。 为取轴序列 的成员。转动副R，棱柱副P，螺旋副H，接触副O是圆柱副C的特例。
[0165] 描述运动链的基本拓扑符号及操作是构成运动链拓扑符号系统的基础，定义如下：
[0166] 【1】运动链由偏序集合(]标识。
[0167] 【2】A[l]为取轴序列A的成员；因轴名l具有唯一的编号对应于A[l]的序号，故A[l]计算复杂度为O(1)。
[0168] 【3】为取轴l的父轴；轴 的计算复杂度为O(1)。计算复杂度O()表示计算过程的操作次数，通常指浮点乘与加的次数。以浮点乘与加的次数表达计算复杂度非常烦琐，故常采用算法循环过程中的主要操作次数；比如：关节位姿、速度、加速度等操作的次数。
[0169] 【4】 为取轴序列 的成员； 计算复杂度为O(1)。
[0170] 【5】llk为取由轴l至轴k的运动链，输出表示为 且 基数记为|llk|。llk执行过程：执行 若 则执行 否则，结束。llk计算复杂度为O(|llk|)。
[0171] 【6】ll为取轴l的子。该操作表示在 中找到成员l的地址k；从而，获得轴l的子A[k]。因 不具有偏序结构，故ll的计算复杂度为
[0172] 【7】lL表示获得由轴l及其子树构成的闭子树，lL为不含l的子树；递归执行ll，计算复杂度为
[0173] 【8】支路、子树及非树弧的增加与删除操作也是必要的组成部分；从而，通过动态Span树及动态图描述可变拓扑结构。在支路llk中，若 则记即 表示在支路中取成员m的子。
[0174] 定义以下表达式或表达形式：
[0175] 轴与杆件具有一一对应性；轴间的属性量 及杆件间的属性量 具有偏序性。
[0176] 约定： 表示属性占位；若属性p或P是关于位置的，则 应理解为坐标系 的原点至F[l]的原点；若属性p或P是关于方向的，则 应理解为坐标系 至F[l]。
[0177] 及 应分别理解为关于时间t的函数 及 且 及 是t0时刻的常数或常数阵列。但是正体的 及 应视为常数或常数阵列。
[0178] 本申请中约定：在运动链符号演算系统中，具有偏序的属性变量或常量，在名称上包含表示偏序的指标；要么包含左上角及右下角指标，要么包含右上角及右下角指标；它们的方向总是由左上角指标至右下角指标，或由右上角指标至右下角指标，本申请中为叙述简便，有时省略方向的描述，即使省略，本领域技术人员通过符号表达式也可以知道，本申请中采用的各参数，对于某种属性符，它们的方向总是由偏序指标的左上角指标至右下角指标，或由右上角指标至右下角指标。例如： 可简述为(表示由k至l)平动矢量； 表示(由k至l的)线位置； 表示(由k至l的)平动矢量；其中：r表示“平动”属性符，其余属性符对应为：属性符φ表示“转动”；属性符Q表示“旋转变换矩阵”；属性符l表示“运动链”；属性符u表示“单位矢量”；属性符ω表示“角速度”；角标为i表示惯性坐标系或大地坐标系；其他角标可以为其他字母，也可以为数字。
[0179] 本申请的符号规范与约定是根据运动链的偏序性、链节是运动链的基本单位这两个原则确定的，反映了运动链的本质特征。链指标表示的是连接关系，右上指标表征参考系。采用这种符号表达简洁、准确，便于交流与书面表达。同时，它们是结构化的符号系统，包含了组成各属性量的要素及关系，便于计算机处理，为计算机自动建模奠定基础。指标的含义需要通过属性符的背景即上下文进行理解；比如：若属性符是平动类型的，则左上角指标表示坐标系的原点及方向；若属性符是转动类型的，则左上角指标表示坐标系的方向。
[0180] (1)lS－杆件l中的点S；而S表示空间中的一点S。
[0181] (2) －杆件k的原点Ok至杆件l的原点Ol的平动矢量；
[0182] 在自然坐标系F[k]下的坐标矢量，即由k至l的坐标矢量；
[0183] (3) －原点Ok至点lS的平动矢量；
[0184] 在F[k]下的坐标矢量；
[0185] (4) －原点Ok至点S的平动矢量；
[0186] 在F[k]下的坐标矢量；
[0188] －运动副 的轴矢量；
[0189] 及 分别在 及F[l]下的坐标矢量； 是轴不变量，为一结构常数；
[0190] 为转动矢量，转动矢量/角矢量 是自由矢量，即该矢量可自由平移；
[0191] (6) －沿轴 的线位置(平动位置)，
[0192] －绕轴 的角位置，即关节角、关节变量，为标量；
[0193] (7)左下角指标为0时，表示机械零位；如：
[0194] －平动轴 的机械零位，
[0195] －转动轴 的机械零位；
[0196] (8)0－三维零矩阵；1－三维单位矩阵；
[0197] (9)约定：“\”表示续行符； 表示属性占位；则
[0198] 幂符 表示 的x次幂；右上角角标∧或表示分隔符；如： 或 为 的x次幂。
[0199] 表示 的转置，表示对集合转置，不对成员执行转置；如：
[0200] 为投影符，表示矢量或二阶张量对参考基的投影矢量或投影序列，即坐标矢量或坐标阵列，投影即是点积运算“·”；如：位置矢量 在坐标系F[k]中的投影矢量记为[0201] 为叉乘符；如： 是轴不变量 的叉乘矩阵；给定任一矢量 的叉乘矩阵为叉乘矩阵是二阶张量。
[0202] 叉乘符运算的优先级高于投影符 的优先级。投影符 的优先级高于成员访问符或 成员访问符 优先级高于幂符
[0203] (10)单位矢量在大地坐标系的投影矢量 单位零位矢量
[0204] (11) －零位时由原点 至原点Ol的平动矢量，且记 表示位置结构参数。
[0205] (12)iQl，相对绝对空间的旋转变换阵；
[0206] (13)以自然坐标轴矢量为参考方向，相对系统零位的角位置或线位置，记为ql，称为自然坐标；关节变量 自然关节坐标为φl；
[0207] (14)对于一给定有序的集合r＝[1,4,3,2]T，记r[x]表示取集合r的第x行元素。常记[x]、[y]、[z]及[w]表示取第1、2、3及4列元素。
[0208] (15)ilj表示由i到j的运动链；llk为取由轴l至轴k的运动链；
[0209] 给定运动链 若n表示笛卡尔直角系，则称 为笛卡尔轴链；若n表示自然参考轴，则称 为自然轴链。
[0210] (16)Rodrigues四元数表达形式：
[0211] 欧拉四元数表达形式：
[0212] 不变量的四元数(也称为轴四元数)表达形式
[0213] 如位置矢量 在笛卡尔三个坐标轴上的投影矢量为 定义 由于左上角指标指明了参考系， 既间接表示了位移矢量 又直接表示了位移坐标矢量，即具有矢量及坐标矢量的双重作用。
[0214] 分块矩阵的高维行列式计算：
[0215] 记<1:n>表示自然数[1:n]的全排列，共有n！个实例。给定属于数域的大小为n×n的矩阵M，其j行i列元素记为 根据行列式定义得
[0216]
[0217] 其中：I[i1,…in]表示排列的逆序个数。式(1)计算复杂度为：n！次n个数积及n！次加法，具有指数计算复杂度，只能适用于维度较小的行列式。对于维度较大的行列式，通常应用Laplace公式进行递规运算，记 为 的伴随矩阵(Adjugate Matrix)，则有
[0218]
[0219] 更简单的算法通常应用高斯消去法或LU分解法，先通过初等变换将矩阵变为三角阵或三角阵的乘积，后计算行列式。上述针对数域的行列式计算方法不适用于高维度的多项式矩阵，需要引入分块矩阵的行列式计算方法。计算矢量多项式(Vector Polynomial)的行列式是一个特定的分块矩阵行列式的计算问题，它在矢量层次上表达了矢量与行列式的内在联系。而分块矩阵行列式计算则从矩阵层次上表达分块矩阵与行列式的内在规律。
[0220] 若给定矢量多项式 其中： 及为3D坐标矢量， 为多项式变量序列；若约定
[0221]
[0222] 则有
[0223]
[0224] 式(3)及式(4)可以推广至n维空间。
[0225] 实施例1
[0226] 给定2个2维行矢量多项式 及一方面，由式(4)得
[0227]
[0228] 另一方面，
[0229]
[0230] 上面的结果验证了式(4)的正确性。
[0231] 给出分块矩阵的行列式计算定理：
[0232] 若记大小为(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行列式与分块矩阵 及 的行列式关系为
[0233]
[0234] 对行列式进行行阶梯化计算原理：
[0235] 对于S×S矩阵，其每一项是关于t1的n阶多项式。计算该矩阵的行列式时，可通过初等行变换将原行列式变为上三角行列式，再将非零的对角线元素相乘，得到行列式的多项式表达式。该式为0，得到t1的所有解。
[0236] 行阶梯化的具体方法为，先对行列式第一列的最高阶次由高到低进行排序，再进行最多(S-1)×n次初等行变换消元，得到第一列只有第一个元素不为0的行列式。再对该行列式第1行及1列的余子式进行初等行变换消元，依次迭代求解。
[0237] 实施例2
[0238] 通过矩阵的初等行变换，得到 的行阶梯矩阵。
[0239] 步骤为：rk代表第k行。得
[0240]
[0241] 则得
[0242] 基于“N进位字”的N阶多项式系统：
[0243] 若n个“n元1阶”多项式幂积 中独立变量重复出现N次，则得到n个“n元N阶”多项式系统 “n元N阶多项式系统”与“n位N进位字” 同构。
[0244]
[0245]
[0246]
[0247]
[0248] n个“n元N阶”多项式系统的Dixon多项式：
[0249] 引入辅助变量[y2,y3,…,yn]，且有
[0250]
[0251] 在多元多重多项式(8)中，用辅助变量Ym的前m个依次替换原变量(Original Variables)Xn中的m个变量，记“|”为替换操作符，得到增广的(Extended)多项式[0252]
[0253]
[0254] 式中右上角标α、α表示幂；
[0255] 由式(6)及式(12)得
[0256]
[0257] 其中：
[0258]
[0259] 定义可分离组合变量 及 如下：
[0260]
[0261] 由式(14)及式(15)可知：替换式 是关于 及 的双重线性型。相应地，用辅助变量替换的多项式系统记为
[0262]
[0263] 给定n个“n元N阶”多项式系统 定义其Dixon多项式为
[0264]
[0265] 由式(17)得
[0266]
[0267] 式(15)中分离变量与文献不同：原变量Xn-1被辅助变量Yn-1替换的次序不同，Dixon多项式也不同。考虑式(13)及式(18)得该多项式的Dixon行列式
[0268]
[0269] 在笛卡尔空间下，由位置矢量或转动矢量构成的行列式表示矢量张成空间的容积(Volume)；在不同笛卡尔空间下具有容积的不变性。其中：
[0270]
[0271] n个“n元N阶”多项式的Dixon行列式的阶次及替换变量项数分别为：
[0272]
[0273] n个“n元N阶”Dixon矩阵：
[0274] 给定n个“n元N阶”多项式系统Fn(Yn-1|Xn-1)，n≥2；存在与消去变量x2,…,xn无关的Dixon矩阵SΘS(x1)，其Dixon多项式 表示为分离变量 及 的双重线性型：
[0275]
[0276] α[l]∈[0,N·(n-l+1)-1],l∈[2:n] (23)
[0277] 为大小为S×S的Dixon矩阵，其第[i][j]成员为单变量x1的N阶多项式：
[0278]
[0279] 其中：
[0280]
[0281] 若
[0282]
[0283] 则有
[0284]
[0285] 考虑式(22)，若 故得
[0286] Det(SΘS(x1))＝0。 (28)
[0287] 称式(28)中“n个n元”为Dixon消元的必要条件，从而获得可行解。若SΘS存在零行或零列向量，则无法建立x1的多项式方程；此时，通过除标量积之外的初等变换，将SΘS变为行阶梯(Row Echelon)矩阵Ech(SΘS)；在计算该矩阵的杻轴(Pivot)的积之后得方阵即在SΘS中选取S′个独立的列向量。
[0288] 任一个n个“n元N阶”多项式系统 的实例(简称多项式)记为其中： 且有 根据
的多项式确定Dixon矩阵、分离变量 及 选取 及 满足
的多项式确定Dixon矩阵、分离变量 及 选取 及 满足
[0289]
[0290] 确定双线性型
[0291]
[0292]
[0293] 其中： 中与 对应的各列线性独立。因 由式(22)及式(25)得
[0294]
[0295] 称其为结式或消去式。式(32)是单变量x1的多项式方程；消去了n-1个未知量；从而，可以获得单变量x1的可行解。若x1同时满足
[0296]
[0297] 则x1为正确解。将已解的x1代入式(34)，因式(32)成立且 任意，故得
[0298]
[0299] 即有
[0300]
[0301] 若有必要条件
[0302]
[0303] 成立，解式(35)，得被消去变量 的解；否则，需要结合式(16)得到全部解。考虑式(25)，因式(22)两边的x1阶次相等，故必有
[0304]
[0305] 若同时满足
[0306]
[0307] 则由式(35)能解得 中n-1个互不相同的组合变量；从而，得到所有独立变量的解。
[0308] 给定n个“n元N阶”多项式 Dixon矩阵计算步骤如下：
[0309] ①确定系统结构。方程数及独立变量数记为n；独立变量记为Xn；多项式复合变量记为 替换变量记为 替换变量数为n-1；大小为S·S的Dixon矩阵记为 其成员系数如式(24)所示，其中：S由式(32)确定；待消去变量为x1。
[0310] ②由式(8)得xα与 对应关系，表达式(11)中 至多有S项。
[0311] ③根据式(19)及Sarrus规则，计算Dixon(Fn(Yn-1|Xn-1))；根据 对应的N进位字运算结果，完成多项式合并。
[0312] ④Dixon矩阵成员如式(32)所示，由式(32)计算Dixon矩阵SΘS的(n+1)·S2个系数。
[0313] ⑤当满足式(37)及式(38)直接解判别准则时，由式(34)及式(35)得全部数值解。
[0314] 实施例3
[0315] 对多项式系统(39)进行Dixon消元。
[0316]
[0317] 步骤为：该式是4个“4元1阶”多项式系统，满足Dixon消元条件。由式(19)及式(22)，得
[0318]
[0319] 其中：
[0320]
[0321] 由式(34)及式(40)得5个解：
[0322]
[0323]
[0324]
[0325] 其中： 不是该方程组的解。将其它解分别代入式(35)。当 时，由式(35)得
[0326]
[0327] 解得：τ3＝1，τ4＝-2。将 τ3及τ4代入式(39)得τ2＝1。同样，可得其他三组解。显然，因变量不满足式(26)，式(40)所示的Dixon矩阵不对称。该例表明Dixon行列式为零对于多重线性多项式系统是充分的。
[0328] 基于轴不变量的定轴转动
[0329] 如图3所示，给定轴矢量 及与其固结的单位矢量 在转动前，对于单位矢量对系统零位轴 的投影矢量为 对系统径向轴 的矩矢量为 径向矢量为
[0330] 轴矢量 相对于杆件 及Ωl或自然坐标系 及F[l]是固定不变的，故称该转动为定轴转动。单位矢量 绕轴 转动 后，转动后的零位矢量 对系统零位轴 的投影矢量为 转动后的零位矢量 对系统径向轴 的矩矢量为轴向分量为 故得具有链指标的Rodrigues矢量方程
[0331]
[0332] 因单位矢量 是任意的且 得具有链指标的Rodrigues转动方程
[0333]
[0334] 若 由式(42)，得 若 即坐标系 与F[l]的方向一致，由式(42)可知：反对称部分 必有 因此，系统零位是自然坐标系 与F[l]
重合的充分必要条件，即初始时刻的自然坐标系方向一致是系统零位定义的前提条件。利用自然坐标系可以很方便地分析多轴系统运动学和动力学。
重合的充分必要条件，即初始时刻的自然坐标系方向一致是系统零位定义的前提条件。利用自然坐标系可以很方便地分析多轴系统运动学和动力学。
[0335]
[0336] 式(43)是关于 和 的多重线性方程，是轴不变量 的二阶多项式。给定自然零位矢量 作为 的零位参考，则 及 分别表示零位矢量及径向矢量。式
(43)即为 对称部分 表示零位轴张量，反对称部分
表示径向轴张量，分别与轴向外积张量 正交，从而确定三维自然轴空间；
式(43)仅含一个正弦及余弦运算、6个积运算及6个和运算，计算复杂度低；同时，通过轴不变量 及关节变量 实现了坐标系及极性的参数化。
(43)即为 对称部分 表示零位轴张量，反对称部分
表示径向轴张量，分别与轴向外积张量 正交，从而确定三维自然轴空间；
式(43)仅含一个正弦及余弦运算、6个积运算及6个和运算，计算复杂度低；同时，通过轴不变量 及关节变量 实现了坐标系及极性的参数化。
[0337] 对于轴链 有
[0338]
[0339] 由式(44)及式(43)得 则 是 及 的多重线性型，其中：l∈ilk。式(43)可表示为
[0340]
[0341] 称(45)为改进的Cayley变换。即有
[0342]
[0343] 由式(46)得规范的位置方程
[0344]
[0345] “居-吉布斯”四元数的确定：
[0346] 对于任意杆件l, 定义与欧拉四元数同构的“居-吉布斯”(Ju-Gibbs)规范四元数：
[0347]
[0348] 其中： 为Gibbs矢量。Gibbs共轭四元数为：
[0349]
[0350] 其中：
[0351]
[0352] 显然， 为 模的平方。因居-吉布斯四元数是四元数，故满足四元数乘法运算[0353]
[0354] 其中：
[0355]
[0356] 由式(52)得
[0357]
[0358] 习惯上，单关节及运动链的期望姿态以规范的Ju-Gibbs四元数(简称规范Ju-Gibbs四元数，即“标部”为1的四元数)表示；但是它们积运算通常是不规范的，即其标部不为1。由式(53)可知：只有给定轴l及 的规范Ju-Gibbs四元数，且两轴正交， 才为规范四元数。
[0359] 由式(53)得
[0360]
[0361] 由四维复数性质得
[0362]
[0363] 记 由式(52)得
[0364]
[0365] 故 为单位Ju-Gibbs四元数。
[0366] 由式(48)至式(50)及式(55)得
[0367]
[0368] 由式(50)、式(54)及式(57)得
[0369]
[0370] 类DCM及性质：
[0371] 对于轴链 规范的姿态方程为：
[0372]
[0373] 由式(59)得
[0374]
[0375] 式中， 为旋转变换矩阵； 表示用辅助变量yl的前l个依次替换原变量τl中的l个变量，记“|”为替换操作符；
[0376] 其中：
[0377]
[0378] 由式(61)可知：iQn及 是关于τk的n重2阶多项式。由式(60)可知：因 与 类似，故称之为类DCM(DCM，方向余弦矩阵)。由式(62)得
[0379]
[0380] 显然，类DCM可以通过Ju-Gibbs四元数表达。因此，式(59)姿态方程及式(47)位置方程是关于Ju-Gibbs四元数的表达式。
[0381] 分块方阵的逆：
[0382] 若给定可逆方阵K、B及C，其中B及C分别为l×l、c×c的方阵；A、D分别为l×c、c×l的矩阵，且
[0383]
[0384] 则有
[0385]
[0386] 基于轴不变量的Dixon行列式计算原理：
[0387] 下面基于轴不变量，提出径向不变量及运动链的Dixon行列式基本性质，为基于轴不变量的机器人逆运动学分析奠定基础。
[0388] 【1】轴不变量
[0389] 首先，轴不变量与坐标轴具有本质区别：坐标轴是具有零位及单位刻度的参考方向，可以描述沿轴向平动的线位置，但不能完整描述绕轴向的角位置，因为坐标轴自身不具有径向参考方向，即不存在表征转动的零位。在实际应用时，需要补充坐标轴的径向参考。坐标轴自身是1D的，3个正交的坐标轴构成3D的笛卡尔标架；轴不变量是3D空间单位参考轴(简称3D参考轴)，具有径向参考零位。“3D参考轴”及其径向参考零位可以确定对应的笛卡尔系。以自然坐标系为基础的轴不变量可以准确地反映运动轴及测量轴的“共轴性”、“极性”与“零位”三个基本属性。
[0390] 其次，轴不变量与欧拉轴具有本质的区别：方向余弦矩阵(DCM)是实矩阵，轴矢量是DCM的特征值1对应的特征矢量，是不变量；固定轴不变量是“3D参考轴”，不仅具有原点及轴向，也有径向参考零位；在自然坐标系下，轴不变量不依赖于相邻固结的自然坐标系，即在相邻固结的自然坐标系下具有不变的自然坐标；轴不变量具有幂零特性等优良的数学操作功能；在自然坐标系统中，通过轴不变量及关节坐标，可以唯一确定DCM及参考极性；没有必要为每一个杆件建立各自的体系，可以极大地简化建模的工作量。
[0391] 同时，以唯一需要定义的笛卡尔直角坐标系为参考，测量轴不变量，可以提高结构参数的测量精度。基于轴不变量的优良操作及属性，可以建立包含拓扑结构、坐标系、极性、结构参量及动力学参量的迭代式的运动学及动力学方程。
[0392] 由式(59)及式(47)可知：多轴系统的姿态及位置方程本质上是多元二阶多项式方程，其逆解本质上归结于多元二阶多项式的消元问题，包含Dixon矩阵及Dixon行列式计算的两个子问题。用式(47)的表达3R机械臂位置方程，是3个“3元2阶”多项式，应用Dixon消元方法计算逆解，有两个替换变量，在计算8×8的Dixon行列式时，最大可能的阶次为16。由式(4)可知：行列式计算是一个排列过程，面临着“组合爆炸”的难题。
[0393] 所有的不在确定的多项式时间内可解的问题称为NP问题。非确定性算法将问题分解为“猜测”与“验证”两个阶段：算法的“猜测”阶段具有非确定性，算法的“验证”阶段具有确定性，通过验证来确定猜测的解是否正确。假如可以在多项式时间内计算出来，就称为多项式非确定性问题。多元多项式的消元通常被认为是NP问题。通常应用 基进行多元多项式的消元，不得不求助于启发式的“猜测”与“验证”来解决问题。
[0394] 【2】径向不变量
[0395] 结构参数 及 是链节l的结构参量，在系统零位时，它们可以通过外部测量得到。如图4所示，零位矢量、径向矢量及轴向矢量是与转动角无关的不变量。其中，零位矢量是特定的径向矢量。
[0396] 任一个矢量可以分解为零位矢量及轴向矢量，故有
[0397]
[0398] 其中：
[0399]
[0400]
[0401] 考虑链节 其D-H参数有
[0402]
[0403] 显然， 是轴l及 的公垂线或公共径向矢量， 是轴l的轴向矢量。由式(65)可知：任一个结构参数矢量 可分解为与坐标系为无关的零位不变量 及轴向不变量 它们的径向矢量记为 结构参数矢量 及轴不变量 唯一确定径向坐标系，具有2个独立维度。
若两个轴向不变量 及 共线，则记为
若两个轴向不变量 及 共线，则记为
[0404]
[0405] 若两个零位不变量 及 与任两个径向不变量 及 共面，则记为
[0406]
[0407] 因此，称式(66)所示的轴向不变量及零位不变量是结构参数矢量对自然轴的分解。
[0408] 由式(69)及式(70)可知：同一个轴的三个径向矢量的行列式为零；同一个轴的任意两个轴向矢量的行列式为零。可以用轴不变量及其导出的不变量来简化Dixon行列式计算。
[0409] 由轴不变量导出的零位矢量、径向矢量及轴向矢量具有以下关系：
[0410]
[0411]
[0412]
[0413] 称式(71)为零位矢量的反转公式；称式(72)为零位矢量与径向矢量的互换公式；称式(73)为径向矢量不变性公式。由式(65)、式(71)至式(73)得
[0414]
[0415]
[0416] 由式(74)得
[0417]
[0418] 因 是 的对称部分的结构常数，故称式(74)为矢量 的对称分解式。因 是的反对称部分的结构常数，故称式(75)为矢量 的反对称分解式。称式(76)为归零等式。
[0419] 【3】运动链Dixon行列式性质
[0420] 定义
[0421]
[0422] 由式(52)得
[0423]
[0424] 其中：
[0425]
[0426] 由式(62)及式(66)得
[0427]
[0428] 由式(79)证得
[0429]
[0430] 式(80)可以将 及 可以转化为关于 的多重线性型。同时，对yl及tl具有对称(轮换)性。由式(67)、式(74)及式(75)得
[0431]
[0432] 式(81)由三个导出的独立结构参量 及一个运动变量tl构成。由式(81)得
[0433]
[0434]
[0435] 由式(80)及式(83)得
[0436]
[0437] 由式(80)及式(84)得
[0438]
[0439] 基于轴不变量的Cayley变换
[0440] 当给定角度 后，其正、余弦及其半角的正、余弦均是常数；为方便表达，记[0441]
[0442] 由式(86)得
[0443]
[0444] 定义
[0445]
[0446] 故有
[0447]
[0449] 基于轴不变量的3R机械臂位置逆解方法
[0450] 给定3R转动链 及期望姿态 轴不变量序列 求关节变量序列 这是3R姿态逆解问题。
[0451] 【1】根据机械臂n元3D矢量位姿方程，获得n个“n元2阶”多项式方程。
[0452] 由式(47)得3R运动学方程
[0453]
[0454] 由式(90)得
[0455]
[0456] 由式(91)得
[0457]
[0458] 若记
[0459]
[0460] 则由式(61)及式得(93)
[0461]
[0462] 由式(92)及式(93)得
[0463]
[0464] 下面，阐述3R机械臂运动学方程的Dixon行列式的结构模型及特点。
[0465] 由式(95)得3R运动学多项式方程
[0466]
[0467] 多项式系统F3(Y2|T2)，根据双线性型行列式通式
[0468]
[0469] 则有
[0470]
[0471] 其中：
[0472]
[0473]
[0474] 由式(18)、式(95)及式(96)得
[0475]
[0476] 由式(22)及式(101)可知式(99)成立。由式(80)及式(93)得
[0477]
[0478]
[0479] 由式(93)、式(102)及式(103)得
[0480]
[0481] 其中：应用式(85)计算 显然，式(104)中的y2阶次β2∈[0:3]及y3阶次β3∈[0:1]。考虑式(101)后三项：
中的y2阶次β2∈[0:3]及y3阶次β3∈[0:1]；
中的y2阶次β2∈[0:2]及y3阶次β3∈[0:1]；
中的y2的阶次β2∈[0:3]及y3的阶次β3∈[0:
1]。由上可知：式(101)中的y2阶次β2∈[0:3]及y3的阶次β3∈[0:1]。故有S＝8。
中的y2阶次β2∈[0:3]及y3阶次β3∈[0:1]；
中的y2阶次β2∈[0:2]及y3阶次β3∈[0:1]；
中的y2的阶次β2∈[0:3]及y3的阶次β3∈[0:
1]。由上可知：式(101)中的y2阶次β2∈[0:3]及y3的阶次β3∈[0:1]。故有S＝8。
[0482] 由式(93)、式(101)至式(104)可知： 中组合变量系数为独立的列向量，故选取的系数来构成方阵 剩余列向量一定与 的各列相关。故式(100)成立。
[0483] 【2】应用“基于轴不变量的Dixon行列式计算”方法，“分块矩阵的高维行列式计算”方法或者“对行列式进行行阶梯化计算”方法简化行列式计算。
[0484] 根据运动链Dixon行列式性质，由式(80)及式(93)得
[0485]
[0486]
[0487] 其中： 分别表示轴2至轴3、轴3至轴3S的零位矢量、径向矢量及轴向矢量。
[0488] 由式(105)得
[0489]
[0490] 由式(106)得
[0491]
[0492] 由式(107)得
[0493]
[0494] 由式(101)得
[0495]
[0496] 将式(108)至式(110)代入式(111)得
[0497]
[0498] 【3】应用n个“n元N阶”多项式的Dixon消元与求解原理完成位姿逆解计算，其中：根据Dixon矩阵的行列式为0，得到一元高阶多项式方程，应用基于友阵的一元高阶多项式方程求解一元高阶多项式方程的解。
[0499] 一元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的解。
[0500] 若多项式p(x)的友阵为
[0501]
[0502] 则由矩阵A的特征向量构成的矩阵为范德蒙德(Vandermonde)矩阵为
[0503]
[0504] 且有
[0505] p(λl)＝|A-λl·1n|＝0。 (115)
[0506] 由式(28)、式(99)及式(100)得
[0507]
[0508] 因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。一般情况下，式(116)是关于τ1的16阶单项式方程。
[0509] 该方法的过程表明：整体与局部、复杂与简单是对立统一的；式(4)将矢量多项式的行列式计算转化为三个矢量的行列式，这一步骤起到了决定性的作用；轴不变量及其导出的不变量都是结构参量，系统方程是关于结构参数的矢量与关节变量(标量)的矢量代数方程。
[0510] 以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 强制运动链伞齿轮滚磨机安装距控制与传动小拖板机构 | 2020-05-18 | 931 |
| 一种基于鱼群算法的机构运动链同构识别设计方法 | 2020-05-18 | 708 |
| 一种基于鱼群算法的机构运动链同构识别设计方法 | 2020-05-19 | 222 |
| 外啮合强制运动链滚磨机 | 2020-05-12 | 506 |
| 内啮合强制运动链滚磨机 | 2020-05-13 | 104 |
| 平面闭式复铰运动链的编号矩阵描述方法 | 2020-05-16 | 536 |
| 强制运动链伞齿轮滚磨机 | 2020-05-11 | 318 |
| 内啮合强制运动链立式滚磨机减速机构 | 2020-05-15 | 586 |
| 用于电力推动车辆的运动链 | 2020-05-12 | 263 |
| 强制运动链伞齿轮变交角可偏置滚磨机滚磨轮夹具 | 2020-05-20 | 570 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
分析报告
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
QQ群二维码
意见反馈
意见反馈