CN108535745A - 一种gnss阵列接收机信号解扩前的欺骗干扰检测方法 - Google Patents

一种gnss阵列接收机信号解扩前的欺骗干扰检测方法 Download PDF

Info

Publication number
CN108535745A
CN108535745A CN201810259165.9A CN201810259165A CN108535745A CN 108535745 A CN108535745 A CN 108535745A CN 201810259165 A CN201810259165 A CN 201810259165A CN 108535745 A CN108535745 A CN 108535745A
Authority
CN
China
Prior art keywords
signal
probability
detection
array
cheating
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
CN201810259165.9A
Other languages
English (en)
Other versions
CN108535745B (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.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN201810259165.9A priority Critical patent/CN108535745B/zh
Publication of CN108535745A publication Critical patent/CN108535745A/zh
Application granted granted Critical
Publication of CN108535745B publication Critical patent/CN108535745B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

本发明属于卫星导航技术领域,具体涉及一种卫星导航信号的欺骗干扰的检测方法。主要包括以下步骤:(S1)信号功率估计,得到阵元接收信号的功率估计结果;(S2)估计不存在欺骗干扰时功率估计结果实部与虚部的概率分布;(S3)给定虚警概率,计算检测门限;(S4)将信号功率估计结果与检测门限进行比较,实现欺骗干扰的检测。本发明检测性能良好且实现简单,不需要捕获和跟踪接收到的卫星信号,运算量小且实施方便,可直接用于GNSS阵列接收机基带采样信号的欺骗干扰检测。

Description

一种GNSS阵列接收机信号解扩前的欺骗干扰检测方法
技术领域
本发明属于卫星导航技术领域,具体涉及一种卫星导航信号的欺骗干扰的检测方法。
背景技术
全球导航卫星系统(Global Navigation Satellite System,GNSS)欺骗干扰隐蔽性强,对GNSS的应用造成了严重的安全威胁。而欺骗干扰信号与真实信号具有相同的信号结构,且功率较小,通常需要在信号解扩后得到信号的相应特征,如信号功率,载波多普勒,载波相位等才能进行检测。而信号解扩就需要捕获和跟踪接收到的所有信号,包括真实信号和欺骗信号,这就需要改变接收机捕获跟踪策略,增加信号跟踪通道数,运算量和硬件复杂度较高。
信号解扩后的欺骗干扰检测方法主要包括:信号功率检测、到达时间检测、到达角度检测、电文校验、定位结果检测等。对于阵列接收机,主要利用信号到达角一致性进行欺骗干扰检测,根据信号相位差实现欺骗干扰检测和判决,主要包括相位单差检测和相位双差检测。其一,通过比较信号到达两个天线的相位差与根据星历计算得到的估计值,进行欺骗干扰判别;其二,基于欺骗信号通常从同一天线发射的特点,不同卫星欺骗信号到达两个天线的相位差相同,故可通过不同卫星信号到达两固定天线的相位差的差值进行欺骗干扰的判别。但要得到信号相位及其差值,需要各个通道分别捕获和跟踪各颗卫星信号,运算量大,耗费资源多。因此,研究在保证干扰检测性能的条件下,检测系统复杂度低,检测计算量小的欺骗干扰检测方法对GNSS阵列接收机的优化设计及欺骗干扰检测性能提升都具有重要意义。
发明内容
本发明为了有效解决GNSS阵列接收机欺骗干扰检测问题,提供了一种信号解扩前的欺骗干扰检测方法。基本思路是:对于阵元间距为半波长的GNSS阵列接收机,导航信号在天线阵元间时延差很小,具有较强的相关性,而不同通道内噪声不相关。因此,利用不同通道的信号与噪声相关性的差异,消除噪声的影响,进行导航信号的功率估计,进而通过功率大小实现欺骗干扰的检测。具体技术方案如下:
一种GNSS阵列接收机信号解扩前的欺骗干扰检测方法,包括以下步骤:
(S1)信号功率估计,
以N元天线阵接收机为例,N为整数,接收阵列信号基带采样结果表示为r(nTs),r(nTs)=[r1(nTs),r2(nTs),...,rN(nTs)]T,ri(nTs)表示第i个阵元接收信号在nTs时刻的采样结果,n表示采样点,Ts为采样间隔;
以第1个阵元对应通道为参考通道,将天线阵第i个阵元,i=2,3,...,N,基带采样信号ri(nTs)与参考通道基带采样信号r1(nTs)的共轭值r1 *(nTs)进行相乘,并累加,得到第i个阵元接收信号的功率估计结果
K为累加点数,“*”表示共轭运算符,进一步将表示为实部和虚部的和:
其中分别为的实部和虚部,j表示虚数单位;
(S2)估计不存在欺骗干扰时的概率分布;
根据卫星分布特性对相位差、可见卫星数和信号功率进行取值,生成仿真数据,按照步骤S1描述方法进行蒙特卡洛仿真,得到不同相位差、可见卫星数和信号功率条件下的值,此处与步骤S1中含义相同,统计得到不存在欺骗干扰条件下的概率分布曲线曲线横坐标为功率,纵坐标为估计结果等于该功率值的概率;
(S3)给定虚警概率,根据所述步骤(S2)中概率分布曲线计算检测门限Th;
设置二元假设为,
H0:不存在欺骗干扰
H1:存在欺骗干扰
虚警概率Pfa定义为不存在欺骗干扰条件下大于门限Th的概率,在给定虚警概率Pfa的情况下,根据概率分布曲线直接读出对应的检测门限Th。
(S4)根据步骤(S1)得到的信号功率估计结果与检测门限Th进行比较,实现欺骗干扰的检测。
采用本发明获得的有益效果是:直接利用不同天线信号和噪声相关性差异,在解扩前估计信号功率,避免高计算量的捕获和跟踪处理,再根据普遍的真实信号可见性,功率分布情况,及检测虚警概率,计算干扰判决门限,进而实现欺骗干扰的检测。本发明检测性能良好且实现简单,不需要捕获和跟踪接收到的卫星信号,运算量小且实施方便,可直接用于GNSS阵列接收机基带采样信号的欺骗干扰检测。
附图说明
图1为本发明所提供的GNSS阵列接收机欺骗干扰检测方法的流程图;
图2为仿真实验中计算的功率估计实部的概率分布示意图;
图3为仿真实验中计算的功率估计虚部的概率分布示意图;
图4为本发明所提供的GNSS阵列接收机欺骗干扰检测方法欺骗干扰检测的ROC曲线。
具体实施方式
下面结合附图和实施例对本发明进行详细说明。
图1是本发明所提供的GNSS阵列接收机欺骗干扰检测方法的流程图。
第一步是阵列信号功率估计,通过对阵列基带采样信号r(nTs)不同阵元信号与参考阵元信号共轭相乘累加得到。
第二步是估计不存在欺骗干扰时的概率分布;
因为存在欺骗干扰和不存在欺骗干扰的情况下,具有不同的概率分布,以此可以进行欺骗干扰的检测。通常接收机在信号解扩前难以获得卫星分布情况和信号功率,且自身位置和方位未知,故通过一般情况下的卫星分布特性和信号的功率分布规律估计无欺骗干扰条件下的功率估计结果的分布情况。可以认为,随机放置的天线阵,信号在两个天线之间的相位差服从-π~π的均匀分布U(-π,π)。以全球定位系统(Global Positioning System,GPS)为例,通常可见卫星数为4~11颗,卫星导航信号到达地面附近的功率基本一致,对应载噪比(Carrier to Noise Ratio,CNR)变化范围为40~50dBHz,认为信号功率在变化范围内均匀分布,可见卫星数也服从4~11的均匀分布。根据上述空间分布特性对相位差、可见卫星数和信号功率大小进行取值,生成仿真数据r(nTs),按第一步描述方法进行蒙特卡洛仿真,统计得到不存在欺骗干扰条件下概率分布曲线的概率分布曲线曲线横坐标为功率,纵坐标为估计结果等于该功率值的概率。
第三步,根据给定虚警概率Pfa计算检测门限Th。
设置二元假设为:
H0:不存在欺骗干扰
H1:存在欺骗干扰
虚警概率Pfa定义为不存在欺骗干扰条件下大于门限Th的概率。在给定虚警概率Pfa的情况下,根据概率分布曲线直接读出对应的检测门限Th。
第四步为欺骗干扰检测,根据信号功率估计结果与检测门限Th进行比较,判决是否存在欺骗干扰。
图2至图3是仿真实验中计算的有无欺骗干扰条件下实部和虚部的概率分布曲线。实施例仿真中认为真实信号个数与欺骗信号个数相同,在4~11内均匀分布,真实信号和欺骗信号在两个天线之间的相位差服从均匀分布U(-π,π),真实信号CNR服从40~50dBHz内的均匀分布,欺骗干扰信号功率比真实信号功率高3dB,图中为按第二步所述,进行106次数值计算的统计结果,并对信号最小功率进行了归一化。
图2中实线为有欺骗干扰条件下实部的概率分布曲线,虚线为无欺骗干扰条件下实部的概率分布曲线,横坐标为有、无欺骗干扰条件下实部的取值范围,纵坐标为对应的概率。
图3中实线为有欺骗干扰条件下实部的概率分布曲线,虚线为无欺骗干扰条件下实部的概率分布曲线。横坐标为有、无欺骗干扰条件下虚部的取值范围,纵坐标为对应的概率。
图4为欺骗干扰检测的受试者工作特征(Receiver Operating Characteristic,ROC)曲线,用于衡量算法的干扰检测性能。横坐标为检测的虚警概率,纵坐标为对应的检测概率。从ROC曲线可以看出,3%的虚警概率即可实现约95%的检测概率,说明本发明提供的GNSS欺骗干扰检测方法性能优良。
以上实施例仅用于帮助理解本发明的方法,应当指出,对于技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以对本发明进行若干修饰和改进,这些修饰和改进也落入本发明权利要求的保护范围内。

Claims (3)

1.一种GNSS阵列接收机信号解扩前的欺骗干扰检测方法,其特征在于包括以下步骤:
(S1)信号功率估计,
以N元天线阵接收机为例,接收阵列信号基带采样结果表示为r(nTs),r(nTs)=[r1(nTs),r2(nTs),...,rN(nTs)]T,ri(nTs)表示第i个阵元接收信号在nTs时刻的采样结果,n表示采样点,Ts为采样间隔;
以第1个阵元对应通道为参考通道,将天线阵第i个阵元,i=2,3,...,N,基带采样信号ri(nTs)与参考通道基带采样信号r1(nTs)的共轭值r1 *(nTs)进行相乘,并累加,得到第i个阵元接收信号的功率估计结果
K为累加点数,进一步将表示为实部和虚部的和:
其中分别为的实部和虚部,j表示虚数单位;
(S2)估计不存在欺骗干扰时的概率分布;
(S3)给定虚警概率,根据所述步骤(S2)中的概率分布,计算检测门限Th;
(S4)根据步骤(S1)得到的信号功率估计结果与检测门限Th进行比较,实现欺骗干扰的检测。
2.如权利要求1所述的一种GNSS阵列接收机信号解扩前的欺骗干扰检测方法,其特征在于,所述步骤(S2)的具体过程为:根据卫星分布特性对相位差、可见卫星数和信号功率进行取值,生成仿真数据,按步骤(S1)的方法进行蒙特卡洛仿真,统计得到不存在欺骗干扰条件下的概率分布曲线曲线横坐标为功率,纵坐标为估计结果等于该功率值的概率。
3.如权利要求1所述的一种GNSS阵列接收机信号解扩前的欺骗干扰检测方法,其特征在于,所述步骤(S3)的具体过程为:
设置二元假设为,
H0:不存在欺骗干扰
H1:存在欺骗干扰
虚警概率Pfa定义为不存在欺骗干扰条件下大于门限Th的概率,在给定虚警概率Pfa的情况下,根据概率分布曲线直接读出对应的检测门限Th。
CN201810259165.9A 2018-03-27 2018-03-27 一种gnss阵列接收机信号解扩前的欺骗干扰检测方法 Active CN108535745B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810259165.9A CN108535745B (zh) 2018-03-27 2018-03-27 一种gnss阵列接收机信号解扩前的欺骗干扰检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810259165.9A CN108535745B (zh) 2018-03-27 2018-03-27 一种gnss阵列接收机信号解扩前的欺骗干扰检测方法

Publications (2)

Publication Number Publication Date
CN108535745A true CN108535745A (zh) 2018-09-14
CN108535745B CN108535745B (zh) 2020-08-14

Family

ID=63485181

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810259165.9A Active CN108535745B (zh) 2018-03-27 2018-03-27 一种gnss阵列接收机信号解扩前的欺骗干扰检测方法

Country Status (1)

Country Link
CN (1) CN108535745B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110471091A (zh) * 2019-08-29 2019-11-19 北京航空航天大学合肥创新研究院 一种基于相关器正交分量的欺骗干扰检测方法
CN112954586A (zh) * 2021-01-29 2021-06-11 北京邮电大学 一种欺骗干扰源定位方法、电子设备及存储介质
CN117452443A (zh) * 2023-12-22 2024-01-26 中国人民解放军国防科技大学 一种基于双门限判决的卫星导航欺骗干扰检测方法及装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140247186A1 (en) * 2012-10-26 2014-09-04 Uti Limited Partnership Global Navigation Satellite System (GNSS) Anti-Interference using Array Processing
CN104536014A (zh) * 2015-01-23 2015-04-22 哈尔滨工业大学 一种基于Relax算法的抗转发式欺骗干扰方法
EP3026462A1 (en) * 2014-11-21 2016-06-01 Akademia Marynarki Wojennej The method of removing interference, and circuit for removing interference in the gnss system
CN105656594A (zh) * 2016-01-12 2016-06-08 杭州电子科技大学 基于信道差异的转发式gnss欺骗干扰检测方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140247186A1 (en) * 2012-10-26 2014-09-04 Uti Limited Partnership Global Navigation Satellite System (GNSS) Anti-Interference using Array Processing
EP3026462A1 (en) * 2014-11-21 2016-06-01 Akademia Marynarki Wojennej The method of removing interference, and circuit for removing interference in the gnss system
CN104536014A (zh) * 2015-01-23 2015-04-22 哈尔滨工业大学 一种基于Relax算法的抗转发式欺骗干扰方法
CN105656594A (zh) * 2016-01-12 2016-06-08 杭州电子科技大学 基于信道差异的转发式gnss欺骗干扰检测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
FEIQIANG CHEN ET AL.: ""Combined algorithm for interference suppression and signal acquisition in GNSS receivers"", 《ELECTRONICS LETTERS》 *
王芝应 等: ""BOC接收机捕获阶段转发式欺骗信号检测算法"", 《全球定位系统》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110471091A (zh) * 2019-08-29 2019-11-19 北京航空航天大学合肥创新研究院 一种基于相关器正交分量的欺骗干扰检测方法
CN110471091B (zh) * 2019-08-29 2021-05-25 北京航空航天大学合肥创新研究院 一种基于相关器正交分量的欺骗干扰检测方法
CN112954586A (zh) * 2021-01-29 2021-06-11 北京邮电大学 一种欺骗干扰源定位方法、电子设备及存储介质
CN117452443A (zh) * 2023-12-22 2024-01-26 中国人民解放军国防科技大学 一种基于双门限判决的卫星导航欺骗干扰检测方法及装置
CN117452443B (zh) * 2023-12-22 2024-03-22 中国人民解放军国防科技大学 一种基于双门限判决的卫星导航欺骗干扰检测方法及装置

Also Published As

Publication number Publication date
CN108535745B (zh) 2020-08-14

Similar Documents

Publication Publication Date Title
Nielsen et al. GNSS spoofing detection for single antenna handheld receivers
Liu et al. On the performance of the cross-correlation detector for passive radar applications
US11714199B2 (en) System and method for detecting spoofing of GNSS signals
CN106125053B (zh) 脉冲多普勒雷达极化抗干扰方法
CN104502899A (zh) 一种自适应的恒虚警率目标检测方法
CN104267379A (zh) 一种基于波形设计的主被动雷达协同抗干扰方法
CN104020451A (zh) 基于聚类的外辐射源雷达目标航迹处理方法
CN103217670B (zh) 一种基于pca的外辐射源微弱信号检测方法
CN103954951A (zh) 基于功率中值和归一化协方差估计的海面目标检测方法
CN104656104A (zh) 基于最大似然估计的卫星导航欺骗信号识别方法及系统
CN108535745A (zh) 一种gnss阵列接收机信号解扩前的欺骗干扰检测方法
CN102944873A (zh) 基于多频点回波幅度逆序统计量的低空目标检测方法
Dehghanian et al. GNSS spoofing detection based on signal power measurements: statistical analysis
CN104977585B (zh) 一种鲁棒的运动声纳目标检测方法
CN110879404A (zh) 一种基于相关峰和残留信号相结合的gnss欺骗干扰检测方法
CN114114192B (zh) 集群目标检测方法
CN103954943A (zh) 网络化雷达系统抗欺骗式干扰方法
CN106249212A (zh) 主瓣压制干扰背景下有源假目标的极化鉴别方法
CN107272030B (zh) 一种基于多个gps卫星信息融合的目标探测方法
De Maio et al. Adaptive radar detection: A Bayesian approach
CN111896928A (zh) 基于有源欺骗式干扰消除的多站雷达目标检测方法
CN112596078A (zh) 一种基于载噪比统计的卫星导航欺骗干扰检测方法
Broumandan et al. A network-based GNSS structural interference detection, classification and source localization
CN111198387A (zh) 一种抗欺骗干扰的空时采样导航定位方法
CN102621535B (zh) 一种高效的协方差矩阵结构估计方法

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