CN113483759B - 一种GNSS、INS、Vision组合导航系统完好性保护级计算方法 - Google Patents

一种GNSS、INS、Vision组合导航系统完好性保护级计算方法 Download PDF

Info

Publication number
CN113483759B
CN113483759B CN202110723664.0A CN202110723664A CN113483759B CN 113483759 B CN113483759 B CN 113483759B CN 202110723664 A CN202110723664 A CN 202110723664A CN 113483759 B CN113483759 B CN 113483759B
Authority
CN
China
Prior art keywords
gnss
ins
navigation system
vision
integrated 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.)
Active
Application number
CN202110723664.0A
Other languages
English (en)
Other versions
CN113483759A (zh
Inventor
姜海涛
李团
施闯
景贵飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN202110723664.0A priority Critical patent/CN113483759B/zh
Publication of CN113483759A publication Critical patent/CN113483759A/zh
Application granted granted Critical
Publication of CN113483759B publication Critical patent/CN113483759B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/20Instruments for performing navigational calculations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/10Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
    • G01C21/12Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
    • G01C21/16Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
    • G01C21/165Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments
    • G01C21/1656Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments with passive imaging devices, e.g. cameras
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42Determining position
    • G01S19/48Determining position by combining or switching between position solutions derived from the satellite radio beacon positioning system and position solutions derived from a further system
    • G01S19/49Determining position by combining or switching between position solutions derived from the satellite radio beacon positioning system and position solutions derived from a further system whereby the further system is an inertial position system, e.g. loosely-coupled
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Navigation (AREA)

Abstract

本发明适用于GNSS、INS、Vision组合导航数据处理技术领域,提供了一种GNSS、INS、Vision组合导航系统完好性保护级计算方法,包括如下步骤:步骤S100:建立GNSS、INS、Vision组合导航系统的误差状态模型和IMU的零偏模型;步骤S200:建立GNSS和视觉的系统观测模型,并分别计算相应的观测矩阵H;步骤S300:构建不同观测模型下的扩展卡尔曼方程;步骤S400:引入GNSS、INS和视觉故障偏差,并计算不同传感器故障模式下的保护级大小;本发明推导出误差状态模型卡尔曼滤波情况下的组合导航系统保护级公式,进一步考虑到GNSS和视觉观测值的故障偏差在实际情况下是未知的,参照LSR算法经典保护级的取值方法,将其放大为最大特征斜率和最小可检测偏差的乘积,得到不同故障模式下系统的保护级。

Description

一种GNSS、INS、Vision组合导航系统完好性保护级计算方法
技术领域
本发明属于GNSS、INS、Vision组合导航数据处理技术领域,尤其涉及一种GNSS、INS、Vision组合导航系统完好性保护级计算方法。
背景技术
随着自动驾驶、机器人和无人机等技术的兴起,城市等复杂GNSS(GlobalNavigation Satellite System)环境下的高精度定位定姿需求不断增长。GNSS/INS(Inertial Navigation System)组合导航能够提供连续、可靠和完备的高精度导航参数,是应用最广的组合导航技术之一。然而,在GNSS复杂环境下,基于低成本MEMS-IMU(Micro-Electro Mechanical-Systems Inertial Measurement Unit)的GNSS/INS组合导航在GNSS信号中断期间导航误差将迅速积累,这无疑降低了系统的可用性。基于相机的视觉导航具有精度高、自主无源等优点,视觉/INS组合可极大降低INS导航误差。但是,视觉/INS组合无法提供全球框架下的导航信息且误差会随距离累积。因此,将GNSS、INS和视觉进行信息融合可实现优势互补,获得更优的导航性能。
然而多源组合导航系统的优秀性能是建立在系统的复杂性之上,而系统的复杂性势必增加系统故障发生的概率。因此,系统的完好性监测对于组合导航系统的可靠性至关重要。系统的完好性主要包括两个部分:保护级和故障检测。目前针对多源组合导航系统的故障检测主要是针对当前信息的“快照法”和基于历史信息的“连续法”,理论相对成熟,而GNSS、INS、Vision组合导航系统的保护级计算却缺乏具体的理论推导,尤其是基于误差状态模型的扩展卡尔曼滤波下GNSS/INS松组合和Vision/INS紧组合系统。
发明内容
本发明提供一种GNSS、INS、Vision组合导航系统完好性保护级计算方法,旨在解决目前针对多源组合导航系统的故障检测主要是针对当前信息的“快照法”和基于历史信息的“连续法”,理论相对成熟,而多源组合导航系统的保护级计算却缺乏具体的理论推导问题。
本发明是这样实现的,一种GNSS、INS、Vision组合导航系统完好性保护级计算方法,包括如下步骤:
步骤S100:建立GNSS、INS、Vision组合导航系统的误差状态模型和IMU的零偏模型;
步骤S200:建立GNSS和视觉的系统观测模型,并分别计算相应的观测矩阵H;
步骤S300:构建不同观测模型下的扩展卡尔曼方程;
步骤S400:引入GNSS、INS和视觉故障偏差,并计算不同传感器故障模式下的保护级大小。
优选的,所述步骤S100中,将位置、速度以及姿态变量的误差量作为组合导航系统的状态量,把IMU的零偏建模为一阶马尔可夫模型。
优选的,所述步骤S200中,GNSS、INS、Vision组合导航系统可分为两个子系统,分别为GNSS/INS的松组合导航系统和Vision/INS紧组合导航系统。
优选的,所述GNSS/INS松组合导航系统其观测值是GNSS的位置信息,观测模型是INS的位置预测值减去GNSS的定位结果。
优选的,所述Vision/INS紧组合导航系统观测值是特征点在当前相机的位置,观测模型是预测的路标点在当前相机的位置减去当前相机实际观测到的特征点的位置。
优选的,所述步骤S300中,需要构建不同的卡尔曼方程用于更新不同观测值到达时系统的状态量。
与现有技术相比,本发明的有益效果是:本发明的一种GNSS、INS、Vision组合导航系统完好性保护级计算方法,推导出误差状态模型卡尔曼滤波情况下的组合导航系统保护级公式,进一步考虑到GNSS和视觉观测值的故障偏差在实际情况下是未知的,参照LSR算法经典保护级的取值方法,将其放大为最大特征斜率和最小可检测偏差的乘积,便可以得到不同故障模式下系统的保护级。
附图说明
图1为本发明的方法步骤示意图;
图2为本发明中的GNSS、INS、Vision组合导航系统完好性监测的整体框架图;
图中:GNSS:全球导航卫星系统;INS:惯性导航系统;IMU:惯性测量单元;PVA:位置、速度和姿态;MEMS:微机电系统;RCTM:残差卡方检验法;AIME:自主完整性外推法;EKF:扩展卡尔曼滤波;HPL:水平保护级;HAL:水平报警限。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
请参阅图1-2,本发明提供一种技术方案:一种GNSS、INS、Vision组合导航系统完好性保护级计算方法,包括如下步骤:
步骤S100:建立GNSS、INS、Vision组合导航系统的误差状态模型和IMU的零偏模型,把位置,速度,姿态等变量的误差量作为组合导航系统的状态量,考虑到IMU的漂移,把IMU的零偏建模为一阶马尔可夫模型;在确定系统模型的基础上,推导出系统的状态方程,为卡尔曼滤波做准备;
步骤S200:建立GNSS和视觉的系统观测模型,并分别计算相应的观测矩阵H,对于GNSS、INS、Vision组合导航系统,可以分为两个子系统,分别是GNSS/INS的松组合导航系统和Vision/INS紧组合导航系统,对于GNSS/INS松组合导航系统,其观测值是GNSS的位置信息,观测模型是INS的位置预测值减去GNSS的定位结果;对于Vision/INS紧组合导航系统,观测值是特征点在当前相机的位置,观测模型是预测的路标点在当前相机的位置减去当前相机实际观测到的特征点的位置,然后,根据前面的系统状态误差方程和系统的观测模型推算出相应的观测矩阵H;
步骤S300:构建不同观测模型下的扩展卡尔曼方程,因为GNSS/INS的松组合系统和Vision/INS紧组合系统的状态量不同,而且两者的观测模型不同,故相关的参数比如观测矩阵,协方差矩阵等都不一样,为此,需要构建不同的卡尔曼方程用于更新不同观测值到达时系统的状态量;
步骤S400:引入GNSS、INS和视觉故障偏差,并计算不同传感器故障模式下的保护级大小,无论是GNSS/INS松组合系统还是Vision/INS紧组合系统,当观测值到达时,扩展卡尔曼的状态更新方程都不是标准的卡尔曼方程。
本实施方式中,惯性导航机械编排(INS Mechanization)是通过连续的INS导航微分方程推导而得到,地心地固坐标系(e)下的导航微分方程可以表示为:
其中,δre,δve,φ分别为载体坐标系原点(IMU测量中心)相对于e坐标系且投影到e坐标系下的位置,速度和姿态矢量;为地球自转角速度矢量在e系下的投影;δge当地重力加速度矢量在e系下的投影;/>是b系相对于i系的旋转角速度在b系得投影;δfb是b系下的比力矢量;/>是b系到e系的旋转矩阵。
组合导航算法中常将IMU的零偏建模为一阶高斯-马尔科夫过程,模型如下:
其中,bg,ba分别是角速度和加速度的零偏误差;wg,wa是两者对应的高斯白噪声;
则GNSS、INS、Vision组合导航系统的误差状态方程表达如下:
其中,是第i个相机的姿态误差和状态误差;且有:
对于GNSS、INS、Vision组合结构中的GNSS/INS松组合导航系统,观测方程是INS的预测位置减去实际的GNSS解算位置,其表达如下:
系统对应的状态转换矩阵Hk表达如下:
其中,是b坐标系下的杆臂向量。
对于GNSS、INS、Vision组合结构中的Vision/INS紧组合导航系统,当一个特征点fj被相机Ci观测到时,fj的观测残差在Ci的投影表达如下:
其中,是相机的实际观测值;/>是相机的预测观测值,且/>的表达如下:
其中,是相机的观测噪声向量;/>是预测的特征点在相机坐标系下的位置。
进一步可得观测残差近似表达如下:
其中,和/>分别是/>相对于δxk和特征点误差/>求的偏导;/>表达如下:
通常,一个静态特征会被多个连续的相机观测到,因此,可以通过将所有的残差叠加在一起来获得该静态特征的完整残差矢量。其表达如下:
对于EKF,残差线性化需要满足如下形式,即残差与误差项成线性化关系,且噪声项为与状态向量无关的零均值的高斯分布:
当多个特征点被多个相机观测到时,式(14)可以重新表达如下:
ro=ATHXδxk+ATn=Hoδxk+no
(15)
其中,A是半正交矩阵(即ATA=I),A的列组成了Hf的左零空间。
对于GNSS、INS、Vision组合结构中的GNSS/INS松组合导航系统,假设k历元前IMU发生故障,当有GNSS观测值时,得到的EKF优化误差状态是:
其中,(')表示受故障干扰的变量;K'k是受故障影响后的增益矩阵rk表达如下:
rk=Hk·δχIMUk (17)
其中,μk是GNSS观测值的高斯白噪声。
则IMU故障模式下,GNSS/INS松组合导航系统的保护级推导步骤如下:
步骤一、公式(16)两边同时减去误差状态量δχIMU,得到(18):
步骤二、公式(18)两边同时乘以Hk,得到(19):
步骤三、把公式(17)带入公式(19),可得公式(20)
步骤四、公式(20)两边同时乘以Lk
其中,Lk的表达如下:
则公式(21)可以重新写成如下的形式:
则IMU故障模式的保护级为:
其中,下标’q’表示的第q个分量,σq表示Lkμk的协方差矩阵第q行对角线元素的平方根,Kmd,imu与GNSS/INS松组合系统IMU故障情况下的完整性风险要求有关。
对于GNSS、INS、Vision组合结构中的GNSS/INS松组合导航系统,假设k历元前GNSS发生故障,则rk可以重新表达如下:
rk=Hk·δxIMUk+fk,GNSS
(25)
其中,fk,GNSS是GNSS观测值的偏差向量。
可以得到:
则GNSS故障模式的保护级为:
PLq,gnss=[(K'k-Lk)rk]q+Kmd,gnssq+[Lkfk,GNSS]q (27)
其中,Kmd,gnss与GNSS/INS松组合系统GNSS故障情况下的完整性风险要求有关。
考虑到式(27)中包含fk,GNSS且未知,此处参考LSR算法经典保护级的取值方法,将Lkfk,GNSS放大为最大特征斜率Max(Slopemq)与最小可检测偏差的乘积。其中特征斜率Slopemq根据fk,GNSS对GNSS故障检测检验统计量Ts和滤波估计的影响进行定义;对于给定的漏检概率,可以从测试统计信息的阈值计算/>
fk,GNSS在滤波估计第q维引起的误差为:
μq=(Lk)qmξb=lqmξb (28)
其中,lqm代表Lk第q行第m列的元素;且有fk,GNSS=[ξb ξb ξb]T.
而GNSS故障检测检验统计量Ts表达如下:
Ts=rk T·Sk -1·rk (29)
fk,GNSS对检测统计量Ts的影响表现在非中心化参数λ上面,λ的表达如下:
λ=ξb 2S_mm (30)
其中,S_mm是矩阵Sk -1=(Hk·Pk,k-1·Hk T+Rk)-1的第m行对角线元素,则Slopemq的定义如下:
则GNSS故障模式的保护级重新表达为:
对于GNSS、INS、Vision组合结构中的Vision/INS紧组合导航系统,假设k历元前相机发生故障偏差fx,visual,当有视觉观测值时,得到的EKF优化误差状态是:
则ro可以重新表达如下:
ro=Hoδxk+no+fo,visual (33)
参考GNSS故障,可以得到:
其中,fo,visual=ATfχ,visual,且Lk,o的表达如下:
考虑到求逆时可能会出现问题,且根据式(15)可得:
因此,Lk,o可以重新表达如下:
式(34)重新写成:
其中
则视觉故障模式下,Vision/INS紧组合系统的保护级表达为:
σq,o表示Bk,on的协方差矩阵第q行对角线元素的平方根;Kmd,cam与Vision/INS紧组合系统相机故障情况下的完整性风险要求有关;对于给定的漏检概率,可以从测试统计信息的阈值计算Slopemq,o的计算如下:
其中,bqm,o表示矩阵Bk,o的q行,m列对应的元素;是第m个路标点观测噪声的标准差。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (2)

1.一种GNSS、INS、Vision组合导航系统完好性保护级计算方法,其特征在于:包括如下步骤:
步骤S100:建立GNSS、INS、Vision组合导航系统的误差状态模型和IMU的零偏模型;
步骤S200:建立GNSS和视觉的系统观测模型,并分别计算相应的观测矩阵H,GNSS、INS、Vision组合导航系统分为两个子系统,分别为GNSS/INS的松组合导航系统和Vision/INS紧组合导航系统,所述GNSS/INS松组合导航系统其观测值是GNSS的位置信息,观测模型是INS的位置预测值减去GNSS的定位结果,所述Vision/INS紧组合导航系统观测值是特征点在当前相机的位置,观测模型是预测的路标点在当前相机的位置减去当前相机实际观测到的特征点的位置;
步骤S300:构建不同观测模型下的扩展卡尔曼方程,需要构建不同的卡尔曼方程用于更新不同观测值到达时系统的状态量;
步骤S400:引入GNSS、INS和视觉故障偏差,并计算不同传感器故障模式下的保护级大小。
2.如权利要求1所述的一种GNSS、INS、Vision组合导航系统完好性保护级计算方法,其特征在于:所述步骤S100中,将位置、速度以及姿态变量的误差量作为组合导航系统的状态量,把IMU的零偏建模为一阶马尔可夫模型。
CN202110723664.0A 2021-06-29 2021-06-29 一种GNSS、INS、Vision组合导航系统完好性保护级计算方法 Active CN113483759B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110723664.0A CN113483759B (zh) 2021-06-29 2021-06-29 一种GNSS、INS、Vision组合导航系统完好性保护级计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110723664.0A CN113483759B (zh) 2021-06-29 2021-06-29 一种GNSS、INS、Vision组合导航系统完好性保护级计算方法

Publications (2)

Publication Number Publication Date
CN113483759A CN113483759A (zh) 2021-10-08
CN113483759B true CN113483759B (zh) 2023-10-17

Family

ID=77936568

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110723664.0A Active CN113483759B (zh) 2021-06-29 2021-06-29 一种GNSS、INS、Vision组合导航系统完好性保护级计算方法

Country Status (1)

Country Link
CN (1) CN113483759B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114136315B (zh) * 2021-11-30 2024-04-16 山东天星北斗信息科技有限公司 一种基于单目视觉辅助惯性组合导航方法及系统
CN115408652B (zh) * 2022-10-31 2023-02-03 北京航空航天大学 一种基于截断误差估计的ekf组合导航保护级确定方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2784445A2 (en) * 2013-03-26 2014-10-01 Honeywell International Inc. Selected aspects of advanced receiver autonomous integrity monitoring application to kalman filter based navigation filter
CN108088470A (zh) * 2018-01-30 2018-05-29 北京大学 一种组合导航慢变斜坡故障完好性监测方法
CN108692739A (zh) * 2017-03-29 2018-10-23 霍尼韦尔国际公司 用于具有异构测量结果的导航系统的完整性监视方法
CN109900300A (zh) * 2019-03-27 2019-06-18 北京航空航天大学 一种用于无人机的组合导航完好性监测系统
CN111060133A (zh) * 2019-12-04 2020-04-24 南京航空航天大学 一种用于城市复杂环境的组合导航完好性监测方法
CN111123304A (zh) * 2019-11-28 2020-05-08 北京航空航天大学 一种视觉导航完好性监测计算方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR3028310A1 (fr) * 2014-11-07 2016-05-13 Thales Sa Procede de dertemination d'un rayon de protection associe a un parametre de navigation d'un systeme hybride de navigation inertielle, et systeme associe
EP3792665A1 (en) * 2019-09-10 2021-03-17 Trimble Inc. Protection level generation methods and systems for applications using navigation satellite system (nss) observations

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2784445A2 (en) * 2013-03-26 2014-10-01 Honeywell International Inc. Selected aspects of advanced receiver autonomous integrity monitoring application to kalman filter based navigation filter
CN108692739A (zh) * 2017-03-29 2018-10-23 霍尼韦尔国际公司 用于具有异构测量结果的导航系统的完整性监视方法
CN108088470A (zh) * 2018-01-30 2018-05-29 北京大学 一种组合导航慢变斜坡故障完好性监测方法
CN109900300A (zh) * 2019-03-27 2019-06-18 北京航空航天大学 一种用于无人机的组合导航完好性监测系统
US10466362B1 (en) * 2019-03-27 2019-11-05 Beihang University Integrated navigation integrity monitoring system for unmanned aerial vehicles
CN111123304A (zh) * 2019-11-28 2020-05-08 北京航空航天大学 一种视觉导航完好性监测计算方法
CN111060133A (zh) * 2019-12-04 2020-04-24 南京航空航天大学 一种用于城市复杂环境的组合导航完好性监测方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
A High-Integrity and Low-Cost Navigation System for Autonomous Vehicles;Bijjahalli Suraj et al.;《IEEE TRANSACTIONS ON INTELLIGENT TRANSPORTATION SYSTEMS》;第22卷(第1期);356-369 *
BDS/GPS组合导航接收机自主完好性监测算法;王尔申;杨福霞;庞涛;曲萍萍;蓝晓宇;姜毅;;北京航空航天大学学报(第04期);35-41 *
Vision-integrated Navigation and Integrity Monitoring for Aircraft Final Approach;Watanabe Yoko et al.;《IFAC PAPERSONLINE》;第53卷(第2期);14730-14735 *
多传感器组合导航完好性研究;史进恒等;《第十一届中国卫星导航年会论文集》;1-6 *
用于精细化无人机管控的GNSS/INS组合导航定位及完好性监测算法;张文宇等;《第十一届中国卫星导航年会论文集》;1-8 *

Also Published As

Publication number Publication date
CN113483759A (zh) 2021-10-08

Similar Documents

Publication Publication Date Title
Youn et al. Combined quaternion-based error state Kalman filtering and smooth variable structure filtering for robust attitude estimation
Ahmed et al. Accurate attitude estimation of a moving land vehicle using low-cost MEMS IMU sensors
US7979231B2 (en) Method and system for estimation of inertial sensor errors in remote inertial measurement unit
CN111102978B (zh) 一种车辆运动状态确定的方法、装置及电子设备
JP5036462B2 (ja) ナビゲーション・システムおよびナビゲーション方法
Nebot et al. Initial calibration and alignment of low‐cost inertial navigation units for land vehicle applications
EP3043148B1 (en) Heading for a hybrid navigation solution based on magnetically calibrated measurements
US8185309B2 (en) Enhanced inertial system performance
US6205401B1 (en) Navigation system for a vehicle, especially a land craft
CN113483759B (zh) 一种GNSS、INS、Vision组合导航系统完好性保护级计算方法
JP7036080B2 (ja) 慣性航法装置
CN108106635A (zh) 惯性卫导组合导航系统的长航时抗干扰姿态航向校准方法
CN112629538A (zh) 基于融合互补滤波和卡尔曼滤波的舰船水平姿态测量方法
JP2012194175A (ja) 姿勢判定方法、位置算出方法及び姿勢判定装置
KR20060048581A (ko) 네비게이션 중복성을 위한 시스템
CN114002725A (zh) 一种车道线辅助定位方法、装置、电子设备及存储介质
CN108981708B (zh) 四旋翼扭矩模型/航向陀螺/磁传感器容错组合导航方法
CN110954102A (zh) 用于机器人定位的磁力计辅助惯性导航系统及方法
US10048074B1 (en) Polar region operating attitude and heading reference system
CN112683269A (zh) 一种附有运动加速度补偿的marg姿态计算方法
CN117053782A (zh) 一种水陆两栖机器人组合导航方法
Youn Magnetic fault–tolerant navigation filter for a UAV
JP2007232444A (ja) 慣性航法装置およびその誤差補正方法
JP7028223B2 (ja) 自己位置推定装置
CN117537814A (zh) 一种矩阵李群上的无迹卡尔曼滤波初始对准方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant