CN112782728B - 一种基于惯性辅助的天线阵欺骗干扰信号检测方法 - Google Patents
一种基于惯性辅助的天线阵欺骗干扰信号检测方法 Download PDFInfo
- Publication number
- CN112782728B CN112782728B CN202110106507.5A CN202110106507A CN112782728B CN 112782728 B CN112782728 B CN 112782728B CN 202110106507 A CN202110106507 A CN 202110106507A CN 112782728 B CN112782728 B CN 112782728B
- Authority
- CN
- China
- Prior art keywords
- satellite
- antenna array
- signal
- detected
- signals
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 65
- 238000012360 testing method Methods 0.000 claims abstract description 28
- 238000000034 method Methods 0.000 claims abstract description 20
- 238000004364 calculation method Methods 0.000 claims description 14
- 239000013598 vector Substances 0.000 claims description 11
- 239000011159 matrix material Substances 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 5
- 230000002452 interceptive effect Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000010276 construction Methods 0.000 claims 1
- 238000005516 engineering process Methods 0.000 abstract description 9
- 238000012986 modification Methods 0.000 description 5
- 230000004048 modification Effects 0.000 description 5
- 230000001629 suppression Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 208000003443 Unconsciousness Diseases 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000009466 transformation Effects 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/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
-
- 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/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
-
- 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)
- Signal Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供一种基于惯性辅助的天线阵欺骗干扰信号检测方法,由惯性辅助解算每一个虚拟真实卫星信号在三元天线阵坐标系中的到达角;计算虚拟真实卫星信号和待检测卫星信号的到达角在三元天线阵不同基线上的载波相位差;基于采用载波相位双差算法构建的欺骗干扰信号的二元假设检验模型及其判决依据,判决每一个待检测卫星信号是否为欺骗干扰信号。充分利用卫星欺骗干扰信号与惯性辅助解算的虚拟真实卫星信号在三元天线坐标系中的到达角差异,通过计算待检测卫星信号与虚拟真实卫星信号在三元天线阵基线上载波相位差的不同来识别欺骗干扰,解决了传统载波相位双差检测技术检测时间较长且对单路卫星欺骗干扰信号检测失效的问题。
Description
技术领域
本发明涉及卫星导航接收***防欺骗技术领域,更具体地,涉及一种基于惯性辅助的天线阵欺骗干扰信号检测方法。
背景技术
近年来,针对卫星导航信号的压制与反压制,欺骗与反欺骗技术得到了充分的研究与发展。卫星导航接收***在恶劣电磁干扰环境中受到的人为干扰主要包括压制式干扰和欺骗式干扰。压制式干扰利用大功率信号对目标用户进行压制,通常会造成卫星接收机环路失锁和定位失败,在战场环境中的隐蔽性较差。欺骗式干扰的功率水平、信号格式以及频谱结构与真实卫星信号相似,旨在迫使卫星导航接收机在毫无意识的情况下误锁在欺骗信号上,产生看似可靠、实则错误的导航定位结果,从而达到控制卫星导航接收机的目的。相比于压制式干扰,欺骗式干扰不容易被对方察觉,危害性更大。因此,有必要针对卫星导航欺骗干扰信号的检测方法展开深入研究。
利用卫星导航接收机的阵列天线对接收到的卫星信号进行到达角检测是当前较为有效的欺骗干扰检测方法之一。常用的到达角检测技术主要包括载波相位单差检测技术和载波相位双差检测技术。其中,载波相位单差检测技术不能直接应用于运动载体;载波相位双差检测技术主要通过遍历接收到的任意两路卫星信号载波相位差是否相等来检测欺骗干扰,虽然解决了运动载体欺骗干扰检测问题,但检测时间较长且对单路卫星欺骗干扰信号检测失效。
发明内容
本发明提供一种克服上述问题或者至少部分地解决上述问题的一种基于惯性辅助的天线阵欺骗干扰信号检测方法,包括:
由惯性辅助解算每一个虚拟真实卫星信号在天线阵坐标系中的到达角;
计算惯性辅助解算的每一个虚拟真实卫星信号的到达角在天线阵不同基线上的载波相位差,以及每一个待检测卫星信号的到达角在天线阵不同基线上的载波相位差;
采用载波相位双差算法构建欺骗干扰信号的二元假设检验模型,并基于二元假设检验模型设定欺骗干扰信号的判决依据;
遍历惯性辅助解算的每一个虚拟真实卫星信号与每一个待检测卫星信号之间的载波相位双差,基于二元假设检验模型及其判决依据,判决每一个待检测卫星信号是否为欺骗干扰信号。
本发明提供的一种基于惯性辅助的天线阵欺骗干扰信号检测方法,充分利用卫星欺骗干扰信号与惯性辅助解算的虚拟真实卫星信号在天线坐标系中的到达角差异,通过计算待检测卫星信号与虚拟真实卫星信号在天线阵基线上载波相位差的不同来识别欺骗干扰,解决了传统载波相位双差检测技术检测时间较长且对单路卫星欺骗干扰信号检测失效的问题。
附图说明
图1为本发明提供的一种基于惯性辅助的天线阵欺骗干扰信号检测方法流程图;
图2为天上各颗卫星与卫星导航接收机天线阵的位置关系示意图;
图3为卫星导航接收机天线阵坐标系示意图;
图4为本发明提供的一种基于惯性辅助的天线阵欺骗干扰信号检测方法的整个流程图。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
针对单天线转发的卫星欺骗干扰场景,传统的对卫星欺骗干扰信号采用载波相位双差算法进行检测的过程为,对于卫星导航接收机天线阵接收到的所有的卫星信号,包括真实卫星信号和欺骗干扰信号,此时不知道哪些是真实卫星信号,哪些是欺骗干扰信号。
由于每一颗真实卫星信号的来向(可以理解为到达角)是不同的,而所有的欺骗干扰信号的来向是相同的。因此,需要将接收到的所有的卫星信号的任意两路的来向进行对比。此时,如果只有一路欺骗干扰信号,则该路欺骗干扰信号的来向与其它所有的真实卫星信号的来向均不同,因此,没有办法识别出该路欺骗干扰信号;只有当存在两路或两路以上的欺骗干扰信号时,这多路欺骗干扰信号的来向与每一路真实卫星信号的来向均不同,但是这多路欺骗干扰信号的来向均相同,因此,判定这多路卫星信号为多路欺骗干扰信号。因此,传统的欺骗干扰信号检测方法不能够检测单路欺骗干扰信号。
需要说明的是,每一路卫星信号在天线阵坐标中的到达角是无法直接测量到的,因此,对任意两路卫星信号的来向(到达角)的比较需要转化为载波相位差的比较。
基于传统欺骗干扰信号检测方法无法对单路欺骗干扰信号检测以及检测时间长的缺陷,本发明提供了一种基于惯性辅助的欺骗干扰信号检测方法,图1为本发明提供的一种基于惯性辅助的天线阵欺骗干扰信号检测方法流程图,如图1所示,方法包括:101、由惯性辅助解算每一个虚拟真实卫星信号在天线阵坐标系中的到达角;102、计算惯性辅助解算的每一个虚拟真实卫星信号的到达角在天线阵不同基线上的载波相位差,以及每一个待检测卫星信号的到达角在天线阵不同基线上的载波相位差;103、采用载波相位双差算法构建欺骗干扰信号的二元假设检验模型,并基于二元假设检验模型设定欺骗干扰信号的判决依据;104、遍历惯性辅助解算的每一个虚拟真实卫星信号与每一个待检测卫星信号之间的载波相位双差,基于二元假设检验模型及其判决依据,判决每一个待检测卫星信号是否为欺骗干扰信号。
可以理解的是,在本发明中,由惯性辅助解算出每一个虚拟真实卫星信号在卫星导航接收机天线阵坐标系中的到达角,该到达角为真实卫星信号的理论值,可作为后续识别卫星欺骗干扰信号的参考。然后,计算惯性辅助解算的每一个虚拟真实卫星信号在天线阵不同基线上的载波相位差,以及通过天线阵测量得到每一个待检测卫星信号在天线阵不同基线上的载波相位差。由于每一个虚拟真实卫星信号在天线阵坐标系中的到达角是由真实卫星在导航坐标系中的俯仰角、方位角通过姿态转换得到,因此,可以通过比较每一个待检测卫星信号在天线阵不同基线上的载波相位差和每一个虚拟真实卫星信号在天线阵不同基线上的载波相位差是否相等,来识别每一个待检测卫星信号是否为卫星欺骗干扰信号。
基于该原理,采用载波相位双差算法构建欺骗干扰信号的二元假设检验模型,并设定二元假设检验模型的判决依据。根据构建的二元假设检验模块及其判决依据,根据每一个虚拟真实卫星信号在天线阵不同基线上的载波相位差和每一个待检测卫星信号在天线阵不同基线上的载波相位差,利用二元假设检验模型及判决依据,对每一个待检测卫星信号进行欺骗干扰信号的识别。遍历所有的虚拟真实卫星信号和所有的待检测卫星信号的载波相位双差,在整个空域范围内检测卫星欺骗干扰信号。
本发明充分利用卫星欺骗干扰信号与惯性辅助解算的虚拟真实卫星信号在天线坐标系中的到达角差异,通过计算待检测卫星信号与虚拟真实卫星信号在天线阵基线上载波相位差的不同来识别欺骗干扰,由于通过惯性辅助解算能够得到每一个真实卫星信号的到达角和载波相位差,在识别欺骗干扰信号时,如果其中一路卫星信号的载波相位差与任何一路真实卫星信号的载波相位差不同时,那么该路卫星信号即为卫星欺骗干扰信号,对于单路欺骗干扰信号,也可以识别,解决了传统载波相位双差检测技术对单路卫星欺骗干扰信号检测失效的问题。
在一种可能的实施例方式中,由惯性辅助解算每一个虚拟真实卫星信号在天线阵坐标系中的到达角包括:利用惯性辅助解算卫星导航接收机天线阵相对于当地导航坐标系的姿态矩阵;通过卫星导航接收机信号处理得到各颗真实卫星在当地导航坐标系下的方位角和俯仰角;基于各颗真实卫星在当地导航坐标系下的方位角和俯仰角以及惯性辅助解算的卫星导航接收机天线阵相对于当地导航坐标系的姿态矩阵,通过坐标转换得到各颗真实卫星在天线阵坐标下的到达角,到达角包括方位角和俯仰角,由于该真实卫星信号的到达角主要由惯性辅助解算得到,并非通过天线阵信号处理检测所得,故称之为惯性辅助解算的虚拟真实卫星信号的到达角。
可以理解的是,天上各颗卫星与接收机天线阵的位置关系示意如图2所示,接收机天线阵的位置为Pu,天上卫星的位置为Psv,PuPsv的单位矢量就是从接收机天线阵到天上卫星的方向余弦矢量。根据真实卫星在当地导航坐标系中的俯仰角α和方位角β的定义可知,将天上第i颗卫星的方向余弦矢量DC,表示在东北天(ENU)导航坐标系中为:
DCENU=[-sinβcosα cosβcosα sinα]T。
利用惯导解算得到天线阵坐标系相对于东北天导航坐标系(当地导航坐标系)的坐标转换矩阵表示为:
其中,由惯导解算的接收机天线阵俯仰角为θ、方位角为横滚角为γ。由此,可得天上卫星在接收机天线阵坐标系中的方向矢量DCb为:
若DCb=[κe κn κu],则由惯性辅助解算的天上第i颗卫星在天线阵坐标系中的俯仰角θi和方位角分别为:
在一种可能的实施例方式中,在计算惯性辅助解算的虚拟真实卫星信号和待检测卫星信号在天线阵不同基线上的载波相位差时,通过研究表明,当卫星导航接收机天线阵元数大于等于三,且天线阵型不为线阵时,阵元间的载波相位差可与卫星入射信号到达角一一对应。为方便分析,以正三角形天线阵为例,设天线阵的基线长度为b,以三个阵元所在的平面为X轴和Y轴所在平面,以A为原点建立直角坐标系。AB基线方向为X轴方向,Z轴指向天顶方向,XYZ轴满足右手定则,如图3所示。
卫星入射信号在天线阵坐标系下的单位方向矢量可以表示为:
阵列天线三个阵元的坐标分别为A=[0,0,0],b=[b,0,0]和相应阵元间基线方向矢量分别为bAB、bAC和bBC。
假设由惯性辅助解算的虚拟真实卫星信号在天线阵坐标系中的到达角为Ψi,则虚拟真实卫星信号i与待检测卫星信号j在天线阵基线AB上的载波相位差可分别表示为:
其中,为待检测卫星信号载波相位测量值的测量噪声,Ψi为由惯性辅助解算的虚拟真实卫星信号在天线阵坐标系中的到达角,当采用高精度的惯导时,该到达角的计算误差相对于天线阵检测实际信号到达角的精度可忽略不计。同理,可得基线AC上的载波相位差如下:
由于AB、AC和BC这三条基线的载波相位差存在相关性,BC基线上的载波相位差可以直接由AB和AC基线上的载波相位差得到,因此BC基线对应的载波相位差对载波相位双差检测没有贡献,不需要单独给出。
由于三元阵为正三角形,所以有:
xB=b,
虚拟真实卫星信号i与待检测卫星信号j在天线阵基线上的载波相位差可以写成:
dΦi=HΛi;
dΦj=HΛj+γj;
其中,
在一种可能的实施例方式中,上述对虚拟真实卫星信号i与待检测卫星信号j在天线阵基线上的载波相位差进行了表示,其中,虚拟真实卫星信号i与待检测卫星信号j的载波相位双差矢量可表示为:
ΔdΦij=dΦi-dΦj=H(Λi-Λj)+(γi-γj)=HRij+ηij;
式中,γi为0,当待检测卫星信号为欺骗干扰信号时,Rij可以认为不为零;当待检测卫星信号为真实卫星信号时,Rij等于零。由此,可利用上式表示的载波相位双差进行欺骗干扰信号的检测。
利用载波相位双差检测的原理,构建欺骗干扰信号检测的二元假设检验模型:
其中,H0/H1的假设条件分别为:
H0:参与载波相位双差检测的待检测卫星信号为欺骗干扰信号;
H1:参与载波相位双差检测的待检测卫星信号为真实卫星信号。
在一种可能的实施例方式,上述构建了欺骗干扰信号检测的二元假设检验模型,本实施例基于构建的二元假设检验模型,设定欺骗干扰信号的判决依据,依据判决依据进行欺骗干扰信号的检测。
具体地,假设三阵元A、B和C输出同一路信号的载波相位测量值误差A、B和C相互独立,且都服从均值为0,方差为σ2的高斯分布。则γAB、γAC服从均值为0,方差为2σ2的高斯分布。那么,ηij可以表示为:
易知ηij服从N(0,C)分布,其中协方差矩阵C可以表示为:
相应的可以求得C-1为:
当待检测信号为卫星真实信号时,载波相位双差矢量的分布满足Rij=ηij~N(0,C);
可见,上式为二维高斯分布,难以直接得到不同假设条件下的概率分布。因此,构造上式高斯随机变量的二次型,可以得到服从自由度为2的分布的检测量TH1,即有:
当待检测信号为卫星欺骗信号时,载波相位双差矢量的分布满足:
Rij=ηij~N(μ,C);
其中,μ为天线阵接收卫星欺骗信号情况下Rij的均值。
此时,构造的统计量TH0服从自由度为2的非中心分布,即有:
其中,
λ=μTC-1μ。
当给定虚警率PFA时,欺骗干扰信号的检测门限可表示为:
其中,为自由度为中心χ2分布右尾概率的反函数。
令检测量表示为r=(Rij)TC-1(Rij),二元假设检测模型的判决依据为:
r>ρth,则判为H0;
r<ρth,则判为H1。
此时,对应的欺骗干扰信号检测概率可表示为:
其中,表示自由度为2、非中心参数为λ的非中心分布χ2右尾概率。
当惯性辅助解算的虚拟真实卫星信号与待检测卫星信号在天线阵坐标系中的到达角差异较小时,很难判别待检测的卫星信号具体为真实卫星信号还是欺骗干扰信号。因此,设定虚警率恒为0.01,欺骗干扰信号最小检测概率PD,min为0.99。将待检测信号欺骗干扰检测概率小于PD,min的空域范围定义为检测盲区。当待检测信号处于检测盲区时,将待检测信号视为真实的卫星信号,并借助卫星接收机的其他欺骗干扰检测手段做进一步处理。其中,可将较小的虚警率和较高的最小检测概率作为本发明方法的输入参数,缩小欺骗干扰信号的检测盲区,最终提升对欺骗干扰信号的检测概率。
在一种可能的实施例方式中,遍历惯性辅助解算的每一个虚拟真实卫星信号与每一个待检测卫星信号之间的载波相位双差,基于二元假设检验模型及其判决依据,判决每一个待检测卫星信号是否为欺骗干扰信号包括:基于任一个虚拟真实卫星信号与每一个待检测卫星信号在天线阵不同基线上的载波相位双差,利用二元假设检验模型及其判决依据,判定相对于任一个虚拟真实卫星信号,所有待检测卫星信号中的欺骗干扰信号;遍历所有的虚拟真实卫星信号与每一个待检测卫星信号,判定出所有待检测卫星信号中的欺骗干扰信号。
可以理解的是,针对每个惯性辅助解算的虚拟真实卫星信号,首先遍历所有的待检测卫星信号,进行欺骗干扰检测,然后直至空域中所有惯性辅助解算的虚拟真实卫星信号遍历完成。
受当前技术条件的限制,卫星信号欺骗干扰源大都采用单一天线发射欺骗干扰信号。针对单天线发射的多星欺骗干扰场景,当待检测信号位于所有检测盲区外,即可判定该待检测信号为欺骗信号。存在欺骗干扰信号位于空域中某一路参考信号周边的检测盲区内,并将其识别为真实卫星信号。此时,该卫星欺骗干扰信号可利用卫星接收机后续的欺骗干扰检测手段做进一步处理,同时不影响其他路的真实卫星信号参与卫导接收机的基带信号处理和导航定位解算,从而实现对单天线发射的多星欺骗干扰信号的检测和抑制。
参见图4,为本发明提供的基于惯性辅助的天线阵欺骗干扰信号检测方法的整体流程图,在对整个空域范围内的待检测卫星信号进行欺骗干扰检测时,对于多个真实卫星信号,比如,有k个真实卫星信号,有n个待检测卫星信号,其中,待检测卫星信号中包括真实卫星信号,那么n>k,通常情况下,n是远大于k。
先从1号虚拟真实卫星信号开始,由惯性辅助解算出1号虚拟真实卫星信号在天线阵坐标系中的到达角,求解1号虚拟真实卫星信号与1号待检测卫星信号在天线阵基线之上的载波相位双差,利用二元假设检验模型及其判决依据,判断1号待检测卫星信号与1号虚拟真实卫星信号之间的r11是否小于ρth,如果小于ρth,则1号待检测卫星信号为真实卫星信号,此时,借助接收机等其它检测手段做进一步的处理;如果不小于ρth,则先将1号待检测卫星信号识别为欺骗干扰信号。然后构建1号虚拟真实卫星信号与2号待检测卫星信号之间的载波相位双差,利用二元假设检验模型及其判决依据,判断1号待检测卫星信号与2号虚拟真实卫星信号之间的r12是否小于ρth,如果小于ρth,则2号待检测卫星信号为真实卫星信号;如果不小于ρth,则先将2号待检测卫星信号识别为欺骗干扰信号;以此类推,根据1号虚拟真实卫星信号的载波相位差,识别n个待检测卫星信号中的真实卫星信号和欺骗干扰信号,理论上,n个待检测卫星信号中仅存在一个对应于1号的真实卫星信号,其余的识别为欺骗干扰信号。同时,经天线阵载波相位双差检测出的真实卫星信号,可结合卫星导航接收机后续的信号处理得到该真实卫星信号所对应的导航卫星PRN(pseudo random noise code,随机噪声码)编号,利用该PRN编号与惯性辅助解算的虚拟真实卫星PRN编号可做进一步欺骗干扰信号识别的校验。
经过一轮的识别,再以2号虚拟真实卫星信号的载波相位差为根据,再次识别n个待检测卫星信号中的真实卫星信号和欺骗干扰信号,理论上,可再次识别出其中的一个待检测卫星信号为真实卫星信号,其余的识别为欺骗干扰信号。遍历k个虚拟真实卫星信号,可以从n个待检测卫星信号中识别出k个真实卫星信号,其余的(n-k)个待检测卫星信号为欺骗干扰信号。
采用该欺骗干扰信号的检测算法,同样的,对于n个待检测卫星信号,算法的复杂度为n*k。而采用传统的对n路卫星信号进行欺骗干扰检测算法需要将任意两路卫星信号的来向进行比较,因此,算法的复杂度为其中,n是远大于k的,因此,本发明提供的欺骗干扰信号检测方法的检测时间相比传统检测时间短。
本发明提供的一种基于惯性辅助的天线阵欺骗干扰信号检测方法,通过计算待检测卫星信号与虚拟真实卫星信号在天线阵基线上载波相位差的不同来识别欺骗干扰,解决了传统载波相位双差检测技术检测时间较长且对单路卫星欺骗干扰信号检测失效的问题;通过设置判决依据中的虚警率和欺骗干扰信号的最小检测概率,可缩小欺骗干扰信号的检测盲区,提升了检测概率,综合检测性能优异。
尽管已描述了本申请的优选实施例,但本领域内的技术人员一旦得知了基本创造概念,则可对这些实施例作出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本申请范围的所有变更和修改。
显然,本领域的技术人员可以对本申请进行各种改动和变型而不脱离本申请的精神和范围。这样,倘若本申请的这些修改和变型属于本申请权利要求及其等同技术的范围之内,则本申请也意图包括这些改动和变型在内。
Claims (5)
1.一种基于惯性辅助的天线阵欺骗干扰信号检测方法,其特征在于,包括:
由惯性辅助解算每一个虚拟真实卫星信号在天线阵坐标系中的到达角;
计算惯性辅助解算的每一个虚拟真实卫星信号的到达角在天线阵不同基线上的载波相位差,以及每一个待检测卫星信号的到达角在天线阵不同基线上的载波相位差;
采用载波相位双差算法构建欺骗干扰信号的二元假设检验模型,并基于二元假设检验模型设定欺骗干扰信号的判决依据;
遍历惯性辅助解算的每一个虚拟真实卫星信号与每一个待检测卫星信号之间的载波相位双差,基于二元假设检验模型及其判决依据,判决每一个待检测卫星信号是否为欺骗干扰信号;
所述计算惯性辅助解算的每一个虚拟真实卫星信号的到达角在天线阵不同基线上的载波相位差,以及每一个待检测卫星信号的到达角在天线阵不同基线上的载波相位差包括:
卫星导航接收机天线阵考虑为正三角形天线阵,设天线阵的基线长度为b,以三个阵元所在的平面为X轴和Y轴所在平面,以A为原点建立直角坐标系,AB基线方向为X轴方向,Z轴指向天顶方向,XYZ轴满足右手定则;
卫星入射信号在天线阵坐标系下的单位方向矢量表示为:
天线阵三个阵元的坐标分别为A=[0,0,0],b=[b,0,0]和相应阵元间基线方向矢量分别为bAB、bAC和bBC,其中,θ为接收机天线阵的俯仰角为,/>为方位角;
虚拟真实卫星信号i与待检测卫星信号j在天线阵基线AB上的载波相位差分别表示为:
其中,为待检测卫星信号载波相位测量值的测量噪声,Ψi为由惯性辅助解算的虚拟真实卫星信号在天线阵坐标系中的到达角;
同理,虚拟真实卫星信号i与待检测卫星信号j在天线阵基线AC上的载波相位差分别表示为:
其中,天线阵为正三角形天线阵,存在:
则虚拟真实卫星信号i与待检测卫星信号j在天线阵基线AB及AC上的载波相位差可以写成:
dΦi=HΛi;
dΦj=HΛj+γj;
其中,
所述采用载波相位双差算法构建欺骗干扰信号的二元假设检验模型包括:
虚拟真实卫星信号i与待检测卫星信号j的载波相位双差矢量能够表示为:
ΔdΦij=dΦi-dΦj=H(Λi-Λj)+(γi-γj)=HRij+ηij;
式中,γi为0,利用载波相位双差检测原理,构建欺骗干扰信号检测的二元假设检验模型为:
其中,H0/H1的假设条件分别为:
H0:参与载波相位双差检测的待检测卫星信号为欺骗干扰信号;
H1:参与载波相位双差检测的待检测卫星信号为真实卫星信号。
2.根据权利要求1所述的欺骗干扰信号检测方法,其特征在于,所述由惯性辅助解算每一个虚拟真实卫星信号在天线阵坐标系中的到达角包括:
利用惯性辅助解算卫星导航接收机天线阵相对于当地导航坐标系的姿态矩阵;
通过卫星导航接收机信号处理得到各颗真实卫星在当地导航坐标系下的方位角和俯仰角;
基于各颗真实卫星在当地导航坐标系下的方位角和俯仰角以及惯性辅助解算的卫星导航接收机天线阵相对于当地导航坐标系的姿态矩阵,通过坐标转换得到每个虚拟真实卫星信号在天线阵坐标下的到达角,所述到达角包括方位角和俯仰角。
3.根据权利要求1所述的欺骗干扰信号检测方法,其特征在于,所述基于二元假设检验模型设定欺骗干扰信号的判决依据包括:
ηij能够表示为:
ηij服从N(0,C)分布,其中协方差矩阵C能够表示为:
相应的,能够求得C-1为:
构造检测量构造的检测量r服从自由度为2的非中心分布,其中,λ=μTC-1μ,μ为天线阵接收欺骗干扰信号情况下Rij的均值;
当给定虚警率PFA时,欺骗干扰信号的检测门限可表示为:
其中,为自由度为中心χ2分布右尾概率的反函数;
则二元假设检验模型的判决依据为:
r>ρth,则判为H0;
r<ρth,则判为H1。
4.根据权利要求3所述的欺骗干扰信号检测方法,其特征在于,还包括通过下式表示对应的欺骗干扰信号检测概率:
其中,表示自由度为2、非中心参数为λ的非中心分布χ2右尾概率。
5.根据权利要求1-4任一项所述的欺骗干扰信号检测方法,其特征在于,所述遍历惯性辅助解算的每一个虚拟真实卫星信号与每一个待检测卫星信号之间的载波相位双差,基于二元假设检验模型及其判决依据,判决每一个待检测卫星信号是否为欺骗干扰信号包括:
基于任一个虚拟真实卫星信号与每一个待检测卫星信号在天线阵不同基线上的载波相位双差,利用二元假设检验模型及其判决依据,判定相对于所述任一个虚拟真实卫星信号,所有待检测卫星信号中的欺骗干扰信号;
遍历所有的虚拟真实卫星信号与每一个待检测卫星信号,判定出所有待检测卫星信号中的欺骗干扰信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110106507.5A CN112782728B (zh) | 2021-01-26 | 2021-01-26 | 一种基于惯性辅助的天线阵欺骗干扰信号检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110106507.5A CN112782728B (zh) | 2021-01-26 | 2021-01-26 | 一种基于惯性辅助的天线阵欺骗干扰信号检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112782728A CN112782728A (zh) | 2021-05-11 |
CN112782728B true CN112782728B (zh) | 2024-03-22 |
Family
ID=75757395
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110106507.5A Active CN112782728B (zh) | 2021-01-26 | 2021-01-26 | 一种基于惯性辅助的天线阵欺骗干扰信号检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112782728B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113640840B (zh) * | 2021-06-30 | 2024-07-02 | 湖南天熠电子科技有限公司 | 一种基于天线阵的非压制gnss欺骗干扰检测抑制方法 |
CN113238257A (zh) * | 2021-07-12 | 2021-08-10 | 航天科工通信技术研究院有限责任公司 | 基于单接收机载波相位差分的gnss欺骗干扰检测方法 |
CN116087994B (zh) * | 2023-04-07 | 2023-07-07 | 中国人民解放军战略支援部队航天工程大学 | 一种基于机器学习的欺骗干扰检测方法 |
CN116859418B (zh) * | 2023-07-05 | 2024-01-26 | 中国人民解放军92728部队 | 基于中继平台监测的卫星导航欺骗干扰检测方法 |
CN117452443B (zh) * | 2023-12-22 | 2024-03-22 | 中国人民解放军国防科技大学 | 一种基于双门限判决的卫星导航欺骗干扰检测方法及装置 |
CN118131277A (zh) * | 2024-05-10 | 2024-06-04 | 厦门大学 | 基于虚拟锚点的gnss欺骗检测*** |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103792550A (zh) * | 2014-02-11 | 2014-05-14 | 哈尔滨工程大学 | 一种基于阵列天线和gps/sins的联合抗干扰方法 |
CN104181572A (zh) * | 2014-05-22 | 2014-12-03 | 南京理工大学 | 一种弹载惯性/卫星紧组合导航方法 |
CN106772455A (zh) * | 2016-12-29 | 2017-05-31 | 南京航空航天大学 | 一种基于惯性信息辅助与参数估计的gnss反欺骗环路跟踪方法 |
CN110244323A (zh) * | 2019-05-24 | 2019-09-17 | 中国科学院光电研究院 | 微、轻型无人机gnss反欺骗***及欺骗信号检测与导航方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
IL169855A (en) * | 2005-07-25 | 2014-05-28 | Elta Systems Ltd | A system and method for locating a receiver location |
-
2021
- 2021-01-26 CN CN202110106507.5A patent/CN112782728B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103792550A (zh) * | 2014-02-11 | 2014-05-14 | 哈尔滨工程大学 | 一种基于阵列天线和gps/sins的联合抗干扰方法 |
CN104181572A (zh) * | 2014-05-22 | 2014-12-03 | 南京理工大学 | 一种弹载惯性/卫星紧组合导航方法 |
CN106772455A (zh) * | 2016-12-29 | 2017-05-31 | 南京航空航天大学 | 一种基于惯性信息辅助与参数估计的gnss反欺骗环路跟踪方法 |
CN110244323A (zh) * | 2019-05-24 | 2019-09-17 | 中国科学院光电研究院 | 微、轻型无人机gnss反欺骗***及欺骗信号检测与导航方法 |
Non-Patent Citations (11)
Title |
---|
GNSS抗欺骗干扰技术研究;耿正霖;聂俊伟;王飞雪;;全球定位***(第04期);全文 * |
Spoofing Detection Using GNSS/INS/Odometer Coupling for Vehicular Navigation;Ali Broumandan;《Sensors》;20180430;全文 * |
利用旋转天线载波相位双差进行欺骗干扰检测技术;张鑫;陈华明;黄仰博;张国柱;欧钢;;武汉大学学报(信息科学版)(第04期);全文 * |
卫星导航欺骗干扰信号仿真与检测关键技术研究;张鑫;《中国博士学位论文全文数据库信息科技辑》;第81-91页 * |
张鑫 ; 陈华明 ; 黄仰博 ; 张国柱 ; 欧钢 ; .利用旋转天线载波相位双差进行欺骗干扰检测技术.武汉大学学报(信息科学版).2016,(第04期),全文. * |
惯导辅助的三元天线阵欺骗干扰检测算法;范广腾;黄仰博;伍微;朱祥维;孙广富;;国防科技大学学报(第02期);第91-95页 * |
惯性信息辅助的卫星导航欺骗检测技术;李四海;《中国惯性技术学报》;20130630;全文 * |
惯性信息辅助的卫星导航欺骗检测技术;李四海;刘洋;张会锁;张晓冬;;中国惯性技术学报(第03期);第336-341页 * |
李四海 ; 刘洋 ; 张会锁 ; 张晓冬 ; .惯性信息辅助的卫星导航欺骗检测技术.中国惯性技术学报.2013,(第03期),第336-341页. * |
耿正霖 ; 聂俊伟 ; 王飞雪 ; .GNSS抗欺骗干扰技术研究.全球定位***.2013,(第04期),全文. * |
采用双天线载波相位差技术的卫星导航接收机抗欺骗方法;黄龙;雍玲;徐博;王飞雪;;国防科技大学学报(第04期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112782728A (zh) | 2021-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112782728B (zh) | 一种基于惯性辅助的天线阵欺骗干扰信号检测方法 | |
Meurer et al. | Robust joint multi-antenna spoofing detection and attitude estimation using direction assisted multiple hypotheses RAIM | |
US8188923B2 (en) | Method of multi-transmitter and multi-path AOA-TDOA location comprising a sub-method for synchronizing and equalizing the receiving stations | |
US20150226834A1 (en) | Evaluating the Position of an Aerial Vehicle | |
Appel et al. | Robust spoofing detection and mitigation based on direction of arrival estimation | |
CN107621645A (zh) | 基于单接收机的欺骗干扰信号检测方法 | |
US11074699B2 (en) | Method for determining a protection radius of a vision-based navigation system | |
CN103592657B (zh) | 一种基于钟差辅助的低可见星下单模raim实现方法 | |
EP2933658B1 (en) | Ground-based system and method to extend the detection of excessive delay gradients using parity corrections | |
US20100271259A1 (en) | Methods and Systems to Diminish False-Alarm Rates in Multi-Hypothesis Signal Detection Through Combinatoric Navigation | |
WO2001081941A1 (en) | Multipath propagation detection and avoidance method and system | |
CN113109843B (zh) | 基于双接收机伪距双差的欺骗信号检测抑制方法及装置 | |
CN108267754B (zh) | 一种基于多接收机间位置相关性的欺骗干扰检测方法 | |
KR102572546B1 (ko) | 단일 주파수 수신기에서 다중 채널에서 발생하는 신호의 차이를 감지하는 장치 및 방법 | |
CN112882068B (zh) | 基于多接收机的gnss抗欺骗干扰方法 | |
KR101957291B1 (ko) | 전자전 지원 시스템(Warfare Support System)의 신호도래방위각 측정 장치 및 그 방법 | |
CN114675305A (zh) | 一种基于压缩感知omp的卫星导航信号干扰抑制方法 | |
CN113640840A (zh) | 一种基于天线阵的非压制gnss欺骗干扰检测抑制方法 | |
Farej et al. | On the Location Accuracy of the Multi-Antenna GPS | |
Zhang et al. | Underwater multi-source DOA tracking using uniform linear array based on improved GM-PHD filter | |
US20210349176A1 (en) | Positioning system and method | |
Zhang et al. | INS-Aided Multi-Antenna GNSS Carrier Phase Double Difference Spoofing Detection | |
CN116774253B (zh) | 一种基于信号到达方向角度差的导航欺骗式干扰检测方法 | |
CN116413675B (zh) | 一种对齐误差下多站雷达协同抗欺骗式干扰方法及装置 | |
CN116087875B (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 |