CN113687397A - 一种紧组合导航转发式欺骗干扰检测方法 - Google Patents

一种紧组合导航转发式欺骗干扰检测方法 Download PDF

Info

Publication number
CN113687397A
CN113687397A CN202111050665.XA CN202111050665A CN113687397A CN 113687397 A CN113687397 A CN 113687397A CN 202111050665 A CN202111050665 A CN 202111050665A CN 113687397 A CN113687397 A CN 113687397A
Authority
CN
China
Prior art keywords
detection
deception
time delay
flag
data
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
Application number
CN202111050665.XA
Other languages
English (en)
Other versions
CN113687397B (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.)
Harbin Engineering University
Original Assignee
Harbin Engineering 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 Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN202111050665.XA priority Critical patent/CN113687397B/zh
Publication of CN113687397A publication Critical patent/CN113687397A/zh
Application granted granted Critical
Publication of CN113687397B publication Critical patent/CN113687397B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/396—Determining accuracy or reliability of position or pseudorange measurements
    • 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/21—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
    • 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)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明公开了一种紧组合导航转发式欺骗干扰检测方法，包括(1)采集GNSS接收机导航参数和SINS导航参数，对数据做预处理；(2)更新残差滑动窗口内数据，根据无欺骗干扰时存储的残差方差数据更新检测阈值；(3)根据前一历元各通道受欺骗情况进入相应检测程序；(4)按照动态时延欺骗信号分段检测的设计思想，依次进行动态时延过程检测和固定时延过程检测；(5)得到欺骗干扰检测结果并对检出欺骗通道所对应的卫星信号予以剔除；(6)更新历元信息，并重复执行(1)到(5)，直至处理完所有历元。本发明实现了动态时延转发式欺骗干扰实时检测，能够检测码相位偏移速率不低于0.1chip/s的动态时延转发式欺骗干扰。

Description

一种紧组合导航转发式欺骗干扰检测方法
技术领域
本发明涉及一种GNSS/SINS紧组合导航转发式欺骗干扰检测，特别是一种基于残差动态时延特征的紧组合导航转发式欺骗干扰检测方法，属于卫星导航欺骗检测技术领域。
背景技术
随着全球导航卫星系统的快速发展，卫星导航定位技术的应用愈加广泛。然而接收机易在转发式欺骗干扰的作用下被诱导到错误的位置，这对接收机的可靠性提出了严峻挑战。
目前，关于转发式欺骗干扰检测的研究主要集中在无惯导的独立接收机导航和有惯导的组合导航两个方面。无惯导时，主要采用卫星信号功率检测、相关峰检测和伪距差异检测等方法，但是这些方法在欺骗信号功率低于真实信号，以及时间延迟缓慢增大的动态时延转发式欺骗干扰场景下失去作用。有惯导时，主要采用紧组合残差χ2检测、紧组合伪距差异检测等方法，但是前者不能指出受欺骗干扰影响的具体卫星，无法针对性地进行欺骗抑制；后者利用接收机与惯导伪距差值或紧组合伪距残差进行欺骗检测，能够检测出固定时延转发式欺骗干扰，但对于动态时延转发式欺骗干扰同样无能为力。
发明内容
针对上述现有技术，本发明要解决的技术问题是提供一种能够实现对动态时延转发式欺骗干扰的实时检测的紧组合导航转发式欺骗干扰检测方法。
为解决上述技术问题，本发明的一种紧组合导航转发式欺骗干扰检测方法，包括以下步骤：
步骤1：采集GNSS接收机导航参数和SINS导航参数，对数据做预处理；
步骤2：更新残差滑动窗口内数据，根据无欺骗干扰时存储的残差方差数据更新检测阈值；
步骤3：对每一颗可见卫星进行判断：残差滑动窗口内的有效数据个数num少于滑动窗口容量m是否成立，若成立，则不具备检测条件，暂不进入检测程序，若全部可见卫星均不具备检测条件，则方法结束；若不成立，则根据该卫星通道前一历元欺骗标志Dec_flag、隔离退出计数器cntr状态进入相应的检测程序；
步骤4：采取动态时延欺骗信号分段检测，依次进行动态时延过程检测和固定时延过程检测；
步骤5：得到欺骗干扰检测结果并对检出欺骗通道所对应的卫星信号予以剔除，得到欺骗干扰检测和组合导航定位结果；
步骤6：更新历元信息，并重复执行步骤1到5，直至处理完所有历元。
本发明还包括：
1.步骤2中更新残差滑动窗口内数据，根据无欺骗干扰时存储的残差方差数据更新检测阈值具体为：
步骤2.1：获得当前历元的伪距和伪距率残差，具体为：
分段线性化处理后的GNSS/SINS紧组合导航系统表达式为：xk＝Φk/k-1xk-1+Γk- 1wk-1和zk＝Hkxk+vk，式中，xk为状态向量；Φk/k-1为一步状态转移矩阵；Γk-1为系统噪声矩阵；zk为量测值向量；Hk为量测矩阵；wk-1和vk为互不相关的零均值白噪声向量，紧组合系统卡尔曼滤波残差ek包括伪距残差ρek和伪距率残差
Figure BDA0003252841770000021
表达式为：ek＝zk-Hkxk/k-1；
步骤2.2：根据无欺骗干扰时存储的残差方差数据更新检测阈值，具体为：根据残差方差离散程度和计算量选择窗口容量，使用未受欺骗干扰时存储的伪距和伪距率残差方差数据，按设定的虚警概率对检测阈值进行实时更新；
步骤2.3：更新残差滑动窗口内数据，具体为：构造容量为m的滑动窗口W＝[ek-m+1,ek-m+2,…,ek]，利用前m-1历元的残差得到滑动均值
Figure BDA0003252841770000022
进而获得滑动方差，每个历元通过移位操作对滑动窗口内数据进行更新。
2.步骤3中若不成立，则根据该卫星通道前一历元欺骗标志Dec_flag、隔离退出计数器cntr状态进入相应的检测程序具体为：
步骤3.1：Dec_flag＝0或1，若Dec_flag为1，则进入隔离检测程序，否则继续判断隔离退出计数器cntr是否为0，若cntr不为0，执行步骤3.2；若cntr为0，执行步骤3.3；
步骤3.2：所述卫星处于隔离退出阶段，执行隔离退出检测程序；
步骤3.3：所述卫星在过去滤波过程内未受到欺骗干扰，执行常规检测程序。
3.步骤4中采取动态时延欺骗信号分段检测，依次进行动态时延过程检测和固定时延过程检测具体为：
步骤4.1：动态时延过程检测，具体为：计算伪距率残差方差，并与伪距率残差方差阈值进行比较，若检测量大于阈值说明存在动态时延过程、存在动态时延转发式欺骗干扰，并将欺骗标志Dec_flag置为1、赋给隔离退出计数器cntr设定的最大计数值，然后结束；若检测量小于阈值说明不存在动态时延过程，执行步骤4.2；
步骤4.2：固定时延过程检测，具体为：计算伪距残差方差，并与伪距残差方差阈值进行比较，若检测量大于阈值说明存在固定时延过程、存在固定时延转发式欺骗干扰，并将欺骗标志Dec_flag置为1、赋给隔离退出计数器cntr设定的最大计数值，然后结束；若检测量小于阈值说明不存在转发式欺骗干扰，执行步骤4.3；
步骤4.3：未检出欺骗干扰时，对隔离退出计数器进行相应操作，具体为：执行常规检测程序，对隔离退出计数器归零，欺骗标志Dec_flag置为0；执行隔离退出检测程序，对隔离退出计数器采取每次减1的操作，直到隔离退出计数器归零，欺骗标志Dec_flag置为0；执行隔离检测程序，赋给隔离退出计数器设定的最大值，欺骗标志Dec_flag置为0。
4.步骤1中采集GNSS接收机导航参数和SINS导航参数，对数据做预处理具体为：
步骤1.1：采集GNSS接收机导航参数，包括：时间标签、卫星PRN号、载波相位、伪距、多普勒频移、信号强度和星历数据；
步骤1.2：采集SINS导航参数，包括比力数据、角速度数据和每次SINS更新的时间标签；
步骤1.3：对采集到的数据做预处理，具体为：首先处理接收机观测数据得到观测信息，包括时间标签、可见卫星PRN号、伪距、伪距率；其次，通过提取星历参数进行卫星位置速度解算得到各颗可见卫星的三维位置、速度、钟差和频漂；最后将SINS导航参数与GNSS接收机导航参数按同一时间标签进行时间同步。
5.步骤5中得到欺骗干扰检测结果并对检出欺骗通道所对应的卫星信号予以剔除，得到欺骗干扰检测和组合导航定位结果具体为：隔离受欺骗干扰卫星通道，并重新构造相应的残差矩阵、量测矩阵、量测噪声协方差矩阵和增益矩阵，组合导航卡尔曼滤波结束、状态更新，并保存当前历元欺骗检测结果。
本发明的有益效果：针对动态时延转发式欺骗干扰检测问题，本发明引入不受外界干扰的捷联式惯性导航系统组成GNSS/SINS紧组合导航系统，并在此基础上通过分析紧组合伪距和伪距率残差的动态时延特征，将一次完整的动态时延转发式欺骗干扰过程分解为“动态时延过程”和“固定时延过程”，提出了基于残差动态时延特征的紧组合导航转发式欺骗干扰检测方法，实现了对动态时延转发式欺骗干扰的实时检测。
本发明利用紧组合伪距和伪距率残差的动态时延特征设计检测量及检测流程，利用残差方差的概率统计特征设计检测阈值，该方法能够有效检测码相位偏移速率不低于0.1chip/s的动态时延转发式欺骗干扰。
附图说明
图1为本发明所述欺骗干扰检测方法整体流程图。
图2为常规检测程序流程图。
图3为隔离检测程序流程图。
图4为隔离退出检测程序流程图。
图5为无转发式欺骗干扰时伪距残差方差变化情况(滑动窗口容量m＝10)。
图6为无转发式欺骗干扰时伪距残差方差变化情况(滑动窗口容量m＝50)。
图7为无转发式欺骗干扰时伪距残差方差变化情况(滑动窗口容量m＝100)。
图8为无转发式欺骗干扰时伪距率残差方差变化情况(滑动窗口容量m＝10)。
图9为无转发式欺骗干扰时伪距率残差方差变化情况(滑动窗口容量m＝50)。
图10为无转发式欺骗干扰时伪距率残差方差变化情况(滑动窗口容量m＝100)。
图11为伪距残差方差概率统计特征。
图12为伪距率残差方差概率统计特征。
图13为动态时延转发式欺骗干扰伪距和伪距率误差变化情况。
图14为PRN1持续动态时延欺骗信号检测结果。
图15为PRN9阶段动态时延欺骗信号检测结果。
图16为PRN2无动态时延欺骗信号检测结果。
具体实施方式
下面结合说明书附图和具体实施方式对本发明做进一步说明。
结合图1至图4，本发明具体实施方式包括以下步骤：
步骤一：采集GNSS接收机导航参数和SINS导航参数，对数据做预处理，具体包括：
1.1采集GNSS接收机导航参数。
1.2采集SINS导航参数。
1.3对采集到的数据做预处理。
步骤二：更新残差滑动窗口内数据，根据无欺骗干扰时存储的残差方差数据更新检测阈值，具体包括：
2.1获得当前历元的伪距和伪距率残差。
2.2根据无欺骗干扰时存储的残差方差数据更新检测阈值
2.3更新残差滑动窗口内数据
步骤三：若残差滑动窗口内的有效数据个数num少于窗口容量m，则不具备检测条件，暂不进入检测程序；否则根据该卫星通道前一历元欺骗标志Dec_flag、隔离退出计数器cntr状态进入相应的检测程序，具体包括：
3.1若残差滑动窗口内的有效数据个数num少于窗口容量m，则暂不进入检测程序；否则根据前一历元各通道欺骗标志Dec_flag、隔离退出计数器cntr状态进入相应的检测程序。
3.2若Dec_flag为1，进入“隔离检测程序”，否则继续判断隔离退出计数器cntr是否为0。
3.3若cntr不为0，说明该卫星处于隔离退出阶段，需执行“隔离退出检测程序”。
3.4若cntr为0，说明该卫星在过去50次滤波过程内未受到欺骗干扰，需执行“常规检测程序”。
步骤四：按照动态时延欺骗信号分段检测的设计思想，依次进行动态时延过程检测和固定时延过程检测，具体包括：
4.1动态时延过程检测。
4.2固定时延过程检测。
4.3未检出欺骗干扰时，三种检测程序需对隔离退出计数器进行相应操作。
步骤五：得到欺骗干扰检测结果并对检出欺骗通道所对应的卫星信号予以剔除，具体包括：
5.1剔除检测出欺骗的通道及相应卫星。
5.2得到欺骗干扰检测和组合导航定位结果。
步骤六：更新历元信息，并重复执行步骤一到五，直至处理完所有历元。
步骤1.1中，采集GNSS接收机导航参数的方法：按RINEX3.02格式读取观测数据(包括时间标签、卫星PRN号、载波相位、伪距、多普勒频移、信号强度等观测值信息)和星历数据。
步骤1.2中，采集SINS输出的比力数据和角速度数据的方法：采集SINS输出的比力数据、角速度数据和每次SINS更新的时间标签。
步骤1.3中，对采集到的数据做预处理的方法：首先对接收机观测数据做简单处理，可直接得到时间标签、可见卫星PRN号、伪距、伪距率等观测信息；其次，通过提取星历参数、进行卫星位置速度解算可以得到各颗可见卫星的三维位置、速度以及钟差、频漂等信息；最后将SINS导航参数与GNSS接收机导航参数按同一时间标签进行时间同步。
步骤2.1中，当前历元的伪距和伪距率残差的获得方法：分段线性化处理后的GNSS/SINS紧组合导航系统表达式为：xk＝Φk/k-1xk-1+Γk-1wk-1和zk＝Hkxk+vk。式中，xk为状态向量；Φk/k-1为一步状态转移矩阵；Γk-1为系统噪声矩阵；zk为量测值向量；Hk为量测矩阵；wk-1和vk为互不相关的零均值白噪声向量。紧组合系统卡尔曼滤波残差ek包括伪距残差ρek和伪距率残差
Figure BDA0003252841770000051
表达式为：ek＝zk-Hkxk/k-1。
步骤2.2中，根据无欺骗干扰时存储的残差方差数据进行检测阈值更新方法：残差方差离散程度随窗口容量增大而减小，但是窗口容量过大会增加计算量。故窗口容量可以取50。在实际应用中需要使用未受欺骗干扰时存储的伪距和伪距率残差方差数据，按设定的虚警概率对检测阈值进行实时更新。
步骤2.3中，残差滑动窗口内数据的更新方法：构造容量为m的滑动窗口W＝[ek-m+1,ek-m+2,…,ek]，利用前m-1历元的残差可以得到滑动均值
Figure BDA0003252841770000052
进而获得滑动方差，每个历元通过移位操作对滑动窗口内数据进行更新。
步骤3.1中，判断是否进入检测程序的方法：判断残差滑动窗口内可用数据数量是否少于窗口容量，若少于则暂不进入检测程序，否则进入检测程序。
步骤3.2中，判断是否进入隔离检测程序的方法：判断欺骗标志是否为1，若为1则进入隔离检测程序，否则需要进一步判断隔离退出计数器是否为0。
步骤3.3中，判断进入隔离退出检测程序的方法：在欺骗标志为0的前提下判断隔离退出计数器是否为0，若不为0则进入隔离退出检测程序。
步骤3.4中，判断进入常规检测程序的方法：在欺骗标志为0的前提下判断隔离退出计数器是否为0，若为0则进入常规检测程序。
步骤4.1中，动态时延过程检测的方法：通过比较伪距和伪距率的动态时延特征可知在动态时延欺骗信号作用下，伪距误差变化缓慢、伪距率误差变化迅速。因此，可对动态时延转发式欺骗干扰进行分段检测，将动态时延转发式欺骗干扰过程等效为“动态时延过程”和“固定时延过程”。在动态时延过程中利用伪距率突变的特征进行检测，在固定时延过程中利用动态时延过程所积累的伪距差异进行检测，从而在不同阶段中分别利用伪距和伪距率的动态时延特征设计欺骗检测算法。
在动态时延过程检测阶段，计算伪距率残差方差，并与伪距率残差方差阈值进行比较，若检测量大于阈值说明存在动态时延过程、存在动态时延转发式欺骗干扰，并将欺骗标志Dec_flag置为1、赋给隔离退出计数器cntr最大计数值50；若检测量小于阈值说明不存在动态时延过程，但还需要进一步做固定时延过程检测以判断是否存在欺骗干扰。
步骤4.2中，固定时延过程检测的方法：在固定时延过程检测阶段，计算伪距残差方差，并与伪距残差方差阈值进行比较。若检测量大于阈值说明存在固定时延过程，说明存在转发式欺骗干扰，并将欺骗标志Dec_flag置为1、赋给隔离退出计数器cntr最大计数值50；若检测量小于阈值说明不存在转发式欺骗干扰。
步骤4.3中，未检出欺骗干扰时，三种检测程序对隔离退出计数器进行相应操作的方法：常规检测程序由于不涉及隔离退出的问题，故对隔离退出计数器直接归零；隔离退出检测程序需要确保刚刚恢复正常的卫星信号在最近50次滤波过程中无欺骗干扰影响，故对其每次采取“减1”的操作，直到隔离退出计数器归零；隔离检测程序未检出欺骗干扰，即判定当前卫星信号已恢复正常，并赋给隔离退出计数器最大计数值50，以便后续执行隔离退出程序。
步骤5.1中，剔除检测出欺骗的通道及相应卫星的方法：隔离受欺骗干扰卫星通道，并重新构造相应的残差矩阵、量测矩阵、量测噪声协方差矩阵和增益矩阵。
步骤5.2中，得到欺骗干扰检测和组合导航定位结果的方法：组合导航卡尔曼滤波结束、状态更新，并保存当前历元欺骗检测结果。
下面结合具体参数给出实施例：
实施例步骤如下：
步骤1：采集GNSS接收机导航参数和SINS导航参数，对数据做预处理。
步骤2：更新残差滑动窗口内数据，根据无欺骗干扰时存储的残差方差数据更新检测阈值。
步骤3：若残差滑动窗口内的有效数据个数num少于窗口容量m，则不具备检测条件，暂不进入检测程序；否则根据该卫星通道前一历元欺骗标志Dec_flag、隔离退出计数器cntr状态进入相应的检测程序。
步骤4：按照动态时延欺骗信号分段检测的设计思想，依次进行动态时延过程检测和固定时延过程检测。
步骤5：得到欺骗干扰检测结果并对检出欺骗通道所对应的卫星信号予以剔除。
步骤6：更新历元信息，并重复执行步骤(1)到(5)，直至处理完所有历元。
本发明以GNSS接收机输出的观测数据和中等精度光纤惯导输出的比力和角速度数据为例，对具体处理过程进行详细说明。数据中的动态时延转发式欺骗干扰为人为添加，该数据具体信息如表1所示：
表1动态时延场景时延参数
Figure BDA0003252841770000071
步骤1中，SINS数据总时长为600s；GNSS接收机数据总时长为590.9s；组合导航起始处理历元为第91个历元。
步骤2中，以PRN1为例，残差滑动窗口分别设置为10、50和100时，无转发式欺骗干扰下的伪距和伪距率残差方差变化情况如图5至图10所示，经综合考虑，残差滑动窗口容量设置为m＝50；图11和图12分别给出了无欺骗干扰影响时，伪距和伪距率残差方差的概率密度和累积分布曲线。当检测概率高于99％时，伪距残差检测阈值可设置为4.65m2，在该阈值下，虚警概率为1.5×10-3；当检测概率高于99％时，伪距率残差检测阈值可设置为0.10m2/s2，在该阈值下，虚警概率为1.2×10-3。
步骤3中，滑动窗口容量预设为m，若卫星序号i大于可见卫星数量n，说明全部可见卫星通道均完成检测，退出检测程序；否则继续执行，并对第i颗卫星的伪距和伪距率残差滑动窗口进行更新。若残差滑动窗口内的有效数据个数num少于窗口容量m，则暂不进入检测程序，否则根据前一历元各通道欺骗标志Dec_flag、隔离退出计数器cntr状态进入相应的检测程序。常规检测程序、隔离检测程序、隔离退出检测程序流程图如图2、图3、图4所示。当检测到欺骗标志Dec_flag＝1时，说明该卫星在上一次卡尔曼滤波过程中已被检出存在欺骗干扰，于是进入“隔离检测程序”进行欺骗干扰检测，否则需要进一步判断其“隔离退出计数器”cntr是否为0。若cntr不为0，则说明该卫星目前正处于隔离退出阶段，需要执行“隔离退出检测程序”；若cntr为0，则说明该卫星至少在过去50次滤波过程内未受到过欺骗干扰，只需执行“常规检测程序”。
步骤4中，以PRN1为例，动态时延转发式欺骗干扰下，伪距和伪距率误差的动态时延特征如图13所示。伪距误差以约14.7m/s的速率缓慢增大至约439.9m后不再继续增加；伪距率误差在2个历元内突变至约14.7m/s，并在伪距误差达到约439.9m后恢复至正常水平。通过比较伪距和伪距率的动态时延特征，可知将动态时延转发式欺骗干扰过程等效为“动态时延过程”和“固定时延过程”可以有效开展欺骗检测算法的设计。检测程序以图2为例。若MV_raw_dot>T1，说明存在动态时延过程，需令Dec_flag＝1和cntr＝50；若小于阈值T1，则需进一步判断是否存在固定时延过程；若MV_raw>T0，则说明存在固定时延过程，需令Dec_flag＝1和cntr＝50；若MV_raw<T0，则表明检测算法未检出欺骗干扰。
步骤5中，在检出欺骗干扰的前提下隔离受欺骗干扰卫星通道，并重新构造相应的残差矩阵、量测矩阵、量测噪声协方差矩阵和增益矩阵；完成状态更新，并保存当前历元欺骗检测结果。
步骤6中，更新历元信息KF_cntr＝KF_cntr+1，循环处理完所有历元。PRN1持续动态时延欺骗信号检测结果如图14所示。PRN9阶段动态时延欺骗信号检测结果如图15所示。PRN2无动态时延欺骗信号检测结果如图16所示。
由动态时延转发式欺骗干扰检测结果可知，本发明可以有效实时检测码相位偏移速率不低于0.1chip/s的动态时延转发式欺骗干扰。
综上，本发明涉及一种GNSS/SINS紧组合导航转发式欺骗干扰检测方法，属于卫星导航欺骗检测技术领域。传统检测方法根据转发式欺骗干扰造成的伪距突变对欺骗信号进行识别，能够检测出固定时延转发式欺骗干扰，但对于动态时延转发式欺骗干扰无能为力。针对动态时延转发式欺骗干扰检测问题，本发明提出了基于残差动态时延特征的紧组合导航转发式欺骗干扰检测方法，根据紧组合残差的动态时延特征选取检测量并设计检测流程，根据残差方差的概率统计特征设计检测阈值，实现了对动态时延转发式欺骗干扰的检测。本发明包括以下步骤：(1)采集GNSS接收机导航参数和SINS导航参数，对数据做预处理；(2)更新残差滑动窗口内数据，根据无欺骗干扰时存储的残差方差数据更新检测阈值；(3)根据前一历元各通道受欺骗情况进入相应的检测程序(4)按照动态时延欺骗信号分段检测的设计思想，依次进行动态时延过程检测和固定时延过程检测(5)得到欺骗干扰检测结果并对检出欺骗通道所对应的卫星信号予以剔除(6)更新历元信息，并重复执行步骤(1)到(5)，直至处理完所有历元。经验证，本发明可以有效实时检测码相位偏移速率不低于0.1chip/s的动态时延转发式欺骗干扰。

Claims (6)

1.一种紧组合导航转发式欺骗干扰检测方法，其特征在于：包括以下步骤：
步骤1：采集GNSS接收机导航参数和SINS导航参数，对数据做预处理；
步骤2：更新残差滑动窗口内数据，根据无欺骗干扰时存储的残差方差数据更新检测阈值；
步骤3：对每一颗可见卫星进行判断：残差滑动窗口内的有效数据个数num少于滑动窗口容量m是否成立，若成立，则不具备检测条件，暂不进入检测程序，若全部可见卫星均不具备检测条件，则方法结束；若不成立，则根据该卫星通道前一历元欺骗标志Dec_flag、隔离退出计数器cntr状态进入相应的检测程序；
步骤4：采取动态时延欺骗信号分段检测，依次进行动态时延过程检测和固定时延过程检测；
步骤5：得到欺骗干扰检测结果并对检出欺骗通道所对应的卫星信号予以剔除，得到欺骗干扰检测和组合导航定位结果；
步骤6：更新历元信息，并重复执行步骤1到5，直至处理完所有历元。
2.根据权利要求1所述的一种紧组合导航转发式欺骗干扰检测方法，其特征在于：步骤2所述更新残差滑动窗口内数据，根据无欺骗干扰时存储的残差方差数据更新检测阈值具体为：
步骤2.1：获得当前历元的伪距和伪距率残差，具体为：
分段线性化处理后的GNSS/SINS紧组合导航系统表达式为：xk＝Φk/k-1xk-1+Γk-1wk-1和zk＝Hkxk+vk，式中，xk为状态向量；Φk/k-1为一步状态转移矩阵；Γk-1为系统噪声矩阵；zk为量测值向量；Hk为量测矩阵；wk-1和vk为互不相关的零均值白噪声向量，紧组合系统卡尔曼滤波残差ek包括伪距残差ρek和伪距率残差
Figure FDA0003252841760000011
表达式为：ek＝zk-Hkxk/k-1；
步骤2.2：根据无欺骗干扰时存储的残差方差数据更新检测阈值，具体为：根据残差方差离散程度和计算量选择窗口容量，使用未受欺骗干扰时存储的伪距和伪距率残差方差数据，按设定的虚警概率对检测阈值进行实时更新；
步骤2.3：更新残差滑动窗口内数据，具体为：构造容量为m的滑动窗口W＝[ek-m+1,ek-m+2,…,ek]，利用前m-1历元的残差得到滑动均值ek-1，进而获得滑动方差，每个历元通过移位操作对滑动窗口内数据进行更新。
3.根据权利要求1或2所述的一种紧组合导航转发式欺骗干扰检测方法，其特征在于：步骤3所述若不成立，则根据该卫星通道前一历元欺骗标志Dec_flag、隔离退出计数器cntr状态进入相应的检测程序具体为：
步骤3.1：Dec_flag＝0或1，若Dec_flag为1，则进入隔离检测程序，否则继续判断隔离退出计数器cntr是否为0，若cntr不为0，执行步骤3.2；若cntr为0，执行步骤3.3；
步骤3.2：所述卫星处于隔离退出阶段，执行隔离退出检测程序；
步骤3.3：所述卫星在过去滤波过程内未受到欺骗干扰，执行常规检测程序。
4.根据权利要求1或2所述的一种紧组合导航转发式欺骗干扰检测方法，其特征在于：步骤4所述采取动态时延欺骗信号分段检测，依次进行动态时延过程检测和固定时延过程检测具体为：
步骤4.1：动态时延过程检测，具体为：计算伪距率残差方差，并与伪距率残差方差阈值进行比较，若检测量大于阈值说明存在动态时延过程、存在动态时延转发式欺骗干扰，并将欺骗标志Dec_flag置为1、赋给隔离退出计数器cntr设定的最大计数值，然后结束；若检测量小于阈值说明不存在动态时延过程，执行步骤4.2；
步骤4.2：固定时延过程检测，具体为：计算伪距残差方差，并与伪距残差方差阈值进行比较，若检测量大于阈值说明存在固定时延过程、存在固定时延转发式欺骗干扰，并将欺骗标志Dec_flag置为1、赋给隔离退出计数器cntr设定的最大计数值，然后结束；若检测量小于阈值说明不存在转发式欺骗干扰，执行步骤4.3；
步骤4.3：未检出欺骗干扰时，对隔离退出计数器进行相应操作，具体为：执行常规检测程序，对隔离退出计数器归零，欺骗标志Dec_flag置为0；执行隔离退出检测程序，对隔离退出计数器采取每次减1的操作，直到隔离退出计数器归零，欺骗标志Dec_flag置为0；执行隔离检测程序，赋给隔离退出计数器设定的最大值，欺骗标志Dec_flag置为0。
5.根据权利要求1或2所述的一种紧组合导航转发式欺骗干扰检测方法，其特征在于：步骤1所述采集GNSS接收机导航参数和SINS导航参数，对数据做预处理具体为：
步骤1.1：采集GNSS接收机导航参数，包括：时间标签、卫星PRN号、载波相位、伪距、多普勒频移、信号强度和星历数据；
步骤1.2：采集SINS导航参数，包括比力数据、角速度数据和每次SINS更新的时间标签；
步骤1.3：对采集到的数据做预处理，具体为：首先处理接收机观测数据得到观测信息，包括时间标签、可见卫星PRN号、伪距、伪距率；其次，通过提取星历参数进行卫星位置速度解算得到各颗可见卫星的三维位置、速度、钟差和频漂；最后将SINS导航参数与GNSS接收机导航参数按同一时间标签进行时间同步。
6.根据权利要求1或2所述的一种紧组合导航转发式欺骗干扰检测方法，其特征在于：步骤5所述得到欺骗干扰检测结果并对检出欺骗通道所对应的卫星信号予以剔除，得到欺骗干扰检测和组合导航定位结果具体为：隔离受欺骗干扰卫星通道，并重新构造相应的残差矩阵、量测矩阵、量测噪声协方差矩阵和增益矩阵，组合导航卡尔曼滤波结束、状态更新，并保存当前历元欺骗检测结果。
CN202111050665.XA 2021-09-08 2021-09-08 一种紧组合导航转发式欺骗干扰检测方法 Active CN113687397B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111050665.XA CN113687397B (zh) 2021-09-08 2021-09-08 一种紧组合导航转发式欺骗干扰检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111050665.XA CN113687397B (zh) 2021-09-08 2021-09-08 一种紧组合导航转发式欺骗干扰检测方法

Publications (2)

Publication Number Publication Date
CN113687397A true CN113687397A (zh) 2021-11-23
CN113687397B CN113687397B (zh) 2023-10-27

Family

ID=78585944

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111050665.XA Active CN113687397B (zh) 2021-09-08 2021-09-08 一种紧组合导航转发式欺骗干扰检测方法

Country Status (1)

Country Link
CN (1) CN113687397B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115042996A (zh) * 2022-04-07 2022-09-13 哈尔滨工程大学 一种多星同时段同区域观测任务规划方法

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080169978A1 (en) * 2006-05-16 2008-07-17 Luke Stephen Powell Detection Of Deception Signal Used To Deceive Geolocation Receiver Of A Satellite Navigation System
EP2746813A1 (en) * 2012-12-21 2014-06-25 Astrium GmbH Detection of spoofing of GNSS navigation signals
US20150070212A1 (en) * 2011-10-27 2015-03-12 Southwest Research Institute Detection of Deceptive Navigation Signal
US20150123846A1 (en) * 2013-11-04 2015-05-07 Electronics And Telecommunications Research Institute Apparatus and method for detecting deception signal in global navigation satellite receiver
KR101677136B1 (ko) * 2015-05-27 2016-11-17 국방과학연구소 단일 암호화 신호원을 이용한 위성항법 기만 검출 시스템 및 방법
CN106842238A (zh) * 2017-01-19 2017-06-13 中国民航大学 基于扩展raim的卫星导航欺骗干扰抑制方法
CN106980130A (zh) * 2017-03-03 2017-07-25 哈尔滨工程大学 一种sins/gnss深组合自适应导航方法
CN110715660A (zh) * 2019-11-21 2020-01-21 中国人民解放军63961部队 基于卡尔曼滤波估计的北斗/惯性组合抗欺骗干扰方法
CN111624626A (zh) * 2020-04-29 2020-09-04 西南电子技术研究所（中国电子科技集团公司第十研究所） 卫星导航系统欺骗干扰信号抑制方法
CN111913195A (zh) * 2020-07-07 2020-11-10 北京自动化控制设备研究所 基于陆基无线电导航信息的gps接收机抗欺骗干扰处理方法
KR102194734B1 (ko) * 2020-03-24 2020-12-23 국방과학연구소 Gps 스푸핑 기만공격을 이용한 안티드론 시스템 및 방법

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080169978A1 (en) * 2006-05-16 2008-07-17 Luke Stephen Powell Detection Of Deception Signal Used To Deceive Geolocation Receiver Of A Satellite Navigation System
US20150070212A1 (en) * 2011-10-27 2015-03-12 Southwest Research Institute Detection of Deceptive Navigation Signal
EP2746813A1 (en) * 2012-12-21 2014-06-25 Astrium GmbH Detection of spoofing of GNSS navigation signals
US20150123846A1 (en) * 2013-11-04 2015-05-07 Electronics And Telecommunications Research Institute Apparatus and method for detecting deception signal in global navigation satellite receiver
KR101677136B1 (ko) * 2015-05-27 2016-11-17 국방과학연구소 단일 암호화 신호원을 이용한 위성항법 기만 검출 시스템 및 방법
CN106842238A (zh) * 2017-01-19 2017-06-13 中国民航大学 基于扩展raim的卫星导航欺骗干扰抑制方法
CN106980130A (zh) * 2017-03-03 2017-07-25 哈尔滨工程大学 一种sins/gnss深组合自适应导航方法
CN110715660A (zh) * 2019-11-21 2020-01-21 中国人民解放军63961部队 基于卡尔曼滤波估计的北斗/惯性组合抗欺骗干扰方法
KR102194734B1 (ko) * 2020-03-24 2020-12-23 국방과학연구소 Gps 스푸핑 기만공격을 이용한 안티드론 시스템 및 방법
CN111624626A (zh) * 2020-04-29 2020-09-04 西南电子技术研究所（中国电子科技集团公司第十研究所） 卫星导航系统欺骗干扰信号抑制方法
CN111913195A (zh) * 2020-07-07 2020-11-10 北京自动化控制设备研究所 基于陆基无线电导航信息的gps接收机抗欺骗干扰处理方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
袁超;康国华;郑峰婴;徐川;张雷;华冰;: "基于加权估计的紧组合导航抗欺骗算法研究", 现代电子技术, no. 19 *
许睿;丁梦羽;孟骞;刘建业;: "MEDLL辅助的GNSS/INS系统欺骗信号辨识方法", 中国惯性技术学报, no. 02 *
赵琳; 李久顺; 程建华: "基于滑动窗口的新息自适应组合导航算法", 系统工程与电子技术, vol. 39, no. 11 *
赵琳; 李久顺; 程建华: "基于滑动窗口的新息自适应组合导航算法", 系统工程与电子技术,, vol. 39, no. 11 *
高帅和; 赵琳: "不同GPS/SINS超紧组合框架的分析与等价性推导", 中国惯性技术学报, vol. 19, no. 5 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115042996A (zh) * 2022-04-07 2022-09-13 哈尔滨工程大学 一种多星同时段同区域观测任务规划方法

Also Published As

Publication number Publication date
CN113687397B (zh) 2023-10-27

Similar Documents

Publication Publication Date Title
CN108241159B (zh) 基于仿生组合导航辅助的抗欺骗式干扰方法
US5917445A (en) GPS multipath detection method and system
CN101629997A (zh) 惯性辅助卫星导航完好性检测装置及检测方法
JP2008506126A5 (zh)
CN111427070B (zh) 一种gnss抗欺骗干扰方法
CN113009520B (zh) 一种卫星导航矢量跟踪环的完好性检测方法
CN115420284B (zh) 一种组合导航系统故障检测与识别方法
CN113687397A (zh) 一种紧组合导航转发式欺骗干扰检测方法
CN114325769B (zh) 一种实时识别及剔除gnss转发欺骗干扰的方法
CN111060935A (zh) Gnss欺骗干扰检测方法
WO2015165908A1 (fr) Procédé et dispositif de contrôle d&#39;intégrité à double niveau de consolidation
CN113109843B (zh) 基于双接收机伪距双差的欺骗信号检测抑制方法及装置
CN117607910B (zh) 一种基于矢量跟踪新息矢量的欺骗检测方法及系统
Mosavi et al. Positioning improvement by combining GPS and GLONASS based on Kalman filter and its application in GPS spoofing situations
CN111913195B (zh) 基于陆基无线电导航信息的gps接收机抗欺骗干扰处理方法
CN112014863A (zh) 基于多地面站数据的gnss功率变化时刻自动分析方法
Marais et al. Multipath and NLOS detection based on the combination of CN0 values and a fish-eye camera
EP1970723A1 (en) Method for the computation of improved SBAS protection levels valid in non-controlled environments
CN115540907A (zh) 一种基于面向星间差分的gps/bds/ins紧组合导航的多故障检测与剔除方法
CN115291253A (zh) 一种基于残差检测的车辆定位完好性监测方法及系统
CN112558112B (zh) 一种gnss状态域慢变斜坡故障完好性监测方法
CN111999750B (zh) 针对杆臂不准的实时单站周跳探测改进方法
van der Merwe et al. Position solution sensitivity for snapshot receivers
CN114814893A (zh) 一种卫星导航信号多址干扰检测方法
CN112882068A (zh) 基于多接收机的gnss抗欺骗干扰方法

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