GNSS接收机抗干扰测试与MUSIC算法实践

1. 项目背景与核心需求

在卫星导航系统测试领域，GNSS模拟器已成为接收机性能验证的关键工具。某高校实验室近期利用GNSS模拟器开展接收机抗干扰测试，这项工作的核心目标是评估接收机在复杂电磁环境下的信号捕获与跟踪能力。不同于常规的静态测试，抗干扰测试需要模拟真实场景中可能遇到的窄带干扰、宽带噪声干扰以及多径效应等多种干扰形态。

从专利文献CN103630910A可以看出，传统功率倒置算法在干扰源附近会产生带状零陷分布，而采用MUSIC空间谱估计算法的新型抗干扰方法能实现更精确的干扰方向定位。这正好契合高校实验室对测试精度的严苛要求——需要在信干比降至20dB时仍能保持有效干扰抑制。

2. 测试系统架构设计

2.1 硬件组成方案

测试系统采用三级架构设计：

1. 信号生成层：使用多通道GNSS模拟器（如Spirent GSS9000），支持GPS L1/L5、Galileo E1/E5等多频点信号生成
2. 干扰注入层：通过矢量信号发生器（如Keysight MXG N5183B）产生：
   • 单音连续波干扰（频率精度±0.01ppm）
   • 线性调频干扰（带宽设置1-10MHz可调）
   • 宽带噪声干扰（噪声密度-110dBm/Hz起）
3. 被测设备层：高校自研的八阵元圆阵接收机，配备FPGA实时处理单元

关键配置技巧：干扰信号功率建议采用阶梯递增方式，从-130dBm开始以5dB为步进，可准确捕捉接收机抗干扰门限。

2.3 空间谱估计实现

基于MUSIC算法的DOA估计流程如下：

matlab复制% 阵列接收数据建模
N = 8; % 阵元数
theta = [-30, 20, 45]; % 真实干扰来向
SNR = 20; % 信噪比(dB)
X = zeros(N, 1000);
for k = 1:length(theta)
    a = exp(-1j*2*pi*(0:N-1)'*sind(theta(k))/2); % 导向矢量
    X = X + a*(randn(1,1000)+1j*randn(1,1000));
end
X = X + 10^(-SNR/20)*(randn(N,1000)+1j*randn(N,1000));

% MUSIC算法