CN113985451A - 一种基于卡尔曼滤波跟踪环路的导航欺骗检测方法和装置 - Google Patents

一种基于卡尔曼滤波跟踪环路的导航欺骗检测方法和装置 Download PDF

Info

Publication number
CN113985451A
CN113985451A CN202111237967.8A CN202111237967A CN113985451A CN 113985451 A CN113985451 A CN 113985451A CN 202111237967 A CN202111237967 A CN 202111237967A CN 113985451 A CN113985451 A CN 113985451A
Authority
CN
China
Prior art keywords
code
navigation
carrier phase
phase variation
code phase
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
CN202111237967.8A
Other languages
English (en)
Other versions
CN113985451B (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.)
Xiangtan University
Original Assignee
Xiangtan 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 Xiangtan University filed Critical Xiangtan University
Priority to CN202111237967.8A priority Critical patent/CN113985451B/zh
Publication of CN113985451A publication Critical patent/CN113985451A/zh
Application granted granted Critical
Publication of CN113985451B publication Critical patent/CN113985451B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/21Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
    • G01S19/215Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service issues related to spoofing

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

本申请涉及一种卫星导航欺骗检测的方法及装置。所述方法包括:获取导航接收机射频前端输出的I/Q两路原始采样数据;对原始采样数据进行码域、频域二维搜索捕获可见星,获得用于跟踪环路的伪随机序列;提取用于滤波器的码相位差异、载波相位偏差和载波变化率;构建基于卡尔曼滤波的码跟踪和载波跟踪环路;建立无导航欺骗信号的零假设和存在导航欺骗信号的备择假设,通过对实时归一化新息平方的区间检测,判决当前接收卫星信号是否包含导航欺骗信号。本方法给出了一个较为完整可靠的导航欺骗信号检测方法,更适合现有单天线接收机硬件场景。

Description

一种基于卡尔曼滤波跟踪环路的导航欺骗检测方法和装置
技术领域
本申请涉及通信技术领域,特别是涉及一种基于卡尔曼滤波跟踪环路的导航欺骗检测方法和装置。
背景技术
全球导航卫星系统(Global Navigation Satellite System,GNSS)在获得广泛应用的同时,也存在一定的安全隐患。由于地面卫星信号功率微弱,民码结构完全公开,极大降低了GNSS的可靠性。早在2012年就有研究人员通过软件定义无线电(SDR)设备完成了多项对民用导航的欺骗干扰实验。随着电子战技术的发展,导航攻击方式将继续增加和多样化,研究利用GNSS信号特征进行导航环境监测的方法和系统具有重要的意义。
无论转发欺骗还是生成欺骗,都存在一定的信号特征差异,虽然已有通过信号功率、信号到达特性、相关峰检测、天线阵列、外部辅助等方法的检测方法,实现了信号能量、多峰、相关峰畸变、信号到达方向、信号到达时间和信号间相关性等检测方法。
但目前的方法存在一些不足:1、依靠欺骗信号时域和频域的对齐过程,对于精度较好的直接对齐检测可行性差。2、需要更改现有接收机的硬件方案,需要支持多天线接收,额外的原始数据传输链路。3、智能算法的应用,需要额外的计算单元完成对数据的计算。4、INS、气压计等通过外部辅助的方式,较大地增大了接收机成本。在此情况下,结合现有接收机硬件方案和数据处理流程,增加对某些特征的检测,同时不引入太高计算复杂度以实现对欺骗检测的方法更具适用性。
发明内容
基于此,有必要针对上述技术问题,提供一种满足低成本、较低计算复杂度、无需额外硬件支持,解决实际应用场景导航欺骗检测问题的特征值检测方法和装置。
一种基于卡尔曼滤波跟踪环路的导航欺骗检测方法,所述方法包括:
所述欺骗检测方法是建立无导航欺骗信号的零假设和存在导航欺骗信号的备择假设,将卡尔曼滤波跟踪环路器的归一化新息平方作为特征量,根据显著水平和其分布函数确定归一化新息平方的置信区间,通过对实时归一化新息平方的域值检测,获得导航信号中是否存在欺骗信号的判断。
所述卡尔曼滤波跟踪环路方法包括以下步骤:
S1,获取接收机I/Q两路中频采样数据SIF
S2,对采样数据SIF进行码域、频域二维搜索,获得可见星C/A码序列、可见星码相位时延和多普勒频移,作为跟踪参数向量初始值;
S3,构建基于载波相位变化量和码相位差异的卡尔曼滤波跟踪环路,得到系统状态方程和观测方程。
S4,通过跟踪参数向量生成本地C/A码序列,将所述1ms采样数据SIF与超前、即时、滞后C/A码序列做互相关运算得到相关器输出积分值,计算得到载波相位变化量以及码相位差异。
S5,根据所述相关器输出积分值更新可见星载噪比,并更新测量噪声协方差矩阵;
S6,将步骤S4所得载波相位变化量和码相位差异作为观测矩阵参数,代入卡尔曼滤波器进行融合计算,获得所述归一化新息平方值,并更新跟踪参数向量,得到最优估计载波相位变化量以及码相位变化量;
S7,通过域值检测判断是否接受原假设,如果是则重复步骤S4至S6,如果否则结束跟踪环路。
在步骤S2中,对采样数据SIF进行码域、频域二维搜索,获得可见星C/A码序列、可见星码相位时延和多普勒频移,作为跟踪参数向量初始值,包括:
根据以下公式进行二维搜索:
Figure BDA0003318094730000031
其中C为卫星对应C/A码序列,根据搜索结果获得间隔大于一个码宽的第一峰值与第二峰值,判断判断第一峰值与第二峰值的比值是否大于预设阈值,常用阈值为2dB,若是则得到可见星序号,其峰值所在对应可见星码相位时延及多普勒频移。若否则搜索其它卫星。
在步骤S3中,构建基于载波相位变化量和码相位差异的卡尔曼滤波跟踪环路,得到系统状态方程和观测方程,包括:
根据载波相位变化量和码相位变化量之间的线性系统关系,可以得到系统状态方程:
Figure BDA0003318094730000032
Figure BDA0003318094730000033
Figure BDA0003318094730000034
Figure BDA0003318094730000035
得到系统状态方程公式:
Xt=FXt-1+wt
其中状态向量X为由码相位差异Δτ、载波相位变化量Δθ、载波频率变化量Δf、载波频率变化率Δα组成的列向量,Xt-1为上一时刻状态,ωt为过程噪声,F为状态转移矩阵:
Figure BDA0003318094730000036
其中β为SIF采样速率与C/A码频率的比值,t为积分时间。
根据观测量与状态向量X的线性关系,可以得到观测方程:
Figure BDA0003318094730000041
Figure BDA0003318094730000042
得到观测方程公式:
Zt=HXt+vt
其中测量向量Zt为由观测码相位差异δτ、观测载波相位变化量δθ组成的列向量,vt为t时刻对应的测量误差,H为将状态向量映射到测量域的转换矩阵:
Figure BDA0003318094730000043
在步骤S4中,通过跟踪参数向量生成本地C/A码序列,将所述1ms采样数据SIF与超前、即时、滞后C/A码序列做互相关运算得到相关器输出积分值,计算得到载波相位变化量以及码相位差异,包括:
通过以下公式计算得到载波相位:
Figure BDA0003318094730000044
其中IP和Qp为即时码相关器输出积分值。
通过以下公式计算得到载波相位误差Δθ
θ=φe(t)-φe(t-1)
通过以下公式计算得到码相位差异:
Figure BDA0003318094730000045
其中d为超前、即时、滞后码的距离,IE和QE、IL和QL分别为超前码和滞后码相关器输出的积分值。
在步骤S5中,根据所述相关器输出积分值更新可见星载噪比,并更新测量噪声协方差矩阵,包括:
根据以下公式计算得到载噪比:
Figure BDA0003318094730000051
Figure BDA0003318094730000052
Figure BDA0003318094730000053
其中IP和Qp为即时码相关器输出积分值,T为积分时间。
根据以下公式计算得到测量噪声协方差:
Figure BDA0003318094730000054
Figure BDA0003318094730000055
其中
Figure BDA0003318094730000059
在步骤S6中,将步骤S4所得载波相位变化量和码相位差异作为观测矩阵参数,代入卡尔曼滤波器进行融合计算,获得所述归一化新息平方值,并更新跟踪参数向量,得到最优估计载波相位变化量以及码相位变化量,包括:
根据系统状态方程获得状态向量预测值
Figure BDA0003318094730000056
更具以下公式预测误差协方差矩阵:
Figure BDA0003318094730000057
其中Q为过程噪声协方差。
根据以下公式计算得到归一化新息平方:
Figure BDA0003318094730000058
其中
Figure BDA0003318094730000061
为测量噪声协方差。
根据以下公式计算得到卡尔曼增益:
Figure BDA0003318094730000062
根据以下公式计算得到状态向量的最优估计:
Figure BDA0003318094730000063
载波相位变化量以及码相位差异的最优估计位于状态向量
Figure BDA0003318094730000064
中。
根据以下公式更新误差协方差矩阵:
Figure BDA0003318094730000065
其中I为行列数均为4行的对角矩阵。
通过域值检测判断是否接受原假设,如果是则重复步骤S4至S6,如果否则结束跟踪环路,包括:
原假设不存在欺骗信号:
H0:E[εt]=m
如果归一化新息平方满足εt∈[r1,r2],则接受原假设。其中置信区间[r1,r2]定义为:
P(εt∈[r1 r2]|H0)=1-α
其中α为显著水平,归一化新息平方服从维度为m的卡方分布,
Figure BDA0003318094730000066
其中m为观测向量X的行数。通过判断每一时刻的归一化新息平方值是否处在置信区间内,判断是否接受原假设。如果超出置信区间则拒绝原假设,认为存在导航欺骗信号。
一种基于卡尔曼滤波跟踪环路的导航欺骗检测装置,所述装置包括:
信号采集模块,连续获取接收机I/Q两路中频采样数据SIF
捕获模块,对采样数据SIF进行码域、频域二维搜索,获得可见星C/A码序列、可见星码相位时延和多普勒频移,作为跟踪参数向量初始值;
相关器模块,所述1ms采样数据SIF与超前、即时、滞后C/A码序列做互相关运算得到相关器输出积分值;
载波相位鉴别器模块,用于基于相关器输出结果计算得到载波相位变化量。
码相位鉴别器模块,用于基于相关器输出结果计算得到码相位差异;
滤波器模块,将所述载波相位变化量和码相位差异作为观测矩阵参数,代入卡尔曼滤波器进行融合计算,获得所述归一化新息平方值,并更新跟踪参数向量,计算得到最优估计载波相位变化量以及码相位变化量;
判决器模块,用于利用归一化新息平方作为指标评估滤波器预测参数与计算参数的变化量,以判断当前导航环境是否存在欺骗信号。
上述基于卡尔曼滤波跟踪环路的导航欺骗检测方法和装置,根据卫星导航欺骗信号基本原理和导航接收机数字信号处理的具体过程,提出了一种适用于接收机跟踪环路状态评估的检测方法,并提取合适的跟踪参数,构建了环路卡尔曼滤波器,建立了适用于接收机跟踪环路的状态转移和观测模型。根据卡尔曼滤波器评估参数归一化新息平方的分布特点,给出了判定导航欺骗信号的假设检验准则。给出了一个完整可靠的基于卡尔曼滤波跟踪环路的导航欺骗检测方法。
附图说明
图1为本发明的基于卡尔曼滤波跟踪环路的导航欺骗检测方法的步骤流程示意图。
图2为本发明提供的卡尔曼滤波跟踪过程中归一化新息平方计算结果图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
在一个实施例中,参见图1至2所示,提供了一种基于卡尔曼滤波跟踪环路的导航欺骗检测方法,首先通过分析导航欺骗的实际过程和导航接收机的数据处理过程,提出了一种卡尔曼滤波跟踪环路的信号跟踪方法;分析了导航欺骗信号对环路滤波器工作状态的影响,研究了不同滤波参数下的检测效果,建立了用于构建完整跟踪环路的参数模型;最后用归一化新息平方作为评价指标,根据显著水平和其分布函数确定归一化新息平方的置信区间,通过对实时归一化新息平方的域值检测,实现判断导航欺骗信号是否存在的问题。具体可以包括以下步骤:
S202,对采样数据SIF进行码域、频域二维搜索,获得可见星C/A码序列、可见星码相位时延和多普勒频移,作为跟踪参数向量初始值;
根据卫星星座及卫星序号,生成一个周期C/A码序列,以采样SIF数据时同样的采样速率对其进行采样,得到采样后的C/A码序列。取与一个周期C/A码序列同样时长下的采样数据SIF,根据以下公式进行码域、频域二维搜索:
Figure BDA0003318094730000081
根据搜索结果获得第一峰值所在的码相位时延、多普勒频移和相关峰值,再在第一峰上下一个码片宽度之外进行二次搜索,获得第二峰的相关峰值。
根据第一相关峰值与第二相关峰值的比值是否大于预设阈值以判断该序号所对应的卫星信号是否存在。
S203,构建基于载波相位变化量和码相位差异的卡尔曼滤波跟踪环路,得到系统状态方程和观测方程。
根据载波相位变化量和码相位变化量之间的线性系统关系,可以得到系统状态方程:
Figure BDA0003318094730000091
Figure BDA0003318094730000092
Figure BDA0003318094730000093
Figure BDA0003318094730000094
得到系统状态方程公式:
Xt=FXt-1+wt
状态转移矩阵F:
Figure BDA0003318094730000095
根据观测量与状态向量X的线性关系,可以得到观测方程:
Figure BDA0003318094730000096
Figure BDA0003318094730000097
观测方程公式:
Zt=HXt+vt
H为将状态向量映射到测量域的转换矩阵:
Figure BDA0003318094730000098
S204,通过跟踪参数向量生成本地C/A码序列,将所述1ms采样数据SIF与超前、即时、滞后C/A码序列做互相关运算得到相关器输出积分值,计算得到载波相位变化量以及码相位差异。
通常超前、滞后码的间距为0.5个码片,即d=0.5,根据步骤S202同样生成采样后的三路C/A码序列,与所述1ms采样数据SIF经过互相关计算后,得各自的相干积分结果IE和QE、IP和QP、IL和QL
通过以下公式计算得到载波相位:
Figure BDA0003318094730000101
根据以下公式计算得到载波相位变化量Δθ
φ=φe(t)-φe(t-1)
根据以下公式计算得到码相位变化量:
Figure BDA0003318094730000102
S205,根据所述相关器输出积分值更新可见星载噪比,并更新测量噪声协方差矩阵;
根据以下公式计算得到载噪比:
Figure BDA0003318094730000103
Figure BDA0003318094730000104
Figure BDA0003318094730000105
根据以下公式计算得到测量噪声协方差:
Figure BDA0003318094730000106
Figure BDA0003318094730000107
其中
Figure BDA0003318094730000117
S206,将步骤S204所得载波相位变化量和码相位差异作为观测矩阵参数,代入卡尔曼滤波器进行融合计算,获得所述归一化新息平方值,并更新跟踪参数向量,得到最优估计载波相位变化量以及码相位变化量,;
根据系统状态方程获得状态向量预测值
Figure BDA0003318094730000111
根据以下公式预测误差协方差矩阵:
Figure BDA0003318094730000112
其中Q为过程噪声协方差,为行列数均为4的对角矩阵。
根据以下公式计算得到归一化新息平方:
εt=(Zt-HXt)T(HPtHT+R)-1(Zt-HXt)
根据以下公式计算得到卡尔曼增益:
Figure BDA0003318094730000113
根据以下公式计算得到状态向量的最优估计:
Figure BDA0003318094730000114
载波相位变化量以及码相位差异的最优估计位于状态向量
Figure BDA0003318094730000115
中。
根据以下公式更新误差协方差矩阵:
Figure BDA0003318094730000116
S207,通过域值检测判断是否接受原假设,如果是则重复步骤S204至S206,如果否则结束跟踪环路;
原假设不存在欺骗信号:
H0:E[εt]=m
如果归一化新息平方满足εt∈[r1,r2],则接受原假设。其中置信区间[r1,r2]定义为:
P(εt∈[r1 r2]|H0)=1-α
通过判断每一时刻的归一化新息平方值是否处在置信区间内,判断是否接受原假设,如果超出置信区间则认为存在导航欺骗信号。
如图2所示,在0~110s的时间段INS值处于置信区间中,在110s时加入了导航欺骗信号,导致INS值超出置信区间。
在一个实施例中提供了一种基于卡尔曼滤波跟踪环路的导航欺骗检测装置,包括:
信号采集模块,连续获取接收机I/Q两路中频采样数据S_IF;
捕获模块,对采样数据S_IF进行码域、频域二维搜索,获得可见星C/A码序列、可见星码相位时延和多普勒频移,作为跟踪参数向量初始值;
相关器模块,所述1ms采样数据S_IF与超前、即时、滞后C/A码序列做互相关运算得到相关器输出积分值;
载波相位鉴别器模块,用于基于相关器输出结果计算得到载波相位变化量。
码相位鉴别器模块,用于基于相关器输出结果计算得到码相位差异;
滤波器模块,将所述载波相位变化量和码相位差异作为观测矩阵参数,代入卡尔曼滤波器进行融合计算,获得所述归一化新息平方值,并更新跟踪参数向量,计算得到最优估计载波相位变化量以及码相位变化量;
判决器模块,用于利用归一化新息平方作为指标评估滤波器预测参数与计算参数的变化量,以判断当前导航环境是否存在欺骗信号。
关于一种基于卡尔曼滤波跟踪环路的导航欺骗检测装置的具体限定可以参见上文中对于一种基于卡尔曼滤波跟踪环路的导航欺骗检测法的限定,在此不再赘述。上述一种基于卡尔曼滤波跟踪环路的导航欺骗检测装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中,也可以以软件形式存储于计算机设备中的存储器中,以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中,提供了一种计算机设备,该计算机设备可以是终端。该计算机设备包括通过系统总线连接的处理器、存储器、通信接口、显示屏和输入装置。其中,该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的通信接口用于与外部的终端进行有线或无线方式的通信,无线方式可通过WIFI、运营商网络、NFC(近场通信)或其他技术实现。该计算机程序被处理器执行时以实现一种基于卡尔曼滤波跟踪环路的导航欺骗检测方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏,该计算机设备的输入装置可以是显示屏上覆盖的触摸层,也可以是计算机设备外壳上设置的按键、轨迹球或触控板,还可以是外接的键盘、触控板或鼠标等。
本领域技术人员可以理解,上述描述仅仅是与本申请方案相关的部分结构的框图,并不构成对本申请方案所应用于其上的计算机设备的限定,具体的计算机设备可以包括比图中所示更多或更少的部件,或者组合某些部件,或者具有不同的部件布置。
在一个实施例中,还提供了一种计算机设备,包括存储器和处理器,存储器中存储有计算机程序,该处理器执行计算机程序时实现上述各方法实施例中的步骤。
在一个实施例中,提供了一种计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现上述各方法实施例中的步骤。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用,均可包括非易失性和易失性存储器中的至少一种。非易失性存储器可包括只读存储器(Read-Only Memory,ROM)、磁带、软盘、闪存或光存储器等。易失性存储器可包括随机存取存储器(Random Access Memory,RAM)或外部高速缓冲存储器。作为说明而非局限,RAM可以是多种形式,比如静态随机存取存储器(Static Random Access Memory,SRAM)或动态随机存取存储器(Dynamic Random Access Memory,DRAM)等。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。

Claims (9)

1.一种基于卡尔曼滤波跟踪环路的导航欺骗检测方法,其特征在于:所述欺骗检测方法是建立无导航欺骗信号的零假设和存在导航欺骗信号的备择假设,将卡尔曼滤波跟踪环路器的归一化新息平方作为特征量,根据显著水平和其分布函数确定归一化新息平方的置信区间,通过对实时归一化新息平方的区间检测,获得导航信号中是否存在欺骗信号的判断。
2.根据权利要求1所述卡尔曼滤波跟踪环路的方法,其特征在于,包括以下步骤:
S1,获取接收机I/Q两路中频采样数据SIF
S2,对采样数据SIF进行码域、频域二维搜索,获得可见星C/A码序列、可见星码相位时延和多普勒频移,作为跟踪参数向量初始值;
S3,构建基于载波相位变化量和码相位差异的卡尔曼滤波跟踪环路,得到系统状态方程和观测方程。
S4,通过跟踪参数向量生成本地C/A码序列,将所述1ms采样数据SIF与超前、即时、滞后C/A码序列做互相关运算得到相关器输出积分值,计算得到载波相位变化量以及码相位差异。
S5,根据所述相关器输出积分值更新可见星载噪比,并更新测量噪声协方差矩阵;
S6,将步骤S4所得载波相位变化量和码相位差异作为观测矩阵参数,代入卡尔曼滤波器进行融合计算,获得所述归一化新息平方值,并更新跟踪参数向量,得到最优估计载波相位变化量以及码相位变化量;
S7,通过所述置信区间判断是否接受原假设,如果是则重复步骤S4至S6,如果否则结束跟踪环路。
3.根据权利要求2所述的方法,其特征在于,所述对采样数据SIF进行码域、频域二维搜索,获得可见星C/A码序列、可见星码相位时延和多普勒频移,作为跟踪参数向量初始值,包括:
根据以下公式进行二维搜索:
Figure FDA0003318094720000021
其中C为卫星对应C/A码序列,根据搜索结果获得间隔大于一个码宽的第一峰值与第二峰值,判断第一峰值与第二峰值的比值是否大于预设阈值,常用阈值为2dB,若是则得到可见星序号,其峰值所在对应可见星码相位时延及多普勒频移。若否则搜索其它卫星。
4.根据权利要求1所述的方法,其特征在于,所述构建基于载波相位变化量和码相位差异的卡尔曼滤波跟踪环路,得到系统状态方程和观测方程,包括:
系统状态方程公式:
Xt=FXt-1+wt
其中状态向量X为由码相位差异Δτ、载波相位变化量Δθ、载波频率变化量Δf、载波频率变化率Δα组成的列向量,Xt-1为上一时刻状态,ωt为过程噪声,F为状态转移矩阵:
Figure FDA0003318094720000022
其中β为SIF采样速率与C/A码频率的比值,t为积分时间。
观测方程公式:
Zt=HXt+vt
其中测量向量Zt为由观测数据提供,为观测码相位差异δτ、观测载波相位变化量δθ组成的列向量,H为将状态向量映射到测量域的转换矩阵:
Figure FDA0003318094720000031
5.根据权利要求2所述的方法,其特征在于,所述通过跟踪参数向量生成本地C/A码序列,将所述1ms采样数据SIF与超前、即时、滞后C/A码序列做互相关运算得到相关器输出积分值,计算得到载波相位变化量以及码相位差异,包括:
通过以下公式计算得到载波相位:
Figure FDA0003318094720000032
其中IP和Qp为即时码相关器输出积分值。
通过以下公式计算得到载波相位变化量Δθ
θ=φe(t)-φe(t-1)
通过以下公式计算得到码相位差异:
Figure FDA0003318094720000033
其中d为超前、即时、滞后码的距离,IE和QE、IL和QL分别为超前码和滞后码相关器输出的积分值。
6.根据权利要求6所述的方法,其特征在于,所述根据所述相关器输出积分值更新可见星载噪比,并更新测量噪声协方差矩阵,包括:
根据以下公式计算得到载噪比:
Figure FDA0003318094720000034
Figure FDA0003318094720000035
Figure FDA0003318094720000036
其中IP和Qp为即时码相关器输出积分值,
Figure FDA0003318094720000041
Figure FDA0003318094720000042
分别为接收信号中信号和总噪声功率的期望和方差,T为积分时间。
根据以下公式计算得到测量噪声协方差:
Figure FDA0003318094720000043
Figure FDA0003318094720000044
其中C/N0|*=10C/N0
7.根据权利要求2所述的方法,其特征在于,所述将步骤S4所得载波相位变化量和码相位差异作为观测矩阵参数,代入卡尔曼滤波器进行融合计算,获得所述归一化新息平方值,并更新跟踪参数向量,得到最优估计载波相位变化量以及码相位变化量,包括:
根据系统状态方程获得状态向量预测值
Figure FDA0003318094720000045
根据以下公式预测误差协方差矩阵:
Figure FDA0003318094720000046
其中Q为过程噪声协方差。
根据以下公式计算得到归一化新息平方:
Figure FDA0003318094720000047
其中
Figure FDA0003318094720000048
为测量噪声协方差。
根据以下公式计算得到卡尔曼增益:
Figure FDA0003318094720000049
根据以下公式计算得到状态向量的最优估计:
Figure FDA00033180947200000410
载波相位变化量以及码相位差异的最优估计位于状态向量
Figure FDA0003318094720000051
中。
根据以下公式更新误差协方差矩阵:
Figure FDA0003318094720000052
其中I为行列数均为4对角矩阵。
8.根据权利要求1所述的方法,其特征在于,所述通过域值检测判断是否接受原假设,如果是则重复步骤S4至S6,如果否则结束跟踪环路,包括:
原假设不存在欺骗信号:
H0:E[εt]=m
如果归一化新息平方满足εt∈[r1,r2],则接受原假设。其中置信区间[r1,r2]定义为:
P(εt∈[r1 r2]|H0)=1-α
其中α为显著水平,归一化新息平方服从维度为m的卡方分布,
Figure FDA0003318094720000053
其中m为观测向量X的行数。通过判断每一时刻的归一化新息平方值是否处在置信区间内,判断是否接受原假设。如果超出置信区间则拒绝原假设,认为存在导航欺骗信号。
9.一种基于卡尔曼滤波跟踪环路的导航欺骗检测装置,其特征在于,所述装置包括:
信号采集模块,连续获取接收机I/Q两路中频采样数据SIF
捕获模块,对采样数据SIF进行码域、频域二维搜索,获得可见星C/A码序列、可见星码相位时延和多普勒频移,作为跟踪参数向量初始值;
相关器模块,所述1ms采样数据SIF与超前、即时、滞后C/A码序列做互相关运算得到相关器输出积分值;
载波相位鉴别器模块,用于基于相关器输出结果计算得到载波相位变化量。
码相位鉴别器模块,用于基于相关器输出结果计算得到码相位差异;
滤波器模块,将所述载波相位变化量和码相位差异作为观测矩阵参数,代入卡尔曼滤波器进行融合计算,获得所述归一化新息平方值,并更新跟踪参数向量,计算得到最优估计载波相位变化量以及码相位变化量;
判决器模块,用于利用归一化新息平方作为指标评估滤波器预测参数与计算参数的变化量,以判断当前导航环境是否存在欺骗信号。
CN202111237967.8A 2021-10-25 2021-10-25 一种基于卡尔曼滤波跟踪环路的导航欺骗检测方法和装置 Active CN113985451B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111237967.8A CN113985451B (zh) 2021-10-25 2021-10-25 一种基于卡尔曼滤波跟踪环路的导航欺骗检测方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111237967.8A CN113985451B (zh) 2021-10-25 2021-10-25 一种基于卡尔曼滤波跟踪环路的导航欺骗检测方法和装置

Publications (2)

Publication Number Publication Date
CN113985451A true CN113985451A (zh) 2022-01-28
CN113985451B CN113985451B (zh) 2022-11-15

Family

ID=79740720

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111237967.8A Active CN113985451B (zh) 2021-10-25 2021-10-25 一种基于卡尔曼滤波跟踪环路的导航欺骗检测方法和装置

Country Status (1)

Country Link
CN (1) CN113985451B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116626716A (zh) * 2023-07-21 2023-08-22 齐鲁空天信息研究院 一种北斗信号的跟踪监测方法
CN117590429A (zh) * 2023-11-06 2024-02-23 中山大学·深圳 一种基于l型阵列的多欺骗信号来波方向检测方法及系统
CN117614789A (zh) * 2024-01-18 2024-02-27 浙江赛思电子科技有限公司 基于类卡尔曼无偏fir滤波器的载波相位跟踪方法及装置

Citations (34)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1031845A2 (en) * 1999-02-23 2000-08-30 Leica Geosystems Inc. Receiver calibration technique for glonass
US6496778B1 (en) * 2000-09-14 2002-12-17 American Gnc Corporation Real-time integrated vehicle positioning method and system with differential GPS
WO2003096061A1 (en) * 2002-05-13 2003-11-20 Honeywell International Inc. Methods and apparatus for radar data processing
US20060161329A1 (en) * 2005-01-14 2006-07-20 Robert Crane System and method for advanced tight coupling of GPS and inertial navigation sensors
EP1793240A1 (en) * 1997-04-21 2007-06-06 Raytheon Company Increased jamming immunity by optimizing processing gain for GPS/INS systems
US20070159392A1 (en) * 2006-01-06 2007-07-12 Honeywell International Inc. Method for improved signal to noise ratio estimation
US20070252754A1 (en) * 2006-04-28 2007-11-01 Alexander Steven B System and method for advanced tight coupling of GPS and navigation based on dead reckoning
DE102008043506A1 (de) * 2008-11-06 2010-05-12 Robert Bosch Gmbh Gerät zur Erkennung von GNSS-Störsendern
CN102183770A (zh) * 2011-05-20 2011-09-14 哈尔滨工程大学 一种抗多径干扰的gps伪随机码跟踪环路及其抗多径干扰方法
EP2397868A1 (en) * 2010-06-15 2011-12-21 The European Union, represented by the European Commission Method of providing an authenticable time-and-location indication
CN102353970A (zh) * 2011-06-10 2012-02-15 北京航空航天大学 一种高抗干扰性能gps/sins组合导航系统及实现方法
CN102540216A (zh) * 2011-11-14 2012-07-04 北京航空航天大学 一种自适应跟踪环路及实现方法
CN103116169A (zh) * 2013-01-20 2013-05-22 哈尔滨工程大学 一种基于矢量跟踪环路的抗干扰方法
CN104793220A (zh) * 2015-05-04 2015-07-22 中国电子科技集团公司第五十四研究所 一种基于多天线的欺骗干扰检测方法
CN105334522A (zh) * 2015-12-07 2016-02-17 北京奇虎科技有限公司 Gps攻击的检测方法及装置
CN105572700A (zh) * 2015-12-21 2016-05-11 北京航空航天大学 一种受干扰通道自检测gps矢量跟踪环路方法
CN106291591A (zh) * 2015-06-23 2017-01-04 霍尼韦尔国际公司 用载波相位和惯性传感器的全球导航卫星系统(gnss)欺骗检测
US20170023365A1 (en) * 2013-09-03 2017-01-26 Litel Instruments System and method for advanced navigation
WO2017067205A1 (zh) * 2015-10-20 2017-04-27 东南大学 一种基于空间域的联合非相干积分矢量跟踪方法
CN106772455A (zh) * 2016-12-29 2017-05-31 南京航空航天大学 一种基于惯性信息辅助与参数估计的gnss反欺骗环路跟踪方法
EP3203267A1 (fr) * 2016-02-05 2017-08-09 Thales Methode de calibrage d'un recepteur de radio-navigation par satellites
CN107765270A (zh) * 2016-08-17 2018-03-06 中国航空工业集团公司西安飞行自动控制研究所 一种基于卡尔曼滤波的卫星导航接收机跟踪环
CN108375779A (zh) * 2018-01-16 2018-08-07 北京理工大学 一种基于Kalman滤波器的导航接收机干扰抑制方法
CN108828628A (zh) * 2018-04-23 2018-11-16 厦门大学 一种欺骗信号检测方法
CN109188469A (zh) * 2018-08-01 2019-01-11 南京航空航天大学 一种gnss信号接收机欺骗信号参数估计方法
US20190129044A1 (en) * 2016-07-19 2019-05-02 Southeast University Cubature Kalman Filtering Method Suitable for High-dimensional GNSS/INS Deep Coupling
CN110531390A (zh) * 2019-07-30 2019-12-03 湘潭大学 多普勒积分法联合电离层残差法检测与修复周跳
EP3588467A1 (en) * 2018-06-21 2020-01-01 Eagle Technology, LLC Satellite automatic identification system (ais) for tracking a plurality of maritime vessels and related methods
WO2020000137A1 (en) * 2018-06-25 2020-01-02 Beijing DIDI Infinity Technology and Development Co., Ltd Integrated sensor calibration in natural scenes
CN111143918A (zh) * 2019-12-06 2020-05-12 湘潭大学 无线电波路径预测方法和装置
CN111650620A (zh) * 2020-05-29 2020-09-11 电子科技大学 一种基于gps导航的轨迹欺骗方法
CN112422227A (zh) * 2020-11-25 2021-02-26 湘潭大学 卫星可信导航干扰资源调度方法和装置
EP3805801A1 (en) * 2019-10-10 2021-04-14 HERE Global B.V. Identifying gnss navigation data as potentially manipulated or as trustworthy at least partially based on an estimated deviation of a second estimate of a satellite state from a first estimate of the satellite state
US20210311206A1 (en) * 2017-10-23 2021-10-07 Deere & Company RTK Vector Phase Locked Loop Architecture

Patent Citations (35)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1793240A1 (en) * 1997-04-21 2007-06-06 Raytheon Company Increased jamming immunity by optimizing processing gain for GPS/INS systems
EP1031845A2 (en) * 1999-02-23 2000-08-30 Leica Geosystems Inc. Receiver calibration technique for glonass
US6496778B1 (en) * 2000-09-14 2002-12-17 American Gnc Corporation Real-time integrated vehicle positioning method and system with differential GPS
WO2003096061A1 (en) * 2002-05-13 2003-11-20 Honeywell International Inc. Methods and apparatus for radar data processing
EP1504282A1 (en) * 2002-05-13 2005-02-09 Honeywell International Inc. Methods and apparatus for radar data processing
US20060161329A1 (en) * 2005-01-14 2006-07-20 Robert Crane System and method for advanced tight coupling of GPS and inertial navigation sensors
US20070159392A1 (en) * 2006-01-06 2007-07-12 Honeywell International Inc. Method for improved signal to noise ratio estimation
US20070252754A1 (en) * 2006-04-28 2007-11-01 Alexander Steven B System and method for advanced tight coupling of GPS and navigation based on dead reckoning
DE102008043506A1 (de) * 2008-11-06 2010-05-12 Robert Bosch Gmbh Gerät zur Erkennung von GNSS-Störsendern
EP2397868A1 (en) * 2010-06-15 2011-12-21 The European Union, represented by the European Commission Method of providing an authenticable time-and-location indication
CN102183770A (zh) * 2011-05-20 2011-09-14 哈尔滨工程大学 一种抗多径干扰的gps伪随机码跟踪环路及其抗多径干扰方法
CN102353970A (zh) * 2011-06-10 2012-02-15 北京航空航天大学 一种高抗干扰性能gps/sins组合导航系统及实现方法
CN102540216A (zh) * 2011-11-14 2012-07-04 北京航空航天大学 一种自适应跟踪环路及实现方法
CN103116169A (zh) * 2013-01-20 2013-05-22 哈尔滨工程大学 一种基于矢量跟踪环路的抗干扰方法
US20170023365A1 (en) * 2013-09-03 2017-01-26 Litel Instruments System and method for advanced navigation
CN104793220A (zh) * 2015-05-04 2015-07-22 中国电子科技集团公司第五十四研究所 一种基于多天线的欺骗干扰检测方法
CN106291591A (zh) * 2015-06-23 2017-01-04 霍尼韦尔国际公司 用载波相位和惯性传感器的全球导航卫星系统(gnss)欺骗检测
WO2017067205A1 (zh) * 2015-10-20 2017-04-27 东南大学 一种基于空间域的联合非相干积分矢量跟踪方法
CN105334522A (zh) * 2015-12-07 2016-02-17 北京奇虎科技有限公司 Gps攻击的检测方法及装置
CN105572700A (zh) * 2015-12-21 2016-05-11 北京航空航天大学 一种受干扰通道自检测gps矢量跟踪环路方法
EP3203267A1 (fr) * 2016-02-05 2017-08-09 Thales Methode de calibrage d'un recepteur de radio-navigation par satellites
US20190129044A1 (en) * 2016-07-19 2019-05-02 Southeast University Cubature Kalman Filtering Method Suitable for High-dimensional GNSS/INS Deep Coupling
CN107765270A (zh) * 2016-08-17 2018-03-06 中国航空工业集团公司西安飞行自动控制研究所 一种基于卡尔曼滤波的卫星导航接收机跟踪环
CN106772455A (zh) * 2016-12-29 2017-05-31 南京航空航天大学 一种基于惯性信息辅助与参数估计的gnss反欺骗环路跟踪方法
US20210311206A1 (en) * 2017-10-23 2021-10-07 Deere & Company RTK Vector Phase Locked Loop Architecture
CN108375779A (zh) * 2018-01-16 2018-08-07 北京理工大学 一种基于Kalman滤波器的导航接收机干扰抑制方法
CN108828628A (zh) * 2018-04-23 2018-11-16 厦门大学 一种欺骗信号检测方法
EP3588467A1 (en) * 2018-06-21 2020-01-01 Eagle Technology, LLC Satellite automatic identification system (ais) for tracking a plurality of maritime vessels and related methods
WO2020000137A1 (en) * 2018-06-25 2020-01-02 Beijing DIDI Infinity Technology and Development Co., Ltd Integrated sensor calibration in natural scenes
CN109188469A (zh) * 2018-08-01 2019-01-11 南京航空航天大学 一种gnss信号接收机欺骗信号参数估计方法
CN110531390A (zh) * 2019-07-30 2019-12-03 湘潭大学 多普勒积分法联合电离层残差法检测与修复周跳
EP3805801A1 (en) * 2019-10-10 2021-04-14 HERE Global B.V. Identifying gnss navigation data as potentially manipulated or as trustworthy at least partially based on an estimated deviation of a second estimate of a satellite state from a first estimate of the satellite state
CN111143918A (zh) * 2019-12-06 2020-05-12 湘潭大学 无线电波路径预测方法和装置
CN111650620A (zh) * 2020-05-29 2020-09-11 电子科技大学 一种基于gps导航的轨迹欺骗方法
CN112422227A (zh) * 2020-11-25 2021-02-26 湘潭大学 卫星可信导航干扰资源调度方法和装置

Non-Patent Citations (14)

* Cited by examiner, † Cited by third party
Title
JAFARNIA JAHROMI A: ""GPS spoofer countermeasure effectiveness based on signal strength, noise power, and C/NO measurements"", 《INTERNATIONAL JOURNAL OF SATELLITE COMMUNICATIONS AND NETWORKING》 *
O"HANLON B W: ""Real-time GPS spoofing detection via correlation of encrypted signals"", 《NAVIGATION》 *
WU Y H: ""Intelligent vehicle safety system based on BeiDou satellite navigation system"", 《IET INTELLIGENT TRANSPORT SYSTEMS》 *
XU Z: ""Analysis of GIS deformation monitoring accuracy of ultra-high voltage substation using Beidou satellite navigation system"", 《ZHEJIANG ELECTRIC POWER》 *
乔纯捷: ""全球导航卫星系统接收机的复信号自适应陷波干扰抑制"", 《国防科技大学学报》 *
任彬彬: ""GNSS调零抗干扰天线的反欺骗性能分析"", 《全球定位系统》 *
刘科: ""GNSS欺骗干扰检测算法与实验验证方法研究"", 《中国优秀博士学位论文全文数据库-信息科技辑》 *
张原: ""基于并行扩展卡尔曼滤波的自主卫星星座导航方法"", 《空间控制技术与应用》 *
李四海等: "惯性信息辅助的卫星导航欺骗检测技术", 《中国惯性技术学报》 *
杨瑛琳: ""基于容积卡尔曼滤波的GNSS抗欺骗式干扰研究"", 《中国优秀硕士学位论文全文数据库-信息科技辑》 *
王雪冬等: "基于Doppler估计的弹载接收机抗欺骗干扰研究", 《导航与控制》 *
许睿: ""欺骗环境下GNSS信号估计与定位修正技术"", 《航空学报》 *
谢党: "《导航对抗中信道模拟器的研究与实现》", 《中国优秀硕士学位论文全文数据库-信息科技辑》 *
顾明星: "".GNSS欺骗与抗欺骗技术研究进展"", 《江苏海洋大学学报》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116626716A (zh) * 2023-07-21 2023-08-22 齐鲁空天信息研究院 一种北斗信号的跟踪监测方法
CN116626716B (zh) * 2023-07-21 2023-09-29 齐鲁空天信息研究院 一种北斗信号的跟踪监测方法
CN117590429A (zh) * 2023-11-06 2024-02-23 中山大学·深圳 一种基于l型阵列的多欺骗信号来波方向检测方法及系统
CN117614789A (zh) * 2024-01-18 2024-02-27 浙江赛思电子科技有限公司 基于类卡尔曼无偏fir滤波器的载波相位跟踪方法及装置
CN117614789B (zh) * 2024-01-18 2024-04-09 浙江赛思电子科技有限公司 基于类卡尔曼无偏fir滤波器的载波相位跟踪方法及装置

Also Published As

Publication number Publication date
CN113985451B (zh) 2022-11-15

Similar Documents

Publication Publication Date Title
CN113985451B (zh) 一种基于卡尔曼滤波跟踪环路的导航欺骗检测方法和装置
US8416133B2 (en) System and method for compensating for faulty measurements
KR100941142B1 (ko) 데이터 통신망을 이용한 위치 추적 시스템 및 그 방법
US8775076B2 (en) Position calculating method and position calculating device
US6459405B1 (en) Satellite-based location system employing knowledge-based sequential signal search strategy
KR20010030317A (ko) 동적 적분 기술을 사용하는 위성에 기초한 위치 설정 시스템
CN104142508A (zh) 用于降低gnss接收器中功耗的方法、装置和系统
JP5636552B2 (ja) 移動受信機による、衛星からの拡散スペクトル信号の取得を最適化する方法
Guo et al. A new FFT acquisition scheme based on partial matched filter in GNSS receivers for harsh environments
CN108931789B (zh) 攻击检测方法、攻击检测器、计算机可读存储介质和终端
ES2613704T3 (es) Dispositivo espaciotemporal multi-antenas multi-correladores para el rechazo de los multi-trayecto de los sistemas de navegación
JP3809153B2 (ja) 衛星測位システム及び衛星測位方法
EP1545019A1 (en) GPS receiver using differential correlation
Shang et al. Detection and mitigation of GNSS spoofing via the pseudorange difference between epochs in a multicorrelator receiver
Xu et al. Estimation of ionospheric total electron content using GNSS observations derived from a smartphone
Tabatabaei et al. Robust adaptive joint tracking of GNSS signal code phases in urban canyons
US20110050496A1 (en) Energy Domain Based Peak Reconstruction Methods And Apparatuses
Andrianarison et al. Efficient and innovative techniques for collective acquisition of weak GNSS signals
CN106291612B (zh) 一种导航卫星星间链路无线信号高性能捕获判决方法
Stienne et al. GNSS dataless signal tracking with a delay semi-open loop and a phase open loop
Li et al. Rapid acquisition assisted by navigation data for inter-satellite links of navigation constellation
CN117269989A (zh) 一种基于ins辅助的GNSS欺骗的检测方法及系统
Siddakatte et al. Enhanced GNSS signal tracking in fading environments using frequency diversity
Jianfeng et al. Low C/N0 carrier tracking loop based on optimal estimation algorithm in GPS software receivers
Gioia et al. Precise time in your pocket: timing performance of Android phones

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