发明名称 微陀螺仪基于角速度估计的自适应模糊滑模控制方法
摘要 本发明提供了一种基于角速度估计的自适应模糊滑模控制方法,并将其应用于三轴微陀螺仪控制系统中,可以实现微陀螺仪系统的轨迹跟踪控制和参数估计。本发明基于微陀螺仪的状态空间方程进行控制器设计,首先,在滑模控制器的基础上设计了一种新型的自适应辨识方法,在线实时更新微陀螺仪的角速度和其它系统参数的估计值;另外利用自适应模糊系统来调节滑模控制切换项中的固定增益,对固定增益的上界值进行模糊逼近,可以减弱滑模控制中出现的抖振现象。应用本发明可以得到良好的跟踪性能,而且对参数变化和外部扰动具有鲁棒性。
申请公布号 CN103336435A 申请公布日期 2013.10.02
申请号 CN201310243814.3 申请日期 2013.06.19
申请人 河海大学常州校区 发明人 辛明缘;费峻涛
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 南京纵横知识产权代理有限公司 32224 代理人 董建林
主权项 1.微陀螺仪基于角速度估计的自适应模糊滑模控制方法,其特征在于,包括以下步骤1)建立微陀螺仪无量纲动态方程;2)建立微陀螺仪的参考模型;3)设计基于角速度估计的自适应模糊滑模控制器,具体为3-1)设计等效控制器,对滑模切换面S求导,令<img file="FDA00003366883700011.GIF" wi="148" he="71" />不考虑不确定性和外加干扰,得到等效控制力u<sub>eq</sub>,u<sub>eq</sub>=K<sup>*T</sup>X-(CB)<sup>-1</sup>CA<sub>m</sub>e其中,<maths num="0001"><![CDATA[<math><mrow><msup><mi>K</mi><mrow><mo>*</mo><mi>T</mi></mrow></msup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msup><msub><mi>&omega;</mi><mi>x</mi></msub><mn>2</mn></msup><mo>-</mo><msup><msub><mi>&omega;</mi><mn>1</mn></msub><mn>2</mn></msup></mtd><mtd><msub><mi>d</mi><mi>xx</mi></msub></mtd><mtd><msub><mi>&omega;</mi><mi>xy</mi></msub></mtd><mtd><msub><mi>d</mi><mi>xy</mi></msub><mo>-</mo><msub><mrow><mn>2</mn><mi>&Omega;</mi></mrow><mi>z</mi></msub></mtd><mtd><msub><mi>&omega;</mi><mi>xz</mi></msub></mtd><mtd><msub><mi>d</mi><mi>xz</mi></msub><mo>+</mo><msub><mrow><mn>2</mn><mi>&Omega;</mi></mrow><mi>y</mi></msub></mtd></mtr><mtr><mtd><msub><mi>&omega;</mi><mi>xy</mi></msub></mtd><mtd><msub><mi>d</mi><mi>xy</mi></msub><mo>+</mo><msub><mrow><mn>2</mn><mi>&Omega;</mi></mrow><mi>z</mi></msub></mtd><mtd><msup><msub><mi>&omega;</mi><mi>y</mi></msub><mn>2</mn></msup><mo>-</mo><msup><msub><mi>&omega;</mi><mn>2</mn></msub><mn>2</mn></msup></mtd><mtd><msub><mi>d</mi><mi>yy</mi></msub></mtd><mtd><msub><mi>&omega;</mi><mi>yz</mi></msub></mtd><mtd><msub><mi>d</mi><mi>yz</mi></msub><mo>-</mo><msub><mrow><mn>2</mn><mi>&Omega;</mi></mrow><mi>x</mi></msub></mtd></mtr><mtr><mtd><msub><mi>&omega;</mi><mi>xz</mi></msub></mtd><mtd><msub><mi>d</mi><mi>xz</mi></msub><mo>-</mo><msub><mrow><mn>2</mn><mi>&Omega;</mi></mrow><mi>y</mi></msub></mtd><mtd><msub><mi>&omega;</mi><mi>yz</mi></msub></mtd><mtd><mrow><mo>(</mo><msub><mi>d</mi><mi>yz</mi></msub><mo>+</mo><msub><mrow><mn>2</mn><mi>&Omega;</mi></mrow><mi>x</mi></msub><mo>)</mo></mrow></mtd><mtd><msup><msub><mi>&omega;</mi><mi>z</mi></msub><mn>2</mn></msup><mo>-</mo><msup><msub><mi>&omega;</mi><mn>3</mn></msub><mn>2</mn></msup></mtd><mtd><msub><mi>d</mi><mi>zz</mi></msub></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>为一个常数矩阵,满足A+BK<sup>*T</sup>=A<sub>m</sub>,X为微陀螺仪的轨迹状态,X<sub>m</sub>为参考模型的轨迹状态e为跟踪误差,满足e=X-X<sub>m</sub>S为滑模切换面函数,满足S=Ce,C为滑模系数A和B为微陀螺仪状态方程的状态参量,A<sub>m</sub>为参考模型状态方程的状态参量;3-2)考虑到所述等效控制力u<sub>eq</sub>中,K<sup>*</sup>未知,用K<sup>*</sup>的估计值K代替K<sup>*</sup>,自适应调节K的参数,并增加切换控制器,实现对不确定性和外加干扰的鲁棒控制得到自适应滑模控制律φ,φ=u<sub>eq</sub>+u<sub>sw</sub>=K<sup>T</sup>X-(CB)<sup>-1</sup>CA<sub>m</sub>e-fsgn(S)其中,u<sub>eq</sub>=K<sup>T</sup>X-(CB)<sup>-1</sup>CA<sub>m</sub>e为等效控制器的输出;u<sub>sw</sub>=-fsgn(S)为切换控制器的输出,实现对不确定性和外加干扰的鲁棒控制;f为固定增益,是系统本身的不确定性和外加干扰总量f<sub>m</sub>的上界值;3-3)基于Lyapunov函数,设计等效控制器中的自适应参数K的自适应律,实现系统参数和角速度的自适应辨识,保证整个闭环系统的全局渐进稳定性,Lyapunov函数V<sub>1</sub>为,<maths num="0002"><![CDATA[<math><mrow><msub><mi>V</mi><mn>1</mn></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>S</mi><mi>T</mi></msup><mi>S</mi><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>tr</mi><mo>[</mo><mover><mi>K</mi><mo>~</mo></mover><msup><mi>M</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><msup><mover><mi>K</mi><mo>~</mo></mover><mi>T</mi></msup><mo>]</mo></mrow></math>]]></maths>自适应律<img file="FDA00003366883700014.GIF" wi="128" he="83" />为,<maths num="0003"><![CDATA[<math><mrow><msup><mover><mi>K</mi><mo>.</mo></mover><mi>T</mi></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><msup><mi>MB</mi><mi>T</mi></msup><msup><mi>C</mi><mi>T</mi></msup><msup><mi>SX</mi><mi>T</mi></msup></mrow></math>]]></maths>其中,估计误差<img file="FDA00003366883700021.GIF" wi="57" he="72" />满足<img file="FDA00003366883700022.GIF" wi="262" he="77" />M为对称正定矩阵;3-4)由于f<sub>m</sub>未知,用自适应模糊系统的输出<img file="FDA00003366883700023.GIF" wi="30" he="71" />逼近切换控制器的固定增益值f,这样,切换控制器的输出变为<img file="FDA00003366883700024.GIF" wi="332" he="87" />则得到自适应模糊滑模控制律<img file="FDA00003366883700025.GIF" wi="38" he="59" />为:<img file="FDA00003366883700026.GIF" wi="968" he="91" />其中,<maths num="0004"><![CDATA[<math><mrow><msub><mover><mi>f</mi><mo>^</mo></mover><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></msubsup><msub><mi>&alpha;</mi><mi>ij</mi></msub><msub><mi>&mu;</mi><mi>A</mi></msub><mrow><mo>(</mo><msub><mi>s</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></msubsup><msub><mi>&mu;</mi><mi>A</mi></msub><mrow><mo>(</mo><msub><mi>s</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>=</mo><msup><msub><mover><mi>a</mi><mo>^</mo></mover><mi>i</mi></msub><mi>T</mi></msup><msub><mi>&xi;</mi><mi>i</mi></msub><mo>,</mo></mrow></math>]]></maths><img file="FDA00003366883700028.GIF" wi="49" he="86" />为<img file="FDA00003366883700029.GIF" wi="45" he="72" />的子变量,<img file="FDA000033668837000210.GIF" wi="48" he="76" />为可变参数,ξ<sub>i</sub>为模糊向量,<maths num="0005"><![CDATA[<math><mrow><msub><mover><mi>a</mi><mo>^</mo></mover><mi>i</mi></msub><mo>=</mo><msup><mrow><mo>[</mo><msub><mi>&alpha;</mi><mrow><mi>i</mi><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>&alpha;</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>&alpha;</mi><mi>im</mi></msub><mo>]</mo></mrow><mi>T</mi></msup><mo>,</mo></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><msub><mi>&xi;</mi><mi>i</mi></msub><mo>=</mo><msup><mrow><mo>[</mo><msub><mi>&xi;</mi><mrow><mi>i</mi><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>&xi;</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>&xi;</mi><mi>im</mi></msub><mo>]</mo></mrow><mi>T</mi></msup><mo>,</mo><mi>S</mi><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>S</mi><mn>1</mn></msub></mtd><mtd><msub><mi>S</mi><mn>2</mn></msub></mtd><mtd><msub><mi>S</mi><mn>3</mn></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>,</mo></mrow></math>]]></maths>μ<sub>A</sub>(s<sub>i</sub>),i=1,2,3为滑模切换面函数的隶属函数,sgn(·)为符号函数,上限m是指隶属函数的个数;3-5)基于Lyapunov函数,设计切换控制器中的可变参数<img file="FDA000033668837000213.GIF" wi="40" he="78" />的自适应律<img file="FDA000033668837000214.GIF" wi="81" he="88" />实现对外部扰动的鲁棒性,保证整个闭环系统的全局渐进稳定性,Lyapunov函数V<sub>2</sub>为,<maths num="0007"><![CDATA[<math><mrow><msub><mi>V</mi><mn>2</mn></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>S</mi><mi>T</mi></msup><mi>S</mi><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>tr</mi><mo>[</mo><mover><mi>K</mi><mo>~</mo></mover><msup><mi>M</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><msup><mover><mi>K</mi><mo>~</mo></mover><mi>T</mi></msup><mo>]</mo><mo>+</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>&eta;</mi></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>3</mn></munderover><msup><msub><mover><mi>a</mi><mo>~</mo></mover><mi>i</mi></msub><mi>T</mi></msup><msub><mover><mi>a</mi><mo>~</mo></mover><mi>i</mi></msub></mrow></math>]]></maths>自适应律<img file="FDA000033668837000216.GIF" wi="50" he="84" />为,<maths num="0008"><![CDATA[<math><mrow><msub><mover><mi>&alpha;</mi><mover><mo>^</mo><mo>.</mo></mover></mover><mi>i</mi></msub><mo>=</mo><mi>&eta;</mi><msub><mi>&lambda;</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub><mo>|</mo><msub><mi>s</mi><mi>i</mi></msub><mo>|</mo><msub><mi>&xi;</mi><mi>i</mi></msub></mrow></math>]]></maths>其中,η,λ<sub>i2</sub>为正常数,可变参数的估计误差<img file="FDA000033668837000218.GIF" wi="48" he="74" />满足<img file="FDA000033668837000219.GIF" wi="272" he="80" />α<sub>i</sub><sup>*</sup>为最优可变参数。
地址 213022 江苏省常州市新北区晋陵北路200号