发明名称 适用于相对阶为1控制系统的平滑非奇异终端滑模控制方法
摘要 适用于相对阶为1控制系统的平滑非奇异终端滑模控制方法。本发明涉及一种非奇异终端滑模控制方法。本发明解决了现有非奇异终端滑模控制方法存在抖振而使控制器无法输出连续平滑的控制信号的问题,和无法应用于相对阶为1的控制系统的问题。若被控系统为相对阶为1的单输入单输出控制系统,则实时获取系统状态微分,设计非奇异终端滑动模态,引入虚拟控制量利用积分作用使得实际输出控制量平滑连续;若为匹配多输入多输出控制系统,则将控制系统的标量实现形式变换为矩阵矢量实现形式;若为非匹配不确定多输入多输出控制系统,则根据能控性指数r做两次非奇异状态变换分解成r个子系统,引入参考模型消除子系统输入通道的时变不确定性,进而设计控制律和辅助控制律。本发明应用于相对阶为1控制系统的平滑非奇异终端滑模控制。
申请公布号 CN104267605B 申请公布日期 2016.10.05
申请号 CN201410578260.7 申请日期 2014.10.25
申请人 哈尔滨工业大学 发明人 王艳敏;曹雨晴;夏红伟;申立群
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 杨立超
主权项 适用于相对阶为1控制系统的平滑非奇异终端滑模控制方法,其特征在于:所述控制方法通过以下步骤实现:步骤一、若被控系统为相对阶为1的单输入单输出控制系统A:<img file="FDA0001019005770000011.GIF" wi="507" he="62" />式中,<img file="FDA0001019005770000012.GIF" wi="30" he="46" />表示单输入单输出控制系统A的状态微分信号,x表示单输入单输出控制系统A的系统状态,u表示单输入单输出控制系统A的控制量,t表示时间,则执行步骤二的控制方法;若被控系统为相对阶为1的多输入多输出控制系统:<img file="FDA0001019005770000013.GIF" wi="693" he="63" />式中,<img file="FDA0001019005770000014.GIF" wi="46" he="55" />表示多输入多输出控制系统状态微分信号,x<sub>s</sub>表示多输入多输出控制系统状态,u<sub>s</sub>表示多输入多输出控制系统控制量,t表示时间,则执行步骤三的控制方法;步骤二、所述单输入单输出控制系统:<img file="FDA0001019005770000015.GIF" wi="254" he="55" />的控制方法具体为:步骤二一,由于单输入单输出控制系统A:<img file="FDA0001019005770000016.GIF" wi="254" he="55" />的相对阶为1,所述系统状态微分信号<img file="FDA0001019005770000017.GIF" wi="33" he="46" />在实际系统中为未知量,则利用高阶滑模鲁棒精确微分器实时获取<img file="FDA0001019005770000018.GIF" wi="59" he="46" /><maths num="0001"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mover><mi>y</mi><mo>&CenterDot;</mo></mover><mo>=</mo><msub><mi>v</mi><mn>0</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>v</mi><mn>0</mn></msub><mo>=</mo><msub><mi>v</mi><mn>1</mn></msub><mo>-</mo><msub><mi>&lambda;</mi><mn>0</mn></msub><mo>|</mo><mi>y</mi><mo>-</mo><mi>x</mi><msup><mo>|</mo><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><mi>s</mi><mi>i</mi><mi>g</mi><mi>n</mi><mrow><mo>(</mo><mi>y</mi><mo>-</mo><mi>x</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mover><mi>v</mi><mo>&CenterDot;</mo></mover><mn>1</mn></msub><mo>=</mo><mo>-</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><mi>s</mi><mi>i</mi><mi>g</mi><mi>n</mi><mrow><mo>(</mo><msub><mi>v</mi><mn>1</mn></msub><mo>-</mo><msub><mi>v</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0001019005770000019.GIF" wi="1950" he="255" /></maths>式中,λ<sub>0</sub>为设计参数一、λ<sub>1</sub>为设计参数二;y,v<sub>0</sub>和v<sub>1</sub>为高阶滑模鲁棒精确微分器状态的中间变量;步骤二二,设计非奇异终端滑动模态s(t),设<img file="FDA00010190057700000110.GIF" wi="322" he="62" />式中,非奇异终端滑动模态s(t)的设计参数c&gt;0;非奇异终端滑动模态s(t)的设计参数p、q为奇数,且满足p&gt;q&gt;0,1&lt;p/q&lt;2;步骤二三,引入虚拟控制量v  (4),使得所述系统状态x相对于虚拟控制量v的相对阶为2,增加了系统的相对阶;基于Lyapunov稳定定理设计虚拟控制量v,保证系统状态到达并维持在预先设计的滑模面s(t)=0上,并对参数摄动和外部扰动具有鲁棒性,且由于积分作用u=∫vdt  (5),使得实际输出控制量u平滑连续;步骤三、所述多输入多输出控制系统:<img file="FDA00010190057700000111.GIF" wi="318" he="55" />的具体控制方法为:步骤三一、若被控相对阶为1的多输入多输出控制系统:<img file="FDA0001019005770000021.GIF" wi="318" he="61" />内的不确定项或扰动项为匹配不确定性,即所述多输入多输出控制系统为匹配多输入多输出控制系统B:<maths num="0002"><math><![CDATA[<mrow><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>A</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>B</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0001019005770000022.GIF" wi="1950" he="63" /></maths>只需将步骤二一至步骤二三中单输入单输出控制系统A的标量实现形式变换为矩阵矢量实现形式即可,而平滑非奇异终端滑模控制方法因其具有的鲁棒性而对匹配不确定性扰动具有不变性;若被控系统为相对阶为1的非匹配不确定多输入多输出控制系统C:<maths num="0003"><math><![CDATA[<mrow><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>A</mi><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><mi>B</mi><mo>+</mo><mi>&Delta;</mi><mi>B</mi><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo></mrow><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>f</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0001019005770000023.GIF" wi="1950" he="71" /></maths>式中,x∈R<sup>n</sup>为系统状态,u∈R<sup>m</sup>为控制量,且1≤m≤n;A∈R<sup>n×n</sup>是已知常数矩阵;B∈R<sup>n×m</sup>是已知常数矩阵,矩阵B的维数为n<sub>1</sub>,n<sub>1</sub>≤m;(A,B)可控,设r为非匹配不确定多输入多输出控制系统C中常数矩阵(A,B)的能控性指数;ΔB(t)∈R<sup>n×m</sup>表示多输入通道内的匹配不确定性,即:ΔB(t)=Bd(t)  (8),式中,d(t)∈R<sup>m×m</sup>时变有界,其范围||d(t)||≤l<sub>d</sub>;f(t)表示非匹配不确定多输入多输出控制系统C的非匹配外部扰动,为R<sup>n</sup>上光滑有界函数,则执行步骤三二的控制方法;步骤三二、所述非匹配不确定多输入多输出控制系统C:<img file="FDA0001019005770000024.GIF" wi="758" he="70" />的具体控制方法为:步骤1,根据非匹配不确定多输入多输出控制系统C的能控性指数r,对非匹配不确定多输入多输出控制系统C做第一次非奇异状态变换:y=F<sub>1</sub>x   (9),式中,y表示第一次非奇异状态变换后的状态,x表示第一次非奇异状态变换前的系统状态,F<sub>1</sub>∈R<sup>n×n</sup>为变换矩阵;则非匹配不确定多输入多输出控制系统C变换为块控标准型控制系统D:<maths num="0004"><math><![CDATA[<mrow><mover><mi>y</mi><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msup><mi>A</mi><mo>&prime;</mo></msup><mi>y</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><msup><mi>B</mi><mo>&prime;</mo></msup><mo>+</mo><msup><mi>&Delta;B</mi><mo>&prime;</mo></msup><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo></mrow><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msup><mi>f</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0001019005770000025.GIF" wi="1950" he="71" /></maths>式中,A′=F<sub>1</sub>AF<sub>1</sub><sup>‑1</sup>表示系数矩阵A经第一次非奇异状态变换后的矩阵;<img file="FDA0001019005770000026.GIF" wi="442" he="87" />表示系数矩阵B经第一次非奇异状态变换后的矩阵,<img file="FDA0001019005770000027.GIF" wi="227" he="63" />表示B′的子矩阵,维数为n<sub>1</sub>;ΔB′(t)=B′d(t)表示多输入通道内的匹配不确定性;f′(t)表示块控标准型控制系统D的非匹配外部扰动,且存在f′(t)=F<sub>1</sub>f(t);步骤2,将步骤1得到的块控标准型控制系统D写为分块形式,相应地块控标准型控制系统D分解为内部子系统一和输入输出子系统一,内部子系统一:<img file="FDA0001019005770000031.GIF" wi="1622" he="126" />输入输出子系统一:<img file="FDA0001019005770000032.GIF" wi="1518" he="127" />控制量u仅出现在输入输出子系统一中;y表示第一次非奇异状态变换后的状态,<img file="FDA0001019005770000033.GIF" wi="317" he="55" /><img file="FDA0001019005770000034.GIF" wi="118" he="54" />为y的分量,维数为n<sub>i</sub>,i=1,…r,n<sub>1</sub>+…+n<sub>r</sub>=n;<img file="FDA0001019005770000035.GIF" wi="523" he="86" /><img file="FDA0001019005770000036.GIF" wi="182" he="63" />为f′的分量,<img file="FDA0001019005770000037.GIF" wi="181" he="63" />为f′的分量;A′<sub>i,j</sub>为内部子系统一的设计矩阵,i=2,…r;步骤3,将步骤2得到的存在状态耦合的内部子系统一做第二次非奇异状态变换:<maths num="0005"><math><![CDATA[<mrow><mi>z</mi><mo>=</mo><msubsup><mi>F</mi><mn>2</mn><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mi>y</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0001019005770000038.GIF" wi="1949" he="62" /></maths>式中,y是第一次非奇异状态变换后的状态,z为第二次非奇异状态变换后的状态,F<sub>2</sub>为变换阵,F<sub>2</sub>∈R<sup>n×n</sup>,<img file="FDA0001019005770000039.GIF" wi="893" he="375" />子矩阵<img file="FDA00010190057700000310.GIF" wi="1154" he="71" />子矩阵<img file="FDA00010190057700000311.GIF" wi="1739" he="126" /><img file="FDA00010190057700000312.GIF" wi="85" he="63" />为B<sub>i,i‑1</sub>的Moore‑Penrose逆,N<sub>i</sub>为变换矩阵F<sub>2</sub>中子矩阵的设计矩阵,相应地块控标准型控制系统D进一步转换为控制系统E:<img file="FDA00010190057700000313.GIF" wi="1078" he="71" />式中,A″=(F<sub>2</sub>)<sup>‑1</sup>A′F<sub>2</sub>表示系数矩阵A′经过第二次非奇异状态变换后的矩阵;<img file="FDA00010190057700000314.GIF" wi="502" he="79" />表示系数矩阵B′经过第二次非奇异状态变换后的矩阵;ΔB″(t)=B″d(t)表示第二次非奇异状态变换后的控制系统E多输入通道内的匹配不确定性;f″(t)=(F<sub>2</sub>)<sup>‑1</sup>f′(t)表示第二次非奇异状态变换后的控制系统E非匹配外部扰动;步骤4,将步骤3得到的控制系统E写成分块形式,相应地控制系统E分解为内部子系统二和输入输出子系统二:内部子系统二:<img file="FDA0001019005770000041.GIF" wi="1622" he="71" />输入输出子系统二:<img file="FDA0001019005770000042.GIF" wi="1518" he="122" />式中,<img file="FDA0001019005770000043.GIF" wi="297" he="77" /><img file="FDA0001019005770000044.GIF" wi="123" he="54" />表示z的分量;<img file="FDA0001019005770000045.GIF" wi="468" he="76" /><img file="FDA0001019005770000046.GIF" wi="221" he="67" />表示f″的分量i=1,…,r,<img file="FDA0001019005770000047.GIF" wi="221" he="63" />表示f″的分量;<img file="FDA0001019005770000048.GIF" wi="156" he="55" />为内部子系统二的设计矩阵,i=2,…r;步骤5,对步骤4得到的输入输出子系统二引入参考模型一来避免步骤4中无法应用传统方法求取B<sub>1,0</sub>d(t)逆矩阵的过程,所述参考模型一为:<maths num="0006"><math><![CDATA[<mrow><mover><mi>&xi;</mi><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>&alpha;</mi><mo>=</mo><mn>1</mn></mrow><mi>r</mi></munderover><msubsup><mi>A</mi><mrow><mn>1</mn><mo>,</mo><mi>&alpha;</mi></mrow><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><msub><mi>z</mi><mi>&alpha;</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&epsiv;B</mi><mrow><mn>1</mn><mo>,</mo><mn>0</mn></mrow></msub><msup><mi>u</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msup><mi>w</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0001019005770000049.GIF" wi="1950" he="127" /></maths>式中,<img file="FDA00010190057700000414.GIF" wi="122" he="54" />表示参考模型一的状态,ε∈R<sup>1</sup>表示参考模型一的设计参数,u′表示设计的参考模型一的控制律,<img file="FDA00010190057700000415.GIF" wi="154" he="51" />表示设计的参考模型一的辅助控制律。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号