主权项 |
一种AUV自主寻底控制方法,其特征在于,包括以下步骤:AUV无动力下潜至预定深度,并抛载下潜压铁;AUV在预定深度悬停,根据姿态传感器的信息调整AUV姿态;AUV从预定深度开始进行垂直寻底操作;AUV到达预定高度后,开始定高航行;所述垂直寻底操作采用自适应滤波算法,具体为:离散系统的状态方程为:<maths num="0001"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>s</mi><mi>k</mi></msub><mo>=</mo><msub><mi>As</mi><mi>k</mi></msub><mo>+</mo><msub><mi>Bu</mi><mi>k</mi></msub><mo>+</mo><msub><mi>w</mi><mi>k</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>z</mi><mi>k</mi></msub><mo>=</mo><msub><mi>Hs</mi><mi>k</mi></msub><mo>+</mo><msub><mi>v</mi><mi>k</mi></msub></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001144789970000011.GIF" wi="1230" he="136" /></maths>初始化:<maths num="0002"><math><![CDATA[<mrow><msub><mi>s</mi><mrow><mn>0</mn><mo>|</mo><mn>0</mn></mrow></msub><mo>=</mo><msub><mover><mi>h</mi><mo>‾</mo></mover><mn>0</mn></msub></mrow>]]></math><img file="FDA0001144789970000012.GIF" wi="166" he="86" /></maths>P<sub>0|0</sub>=10<sup>‑6</sup>*I预测:<maths num="0003"><math><![CDATA[<mrow><msub><mover><mi>s</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><mi>A</mi><msub><mover><mi>s</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>-</mo><mn>1</mn><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>Bu</mi><mi>k</mi></msub></mrow>]]></math><img file="FDA0001144789970000013.GIF" wi="438" he="87" /></maths>P<sub>k|k‑1</sub>=AP<sub>k‑1|k‑1</sub>A<sup>T</sup>+BQ<sub>k‑1</sub>B<sup>T</sup>更新:K<sub>k</sub>=P<sub>k|k‑1</sub>H<sup>T</sup>(HP<sub>k|k‑1</sub>H<sup>T</sup>+R<sub>k</sub>)<sup>‑1</sup><maths num="0004"><math><![CDATA[<mrow><msub><mover><mi>s</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>|</mo><mi>k</mi></mrow></msub><mo>=</mo><msub><mover><mi>s</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>K</mi><mi>k</mi></msub><mrow><mo>[</mo><mrow><msub><mi>z</mi><mi>k</mi></msub><mo>-</mo><mi>H</mi><msub><mover><mi>s</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow><mo>]</mo></mrow></mrow>]]></math><img file="FDA0001144789970000014.GIF" wi="614" he="126" /></maths>P<sub>k|k</sub>=(I‑K<sub>k</sub>)HP<sub>k|k‑1</sub>其中,s<sub>k</sub>=h<sub>k</sub>,表示k时刻载体所在的对底高度,其初始化由深度信息结合实地海深估算得到;u<sub>k</sub>=V<sub>k</sub>Δt,V<sub>k</sub>表示k时刻载体的垂向速度,由加速度积分得到,Δt表示速度间隔时间;Q<sub>k</sub>=diag(w<sub>k</sub><sup>2</sup>),R<sub>k</sub>=diag(v<sub>k</sub><sup>2</sup>),I是单位矩阵;z<sub>k</sub>是k时刻多普勒的测底高度值;w<sub>k</sub>是过程驱动噪声,v<sub>k</sub>为观测噪声,Q<sub>k</sub>、R<sub>k</sub>分别为过程噪声和观测噪声的协方差矩阵。 |