惯性导航系统/测距机系统辅助的机载全球导航卫星系统欺骗式干扰自主检测算法
2022-09-22赵崴震钟伦珑
卢 丹 赵崴震 钟伦珑
(中国民航大学天津市智能信号与图像处理重点实验室 天津 300300)
1 引言
全球导航卫星系统(Global Navigation Satellite System, GNSS)的定位精度高和实时性好等优点使其在民航航行新技术中扮演着重要角色[1]。然而,GNSS信号功率低和民用粗码结构公开等因素导致其极易受到欺骗式干扰的影响[2,3]。欺骗式干扰使用假冒的卫星信号,致使机载接收机生成错误的定位结果,这严重影响了飞机的飞行安全。软件无线电技术(Software Defined Radio, SDR)与硬件平台的发展使得个人生成模拟GNSS信号将更加容易[4,5],未来蓄意破坏者对目标用户进行有意干扰的频率将大大增加。GNSS欺骗逐渐成为不容忽视的威胁。
针对欺骗式干扰,主要有信号功率监视[6]、密码认证[7-9]、天线阵列处理[10,11]和使用辅助导航信息[12]等检测方式。惯性导航系统(Inertial Navigation System, INS)是机载自主式导航设备,不受无线电信号的干扰,与GNSS有较高的互补性,是检测机载GNSS欺骗最常用的辅助导航工具[13,14]。基于卡尔曼滤波新息序列的欺骗检测算法,不需要对基础设施进行修改,也不需要增加额外的设备,利用飞机上已有的INS/GNSS数据即可实现检测,便于飞行过程中使用[15]。欺骗式干扰会导致INS/GNSS系统量测值发生幅值较大的阶跃变化,此时新息序列的均值不为0,通过该新息序列构造的卡方检验统计量会超过阈值,从而检测到欺骗式干扰的存在[16]。也就是说,卡方检验算法有意义的条件是构成卡方检验统计量的新息序列的均值能够在存在欺骗式干扰的情况下始终不为0。然而,当量测值突变幅度保持不变时,卡尔曼滤波的估计值和一步预测值会由于欺骗量测值的修正而逐渐和量测值“一致”起来,新息序列均值重新归0,检验统计量重新回到阈值之下,卡方检测算法失效,即使后面存在欺骗,检测算法也无法通过统计量与阈值的关系进行判断。因此,INS/GNSS系统无法依据传统卡方检测算法来判断给出的定位信息是否正确。另外,以上使用辅助导航信息进行欺骗检测的方法大多局限于INS/GNSS组合方式,并没有考虑到机载多导航源这一特点。
机载测距机(Distance Measuring Equipment,DME)是一种由地面应答机和机载询问机两部分组成的无线电测距导航系统。该系统通过测量脉冲信号在两者之间的传播时间来获取斜距信息,进而确定飞机的3维位置[17]。GNSS欺骗信号的功率比真实GNSS信号功率略高且在接收机热噪声下,DME信号的功率远高于接收机热噪声,因此DME系统不会受到欺骗式干扰的影响[18]。本文使用联邦卡尔曼滤波进行INS,GNSS和DME的数据融合,引入融合后的INS/DME数据参与构建欺骗检验统计量。当机载GNSS欺骗达到一定程度时,本文算法能够给出较为准确的欺骗持续时间。本文进一步将所设计的算法应用于组合导航系统中,得到了正确的定位信息,最后通过仿真实验对所提算法的有效性进行了验证。
2 组合导航模型与问题描述
本节分别给出了松组合方式下INS/GNSS系统和INS/DME系统的卡尔曼滤波的状态方程和量测方程,并对欺骗式干扰环境下的INS/GNSS量测方程进行量化分析,简要说明了传统检测算法的失效原因。
2.1 状态方程与量测方程的构建
使用导航系原点固定的INS直接输出北、东和下3个方向的位置信息,以INS的误差参数作为状态量,得到卡尔曼滤波的状态方程为
2.2 INS/GNSS量测方程分析
变,而此时的统计量会再次超过阈值。因此,当统计量第2次超过阈值后,不一定就是由于欺骗结束而导致,进而说明,传统算法无法依据统计量超过阈值的次数来判断欺骗是否结束。也就是说,传统算法存在检测不到欺骗持续时间的局限性。
3 欺骗检测算法
本文主要由传感器模块、欺骗检测模块和定位解算模块构成,如图1所示。传感器模块提供检测所必需的INS, GNSS和DME位置信息。欺骗检测模块给出了基于重构新息序列的有限记忆卡方检测算法的构造过程并证明了该方法的可行性。在此基础上,本文进一步利用DME数据,给出定位解算模块。该模块能够在存在GNSS欺骗的情况下提供正确的飞行引导。下面对新的检测算法进行详细介绍。
图1 新的检测算法原理图
3.1 重构新息序列
在式(13)的基础上进一步构造检验统计量,用于实现对欺骗开始和结束时间的较为准确的获取。下面给出有限记忆检验统计量的具体构造过程和相关的判决准则。
3.2 构建有限记忆卡方检验统计量
反之,则不存在欺骗式干扰。使用基于重构新息序列的有限记忆卡方检验统计量能够较为准确地判断欺骗的开始和结束,为使用GNSS还是DME进行准确的飞行引导提供了判断依据。
3.3 新的检测算法有效性证明
根据2.2节的分析可知,新息序列的均值是否为0,直接决定了后续构造的检验统计量能否检测到欺骗。下面给出新的检测算法的有效性证明。
将式(12)代入式(11),对新的量测值进行化简,得到
4 定位算法
5 仿真实验
本节通过仿真来验证新算法的有效性。模拟一段时间为180 s的直线飞行轨迹,初始位置为北纬40°,东经116°,高度1000 m。该实验使用GNSS模拟器生成GNSS信号,然后用软件接收机接收得到GNSS的位置信息。按照实际应用中对INS传感器精度的要求,设置陀螺仪随机噪声和加速度计随机噪声分别为0.07 °/h2和3 mg,设置陀螺仪和加速度计的1阶马尔可夫过程噪声分别为0.08 °/h2和12 mg,设置马尔可夫过程的相关时间为3600 s。根据真实DME测量精度,设置其测量误差为180 m。
本文设置欺骗的开始时间和结束时间分别为第40 s和第141 s。默认使用的欺骗信号在北、东、地3个方向上引起的偏差相同,以东方向位置误差为例,如图2所示。由图2可知,如果接收机接收真实GNSS信号,则接收机解算出来的位置误差均值始终为0。如果软件接收机从第40 s开始接收欺骗信号,可以看到,在第40~140 s这一时间段内,接收机解算的位置误差均值为300 m。由此说明,接收机接收欺骗信号解算的位置与接收真实GNSS信号解算位置存在300 m的偏差。
使用传统瞬时卡方算法对图2的欺骗场景进行检测,如图3所示。由图3可以看到,传统算法在第40 s超过阈值而报警,但在第53 s之后,传统算法的检验统计量重新回到阈值之下。第53~140 s这一时间段内,虽然存在GNSS欺骗但不会出现警报,直到第141 s欺骗结束引起量测值第2次发生突变,检验统计量再次超过阈值而告警。由此不难得出,传统算法给出的统计量与阈值的关系不能用来判断欺骗是否结束。
图2 接收机有、无欺骗的东向位置误差
图3 基于传统新息序列的瞬时卡方检测算法
在相同的实验条件下,使用本文设计的算法对图2的欺骗场景进行检测,如图4所示。滑动窗口的长度N设置为10,有限记忆卡方检验统计量在第41 s超过阈值而发出警报,表明该算法检测到了欺骗式干扰的存在。在第149 s检验统计量重新回到阈值之下,表明本文算法检测到了GNSS欺骗结束。由此可知,新的检测算法可以依靠统计量与阈值的关系来判断欺骗的开始和结束。
图4 基于重构新息序列的有限记忆卡方检测算法
图5分别给出了INS/GNSS系统位置误差、INS/DME系统定位误差和本文算法的误差曲线,其中3个子图分别表示北、东、地3个方向。由图可知,在没有欺骗式干扰时,INS/GNSS系统的误差曲线均值为0。第40~141 s这一时间段内,GNSS欺骗导致INS/GNSS系统的误差曲线被带到0均值以外。无论是否存在GNSS欺骗,INS/DME系统误差曲线始终保持0均值。然而,DME系统本身的定位误差远远大于GNSS,从而导致INS/DME的定位误差大于INS/GNSS。当本文所提的检测算法没有超过阈值时,使用INS/GNSS/DME作为组合导航系统的定位解。当检测到GNSS欺骗后,本文的定位算法立刻停止GNSS参与系统定位,因此红色的误差曲线能够始终保持0均值的统计规律。在GNSS系统隔离期间,本文的定位算法实质上就是INS/DME,组合导航系统的定位精度会发生下降。然而,无论是否存在欺骗式干扰,本文给出的定位方案都优于单一使用INS/GNSS或INS/DME。
图5 组合导航系统位置误差
下面分析不同DME测量精度下,欺骗引起的位置偏差对新算法检测概率的影响。接收机分别接收不同的欺骗信号,使其解算位置分别偏离真实GNSS位置20 m, 25 m, 30 m直到100 m。在存在GNSS欺骗的时间段内,如果出现检验统计量超过阈值的情况,则视该次实验检测到GNSS欺骗。进行1000次蒙特卡罗仿真,统计传统算法和本文算法检测到欺骗的次数,如图6所示。由图6可知,DME的定位精度越高,能够检测到的位置偏差量越小,检测概率越高。如果DME的测量误差为20 m,当3个方向都偏离20 m时,本文算法的检测概率就达到100%。相同实验条件下,传统算法需要当3个方向都偏离70 m时检测概率才能达到100%。由此说明,如果未来DME的定位精度提高或者存在其他定位精度更高的辅助导航源,所提算法都将有比现在更好的检测效果。
图6 不同DME测量误差下本文算法的检测概率
上面实验提到,新算法的滑动窗口长度设置为10,下面给出N= 10 的选取依据。N的选取需要综合考虑欺骗引起的位置偏差量、漏警概率和时间延迟量3个因素。接收机分别接收不同的欺骗信号,使其解算位置分别偏离真实GNSS位置200 m, 250 m,300 m直到500 m。在存在GNSS欺骗的时间段内,如果检验统计量出现低于阈值的情况,则视该次实验出现漏警。进行1000次蒙特卡罗实验,对位置偏差量和滑动窗口长度共同影响下的新算法的漏警概率进行统计,如图7所示。由图7可知,当本文算法的滑动窗口长度N设置为1时,其漏警概率在位置偏差为500 m时才能达到0%。相同实验条件下,当本文算法的滑动窗口长度N设置为10时,其漏警概率在加入250 m的GNSS欺骗时就能够达到0%。对图7进一步分析可知,如果GNSS欺骗为250 m及以上时,新算法的N≥10就可以保证漏警概率为0%。而对于200 m及以下的GNSS欺骗,需满足新算法的N≥20才可以达到0%的漏警概率。
图7 不同滑窗长度下本文算法的漏警概率
进一步分析不同程度的欺骗下,滑动窗口长度对新算法检测延迟的影响,如图8所示。由图8可知,当N/=1时,滑动窗口的存在会造成本文算法在欺骗检测上出现延迟。对于任意一个GNSS欺骗,滑动窗口的长度N越大,两幅子图的时间延迟量越大,所设计的检测算法给出的欺骗持续时间越不准确。假如欺骗式干扰至少会引起250 m的位置偏差,在满足漏警概率为0%的条件下,应让N尽量小。综合欺骗引起的位置偏差量、漏警概率和时间延迟量3个因素,本文选取滑动窗口的长度N为10。
图8 不同滑窗长度下本文算法的检测延迟
6 结束语
基于传统新息序列的卡方检测算法无法给出欺骗的持续时间,从而导致INS/GNSS系统无法根据该算法的检测结果判断其定位解算是否正确。本文提出一种基于重构新息序列的有限记忆卡方检测算法。本算法使用INS,GNSS和DME的数据构造一种不参与卡尔曼滤波的新息序列,然后将该新息序列构造成有限记忆卡方检验统计量,从而实现对欺骗式干扰的检测。仿真表明,当欺骗式干扰造成250 m及以上的位置偏差时,所设计的算法能够实现对欺骗持续时间较为准确的获取。本文进一步将所提算法应用于INS/GNSS/DME系统中,给出了正确的飞行引导方案。由检测概率的仿真图可知,如果未来DME的定位精度提高,或者存在其他定位精度更高的导航源,本文算法都将有比现在更好的检测效果和定位效果。