发明名称 一种自主水下航行器动态抗饱和俯仰角控制方法
摘要 本发明提出一种自主水下航行器动态抗饱和俯仰角控制方法,方法中基于滑模控制设计了自适应律,在线估计系统参数及未知干扰的上界,不断地通过当前系统的运行状态与期望状态相比较,通过自适应律来改变控制器,使系统能够达到预期的性能指标。同时方法设计了抗饱和动态补偿器,用辅助变量w来减缓输入饱和对系统的影响,以保证控制系统有良好的性能。
申请公布号 CN104793625B 申请公布日期 2017.05.03
申请号 CN201510184979.7 申请日期 2015.04.17
申请人 西北工业大学 发明人 崔荣鑫;严卫生;张欣
分类号 G05D1/08(2006.01)I 主分类号 G05D1/08(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 陈星
主权项 一种自主水下航行器动态抗饱和俯仰角控制方法,其特征在于:包括以下步骤:步骤1:忽略水下航行器俯仰运动和横滚运动的耦合,建立水下航行器俯仰运动模型:<maths num="0001"><math><![CDATA[<mrow><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mo>=</mo><msub><mi>&omega;</mi><mi>z</mi></msub></mrow>]]></math><img file="FDA0001174592170000011.GIF" wi="126" he="69" /></maths><maths num="0002"><math><![CDATA[<mrow><msub><mover><mi>&omega;</mi><mo>&CenterDot;</mo></mover><mi>z</mi></msub><mo>=</mo><msub><mi>c</mi><mn>1</mn></msub><msub><mi>&omega;</mi><mi>z</mi></msub><mo>+</mo><msub><mi>c</mi><mn>2</mn></msub><mi>&alpha;</mi><mo>+</mo><msub><mi>c</mi><mn>3</mn></msub><mi>u</mi><mo>+</mo><mi>D</mi></mrow>]]></math><img file="FDA0001174592170000012.GIF" wi="486" he="63" /></maths><maths num="0003"><math><![CDATA[<mrow><mover><mi>&alpha;</mi><mo>&CenterDot;</mo></mover><mo>=</mo><msub><mi>c</mi><mn>4</mn></msub><msub><mi>&omega;</mi><mi>z</mi></msub><mo>+</mo><msub><mi>c</mi><mn>5</mn></msub><mi>&alpha;</mi><mo>+</mo><msub><mi>c</mi><mn>6</mn></msub><mi>u</mi></mrow>]]></math><img file="FDA0001174592170000013.GIF" wi="358" he="55" /></maths>y=θ其中θ是俯仰角,ω<sub>z</sub>是俯仰角速度,α是攻角,u是控制输入,D表示非线性项和外界未知干扰,y是输出,c<sub>1</sub>,c<sub>2</sub>是未知的需要估计的模型参数,c<sub>3</sub>是已知的控制增益,c<sub>4</sub>,c<sub>5</sub>,c<sub>6</sub>为设定的模型参数;控制输入u取决于控制器输出的舵角δ<sub>e</sub>:<maths num="0004"><math><![CDATA[<mrow><mi>u</mi><mo>=</mo><mi>s</mi><mi>a</mi><mi>t</mi><mrow><mo>(</mo><msub><mi>&delta;</mi><mi>e</mi></msub><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>u</mi><mi>m</mi></msub><mo>,</mo></mrow></mtd><mtd><mrow><msub><mi>&delta;</mi><mi>e</mi></msub><mo>&gt;</mo><msub><mi>u</mi><mi>m</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&delta;</mi><mi>e</mi></msub><mo>,</mo></mrow></mtd><mtd><mrow><mo>-</mo><msub><mi>u</mi><mi>m</mi></msub><mo>&lt;</mo><msub><mi>&delta;</mi><mi>e</mi></msub><mo>&lt;</mo><msub><mi>u</mi><mi>m</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><msub><mi>u</mi><mi>m</mi></msub><mo>,</mo></mrow></mtd><mtd><mrow><msub><mi>&delta;</mi><mi>e</mi></msub><mo>&lt;</mo><msub><mi>u</mi><mi>m</mi></msub></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001174592170000014.GIF" wi="698" he="223" /></maths>其中u<sub>m</sub>为控制器幅值;步骤2:对步骤1设计得到的水下航行器俯仰运动模型,采用以下控制律和自适应律进行控制;所述控制律为:<maths num="0005"><math><![CDATA[<mrow><msub><mi>&delta;</mi><mi>e</mi></msub><mo>=</mo><mfrac><mn>1</mn><msub><mi>c</mi><mn>3</mn></msub></mfrac><mo>&lsqb;</mo><mo>-</mo><mi>&tau;</mi><mrow><mo>(</mo><mi>s</mi><mo>-</mo><mi>w</mi><mo>)</mo></mrow><mo>-</mo><mi>k</mi><mo>*</mo><mi>s</mi><mi>g</mi><mi>n</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>-</mo><msup><mover><mi>h</mi><mo>^</mo></mover><mi>T</mi></msup><mi>&phi;</mi><mo>-</mo><mover><mi>&rho;</mi><mo>^</mo></mover><mo>*</mo><mi>s</mi><mi>g</mi><mi>n</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>+</mo><msub><mover><mi>&theta;</mi><mo>&CenterDot;&CenterDot;</mo></mover><mi>d</mi></msub><mo>-</mo><mi>&lambda;</mi><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>z</mi></msub><mo>-</mo><msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mi>d</mi></msub><mo>)</mo></mrow><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0001174592170000015.GIF" wi="1310" he="127" /></maths>其中s为滑模面,取<img file="FDA0001174592170000016.GIF" wi="218" he="46" />e为俯仰角误差e=θ‑θ<sub>d</sub>,θ<sub>d</sub>为期望俯仰角,λ为满足Hurwitz条件的正常数;sgn(s)表示s的符号;参数τ和k由设定的滑模面s的趋近律得到,滑模面s的趋近律为指数趋近律<img file="FDA0001174592170000017.GIF" wi="403" he="55" />且τ>0,k>0;w为设计的动态抗饱和补偿器状态,设计的动态抗饱和补偿器为:<maths num="0006"><math><![CDATA[<mrow><mover><mi>w</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mo>-</mo><mover><mi>k</mi><mo>&OverBar;</mo></mover><mi>w</mi><mo>-</mo><msup><mfrac><mrow><mo>|</mo><msub><mi>c</mi><mn>3</mn></msub><msup><mi>s</mi><mi>T</mi></msup><mi>&Delta;</mi><mi>u</mi><mo>|</mo><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo>(</mo><mi>&Delta;</mi><mi>u</mi><mo>)</mo></mrow></mrow><mrow><mo>|</mo><mo>|</mo><mi>w</mi><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup></mrow></mfrac><mn>2</mn></msup><mi>w</mi><mo>+</mo><mi>&Delta;</mi><mi>u</mi><mo>,</mo></mrow></mtd><mtd><mrow><mo>|</mo><mo>|</mo><mi>w</mi><mo>|</mo><mo>|</mo><mo>&GreaterEqual;</mo><mi>i</mi></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><mo>|</mo><mo>|</mo><mi>w</mi><mo>|</mo><mo>|</mo><mo>&lt;</mo><mi>i</mi></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001174592170000018.GIF" wi="949" he="303" /></maths>其中Δu=u‑δ<sub>e</sub>,s<sup>T</sup>表示s的转置,<img file="FDA0001174592170000019.GIF" wi="38" he="55" />表示给定的抗饱和参数,且满足<img file="FDA00011745921700000110.GIF" wi="290" he="111" />||·||表示范数,ι取值范围为[0,0.001];所述自适应律为:<maths num="0007"><math><![CDATA[<mrow><mover><mover><mi>h</mi><mo>^</mo></mover><mo>&CenterDot;</mo></mover><mo>=</mo><mo>-</mo><msub><mi>&epsiv;</mi><mn>1</mn></msub><mover><mi>h</mi><mo>^</mo></mover><mo>+</mo><msub><mi>k</mi><mn>1</mn></msub><mi>&phi;</mi><mi>s</mi></mrow>]]></math><img file="FDA0001174592170000021.GIF" wi="302" he="85" /></maths><maths num="0008"><math><![CDATA[<mrow><mover><mover><mi>&rho;</mi><mo>^</mo></mover><mo>&CenterDot;</mo></mover><mo>=</mo><mo>-</mo><msub><mi>&epsiv;</mi><mn>2</mn></msub><mover><mi>&rho;</mi><mo>^</mo></mover><mo>+</mo><msub><mi>k</mi><mn>2</mn></msub><mo>|</mo><mi>s</mi><mo>|</mo></mrow>]]></math><img file="FDA0001174592170000022.GIF" wi="325" he="78" /></maths>其中,h<sup>T</sup>=[c<sub>1</sub>,c<sub>2</sub>],φ=[ω<sub>z</sub>,α]<sup>T</sup>,<img file="FDA0001174592170000023.GIF" wi="35" he="63" />为h的估计值;ρ为非线性项和外界未知干扰D的上界,<img file="FDA0001174592170000024.GIF" wi="38" he="55" />为ρ的估计值,ε<sub>1</sub>,ε<sub>2</sub>,k<sub>1</sub>,k<sub>2</sub>均为给定的自适应律设计参数。
地址 710072 陕西省西安市友谊西路127号