CN113687397B - 一种紧组合导航转发式欺骗干扰检测方法 - Google Patents
一种紧组合导航转发式欺骗干扰检测方法 Download PDFInfo
- Publication number
- CN113687397B CN113687397B CN202111050665.XA CN202111050665A CN113687397B CN 113687397 B CN113687397 B CN 113687397B CN 202111050665 A CN202111050665 A CN 202111050665A CN 113687397 B CN113687397 B CN 113687397B
- Authority
- CN
- China
- Prior art keywords
- detection
- deception
- time delay
- residual
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 113
- 238000001514 detection method Methods 0.000 claims abstract description 130
- 238000007781 pre-processing Methods 0.000 claims abstract description 12
- 238000002955 isolation Methods 0.000 claims description 42
- 239000011159 matrix material Substances 0.000 claims description 21
- 238000001914 filtration Methods 0.000 claims description 12
- 238000005259 measurement Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 238000013461 design Methods 0.000 abstract description 7
- 238000011897 real-time detection Methods 0.000 abstract description 3
- 230000000694 effects Effects 0.000 description 3
- 101100191408 Arabidopsis thaliana PRN2 gene Proteins 0.000 description 2
- 230000035772 mutation Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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和伪距率残差表达式为：ek＝zk-Hkxk/k-1；
步骤2.2：根据无欺骗干扰时存储的残差方差数据更新检测阈值，具体为：根据残差方差离散程度和计算量选择窗口容量，使用未受欺骗干扰时存储的伪距和伪距率残差方差数据，按设定的虚警概率对检测阈值进行实时更新；
步骤2.3：更新残差滑动窗口内数据，具体为：构造容量为m的滑动窗口W＝[ek-m+1,ek-m+2,…,ek]，利用前m-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：所述卫星在过去滤波过程内未受到欺骗干扰，执行常规检测程序。
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和伪距率残差表达式为：ek＝zk-Hkxk/k-1。
步骤2.2中，根据无欺骗干扰时存储的残差方差数据进行检测阈值更新方法：残差方差离散程度随窗口容量增大而减小，但是窗口容量过大会增加计算量。故窗口容量可以取50。在实际应用中需要使用未受欺骗干扰时存储的伪距和伪距率残差方差数据，按设定的虚警概率对检测阈值进行实时更新。
步骤2.3中，残差滑动窗口内数据的更新方法：构造容量为m的滑动窗口W＝[ek-m+1,ek-m+2,…,ek]，利用前m-1历元的残差可以得到滑动均值进而获得滑动方差，每个历元通过移位操作对滑动窗口内数据进行更新。
步骤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动态时延场景时延参数
步骤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 (3)
1.一种紧组合导航转发式欺骗干扰检测方法，其特征在于：包括以下步骤：
步骤1：采集GNSS接收机导航参数和SINS导航参数，对数据做预处理；
步骤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和伪距率残差表达式为：ek＝zk-Hkxk/k-1；
步骤2.2：根据无欺骗干扰时存储的残差方差数据更新检测阈值，具体为：根据残差方差离散程度和计算量选择窗口容量，使用未受欺骗干扰时存储的伪距和伪距率残差方差数据，按设定的虚警概率对检测阈值进行实时更新；
步骤2.3：更新残差滑动窗口内数据，具体为：构造容量为m的滑动窗口W＝[ek-m+1,ek-m+2,…,ek]，利用前m-1历元的残差得到滑动均值进而获得滑动方差，每个历元通过移位操作对滑动窗口内数据进行更新；
步骤3：对每一颗可见卫星进行判断：残差滑动窗口内的有效数据个数num少于滑动窗口容量m是否成立，若成立，则不具备检测条件，暂不进入检测程序，若全部可见卫星均不具备检测条件，则方法结束；若不成立，则根据卫星通道前一历元欺骗标志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：采取动态时延欺骗信号分段检测，依次进行动态时延过程检测和固定时延过程检测，具体为：
步骤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：得到欺骗干扰检测结果并对检出欺骗通道所对应的卫星信号予以剔除，得到欺骗干扰检测和组合导航定位结果；
步骤6：更新历元信息，并重复执行步骤1到5，直至处理完所有历元。
2.根据权利要求1所述的一种紧组合导航转发式欺骗干扰检测方法，其特征在于：步骤1所述采集GNSS接收机导航参数和SINS导航参数，对数据做预处理具体为：
步骤1.1：采集GNSS接收机导航参数，包括：时间标签、卫星PRN号、载波相位、伪距、多普勒频移、信号强度和星历数据；
步骤1.2：采集SINS导航参数，包括比力数据、角速度数据和每次SINS更新的时间标签；
步骤1.3：对采集到的数据做预处理，具体为：首先处理接收机观测数据得到观测信息，包括时间标签、可见卫星PRN号、伪距、伪距率；其次，通过提取星历参数进行卫星位置速度解算得到各颗可见卫星的三维位置、速度、钟差和频漂；最后将SINS导航参数与GNSS接收机导航参数按同一时间标签进行时间同步。
3.根据权利要求1所述的一种紧组合导航转发式欺骗干扰检测方法，其特征在于：步骤5所述得到欺骗干扰检测结果并对检出欺骗通道所对应的卫星信号予以剔除，得到欺骗干扰检测和组合导航定位结果具体为：隔离受欺骗干扰卫星通道，并重新构造相应的残差矩阵、量测矩阵、量测噪声协方差矩阵和增益矩阵，组合导航卡尔曼滤波结束、状态更新，并保存当前历元欺骗检测结果。
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 CN113687397A (zh) | 2021-11-23 |
| CN113687397B true 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) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2746813A1 (en) * | 2012-12-21 | 2014-06-25 | Astrium GmbH | Detection of spoofing of GNSS navigation signals |
| 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 스푸핑 기만공격을 이용한 안티드론 시스템 및 방법 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7623068B2 (en) * | 2006-05-16 | 2009-11-24 | Southwest Research Institute | Detection of deception signal used to deceive geolocation receiver of a satellite navigation system |
| US9366762B2 (en) * | 2011-10-27 | 2016-06-14 | Southwest Research Institute | Detection of deceptive navigation signal |
| KR101930354B1 (ko) * | 2013-11-04 | 2018-12-18 | 한국전자통신연구원 | 위성항법 수신기에서의 기만 신호 검출 장치 및 방법 |
-
2021
- 2021-09-08 CN CN202111050665.XA patent/CN113687397B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2746813A1 (en) * | 2012-12-21 | 2014-06-25 | Astrium GmbH | Detection of spoofing of GNSS navigation signals |
| 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 (6)
| Title |
|---|
| MEDLL辅助的GNSS/INS系统欺骗信号辨识方法;许睿;丁梦羽;孟骞;刘建业;;中国惯性技术学报(第02期);全文 * |
| 不同GPS/SINS超紧组合框架的分析与等价性推导;高帅和; 赵琳;中国惯性技术学报;第19卷(第5期);全文 * |
| 基于加权估计的紧组合导航抗欺骗算法研究;袁超;康国华;郑峰婴;徐川;张雷;华冰;;现代电子技术(第19期);全文 * |
| 基于滑动窗口的新息自适应组合导航算法;赵琳; 李久顺; 程建华;系统工程与电子技术;第39卷(第11期);全文 * |
| 袁超 ; 康国华 ; 郑峰婴 ; 徐川 ; 张雷 ; 华冰 ; .基于加权估计的紧组合导航抗欺骗算法研究.现代电子技术.2013,(第19期),全文. * |
| 赵琳 ; 李久顺 ; 程建华. 基于滑动窗口的新息自适应组合导航算法.系统工程与电子技术.2017,第39卷(第11期),全文. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113687397A (zh) | 2021-11-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108241159B (zh) | 基于仿生组合导航辅助的抗欺骗式干扰方法 | |
| US5917445A (en) | GPS multipath detection method and system | |
| CN107505642B (zh) | 一种ins辅助的实时bds单频周跳探测方法 | |
| JP5425039B2 (ja) | 衛星信号判定装置及びプログラム | |
| KR20140138027A (ko) | 다중 모드 내비게이션을 위한 수신기 및 방법 | |
| CN108974054B (zh) | 无缝列车定位方法及其系统 | |
| CN110673168B (zh) | 异步多用户联合欺骗信号检测方法及装置 | |
| CN113687397B (zh) | 一种紧组合导航转发式欺骗干扰检测方法 | |
| EP2208084A1 (en) | Survivor locating method and apparatus using search and rescue beacon equipped with navigation chipset | |
| CN111913195B (zh) | 基于陆基无线电导航信息的gps接收机抗欺骗干扰处理方法 | |
| CN112014863A (zh) | 基于多地面站数据的gnss功率变化时刻自动分析方法 | |
| WO1998049577A2 (en) | Methods for gyro bias estimation using gps | |
| CN116859415A (zh) | 一种快速稳健高精度的多故障卫星识别及定位方法 | |
| Sparks et al. | Extreme ionospheric storms and their impact on WAAS | |
| Konno et al. | Analysis of ionosphere gradient using Japan GEONET data | |
| EP2418513A1 (en) | Computing of robust and improved signal-in-space accuracy parameters in a regional or global navigation satellite system | |
| CN111999750B (zh) | 针对杆臂不准的实时单站周跳探测改进方法 | |
| van der Merwe et al. | Position solution sensitivity for snapshot receivers | |
| CN113534200A (zh) | 一种多模多频卫星导航授时接收机抗欺骗方法 | |
| Zair et al. | Coupling outlier detection with particle filter for GPS-based localization | |
| Yang et al. | Detection and exclusion of incipient fault for gnss-based train positioning under non-gaussian assumption | |
| Lee et al. | Cycle Slip Detection to Provide Stable Precision Positioning under Severe Ionospheric Conditions | |
| CN117538914A (zh) | 惯性辅助城市环境下的gnss多粗差探测方法 | |
| CN114325769B (zh) | 一种实时识别及剔除gnss转发欺骗干扰的方法 | |
| Giremus et al. | A GLR algorithm to detect and exclude up to two simultaneous range failures in a GPS/Galileo/IRS case |
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 |