[发明专利]一种机械臂控制高速矩阵级联乘法计算实现方法在审
申请号: | 202211216732.5 | 申请日: | 2022-09-30 |
公开(公告)号: | CN115816439A | 公开(公告)日: | 2023-03-21 |
发明(设计)人: | 禹霁阳;汪路元;李欣;庞亚龙;刘伟伟;胡玉茜;李瑞军;范延芳;田苗苗;张思博;辛鹏飞;朱超 | 申请(专利权)人: | 北京空间飞行器总体设计部 |
主分类号: | B25J9/16 | 分类号: | B25J9/16;G06F17/16 |
代理公司: | 中国航天科技专利中心 11009 | 代理人: | 杨春颖 |
地址: | 100094 *** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 机械 控制 高速 矩阵 级联 乘法 计算 实现 方法 | ||
1.一种机械臂控制高速矩阵级联乘法计算实现方法,其特征在于,包括:
确定任意相邻的两个矩阵和其中,Mk和Nk表示第k个矩阵的行数和列数,Mk-1和Nk-1表示第k-1个矩阵的行数和列数,k=0,1,2,...,S-1,S≥1且S∈Z+;
对于矩阵和进行矩阵稀疏特征的等效降阶处理,并确定等效降阶后矩阵阶数;
根据确定的等效降阶后矩阵阶数,确定矩阵乘法最优次序;
按照确定的矩阵乘法最优次序,将级联矩阵乘法转换为并行点积计算,得到计算结果并输出。
2.根据权利要求1所述的机械臂控制高速矩阵级联乘法计算实现方法,其特征在于,矩阵表示用于连乘的S个矩阵Γ(k,k+1),Γ(k,k+1)表示S自由度空间机械臂相邻自由度关节坐标系之间的齐次变换矩阵。
3.根据权利要求2所述的机械臂控制高速矩阵级联乘法计算实现方法,其特征在于,Γ(k,k+1)表示如下:
其中,表示两个相邻关节坐标系之间的平移变化,ξ(k,k+1)表示两个相邻关节坐标系之间的旋转变化。
4.根据权利要求3所述的机械臂控制高速矩阵级联乘法计算实现方法,其特征在于,
其中,θk+1表示机械臂第k+1个关节绕其转轴的旋转角度。
5.根据权利要求1所述的机械臂控制高速矩阵级联乘法计算实现方法,其特征在于,S自由度空间机械臂整臂末端位姿齐次变换矩阵Γ(0,S)和矩阵连乘结果YM×N分别表示如下:
6.根据权利要求1所述的机械臂控制高速矩阵级联乘法计算实现方法,其特征在于,对于矩阵和进行矩阵稀疏特征的等效降阶处理,并确定等效降阶后矩阵阶数,包括:
对矩阵和进行四分块,得到矩阵对应的四分块矩阵和矩阵对应的四分块矩阵;
对各分块矩阵进行分块乘法计算,并对各被乘矩阵进行零行向量个数检测,对各乘矩阵进行零列向量个数检测;
根据各被乘矩阵零行向量个数检测结果和各乘矩阵零列向量个数检测结果,确定等效降阶后矩阵阶数。
7.根据权利要求6所述的机械臂控制高速矩阵级联乘法计算实现方法,其特征在于,通过如下公式对矩阵和进行四分块:
其中,表示矩阵中0至(Mk/2)-1行、0至(Nk/2)-1列的数据,表示矩阵中0至(Mk/2)-1行、(Nk/2)至Nk-1列的数据,表示矩阵中(Mk/2)至Mk-1行、0至(Nk/2)-1列的数据,表示矩阵中(Mk/2)至Mk-1行、(Nk/2)至Nk-1列的数据;表示矩阵中0至(Mk-1/2)-1行、0至(Nk-1/2)-1列的数据,表示矩阵中0至(Mk-1/2)-1行、(Nk-1/2)至Nk-1-1列的数据,表示矩阵中(Mk-1/2)至Mk-1-1行、0至(Nk-1/2)-1列的数据,表示矩阵中(Mk-1/2)至Mk-1-1行、(Nk-1/2)至Nk-1-1列的数据。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于北京空间飞行器总体设计部,未经北京空间飞行器总体设计部许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202211216732.5/1.html,转载请声明来源钻瓜专利网。