发明名称 微陀螺仪的自适应反演非奇异终端滑模控制方法
摘要 本发明提出了一种微陀螺仪的自适应反演非奇异终端滑模控制方法,将反演法与终端滑模相结合,根据Lyapunov稳定性理论设计反演终端滑模控制律,使系统状态能够在很短的有限时间内收敛到平衡点,同时利用自适应控制辨识出微陀螺仪的角速度和其它系统参数,进一步考虑到现有终端滑模控制存在的奇异性问题,引入非奇异终端滑模控制,实现了系统的全局非奇异控制,控制器响应速度完全可以和传统终端滑模相媲美,具有重要的理论和实际应用价值。本发明在保证收敛速度和跟踪性能的同时,对外界干扰具有较强的鲁棒性和自适应能力。
申请公布号 CN103728882A 申请公布日期 2014.04.16
申请号 CN201410006856.X 申请日期 2014.01.07
申请人 河海大学常州校区 发明人 严维锋;费峻涛
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 南京纵横知识产权代理有限公司 32224 代理人 董建林
主权项 1.微陀螺仪的自适应反演非奇异终端滑模控制方法,其特征在于,包括以下步骤:1)构建微陀螺仪系统的数学模型为:<maths num="0001"><![CDATA[<math><mrow><mover><mi>q</mi><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover><mo>+</mo><mrow><mo>(</mo><mi>D</mi><mo>+</mo><mn>2</mn><mi>&Omega;</mi><mo>)</mo></mrow><mrow><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mo>+</mo><mi>Kq</mi><mo>=</mo><mi>u</mi><mo>+</mo><mi>f</mi></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,q为微陀螺仪的质量块在驱动轴和感测轴两轴的位置向量,为微陀螺仪系统的输出;u为微陀螺仪的控制输入;D为阻尼矩阵;K包含了两轴的固有频率和耦合的刚度系数;Ω为角速率矩阵;f为系统的参数不确定性和外部干扰;2)构建自适应反演非奇异终端滑模控制器,得到自适应反演非奇异终端滑模控制律:2-1)定义变量X<sub>1</sub>和X<sub>2</sub>,令:<img file="FDA0000454131740000012.GIF" wi="331" he="76" />基于反演设计理论,将微陀螺仪的动力学方程式(3)变化为如下形式:<maths num="0002"><![CDATA[<math><mfenced open='{' close='' separators=''><mtable><mtr><mtd><msub><mover><mi>X</mi><mo>&CenterDot;</mo></mover><mn>1</mn></msub><mo>=</mo><msub><mi>X</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mover><mi>X</mi><mo>&CenterDot;</mo></mover><mn>2</mn></msub><mo>=</mo><mo>-</mo><mrow><mo>(</mo><mi>D</mi><mo>+</mo><mn>2</mn><mi>&Omega;</mi><mo>)</mo></mrow><msub><mi>X</mi><mn>2</mn></msub><mo>-</mo><mi>K</mi><msub><mi>X</mi><mn>1</mn></msub><mo>+</mo><mi>u</mi><mo>+</mo><mi>f</mi></mtd></mtr></mtable><mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></mfenced></math>]]></maths>2-2)定义误差变量e<sub>2</sub>为:e<sub>2</sub>=X<sub>2</sub>-α<sub>1</sub>    (11)其中q<sub>r</sub>为质量块沿两轴的理想位置输出向量,α<sub>1</sub>为虚拟控制量,且有:<maths num="0003"><![CDATA[<math><mrow><msub><mi>&alpha;</mi><mn>1</mn></msub><mo>=</mo><mo>-</mo><msub><mi>c</mi><mn>1</mn></msub><msub><mi>e</mi><mn>1</mn></msub><mo>+</mo><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mi>r</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中c<sub>1</sub>为误差系数,是非零正常数,e<sub>1</sub>=X<sub>1</sub>-q<sub>r</sub>;2-3)定义非奇异终端滑模面s<sub>c</sub>为:<maths num="0004"><![CDATA[<math><mrow><msub><mi>s</mi><mi>c</mi></msub><mo>=</mo><msub><mi>e</mi><mn>1</mn></msub><mo>+</mo><mfrac><mn>1</mn><mi>&beta;</mi></mfrac><msubsup><mi>e</mi><mn>2</mn><mrow><msub><mi>p</mi><mn>1</mn></msub><mo>/</mo><msub><mi>p</mi><mn>2</mn></msub></mrow></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>16</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,β,p<sub>1</sub>,p<sub>2</sub>为滑模面参数,满足:β>0,p<sub>1</sub>,p<sub>2</sub>为奇数,且1<p<sub>1</sub>/p<sub>2</sub><2;2-4)对于所述微陀螺仪系统,采用式(16)所述的滑模面,设计反演非奇异终端滑模控制律φ由四个控制律u<sub>0</sub>,u<sub>1</sub>,u<sub>2</sub>,u<sub>3</sub>组成:φ=u<sub>0</sub>+u<sub>1</sub>+u<sub>2</sub>+u<sub>3</sub>    (19)其中,<maths num="0005"><![CDATA[<math><mrow><msub><mi>u</mi><mn>0</mn></msub><mo>=</mo><mrow><mo>(</mo><mi>D</mi><mo>+</mo><mn>2</mn><mi>&Omega;</mi><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>e</mi><mn>2</mn></msub><mo>+</mo><msub><mi>&alpha;</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><mi>K</mi><mrow><mo>(</mo><msub><mi>e</mi><mn>1</mn></msub><mo>+</mo><msub><mi>q</mi><mi>r</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mover><mi>&alpha;</mi><mo>&CenterDot;</mo></mover><mn>1</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>20</mn><mo>)</mo></mrow></mrow></math>]]></maths>D,K,Ω分别为微陀螺仪的三个参数矩阵;<maths num="0006"><![CDATA[<math><mrow><msub><mi>u</mi><mn>1</mn></msub><mo>=</mo><mo>-</mo><mi>&beta;</mi><mfrac><msub><mi>p</mi><mn>2</mn></msub><msub><mi>p</mi><mn>1</mn></msub></mfrac><mi>diag</mi><mrow><mo>(</mo><msup><msub><mi>e</mi><mn>2</mn></msub><mrow><mn>1</mn><mo>-</mo><msub><mi>p</mi><mn>1</mn></msub><mo>/</mo><msub><mi>p</mi><mn>2</mn></msub></mrow></msup><mo>)</mo></mrow><mover><msub><mi>e</mi><mn>1</mn></msub><mo>&CenterDot;</mo></mover><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>21</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><msub><mi>u</mi><mn>2</mn></msub><mo>=</mo><mo>-</mo><mi>&beta;</mi><mfrac><msub><mi>p</mi><mn>2</mn></msub><msub><mi>p</mi><mn>1</mn></msub></mfrac><mi>diag</mi><mrow><mo>(</mo><msup><msub><mi>e</mi><mn>2</mn></msub><mrow><mn>1</mn><mo>-</mo><msub><mi>p</mi><mn>1</mn></msub><mo>/</mo><msub><mi>p</mi><mn>2</mn></msub></mrow></msup><mo>)</mo></mrow><mfrac><msub><mi>s</mi><mi>c</mi></msub><msup><mrow><msub><mrow><mo>|</mo><mo>|</mo><mi>s</mi></mrow><mi>c</mi></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup></mfrac><mrow><mo>(</mo><msup><msub><mi>e</mi><mn>1</mn></msub><mi>T</mi></msup><msub><mi>e</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>22</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><msub><mi>u</mi><mn>3</mn></msub><mo>=</mo><mo>-</mo><mi>&rho;</mi><mfrac><msup><mrow><mo>[</mo><msup><msub><mi>s</mi><mi>c</mi></msub><mi>T</mi></msup><mfrac><mn>1</mn><mi>&beta;</mi></mfrac><mi>diag</mi><mrow><mo>(</mo><msup><msub><mi>e</mi><mn>2</mn></msub><mrow><msub><mi>p</mi><mn>1</mn></msub><mo>/</mo><msub><mi>p</mi><mn>2</mn></msub><mo>-</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><mo>]</mo></mrow><mi>T</mi></msup><msup><mrow><mo>|</mo><mo>|</mo><msup><msub><mi>s</mi><mi>c</mi></msub><mi>T</mi></msup><mfrac><mn>1</mn><mi>&beta;</mi></mfrac><mi>diag</mi><mrow><mo>(</mo><msup><msub><mi>e</mi><mn>2</mn></msub><mrow><msub><mi>p</mi><mn>1</mn></msub><mo>/</mo><msub><mi>p</mi><mn>2</mn></msub><mo>-</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup></mfrac><mo>|</mo><mo>|</mo><msub><mi>s</mi><mi>c</mi></msub><mo>|</mo><mo>|</mo><mo>|</mo><mo>|</mo><mfrac><mn>1</mn><mi>&beta;</mi></mfrac><mi>diag</mi><mrow><mo>(</mo><msup><msub><mi>e</mi><mn>2</mn></msub><mrow><msub><mi>p</mi><mn>1</mn></msub><mo>/</mo><msub><mi>p</mi><mn>2</mn></msub><mo>-</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><mo>|</mo><mo>|</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>23</mn><mo>)</mo></mrow></mrow></math>]]></maths>ρ为反演非奇异终端滑模控制参数;2-5)由于微陀螺仪的三个参数矩阵D,K,Ω未知,根据自适应控制理论,用估计值<img file="FDA0000454131740000025.GIF" wi="247" he="68" />分别替代式(20)中的参数矩阵D,K,Ω,并设计三个估计值的自适应算法,在线实时更新估计值,则式(20)所述的控制律u<sub>0</sub>调整为u′<sub>0</sub>:<maths num="0009"><![CDATA[<math><mrow><msup><msub><mi>u</mi><mn>0</mn></msub><mo>&prime;</mo></msup><mo>=</mo><mrow><mo>(</mo><mover><mi>D</mi><mo>^</mo></mover><mo>+</mo><mn>2</mn><mover><mi>&Omega;</mi><mo>^</mo></mover><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>e</mi><mn>2</mn></msub><mo>+</mo><msub><mi>&alpha;</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><mover><mi>K</mi><mo>^</mo></mover><mrow><mo>(</mo><msub><mi>e</mi><mn>1</mn></msub><mo>+</mo><msub><mi>q</mi><mi>r</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mover><mi>&alpha;</mi><mo>&CenterDot;</mo></mover><mn>1</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>25</mn><mo>)</mo></mrow></mrow></math>]]></maths>2-6)用所述步骤2-5)调整后的控制律u'<sub>0</sub>代替步骤2-4)中的控制律u<sub>0</sub>,带入到式(19)中,得到自适应反演非奇异终端滑模控制器的控制律φ':φ'=u'<sub>0</sub>+u<sub>1</sub>+u<sub>2</sub>+u<sub>3</sub>    (26)2-7)将自适应反演非奇异终端滑模控制器的控制律φ'作为微陀螺仪系统的控制输入u,带入微陀螺仪系统的数学模型中,实现对微陀螺仪系统的跟踪控制。
地址 213022 江苏省常州市新北区晋陵北路200号