[发明专利]机械臂伺服系统的有限时间协同控制方法无效

专利信息
申请号: 201410398181.8 申请日: 2014-08-14
公开(公告)号: CN104216284A 公开(公告)日: 2014-12-17
发明(设计)人: 陈强;汤筱晴;翟双坡 申请(专利权)人: 浙江工业大学
主分类号: G05B13/04 分类号: G05B13/04
代理公司: 杭州天正专利事务所有限公司 33201 代理人: 王兵;黄美娟
地址: 310014 浙*** 国省代码: 浙江;33
权利要求书: 查看更多 说明书: 查看更多
摘要: 机械臂伺服系统的有限时间协同控制方法,包括:建立机械臂伺服系统的动态模型,初始化系统状态、采样时间以及相关控制参数;根据微分中值定理,将系统中的非线性输入死区线性近似为一个简单的时变系统,推导出带有未知死区的机械臂伺服系统模型;计算控制系统的跟踪误差、非奇异终端滑模面及其一阶导数;基于带有未知死区的机械臂伺服系统模型,选择神经网络逼近未知函数,并根据系统跟踪误差、非奇异终端滑模面,设计神经网络有限时间协同控制器,更新神经网络权值矩阵;该方法不仅能够避免非线性死区的附加逆补偿,减小滑模高频抖振问题,而且可以实现机械臂伺服系统的有限时间快速跟踪。
搜索关键词: 机械 伺服系统 有限 时间 协同 控制 方法
【主权项】:
机械臂伺服系统的有限时间协同控制方法,包含以下步骤:步骤1,建立如式(1)所示机械臂伺服系统的动态模型,初始化系统状态、采样时间以及相关控制参数;<mrow><mi>M</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mover><mi>x</mi><mrow><mo>.</mo><mo>.</mo></mrow></mover><mo>+</mo><msub><mi>V</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mover><mi>x</mi><mo>.</mo></mover><mo>)</mo></mrow><mover><mi>x</mi><mo>.</mo></mover><mo>+</mo><mi>G</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>+</mo><mi>F</mi><mrow><mo>(</mo><mover><mi>x</mi><mo>.</mo></mover><mo>)</mo></mrow><mo>=</mo><mi>&tau;</mi><mo>-</mo><msub><mi>&tau;</mi><mi>d</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>其中分别为位置,速度和加速度;M(x)∈Rn×n为对称正定矩阵;为科式力;G(x)∈Rn为重力;表示摩擦力;τd∈Rn×1为外部扰动;τ∈Rn×1为死区输出值,表示为:<mrow><mi>&tau;</mi><mo>=</mo><mi>D</mi><mrow><mo>(</mo><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>g</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow></mtd><mtd><mi>if</mi></mtd><mtd><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&GreaterEqual;</mo><msub><mi>b</mi><mi>r</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>if</mi></mtd><mtd><msub><mi>b</mi><mi>l</mi></msub><mo>&lt;</mo><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&lt;</mo><msub><mi>b</mi><mi>r</mi></msub></mtd></mtr><mtr><mtd><msub><mi>g</mi><mi>l</mi></msub><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow></mtd><mtd><mi>if</mi></mtd><mtd><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&le;</mo><msub><mi>b</mi><mi>l</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>其中u(t)∈R是实际控制信号;gl(u),gr(u)为未知非线性函数;bl和br为死区未知宽度参数,并且满足bl<0,br>0;步骤2,根据微分中值定理,将系统中的非线性输入死区线性近似为一个简单的时变系统,推导出带有未知死区的机械臂伺服系统模型;2.1根据微分中值定理,存在ξl∈(‑∞,bl)和ξr∈(br,+∞)使gl(u)=g′l(ξ′l)(u‑bl)(3)其中ξ′l∈(‑∞,bl],g′ll)为函数gl(u)在ξ′l处的倒数;gr(u)=g′r(ξ′r)(u‑br)(4)其中ξ′r∈[bl,+∞),g′rr)为函数gr(u)在ξ′r处的倒数;根据式(3)和式(4),可将式(2)改写为其中<mrow><mi>&rho;</mi><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mo>-</mo><msubsup><mi>g</mi><mi>r</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><msub><mi>&xi;</mi><mi>r</mi></msub><mo>)</mo></mrow><msub><mi>b</mi><mi>r</mi></msub></mtd><mtd><mi>if</mi></mtd><mtd><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&GreaterEqual;</mo><msub><mi>b</mi><mi>r</mi></msub></mtd></mtr><mtr><mtd><mo>-</mo><mo>[</mo><msubsup><mi>g</mi><mi>l</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><msub><mi>&xi;</mi><mi>l</mi></msub><mo>)</mo></mrow><mo>+</mo><msubsup><mi>g</mi><mi>r</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><msub><mi>&xi;</mi><mi>r</mi></msub><mo>)</mo></mrow><mo>]</mo><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd><mtd><mi>if</mi></mtd><mtd><msub><mi>b</mi><mi>l</mi></msub><mo>&lt;</mo><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&lt;</mo><msub><mi>b</mi><mi>r</mi></msub></mtd></mtr><mtr><mtd><mo>-</mo><msubsup><mi>g</mi><mi>l</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><msub><mi>&xi;</mi><mi>l</mi></msub><mo>)</mo></mrow></mtd><mtd><mi>if</mi></mtd><mtd><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&le;</mo><msub><mi>b</mi><mi>l</mi></msub></mtd></mtr></mtable></mfenced><mo>.</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>其中ξl∈(‑∞,bl],ξr∈[bl,+∞),并且2.2由式(1)和式(5)可得带有未知死区的机械臂伺服系统模型为:其中|ρ(u)|≤ρN,ρN是未知正常数,满足ρN=(gr1+gl1)max{br,bl};步骤3,计算控制系统的跟踪误差、非奇异终端滑模面及其一阶导数;3.1定义控制系统的跟踪误差为e(t)=xd‑x    (11)其中xd为二阶可导期望轨迹;3.2在设计协同控制器的过程中,定义协同流型为:<mrow><mi>s</mi><mo>=</mo><mover><mi>e</mi><mo>.</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>&lambda;e</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>其中λ为影响系统收敛的正常数;根据如式(13)所示的非奇异终端滑模面,设计有限时间协同控制器,使系统误差快速趋向于式(12)所示的协同流型;<mrow><mi>q</mi><msup><mover><mi>s</mi><mo>.</mo></mover><mrow><mi>p</mi><mo>/</mo><mi>r</mi></mrow></msup><mo>+</mo><mi>s</mi><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>其中q为正定常数,p>0,r>0,1<p/r<2;3.3对式(12)求导,结合式(6)和式(11)可以推导出其中非线性未知函数h为<mrow><mi>h</mi><mo>=</mo><mi>M</mi><mrow><mo>(</mo><msub><mover><mi>x</mi><mrow><mo>.</mo><mo>.</mo></mrow></mover><mi>d</mi></msub><mo>+</mo><mi>&lambda;</mi><mover><mi>e</mi><mo>.</mo></mover><mo>)</mo></mrow><mo>+</mo><msub><mi>V</mi><mi>m</mi></msub><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>.</mo></mover><mi>d</mi></msub><mo>+</mo><mi>&lambda;e</mi><mo>)</mo></mrow><mo>+</mo><mi>G</mi><mo>+</mo><mi>F</mi><mo>-</mo><mi>&rho;</mi><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow><mo>.</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mrow>步骤4,基于带有未知死区的机械臂伺服系统模型,选择神经网络逼近未知函数,并根据系统跟踪误差、非奇异终端滑模面,设计神经网络有限时间协同控制器,更新神经网络权值矩阵;4.1根据式(13)和式(14),控制器的表达式为<mrow><mi>u</mi><mo>=</mo><msubsup><mover><mi>W</mi><mo>^</mo></mover><mi>a</mi><mi>T</mi></msubsup><mi>&phi;</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>+</mo><mi>k</mi><msup><mi>s</mi><mrow><mi>r</mi><mo>/</mo><mi>p</mi></mrow></msup><mo>+</mo><mi>&delta;sgn</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>16</mn><mo>)</mo></mrow></mrow>其中为理想权重的估计值,为神经网络基函数,x为神经网络输入向量,k=q‑1是控制参数,δ是一个正的常数,满足εN为一个正常数,表示神经网络逼近误差的上限,τN为外部扰动的上限常数,为神经网络权值估计误差;4.2设计神经网络的权值更新律为:<mrow><msub><mover><mi>W</mi><mover><mo>^</mo><mo>.</mo></mover></mover><mi>a</mi></msub><mo>=</mo><msub><mi>K</mi><mi>C</mi></msub><mi>&phi;</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><msup><mi>s</mi><mi>T</mi></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>17</mn><mo>)</mo></mrow></mrow>其中KC为正定对角矩阵,φ(x)通常取为以下高斯函数:<mrow><mi>&phi;</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><msup><mrow><mo>|</mo><mi>x</mi><mo>-</mo><mi>c</mi><mo>|</mo></mrow><mn>2</mn></msup><mrow><mn>2</mn><msup><mi>b</mi><mn>2</mn></msup></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>18</mn><mo>)</mo></mrow></mrow>其中c=[c1,c2,...,cn]T是高斯函数的中心;b是高斯函数的宽度;由式(18)可以推出0<φ(x)≤1;步骤5,设计李雅普诺夫函数则可以证明闭环控制系统中的所有信号均是一致有界的;同时,系统误差e可以在有限时间内收敛至平衡点e=0。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于浙江工业大学,未经浙江工业大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/patent/201410398181.8/,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top