精确时延补偿的自适应旁瓣相消射电抗干扰方法研究*1

黄　达，王　壮，程　翥，董　亮，苏　瑶

(1. 国防科技大学电子科学与工程学院ATR国防科技重点实验室，湖南 长沙　410073；2. 国防科技大学电子科学与工程学院电子工程所，湖南 长沙　410073；3. 中国科学院云南天文台，云南 昆明　650011；4. 西安电子科技大学，陕西 西安　710071)



精确时延补偿的自适应旁瓣相消射电抗干扰方法研究*1

黄达1，王壮1，程翥2，董亮3，苏瑶4

(1. 国防科技大学电子科学与工程学院ATR国防科技重点实验室，湖南 长沙410073；2. 国防科技大学电子科学与工程学院电子工程所，湖南 长沙410073；3. 中国科学院云南天文台，云南 昆明650011；4. 西安电子科技大学，陕西 西安710071)

摘要：导航信号作为射电干扰信号影响L波段射电观测时，与雷达、通信领域的干扰不同，其功率通常小于噪声功率，使得利用干扰信号相关性的传统抗干扰方法效果不佳，且缺乏有效的评价指标。针对以上问题，依据射电信号处理流程构建了射电抗干扰系统框架，并从射电观测数据应用出发，给出了干扰抑制度和信号损失度两个指标。同时提出在精确补偿时延后利用抑制干扰子空间噪声分量的自适应旁瓣相消方法来消除导航信号的影响；仿真结果表明，本方法能有效抑制导航信号对射电观测的影响，且性能优于传统基于特征子空间的自适应旁瓣相消方法。

关键词：射电天文抗干扰；自适应旁瓣相消；导航信号；可变分数时延滤波器

随着射电天文的发展，国际通讯联盟所划分的射电观测频段已不能满足射电天文研究的需求，观测保护频段以外的射电源，常常受到干扰信号的影响，极大限制了射电源的观测质量。导航信号是L频段射电观测的常见干扰，根据全球定位系统和北斗导航卫星的ICD文件可知①-②，全球定位系统导航信号的主要频率是L1：1 575.42 ± 1.023 MHz、L2：1 227.6 ± 10.23 MHz等，北斗导航信号的主要频率是：B1：1 556.098 ± 2.046 MHz、B3：1 268.52 ± 10.23 MHz等，而这些频率与红移的分子谱线、脉冲星等射电源的观测频段重合。同时，根据ICD文件中导航信号的最小电平可以换算其功率通量谱密度为-196～-201 dBw/(m2/Hz)，而在同样带宽下1 mJy的射电信号功率通量谱密度为-290 dBw/(m2/Hz)，因此在接收时导航信号比射电信号强89～91 dB。此外，天线接收系统对接收信号也有影响，由于导航信号往往从旁瓣进入射电接收机，经过放大、滤波等过程其功率通常小于噪声功率但大于射电信号功率。

通过以上分析可知，如何抑制射电观测中弱于噪声的干扰信号是本文讨论的主要问题与难点。对于导航信号这类弱于噪声的干扰信号而言，完全获取信号的先验信息难度较大，但可利用信号本身的特点来抗干扰，如空间特性和相关特性。同时，抗干扰方法的实施难免对射电信号处理有一定影响，因此抗干扰系统的建立既需要有效利用干扰信号本身的特点，同时也需尽量避免对射电信号处理的影响。进一步来说，有效地利用干扰信号的空间特性需要基于多通道接收的方式，文[1-2]利用干涉仪和相控阵天线自身多通道的特性提取干扰信号的子空间，通过正交投影的方式消除干扰信号，然而由于信号空间与干扰空间不垂直，导致对观测信号影响较大。为了利用干扰信号的相关性，借助辅助天线是常用方式。文[3-4]通过辅助天线估计更为精准的干扰子空间，并运用如正交子空间、互子空间投影和多阵列旁瓣相消等方式抑制干扰信号。对于大型单射电天线而言，要利用干扰信号的空间特性和相关性，则需要借助多个辅助天线。文[5]从雷达信号处理角度出发，基于特征子空间的自适应旁瓣相消方法，利用辅助阵列中强干扰信号的特征子空间与噪声子空间良好的垂直性，消除了接收信号中的强干扰信号。对于抑制射电观测中弱于噪声的导航信号而言，需根据射电信号处理流程重新构建抗干扰系统框架，并且改进自适应旁瓣相消方法有效抑制导航信号。同时考虑到射电观测数据应用，应给出新的评价指标来评估干扰抑制程度和信号损失程度。此外，实际情况中射电天线(主天线)口径通常较大，使架设辅助阵列离主天线较远，因此由于空间位置产生的信号时延导致导航信号相关性减少，需要考虑时延补偿的问题。

综上所述，本文首先依据射电信号处理流程构建了基于辅助天线的射电抗干扰系统框架，并从射电观测数据应用出发，给出了干扰抑制度和信号损失度两个指标。同时本文在精确补偿主天线接收信号时延后，利用抑制干扰子空间噪声分量的自适应旁瓣相消方法来抑制导航信号的影响。最后利用蒙特卡洛仿真方式，讨论了不同主天线和辅助阵列间距(即主、辅天线距离)、干扰波达角度、观测时间、辅助阵元口径和阵元数量对抗干扰性能的影响。

1基于辅助阵列的射电抗干扰系统框架构建

在大型单射电天线中，利用导航信号的空间特性和相关性抑制干扰，需要借助多通道的辅助天线实现，即辅助阵列。自适应旁瓣相消作为常用的抗干扰方法之一，本质是利用辅助阵列与主天线中干扰信号的相关性和干扰的空间特性获得辅助阵列的最优权矢量，使合成的天线方向图在干扰方向上产生零点，达到干扰抑制的目的。此外，为了降低抗干扰系统对射电信号处理的影响，抗干扰系统应独立于射电信号处理流程之外。根据以上分析，建立基于辅助阵列的射电抗干扰系统框架，如图1。

图1　基于辅助天线的射电抗干扰系统框架

(1)指标构建：与雷达、通信领域不同，射电天文抗干扰中干扰信号通常小于噪声信号。因此有必要着眼于射电天文数据应用提出评价指标来评估干扰消除程度和信号损失程度。

(2)时延补偿：由于射电天文望远镜口径通常较大，利用辅助天线抗干扰时，其架设位置离射电天线(主天线)位置较远，因此需要精确补偿信号时延，提高导航信号的相关性。

(3)自适应旁瓣相消：完全已知所有导航卫星的先验信息难度较大，但其运行轨道相对固定易于预测，因此可利用导航信号良好的空间特性消除影响。自适应旁瓣相消方法中基于特征子空间类的方法可以在干扰较强时，有效避免估计误差的影响，获得辅助阵列的最优权矢量。然而导航信号功率通常小于噪声功率，因此需要改进此类方法。

2精确时延补偿的自适应旁瓣相消射电抗干扰方法

2.1评价指标

在射电抗干扰中，干扰信号功率通常小于噪声功率，传统的性能指标如相消比等，无法有效地刻画干扰和信号分量的变化程度。为解决此问题，本文从射电观测数据应用出发，以辅助天线接收信号以及先验的射电信号为标尺，分别衡量了抗干扰前后主天线接收信号中干扰信号与射电信号的相关程度，构建了干扰抑制度GI和信号损失度GS表征干扰抑制程度与信号损失程度，其定义分别为：

干扰抑制度GI：

(1)

同理，信号损失度GS：

(2)

2.2时延补偿

为了获取良好的滤波器通带特性，本文采用基于非迭代的加权最小二乘的可变分数时延全通滤波器实现分数时延补偿[8]，滤波器结构基于farrow结构[9]，该结构基于不同准则逼近理想的滤波器获得滤波器系数后，利用高阶分数时延的多项式近似得到滤波器系数。此方式即可将分数时延与滤波器系数分离，使整个滤波器的传递函数简化为分数时延对多个子滤波器的加权和，因此该结构大大减少了求解滤波器系数的复杂度。加权最小二乘滤波方法求解滤波器系数的主要思想是在一定通带范围和时延因子(分数时延与采样频率的乘积)范围内，使可变分数时延滤波器的频率响应与理想的分数时延频率响应函数误差最小，即

(3)

其中，B是由滤波器系数组成的矩阵；W(ω,p)是非负的权值函数；eH为可变分数时延滤波器与理想分数时延频率响应的误差函数；p为除去整数部分的分数时延因子。不难发现(3)式中存在积分过程，利用泰勒级数展开可简化积分过程，即通过非迭代的方式，快速获得滤波器系数。图2(a)和图2(b)分别是分数时延因子p变化范围为-0.5～0.5，归一化频率ω变化范围为0～0.9π时，基于加权最小二乘滤波的可变分数时延全通滤波器输出的相位时延和与理想的分数时延频率响应函数误差结果。

图2WLS可变分数时延滤波器结果

Fig.2Results of the WLS-VFD filter

通过结果可得，基于非迭代加权最小二乘滤波的可变分数时延滤波器能在宽频带范围下可实现精确的分数时延。

2.3自适应旁瓣相消

自适应旁瓣相消方法的核心问题是求解最优权值矢量。基于特征子空间的自适应旁瓣相消方法是利用干扰子空间与噪声子空间的正交特性，将权值矢量向干扰子空间投影得到最优权值矢量，可见准确估计干扰子空间十分重要。辅助阵列的接收信号主要是干扰信号和噪声信号，可通过估计辅助阵列的协方差矩阵进而获得干扰子空间和噪声子空间。假设存在p个互不相关的窄带干扰，协方差矩阵R^x特征分解如下式：

(4)

其中，R^x是由k(k≥p)个辅助阵元数据长度为N的接收信号估计的协方差矩阵；DI和DN是分别由前p个大特征值和后k-p个小特征值构成的对角阵；UI和UN分别是DI和DN对应特征向量构成的矩阵，其列矢量分别张成的空间为干扰子空间和噪声子空间。

(1)首先根据导航信号入射方向(即波达角)和主天线与辅助阵列参考阵元间的空间位置估计信号时延τ；

(2)根据信号时延τ和采样频率计算时延因子I+p，通过补零对整数时延I补偿，通过2.2中的可变分数时延滤波器对分数时延p补偿；

(3)通过有限长度的快拍数估计辅助阵列的协方差矩阵和主天线与辅助阵列的相关矢量，即R^x、R^xy；

(5)

(6)

通过以上方式可将干扰子空间的噪声分量消除一部分，使干扰子空间与噪声子空间有更好的正交性，进而获得更优的抗干扰性能。

3仿真结果分析

3.1参数设置

本文以FAST 30 m模型脉冲星天线为主天线模型[10]。观测不同的射电信号，其信号形式有较大差异，无法用统一的形式表示射电信号，因此本文用功率谱流量密度为1 mJy的随机信号代替射电信号；辅助天线是由无方向性阵元构成的均匀线阵，且阵元间距设置为0.09 m。以调制C/A码B3频点(1 268.52 MHz)的北斗导航信号为干扰信号。

为了定量地对比时延补偿对本文方法的影响和本方法较传统的基于特征子空间自适应旁瓣相消方法性能的影响，分别对比了在时延补偿和时延未补偿情况下的本文方法与基于特征子空间的自适应旁瓣相消方法，并讨论了不同主天线和辅助阵列间距(即主、辅天线距离)、干扰波达角度、观测时间、辅助阵元口径和阵元数量对抗干扰性能的影响，参数设置如表1。此外，本文利用蒙特卡洛模拟方式消除随机误差对算法评价可靠性的影响。根据大数定理可知，在置信度α=0.03，方差为0.074 3(辅助天线噪声标准差)，蒙特卡洛次数为200次时，蒙特卡洛误差为ε=0.015 8，因此蒙特卡洛次数为200次时可以满足比较算法性能的要求。

表1　仿真参数对比

3.2结果分析

3.2.1主天线与辅助阵列距离的影响

主天线与辅助阵列中参考阵元的距离远近导致接收信号的时延不同进而影响导航信号的相关性。因此本文讨论了主、辅天线距离对方法性能的影响，如图3。

结合表1和图3(a)可知，随着主天线与辅助阵列距离的增加，时延补偿的本文方法干扰抑制度比时延补偿和未补偿的自适应旁瓣相消方法高15 dB左右，同时，时延未补偿的本文方法干扰抑制度，由于时延的影响随主、辅天线距离增加而下降；另外，无论时延是否补偿，本文方法和基于特征子空间的自适应旁瓣相消方法，随着主天线与辅助阵列距离增加，信号损失度相差不大，约1 dB，如图3(b)。

图3不同主天线与辅助阵列距离的方法性能对比

Fig.3Performance comparison with different distance between the main-antenna and the auxiliary antennas as the variable

3.2.2波达方向和阵元口径的影响

不同来向的导航信号由于方向图的影响，功率差异很大，导致干扰信号(导航信号)在主天线中干噪比差异也很大；另一方面，阵元口径影响干扰信号的接收功率。因此本文讨论了不同信号波达角θ(与主瓣轴的夹角)和阵元口径对方法性能的影响，如图4。

图4不同波达角和辅助阵元口径的方法性能对比

Fig.4Performance comparison with direction of arrival as the variable and aperture of the auxiliary antennas as the variable

结合表1和图4(a)可得，随着导航信号波达角远离主瓣，其信号功率也随方向图变化而减小。时延补偿和时延未补偿的本文方法干扰抑制度也随之减小，但是时延补偿的本文方法比自适应旁瓣相消方法干扰抑制度高7～20 dB，而时延未补偿的本文方法仅在信号波达角小于2.5°时比自适应旁瓣相消方法干扰抑制度高；图4(b)中导航信号波达角在1.5°范围内信号损失度比旁瓣方向上更大，是因为合成的方向图使主瓣零陷，零陷程度越深，信号的损失越大，因此无论时延是否补偿，本文方法信号损失度比自适应旁瓣相消方法高约3 dB。当导航信号从旁瓣甚至远旁瓣进入接收机时，旁瓣出现零陷，因此信号损失度GS逐渐趋于0。

结合表1和图4(c)可得，随着辅助阵元口径增加(等价于干噪比的增加)，时延补偿的本文方法和自适应旁瓣相消方法干扰抑制度GI均不断提高，由于口径的增大使得干扰信号接收功率增大。时延补偿的自适应旁瓣相消方法在辅助阵元口径大于1.4 m时，干扰抑制度超过时延未补偿的本文方法，但其干扰抑制度比时延补偿的本文方法的GI低7～10 dB。随着阵元口径的增加，本文方法与自适应旁瓣相消方法的信号损失度相差2 dB内，如图4(d)。

3.2.3观测时间和阵元数量的影响

观测时间是指参与估计R^x、R^xy运算的快拍数量(即接收信号数据长度)N所对应的时间，不同的观测时间影响相关矩阵的估计精度，而阵元数量影响合成方向图在干扰来向零点的零陷程度和可消除干扰个数，因此本文讨论了不同观测时间和阵元数量对方法性能的影响，如图5。

结合表1和图5(a)可得，随着观测时间的增加，时延补偿的本文方法干扰抑制性能不断提高，干扰抑制度比自适应旁瓣相消方法高10～18 dB，且比时延未补偿的本文方法干扰抑制度高10～12 dB。相比自适应旁瓣相消方法，时延补偿的本文方法的信号损失度GS高约1 dB，如图5(b)。

图5不同观测时间和阵元数量的方法性能对比

Fig.5Performance comparison between different observation time and between different numbers of auxiliary antenna

结合表1和图5(c)可知，随着阵元数的增加，两种抗干扰方法的干扰抑制度都会随之增加，但时延补偿的本文方法比自适应旁瓣相消方法干扰抑制度GI高约14 dB，且比时延未补偿的本方法干扰抑制度高约4 dB。此外，无论时延是否补偿，本文方法的信号损失度GS整体比自适应旁瓣相消方法高约1 dB，如图5(d)。

4结论

本文针对导航信号干扰L波段射电观测这一实际情况，构建了基于辅助阵列的抗干扰系统框架，并从射电数据应用出发，提出了干扰抑制度和信号损失度两个指标。利用导航信号良好的空间特性和相关特性，改进了传统的基于特征子空间的自适应旁瓣相消方法，使其在干扰功率小于噪声功率时依然可有效地消除干扰。此外，本文考虑到由于射电天线与辅助阵列空间相对位置造成的时延导致导航信号相关性降低问题。利用基于非迭代的加权最小二乘的可变分数时延全通滤波器实现了时延的精确补偿，提高了导航信号的相关性。最后基于所提指标，通过蒙特卡洛仿真方式，对比了本文方法与自适应旁瓣相消方法，结果表明本文方法的干扰抑制度提高了10～20 dB，但信号损失度增加约1 dB。

参考文献：

[1]Hellbourg G, Trainini T, Weber R, et al. RFI subspace estimation techniques for new generation radio telescopes[C]// Signal Processing Conference (EUSIPCO), 2012 Proceedings of the 20th European. 2012: 200-204.

[2]Sardarabadi A M, van der Veen A J. Subspace estimation using factor analysis[C]// Sensor Array and Multichannel Signal Processing Workshop (SAM),2012 IEEE 7th.2012:477-480.

[3]Hellbourg G, Chippendale A P, Kesteven M J, et al. Reference antenna-based subspace tracking for RFI mitigation in radio astronomy[C] // 2014 IEEE Global Conference on Signal and Information Processing (GlobalSIP). 2014: 1286-1290.

[4]Jeffs B D, Li L, Warnick K F. Auxiliary antenna-assisted interference mitigation for radio astronomy arrays[J]. IEEE Transactions on Signal Processing, 2005, 53(2): 439-451.

[5] 赵永波, 王志慧, 张守宏. 一种基于特征空间的自适应天线旁瓣相消算法[J]. 航空计算技术, 2004, 34(1): 5-8.

Zhao Yongbo, Wang Zhihui, Zhang Shouhong. An eigenspace-based algorithm for adaptive antenna sidelobe cancellation[J]. Aeronautical Computer Technique, 2004, 34(1): 5-8.

[6]秦顺友, 李光, 张立军, 等. 65米射电望远镜天线测量技术[C]// 中国电子科技集团公司第五十四研究所. 2013年全国天线年会. 2013: 1-61.

[7]张洪波, 毛佩锋, 汪敏, 等. 40 m口径射电望远镜[J]. 天文研究与技术——国家天文台台刊, 2008, 5(2): 187-191.

Zhang Hongbo, Mao Peifeng, Wang Min, et al. The 40m radio telescope[J]. Astronomical Research & Technology——Publications of National Astronomical Observatories of China, 2008, 5(2): 187-191.

[8]Deng Tianbo. Noniterative WLS design of allpass variable fractional-delay digital filters[J]. IEEE Transactions on Circuits and Systems I: Regular Papers, 2006, 53(2): 358-371.

[9]Farrow C W. A continuously variable digital delay element[J]. IEEE International Symposium on Circuits and Systems, 1988, 3: 2641-2645.

[10] 高健健, 曹洋, 甘恒谦, 等. FAST 30 m模型L-波段馈源的设计、制作和性能测试[J]. 天文研究与技术——国家天文台台刊, 2007, 4(1): 48-52.

Gao Jianjian, Cao Yang, Gan Henqian, et al. Design、manufacture and test of the L-band feed for the FAST 30m demonstrator[J]. Astronomical Research & Technology——Publications of National Astronomical Observatories of China, 2007, 4(1): 48-52.

*基金项目：国家自然科学基金 (11173068, 11303094) 资助.

收稿日期：2015-07-22；

修订日期：2015-08-14

作者简介：黄达，男，硕士. 研究方向：射电天文抗干扰算法. Email: huangda_nudt@sina.com

中图分类号：P161

文献标识码：A

文章编号:1672-7673(2016)03-0284-09

Mitigation of the Radio Frequency Interference Using an Adaptive Side-Lobe Cancellation Based on Accurate Delay Compensation

Huang Da1，Wang Zhuang1，Cheng Zhu2，Dong Liang3，Su Yao4

(1. ATR Key Lab, School of Electronic Science and Engineering, National University of Defense Technology, Changsha 410073, China; 2. Research Institution of Electronic Engineering, School of Electronic Science and Engineering, National University of Defense Technology, Changsha 410073, China; 3. Yunnan Observatories, Chinese Academy of Sciences, Kunming 650011, China; 4. School of Electronic Engineering Xidian University, Xi′an 710071, China)

Abstract:When navigation signals, as radio interferences, affect observation in L-band, they are different from other types of interferences in communication field or radar.Using the correlation interference, traditional methods are ineffective and they lack effective evaluation indexes, too, because the power of navigation signal is normally weaker than noise power. To solve the above problems, this paper constructs a system framework of radio astronomy interference based on radio signal processing, and raises two indexes: interference suppression and signal loss based on the application of the observation data. After accurate delay compensation, this paper proposes a method to eliminate the effects of the navigation signals, using the adaptive side-lobe cancellation that suppresses the noise component of interference subspace. From the results of the simulation, conclusion can be drawn that the proposed method can effectively suppress the navigation signal in the observations, and the performance is better than that of traditional adaptive side-lobe cancellation based in feature subspace.

Key words:Radio astronomy interference; Adaptive side-lobe cancellation; Navigation signal; Variable fractional delay filter

CN 53-1189/PISSN 1672-7673

①http://www.gps.gov/technical/icwg/ICD-GPS-870A.pdf

②http://www.beidou.gov.cn/attach/2012/12/27/201212273da29c5eb8274deb8cd2b178228ba2bd.pdf