徐文璞，陈奇东，黄晓，刘睿，刘杨

(中国电波传播研究所,山东 青岛 266107)

0 引言

GNSS 经过多年的发展，被广泛应用在民用及军用领域，由于其采用的技术手段和工作原理逐渐为人知悉，近年来对GNSS 导航对抗的研究被广泛关注.在军用领域普遍采用具备自适应抗干扰调零天线的接收机，以提高GNSS 抗干扰能力.针对GNSS 导航对抗场景，自适应抗干扰天线的抗干扰算法已有很多研究，在更好的硬件平台支撑下，抗干扰算法不断完善与优化，抗干扰性能更好[1-2].而针对抗干扰天线的干扰方法的研究相对欠缺，采用分布式布站的方式，作为自适应抗干扰天线的干扰策略较为普遍[3-6]，文献[7]研究了GNSS 信号体制本身固有的抗干扰性能.

基于自适应调零天线的抗干扰算法是以非平稳干扰信号为基础进行设计，因此提出采用非平稳信号进行GNSS 干扰，并对多种非平稳干扰信号与传统宽带干扰的干扰效果进行仿真，验证了采用非平稳干扰信号对抗GNSS 抗干扰天线的有效性.

1 针对自适应抗干扰调零天线的干扰策略

GNSS 的干扰以及抗干扰不断进步[8-13]，最新的干扰及抗干扰算法逐渐装备于军队中，在美军的GPS 现代化规划中，已将自适应调零天线应用于GPS 接收机中，从而提高接收机干扰容限.作为一种显著提高GPS 接收机抗干扰能力的手段，自适应调零天线可以对接收到的干扰信号形成空间零陷，从而最大限度地抑制该方向的干扰.

自适应调零天线通过对多个阵元设计特定的布局，对多个阵元分别设计信号处理通道，通过对多路信号赋权值，利用滤波的方式滤除干扰信号，获取期望信号.自适应调零天线的组成框如图1 所示.

图1 自适应调零天线的组成框

在进行信号的分析处理时，可根据信号的平稳特征，将信号分为平稳信号以及非平稳信号.对于平稳信号时间序列，其均值和方差是与时间无关的常数，协方差与时间无关，而与时间间隔有关.利用样本时间序列所得到的拟合曲线，可在未来的一段时间内反映时间序列的均值、方差、协方差等信息.对于非平稳信号时间序列，样本时间序列的本质特质只存在于当前发生的时间段内，其均值、方差、协方差不是常数，通过对该样本时间序列进行拟合的曲线，不足以反映该信号在未来时间段的特征信息.

对于统计特性随时间呈周期变化的信号，称为循环平稳随机过程，根据该信号的具体特征，可将该随机过程细分为严格平稳和广义平稳两种.

1)严格平稳过程

严格平稳随机过程的典型特点为：该信号的n个随机变量以及其对应的的联 合分布 函数的 时延均相同.

2)广义平稳过程

当某个随机信号 {X(t),t∈T}，满足以下特点时，该信号为广义平稳随机过程.

式中，R为相关计算函数.

对该信号进行自相关计算，相应的自相关函数对于时延 τ=0 共轭对称，可表示为

将该信号输入到线性时不变系统，输出仍是平稳随机信号，对应的时频函数可表示为：

式中：Y表示时频响应；H(f) 是冲激响应h(t) 的傅里叶变换.

目前绝大多数调零天线的抗干扰算法，如空域、空时以及空频自适应处理算法，都是基于干扰信号广义平稳这一假设，即干扰信号的均值为常数，其协方差矩阵与时间起点无关，只与时间间隔有关.因此，如果采用非平稳干扰信号体制，其均值或协方差是时变的，通过单纯的频率分析(傅里叶变换)不具备较好的时频分辨能力，也无法精确地对干扰信号进行定位.此外，干扰信号的非平稳特性打破了天线阵列抗干扰算法设计的前提(干扰信号广义平稳)，对此类抗干扰算法的性能产生较大影响.因此，可采用非平稳干扰制式，对此条件下的干扰效果进行分析评估，验证非平稳干扰信号对GNSS 抗干扰天线的干扰效果.

2 非平稳干扰制式的仿真及干扰效果验证

空频抗干扰方法因其抗干扰性能较优，抗干扰算法的计算复杂度低，目前在军事和民用领域使用最为广泛，因此本仿真选择的GNSS 调零天线的抗干扰方式为空频抗干扰方法.通过仿真维纳过程、分段平稳信号、白噪声(white noise，WN)中时变正弦组合、均值具有趋向性的非平稳随机信号、周期平稳随机信号等非平稳信号的干扰效果，验证非平稳干扰信号应用于GNSS 干扰策略的可行性.

2.1 仿真条件

1)目前性能抗干扰性能最优的自适应抗干扰调零天线的阵元数量为7个，因此本仿真选择的阵列排布方式为七阵元圆阵，阵列排布如图2 所示.

图2 七阵元抗干扰调零天线阵列分布

2)仿真的卫星方向及干扰方向如图3 所示，卫星方向设置为方位角40°，俯仰角50°，扩频干扰和非平稳干扰方向均设置为方位角50°，俯仰角40°.

图3 七阵元抗干扰调零天线阵列分布

3)仿真的干扰信号频率覆盖1 100 MHz～1 600 MHz，该频段覆盖GNSS 的各卫星导航系统，从而验证干扰制式对GNSS 不同频率的干扰效果.

4)空频抗干扰方法对传统宽带扩频信号在干扰方向(方位角50°，俯仰角40°)的衰减增益随频率变化的情况如图4 所示.

图4 扩频信号干扰效果图

2.2 维纳过程干扰效果验证

将WN 输入到一个积分器，其输出即为维纳过程.维纳过程作为一种非平稳正态随机过程，若用x(t) 来表示它，有

式中：rxx为求自 相关函数；α 为输入WN 的方差.若 ω(t) 为平稳正态WN，则维纳过程x(t)为

用维纳过程产生的非平稳信号进行仿真分析，其时域波形如图5 所示.

图5 维纳过程非平稳信号时域仿真

七阵元抗干扰调零天线对维纳过程非平稳信号干扰的衰减增益随频率的变化情况如图6 所示.

图6 维纳过程非平稳信号干扰效果

通过与扩频信号干扰效果类比分析，可以非常明显的观察到采用维纳过程作为干扰信号，对空频抗干扰的干扰效果比扩频干扰的干扰效果好了10 dB 以上，说明空频抗干扰对维纳过程非平稳干扰抑制能力较差.

2.3 分段平稳信号干扰效果验证

对于分段平稳随机信号，在完整的观测周期中，该信号的统计特征为非平稳的，但是在多个观测时刻，该信号会有统计特征的突变，如果将信号以突变时刻为分割点，切分为多个时间段进行分析，每段时间内信号的统计特征则保持着平稳特性.

对各段平稳随机信号可采用不同的平稳随机信号模型来描述其统计特性，通常采用自回归模型(autoregressive model，AR)进行信号仿真，将WN 输入到线性时不变稳定系统，则输出即为分段平稳随机信号.

将AR(P)模型的x(n) 表示为

式中，ω(n) 为零均值与方差的正态WN.

分别对p=2、p=3 以及它们的组合信号进行仿真分析，得到的分段平稳信号如图7 所示.

图7 分段平稳信号时域仿真

七阵元抗干扰调零天线对分段平稳信号干扰的衰减增益随频率的变化情况如图8 所示.

图8 分段平稳信号干扰效果

通过与扩频信号干扰效果类比分析，采用分段平稳信号作为干扰信号，对空频抗干扰的干扰效果比扩频干扰的干扰效果低了15 dB 以上，说明空频抗干扰对分段平稳信号的抑制能力较好.

2.4 时变正弦信号干扰效果验证

时变正弦组合是指恒定或时变幅值与时变频率的正弦波，其中恒定幅值与线性时变频率的正弦波为通常的线性调频(linear frequency modulation，LMF)信号.

2.4.1 固定幅值与时变频率的正弦波干扰仿真

为研究方便起见，令固定幅值为1 与取样间隔Δt=1，设时变频率的正弦波为

经过推导后可得：

对固定幅值与时变频率的正弦波模型进行仿真分析，固定幅值与时变频率的正弦波模型的非平稳信号如图9 所示.

图9 固定幅值与时变频率的正弦波信号时域仿真

七阵元抗干扰调零天线对固定幅值与时变频率的正弦波信号干扰的衰减增益随频率的变化情况如图10 所示.

图10 固定幅值与时变频率的正弦波信号干扰效果

通过与扩频信号干扰效果类比分析，采用固定幅值与时变频率的正弦波信号作为干扰信号，在大部分频段对空频抗干扰的干扰效果比扩频干扰的干扰效果提升约15 dB，但是存在少部分频段与扩频干扰有相似的干扰效果.

进一步考虑固定幅值与时变频率的正弦波信号间断出现的干扰方式，仿真结果图如图11～12 所示.

图11 固定幅值与时变频率的间断正弦波信号时域仿真

图12 固定幅值与时变频率的间断正弦波信号干扰效果

可以发现固定幅值与时变频率的正弦波间断信号比干扰信号连续出现时的干扰能力更好一点，但还是在部分频段存在被抗干扰的风险.

2.4.2 时变幅值与时变频率正弦波干扰仿真

设具有时变幅值与时变频率的正弦波为

对时变幅值与时变频率的正弦波模型进行仿真分析，时变幅值与时变频率的正弦波模型的非平稳信号如图13 所示.

图13 时变幅值与时变频率的正弦波信号时域仿真

七阵元抗干扰调零天线对时变幅值与时变频率的正弦波信号干扰的衰减增益随频率的变化情况如图14 所示.

图14 时变幅值与时变频率的正弦波信号干扰效果

通过与扩频信号干扰效果类比分析，采用时变幅值与时变频率的正弦波信号作为干扰信号，对空频抗干扰的干扰效果比扩频干扰的干扰效果提升约15 dB，说明空频抗干扰对时变幅值与时变频率的正弦波信号的抑制能力较差.

2.5 均值具有趋向性的非平稳随机信号干扰效果验证

均值具有趋向性的非平稳随机信号的特点是它的均值是在某一个趋向性曲线的周围随机波动的.该信号为线性时不变系统，当输入为WN时，输出为均值具有趋向性的非平稳信号.

设线性时不变离散系统的脉冲响应为h(k)，它有n个极点，若输入为WN时，其输出为

对均值具有趋向性的非平稳随机信号模型进行仿真分析，固定幅值与时变频率的正弦波模型的非平稳信号如图15 所示.

图15 均值具有趋向性的非平稳随机信号时域仿真

七阵元抗干扰调零天线对均值具有趋向性的非平稳随机信号干扰的衰减增益随频率的变化情况如图16 所示.

图16 均值具有趋向性的非平稳随机信号干扰效果图

通过与扩频信号干扰效果类比分析，采用均值具有趋向性的非平稳随机信作为干扰信号，对空频抗干扰的干扰效果比扩频干扰的干扰效果提升15 dB，说明空频抗干扰对均值具有趋向性的非平稳随机信号非平稳干扰抑制能力很差.

2.6 周期平稳随机信号干扰效果验证

周期平稳随机信号的典型特征为：其统计特征随时间周期性变化或近似周期性变化，该信号作为特殊的非平稳随机信号，可通过线性周期时变系统产生.

对于2 阶非平稳信号 {x(t)(t∈(-∞,∞))}，其均值和自相关函数均呈现以T 为周期的周期性变化，即

式中：mx(t)为均值函数；rxx(t)为自相关函数.

对该周期平稳随机信号进行仿真建模，周期平稳随机模型的非平稳信号如图17 所示.

图17 周期平稳随机信号时域仿真

七阵元抗干扰调零天线对周期平稳随机信号干扰的衰减增益随频率的变化情况如图18 所示.

图18 周期平稳随机信号干扰效果

通过与扩频信号干扰效果进行类比分析，采用均值具有趋向性的非平稳随机信作为干扰信号，对空频抗干扰的干扰效果与扩频干扰的干扰效果相似，只在部分频段比扩频干扰的干扰效果好.

3 结束语

通过对多种典型非平稳干扰信号进行仿真，并对采用空频抗干扰方法的GNSS 七阵元调零天线进行干扰效果验证.对比非平稳干扰信号与扩频信号的干扰性能，优选出比扩频信号干扰效果更好的非平稳干扰信号如表1 所示.

表1 非平稳信号干扰效果分析

由表1 可知，时变幅值与时变频率的正弦波信号，在仿真的GNSS 频段中，都表现出较好的干扰性能.因此，采用时变幅值与时变频率的正弦波作为干扰信号制式，可以有效对抗GNSS 抗干扰天线.