[发明专利]一种单轴旋转的捷联惯导系统初始对准方法有效
申请号: | 200910071734.8 | 申请日: | 2009-04-08 |
公开(公告)号: | CN101514900A | 公开(公告)日: | 2009-08-26 |
发明(设计)人: | 孙枫;孙伟;李国强;薛媛媛;奔粤阳;高伟;赵彦雷;吕玉红;李雪峰 | 申请(专利权)人: | 哈尔滨工程大学 |
主分类号: | G01C21/16 | 分类号: | G01C21/16;G01C25/00 |
代理公司: | 暂无信息 | 代理人: | 暂无信息 |
地址: | 150001黑龙江省哈尔滨市南岗区南通*** | 国省代码: | 黑龙江;23 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明提供的是一种单轴旋转的捷联惯导系统初始对准方法。对于载体静止状态下的捷联惯性导航系统,在其采集陀螺仪输出和加速度计输出信息完成粗对准的基础之上,建立载体坐标系和计算地理坐标系之间的转换矩阵;建立以速度误差为状态变量的卡尔曼滤波状态方程及速度误差为量测量的量测方程;通过卡尔曼滤波技术估计出载体失准角并反馈到系统中完成系统的初始对准。本发明能克服地理坐标系等效陀螺漂移对方位失准角估算精度的影响,提高对准精度。 | ||
搜索关键词: | 一种 旋转 捷联惯导 系统 初始 对准 方法 | ||
【主权项】:
1、一种单轴旋转的捷联惯导系统初始对准方法,其特征在于包括以下步骤:(1)利用全球定位系统GPS确定载体的初始位置参数,并装订至导航计算机中;(2)采集光纤陀螺捷联惯性导航系统的光纤陀螺仪和石英加速度计输出的数据;(3)对采集到的光纤陀螺仪和石英加速度计的数据进行处理,采用解析法来完成系统的粗对准,初步确定载体的姿态信息;(4)惯性测量单元绕载体坐标系oyb轴正向旋转45度,确定IMU坐标系与载体坐标系之间的初始相对位置;(5)确定两坐标系相对初始位置关系后,惯性测量单元绕载体坐标系方位轴ozb正向以角速度ω=6°/s连续转动;(6)建立以速度误差为状态变量的卡尔曼滤波状态方程及速度误差为量测量的量测方程;1)建立卡尔曼滤波的状态方程:用一阶线性微分方程来描述旋转捷联惯导系统的状态误差:X · ( t ) = A ( t ) X ( t ) + B ( t ) W ( t ) ]]> 其中X(t)为t时刻系统的状态向量;A(t)和B(t)分别为系统的状态矩阵和噪声矩阵;W(t)为系统噪声向量;系统的状态向量为:系统的白噪声向量为:W=[ax ay ωz ωy ωz 0 0 0 0 0]T其中δVe、δVn分别表示东向、北向的速度误差;分别为IMU坐标系oxs、oys轴加速度计零偏;εx、εy、εz分别为IMU坐标系oxs、oys、ozs轴陀螺的常值漂移;ax、ay分别为IMU坐标系oxs、oys轴加速度计的白噪声误差;ωx、ωy、ωz分别为IMU坐标系oxs、oys、ozs轴陀螺的白噪声误差;系统的状态转移矩阵为:A ( t ) = F 5 × 5 T 5 × 5 0 5 × 5 0 5 × 5 ]]> 其中:F 5 × 5 = F 2 × 2 F 2 × 3 F 3 × 2 F 3 × 3 ]]>F 2 × 2 = V n tan L R n 2 ω ie sin L + V e tan L R n - 2 ( ω ie sin L + V e tan L R n ) 0 ]]>F 2 × 3 = 0 - f u f n f u 0 - f e ]]>F 3 × 2 = 0 - 1 R m 1 R n 0 tan L R n 0 ]]>F 3 × 3 = 0 ω ie sin L + V e tan L R n - ( ω ie cos L + V e R n ) - ( ω ie sin L + V e tan L R n ) 0 - V n R m ω ie cos L + V e R n V n R m 0 ]]> Ve、Vn分别表示东向、北向的速度;ωie表示地球自转角速度;Rm、Rn分别表示地球子午、卯酉曲率半径;L表示当地纬度;fe、fn、fu分别表示为导航坐标系下东向、北向、天向的比力;令捷联矩阵Tsn′为:T s n ′ = T b n ′ T s b = T 11 T 12 T 13 T 21 T 22 T 23 T 31 T 32 T 33 ]]> 则T 5 × 5 = T 11 T 12 0 0 0 T 21 T 22 0 0 0 0 0 T 11 T 12 T 13 0 0 T 21 T 22 T 23 0 0 T 31 T 32 T 33 ; ]]> 2)建立卡尔曼滤波的量测方程:用一阶线性微分方程来描述旋转捷联惯导系统的量测方程如下:Z(t)=H(t)X(t)+v(t)其中:Z(t)表示t时刻系统的量测向量;H(t)表示系统的量测矩阵;v(t)表示系统的量测噪声;系统量测矩阵为:H ( t ) = 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 ]]> 量测量为水平速度误差:v(t)=[ve vn]T;(7)利用步骤(6)估计出平台失准角,并在精对准结束时用它来修正系统的捷联姿态矩阵,完成精确的初始对准。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于哈尔滨工程大学,未经哈尔滨工程大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/200910071734.8/,转载请声明来源钻瓜专利网。