CN113985451A - 一种基于卡尔曼滤波跟踪环路的导航欺骗检测方法和装置 - Google Patents
一种基于卡尔曼滤波跟踪环路的导航欺骗检测方法和装置 Download PDFInfo
- 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
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/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
- G01S19/215—Interference 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码序列、可见星码相位时延和多普勒频移,作为跟踪参数向量初始值,包括:
根据以下公式进行二维搜索:
其中C为卫星对应C/A码序列,根据搜索结果获得间隔大于一个码宽的第一峰值与第二峰值,判断判断第一峰值与第二峰值的比值是否大于预设阈值,常用阈值为2dB,若是则得到可见星序号,其峰值所在对应可见星码相位时延及多普勒频移。若否则搜索其它卫星。
在步骤S3中,构建基于载波相位变化量和码相位差异的卡尔曼滤波跟踪环路,得到系统状态方程和观测方程,包括:
根据载波相位变化量和码相位变化量之间的线性系统关系,可以得到系统状态方程:
得到系统状态方程公式:
Xt=FXt-1+wt
其中状态向量X为由码相位差异Δτ、载波相位变化量Δθ、载波频率变化量Δf、载波频率变化率Δα组成的列向量,Xt-1为上一时刻状态,ωt为过程噪声,F为状态转移矩阵:
其中β为SIF采样速率与C/A码频率的比值,t为积分时间。
根据观测量与状态向量X的线性关系,可以得到观测方程:
得到观测方程公式:
Zt=HXt+vt
其中测量向量Zt为由观测码相位差异δτ、观测载波相位变化量δθ组成的列向量,vt为t时刻对应的测量误差,H为将状态向量映射到测量域的转换矩阵:
在步骤S4中,通过跟踪参数向量生成本地C/A码序列,将所述1ms采样数据SIF与超前、即时、滞后C/A码序列做互相关运算得到相关器输出积分值,计算得到载波相位变化量以及码相位差异,包括:
通过以下公式计算得到载波相位:
其中IP和Qp为即时码相关器输出积分值。
通过以下公式计算得到载波相位误差Δθ:
△θ=φe(t)-φe(t-1)
通过以下公式计算得到码相位差异:
其中d为超前、即时、滞后码的距离,IE和QE、IL和QL分别为超前码和滞后码相关器输出的积分值。
在步骤S5中,根据所述相关器输出积分值更新可见星载噪比,并更新测量噪声协方差矩阵,包括:
根据以下公式计算得到载噪比:
其中IP和Qp为即时码相关器输出积分值,T为积分时间。
根据以下公式计算得到测量噪声协方差:
在步骤S6中,将步骤S4所得载波相位变化量和码相位差异作为观测矩阵参数,代入卡尔曼滤波器进行融合计算,获得所述归一化新息平方值,并更新跟踪参数向量,得到最优估计载波相位变化量以及码相位变化量,包括:
其中Q为过程噪声协方差。
根据以下公式计算得到归一化新息平方:
根据以下公式计算得到卡尔曼增益:
根据以下公式计算得到状态向量的最优估计:
根据以下公式更新误差协方差矩阵:
其中I为行列数均为4行的对角矩阵。
通过域值检测判断是否接受原假设,如果是则重复步骤S4至S6,如果否则结束跟踪环路,包括:
原假设不存在欺骗信号:
H0:E[εt]=m
如果归一化新息平方满足εt∈[r1,r2],则接受原假设。其中置信区间[r1,r2]定义为:
P(εt∈[r1 r2]|H0)=1-α
其中α为显著水平,归一化新息平方服从维度为m的卡方分布,其中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,根据以下公式进行码域、频域二维搜索:
根据搜索结果获得第一峰值所在的码相位时延、多普勒频移和相关峰值,再在第一峰上下一个码片宽度之外进行二次搜索,获得第二峰的相关峰值。
根据第一相关峰值与第二相关峰值的比值是否大于预设阈值以判断该序号所对应的卫星信号是否存在。
S203,构建基于载波相位变化量和码相位差异的卡尔曼滤波跟踪环路,得到系统状态方程和观测方程。
根据载波相位变化量和码相位变化量之间的线性系统关系,可以得到系统状态方程:
得到系统状态方程公式:
Xt=FXt-1+wt
状态转移矩阵F:
根据观测量与状态向量X的线性关系,可以得到观测方程:
观测方程公式:
Zt=HXt+vt
H为将状态向量映射到测量域的转换矩阵:
S204,通过跟踪参数向量生成本地C/A码序列,将所述1ms采样数据SIF与超前、即时、滞后C/A码序列做互相关运算得到相关器输出积分值,计算得到载波相位变化量以及码相位差异。
通常超前、滞后码的间距为0.5个码片,即d=0.5,根据步骤S202同样生成采样后的三路C/A码序列,与所述1ms采样数据SIF经过互相关计算后,得各自的相干积分结果IE和QE、IP和QP、IL和QL。
通过以下公式计算得到载波相位:
根据以下公式计算得到载波相位变化量Δθ:
△φ=φe(t)-φe(t-1)
根据以下公式计算得到码相位变化量:
S205,根据所述相关器输出积分值更新可见星载噪比,并更新测量噪声协方差矩阵;
根据以下公式计算得到载噪比:
根据以下公式计算得到测量噪声协方差:
S206,将步骤S204所得载波相位变化量和码相位差异作为观测矩阵参数,代入卡尔曼滤波器进行融合计算,获得所述归一化新息平方值,并更新跟踪参数向量,得到最优估计载波相位变化量以及码相位变化量,;
其中Q为过程噪声协方差,为行列数均为4的对角矩阵。
根据以下公式计算得到归一化新息平方:
εt=(Zt-HXt)T(HPtHT+R)-1(Zt-HXt)
根据以下公式计算得到卡尔曼增益:
根据以下公式计算得到状态向量的最优估计:
根据以下公式更新误差协方差矩阵:
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,如果否则结束跟踪环路。
9.一种基于卡尔曼滤波跟踪环路的导航欺骗检测装置,其特征在于,所述装置包括:
信号采集模块,连续获取接收机I/Q两路中频采样数据SIF;
捕获模块,对采样数据SIF进行码域、频域二维搜索,获得可见星C/A码序列、可见星码相位时延和多普勒频移,作为跟踪参数向量初始值;
相关器模块,所述1ms采样数据SIF与超前、即时、滞后C/A码序列做互相关运算得到相关器输出积分值;
载波相位鉴别器模块,用于基于相关器输出结果计算得到载波相位变化量。
码相位鉴别器模块,用于基于相关器输出结果计算得到码相位差异;
滤波器模块,将所述载波相位变化量和码相位差异作为观测矩阵参数,代入卡尔曼滤波器进行融合计算,获得所述归一化新息平方值,并更新跟踪参数向量,计算得到最优估计载波相位变化量以及码相位变化量;
判决器模块,用于利用归一化新息平方作为指标评估滤波器预测参数与计算参数的变化量,以判断当前导航环境是否存在欺骗信号。
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)
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)
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 |
-
2021
- 2021-10-25 CN CN202111237967.8A patent/CN113985451B/zh active Active
Patent Citations (35)
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 (5)
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 |