发明名称 基于滑模补偿的微陀螺仪鲁棒神经网络控制系统及方法
摘要 本发明公开了一种基于滑模补偿的微陀螺仪鲁棒神经网络控制系统及方法,所述控制系统包括给定轨迹生成模块,滑模面定义模块,神经网络控制器、权值自适应机制模块,滑模补偿器,微陀螺仪系统,比例微分控制模块,第一加法器和第二加法器,基于该控制系统,建立基于滑模面的微陀螺仪动力学模型,设计控制器结构,设计RBF网络权值的更新算法,以实现微陀螺仪轨迹追踪。本发明的控制方法能够在线补偿微陀螺仪的未知动态特性以及噪声干扰的影响,使得微陀螺仪的振动轨迹完全跟踪上参考轨迹,同时提高了系统抗干扰的鲁棒性和可靠性,网络权值的更新算法基于Lyapunov稳定性理论设计,保证闭环系统的稳定性,本发明为微陀螺仪应用范围的扩展提供了有力基础。
申请公布号 CN103116275B 申请公布日期 2016.04.06
申请号 CN201310065126.2 申请日期 2013.03.01
申请人 河海大学常州校区 发明人 杨玉正;费峻涛
分类号 G05B13/02(2006.01)I 主分类号 G05B13/02(2006.01)I
代理机构 南京纵横知识产权代理有限公司 32224 代理人 董建林
主权项 一种基于滑模补偿的微陀螺仪鲁棒神经网络控制系统的控制方法,其特征在于:所述控制系统包括:给定轨迹生成模块101,用于输出微陀螺仪两轴振动的参考轨迹;滑模面定义模块102,用于接收追踪误差,并产生一个滑模面信号输出;神经网络控制器103,用于接收参考轨迹和追踪误差信号,并产生神经网络控制器输出;所述神经网络控制器103选用的神经网络结构为RBF神经网络,它包含三层结构:输入层,隐层和输出层,输入层接受系统中的可测量信号输入,隐层采用高斯基函数计算非线性映射后的输出,输出层通过加权各隐层节点的输出得到整个RBF神经网络的输出;权值自适应机制模块104,用于接收滑模面信号,并产生神经网络权值更新算法;常规比例微分控制模块105,用于接收滑模面信号,并产生比例微分控制输出;滑模补偿器106,用于接收滑模面信号,并产生滑模补偿信号输出;微陀螺仪系统107,被控对象数学模型,考虑了外界干扰的影响,输出振动轨迹的位置和速度信号;第一加法器108,用于把参考轨迹与微陀螺仪的位置和速度输出相减,并产生追踪误差输出;第二加法器109,用于接收神经网络输出信号,比例微分控制输出信号和滑模补偿信号,产生微陀螺仪的控制输入;所述控制方法,包括以下步骤,1)建立基于滑模面的微陀螺仪动力学模型;具体为:1‑1)考虑到制造误差和外界干扰作用,两轴微机械陀螺仪的动力学方程的向量形式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>M</mi><mover><mi>q</mi><mo>&CenterDot;&CenterDot;</mo></mover><mo>+</mo><mi>D</mi><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mo>+</mo><mi>K</mi><mi>q</mi><mo>=</mo><msub><mi>&tau;</mi><mi>u</mi></msub><mo>-</mo><mn>2</mn><mi>&Omega;</mi><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mo>+</mo><msub><mi>&tau;</mi><mi>d</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000913363830000011.GIF" wi="1420" he="71" /></maths>式中,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mi>q</mi><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><mi>x</mi></mtd></mtr><mtr><mtd><mi>y</mi></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></mtd><mtd><mrow><msub><mi>&tau;</mi><mi>u</mi></msub><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>u</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><msub><mi>u</mi><mi>y</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></mtd><mtd><mrow><msub><mi>&tau;</mi><mi>d</mi></msub><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>d</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><msub><mi>d</mi><mi>y</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></mtd><mtd><mrow><mi>M</mi><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><mi>m</mi></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>m</mi></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></mtd><mtd><mrow><mi>D</mi><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>d</mi><mrow><mi>x</mi><mi>x</mi></mrow></msub></mtd><mtd><msub><mi>d</mi><mrow><mi>x</mi><mi>y</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>d</mi><mrow><mi>x</mi><mi>y</mi></mrow></msub></mtd><mtd><msub><mi>d</mi><mrow><mi>y</mi><mi>y</mi></mrow></msub></mtd></mtr></mtable></mfenced></mrow></mtd></mtr></mtable><mo>,</mo></mrow>]]></math><img file="FDA0000913363830000012.GIF" wi="1445" he="174" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open = '' close = ''><mtable><mtr><mtd><mrow><mi>K</mi><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>k</mi><mrow><mi>x</mi><mi>x</mi></mrow></msub></mtd><mtd><msub><mi>k</mi><mrow><mi>x</mi><mi>y</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>k</mi><mrow><mi>x</mi><mi>y</mi></mrow></msub></mtd><mtd><msub><mi>k</mi><mrow><mi>y</mi><mi>y</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></mtd><mtd><mrow><mi>&Omega;</mi><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><msub><mi>m&Omega;</mi><mi>z</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>m&Omega;</mi><mi>z</mi></msub></mrow></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000913363830000013.GIF" wi="779" he="173" /></maths>m为质量块的质量;x,y分别为质量块沿驱动轴和感测轴的位置;d<sub>xx</sub>,d<sub>xy</sub>,d<sub>yy</sub>为微陀螺仪的阻尼系数,k<sub>xx</sub>,k<sub>xy</sub>,k<sub>yy</sub>为微陀螺仪的弹簧系数;Ω<sub>z</sub>是微陀螺仪工作环境中的角速度;u<sub>x</sub>,u<sub>y</sub>是控制输入;d<sub>x</sub>,d<sub>y</sub>是外界干扰作用;1‑2)通过给定轨迹生成模块101,输出质量块两轴振动的参考轨迹q<sub>d</sub>;1‑3)通过第一加法器108,将参考轨迹q<sub>d</sub>与微陀螺仪系统107的位置和速度输出相减,并产生追踪误差函数e(t);1‑4)通过滑模面定义模块102,输出滑模面函数s(t),<img file="FDA0000913363830000021.GIF" wi="366" he="70" />其中,t为时间,Λ=Λ<sup>T</sup>>0,为滑模面系数;1‑5)结合式(2),得到基于滑模面的微陀螺仪的动力学方程<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>M</mi><mover><mi>s</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mo>-</mo><mi>D</mi><mi>s</mi><mo>-</mo><msub><mi>&tau;</mi><mi>u</mi></msub><mo>+</mo><mi>f</mi><mrow><mo>(</mo><mi>g</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&tau;</mi><mi>d</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000913363830000029.GIF" wi="1637" he="67" /></maths>式中,<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>f</mi><mrow><mo>(</mo><mi>g</mi><mo>)</mo></mrow><mo>=</mo><mi>M</mi><mrow><mo>(</mo><msub><mover><mi>q</mi><mo>&CenterDot;&CenterDot;</mo></mover><mi>d</mi></msub><mo>+</mo><mi>&Lambda;</mi><mover><mi>e</mi><mo>&CenterDot;</mo></mover><mo>)</mo></mrow><mo>+</mo><mi>D</mi><mrow><mo>(</mo><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mi>d</mi></msub><mo>+</mo><mi>&Lambda;</mi><mi>e</mi><mo>)</mo></mrow><mo>+</mo><mi>K</mi><mi>q</mi><mo>+</mo><mn>2</mn><mi>&Omega;</mi><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000913363830000022.GIF" wi="1502" he="70" /></maths>f(g)表示未知的微陀螺仪函数,<img file="FDA0000913363830000023.GIF" wi="175" he="71" /><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>g</mi><mo>=</mo><msup><mfenced open = '[' close = ']'><mtable><mtr><mtd><msup><mi>e</mi><mi>T</mi></msup></mtd><mtd><msup><mover><mi>e</mi><mo>&CenterDot;</mo></mover><mi>T</mi></msup></mtd><mtd><msubsup><mi>q</mi><mi>d</mi><mi>T</mi></msubsup></mtd><mtd><msubsup><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mi>d</mi><mi>T</mi></msubsup></mtd><mtd><msubsup><mover><mi>q</mi><mo>&CenterDot;&CenterDot;</mo></mover><mi>d</mi><mi>T</mi></msubsup></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>,</mo></mrow>]]></math><img file="FDA0000913363830000024.GIF" wi="619" he="103" /></maths>为可检测到的信号,作为神经网络的输入;2)设计控制器结构;具体为:2‑1)神经网络控制器103接收参考轨迹和追踪误差信号,产生神经网络控制器输出信号<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mover><mi>f</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>g</mi><mo>)</mo></mrow><mo>,</mo><mover><mi>f</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>g</mi><mo>)</mo></mrow><mo>=</mo><msup><mover><mi>W</mi><mo>^</mo></mover><mi>T</mi></msup><mi>&phi;</mi><mrow><mo>(</mo><mi>g</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000913363830000025.GIF" wi="510" he="86" /></maths>2‑2)常规比例微分控制模块105接收滑模面信号,输出比例微分控制输出信号K<sub>v</sub>s,<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>K</mi><mi>v</mi></msub><mi>s</mi><mo>=</mo><msub><mi>K</mi><mi>v</mi></msub><mrow><mo>(</mo><mover><mi>e</mi><mo>&CenterDot;</mo></mover><mo>+</mo><mi>&Lambda;</mi><mi>e</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000913363830000026.GIF" wi="381" he="70" /></maths>2‑3)滑模补偿器106接收滑模面信号,输出滑模补偿信号(ε<sub>N</sub>+b<sub>d</sub>)sgn(s);2‑4)第二加法器109接收上述三种信号,产生微陀螺仪的控制输入τ<sub>u</sub>,<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mi>&tau;</mi><mi>u</mi></msub><mo>=</mo><msup><mover><mi>W</mi><mo>^</mo></mover><mi>T</mi></msup><mi>&phi;</mi><mrow><mo>(</mo><mi>g</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>K</mi><mi>v</mi></msub><mi>s</mi><mo>+</mo><mrow><mo>(</mo><msub><mi>&epsiv;</mi><mi>N</mi></msub><mo>+</mo><msub><mi>b</mi><mi>d</mi></msub><mo>)</mo></mrow><mi>sgn</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000913363830000027.GIF" wi="1636" he="85" /></maths>2‑5)将控制输入式(15)带入动力学方程式(5)中,得到闭环动力学方程为:<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><mi>M</mi><mover><mi>s</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mo>-</mo><mrow><mo>(</mo><msub><mi>K</mi><mi>v</mi></msub><mo>+</mo><mi>D</mi><mo>)</mo></mrow><mi>s</mi><mo>+</mo><msup><mover><mi>W</mi><mo>~</mo></mover><mi>T</mi></msup><mi>&phi;</mi><mrow><mo>(</mo><mi>g</mi><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><mi>&epsiv;</mi><mo>(</mo><mi>g</mi><mo>)</mo><mo>-</mo><msub><mi>&tau;</mi><mi>d</mi></msub><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><msub><mi>&epsiv;</mi><mi>N</mi></msub><mo>+</mo><msub><mi>b</mi><mi>d</mi></msub><mo>)</mo></mrow><mi>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>]]></math><img file="FDA0000913363830000028.GIF" wi="1612" he="79" /></maths>其中,<img file="FDA0000913363830000031.GIF" wi="260" he="70" />为权值估计误差,W<sup>*</sup>为最优权值,<img file="FDA0000913363830000032.GIF" wi="60" he="70" />为实时权值,ε(g)=f(g)‑W<sup>*T</sup>φ(g),为RBF神经网络逼近误差函数,φ(g)=[φ<sub>j</sub>(g)],φ<sub>j</sub>(g)为各隐层节点的输出向量,j=1,2,…n<sub>2</sub>,n<sub>2</sub>表示隐层节点个数,K<sub>v</sub>为状态反馈增益,<img file="FDA0000913363830000033.GIF" wi="286" he="78" />ε<sub>N</sub>为神经网络逼近误差的上界,b<sub>d</sub>为外界干扰的上界,sgn()是符号函数;3)设计RBF网络权值的更新算法,采用基于Lyapunov稳定性理论设计RBF网络权值的更新算法,具体为通过权值自适应机制模块104输出神经网络权值的更新算法:<img file="FDA0000913363830000034.GIF" wi="293" he="85" />将更新算法应用于RBF神经网络,确保追踪滑模面s(t)收敛到零,其中,F=F<sup>T</sup>>0,为权值自适应增益矩阵;所述Lyapunov候选函数L为:<maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><mi>L</mi><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>s</mi><mi>T</mi></msup><mi>M</mi><mi>s</mi><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>t</mi><mi>r</mi><mrow><mo>(</mo><msup><mover><mi>W</mi><mo>~</mo></mover><mi>T</mi></msup><msup><mi>F</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mover><mi>W</mi><mo>~</mo></mover><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000913363830000035.GIF" wi="574" he="119" /></maths>其中,<maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><mi>M</mi><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><mi>m</mi></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>m</mi></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000913363830000036.GIF" wi="308" he="151" /></maths>m为质量块的质量。
地址 213022 江苏省常州市新北区晋陵北路200号