分布式间歇干扰下基于SMI的GNSS空时自适应处理器性能分析

王解，刘文祥，陈飞强，欧钢; WANG Jie; LIU Wenxiang; CHEN Feiqiang; OU Gang

分布式间歇干扰下基于SMI的GNSS空时自适应处理器性能分析^*

doi: 10.11887/j.cn.202306013

王解^1,2 ，刘文祥¹ ，陈飞强¹ ，欧钢¹

1. 国防科技大学电子科学学院, 湖南长沙 410073

2. 国防科技大学电子对抗学院, 安徽合肥 230037

基金项目: 国家自然科学基金资助项目（42204039）

详细信息

作者简介

王解（1989—），男，安徽安庆人，讲师，博士，E-mail: wangjiefjj@163.com

中图分类号: P228

文献标识码: A

文章编号: 1001-2486(2023)06-090-10

Performance analysis of SMI based GNSS space-time adaptive processing under distributed intermittent interferences

WANG Jie^1,2 ， LIU Wenxiang¹ ， CHEN Feiqiang¹ ， OU Gang¹

1. College of Electronic Science and Technology, National University of Defense Technology, Changsha 410073 , China

2. College of Electronic Engineering, National University of Defense Technology, Hefei 230037 , China

摘要

针对全球导航卫星系统接收机面临的分布式间歇干扰威胁，全面准确地分析了基于采样矩阵直接求逆方法的空时自适应处理（space-time adaptive processing，STAP）的抗干扰性能。在建立性能分析模型的基础上，根据间歇干扰造成的采样协方差矩阵不匹配情况，通过分类评估、梳理归纳、理论推导全面分析了不同情况下STAP的性能。分析结果表明在采样长度小于闪烁周期且使用预采用处理方法时，空时自适应处理器会出现漏采样，干扰抑制性能会急剧下降，漏采样出现的频率为各干扰闪烁频率之和。数值和仿真分析验证了理论分析的结果，在此基础上讨论了考虑分布式间歇干扰时STAP的优化设计。

关键词

分布式间歇干扰 / 直接矩阵求逆 / GNSS / 空时自适应处理 / 性能分析

Abstract

Considering the distributed intermittent interference threat faced by global navigation satellite system receiver, the anti-jamming performance of STAP(space-time adaptive processing) based on sample matrix inversion method was comprehensively and accurately analyzed. Based on the performance analysis model, according to the mismatch of sampling covariance matrix caused by intermittent interferences, the performance of STAP under different conditions was comprehensively analyzed through classification evaluation, combing induction and theoretical derivation. The analysis results show that when the sampling length is less than the scintillation period and the pre-sample method is used, the missed sampling of space-time adaptive processor will occur and the interference suppression performance will decrease sharply, and the frequency of missed sampling is the sum of the scintillation frequencies of each interference. Numerical and simulation analyses verify the results of analysis, and based on this, the optimal design of STAP in the face of distributed intermittent interference is discussed.

Keywords

distributed intermittent interferences / sample matrix inversion / GNSS / space-time adaptive processing / performance analysis

1 信号模型和阵列接收机处理流程 1.1 分布式间歇干扰模型 1.2 阵列导航接收机抗干扰处理流程 2 分布式间歇干扰影响评估 2.1 分布式间歇干扰的影响机理和关键因素 2.2 干扰影响的分类评估 2.2.1 采样长度大于闪烁周期 2.2.2 采样长度小于闪烁周期 3 自适应处理器性能理论分析 4 数值分析和仿真验证 5 讨论 6 结论

面对日益严重的干扰威胁，抗干扰全球导航卫星系统（global navigation satellite system，GNSS）接收机技术不断得到发展并被应用于不同领域，其中，利用天线阵抗干扰技术被认为最为有效，能够有效抑制多个窄带和宽带干扰信号。在干扰信号连续平稳的条件下，天线阵抗干扰处理表现出较好的干扰抑制性能。然而文献^[1-4]的研究工作表明当前GNSS面临的干扰威胁多种多样并在不断地上演着，非平稳和非连续干扰是所面临威胁的一部分。目前很少有文献分析非连续信号对天线阵接收机的影响。文献^[5]首先提醒关注非平稳干扰对天线阵卫星接收机的干扰。文献^[6]的研究表明天线阵技术不能很好应对非连续的脉冲干扰，必须进一步使用脉冲消隐技术才能应对该干扰场景。这些文献虽然关注了导航频段的非平稳或非连续信号干扰，但是并没有针对非平稳干扰对接收机性能影响展开全面分析。本文针对分布式间歇干扰对阵列GNSS接收机抗干扰处理环节的影响进行精确分析。

分布式间歇干扰是指多个干扰信号在空间上由不同方向入射，在时间上断续存在，其中的单个干扰呈现脉冲干扰特性。分布式间歇干扰之所以对基于递推算法的阵列处理影响较大，是因为信号突变使得滤波器收敛速度下降甚至无法收敛^[7]，但是对于使用矩阵直接求逆（sample matrix inversion，SMI）算法的抗干扰处理，其收敛速度和采样频率以及长度有关，此时分布式间歇干扰将导致采样协方差矩阵和真实协方差矩阵不匹配。在通信领域也存在由通道间干扰形成的间歇式分布干扰，研究表明干扰将导致基于SMI的阵列处理性能下降^[8-10]，并给出了优化采样的方法。文献^[11-12]分析了在杂波环境下，雷达中空时自适应处理（space-time adaptive processing，STAP）方法的性能，文献^[13-15]也分析了在特殊的样本不匹配情况下，阵列雷达处理中信噪比损失概率分布和检测概率。但是GNSS 阵列处理与通信和雷达阵列处理有明显的差别，例如，雷达处理中采样协方差矩阵只能通过相邻距离单元采样数据计算；GNSS应用中采样协方差矩阵的计算则有多种选择方式，这些差别让分布式间歇干扰对GNSS阵列处理的影响明显不同。

在分析分布式间歇干扰的影响时，本文假设干扰的总个数小于阵列的自由度。分布式间歇干扰可由无意干扰或有意干扰形成，而大量无意干扰形成分布式间歇干扰场景比较少见，且有意干扰追求用尽可能少的干扰个数实现最优干扰，因此该假设是符合实际的。并且文章将阵列处理模块作为分析的目标，主要关注卫星信号相关前的干扰抑制性能，将输出功率和信噪比作为性能评估指标。

文章首先介绍分布式间歇干扰的模型并简要介绍GNSS阵列处理的一般流程，然后在分析分布式间歇干扰影响机理和关键因素的基础上，对干扰的影响进行分类评估，在建立分类的不匹配模型之后进行了归纳并给出理论分析结果，并通过数值仿真进行验证，最后提出抗干扰处理设计和优化意见。

1 信号模型和阵列接收机处理流程

1.1 分布式间歇干扰模型

电磁环境中存在大量间歇干扰，在导航频段内，来自塔康导航系统的脉冲信号就属于间歇干扰^[6]，若间歇干扰由不同方向入射就会构成分布式间歇干扰场景。假设各干扰信号为矩形脉冲调制信号，其载波及带宽和目标卫星信号一致，并假设干扰信号为上升时间和下降时间均为0的理想开关信号。其时域波形描述为：

j (t) = c (t) c o s (2 π f_{0} t) W (t)

(1)

式中：c（t）为宽带噪声调制信号；f₀为载波频率；W（t）为方波，其表达式为

(2)

其中，A为信号的幅值，τ为有效脉冲宽度，T为脉冲周期。

在进行影响分析时，假定干扰相对接收机的空间分布是固定的，在时间上，各干扰信号的到达时间是随机的。给出其时间和空间（到达角：俯仰角θ，方位角φ）的示意图如图1所示。

图1空间分布式间歇干扰信号示意图

Fig.1Schematic diagram of spatially distributed intermittent interferences

1.2 阵列导航接收机抗干扰处理流程

在信号处理流程上，基于数字信号处理的阵列接收机与单天线接收机相比增加了抗干扰处理环节，如图2所示。其信号处理主要包括了射频前端、AD转换、抗干扰处理、通道处理等环节，间歇干扰对上述各环节都会产生一定的影响，文章主要关注间歇干扰对抗干扰处理环节的影响。

图2阵列接收机信号的处理流程

Fig.2Signal processing flow of array receiver

基于SMI方法的导航接收机STAP抗干扰处理中，利用不同的优化准则求解抗干扰滤波器系数，以最小方差无失真响应（minimum variance distortionless response，MVDR）准则为例，其通过使输出功率为0，并约束卫星信号增益为1，求解滤波器系数：

w_{o p t} = \frac{R_{y}^{- 1} a (θ_{0})}{a^{H} (θ_{0}) R_{y}^{- 1} a (θ_{0})}

(3)

式中：a（θ₀）是卫星信号的导向矢量，由信号入射角度θ₀和阵型计算得来；R_y是真实的噪声和干扰信号的协方差矩阵，一般是无法准确获取的。由于卫星信号功率在噪声功率之下，因此阵列接收机在处理过程中，认为采样信号为噪声和干扰信号，在采样信号为平稳信号的假设下，将采样信号的协方差矩阵

{\hat{R}}_{y}

代替真实协方差矩阵进行解算。

假设采样长度为T₀时，采样协方差矩阵可以表示为：

\begin{matrix} {\hat{R}}_{y} = \frac{1}{T_{0}} \int_{t}^{t + T_{0}} [(s (t) + s_{j} (t) + n (t)) \cdot \\ {(s (t) + s_{j} (t) + n (t))}^{H}] d t \end{matrix}

(4)

式中：s（t）为天线阵接收的理想卫星信号矢量，

s （ t ） = a (θ_{0}) s （ t ）; s_{j} （ t ） = \sum_{n = 1}^{N} s_{j n} （ t ）

为天线阵接收的干扰信号矢量，

s_{j n} （ t ） = a (θ_{j n}) s_{j n} （ t ） ， a (θ_{j n})

为第n个干扰信号的导向矢量；n（t）为噪声信号。由于在卫星导航应用中，卫星信号功率远小于噪声功率，并假设各干扰信号不相关，因此式（4）可写为：

{\hat{R}}_{y} \approx R_{j} + R_{n} = R_{j} + σ^{2} I

(5)

式中，σ² 为噪声功率，R_j表示干扰信号的采样协方差矩阵，R_n为噪声信号协方差矩阵。

2 分布式间歇干扰影响评估

2.1 分布式间歇干扰的影响机理和关键因素

在分布式间歇干扰场景下，干扰信号平稳假设不再成立，使用训练样本计算的协方差矩阵和实际处理信号的协方差矩阵产生了不匹配，造成信噪比损失^[13]，可表示为：

\begin{matrix} ρ ({\hat{R}}_{y}) = \frac{{[a^{H} (θ_{0}) {\hat{R}}_{y}^{- 1} a (θ_{0})]}^{2}}{a^{H} (θ_{0}) R_{y}^{- 1} a (θ_{0})} \cdot \\ \frac{1}{a^{H} (θ_{0}) {\hat{R}}_{y}^{- 1} R_{y} {\hat{R}}_{y}^{- 1} a (θ_{0})} . \end{matrix}

(6)

式中，

ρ ({\hat{R}}_{y})

是小于1的正实数，从抗干扰的角度来看，需要使信噪比损失最小，从而达到好的抗干扰效果。

根据信噪比损失

ρ ({\hat{R}}_{y})

的表达式，可以看出

{\hat{R}}_{y}

的不匹配程度决定了分布式间歇干扰的干扰效果。在忽略信号功率和到达角分布的情况下，分析分布式间歇干扰的模型和STAP过程。抗干扰处理过程中采样矩阵计算方式对不匹配程度影响较大，当被处理信号和采样信号是相同的数据集时，称之为即时采样计算，当采样信号集为处理信号集的前一数据块时，称之为预采样计算，两种采样处理方式面对分布式间歇干扰的表现差别较大。从采样信号协方差矩阵的计算式（4）来看，间歇干扰的闪烁周期和信号采样长度对抗干扰处理效果影响也较大。

2.2 干扰影响的分类评估

在分析影响干扰效果的关键因素之后，首先推导了不同场景下训练样本计算协方差矩阵和实际处理信号的真实协方差矩阵不匹配情况，然后总结了间歇干扰下采样协方差矩阵不匹配的几种表现形式，理论分析了不匹配可能造成的信噪比损失，并对出现的概率和影响作了定性分析。

根据上文关键影响因素的分析，下面根据采样长度和处理方式进行具体分析。

2.2.1 采样长度大于闪烁周期

在采样长度大于闪烁周期的情况下，无论是预采样处理还是即时采样处理都有：

\begin{matrix} {\hat{R}}_{y} \approx R_{n} + \sum_{n = 1}^{N} d_{n} p_{j n} a (θ_{j n}) a^{H} (θ_{j n}) \\ = R_{n} + \sum_{n = 1}^{N} d_{n} R_{j n} \end{matrix}

(7)

式中：p_j_n为干扰n的功率；d_n为第n个干扰的占空比，为简化分析，假设所有干扰占空比相同；

R_{j n} = p_{j n} a (θ_{j n}) a^{H} (θ_{j n})

。

使用上述采样矩阵进行抗干扰处理后，由于不同时刻天线阵计算协方差矩阵与真实协方差矩阵的不同，因此会在不同时刻产生不同的信噪比损失，而信噪比损失由不同时刻信号的开关决定。

处理信号的真实协方差矩阵为：

R_{y} = R_{n} + \sum_{n \in {X}} R_{j n}

(8)

式中，{X} 表示信号为开状态的干扰集合。特殊地，干扰信号全部为开或者全部为关（如图3中的数据段①和④）。

2.2.2 采样长度小于闪烁周期

当采样长度远小于闪烁周期时，需要区分不同采样处理方式，如图4所示。当处理方式是预采样处理时，会出现采样信号和处理信号的不匹配，可分为两种情况。

一种不匹配是采样信号不包含处理数据中的一个或者几个干扰信号，导致无法抑制未被采样的干扰信号（如图4所示，若使用数据集a计算的结果处理数据集b，则干扰1无法被抑制），此时采样信号计算协方差矩阵为：

图3采样长度远大于闪烁周期时采样数据集和处理数据集示意图

Fig.3Sampling data set and processing data set when the sampling length is much larger than the scintillation period

图4采样长度远小于闪烁周期时采样数据集和处理数据集示意图

Fig.4Sampling data set and processing data set when the sampling length is much less than the scintillation period

\begin{matrix} {\hat{R}}_{y} \approx R_{n} + \sum_{n \in {X}} p_{j n} a (θ_{j n}) a^{H} (θ_{j n}) \\ = R_{n} + \sum_{n \in {X}} R_{j n} \end{matrix}

(9)

处理信号的真实协方差矩阵为：

R_{y} = R_{n} + \sum_{n \in {X + Y}} R_{j n}

(10)

式中，{Y}表示处理信号中未被采样的干扰信号。

另一种不匹配是采样信号中存在某个或某些干扰信号，处理信号中没有相应的干扰信号（如图4所示，若使用数据集c计算的协方差矩阵的结果处理数据集d，则处理数据中不包含干扰2），此时采样信号计算协方差矩阵为：

{\hat{R}}_{y} \approx R_{n} + \sum_{n \in {X + Z ∣} p_{j n} a (θ_{j n}) a^{H} (θ_{j n}) = R_{n} + \sum_{n \in {X + Z}} R_{j n}

(11)

式中，{Z}表示被采样但处理信号中不包含的干扰信号。处理信号的真实协方差矩阵为：

\begin{matrix} R_{y} = R_{n} + \sum_{n \in {X}} p_{j n} a (θ_{j n}) a^{H} (θ_{j n}) \\ = R_{n} + \sum_{n \in {X}} R_{j n} \end{matrix}

(12)

两种不匹配的情况都会产生信噪比损失，发生的频率为各干扰源闪烁频率之和。

使用即时采样时，只有在处理数据段中包含干扰信号跳变过程才会造成协方差矩阵不匹配，在处理该数据段时采样信号为：

\begin{matrix} {\hat{R}}_{y} \approx R_{n} + \sum_{n \in {X}} p_{j n} a (θ_{j n}) a^{H} (θ_{j n}) + \\ d_{p} p_{j y} a (θ_{j y}) a^{H} (θ_{j y}) \\ = R_{n} + \sum_{n \in {X}} R_{j n} + d_{p} R_{j y} \end{matrix}

(13)

式中，下标jy表示采样集中包含跳变的干扰，d_p为被采样长度占整个闪烁周期的比例。若闪烁周期很小，则在采样长度远小于闪烁周期情况下对信号接收和处理的影响几乎可以忽略不计。

通过以上梳理，全面分析了分布式间歇干扰导致的采样数据计算协方差矩阵和处理信号真实协方差矩阵的不匹配情况。根据分析的结果，在进行干扰影响的评估时，区分采样长度和间歇干扰的闪烁周期对比的两种情况。若采样长度远大于干扰闪烁周期，则计算闪烁周期内信号处理的平均信噪比损失；若采样长度远小于干扰闪烁周期，则分别计算不同情况的信噪比损失以及发生的频率。表1总结了几种不匹配情况和影响评估方法。

进一步，可将分布式间歇干扰造成的计算协方差矩阵和处理信号的真实协方差矩阵的不匹配情况归纳为对干扰信号采样不充分、漏采样和过采样，漏采样是采样不充分的特殊情况。

采样不充分时，采样信号和被处理信号的协方差矩阵的关系为：

R_{y} = {\hat{R}}_{y} + \sum_{n \in {Y}} (1 - r_{j n}) p_{j n} a (θ_{j n}) a^{H} (θ_{j n})

(14)

式中：{Y} 表示未被充分采样的干扰信号；r_j_n表示信号采样计算功率和实际功率的比值，r_j_n=0是表示该信号完全未被采样，属于漏采样。表1中第4，5种情况均属于采样不充分，其中第5种情况下r_j_n=0，属于漏采样。

表1协方差矩阵不匹配情况分类及影响评估方法

Table1Classification of covariance matrix mismatch and its impact assessment method

过采样时，采样信号和实际信号的协方差矩阵的关系为：

R_{y} = {\hat{R}}_{y} - \sum_{n \in {Z}} r_{j n} p_{j n} a (θ_{j n}) a^{H} (θ_{j n})

(15)

式中：{Z} 表示过采样的干扰信号；特殊地，r_j_n=1表示采样数据段内干扰信号始终存在，但处理数据段没有该干扰信号，称之为完全过采样。表1第3，6种情况属于过采样。

表1中的第1，2种情况下，由于采样长度远大于闪烁周期，一个采样长度内所有干扰信号都被采样，但是闪烁干扰占空比不为1，导致采样协方差矩阵不能反映真实瞬时功率，因此在干扰信号开关的不同时刻，其协方差矩阵不匹配既有采样不充分也有过采样。

尽管上文对分布式间歇干扰产生的协方差矩阵不匹配进行了分类和精确推导，但仍然不能得出干扰对信噪比产生的损失情况，接下来通过理论推导了各类情况下自适应处理器的性能，并通过设置典型场景，利用数值和仿真精确评估间歇干扰的影响。

3 自适应处理器性能理论分析

第一类情况：采样不充分。协方差矩阵的不匹配表示为：

R_{y} = {\hat{R}}_{y} + k_{1} q q^{H}

(16)

式中，q是被部分采样的干扰信号Q的单位导向矢量，k₁表示信号漏采样导致的功率损失大小。

由于干扰信号Q已经被部分采样，因此有：

R_{y} = {\hat{R}}_{y} + k_{1} q q^{H} = (R_{x} + k_{0} q q^{H}) + k_{1} q q^{H}

(17)

{\hat{R}}_{y} = R_{x} + k_{0} q q^{H}

(18)

其中，k₀表示采样计算所得的干扰信号Q的功率，R_x表示除干扰信号Q之外其他干扰信号和噪声信号的协方差矩阵。

由于干扰信号功率远大于噪声功率，因此即使是部分被采样，依然可以假定

k_{0} ≫ σ^{2}

。首先对

{\hat{R}}_{y}

进行特征值分解，在干扰源数量P小于阵元个数的情况下：

\begin{matrix} {\hat{R}}_{y} = \sum_{m = 1}^{N} λ_{m} e_{m} e_{m}^{H} \\ \approx \sum_{m = 1}^{P} λ_{m} e_{m} e_{m}^{H} + σ^{2} \sum_{m = P + 1}^{N} e_{m} e_{m}^{H} \\ = U_{j} Σ_{j} U_{j}^{H} + σ^{2} \sum_{m = P + 1}^{N} e_{m} e_{m}^{H} \end{matrix}

(19)

式中：Σ_j=diag{λ₁，λ₂，···，λ_P}是一个P维对角矩阵，其对角线上的元素对应

{\hat{R}}_{y}

的P个最大的特征值；U_j=[e₁，e₂，···，e_P] 是一个n×P维矩阵，其中的列向量是

{\hat{R}}_{y}

的P个最大的特征值对应的特征向量，这些特征向量张成的子空间也就是干扰子空间，等价于P个干扰对应的导向矢量a₁，a₂，···，a_P 张成的子空间，即span{e₁，e₂，···，e_P}=span{a₁，a₂，···，a_P}，e_P₊₁，e_P₊₂，···，e_N是

{\hat{R}}_{y}

的N-P个较小特征值对应的特征向量，这些向量张成了与干扰子空间正交的噪声子空间。由于λ_m远小于σ²，

{\hat{R}}_{y}^{- 1} = \sum_{m = 1}^{n} \frac{1}{λ_{m}} e_{m} e_{m}^{H} \approx \frac{1}{σ^{2}} \sum_{m = P + 1}^{n} e_{m} e_{m}^{H}

(20)

显然干扰信号Q的导向矢量q属于干扰子空间，因此，

(21)

根据矩阵求逆引理和式（20）、式（21）得

\begin{matrix} R_{y}^{- 1} = {({\hat{R}}_{y} + k_{1} q q^{H})}^{- 1} \\ = {\hat{R}}_{y}^{- 1} - \frac{k_{1}}{1 + k_{1} q^{H} {\hat{R}}_{y}^{- 1} q} {\hat{R}}_{y}^{- 1} q q^{H} {\hat{R}}_{y}^{- 1} \\ \approx {\hat{R}}_{y}^{- 1} \end{matrix}

(22)

式（22）代入式（6）可知在第一类情况下ρ_a1=1，信噪比基本无损失。

第二类情况: 漏采样。协方差矩阵的不匹配表示为：

R_{y} = {\hat{R}}_{y} + k_{2} v v^{H}

(23)

式中，v是完全未被采样的干扰信号V的单位导向矢量，k₂表示干扰信号V的功率。

由于干扰信号V完全未被采样，因此，式（22）在该情况下不成立。此时被漏采样的干扰方向未被置零，会导致处理后的输出中包含该干扰信号，信噪比损失ρ_a2取决于干扰信号导向矢量、功率等因素。

第三类情况: 过采样。协方差矩阵的不匹配表示为：

R_{y} = {\hat{R}}_{y} - b u u^{H}

(24)

式中，u为被过采样干扰信号U的单位导向矢量，b表示干扰信号U的功率。此时u属于

{\hat{R}}_{y}

对应的干扰子空间，根据式（22）的推导过程，同样有

R_{y}^{- 1} \approx {\hat{R}}_{y}^{- 1}

，因此在过采样下信噪比损失较小，特殊地，在完全过采样时，当干扰信号U的入射方向和卫星信号方向接近或一致时，此方向信号增益被置零，将导致

{\hat{R}}_{y}^{- 1} a (θ_{0}) \to 0

，因此在该方向会造成较大的信噪比损失。

对照表1，当采样长度远大于闪烁周期时，不同时刻干扰信号既有采样不充分又有过采样，通过上述推导可知该情况下信噪比基本无损失；当采样长度远小于闪烁周期且即时采样处理时，存在过采样和采样不充分，此时信噪比也基本无损失；当采样长度远小于闪烁周期且预采样处理时，若采样长度结束后处理信号段存在干扰信号关，则出现干扰信号过采样，此时信噪比基本无损失，但是若采样长度结束后处理信号段存在干扰信号开，则出现干扰信号完全未被采样，此时会有明显的信噪比损失。

综上，通过分类、归纳和理论分析表明闪烁干扰导致的采样计算协方差矩阵和真实协方差不匹配只有在一种情况下才会造成明显的信噪比损失，即抗干扰处理器采样长度远小于闪烁周期并使用预采样处理，在采样结束后处理信号段存在干扰信号开时发生，处理该数据段时干扰信号不能完全被抑制。该情况出现的频率为干扰信号闪烁频率之和，此时信噪比损失与干扰信号导向矢量、功率、干扰分布等因素相关，下文将通过数值仿真进行定性和定量分析。

4 数值分析和仿真验证

文章分析时使用四阵元中心圆阵抗干扰接收机，抗干扰处理过程中时域抽头为10个，采样率为40 MHz。设置3个均为宽带脉冲调制的干扰信号，周期为500 μs，占空比为50%，干扰到达角分别为[1°，30°]、[5°，60°]和[10°，120°]，干扰功率为-68 dBm，卫星信号功率为-128 dBm，噪声功率为-100 dBm，各脉冲上升沿到达时间不同，如图1所示。

在以上设置的基础上，首先通过数值分析第3节中不同情况下的信噪比损失。

第一类情况：以图3中的数据段①为例，此时3个干扰信号状态全为开，采样长度远大于脉冲信号周期，干扰信号功率被平均，协方差矩阵不匹配体现为3个干扰信号采样不充分。此时不同方向的信噪比损失如图5所示。结果表明此时信噪比损失小于0.5 dB，证明ρ_a1≈1，与理论分析结果一致。

图5采样不充分时的不同方向信噪比损失

Fig.5Loss of signal-to-noise ratio in different directions with insufficient sampling

第二类情况：以图4中数据集a和数据集b为例，使用预采样处理，数据集b中的干扰1漏采样。此时不同方向的信噪比损失如图6所示，干扰1未被抑制导致在部分方向上信噪比损失达到25 dB以上，验证了理论分析的结果。

图6漏采样时不同方向的信噪比损失

Fig.6Loss of signal-to-noise ratio in different directions with missed sampling

第三类情况：场景1，以图3中的数据集④为例，干扰信号为关状态，采样矩阵中则包含3个干扰信号，不同方向的信噪比损失如图7所示，此时，在3个干扰的方向信噪比有明显损失，但是在其他方向损失5 dB左右；场景2，若以图4中数据集c和数据集d为例，数据集c作为采样数据包含干扰1，数据集d作为处理数据不包含干扰1，不同方向的信噪比损失如图8所示，此时，同样只有在干扰1的方向信噪比有明显损失，其他方向损失较小。以上均为完全过采样（处理信号中完全不包括干扰信号），非完全过采样的情况和采样不充分的情况类似，信噪比基本无损失。

图7过采样时不同方向信噪比损失（场景1）

Fig.7Loss of signal-to-noise ratio in different directions with redundant sampling (scenario 1)

图8过采样时不同方向信噪比损失（场景2）

Fig.8Loss of signal-to-noise ratio in different directions with redundant sampling (scenario 2)

以上数值分析的结果和理论分析的结果一致，只有在漏采样时，信噪比有较大的损失，同时在完全过采样时，干扰信号来向的信噪比损失较大，这是由于该方向的卫星信号也被抑制。用抗干扰处理仿真系统也可验证上述结论，在上述参数设置的基础上，设置卫星信号的来向为[80°，220°]，仿真时长20 ms。进行不同采样长度和采样方式下的仿真处理，包括采样长度为1 ms的即时处理和预处理，采样长度为100 μs的即时处理和预处理。不同设置下的抗干扰处理后输出功率和卫星信号捕获处理相关积分结果如图9~12所示。

根据图9和图10的仿真结果可知，当采样长度为1 ms时，不论是预采样处理还是即时采样处理都能较好抑制干扰，能够完成卫星信号捕获。

图9采样长度为1 ms的即时采样处理结果

Fig.9Results of immediate sampling processing with 1 ms sampling length

图10采样长度为1 ms的预采样处理结果

Fig.10Results of pre-sampling processing with 1 ms sampling length

在采样长度为100 μs时，即时采样处理也能够很好地抑制干扰并取得很好的捕获结果（如图11所示）。但是在预采样处理时，由于有漏采样出现导致漏采样干扰不能被抑制，抗干扰性能急剧下降，导致信号捕获处理不成功（如图12所示）。上述仿真结果与理论分析的结果一致。

图11采样长度为100 μs的即时采样处理结果

Fig.11Results of immediate sampling processing with 100 μs sampling length

图12采样长度为100 μs的预采样处理结果

Fig.12Results of pre-sampling processing with 100 μs sampling length

5 讨论

通过数值和仿真分析较好地验证了分类分析、归纳、理论推导的准确性，分布式间歇干扰对基于SMI的抗干扰处理的干扰抑制性能需要区分不同的处理方式和采样长度。只有当采样长度小于闪烁周期且使用预采样处理时，面对分布式间歇干扰，抗干扰处理器性能才会明显下降。从抗干扰角度看，在干扰闪烁周期不确定的情况下如何选择最优的抗干扰处理方法尤为重要^[16-17]。

若选择较长的采样长度，则需要使用预采样处理的方式，因为此时若选择即时采样处理则会导致较大的处理延迟，更重要的是所需要的存储空间会大大增加。采用预采样处理的风险在于不能保证采用长度大于闪烁周期，避免漏采样，特别是不同干扰闪烁周期不一致的情况下。且在高动态情况下长采样长度还会导致干扰在空间扩展，消耗空间自由度，降低抗干扰性能。一种改善方法是先估计干扰闪烁周期再自适应调整采样长度。

在不能确定闪烁周期的情况下，选择小的采样长度并且使用即时采样处理方法可以获得较好的抗干扰效果。从理论上来说，在假设信号平稳的前提下，采样协方差矩阵的计算样本点数大于2M-3 就可以使协方差矩阵不匹配造成的信噪比损失小于3 dB，其中M表示协方差矩阵维数。但是即时采样处理会带来信号处理的延迟，并且需要额外的内存来预存采样信号，在接收机对处理延迟不敏感应用中可以采用即时采样处理。若脉冲干扰导致信号平稳的假设不成立，采样协方差矩阵计算样本数与信噪比损失之间的关系将发生变化^[18]，最小计算样本数需要增加。

综上，面对分布式间歇干扰，采样长度和处理方式有两种选择：一是使用预采样处理提高采样长度；二是使用即时采样处理且尽量小的采样长度。但是这种选择在仅评估干扰抑制能力时是较好的选择，如果考虑相关后的性能，则需要进行进一步的评估，因为面对干扰信号的闪烁，最优权值不断变化，会导致不同的信号失真，造成负面的影响，这些影响如何需要进一步的分析评估^[19]。因此，在改善参数设置的基础上，如何结合其他的抗干扰方法消除闪烁干扰仍有待研究。

6 结论

文章全面精确分析了在不同情况下，基于SMI的STAP抗干扰面对分布式间歇干扰的抑制性能，将分布式间歇干扰造成的协方差矩阵不匹配归纳分类，理论分析表明在采样长度小于闪烁周期且使用预采样处理方法时，会出现漏采样，干扰抑制性能会急剧下降。漏采样出现的频率为各干扰闪烁频率之和，数值和仿真分析验证了理论分析的结果，在此基础上讨论了面对分布闪烁干扰如何设置抗干扰参数，以及可能面临的挑战，为应对分布式间歇干扰提供了重要参考。

图1空间分布式间歇干扰信号示意图

Fig.1Schematic diagram of spatially distributed intermittent interferences

下载: 全尺寸图片

图2阵列接收机信号的处理流程

Fig.2Signal processing flow of array receiver

下载: 全尺寸图片

图3采样长度远大于闪烁周期时采样数据集和处理数据集示意图

Fig.3Sampling data set and processing data set when the sampling length is much larger than the scintillation period

下载: 全尺寸图片

图4采样长度远小于闪烁周期时采样数据集和处理数据集示意图

Fig.4Sampling data set and processing data set when the sampling length is much less than the scintillation period

下载: 全尺寸图片

图5采样不充分时的不同方向信噪比损失

Fig.5Loss of signal-to-noise ratio in different directions with insufficient sampling

下载: 全尺寸图片

图6漏采样时不同方向的信噪比损失

Fig.6Loss of signal-to-noise ratio in different directions with missed sampling

下载: 全尺寸图片

图7过采样时不同方向信噪比损失（场景1）

Fig.7Loss of signal-to-noise ratio in different directions with redundant sampling (scenario 1)

下载: 全尺寸图片

图8过采样时不同方向信噪比损失（场景2）

Fig.8Loss of signal-to-noise ratio in different directions with redundant sampling (scenario 2)

下载: 全尺寸图片

图9采样长度为1 ms的即时采样处理结果

Fig.9Results of immediate sampling processing with 1 ms sampling length

下载: 全尺寸图片

图10采样长度为1 ms的预采样处理结果

Fig.10Results of pre-sampling processing with 1 ms sampling length

下载: 全尺寸图片

图11采样长度为100 μs的即时采样处理结果

Fig.11Results of immediate sampling processing with 100 μs sampling length

下载: 全尺寸图片

图12采样长度为100 μs的预采样处理结果

Fig.12Results of pre-sampling processing with 100 μs sampling length

下载: 全尺寸图片

表1协方差矩阵不匹配情况分类及影响评估方法

Table1Classification of covariance matrix mismatch and its impact assessment method

下载: 全尺寸图片

图1空间分布式间歇干扰信号示意图

Fig.1Schematic diagram of spatially distributed intermittent interferences

图2阵列接收机信号的处理流程

Fig.2Signal processing flow of array receiver

图3采样长度远大于闪烁周期时采样数据集和处理数据集示意图

Fig.3Sampling data set and processing data set when the sampling length is much larger than the scintillation period

图4采样长度远小于闪烁周期时采样数据集和处理数据集示意图

Fig.4Sampling data set and processing data set when the sampling length is much less than the scintillation period

图5采样不充分时的不同方向信噪比损失

Fig.5Loss of signal-to-noise ratio in different directions with insufficient sampling

图6漏采样时不同方向的信噪比损失

Fig.6Loss of signal-to-noise ratio in different directions with missed sampling

图7过采样时不同方向信噪比损失（场景1）

Fig.7Loss of signal-to-noise ratio in different directions with redundant sampling (scenario 1)

图8过采样时不同方向信噪比损失（场景2）

Fig.8Loss of signal-to-noise ratio in different directions with redundant sampling (scenario 2)

图9采样长度为1 ms的即时采样处理结果

Fig.9Results of immediate sampling processing with 1 ms sampling length

图10采样长度为1 ms的预采样处理结果

Fig.10Results of pre-sampling processing with 1 ms sampling length

图11采样长度为100 μs的即时采样处理结果

Fig.11Results of immediate sampling processing with 100 μs sampling length

图12采样长度为100 μs的预采样处理结果

Fig.12Results of pre-sampling processing with 100 μs sampling length

表1协方差矩阵不匹配情况分类及影响评估方法

Table1Classification of covariance matrix mismatch and its impact assessment method

图1空间分布式间歇干扰信号示意图

Fig.1Schematic diagram of spatially distributed intermittent interferences

图2阵列接收机信号的处理流程

Fig.2Signal processing flow of array receiver

图3采样长度远大于闪烁周期时采样数据集和处理数据集示意图

Fig.3Sampling data set and processing data set when the sampling length is much larger than the scintillation period

图4采样长度远小于闪烁周期时采样数据集和处理数据集示意图

Fig.4Sampling data set and processing data set when the sampling length is much less than the scintillation period

图5采样不充分时的不同方向信噪比损失

Fig.5Loss of signal-to-noise ratio in different directions with insufficient sampling

图6漏采样时不同方向的信噪比损失

Fig.6Loss of signal-to-noise ratio in different directions with missed sampling

图7过采样时不同方向信噪比损失（场景1）

Fig.7Loss of signal-to-noise ratio in different directions with redundant sampling (scenario 1)

图8过采样时不同方向信噪比损失（场景2）

Fig.8Loss of signal-to-noise ratio in different directions with redundant sampling (scenario 2)

图9采样长度为1 ms的即时采样处理结果

Fig.9Results of immediate sampling processing with 1 ms sampling length

图10采样长度为1 ms的预采样处理结果

Fig.10Results of pre-sampling processing with 1 ms sampling length

图11采样长度为100 μs的即时采样处理结果

Fig.11Results of immediate sampling processing with 100 μs sampling length

图12采样长度为100 μs的预采样处理结果

Fig.12Results of pre-sampling processing with 100 μs sampling length

表1协方差矩阵不匹配情况分类及影响评估方法

Table1Classification of covariance matrix mismatch and its impact assessment method

MIRALLES D, MOGHADAM M S, AKOS D M. GNSS threat monitoring and reporting with the android raw GNSS measurements and STRIKE3[C]//Proceedings of the 32nd International Technical Meeting of the Satellite Division of the Institute of Navigation,2019:275-289.

PATTINSON M, DUMVILLE M, YING Y Q,et al. Standardization of GNSS threat reporting and receiver testing through international knowledge exchange,experimentation and exploitation[STRIKE3][R/OL].(2020-02-07)[2021-08-21].https://www.gps.gov/cgsic/meetings/2017/dumville.pdf.

LIU F L, WANG Y D, YUE L,et al. Jamming method based on optimal power difference for LMS-GPS receiver[J]. Progress in Electromagnetics Research M,2019,77:167-175.

WANG J, LIU W X, CHEN F Q,et al. GNSS array receiver faced with overloaded interferences:anti-jamming performance and the incident directions of interferences[J]. Journal of Systems Engineering and Electronics,2023,34(2):335-341.

KERR T H,Ⅲ. Vulnerability of recent GPS adaptive antenna processing(and all STAP/SLC)to statistically nonstationary jammer threats[C]//Proceedings of International Symposium on Optical Science and Technology,2001.

KONOVALTSEV A, DE LORENZO D S, HORNBOSTEL A,et al. Mitigation of continuous and pulsed radio interference with GNSS antenna arrays[C]//Proceedings of 21st International Technical Meeting of the Satellite Division of the Institute of Navigation,2008.

朱林, 方胜良, 吴付祥. 分布式闪烁干扰下卫星数据链效能评估与仿真[J]. 信息工程大学学报,2014,15(6):697-701,713. ZHU L, FANG S L, WU F X. Evaluation and simulation of the efficiency of satellite data links in the presence of distributing blinking jamming[J]. Journal of Information Engineering University,2014,15(6):697-701,713.(in Chinese)

SIDDIQUI F A, SRENG V, DANILO-LEMOINE F,et al. Suppression of intermittent interference using smart antenna with distributed training scheme[C]//Proceedings of 11th International Bhurban Conference on Applied Sciences & Technology Islamabad,2014:425-429.

BESSON O, STOICA P, KAMIYA Y. Direction finding in the presence of an intermittent interference[J]. IEEE Transactions on Signal Processing,2002,50(7):1554-1564.

SIDDIQUI F A. Suppression of intermittent interference using smart antennas[D]. Canada: Carleton University,2003.

MCDONALD K F, BLUM R S. Exact performance of STAP algorithms with mismatched steering and clutter statistics[J]. IEEE Transactions on Signal Processing,2000,48(10):2750-2763.

BLUM R S, MCDONALD K F. Analysis of STAP algorithms for cases with mismatched steering and clutter statistics[J]. IEEE Transactions on Signal Processing,2000,48(2):301-310.

BESSON O. Analysis of the SNR loss distribution with covariance mismatched training samples[J]. IEEE Transactions on Signal Processing,2020,68:5759-5768.

RAGHAVAN R S. False alarm analysis of the AMF algorithm for mismatched training[J]. IEEE Transactions on Signal Processing,2018,67(1):83-96.

BESSON O. Detection in the presence of surprise or undernulled interference[J]. IEEE Signal Processing Letters,2007,14(5):352-354.

RABIDEAU D J, STEINHARDT A O. Improving the performance of adaptive arrays in nonstationary environments through data-adaptive training[C]//Proceedings of Conference Record of the Thirtieth Asilomar Conference on Signals, Systems and Computers,2002:75-79.

AKDEMIR Ş B, CANDAN Ç. Study of angle of arrival estimation performance of MVDR and Capon methods under intermittent interference[C]//Proceedings of 26th Signal Processing and Communications Applications Conference(SIU),2018.

WANG J, OU G, LIU W X,et al. Performance analysis of SMI filter for antenna array receiver in pulse interference environment[J]. Frontiers in Physics,2023,11:1125431.

WANG J, LIU W X, OU G,et al. Channel scintillations of array global navigation satellite system receiver under distributed intermittent interferences[J]. IET Radar, Sonar & Navigation,2023,17(2):227-235.

Home

1 信号模型和阵列接收机处理流程

2 分布式间歇干扰影响评估

3 自适应处理器性能理论分析

4 数值分析和仿真验证

5 讨论

6 结论