-
0 引言
-
全球导航卫星系统(Global Navigation Satellite System,GNSS)因其可为各类军民载体提供全天候高精度的定位、测速和授时(Position Velocity Time,PVT)服务而得到广泛应用.由于卫星导航信号的脆弱性,GNSS接收机极易受到恶意人为干扰,使得平台导航系统失效[1].一般GNSS调零天线针对压制式干扰,但从近几年针对无人武器平台的GNSS干扰热门事件来看,欺骗式干扰的威胁越来越大,这是因为欺骗式干扰信号不易被GNSS接收机察觉,并且随着软件无线电技术的发展,欺骗干扰设备所需成本也越来越低.转发式欺骗干扰技术上容易实现,成本低,是目前最主要的欺骗干扰方式.抗欺骗式干扰技术可分为欺骗干扰检测和欺骗干扰抑制:欺骗干扰检测主要实现对接收信号中是否存在欺骗干扰信号进行检测[2]; 欺骗干扰抑制主要实现消除欺骗干扰信号对GNSS系统的影响[3].对欺骗干扰抑制方法的研究相对其检测方法要少很多,特别是不改变GNSS编码结构下的接收机自主抑制方法尤其困难.目前研究较多的欺骗抑制方法有残留信号检测算法,其根据欺骗干扰源无法将真实卫星信号完全抵消,从GNSS接收机接收信号的缓冲样本中将系统重构建的欺骗信号减去,即可实现对欺骗干扰的消除[4],该方法对欺骗干扰抑制的效果属于中等水平,需要额外的存储空间和专门用于重构欺骗信号的通道,并且在欺骗信号功率明显大于真实卫星信号时,因检测不到真实信号而失效.接收机自主完整性监测(Receiver Autonomous Integrity Monitoring,RAIM)是另一种欺骗干扰抑制方法,其把欺骗干扰信号当成故障信号,通过比较伪距测量值而剔除因存在欺骗信号而导致的异常值,从而实现对欺骗干扰的抑制[5],该方法成本低,实现简单,但欺骗信号只能是一个.最可靠的是空域零陷控制算法,由于欺骗干扰大多采用单天线发射,因此信号到达方向(Direction of Arrival,DOA)是区分真实信号和欺骗信号的可靠因素,在判断出欺骗信号来向的基础上[6],可通过零陷形成技术降低欺骗信号方向上的增益,实现对欺骗干扰的抑制[7].该方法虽然能达到良好的抑制效果,但是由于需要进行矩阵求逆运算,算法复杂度较高.针对上述问题,本文提出一种基于调零天线的欺骗干扰抑制算法,能弥补残留信号检测算法和RAIM算法的不足,相对于空域零陷控制算法,运算复杂度更低.
-
本文首先分析欺骗干扰信号模型,然后提出基于调零天线的欺骗干扰抑制算法,最后通过MATLAB仿真验证该算法的可行性.
-
1 欺骗式干扰信号模型
-
GNSS接收机的基本结构如图1所示,整体分为天线、射频前端、信号处理、定位4个部分.通常,GNSS接收机射频前端输出为中频信号,中频信号转换为数字信号后,通过各种信号处理算法实现不同功能,包括了干扰信号检测和抑制等[8].而干扰抑制后的信号,则进入接收机,实现载体定位和授时功能[9].
-
GNSS接收机天线为阵列天线.在存在欺骗干扰的情况下,其单元接收信号包括了欺骗干扰、卫星信号、噪声3个部分,可以表示为
-
其中: 示第i个卫星的真实信号,M表示真实信号的个数; 表示第j个卫星的欺骗式干扰信号,J表示欺骗信号的个数; n(t)为噪声信号.由于GNSS卫星采用了直接序列扩频技术,所以卫星i的信号在接收机处的信号可表示为
-
其中:PTi为接收机接收的卫星i信号的功率; CTi(t)为接收机接收的卫星i信号的伪随机噪声(Pseudo Random Noise,PRN)码; DTi(t)为接收机接收的卫星i信号的导航数据码; τTi为接收机接收的卫星i信号的码相位; f0为中心频率; φTi为载波相移.
-
由于转发式欺骗干扰信号是对真实卫星信号的延时和功率放大,即欺骗信号和真实卫星信号具有相同的信号结构,因此,卫星i对应的欺骗信号可以表示为
-
其中:Δt即为欺骗干扰相对于真实卫星信号的延时; k为欺骗信号对真实卫星信号功率放大系数.
-
根据图1所示的接收机结构,欺骗式干扰的检测算法在信号处理阶段完成,输入信号为数字中频信号.另外,天线采用均匀面阵,但为了叙述简便,假设阵列天线为L+1个阵元构成的均匀线型阵列,信号处理模块把中频模拟信号经过ADC转换为中频数字信号x(k),其中x(k)为L+1维列矢量,每个元素对应一个阵元接收信号,k为采样序号.根据阵列天线基本理论,接收信号表示为
-
式中,v(θ)表示入射角度为θ的信号或者干扰的方向矢量,n(k)为噪声矢量.
-
2 基于调零天线的欺骗干扰抑制方法实现
-
2.1 算法结构
-
基于调零天线的欺骗干扰抑制的目的是在欺骗干扰的来波方向上形成零陷来实现对欺骗干扰的抑制.经典文献[10]指出,基于功率倒置(Power Inversion,PI)算法的调零天线能在强干扰处产生零陷,而对微弱的欺骗信号是不适用的.由于PI算法的调零天线实现简单可靠,因此,本文所提算法的思想是人为构建大功率信号,模拟其从欺骗信号DOA入射,以使调零天线能在该方向生成零陷,从而实现对该方向信号的抑制,即欺骗干扰信号的抑制.欺骗式干扰的入射角度需要通过DOA估计方法得到,比如阵列信号处理中的自适应滤波算法[11].
-
图1 GNSS接收机的通用内部结构
-
Fig.1 General internal structure of GNSS receiver
-
基于调零天线的欺骗干扰抑制算法结构如图2所示,主要包含压制式干扰生成模块、信号求和模块和调零天线模块3部分,其中压制式干扰生成模块和信号求和模块主要是完成调零天线输入信号的产生.基于调零天线的欺骗干扰抑制算法结构具体实现步骤如下:
-
1)压制式干扰生成模块根据欺骗干扰的DOA信息,人工模拟产生J个大功率的带通信号,其中J是接收信号中的欺骗干扰个数;
-
2)信号求和模块将卫星导航接收阵列接收信号和压制式干扰模块产生的人工干扰求和作为调零天线模块的输入;
-
3)调零天线模块对其输入信号进行权值求解和加权求和,最终将零陷对准欺骗干扰信号,从而实现对欺骗干扰的抑制.
-
压制式干扰生成模块和信号求和模块是对调零天线模块输入信号的生成.
-
图2 基于调零天线的欺骗干扰抑制算法结构
-
Fig.2 Structure of anti-spoofing algorithm based on nulling antenna
-
2.2 人工干扰信号生成
-
根据图2所示结构,调零天线新的输入就是真实输入和压制干扰生成模块输出二者的叠加.压制式干扰生成模块根据欺骗式干扰入射角度θj,人工产生任意J个不相干的大功率带通信号:
-
其中,aj(k)为高斯分布随机序列,ωc为卫星信号的中频载波频率,Ts为采样周期,v(θj)为第j个欺骗干扰信号的方向矢量,j=1,···,J.将上述J个人工生成的压制式干扰信号求和得到
-
其中,s(k)是L+1维列矢量.
-
信号求和模块将卫星导航接收机天线阵面的接收信号x(k)和压制式干扰模块产生的人工干扰s(k)进行求和:
-
以上求和结果y(k)作为调零天线的输入.显然,y(k)也是L+1维列矢量.
-
2.3 调零天线结构
-
同样采用L+1个阵元构成的均匀线阵,阵元间距为载波波长的一半的阵列结构对调零天线进行介绍.可用图3表示调零天线的原理结构组成,其主要采用功率倒置算法实现抗干扰功能.调零天线模块输入信号为
-
其中,输入信号包含真实卫星信号、欺骗干扰信号和噪声信号.为了将其形式化表达,将阵元0对应信号y0(k)作为调零天线参考阵元接收到的信号d(k),假设其接收真实卫星信号为sT(k),空间存在J个欺骗干扰信号,噪声信号表示为n0(k),则参考信号d(k)可用下式表示:
-
其余L个阵元对应信号作为辅助阵元信号,其接收信号矢量可写为
-
式中:va(θm),va(θj)分别为真实卫星信号和欺骗干扰信号的方向矢量; n(k)为辅助阵元阵列接收到的噪声信号.
-
图3 调零天线结构
-
Fig.3 Structure of nulling antenna
-
由于真实卫星信号本身强度就相对较弱,再经过扩频调制之后,其信号强度甚至会远低于噪声信号的强度.那么,式(9)可近似表示成
-
式(10)可写为
-
经过自适应滤波之后,非参考阵元阵列的权矢量可表示成
-
经过加权处理后的调零天线输出信号为
-
阵列方向图可以表示为
-
式(15)中的权值和方向矢量包含了第0个阵元的信息,所以
-
2.4 权值求解方法
-
调零天线图3中的权值更新模块通常采用最小均方(Least Mean Square,LMS)算法计算权值.LMS算法的工作原理是使参考阵元和非参考阵元阵列接收信号的均方误差最小化,并通过当前时间的权矢量与调零天线输入迭代计算后一时刻的权矢量值,重复以上运算直到权矢量收敛[12].该算法具体实现步骤如下:
-
1)设置初始取值为0的权矢量wa(0)与迭代步长μ,0<μ<2/λmax,其中λmax为调零天线输入信号协方差矩阵的最大特征值;
-
2)计算此时非参考阵元阵列输出信号;
-
3)利用参考阵元接收信号,计算阵列输出信号z(k)=d(k)-u(k);
-
4)依据快拍数k,对权矢量进行更新;
-
5)回退执行2),直到权矢量收敛;
-
6)权矢量收敛后,将最后的z(k)作为调零天线输出信号.
-
3 算法性能分析
-
空域零陷控制算法和本文提出的基于调零天线的欺骗干扰抑制算法都能达到较好的欺骗抑制效果,但是二者算法复杂度存在区别.
-
空域零陷控制算法满足以下条件即可实现对欺骗干扰的抑制.
-
其中:v0,vi分别为真实信号和欺骗信号的方向矢量; K为接收信号中存在的欺骗信号来向个数.若矩阵
-
则式(17)可写为
-
因此,求得最优权矢量表达式为
-
由式(20)可知,空域零陷控制算法对权矢量的求解涉及到矩阵求逆T-1运算,而通常仅矩阵求逆算法的复乘计算量就能达到L4的数量级,其中L是阵元数.而采用LMS算法计算权值的调零天线只需对权值简单地进行迭代即可确定最优权矢量值,复乘计算量最高达到L3的数量级.当然,两种算法的具体复杂度还跟采样快拍有关.但从以上推导可以得出:基于调零天线的欺骗干扰抑制算法复杂度比空域零陷控制算法复杂度有所降低.
-
4 仿真分析
-
下面对本文提出的欺骗干扰抑制方法进行仿真分析,假设接收机阵列为8阵元均匀线阵,阵元间距为0.5倍波长,仿真假设接收机接收信号包括了一个卫星信号,信号来向分别为30°,两个转发式欺骗干扰信号,信号来向为-40°和60°,以及噪声信号.为了简单,GNSS接收机仅仅研究L1频段的GPS信号,PRN码为长度是1 023的GOLD码,码速率为1.023 MHz,接收信号通过射频前端转换为中频信号,中频频率4.092 MHz,采样频率为37.851 MHz,信号信噪比SNR=-20 dB,干噪比INR=-19 dB.噪声为功率为1的高斯白噪声.
-
用大功率带通信号模拟两个压制式干扰,其中a(k)为均值为0,功率为INR=10 dB的高斯分布随机序列.两个人为干扰的区别是基带信号a(k)独立生成,彼此不相关.其时域波形和频谱图分别如图4、5所示,两个人为干扰时域波形不同,但频谱是相同的.
-
然后,根据式(7)和式(14)对GNSS接收机接收信号进行调零天线抗干扰处理.采用LMS算法计算最优权矢量,图6为根据式(15)计算的调零天线方向图.当不增加人为干扰时,由于欺骗信号功率太弱,无法形成零陷,方向图是水平实线.当采用本文方法后,在-40°和60°两个干扰信号方向,都能够形成零陷.
-
图4 带通信号时域波形
-
Fig.4 Time domain waveform of bandpass signal
-
图5 带通信号频谱
-
Fig.5 Spectrum of bandpass signal
-
欺骗干扰通常和真实卫星信号功率相当,都较弱.调零天线主要目的是减小欺骗干扰功率,保持真实信号功率不变,这些都可以从干扰和真实信号DOA对应的方向图上看出.为此,在相同假设条件下改变阵元数量和人为干扰功率,对真实信号、两个干扰信号天线增益进行仿真,如图7、图8所示.仿真结果说明,真实信号功率和天线数量、人为干扰功率大小没有显著关系.图7显示,欺骗式干扰抑制度在天线数量满足自由度要求后,维持在-55 dB左右,不随天线数量增加而增加.图8显示,干扰功率对干扰抑制度影响显著,干扰功率越大,抑制越多.
-
抗干扰模块在抑制干扰时,不能对真实信号产生较大影响,至少要保证真实信号能够被GNSS接收机中PN码检测器正常检测到相关峰.根据扩频通信接收中相关峰检测方法,对抗干扰前后信号进行分析.图9为未进行欺骗干扰抑制接收机捕获的相关峰,仿真中可以看到接收机捕获到3个相关峰,说明接收信号中存在2个欺骗干扰信号.干扰抑制后,同样进行相关峰检测,如图10所示,则只有一个相关峰,说明干扰信号被有效抑制,并且真实卫星信号得到了有效保护.
-
图6 两种情况调零天线方向图
-
Fig.6 Directional diagram of nulling antenna in two cases
-
图7 天线增益与天线数量关系
-
Fig.7 Relationship between antenna gain and antenna number
-
5 结论
-
通过阵列天线接收信号在空域的处理,可以确定干扰信号角度,然后人为把欺骗式干扰转换为压制式干扰,再采用传统调零天线算法在欺骗干扰方向形成天线方向图零陷,从而实现对欺骗干扰抑制.通过MATLAB仿真实验验证了算法的可行性,对于GNSS抗干扰天线工程设计具有参考价值.
-
图8 天线增益与干扰功率关系
-
Fig.8 Relationship between antenna gain and interference power
-
图9 干扰抑制前接收机捕获相关峰
-
Fig.9 Correlation peak captured by receiver before interference suppression
-
图10 干扰抑制后接收机捕获相关峰
-
Fig.10 Correlation peak captured by receiver after interference suppression
-
参考文献
-
[1] Lu Q,Feng X Z,Zhou C.A detection and weakening method for GNSS time-synchronization attacks[J].IEEE Sensors Journal,2021,21(17):19069-19077
-
[2] Li J Z,Zhu X W,Ouyang M J,et al.Research on multi-peak detection of small delay spoofing signal[J].IEEE Access,8:151777-151787
-
[3] Feng W K,Friedt J M,Goavec-Merou G,et al.Software-defined radio implemented GPS spoofing and its computationally efficient detection and suppression[J].IEEE Aerospace and Electronic Systems Magazine,2021,36(3):36-52
-
[4] Wesson K,Shepard D,Humphreys T,et al.Straight talk on anti-spoofing securing the future of PNT[J].GPS World,2012,3(4):65-70
-
[5] Sun Y,Fu L.A new threat for pseudorange-based RAIM:adversarial attacks on GNSS positioning[J].IEEE Access,2019,7:126051-126058
-
[6] 孙凤林,尹继亮,李凤,等.GNSS欺骗干扰的空时联合检测算法[J].计算机仿真,2022,39(5):60-65;SUN Fenglin,YIN Jiliang,LI Feng,et al.Space-time joint detection algorithm for GNSS spoofing interference[J].Computer Simulation,2022,39(5):60-65
-
[7] Hu Y F,Bian S F,Li B,et al.A novel array-based spoofing and jamming suppression method for GNSS receiver[J].IEEE Sensors Journal,2018,18(7):2952-2958
-
[8] 刘丁浩,吕晶,马蕊,等.卫星导航系统欺骗与抗欺骗技术研究与展望[J].通信技术,2017,50(5):837-843;LIU Dinghao,LÜ Jing,MA Rui,et al.The research and prospect of spoofing and anti-spoofing technology in the satellite navigation system[J].Communications Technology,2017,50(5):837-843
-
[9] Rychlicki M,Kasprzyk Z,Rosiński A.Analysis of accuracy and reliability of different types of GPS receivers[J].Sensors(Basel,Switzerland),2020,20(22):6498-6511
-
[10] Compton R T.The power-inversion adaptive array:concept and performance[J].IEEE Transactions on Aerospace and Electronic Systems,1979,15(6):803-814
-
[11] Zeng H,Ahmad Z,Zhou J W,et al.DOA estimation algorithm based on adaptive filtering in spatial domain[J].China Communications,2016,13(12):49-58
-
[12] Meng D W,Feng Z M,Lu M Q.Anti-jamming with adaptive arrays utilizing power inversion algorithm[J].Tsinghua Science & Technology,2008,13(6):796-799
-
摘要
欺骗式干扰由于实现成本低、干扰能力强,成为GNSS中主流干扰来源,但现有GNSS抗干扰天线仅仅针对压制式干扰.根据欺骗式干扰检测获得的干扰到达方向信息,采用人为在该角度增加压制式干扰的方法,现有GNSS调零天线就可以在欺骗干扰角度形成方向图零陷,实现欺骗干扰抑制.通过对天线方向图和接收机扩频码同步仿真,验证了新方法的有效性.
Abstract
Nowadays spoofing is the dominant threats for GNSS receiver due to its low cost and high efficiency.However,the nulling antenna,which is widely used in GNSS receiver,can only mitigate the interference via high power.According to the DOA information of spoofing,some pseudo interferences are added to the received signal and imping on the antenna array from the same DOAs as the spoofing.Then the following traditional nulling antenna can suppress the spoofing since the beam pattern generate nulls at the DOAs of spoofing.Finally,the simulations of beam pattern and PN code synchronization illustrate the performance of the proposed method.