发明名称 一种工业机械手抗干扰鲁棒自适应PID控制方法
摘要 本发明涉及一种工业机械手抗干扰鲁棒自适应PID控制方法,其特征在于:鲁棒自适应PID控制器包括PID控制项、自适应控制项、鲁棒控制项,PID控制项与鲁棒控制项的输入均为机械手的位置与速度误差,自适应控制项的输入为机械手不确定动力学参数的估计值,PID控制项、自适应控制项、鲁棒控制项的输出经过累加器叠加,实现鲁棒自适应PID控制。
申请公布号 CN103901776B 申请公布日期 2016.08.17
申请号 CN201410136080.3 申请日期 2014.04.04
申请人 哈尔滨工程大学 发明人 徐健;乔磊;汪慢;周佳加;张耕实
分类号 G05B13/04(2006.01)I;B25J9/18(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 代理人
主权项 一种工业机械手抗干扰鲁棒自适应PID控制方法,其特征在于:鲁棒自适应PID控制器包括PID控制项、自适应控制项、鲁棒控制项,PID控制项与鲁棒控制项的输入均为机械手的位置与速度误差,自适应控制项的输入为机械手不确定动力学参数的估计值,PID控制项、自适应控制项、鲁棒控制项的输出经过累加器叠加,实现鲁棒自适应PID控制;扰动信号上确界已知时的鲁棒自适应PID控制器具体通过如下方法设计,1)被控对象的数学描述:对于一个具有N个转动关节的机械手系统,其动力学可由以下二阶非线性微分方程描述:<maths num="0001"><math><![CDATA[<mrow><mi>M</mi><mrow><mo>(</mo><mi>q</mi><mo>)</mo></mrow><mover><mi>q</mi><mo>&CenterDot;&CenterDot;</mo></mover><mo>+</mo><mi>C</mi><mrow><mo>(</mo><mi>q</mi><mo>,</mo><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mo>)</mo></mrow><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mo>+</mo><mi>G</mi><mrow><mo>(</mo><mi>q</mi><mo>)</mo></mrow><mo>+</mo><mi>u</mi><mo>=</mo><mi>&tau;</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000949420870000011.GIF" wi="1346" he="55" /></maths>式(1)中,q∈R<sup>n</sup>——机械手关节角位移向量;M(q)∈R<sup>n×n</sup>——机械手的惯性矩阵;<img file="FDA0000949420870000012.GIF" wi="277" he="63" />——机械手的离心力和哥式力矩阵;G(q)∈R<sup>n</sup>——机械手的重力项;τ∈R<sup>n</sup>——机械手的外部控制力矩;u∈R<sup>n</sup>——机械手所受的各种外界扰动项;2)被控对象的动力学特性:①M(q)是正定对称矩阵;②<img file="FDA0000949420870000013.GIF" wi="331" he="71" />为斜对称矩阵;③存在一个依赖于机械手参数的参数向量,使得M(q),<img file="FDA0000949420870000014.GIF" wi="189" he="54" />G(q)满足如下线性关系:<maths num="0002"><math><![CDATA[<mrow><mi>M</mi><mrow><mo>(</mo><mi>q</mi><mo>)</mo></mrow><msub><mover><mi>q</mi><mo>&CenterDot;&CenterDot;</mo></mover><mi>k</mi></msub><mo>+</mo><mi>C</mi><mrow><mo>(</mo><mi>q</mi><mo>,</mo><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mo>)</mo></mrow><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mi>k</mi></msub><mo>+</mo><mi>G</mi><mrow><mo>(</mo><mi>q</mi><mo>)</mo></mrow><mo>=</mo><mi>&Psi;</mi><mrow><mo>(</mo><mi>q</mi><mo>,</mo><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mo>,</mo><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mi>k</mi></msub><mo>,</mo><msub><mover><mi>q</mi><mo>&CenterDot;&CenterDot;</mo></mover><mi>k</mi></msub><mo>)</mo></mrow><mi>P</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000949420870000015.GIF" wi="1349" he="62" /></maths>式(2)中,<img file="FDA0000949420870000021.GIF" wi="419" he="63" />为已知关节变量函数回归矩阵,即为机械手广义坐标及其各阶导数的已知函数矩阵;P∈R<sup>m</sup>为描述机械手质量特性的未知定常参数向量;3)扰动信号上确界已知时的机械手鲁棒自适应PID控制器设计:针对式(1)描述的机械手系统,在扰动信号上确界已知时,鲁棒自适应PID控制器设计如下:<maths num="0003"><math><![CDATA[<mrow><mi>&tau;</mi><mo>=</mo><mo>-</mo><msub><mi>K</mi><mi>P</mi></msub><mi>e</mi><mo>-</mo><msub><mi>K</mi><mi>D</mi></msub><mover><mi>e</mi><mo>&CenterDot;</mo></mover><mo>-</mo><msub><mi>K</mi><mi>I</mi></msub><mrow><mo>(</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mi>t</mi></msubsup><mrow><mi>e</mi><mi>d</mi><mi>t</mi></mrow><mo>)</mo></mrow><mo>+</mo><mi>&Psi;</mi><mrow><mo>(</mo><mi>q</mi><mo>,</mo><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mo>,</mo><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mi>k</mi></msub><mo>,</mo><msub><mover><mi>q</mi><mo>&CenterDot;&CenterDot;</mo></mover><mi>k</mi></msub><mo>)</mo></mrow><mover><mi>P</mi><mo>^</mo></mover><mo>+</mo><mi>V</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000949420870000022.GIF" wi="1465" he="94" /></maths>V=[v<sub>1</sub>, … v<sub>n</sub>]<sup>T</sup>    (4)<maths num="0004"><math><![CDATA[<mrow><msub><mi>v</mi><mi>i</mi></msub><mo>=</mo><mo>-</mo><mrow><mo>(</mo><msub><mi>b</mi><mn>1</mn></msub><mo>+</mo><msub><mi>b</mi><mn>2</mn></msub><mo>|</mo><mo>|</mo><mi>e</mi><mo>|</mo><mo>|</mo><mo>+</mo><msub><mi>b</mi><mn>3</mn></msub><mo>|</mo><mo>|</mo><mover><mi>e</mi><mo>&CenterDot;</mo></mover><mo>|</mo><mo>|</mo><mo>+</mo><msub><mi>b</mi><mn>4</mn></msub><mo>|</mo><mo>|</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mi>t</mi></msubsup><mi>e</mi><mi>d</mi><mi>t</mi><mo>|</mo><mo>|</mo><mo>)</mo></mrow><mi>s</mi><mi>g</mi><mi>n</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000949420870000023.GIF" wi="1370" he="116" /></maths><img file="FDA0000949420870000024.GIF" wi="44" he="62" />为P的估计值,取<img file="FDA0000949420870000025.GIF" wi="43" he="55" />的估计率为:<maths num="0005"><math><![CDATA[<mrow><mover><mover><mi>P</mi><mo>^</mo></mover><mo>&CenterDot;</mo></mover><mo>=</mo><mo>-</mo><msup><mi>&Phi;&Psi;</mi><mi>T</mi></msup><mrow><mo>(</mo><mi>q</mi><mo>,</mo><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mo>,</mo><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mi>k</mi></msub><mo>,</mo><msub><mover><mi>q</mi><mo>&CenterDot;&CenterDot;</mo></mover><mi>k</mi></msub><mo>)</mo></mrow><mi>x</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000949420870000026.GIF" wi="1162" he="88" /></maths>取控制增益矩阵为:K<sub>P</sub>=diag[k<sub>P1</sub>, …, k<sub>Pn</sub>]   (7)K<sub>D</sub>=diag[k<sub>D1</sub>, …, k<sub>Dn</sub>]    (8)K<sub>I</sub>=diag[k<sub>I1</sub>, …, k<sub>In</sub>]    (9)其中,e=q‑q<sub>d</sub>为机械手关节角位置误差,<img file="FDA0000949420870000027.GIF" wi="189" he="55" />为机械手关节角速度误差,q<sub>d</sub>为期望的关节角度;<img file="FDA0000949420870000028.GIF" wi="822" he="94" />γ为常数,且γ>0;b<sub>1</sub>,b<sub>2</sub>,b<sub>3</sub>和b<sub>4</sub>均为大于零的常数;矩阵Φ∈R<sup>m×m</sup>为正定对称阵,k<sub>Pi</sub>,k<sub>Di</sub>,k<sub>Ii</sub>均为正常数,且k<sub>Di</sub>=k<sub>Ii</sub>,i=1,2,…,n。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室