CN113819911A - Gnss失锁下基于自适应容错滤波的导航方法 - Google Patents
Gnss失锁下基于自适应容错滤波的导航方法 Download PDFInfo
- Publication number
- CN113819911A CN113819911A CN202111067791.6A CN202111067791A CN113819911A CN 113819911 A CN113819911 A CN 113819911A CN 202111067791 A CN202111067791 A CN 202111067791A CN 113819911 A CN113819911 A CN 113819911A
- Authority
- CN
- China
- Prior art keywords
- filtering
- matrix
- gnss
- fault
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/26—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network
- G01C21/28—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network with correlation of data from several navigational instruments
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/10—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
- G01C21/12—Navigation; 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/16—Navigation; 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/165—Navigation; 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining 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/42—Determining position
- G01S19/45—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
- G01S19/47—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement the supplementary measurement being an inertial measurement, e.g. tightly coupled inertial
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)
- Signal Processing (AREA)
- Navigation (AREA)
Abstract
本发明公开一种GNSS失锁下基于自适应容错滤波的导航方法,解决SINS/GNSS组合导航在GNSS失锁下使用单纯SINS定位误差越来越大的问题。包括:根据状态转移矩阵与状态变量,建立组合导航的滤波状态方程;根据所述滤波状态方程,卫星有效时结合惯导输出的位置、速度与GNSS接收机输出的位置、速度作差作为滤波器观测值,卫星无效时采用车辆运动学约束构造滤波器观测值,建立滤波观测方程;卫星无效时根据所述滤波观测方程构造用于故障判断的χ2统计量;根据故障判断结果分别计算自适应容错滤波的渐消因子矩阵和加权因子矩阵;根据所述渐消因子矩阵和加权因子矩阵进行自适应容错滤波,结合k时刻的位置、速度信息,得到k时刻的最终导航信息。
Description
技术领域
本发明属于组合导航定位领域,具体涉及一种GNSS失锁下基于自适应容错滤波的导航方法。
背景技术
SINS(Strapdown Inertial Navigation System,捷联惯性导航系统)与GNSS(Global Navigation Satellite System,全球导航卫星系统)良好的互补性,使其成为车载组合导航常用的导航定位方式。GNSS能够在地球表面绝大部分地区准确确定车辆的位置,但在城市环境中,由于树木、建筑物的遮挡,或车辆行驶在隧道、地下停车场等区域时,GNSS由于不能接收足够多的卫星信号会出现失锁情况,仅靠SINS进行导航定位,其定位误差会随时间增加逐渐累积。
目前GNSS失锁引起的导航定位精度降低的问题受到了广泛的关注。现有技术中增加额外传感器辅助SINS能够在一定程度上抑制定位误差的累积。例如,里程计能够提供载体的前向速度,将其与SINS组合能够有效减小导航参数的累积误差。相机的引入可以增加图像的观测信息,能够有效提高系统导航精度。近年来随着AI(人工智能)算法的日渐成熟,利用AI算法辅助SINS可以增强组合导航系统在卫星失效时的性能。但是AI算法计算量较大,在实际使用过程中存在一定的局限性。运动学约束利用车辆在正常行驶过程中横向和竖向速度近似为0的约束条件修正SINS的导航误差,可以在不增加额外传感器的情况下提高导航精度,在低成本组合导航系统中应用广泛。
卡尔曼滤波是多传感器数据融合的经典算法,它要求系统的数学模型和噪声统计特性是先验已知的。但是在实际工程应用中,它们大多是未知的。卡尔曼滤波器使用不准确的系统模型或噪声统计特性会产生较大的状态估计误差,甚至使滤波发散。渐消滤波作为自适应滤波,通过构造渐消因子对一步预测均方误差阵进行修正,从而提高滤波器的调节能力,它通过开窗法确定滤波残差协方差的估计值,对滑动窗口内的数据进行算术平均,不能突出新近观测信息的作用。
车辆在GNSS失锁场景下行驶可能会出现侧滑、跳跃的情况,比如驶过隧道中的积水、减速带、快速转弯,此时车辆不满足运动学的约束条件,故障引入的错误信息会影响滤波效果,甚至发散,影响系统的稳定性。因此,对故障状态进行及时准确的判断,并采取相应的应对措施对于提高系统的导航精度和稳定性十分重要。
发明内容
本发明公开一种GNSS失锁下基于自适应容错滤波的导航方法,解决SINS/GNSS组合导航在GNSS失锁下使用单纯SINS定位误差越来越大的问题。
本发明通过以下技术方案实现。
一种GNSS失锁下基于自适应容错滤波的导航方法,包括:
根据状态转移矩阵与状态变量,建立组合导航的滤波状态方程;
根据所述滤波状态方程,卫星有效时结合惯导输出的位置、速度与GNSS接收机输出的位置、速度作差作为滤波器观测值,卫星无效时采用车辆运动学约束构造滤波器观测值,建立滤波观测方程;
卫星无效时根据所述滤波观测方程构造用于故障判断的χ2统计量;
根据故障判断结果分别计算自适应容错滤波的渐消因子矩阵和加权因子矩阵;
根据所述渐消因子矩阵和加权因子矩阵进行自适应容错滤波,结合k时刻的位置、速度信息,得到k时刻的最终导航信息。
本发明的有益效果:
1、在GNSS失锁场景下,本发明在运动学约束辅助SINS的基础上提出一种自适应容错滤波方法,通过构造χ2统计量来检测系统是否出现故障,设计故障判断准则进行相应处理,在无故障情况下通过限定记忆指数加权的渐消滤波增加对观测信息的利用权重,在发生故障情况下通过残差加权减弱错误的观测信息对滤波更新的影响。
2、传统渐消滤波通过开窗法确定滤波残差协方差的估计值,对滑动窗口内的数据进行算术平均,不能突出新近观测信息的作用,本发明对滑动窗口内的数据进行限定记忆指数加权,并推导给出渐消因子的确定方法,提高了滤波器对新近观测数据的利用权重。
3、本发明设计的故障判断准则对观测故障较敏感,具有告警延迟量小和计算简便的特点。
4、本发明提出的自适应容错滤波算法能有效降低统计特性不准确及异常观测值对系统的不利影响,取得良好的定位效果,具有重要的工程实际意义。
附图说明
图1为本发明GNSS失锁下基于自适应容错滤波的导航方法流程图。
具体实施方式
下面结合附图对本发明作详细说明。
如图1所示,本发明的GNSS失锁下基于自适应容错滤波的导航方法,具体包括:
步骤一、根据状态转移矩阵与状态变量,建立组合导航的滤波状态方程;
Xk=Φk,k-1Xk-1+Γk,k-1Wk-1;
其中,Φk,k-1为状态转移矩阵,Γk,k-1为系统噪声输入矩阵,Xk为k时刻的状态变量,Wk-1为k-1时刻的系统噪声矩阵;
在具体实施时,系统的导航坐标系n系选取为当地地理坐标系,坐标轴分别指向东、北、天方向,车体坐标系m系x、y、z轴分别指向车体的右、前、上方向;惯导系统在载车安装时不能保证惯导坐标系b系和m系重合,存在安装偏差角[αθ αγ αψ]T,所以b系和m系存在转换矩阵不考虑横滚安装偏差角,的误差方程为
式中,δα=[δαθ 0 δαψ]T,“×”表示向量叉乘。δα表示安装误差角向量,δαθ表示俯仰安装误差角,δαψ表示航向安装误差角;
其中,δL,δλ,δh分别为惯导系统的纬度误差、经度误差和高度误差;δvE,δvN,δvU分别为东向、北向、天向速度误差;φE、φN、φU分别为东、北、天三个方向的数学平台误差角;分别为加速度计沿惯导坐标系右、前、上三个方向的随机常值偏置;εx,εy,εz分别为陀螺仪沿惯导坐标系右、前、上三个方向的随机常值漂移;安装误差角δαθ、δαψ为随机常数。
步骤二、根据所述滤波状态方程,当卫星信号有效时结合惯导输出的位置、速度与GNSS接收机输出的位置、速度作差,作为滤波器观测值建立滤波观测方程;
Z1=H1X+V1
其中,观测向量Z1=[δL δλ δh δvE δvN δvU]T,δL,δλ,δh为惯导输出的位置与GNSS接收机输出的位置之差;δvE,δvN,δvU为惯导输出的速度与GNSS接收机输出的速度之差;V1为观测噪声,H1为观测矩阵,表达式为
I3×3表示3维单位矩阵,03×3表示3维零矩阵,03×11表示3×11维零矩阵;
本实施例中,当卫星信号无效时,利用运动学约束条件构造滤波器观测方程:
Z2=H2X+V2
M1(1,*)表示取M1的第1行,其余符号含义类似,V2表示观测噪声;
步骤三、当卫星信号无效时,根据所述滤波观测方程构造用于故障判断的χ2统计量;具体为:
本实施例中,设置信度水平α对应的检验阈值εi,
具体实施时,上述故障判断结果可以参照下表1进行判断:
表1故障判断准则
步骤四、根据故障判断结果分别计算自适应容错滤波的渐消因子矩阵和加权因子矩阵;
本实施例中,所述自适应容错滤波算法具体如下:
其中渐消因子矩阵为Sk=diag(s1,s2,...,sn),加权因子矩阵为ΛK=diag(α1,α2,...,αm),n为系统状态维数,m为观测量维数;
①当系统未发生故障时,Λk=I,所述渐消因子矩阵的计算方法如下:
步骤1:计算k时刻的滤波残差协方差的估计值为
其中,N为窗口大小,b为指数因子,0<b<1。在具体实施时,异常的残差数据一般不记录在窗口内。
步骤2:根据k时刻的滤波残差协方差表达式
Hk=[02×3 D1,2×3 D2,2×3 02×6 D3,2×2]2×17,
Sk=diag(1,1,1,s1,s1,s1,s2,s2,s2,1,1,1,1,1,1,s2,s2),
J1为Jk的对角线第4-6个元素组成的对角阵,J2为Jk的对角线第7-9个元素组成的对角阵,J3为Jk的对角线第16-17个元素组成的对角阵。
通过推导计算
解算得
②当系统观测发生故障时,Sk=I,下面给出加权因子矩阵的计算方法:
其中αi是Λk的第i个元素,Ck,i是Ck的第i个对角元素。
步骤五、根据所述渐消因子矩阵和加权因子矩阵进行自适应容错滤波,结合k时刻的位置、速度信息,得到k时刻的最终导航信息。
综上所述,以上仅为本发明的较佳实例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种GNSS失锁下基于自适应容错滤波的导航方法,其特征在于,包括:
根据状态转移矩阵与状态变量,建立组合导航的滤波状态方程;
根据所述滤波状态方程,卫星有效时结合惯导输出的位置、速度与GNSS接收机输出的位置、速度作差作为滤波器观测值,卫星无效时采用车辆运动学约束构造滤波器观测值,建立滤波观测方程;
当卫星信号无效时,根据所述滤波观测方程构造用于故障判断的χ2统计量;
根据故障判断结果分别计算自适应容错滤波的渐消因子矩阵和加权因子矩阵;
根据所述渐消因子矩阵和加权因子矩阵进行自适应容错滤波,结合k时刻的位置、速度信息,得到k时刻的最终导航信息。
2.如权利要求1所述的一种GNSS失锁下基于自适应容错滤波的导航方法,其特征在于,所述组合导航的滤波状态方程具体为:
Xk=Φk,k-1Xk-1+Γk,k-1Wk-1;
其中,Φk,k-1为状态转移矩阵,Γk,k-1为系统噪声输入矩阵,Xk为k时刻的状态变量,Wk-1为k-1时刻的系统噪声矩阵。
10.如权利要求9所述的一种GNSS失锁下基于自适应容错滤波的导航方法,其特征在于,所述渐消因子矩阵的计算方法具体如下:
①当系统未发生故障时,Λk=I,
步骤1:计算k时刻的滤波残差协方差的估计值为
其中,N为窗口大小,b为指数因子,0<b<1;
步骤2:根据k时刻的滤波残差协方差表达式
Hk=[02×3 D1,2×3 D2,2×3 02×6 D3,2×2]2×17,
Sk=diag(1,1,1,s1,s1,s1,s2,s2,s2,1,1,1,1,1,1,s2,s2)
J1为Jk的对角线第4-6个元素组成的对角阵,J2为Jk的对角线第7-9个元素组成的对角阵,J3为Jk的对角线第16-17个元素组成的对角阵;
通过推导计算
解算得
②当系统观测发生故障时,Sk=I,下面给出加权因子矩阵的计算方法:
其中αi是Λk的第i个元素,Ck,i是Ck的第i个对角元素。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111067791.6A CN113819911A (zh) | 2021-09-13 | 2021-09-13 | Gnss失锁下基于自适应容错滤波的导航方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111067791.6A CN113819911A (zh) | 2021-09-13 | 2021-09-13 | Gnss失锁下基于自适应容错滤波的导航方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113819911A true CN113819911A (zh) | 2021-12-21 |
Family
ID=78914401
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111067791.6A Pending CN113819911A (zh) | 2021-09-13 | 2021-09-13 | Gnss失锁下基于自适应容错滤波的导航方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113819911A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117109571A (zh) * | 2023-10-25 | 2023-11-24 | 北京控制工程研究所 | 导航误差快速收敛方法、装置、电子设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060074558A1 (en) * | 2003-11-26 | 2006-04-06 | Williamson Walton R | Fault-tolerant system, apparatus and method |
CN111189441A (zh) * | 2020-01-10 | 2020-05-22 | 山东大学 | 一种多源自适应容错联邦滤波组合导航系统及导航方法 |
-
2021
- 2021-09-13 CN CN202111067791.6A patent/CN113819911A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060074558A1 (en) * | 2003-11-26 | 2006-04-06 | Williamson Walton R | Fault-tolerant system, apparatus and method |
CN111189441A (zh) * | 2020-01-10 | 2020-05-22 | 山东大学 | 一种多源自适应容错联邦滤波组合导航系统及导航方法 |
Non-Patent Citations (2)
Title |
---|
胡杰等: "基于速度约束与模糊自适应滤波的车载组合导航", 兵工学报, vol. 41, no. 2, pages 231 - 238 * |
麦珍珍: "基于联邦滤波的多传感器组合导航算法研究", 中国优秀硕士学位论文全文库 信息科技辑, no. 02, pages 43 - 55 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117109571A (zh) * | 2023-10-25 | 2023-11-24 | 北京控制工程研究所 | 导航误差快速收敛方法、装置、电子设备及存储介质 |
CN117109571B (zh) * | 2023-10-25 | 2024-01-23 | 北京控制工程研究所 | 导航误差快速收敛方法、装置、电子设备及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111780755B (zh) | 一种基于因子图和可观测度分析的多源融合导航方法 | |
CN110823217B (zh) | 一种基于自适应联邦强跟踪滤波的组合导航容错方法 | |
CN112254718B (zh) | 一种运动约束辅助的基于改进Sage-Husa自适应滤波的水下组合导航方法 | |
CN109724599B (zh) | 一种抗野值的鲁棒卡尔曼滤波sins/dvl组合导航方法 | |
CN109459019B (zh) | 一种基于级联自适应鲁棒联邦滤波的车载导航计算方法 | |
CN111928846B (zh) | 一种基于联邦滤波的多源融合即插即用组合导航方法 | |
CN106679693A (zh) | 一种基于故障检测的矢量信息分配自适应联邦滤波方法 | |
CN112505737B (zh) | 一种gnss/ins组合导航方法 | |
CN110231636B (zh) | Gps与bds双模卫星导航系统的自适应无迹卡尔曼滤波方法 | |
US8898013B2 (en) | Navigation device and process integrating several hybrid inertial navigation systems | |
CN112432644A (zh) | 基于鲁棒自适应无迹卡尔曼滤波的无人艇组合导航方法 | |
CN113984054A (zh) | 基于信息异常检测的改进Sage-Husa自适应融合滤波方法及多源信息融合设备 | |
CN111044051A (zh) | 一种复合翼无人机容错组合导航方法 | |
CN113819911A (zh) | Gnss失锁下基于自适应容错滤波的导航方法 | |
CN111964695A (zh) | 一种无人机组合导航系统故障检测与隔离方法 | |
CN115143954A (zh) | 一种基于多源信息融合的无人车导航方法 | |
CN113008229B (zh) | 一种基于低成本车载传感器的分布式自主组合导航方法 | |
CN116358566B (zh) | 一种基于抗差自适应因子的粗差探测组合导航方法 | |
CN112697154A (zh) | 一种基于矢量分配的自适应多源融合导航方法 | |
CN116224407B (zh) | 一种gnss和ins组合导航定位方法及系统 | |
CN115291253B (zh) | 一种基于残差检测的车辆定位完好性监测方法及系统 | |
CN115014321B (zh) | 一种基于自适应鲁棒滤波的仿生偏振多源融合定向方法 | |
CN115667845A (zh) | 用于移动载体的导航辅助方法 | |
CN111854728B (zh) | 一种基于广义相对熵的容错滤波方法 | |
Eldesoky et al. | Real time localization solution for land vehicle application using low-cost integrated sensors with GPS |
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 |