[发明专利]一种基于无迹卡尔曼滤波的极区环境下舰船大方位失准角传递对准方法在审

专利信息
申请号: 201610265830.6 申请日: 2016-04-26
公开(公告)号: CN105783943A 公开(公告)日: 2016-07-20
发明(设计)人: 程建华;王通达;范晓亮;陈岱岱;管冬雪;康瑛瑶;蒋国桉;董楠楠;于东伟;张为 申请(专利权)人: 哈尔滨工程大学
主分类号: G01C25/00 分类号: G01C25/00
代理公司: 暂无信息 代理人: 暂无信息
地址: 150001 黑龙江省哈尔滨市南岗区*** 国省代码: 黑龙江;23
权利要求书: 查看更多 说明书: 查看更多
摘要: 发明公开了一种基于无迹卡尔曼滤波的极区环境下舰船大方位失准角传递对准方法。步骤一:完成子惯导系统的启动、预热准备,子惯导系统利用主惯导系统发送的导航参数完成一次装订;步骤二:子惯导系统进行惯导解算,同步采集主、子惯导系统在格网系下输出的速度和姿态信息,得到速度和姿态误差来构成量测量Z;步骤三:依据格网系下的导航力学编排,结合格网导航误差方程,采用“速度+姿态”的匹配方式,建立格网系下的系统状态方程和量测方程;步骤四:进行无迹卡尔曼滤波解算,估算出子惯导系统的姿态失准角、速度的状态估算值,完成传递对准。本发明解决了极区环境下大方位失准角的传递对准问题,提高了极区的传递对准精度和适用性。
搜索关键词: 一种 基于 卡尔 滤波 环境 舰船 方位 失准 传递 对准 方法
【主权项】:
一种基于无迹卡尔曼滤波的极区环境下舰船大方位失准角传递对准方法,其特征在于:包括以下步骤:步骤一:完成子惯导系统的启动、预热准备,子惯导系统利用主惯导系统发送的导航参数完成一次装订,所述的导航参数包括速度、姿态矩阵和位置信息;步骤二:子惯导系统进行惯导解算,同步采集主、子惯导系统在格网系下输出的速度和姿态信息,得到速度和姿态误差来构成量测量Z:Z=[δVGφG]T其中:δVG是格网系下的速度误差值;φG是格网系下的姿态失准角;步骤三:在极区舰船存在大方位失准角的情况下,依据格网系下的导航力学编排,结合格网导航误差方程,采用“速度+姿态”的匹配方式,建立格网系下的系统状态方程和量测方程;状态变量为:系统状态方程为:<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>其中,是速度误差的导数,是姿态失准角的导数,是加速度常值漂移的导数,是陀螺常值漂移的导数,是安装失调角的导数,是子惯导的计算姿态矩阵,是子惯导载体坐标系s到格网坐标系G间的转换矩阵,是主惯导载体坐标系m到子惯导计算载体坐标系s′的转换矩阵,是主惯导载体坐标系m到子惯导载体坐标系s的转换矩阵,fs是加速度计量测量,是格网系下的地球自转速率,是格网系下的位置速率,δVG是格网系下的速度误差,是格网系下地球自转角速率的误差,是格网系下位置速率误差,VG是格网系下的速度,为加速度常值偏移,为加速度计随机偏移,为主惯导相对于格网系的角速率在载体坐标系s下的投影,为陀螺仪常值漂移,为陀螺仪随机漂移;量测方程为:<mrow><mover><mi>Z</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mi>H</mi><mi>X</mi><mo>+</mo><mi>v</mi></mrow>其中,Z是量测量,H是量测矩阵,v是量测噪声,其中H如下所示:<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>步骤四:利用设计的格网坐标系下状态方程、量测方程和构造的量测量,进行无迹卡尔曼滤波解算,估算出子惯导系统的姿态失准角、速度的状态估算值,完成传递对准。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于哈尔滨工程大学,未经哈尔滨工程大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/patent/201610265830.6/,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top