发明名称 一种基于无迹卡尔曼滤波的极区环境下舰船大方位失准角传递对准方法
摘要 本发明公开了一种基于无迹卡尔曼滤波的极区环境下舰船大方位失准角传递对准方法。步骤一:完成子惯导系统的启动、预热准备,子惯导系统利用主惯导系统发送的导航参数完成一次装订;步骤二:子惯导系统进行惯导解算,同步采集主、子惯导系统在格网系下输出的速度和姿态信息,得到速度和姿态误差来构成量测量Z;步骤三:依据格网系下的导航力学编排,结合格网导航误差方程,采用“速度+姿态”的匹配方式,建立格网系下的系统状态方程和量测方程;步骤四:进行无迹卡尔曼滤波解算,估算出子惯导系统的姿态失准角、速度的状态估算值,完成传递对准。本发明解决了极区环境下大方位失准角的传递对准问题,提高了极区的传递对准精度和适用性。
申请公布号 CN105783943A 申请公布日期 2016.07.20
申请号 CN201610265830.6 申请日期 2016.04.26
申请人 哈尔滨工程大学 发明人 程建华;王通达;范晓亮;陈岱岱;管冬雪;康瑛瑶;蒋国桉;董楠楠;于东伟;张为
分类号 G01C25/00(2006.01)I 主分类号 G01C25/00(2006.01)I
代理机构 代理人
主权项 一种基于无迹卡尔曼滤波的极区环境下舰船大方位失准角传递对准方法,其特征在于:包括以下步骤:步骤一:完成子惯导系统的启动、预热准备,子惯导系统利用主惯导系统发送的导航参数完成一次装订,所述的导航参数包括速度、姿态矩阵和位置信息;步骤二:子惯导系统进行惯导解算,同步采集主、子惯导系统在格网系下输出的速度和姿态信息,得到速度和姿态误差来构成量测量Z:Z=[δV<sup>G</sup>φ<sup>G</sup>]<sup>T</sup>其中:δV<sup>G</sup>是格网系下的速度误差值;φ<sup>G</sup>是格网系下的姿态失准角;步骤三:在极区舰船存在大方位失准角的情况下,依据格网系下的导航力学编排,结合格网导航误差方程,采用“速度+姿态”的匹配方式,建立格网系下的系统状态方程和量测方程;状态变量为:<img file="FDA0000974984440000011.GIF" wi="637" he="94" />系统状态方程为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>&delta;</mi><msup><mover><mi>V</mi><mo>&CenterDot;</mo></mover><mi>G</mi></msup><mo>=</mo><mrow><mo>(</mo><msubsup><mi>C</mi><msup><mi>s</mi><mo>&prime;</mo></msup><mi>G</mi></msubsup><mo>-</mo><msubsup><mi>C</mi><mi>s</mi><mi>G</mi></msubsup><mo>)</mo></mrow><msup><mi>f</mi><mi>s</mi></msup><mo>-</mo><mrow><mo>(</mo><mn>2</mn><msubsup><mi>&omega;</mi><mrow><mi>i</mi><mi>e</mi></mrow><mi>G</mi></msubsup><mo>+</mo><msubsup><mi>&omega;</mi><mrow><mi>e</mi><mi>G</mi></mrow><mi>G</mi></msubsup><mo>)</mo></mrow><mo>&times;</mo><msup><mi>&delta;V</mi><mi>G</mi></msup></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mrow><mo>(</mo><mn>2</mn><msubsup><mi>&delta;&omega;</mi><mrow><mi>i</mi><mi>e</mi></mrow><mi>G</mi></msubsup><mo>+</mo><msubsup><mi>&delta;&omega;</mi><mrow><mi>e</mi><mi>G</mi></mrow><mi>G</mi></msubsup><mo>)</mo></mrow><mo>&times;</mo><msup><mi>V</mi><mi>G</mi></msup><mo>+</mo><msubsup><mi>C</mi><msup><mi>s</mi><mo>&prime;</mo></msup><mi>G</mi></msubsup><msubsup><mo>&dtri;</mo><mi>s</mi><mi>s</mi></msubsup><mo>+</mo><msubsup><mi>C</mi><msup><mi>s</mi><mo>&prime;</mo></msup><mi>G</mi></msubsup><msubsup><mo>&dtri;</mo><mi>w</mi><mi>s</mi></msubsup></mrow></mtd></mtr><mtr><mtd><mrow><msup><mover><mi>&phi;</mi><mo>&CenterDot;</mo></mover><mi>G</mi></msup><mo>=</mo><mrow><mo>(</mo><msubsup><mi>C</mi><mi>m</mi><mrow><msup><mi>s</mi><mo>&prime;</mo></msup></mrow></msubsup><mo>-</mo><msubsup><mi>C</mi><mi>m</mi><mi>s</mi></msubsup><mo>)</mo></mrow><msubsup><mi>&omega;</mi><mrow><mi>G</mi><mi>m</mi></mrow><mi>m</mi></msubsup><mo>+</mo><msub><mi>&omega;</mi><mi>f</mi></msub><mo>-</mo><msubsup><mi>C</mi><mi>s</mi><mi>G</mi></msubsup><msubsup><mi>&epsiv;</mi><mi>s</mi><mi>s</mi></msubsup><mo>-</mo><msubsup><mi>C</mi><mi>s</mi><mi>G</mi></msubsup><msubsup><mi>&epsiv;</mi><mi>w</mi><mi>s</mi></msubsup></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mover><mo>&dtri;</mo><mo>&CenterDot;</mo></mover><mi>s</mi><mi>s</mi></msubsup><mo>=</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mover><mi>&epsiv;</mi><mo>&CenterDot;</mo></mover><mi>s</mi><mi>s</mi></msubsup><mo>=</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mover><mi>&mu;</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000974984440000012.GIF" wi="892" he="519" /></maths>其中,<img file="FDA0000974984440000013.GIF" wi="101" he="55" />是速度误差的导数,<img file="FDA0000974984440000014.GIF" wi="61" he="70" />是姿态失准角的导数,<img file="FDA0000974984440000015.GIF" wi="61" he="70" />是加速度常值漂移的导数,<img file="FDA0000974984440000016.GIF" wi="46" he="63" />是陀螺常值漂移的导数,<img file="FDA0000974984440000017.GIF" wi="38" he="55" />是安装失调角的导数,<img file="FDA0000974984440000018.GIF" wi="62" he="69" />是子惯导的计算姿态矩阵,<img file="FDA0000974984440000019.GIF" wi="62" he="63" />是子惯导载体坐标系s到格网坐标系G间的转换矩阵,<img file="FDA00009749844400000110.GIF" wi="70" he="71" />是主惯导载体坐标系m到子惯导计算载体坐标系s′的转换矩阵,<img file="FDA00009749844400000111.GIF" wi="70" he="63" />是主惯导载体坐标系m到子惯导载体坐标系s的转换矩阵,f<sup>s</sup>是加速度计量测量,<img file="FDA00009749844400000112.GIF" wi="65" he="70" />是格网系下的地球自转速率,<img file="FDA00009749844400000113.GIF" wi="78" he="69" />是格网系下的位置速率,δV<sup>G</sup>是格网系下的速度误差,<img file="FDA00009749844400000114.GIF" wi="97" he="70" />是格网系下地球自转角速率的误差,<img file="FDA00009749844400000115.GIF" wi="99" he="63" />是格网系下位置速率误差,V<sup>G</sup>是格网系下的速度,<img file="FDA00009749844400000116.GIF" wi="62" he="63" />为加速度常值偏移,<img file="FDA00009749844400000117.GIF" wi="70" he="62" />为加速度计随机偏移,<img file="FDA00009749844400000118.GIF" wi="85" he="63" />为主惯导相对于格网系的角速率在载体坐标系s下的投影,<img file="FDA00009749844400000119.GIF" wi="54" he="61" />为陀螺仪常值漂移,<img file="FDA00009749844400000120.GIF" wi="54" he="62" />为陀螺仪随机漂移;量测方程为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><mi>Z</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mi>H</mi><mi>X</mi><mo>+</mo><mi>v</mi></mrow>]]></math><img file="FDA0000974984440000021.GIF" wi="245" he="55" /></maths>其中,Z是量测量,H是量测矩阵,v是量测噪声,其中H如下所示:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>H</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>I</mi><mrow><mn>2</mn><mo>&times;</mo><mn>2</mn></mrow></msub></mtd><mtd><msub><mn>0</mn><mrow><mn>2</mn><mo>&times;</mo><mn>2</mn></mrow></msub></mtd><mtd><msub><mn>0</mn><mrow><mn>2</mn><mo>&times;</mo><mn>2</mn></mrow></msub></mtd><mtd><msub><mn>0</mn><mrow><mn>2</mn><mo>&times;</mo><mn>2</mn></mrow></msub></mtd><mtd><msub><mn>0</mn><mrow><mn>2</mn><mo>&times;</mo><mn>2</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mn>0</mn><mrow><mn>3</mn><mo>&times;</mo><mn>3</mn></mrow></msub></mtd><mtd><msub><mi>I</mi><mrow><mn>3</mn><mo>&times;</mo><mn>3</mn></mrow></msub></mtd><mtd><msub><mn>0</mn><mrow><mn>3</mn><mo>&times;</mo><mn>3</mn></mrow></msub></mtd><mtd><msub><mn>0</mn><mrow><mn>3</mn><mo>&times;</mo><mn>3</mn></mrow></msub></mtd><mtd><mrow><mo>-</mo><msubsup><mi>C</mi><msup><mi>s</mi><mo>&prime;</mo></msup><mi>G</mi></msubsup></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000974984440000022.GIF" wi="758" he="143" /></maths>步骤四:利用设计的格网坐标系下状态方程、量测方程和构造的量测量,进行无迹卡尔曼滤波解算,估算出子惯导系统的姿态失准角、速度的状态估算值,完成传递对准。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室