[发明专利]基于加速度计输出增量的捷联惯性导航方法有效

专利信息
申请号: 201410076442.4 申请日: 2014-03-04
公开(公告)号: CN103954282A 公开(公告)日: 2014-07-30
发明(设计)人: 奔粤阳;杨晓龙;夏健钟;李敬春;阮双双;杨祥龙;刘适;尹冬寒;李明宇;赵维珩 申请(专利权)人: 哈尔滨工程大学
主分类号: G01C21/16 分类号: G01C21/16
代理公司: 暂无信息 代理人: 暂无信息
地址: 150001 黑龙江省哈尔滨市南岗区*** 国省代码: 黑龙江;23
权利要求书: 查看更多 说明书: 查看更多
摘要: 发明提供的是一种基于加速度计输出增量的捷联惯性导航方法。首先利用陀螺仪的输出计算出船舶的姿态信息和载体坐标系到地理坐标系之间的转换矩阵,然后利用解算得到的和载体坐标系到地理坐标系之间的转换矩阵将加速度计输出比力信息转换到地理坐标系,在地理坐标系上利用当前时刻的比力信息和上一时刻的比力信息得到当前时刻的加速度计输出增量,利用增量信息解算出舰船的加速度信息,然后利用计算得到的加速度信息解算出舰船的速度信息和位置信息。由于在计算加速度计输出增量时,利用当前时刻与上一时刻做差,将零位误差减掉,从而减小加速度计零位误差对系统的影响达到提高系统定位精度的目的。
搜索关键词: 基于 加速度计 输出 增量 惯性 导航 方法
【主权项】:
一种基于加速度计输出增量的捷联惯性导航方法,其特征是:步骤一:舰船上的捷联惯性导航设备启动后,对导航设备进行冷启动;步骤二:在导航系统开始工作之后,首先利用陀螺仪的输出进行捷联惯性导航系统姿态解算;步骤三:利用步骤二的解算姿态,将相邻两次的加速度计输出值fb(k‑1)和fb(k)从载体坐标系转换到导航坐标系上,得到fn(k‑1)和fn(k);步骤四:由步骤三得到的导航坐标系上的比力信息,得到相邻两时刻的导航解算速度方程为:<mrow><msup><mover><mi>v</mi><mo>&CenterDot;</mo></mover><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msup><mi>f</mi><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><mn>2</mn><msubsup><mi>&omega;</mi><mi>ie</mi><mi>n</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><msubsup><mi>&omega;</mi><mi>en</mi><mi>n</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>)</mo></mrow><msup><mi>v</mi><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><msup><mover><mi>v</mi><mo>&CenterDot;</mo></mover><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msup><mi>f</mi><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><msubsup><mrow><mn>2</mn><mi>&omega;</mi></mrow><mi>ie</mi><mi>n</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msubsup><mi>&omega;</mi><mi>en</mi><mi>n</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>&times;</mo><msup><mi>v</mi><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow>其中:分别为k‑1时刻和k时刻地球自转角速度在导航坐标系上的投影;分别为k‑1时刻和k时刻载体相对于地球的转动角速率在导航坐标系上的投影;vn(k‑1),、vn(k)分别为k‑1时刻和k时刻载体速度在导航坐标系上的投影;步骤五:由k‑1时刻到k时刻,的变化都是小量,忽略它们的变化,由步骤四中的式<mrow><msup><mover><mi>v</mi><mo>&CenterDot;</mo></mover><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msup><mi>f</mi><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><mn>2</mn><msubsup><mi>&omega;</mi><mi>ie</mi><mi>n</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><msubsup><mi>&omega;</mi><mi>en</mi><mi>n</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>)</mo></mrow><msup><mi>v</mi><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow>和式<mrow><msup><mover><mi>v</mi><mo>&CenterDot;</mo></mover><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msup><mi>f</mi><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><msubsup><mrow><mn>2</mn><mi>&omega;</mi></mrow><mi>ie</mi><mi>n</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msubsup><mi>&omega;</mi><mi>en</mi><mi>n</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>&times;</mo><msup><mi>v</mi><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow>做差得到:<mrow><mrow><mo>(</mo><msup><mover><mi>v</mi><mo>&CenterDot;</mo></mover><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msup><mover><mi>v</mi><mo>&CenterDot;</mo></mover><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>)</mo></mrow><mo>/</mo><mi>H</mi><mo>=</mo><mrow><mo>(</mo><msup><mi>f</mi><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msup><mi>f</mi><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>)</mo></mrow><mo>/</mo><mi>H</mi><mo>+</mo><mrow><mo>(</mo><msubsup><mrow><mn>2</mn><mi>&omega;</mi></mrow><mi>ie</mi><mi>n</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msubsup><mi>&omega;</mi><mi>en</mi><mi>n</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>&times;</mo><mrow><mo>(</mo><msup><mi>v</mi><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msup><mi>v</mi><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>)</mo></mrow><mo>/</mo><mi>H</mi></mrow>也即:<mrow><msup><mover><mi>a</mi><mo>&CenterDot;</mo></mover><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><msup><mi>f</mi><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msup><mi>f</mi><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>)</mo></mrow><mo>/</mo><mi>H</mi><mo>+</mo><mrow><mo>(</mo><msubsup><mrow><mn>2</mn><mi>&omega;</mi></mrow><mi>ie</mi><mi>n</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msubsup><mi>&omega;</mi><mi>en</mi><mi>n</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>&times;</mo><msup><mover><mi>v</mi><mo>&CenterDot;</mo></mover><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow>其中,由k‑1时刻到k时刻时间长度为H,解微分方程式<mrow><msup><mover><mi>a</mi><mo>&CenterDot;</mo></mover><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><msup><mi>f</mi><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msup><mi>f</mi><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>)</mo></mrow><mo>/</mo><mi>H</mi><mo>+</mo><mrow><mo>(</mo><msubsup><mrow><mn>2</mn><mi>&omega;</mi></mrow><mi>ie</mi><mi>n</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msubsup><mi>&omega;</mi><mi>en</mi><mi>n</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>&times;</mo><msup><mover><mi>v</mi><mo>&CenterDot;</mo></mover><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow>得到载体在k时刻的加速度信息;步骤六:在得到载体的加速度信息之后,对其进行一次积分运算,即:<mrow><msup><mi>v</mi><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Integral;</mo><msub><mi>t</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><msub><mi>t</mi><mi>k</mi></msub></munderover><msup><mi>a</mi><mi>n</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mi>dt</mi></mrow>得到载体在k时刻的速度;步骤七:在步骤六得到载体的速度信息之后,对速度进行一次积分,得到载体的位置信息,完成导航系统的导航定位。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

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

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

×

专利文献下载

说明:

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

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

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

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

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

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

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

钻瓜专利网在线咨询

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

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