CN101038170A - 一种在线估计惯性/卫星组合导航***数据不同步时间的方法 - Google Patents
一种在线估计惯性/卫星组合导航***数据不同步时间的方法 Download PDFInfo
- Publication number
- CN101038170A CN101038170A CNA2007100636394A CN200710063639A CN101038170A CN 101038170 A CN101038170 A CN 101038170A CN A2007100636394 A CNA2007100636394 A CN A2007100636394A CN 200710063639 A CN200710063639 A CN 200710063639A CN 101038170 A CN101038170 A CN 101038170A
- Authority
- CN
- China
- Prior art keywords
- navigation system
- asynchronous time
- error
- satellite
- inertial navigation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Navigation (AREA)
Abstract
一种在线估计惯性/卫星组合导航***数据不同步时间的方法,将惯导***与卫星导航***之间的数据不同步时间作为估计器的状态,因此状态方程包含两部分,一部分是惯导***的误差状态,另一部分是两个子***不同步时间的误差模型;量测量采用两个***输出位置之差,由于两个***的数据不是同一个时刻的更新数据,因而量测模型中除了两个子***相应的测量误差外,还包含了由于数据不同步时间造成的位置差;根据状态方程和量测方程用滤波的方法计算得到状态量的最优估计,包括两个***数据不同步时间。本发明能够实时估计出惯导***和卫星导航***之间的数据不同步时间,有效提高组合导航***的精度,并且简单实用,不需要增加任何硬件负担。
Description
技术领域
本发明涉及一种在线估计惯性/卫星组合导航***数据不同步时间的方法,属于导航定位技术领域。
背景技术
惯性导航***(INS)自主性和隐蔽性好,能连续提供多种较高精度的导航参数(位置、速度、姿态、航向等)的输出信息,频带宽,但其误差(尤其是位置误差)随时间积累,不能长时间单独承担高精度导航的任务。而卫星导航***的定位和测速精度高,且基本不受地域、时间限制。但当载体作大机动飞行或有地形遮蔽时,卫星导航信息有可能中断,或动态误差过大,不能使用;另外卫星接收机数据的更新频率比较低,难以满足实时控制的要求。由于惯导***和卫星导航***具有互补的特点,惯性/卫星组合导航***成为目前最理想的导航***,在各种领域得到了广泛的应用。
组合导航***就是以最优的方式融合来自各子***的数据,综合它们各自的优势,并提供比任何单一子***更精确、更可靠的导航输出。在进行数据融合设计时,首先必须保证用于融合的导航数据在时间上是一致的,组合导航的设计才具有实际意义。实际上,组合导航***中的各个子***往往具有不同的数据更新率,时标漂移以及计算和通信的时延也是客观存在的,这些都会导致子***的数据不同步。不同步误差对导航***性能的影响是非常显著的,特别是在动态应用环境中。所以在组合导航***的设计中,研究有效的数据时间同步技术具有非常重要的意义。
目前解决数据时间同步问题的方法主要是利用硬件或软件定时器来获取惯导***和卫星导航***的不同步时间,再利用高阶保持器得到同步点上的外推数据值,从而实现数据同步。这种方法不仅增加了***的成本和复杂度,而且在同步过程中需要占用较多的CPU时间。
发明内容
本发明的技术问题是:针对现有技术的不足,提供一种简单实用的在线估计组合导航***中惯导***和卫星导航***之间数据不同步时间的方法。
本发明的技术解决方案:一种在线估计惯性/卫星组合导航***数据不同步时间的方法,其特点在于:
(1)分别建立相应的惯导***误差模型,及惯导***和卫星导航***之间不同步时间的误差模型;
(2)选取所述的惯导***和所述的惯导***和卫星导航***之间不同步时间误差作为估计器的状态变量;
(3)采用所述的惯导***给出的位置信息和所述的卫星***给出的位置信息的差值作为量测量,建立相应的量测方程;
(4)采用滤波方法估计所述的惯导误差和所述的惯导***和卫星导航***之间不同步时间;
(5)最后对所述的惯导***的输出进行校正,以提高组合导航***的精度。
本发明的原理:在建立了惯导和卫星***的数据不同步时间误差模型的基础上,采用估计器对其进行在线估计,估计器的状态量包括了惯导***的误差及两个***数据不同步时间误差;量测量为惯导***给出的位置和由卫星导航***得到的位置之差;根据状态方程和量测方程用滤波的方法计算得到对输入状态量的最优估计。
本发明与现有技术相比具有如下优点:本发明利用估计器在组合过程中在线估计惯导***和卫星导航***之间数据的不同步时间,能够大大提高惯性/卫星组合导航***的精度,简单实用,而且不需要增加任何硬件负担。
附图说明
图1为实施本发明的在线估计***的模型;
图2为实施本发明的飞行轨迹(经-纬-高度曲线)
图3为实施本发明后不同步时间的估计曲线(从起飞开始估计);
图4为实施本发明后位置估计误差曲线(从起飞开始估计);
图5为实施本发明后不同步时间的估计曲线(从拐弯开始估计);
图6为实施本发明后位置估计误差曲线(从拐弯开始估计)。
具体实施方式
下面以一架飞机的飞行过程为实例来阐述本发明的具体实施过程。
图1是实施本发明的在线估计***的模型,该***的核心部分是构建估计器的状态方程和量测方程。本发明将惯导***与卫星导航***之间的数据不同步时间作为估计器的状态,在估计惯导***位置、速度误差等状态的同时,也将它估计出来。因此状态方程由两部分组成,一部分是惯导***的误差状态,另一部分是两个子***不同步时间的误差模型。
惯导***的主要误差包括速度误差、姿态误差、位置误差以及惯性器件的误差。选择东北天地理坐标系作为建立动力学方程的基准坐标系,设δVE,δVN,δVU分别代表沿东、北、天方向的速度误差;φE,φN,φU为平台误差角;δL,δλ,δh分别代表纬度、经度及高度误差,εx,εy,εz为陀螺漂移,x,y,z为加速度计零位误差,并假设沿机体坐标系三个轴向安装的陀螺和加速度计的误差模型相同,均为一阶马尔可夫过程,则惯导的误差方程表示为:
(1)速度误差方程
(2)平台角误差方程
(3)位置误差方程
(4)陀螺漂移误差方程
(5)加速度计漂移误差方程
写成矩阵形式如下:
其中XI(t)是***状态向量,WI(t)是***噪声向量,FI(t)是***矩阵,GI(t)是***噪声矩阵。
两个子***不同步时间误差模型的建立步骤或方法如下:
(1)对不同步时间建立数学模型:
惯导***和卫星接收机的数据更新率是不同的,一般来说,惯导数据更新率是50~100Hz,卫星的数据更新频率是1Hz。尽管从理论上说,这些子***的信息更新频率是固定的,但是在实际的导航过程中,由于温度特性等因素影响,惯导***的频标会出现漂移。GPS卫星***准确地在秒同步脉冲(1PPS,每秒一个脉冲)到来时进行伪矩的测量,但是必须注意,导航信息更新时刻并不是数据输出时刻。在每个导航子***中,导航信息更新完成后,还要经过计算和通信才能最终输出导航数据。一般来说,由于频标漂移而导致的数据更新周期的变化是缓慢而细微的,而卫星输出信息相对惯导的滞后时间是主要的不同步时间。基于以上的分析,惯导***与卫星导航***之间的不同步时间包含两部分,一是两个***时间基准之间的差值,二是时钟相对钟漂的大小,前者可表示为一个随机常数,而后者可用一阶马尔可夫过程来描述:
τa=τc+τr+wa
式中τc为随机常数;τr为一阶马尔可夫过程;wa为白噪声。
(2)对所建立的不同步时间的数学模型求导,即可得到相应的误差模型:
式中βτ为相关时间。
不同步时间误差模型的矩阵形式为:
估计器的量测方程的建立方法如下:
(1)采用两个***输出的位置之差作为量测量;
(2)从运动学角度考虑两个***间不同步时间对位置差的影响,并将其作为量测模型的一部分。当载体的速度发生变化时(方向或大小),不同步时间的可观测性越好;
(3)综合两个***各自的测量误差,建立完整的测量方程。
根据建立的状态方程和量测方程采用滤波的方法计算得到各个状态量的最优估计,包括两个***数据不同步时间。
在得到位置、速度误差的估计值后,对惯导***的相应输出进行校正,并以校正后的信息作为组合导航***的输出,从而提高了组合导航***的精度。
从以上分析可得,只要飞机有机动飞行过程就能够估计出两个***间的数据不同步时间。由于飞机在正常飞行时,都有起飞的飞行过程,而飞机的起飞过程中一般都包含有加速滑跑、爬升、拉平、拐弯飞行等必需的飞行动作,因此,可以利用飞机的飞行过程完成两个子***之间的数据时间同步。考虑到惯导***一般在地面就能够通电工作,而卫星***可能在飞机起飞时就能够工作,也可能由于遮挡等原因,需要飞机飞到一定高度方能正常工作,因而我们也分别考虑了两种情况。第一种情况是卫星***在飞机起飞时就能够正常工作,因而可以利用整个飞行过程,包括加速滑跑、爬升、拉平、拐弯等飞行动作,来估计两个***数据之间的不同步时间。
另一种情况是当飞机飞到一定高度后,卫星***才能够正常工作,这时飞机只需做一个加速飞行或拐弯的飞行动作,一样可以估计出不同步时间。由于加速飞行的过程与起飞时的加速滑跑类似,因此仿真了飞机作拐弯动作时本发明方法对不同步时间、组合***位置的估计情况。
在仿真中,均假设在机体坐标系三轴使用性能相同的陀螺仪和加速度计,误差特性一致,均设为有色噪声,其中陀螺随机常值漂移为0.1°/hr,一阶马尔可夫漂移为0.1°/hr,相关时间为3600s,白噪声为0.01°/hr;加速度计零偏为10-3g,一阶马尔可夫过程为10-3g,相关时间为3600s,白噪声标准差为10-3g。而卫星***的位置测量误差为20m。
图2仿真了一个飞机飞行过程的轨迹曲线,飞机从静止开始作加速滑跑,拉高起飞,爬升,到一定的高度后拉平,之后拐弯、平飞至规定的航线。对本发明方法的验证仿真就采用了图2所示的飞行轨迹。
图3是在飞机起飞时本发明方法对不同步时间的估计结果,图4是在飞机起飞时本发明方法对组合***位置的估计误差,以及相同仿真条件下当不估计不同步时间时,组合***的位置估计误差。从图3和图4可以看到,实施本发明10秒左右,即可估计出不同步时间并且估计误差为0.05秒(1σ),而位置估计误差由不估计不同步时间时的150米左右降到了15米左右,精度提高近一个数量级。
图5是在飞机作拐弯动作时本发明方法对不同步时间的估计结果,图6是在飞机作拐弯动作时本发明方法对组合***位置的估计误差,以及相同仿真条件下当不估计不同步时间时,组合***的位置估计误差。从图5和图6可以看到,实施本发明6~7秒左右,即可估计出不同步时间并且估计误差为0.05秒(1σ),而位置估计误差由不估计不同步时间时的100多米降到了10米左右,精度同样提高近一个数量级。
从两种情况的仿真结果来看,利用本发明提出的方法,均能够很快估计出两个***之间的数据不同步时间(大约在10秒左右),且估计误差为0.05秒(1σ)。另外,在估计不同步时间的同时,数据融合的精度与不估计不同步时间的情况相比,精度提高近一个数量级,证明这种方法是非常有效的。
最后所应说明的是:以上实施实例仅用以说明而非限制本发明的技术方案,所有的不脱离本发明的精神和范围的修改或局部替换,均应涵盖在本发明的权利要求范围当中。
Claims (3)
1、一种在线估计惯性/卫星组合导航***数据不同步时间的方法,其特征在于:
(1)分别建立相应的惯导***误差模型,及惯导***和卫星导航***之间不同步时间的误差模型;
(2)选取所述的惯导***和所述的惯导***和卫星导航***之间不同步时间误差作为估计器的状态变量;
(3)采用所述的惯导***给出的位置信息和所述的卫星***给出的位置信息的差值作为量测量,建立相应的量测方程;
(4)采用滤波方法估计所述的惯导误差和所述的惯导***和卫星导航***之间不同步时间;
(5)最后对所述的惯导***的输出进行校正,以提高组合导航***的精度。
2、根据权利要求1所述的一种在线估计惯性/卫星组合导航***数据不同步时间的方法,其特征还在于:所述的惯导***和卫星导航***之间不同步时间误差模型的建立方法如下:
(1)对不同步时间建立数学模型:不同步时间包含两部分,一是两个***时间基准之间的差值,二是时钟相对钟漂的大小,前者可表示为一个随机常数,而后者可用一阶马尔可夫过程来描述;
(2)对所建立的不同步时间的数学模型求导,即可得到相应的误差模型。
3、根据权利要求1所述的一种在线估计惯性/卫星组合导航***数据不同步时间的方法,其特征还在于:所述的量测方程的建立方法如下:
(1)采用两个***输出的位置之差作为量测量;
(2)从运动学角度考虑两个***间不同步时间对位置差的影响,并将其作为量测模型的一部分;
(3)综合两个***各自的测量误差,建立完整的测量方程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007100636394A CN101038170B (zh) | 2007-02-07 | 2007-02-07 | 一种在线估计惯性/卫星组合导航***数据不同步时间的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007100636394A CN101038170B (zh) | 2007-02-07 | 2007-02-07 | 一种在线估计惯性/卫星组合导航***数据不同步时间的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101038170A true CN101038170A (zh) | 2007-09-19 |
CN101038170B CN101038170B (zh) | 2010-06-30 |
Family
ID=38889230
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2007100636394A Expired - Fee Related CN101038170B (zh) | 2007-02-07 | 2007-02-07 | 一种在线估计惯性/卫星组合导航***数据不同步时间的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101038170B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101793529A (zh) * | 2010-03-03 | 2010-08-04 | 北京航空航天大学 | 一种双伪卫星辅助惯性导航***位置校准方法 |
CN102589570A (zh) * | 2012-01-17 | 2012-07-18 | 北京理工大学 | 一种船用惯性导航***单点海上校准方法 |
CN106767788A (zh) * | 2017-01-04 | 2017-05-31 | 北京航天自动控制研究所 | 一种组合导航方法和*** |
CN109724598A (zh) * | 2019-03-08 | 2019-05-07 | 哈尔滨工程大学 | 一种gnss/ins松组合时延误差的估计及补偿方法 |
CN111256691A (zh) * | 2020-02-17 | 2020-06-09 | 苏州芯智谷智能科技有限公司 | 基于gnss/mems惯性组合芯片的组网硬件时间基准建立方法 |
CN113311463A (zh) * | 2020-02-26 | 2021-08-27 | 北京三快在线科技有限公司 | Gps延迟时间在线补偿方法、装置、电子设备和存储介质 |
CN115597571A (zh) * | 2022-12-15 | 2023-01-13 | 西南应用磁学研究所(中国电子科技集团公司第九研究所)(Cn) | 电子罗盘传感器误差及安装误差快速标定与补偿方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1372127A (zh) * | 2001-01-23 | 2002-10-02 | 林清芳 | 改进的定位和数据集成方法及其*** |
CN1139786C (zh) * | 2002-02-06 | 2004-02-25 | 何秀凤 | 自主式定位定向导航仪 |
CN1851407A (zh) * | 2006-05-30 | 2006-10-25 | 梁保山 | 高精度微机电组合惯导装置 |
-
2007
- 2007-02-07 CN CN2007100636394A patent/CN101038170B/zh not_active Expired - Fee Related
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101793529A (zh) * | 2010-03-03 | 2010-08-04 | 北京航空航天大学 | 一种双伪卫星辅助惯性导航***位置校准方法 |
CN101793529B (zh) * | 2010-03-03 | 2012-05-23 | 北京航空航天大学 | 一种双伪卫星辅助惯性导航***位置校准方法 |
CN102589570A (zh) * | 2012-01-17 | 2012-07-18 | 北京理工大学 | 一种船用惯性导航***单点海上校准方法 |
CN106767788A (zh) * | 2017-01-04 | 2017-05-31 | 北京航天自动控制研究所 | 一种组合导航方法和*** |
CN106767788B (zh) * | 2017-01-04 | 2019-07-19 | 北京航天自动控制研究所 | 一种组合导航方法和*** |
CN109724598A (zh) * | 2019-03-08 | 2019-05-07 | 哈尔滨工程大学 | 一种gnss/ins松组合时延误差的估计及补偿方法 |
CN111256691A (zh) * | 2020-02-17 | 2020-06-09 | 苏州芯智谷智能科技有限公司 | 基于gnss/mems惯性组合芯片的组网硬件时间基准建立方法 |
CN113311463A (zh) * | 2020-02-26 | 2021-08-27 | 北京三快在线科技有限公司 | Gps延迟时间在线补偿方法、装置、电子设备和存储介质 |
CN115597571A (zh) * | 2022-12-15 | 2023-01-13 | 西南应用磁学研究所(中国电子科技集团公司第九研究所)(Cn) | 电子罗盘传感器误差及安装误差快速标定与补偿方法 |
CN115597571B (zh) * | 2022-12-15 | 2023-03-28 | 西南应用磁学研究所(中国电子科技集团公司第九研究所) | 电子罗盘传感器误差及安装误差快速标定与补偿方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101038170B (zh) | 2010-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108535755B (zh) | 基于mems的gnss/imu车载实时组合导航方法 | |
CN108226980B (zh) | 基于惯性测量单元的差分gnss与ins自适应紧耦合导航方法 | |
CN101858748B (zh) | 高空长航无人机的多传感器容错自主导航方法 | |
CN107588769B (zh) | 一种车载捷联惯导、里程计及高程计组合导航方法 | |
CN101256080B (zh) | 卫星/惯性组合导航***的空中对准方法 | |
CN110487301A (zh) | 一种雷达辅助机载捷联惯性导航***初始对准方法 | |
CN104655131B (zh) | 基于istssrckf的惯性导航初始对准方法 | |
CN101038170A (zh) | 一种在线估计惯性/卫星组合导航***数据不同步时间的方法 | |
CN109931955B (zh) | 基于状态相关李群滤波的捷联惯性导航***初始对准方法 | |
CN103363992B (zh) | 基于梯度下降的四旋翼无人机姿态航向参考***解算方法 | |
CN103557871A (zh) | 一种浮空飞行器捷联惯导空中初始对准方法 | |
CN104833352A (zh) | 多介质复杂环境下高精度视觉/惯性组合导航方法 | |
CN108759824A (zh) | 高精度定位导航***及方法 | |
CN109708663B (zh) | 基于空天飞机sins辅助的星敏感器在线标定方法 | |
CN104457748A (zh) | 一种嵌入式瞄准吊舱测姿***及其传递对准方法 | |
CN102645223B (zh) | 一种基于比力观测的捷联惯导真空滤波修正方法 | |
CN104165641A (zh) | 一种基于捷联惯导/激光测速仪组合导航***的里程计标定方法 | |
CN103674034A (zh) | 多波束测速测距修正的鲁棒导航方法 | |
CN103630136A (zh) | 冗余传感器配置下基于三级滤波的导航参数最优融合方法 | |
CN103900566B (zh) | 一种消除地球自转角速度对旋转调制型捷联惯导***精度影响的方法 | |
CN105241456A (zh) | 巡飞弹高精度组合导航方法 | |
CN103674064A (zh) | 捷联惯性导航***的初始标定方法 | |
Liu et al. | Interacting multiple model UAV navigation algorithm based on a robust cubature Kalman filter | |
CN101929862A (zh) | 基于卡尔曼滤波的惯性导航***初始姿态确定方法 | |
CN102707080A (zh) | 一种星敏感器模拟捷联惯导陀螺的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100630 Termination date: 20170207 |