[发明专利]基于全向移动机器人运动学建模的误差模型预测控制方法有效
申请号: | 201910142417.4 | 申请日: | 2019-02-26 |
公开(公告)号: | CN109885052B | 公开(公告)日: | 2022-03-25 |
发明(设计)人: | 黄林青;魏武;王栋梁;孙金权;周方华;罗永恒 | 申请(专利权)人: | 华南理工大学 |
主分类号: | G05D1/02 | 分类号: | G05D1/02;G05D1/00 |
代理公司: | 广州粤高专利商标代理有限公司 44102 | 代理人: | 何淑珍;黄海波 |
地址: | 510640 广*** | 国省代码: | 广东;44 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 基于 全向 移动 机器人 运动学 建模 误差 模型 预测 控制 方法 | ||
1.基于全向移动机器人运动学建模的误差模型预测控制方法,其特征在于,所述方法包括以下步骤:
S11、建立FM-OMR四个麦克纳姆轮之间的速度约束运动学模型,具体包括:
S111、首先对FM-OMR的结构进行描述;
S112、其次对FM-OMR的速度约束进行描述:
将车轮的转速表示为车轮的平移速度表示为而Vir,i=1,2,3,4,是自由移动辊的切向速度,将机器人框架中轮子的合成速度表示为:
V1x=V1m+V1rcosα1;V1y=V1rsinα1
V2x=V2m+V2rcosα2;V2y=V2rsinα2
V3x=V3m+V3rcosα3;V3y=V3rsinα3
V4x=V4m+V4rcosα4;V4y=V4rsinα4
其中α1=α3=45,α2=α4=-45,V1m、V2m、V3m、V4m分别表示四个麦克纳姆轮的平移速度,V1y、V2y、V3y、V4y表示分别表示四个麦克纳姆轮沿y轴的合计速度,此外,由于轮速与机器人本体的速度是刚性连接的,四个轮子速度V1x、V2x、V3x、V4x与机器人本体的速度关系建立为:
V1x=vx-l1w;V1y=vy+l2w
V2x=vx+l1w;V2y=vy+l2w
V3x=vx-l1w;V3y=vy-l2w
V4x=vx+l1w;V4y=vy-l2w
因此,获得FM-OMR机体的速度以及四个轮子的角速度的关系如下:
其中vx、vy是机体沿着x、y方向的速度,ω是机体沿着z方向的角速度,其中R是如下雅可比矩阵:
其中l1、l2表示车轮到质心的纵向和横向距离;Rw表示麦克纳姆轮的半径;
其中矩阵R不是方阵,因此没有逆矩阵,然而,广义逆矩阵表示为:
因此得到:
由上面获得:
从上式得到以下等式:
ω1+ω2=ω3+ω4
调整四个麦克纳姆轮各自的速度,实现在给定的平面内期望的全向运动,每个轮子的约束产生作用于整个机体,四个麦克纳姆轮之间的速度约束关系由ω1+ω2=ω3+ω4表示,车身速度与四个轮子角速度之间的关系已经得出,机器人本体的速度表示如下:
其中,
S113、最后构建FM-OMR的运动学模型:
从上述式子,得出非完整FM-OMR的运动学模型如下:
S12、建立FM-OMR的跟踪误差运动学模型,具体过程为:
S121、确定领导者Rr以及跟随者R的物理参数,其中对于领导者小车,有:
其中,l3、l4表示车轮到质心的纵向和横向距离;
S122、根据领导者以及跟随者在世界坐标系下的位置和角度计算Rr在R坐标系下的相对位置和角度(xe,ye,θe):
其中,(x,y,θ)以及(xr,yr,θr)是分别都是对应于R和Rr的惯性协同坐标系;
S123、根据所得的(xe,ye,θe)计算其导数,建立其导数与输入[u1,u2,u3]之间的状态方程,所述导数表示为:
其中,
a=w1r+w2r+w3r+w4r
b=w1r-w2r-w3r+w4r
c=w1+w2+w3+w4
d=w1-w2-w3+w4
其中,a、b是Rr四个车轮角速度的线性组合,c、d是R四个车轮角速度的线性组合;
则综上四个麦克纳姆轮之间具有速度约束的轨迹跟踪误差运动学模型为:
其中,
u3=wr-w;
S124、通过控制输入[u1,u2,u3],误差状态量[xe,ye,θe]收敛到原点;
S13、针对FM-OMR的轨迹跟踪问题,针对跟踪误差运动学模型,设计结合速度约束方程的误差模型预测控制器;设计结合速度约束方程的误差模型预测控制器时,连续时变系统用离散时间模型的形式重新表示,每个特征点的离散时间状态空间模型表示为:
ze(k+1)=L(k)ze(k)+M(k)u(k)
其中,表示输入向量,表示状态量,参数n和m表示状态和输入变量的数量;根据采样定理,给定连续时间系统用采样周期T离散化,因此,上式写成如下形式:
所述状态空间模型中,为了导出给定预测层Np的最优控制序列,损失函数表示为:
其中,Δu(k)表示输入向量的增量,有Δu(k)=u(k)-u(k-1),ze(k+j|k)和Δu(k+j|k)表示基于当前时间k,在时间k+j,上的预测变量,此外,Nc代表控制层,Np代表预测层,最后,Q和P表示适当的加权矩阵,T为采样周期;
从以上式子得出:
定义如下形式的预测向量:
其中,Np=3,Nc=2,且对于所有k≥0时刻,约束表示如下:
其中状态量ze的最大值和最小值输入向量的最大值和最小值输入向量的增量Δu(k)的最小值和最大值均为某一给定的常数;
预测状态表示为如下形式:
最终,基于约束条件下代价函数的EMPC算法通过解以下约束来对问题进行优化:
这取决于:
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于华南理工大学,未经华南理工大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201910142417.4/1.html,转载请声明来源钻瓜专利网。