一种抑制GNSS信息异常的滤波增益动态调整方法技术

本发明专利技术公开了一种抑制GNSS信息异常的滤波增益动态调整方法，在标准卡尔曼滤波算法的基础上，通过新息向量构造符合χ2分布的指标量并进行χ2检验，对观测异常值进行检测。根据χ2检验的参数构造比例因子，使用比例因子对相应的滤波增益矩阵进行降低，达到抑制异常观测值对系统影响的目的。本发明专利技术对正常噪声无影响，对异常值无错检、漏检，且对单独时刻异常值和持续异常值均有较高效作用，可使系统性能提高90％以上，将异常值的作用时间降5倍左右，大大提高了组合导航系统的性能。

【技术实现步骤摘要】

本专利技术涉及INS/GNSS组合导航系统，特别是涉及一种抑制GNSS信息异常的滤波增益动态调整方法。
技术介绍
惯性导航系统(Inertial navigation system，INS)/全球导航卫星系统(Global navigation satellite system，GNSS)组合导航系统中INS数据信息更新率高、噪声低，在短时间内有较高精度，但其系统原理决定了导航误差会随时间增加而积累，GNSS误差与时间无关，可全天候、全球范围提供高精度的导航信息，但其信息更新率较低。INS/GNSS系统经过滤波方法进行信息融合，可实现优势互补，得到更加精确的导航结果。卡尔曼滤波器(Kalman filter，KF)对于过程噪声和观测噪声为Gauss白噪声序列的线性系统，其滤波结果在无偏、一致和渐进有效的意义下是最优的。GNSS信号易受外界环境干扰的影响，尤其是在高楼林立的城市或是隧道峡谷等对信号遮挡严重的区域或者多径效应明显的环境中，GNSS会出现较大的观测异常值。此时观测值密度会表现为拖尾分布，观测噪声不再服从Gauss分布，不符合卡尔曼滤波的基本假设，滤波性能下降，无法满足需求。在卡尔曼滤波的基础上，相继出现了粒子滤波、自适应滤波和H∞滤波等抗差性较好的滤波方法。粒子滤波可适用于非线性系统，抗差性能较好，但其计算量大、实时性较差，在组合导航领域应用较少；自适应滤波方法包括LMS和RLS方法，其跟踪性能较差；H∞滤波依据滤波器性能指标分析误差信息，能量最小，其牺牲平均估计精度来换取抗差性能。
技术实现思路
专利技术目的：本专利技术的目的是提供一种能够抑制INS/GNSS组合导航系统中GNSS信息异常的滤波增益动态调整方法。技术方案：为达到此目的，本专利技术采用以下技术方案：本专利技术所述的抑制GNSS信息异常的滤波增益动态调整方法，包括以下步骤：S1：对标准卡尔曼滤波算法中的新息向量ik构建符合χ2分布的指标量I进行χ2检验，原假设H0:I～χ2(t)，指标量I如式(1)所示： I = i k T ( H k P k , k - 1 H k T + R k ) - 1 i k - - - ( 1 ) ]]>式(1)中，ik为新息向量，Hk为观测矩阵，Pk,k-1为由k-1时刻到k时刻的预测误差方差矩阵，Rk为观测噪声协方差矩阵；S2：选取α为显著性水平；S3：设通过查表方法得到α对应的上分位点为阈值；S4：对各个时刻的指标量I进行检测：如果连续n个时刻的n≥预设时间，则放弃GNSS，仅采用INS导航，进行步骤S7；否则，则进行步骤S5；S5：构建比例因子μ： μ = I k χ α 2 - - - ( 2 ) ]]>式(2)中，Ik为新息向量矩阵；S6：用代替Kk重新进行滤波计算，得到对异常值抑制后的滤波结果，如式(3)所示： K ‾ k = K k / μ - - - ( 3 ) ]]>式(3)中，Kk是滤波器增益矩阵；S7：结束。有益效果：本专利技术在标准卡尔曼滤波算法的基础上提出了一种检测和抑制异常观测值影响的方法，对正常噪声无影响，对异常值无错检、漏检，且对单独时刻异常值和持续异常值均有较高效作用，可使系统性能提高90％以上，将异常值的作用时间降5倍左右，大大提高了组合导航系统的性能。具体实施方式下面结合具体实施方式对本专利技术的技术方案作进一步的介绍。本专利技术在标准卡尔曼滤波算法的基础上，公开了一种抑制GNSS信息异常的滤波增益动态调整方法。下面先介绍一下标准卡尔曼滤波算法。由于实际应用中噪声的不确定性，卡尔曼滤波算法采用了先预测再修正的思想将系统的状态方程和观测方程相结合，最终得到可实现状态最优估计的方程组。假设有动态系统模型： ...

【技术保护点】
一种抑制GNSS信息异常的滤波增益动态调整方法，其特征在于：包括以下步骤：S1：对标准卡尔曼滤波算法中的新息向量ik构建符合χ2分布的指标量I进行χ2检验，原假设H0:I～χ2(t)，指标量I如式(1)所示：I=ikT(HkPk,k-1HkT+Rk)-1ik---(1)]]>式(1)中，ik为新息向量，Hk为观测矩阵，Pk,k‑1为由k‑1时刻到k时刻的预测误差方差矩阵，Rk为观测噪声协方差矩阵；S2：选取α为显著性水平；S3：设通过查表方法得到α对应的上分位点为阈值；S4：对各个时刻的指标量I进行检测：如果连续n个时刻的n≥预设时间，则放弃GNSS，仅采用INS导航，进行步骤S7；否则，则进行步骤S5；S5：构建比例因子μ：μ=Ikχα2---(2)]]>式(2)中，Ik为新息向量矩阵；S6：用代替Kk重新进行滤波计算，得到对异常值抑制后的滤波结果，如式(3)所示：K‾k=Kk/μ---(3)]]>式(3)中，Kk是滤波器增益矩阵；S7：结束。

【技术特征摘要】
1.一种抑制GNSS信息异常的滤波增益动态调整方法，其特征在于：包括以下步骤：S1：对标准卡尔曼滤波算法中的新息向量ik构建符合χ2分布的指标量I进行χ2检验，原假设H0:I～χ2(t)，指标量I如式(1)所示： I = i k T ( H k P k , k - 1 H k T + R k ) - 1 i k - - - ( 1 ) ]]>式(1)中，ik为新...

【专利技术属性】
技术研发人员：王立辉，乔楠，余乐，张月新，
申请(专利权)人：东南大学，
类型：发明
国别省市：江苏;32