CN108196270B - 一种双基线欺骗干扰检测算法的性能分析方法 - Google Patents

一种双基线欺骗干扰检测算法的性能分析方法 Download PDF

Info

Publication number
CN108196270B
CN108196270B CN201711494740.5A CN201711494740A CN108196270B CN 108196270 B CN108196270 B CN 108196270B CN 201711494740 A CN201711494740 A CN 201711494740A CN 108196270 B CN108196270 B CN 108196270B
Authority
CN
China
Prior art keywords
baseline
double
deception
function
performance
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
Application number
CN201711494740.5A
Other languages
English (en)
Other versions
CN108196270A (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.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and 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 Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201711494740.5A priority Critical patent/CN108196270B/zh
Publication of CN108196270A publication Critical patent/CN108196270A/zh
Application granted granted Critical
Publication of CN108196270B publication Critical patent/CN108196270B/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)
  • Radar Systems Or Details Thereof (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明公开了一种双基线欺骗干扰检测算法性能分析方法,首先将两条基线S1、S2长度假设为高斯分布;然后建立欺骗检测的判断函数,使欺骗与未欺骗结果不同;接着根据判断函数,通过假设检验,建立性能分析数学模型;最后通过模型得到算法检测概率,确定双基线欺骗干扰测试方法性能。本发明在设定漏检率Pfa的情况下,通过判断Pd大小,就可以判断双基线欺骗干扰检测算法的性能,分析方法简单,分析速度快,为不同环境下双基线欺骗干扰检测算法的参数设置提供了理论依据。

Description

一种双基线欺骗干扰检测算法的性能分析方法
技术领域
本发明涉及干扰欺骗方法,特别是涉及一种双基线欺骗干扰检测算法的性能分析方法。
背景技术
目前,卫星信号阻塞和欺骗阻碍了它在某些特定区域的应用。由于信号阻塞和欺骗情况的存在会导致信号出现丢失或错误的情况,这对于使用卫星定位的系统很不利。相对于信号阻塞,信号欺骗对卫星定位导航的影响更大。
目前已经存在了各种各样的欺骗检测算法,从简单的观测卫星信号功率大小到寻找相关器输出峰值的位置,再到利用非无线信号源辅助进行干扰检测(例如惯导测量单元)和利用双基线进行欺骗干扰检测。双基线欺骗干扰算法与传统欺骗检测方法,检测过程容易,检测结果更可靠。然而影响双基线欺骗干扰检测算法性能因素有基线长度值、相关系数,不同的参数设置,适用于不同的情况。为了得到最佳的检测效果,有必要对双基线欺骗干扰检测算法进行性能分析,选择最适合的参数设置。
发明内容
本发明的目的是提供一种基于双基线长度的欺骗干扰检测算法的性能分析方法,为后续提高此欺骗干扰信号检测算法做铺垫,以提高导航定位精度。
实现本发明目的的技术解决方案为:一种双基线欺骗干扰检测算法的性能分析方法,包括如下步骤:
步骤1、条件假设:将两条基线S1、S2长度假设为高斯分布;
步骤2:判断函数:建立欺骗检测的判断函数,使欺骗与未欺骗结果不同;
步骤3、性能分析模型:根据判断函数,通过假设检验,建立性能分析数学模型;
步骤4、分析性能:通过模型得到算法检测概率,确定双基线欺骗干扰检测算法性能分析方法。
本发明与现有技术相比,其显著优点为:本发明在设定漏检率Pfa的情况下,通过判断Pd大小,就可以判断双基线欺骗干扰检测算法的性能,分析方法简单,分析速度快,为不同环境下双基线欺骗干扰检测算法的参数设置提供了理论依据。
附图说明
图1是四个接收机的放置方式图。
图2是双基线欺骗干扰检测算法的性能分析方法流程图。
图3是d=2σ,d=4σ,d=6σ三种情况下判断函数1的分析结果图。
图4是ρ=-0.6,ρ=-0.4,ρ=0,ρ=0.4,ρ=0.6五种情况下判断函数1的分析结果图。
图5是d=2σ,d=4σ,d=6σ三种情况下判断函数2的分析结果图。
图6是d=4σ时两种判断函数的分析结果对比图。
图7是d=6σ时两种判断函数的分析结果对比图。
具体实施方案
下面结合附图和具体实施例,进一步说明本发明方案。
一种双基线欺骗干扰检测算法的性能分析方法,包括如下步骤:
步骤1、条件假设:假设双基线S1,S2均属于高斯分布,其具体表示如下:
Figure GDA0003163126620000021
其中,dS1,dS2分别为双基线的数学期望,
Figure GDA0003163126620000022
为基线的方差,N(·)表示高斯分布。
步骤2、判断函数:建立欺骗检测中需要使用到的判断函数,只要得到的欺骗与未欺骗结果不同即可,实验中设置的判断函数具体如下:
(1)判断函数1设置为dS1+S2=S1+S2;由于基线S1,S2为高斯分布,可以得到判断函数
Figure GDA0003163126620000023
其中,
Figure GDA0003163126620000024
为判断函数的方差。ρ为相关系数。
(2)判断函数2设置为最大最小模型函数dmax min,当未被欺骗时,判断函数设为dmin;当受到欺骗时,判断函数设为dmax。由于基线S1,S2为高斯分布,可以得到判断函数dmax min的分布。
步骤3、性能分析模型:根据步骤2中所述的两种判断函数,分别通过假设检验建立其数学模型,具体如下:
(1)当判断函数设为dS1+S2=S1+S2时,模型建立步骤如下:
1)假设无欺骗干扰为假设检验的H0条件,表示如下:
Figure GDA0003163126620000031
其中,dns代表无欺骗情况下的判断函数。
由(2)式可以拓展出此情况下的漏检概率Pfa为:
Figure GDA0003163126620000032
根据高斯分布的特点,通过公式(4)可以计算出阈值α
Figure GDA0003163126620000033
Figure GDA0003163126620000034
2)假设存在欺骗为假设检验中的H1条件,表示如下:
Figure GDA0003163126620000035
其中,ds代表被欺骗情况下的判断函数。
由(6)式可以拓展出此情况下的检测概率Pd为:
Figure GDA0003163126620000036
将由式(5)计算得到的阈值α代入式(7),得到检测概率公式(8):
Figure GDA0003163126620000037
假设nσ=nσS1=nσS2=dS1=dS2,则公式(8)可以简化如下:
Figure GDA0003163126620000038
(2)当判断函数设为最大最小模型时,模型建立步骤如下:
1)同样假设无欺骗为假设检验中的H0情况,表示如下:
H0:min(S1,S2)>α
Figure GDA0003163126620000041
由此可以计算出其漏检概率Pfa为:
Pfa=P(min(S1,S2)<α|H0)=1-(1-P(S1<α)·P(S2<α)) (11)
Figure GDA0003163126620000042
Figure GDA0003163126620000043
2)同样假设被欺骗为假设检验中的H1情况,表示如下:
Figure GDA0003163126620000044
由上式可以得出检测概率的公式为:
Pd=P(max(S1,S2)<α|H1)=(P(S1<α))·(P(S2<α)) (15)
Figure GDA0003163126620000045
假设nσS1=nσS2=nσ=dS1=dS2,代入式(13)可得:
Figure GDA0003163126620000046
由上式可以推导出阈值α:
Figure GDA0003163126620000047
将式(18)代入式(16),可以得到:
Figure GDA0003163126620000048
步骤4、分析性能:通过模型得到算法检测概率,确定双基线欺骗干扰检测算法的性能分析方法,在漏检率Pfa一样的情况下,检测率Pd值越大,说明干扰检测算法性能越好。
基于本发明设计的双基线欺骗干扰检测算法的性能分析方法,下面分别显示当基线长度和相关系数不同时的性能分析结果的比较。具体如下:
(1)在基线长度d=2σ,d=4σ,d=6σ条件下,判断函数为1时的性能分析图如图3所示。图3显示出不同基线长度下,算法的检测性能对比。
(2)在相关系数ρ=-0.6,ρ=-0.4,ρ=0,ρ=0.4,ρ=0.6条件下,判断函数1时的性能分析图如图4所示。图4显示出不同相关系数下,算法的检测性能对比。
(3)在基线长度d=2σ,d=4σ,d=6σ条件下,判断函数2时的性能分析图如图5所示。图5显示出不同基线长度下,算法的检测性能对比。
(4)图6和7显示出当基线长度d=4σ,d=6σ,相关系数ρ=0的条件下,选用判断函数1和判断函数2性能分析结果图,可以看出,判断的结果大致相同,说明本发明方法分析方法有效。

Claims (3)

1.一种双基线欺骗干扰检测算法的性能分析方法,其特征在于,包括如下步骤:
步骤1、条件假设:将两条基线S1、S2长度假设为高斯分布;
步骤2:判断函数:建立欺骗检测的判断函数,使欺骗与未欺骗结果不同;
步骤3、性能分析模型:根据判断函数,通过假设检验,建立性能分析数学模型;
步骤4、分析性能:通过模型得到算法检测概率,确定双基线欺骗干扰测试方法性能;
步骤1中基线S1、S2长度具体表示如下:
Figure FDA0003163126610000011
其中,dS1,dS2分别为双基线的数学期望,
Figure FDA0003163126610000012
为基线的方差,N(·)表示高斯分布。
步骤2的判断函数为:
Figure FDA0003163126610000013
其中,
Figure FDA0003163126610000014
为判断函数的方差,ρ为相关系数,dS1,dS2分别为双基线的数学期望,
Figure FDA0003163126610000015
为基线的方差,N(·)表示高斯分布。
则步骤3建立性能分析模型步骤为:
步骤3.1、假设无欺骗干扰为假设检验的H0条件,表示如下:
Figure FDA0003163126610000016
其中,dns代表无欺骗情况下的判断函数;
步骤3.2、设置漏检概率Pfa,计算出阈值α:
Figure FDA0003163126610000017
步骤3.3、假设存在欺骗为假设检验中的H1条件,表示如下:
Figure FDA0003163126610000018
其中,ds代表被欺骗情况下的判断函数;
步骤3.4、根据阈值α,计算检测概率Pd
Figure FDA0003163126610000021
2.根据权利要求1所述的双基线欺骗干扰检测算法的性能分析方法,其特征在于,步骤2的判断函数为最大最小模型函数dmaxmin,当未被欺骗时,判断函数设为dmin;当受到欺骗时,判断函数设为dmax
3.根据权利要求2所述的双基线欺骗干扰检测算法的性能分析方法,其特征在于,步骤3建立性能分析模型步骤为:
步骤3.1、假设无欺骗为假设检验中的H0情况,表示如下:
Figure FDA0003163126610000022
步骤3.2、设置漏检概率Pfa,计算出阈值α:
Figure FDA0003163126610000023
步骤3.3、假设存在欺骗为假设检验中的H1条件,表示如下:
Figure FDA0003163126610000024
步骤3.4、根据阈值α,计算检测概率Pd
Figure FDA0003163126610000025
CN201711494740.5A 2017-12-31 2017-12-31 一种双基线欺骗干扰检测算法的性能分析方法 Active CN108196270B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711494740.5A CN108196270B (zh) 2017-12-31 2017-12-31 一种双基线欺骗干扰检测算法的性能分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711494740.5A CN108196270B (zh) 2017-12-31 2017-12-31 一种双基线欺骗干扰检测算法的性能分析方法

Publications (2)

Publication Number Publication Date
CN108196270A CN108196270A (zh) 2018-06-22
CN108196270B true CN108196270B (zh) 2021-09-21

Family

ID=62587513

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711494740.5A Active CN108196270B (zh) 2017-12-31 2017-12-31 一种双基线欺骗干扰检测算法的性能分析方法

Country Status (1)

Country Link
CN (1) CN108196270B (zh)

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7219032B2 (en) * 2002-04-20 2007-05-15 John Louis Spiesberger Estimation algorithms and location techniques
KR20070016394A (ko) * 2005-08-03 2007-02-08 삼성탈레스 주식회사 표적 추적에서의 탐지 영역 데이터를 이용한 기만 재밍탐지 방법
CN101349747A (zh) * 2007-07-20 2009-01-21 汤子跃 一种雷达抗有源欺骗干扰性能定量评估方法
US9366762B2 (en) * 2011-10-27 2016-06-14 Southwest Research Institute Detection of deceptive navigation signal
CN102538793B (zh) * 2011-12-23 2015-03-18 北京控制工程研究所 一种双基线非合作目标双目测量系统
JP2013251808A (ja) * 2012-06-01 2013-12-12 Lighthouse Technology & Consulting Co Ltd 衛星測位信号生成装置及び方法、生成プログラム、並びに、衛星測位信号受信装置及び方法
CN103454650B (zh) * 2013-08-20 2015-06-24 北京航空航天大学 一种视觉辅助卫星完好性监测方法
CN104181552B (zh) * 2014-08-21 2017-07-25 武汉大学 一种动态gnss接收机的抗干扰正态零陷加宽的方法
CN104793220A (zh) * 2015-05-04 2015-07-22 中国电子科技集团公司第五十四研究所 一种基于多天线的欺骗干扰检测方法
CN105656594B (zh) * 2016-01-12 2019-04-30 杭州电子科技大学 基于信道差异的转发式gnss欺骗干扰检测方法
CN105607092B (zh) * 2016-01-27 2017-10-31 中国人民解放军国防科学技术大学 基于tdoa和功率测量值的gnss欺骗干扰定位方法

Also Published As

Publication number Publication date
CN108196270A (zh) 2018-06-22

Similar Documents

Publication Publication Date Title
JP5985372B2 (ja) 目標検出装置及び目標検出方法
KR101271430B1 (ko) 무선 통신 시스템에서의 수신 신호 검출 방법 및 장치
KR20150131779A (ko) 순차통계 일정 오경보율 검파의 처리 속도 향상방법
US11303485B2 (en) Signal acquisition method and device
JP2015507866A (ja) 無線信号のなりすましを検出するシステムおよび方法
TW200925632A (en) GNSS satellite signal handling method and correlator implementing the same
CN108881084B (zh) 一种基于gp分布的bpsk/qpsk信号识别方法
CN104253659B (zh) 一种频谱检测方法及其装置
CN113219504B (zh) 定位信息确定方法及装置
CN109919957B (zh) 一种基于动态视觉传感器的角点检测方法
CN108196270B (zh) 一种双基线欺骗干扰检测算法的性能分析方法
CN105656511B (zh) 一种适应于有频偏和低信噪比环境下的差分相关捕获方法
WO2021196765A1 (zh) 到达时间确定方法、装置、终端设备及存储介质
JP2010117332A (ja) バイスタティックレーダ装置
CN108718223B (zh) 一种非合作信号的盲频谱感知方法
CN105553619B (zh) 基于调制识别辅助的突发信号检测方法
EP3255449A1 (en) Acquisition method and device of time of arrival for positioning mobile terminal
CN105959035B (zh) 一种直扩信号截获检测方法
CN1904639A (zh) 一种采用分形特征值的信号检测方法
CN110110341B (zh) 一种判决门限自动推荐的脉冲检测方法
RU2572219C1 (ru) Способ обработки сигнала шумоизлучения объекта
US20200257878A1 (en) Systems and methods for distance independent differential signature detection
RU2634786C1 (ru) Способ определения маневра шумящего объекта
CN110939875B (zh) 基于谱纹分析的管道泄露检测与定位方法及其系统
RU2570430C1 (ru) Способ классификации шумящих объектов

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