[发明专利]一种基于双IMU/DGPS组合的相对姿态测量实时动态滤波方法有效
申请号: | 201110385043.2 | 申请日: | 2011-11-28 |
公开(公告)号: | CN102506857A | 公开(公告)日: | 2012-06-20 |
发明(设计)人: | 芦佳振;张春熹;李婕;李保国 | 申请(专利权)人: | 北京航空航天大学 |
主分类号: | G01C21/10 | 分类号: | G01C21/10;G01S19/41 |
代理公司: | 北京永创新实专利事务所 11121 | 代理人: | 周长琪 |
地址: | 100191*** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明公开基于双IMU/DGPS组合的相对姿态测量实时动态滤波方法,采用双光纤捷联惯导系统,通过捷联惯导实时解算,得到主惯导与从惯导系统的导航信息;判断DGPS信息是否更新,产生两种情况:当更新时,主惯导与从惯导系统进行滤波修正,建立组合导航滤波器的量测方程;当没有更新时,使用主惯导系统对从惯导系统进行滤波修正,建立组合导航滤波器的量测方程;两种情况得到的组合滤波方程通过离散化,建立离散型卡尔曼滤波器的递推方程,解算得到主惯导与从惯导系统的俯仰、横滚和航向角;随后进行相对姿态矩阵解算,得到主惯导与从惯导相对姿态角的主值。本发明提高了导航系统的稳定性,且可实时输出导航系统在测量过程中的速度、位置、姿态信息,测量范围广。 | ||
搜索关键词: | 一种 基于 imu dgps 组合 相对 姿态 测量 实时 动态 滤波 方法 | ||
【主权项】:
1.一种基于双IMU/DGPS组合的相对姿态测量实时动态滤波方法,其特征在于:采用主惯导系统与从惯导系统两个光纤捷联惯导系统与一个DGPS进行组合,具体滤波方法由下述步骤来完成:步骤1:主惯导系统和从惯导系统根据各自的惯性测量单元测得的载体加速度数据和角速度数据,分别通过捷联惯导实时解算,得到主惯导系统和从惯导系统的导航信息;步骤2:判断DGPS信息是否有更新;当DGPS信息更新时,进入步骤3;当DGPS信息丢失或未更新时,则进入步骤4。步骤3:根据由DGPS获取的载体的位置、速度等信息,与步骤一中得到的主惯导系统与从惯导系统的导航信息,分别对主惯导系统与从惯导系统进行滤波修正,建立组合导航滤波器的量测方程,随后进入步骤5;步骤4:使用主惯导系统的速度、位置信息对从惯导系统进行滤波修正,建立双惯导组合滤波器方程,随后进入步骤5;步骤5:将惯导系统组合滤波器量测方程和状态方程离散化,建立离散型卡尔曼滤波器的递推方程;步骤6:通过步骤5中建立的离散型卡尔曼滤波器的递推方程,进行Kalman滤波解算,得到主惯导系统与从惯导系统的俯仰、横滚和航向角;步骤7:根据步骤6得到的主惯导系统的俯仰、横滚和航向角ψAE、ψAN和ψAU,与从惯导系统的俯仰、横滚和航向角ψBE、ψBM和ψBU,建立由地理坐标系到主惯导系统载体坐标系的转换矩阵为
由地理坐标系到从惯导系统载体坐标系的转换矩阵为
则:C A n = cos ψ AN cos ψ AU - sin ψ AN sin ψ AE sin ψ AU - cos ψ AE sin ψ AU sin ψ AN cos ψ AU + cos ψ AN sin ψ AE sin ψ AU cos ψ AN sin ψ AU + sin ψ AN sin ψ AE cos ψ AU cos ψ AE cos ψ AU sin ψ AN sin ψ AU - cos ψ AN sin ψ AE cos ψ AU - sin ψ AN cos ψ AE sin ψ AE cos ψ AN cos ψ AE ]]>C B n = cos ψ BN cos ψ BU - sin ψ BN sin ψ BE sin ψ BU - cos ψ BE sin ψ BU sin ψ BN cos ψ BU + cos ψ BN sin ψ BE sin ψ BU cos ψ BN sin ψ BU + sin ψ BN sin ψ BE cos ψ BU cos ψ BE cos ψ BU sin ψ BN sin ψ BU - cos ψ BN sin ψ BE cos ψ BU - sin ψ BN cos ψ BE sin ψ BE cos ψ BN cos ψ BE ]]> 得到由主惯导系统载体坐标系到从惯导系统载体坐标系的坐标转换矩阵
为:C A B = C n B · C A n = ( C B n ) T · C A n ]]> 则主惯导系统与从惯导相对姿态角的主值为:ψ EAB 0 = arctan C A B ( 3,2 ) 1 - ( C A B ( 3,2 ) ) 2 ]]>ψ NAB 0 = arctan ( - C A B ( 3 , 1 ) C A B ( 3,3 ) ) ]]>ψ UAB 0 = arctan ( - C A B ( 1,2 ) C A B ( 2,2 ) ) ]]> 步骤8:主惯导系统与从惯导系统间的相对姿态置信度估计。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于北京航空航天大学,未经北京航空航天大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201110385043.2/,转载请声明来源钻瓜专利网。