发明名称 一种多机器人系统有限时间鲁棒协同跟踪控制方法
摘要 一种多机器人系统有限时间鲁棒协同跟踪控制方法,涉及多机器人系统的控制方法。为了解决现有的多机器人控制系统控制方法的鲁棒性较差的问题和多机器人系统的整体通讯负担过重的问题。本发明首先建立多机器人系统中跟随机器人的动力学模型<img file="DDA0000793306960000011.GIF" wi="984" he="91" /></maths>动力学模型可线性化为:<img file="DDA0000793306960000012.GIF" wi="1085" he="86" /></maths>定义变量q<sub>ri</sub>、z<sub>1i</sub>、z<sub>2i</sub>,结合虚拟控制器α<sub>1i</sub>得到<img file="DDA0000793306960000013.GIF" wi="484" he="84" /></maths>设计分布式控制律<img file="DDA0000793306960000014.GIF" wi="837" he="139" /></maths>和线性化参数自适应律<img file="DDA0000793306960000015.GIF" wi="286" he="101" />实现每个跟随机器人在有限时间内追随具有动态时变轨迹的领航机器人且跟踪误差有界,完成多机器人系统有限时间跟踪控制。本发明适用于多机器人系统的控制领域。
申请公布号 CN105068427A 申请公布日期 2015.11.18
申请号 CN201510547586.8 申请日期 2015.08.31
申请人 哈尔滨工业大学 发明人 姚俊羽;孙延超;李传江;王鹏宇;马广富;李东禹
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 杨立超
主权项 一种多机器人系统有限时间鲁棒协同跟踪控制方法,其特征在于包括以下步骤:步骤1、建立含有1个领航机器人和n个跟随机器人的多机器人系统,领航机器人标记为0,跟随机器人记为v<sub>L</sub>={1,2,…n};跟随机器人i∈v<sub>L</sub>的动力学模型可由如下Euler‑Lagrange方程描述:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>M</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>q</mi><mi>i</mi></msub><mo>)</mo></mrow><msub><mover><mi>q</mi><mo>&CenterDot;&CenterDot;</mo></mover><mi>i</mi></msub><mo>+</mo><msub><mi>C</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>q</mi><mi>i</mi></msub><mo>,</mo><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mi>i</mi></msub><mo>)</mo></mrow><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mi>i</mi></msub><mo>+</mo><msub><mi>g</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>q</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>&tau;</mi><mi>i</mi></msub><mo>,</mo><mi>i</mi><mo>&Element;</mo><msub><mi>v</mi><mi>L</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000793306930000011.GIF" wi="1454" he="87" /></maths>其中,q<sub>i</sub>∈R<sup>n</sup>为广义坐标;<img file="FDA0000793306930000012.GIF" wi="158" he="74" />为广义速度,<img file="FDA0000793306930000013.GIF" wi="158" he="74" />为广义加速度;M<sub>i</sub>(q<sub>i</sub>)∈R<sup>n×n</sup>为惯量矩阵,是对称正定的;<img file="FDA0000793306930000014.GIF" wi="343" he="77" />为Coriolis力/偏心力;g<sub>i</sub>(q<sub>i</sub>)∈R<sup>n</sup>为广义有势力;τ<sub>i</sub>∈R<sup>n</sup>为作用于跟随机器人i的广义控制力矢量;n是多机器人系统中机器人的维数;对于任意向量x∈R<sup>n</sup>,y∈R<sup>n</sup>,式(1)所示的动力学模型线性化为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>M</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>q</mi><mi>i</mi></msub><mo>)</mo></mrow><mi>x</mi><mo>+</mo><msub><mi>C</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>q</mi><mi>i</mi></msub><mo>,</mo><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mi>i</mi></msub><mo>)</mo></mrow><mi>y</mi><mo>+</mo><msub><mi>g</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>q</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><mo>-</mo><msub><mi>Y</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>q</mi><mi>i</mi></msub><mo>,</mo><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mi>i</mi></msub><mo>,</mo><mi>x</mi><mo>,</mo><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mo>)</mo></mrow><msub><mi>&Theta;</mi><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000793306930000015.GIF" wi="1436" he="85" /></maths>其中,Y<sub>i</sub>为回归矩阵,Θ<sub>i</sub>为包含跟随机器人i物理参数的常值向量;步骤2、根据多机器人系统的通讯拓扑,计算多机器人系统的有向图图论中的加权邻接矩阵A和Laplacian矩阵;步骤3、定义以下变量:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>q</mi><mrow><mi>r</mi><mi>i</mi></mrow></msub><mo>=</mo><mo>-</mo><mi>a</mi><munderover><mo>&Sigma;</mo><mrow><msup><mi>j</mi><mo>*</mo></msup><mo>=</mo><mn>0</mn></mrow><mi>n</mi></munderover><msub><mi>a</mi><mrow><msup><mi>ij</mi><mo>*</mo></msup></mrow></msub><mrow><mo>(</mo><msub><mi>q</mi><mi>i</mi></msub><mo>-</mo><msub><mi>q</mi><msup><mi>j</mi><mo>*</mo></msup></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000793306930000016.GIF" wi="1346" he="160" /></maths>z<sub>1i</sub>=q<sub>i</sub>‑q<sub>ri</sub>     (4)<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>z</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mo>=</mo><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mi>i</mi></msub><mo>-</mo><msub><mi>&alpha;</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000793306930000017.GIF" wi="1217" he="70" /></maths>其中,q<sub>ri</sub>为辅助变量,z<sub>1i</sub>为跟踪误差,z<sub>2i</sub>为构造的误差,a为一正常数,<img file="FDA0000793306930000019.GIF" wi="70" he="71" />为加权邻接矩阵A的元素,元素<img file="FDA00007933069300000110.GIF" wi="70" he="71" />表示跟随机器人i和机器人j<sup>*</sup>之间信息传递关系;i和j表示参数针对的是跟随机器人,i和j分别属于v<sub>L</sub>={1,2,…n};i<sup>*</sup>和j<sup>*</sup>表示参数针对的是所有机器人,i<sup>*</sup>和j<sup>*</sup>分别属于v={0,1,2,…,n};α<sub>1i</sub>为虚拟控制器,其形式如式(6)所示:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>&alpha;</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mo>=</mo><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mrow><mi>r</mi><mi>i</mi></mrow></msub><mo>-</mo><msub><mi>K</mi><mn>1</mn></msub><mi>S</mi><mi>i</mi><mi>g</mi><msup><mrow><mo>(</mo><msub><mi>z</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mo>)</mo></mrow><mi>&alpha;</mi></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000793306930000018.GIF" wi="1282" he="85" /></maths>其中,K<sub>1</sub>为参数矩阵,K<sub>1</sub>=diag(K<sub>11</sub>,K<sub>12</sub>,…,K<sub>1n</sub>);K<sub>11</sub>,K<sub>12</sub>,…,K<sub>1n</sub>为K<sub>1</sub>中的元素;α∈(0,1);<img file="FDA0000793306930000021.GIF" wi="66" he="75" />为q<sub>ri</sub>的一阶导数;对变量z<sub>1i</sub>求导后,再引入如公式(6)所示的虚拟控制器,可得<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mrow><mn>1</mn><mi>i</mi></mrow></msub><mo>=</mo><msub><mi>z</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mo>-</mo><msub><mi>K</mi><mn>1</mn></msub><mi>S</mi><mi>i</mi><mi>g</mi><msup><mrow><mo>(</mo><msub><mi>z</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mo>)</mo></mrow><mi>&alpha;</mi></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000793306930000022.GIF" wi="1275" he="85" /></maths>其中,函数Sig(·)<sup>α</sup>定义为Sig(ξ)<sup>α</sup>=[|ξ<sub>1</sub>|<sup>α</sup>sgn(ξ<sub>1</sub>),…,|ξ<sub>n</sub>|<sup>α</sup>sgn(ξ<sub>n</sub>)]<sup>T</sup>其中,ξ=[ξ<sub>1</sub>,ξ<sub>2</sub>,…,ξ<sub>n</sub>]<sup>T</sup>∈R<sup>n</sup>,sgn(·)为符号函数;步骤4、设计分布式跟踪控制律和线性化参数自适应律,完成多机器人系统有限时间跟踪控制。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号