前视多通道SAR自适应鉴别抑制欺骗干扰

孟智超 卢景月 张磊

引用本文:
Citation:

前视多通道SAR自适应鉴别抑制欺骗干扰

    作者简介: 孟智超(1997–),男,陕西宝鸡人,博士生。2018年6月在西安电子科技大学电子工程学院获得学士学位,现为西安电子科技大学电子工程学院博士研究生。主要研究方向为合成孔径雷达干扰与抗干扰;卢景月(1994–),男,河北沧州人,博士生。2016年在中国海洋大学信息科学与工程学院获得学士学位,现为西安电子科技大学电子工程学院博士研究生。主要研究方向为合成孔径雷达成像;张 磊(1984–),男,浙江金华人,博士,副教授。2012年在西安电子科技大学电子工程学院获得博士学位,现担任西安电子科技大学电子工程学院副教授。研究方向为雷达成像、目标识别。E-mail: leizhang@xidian.edu.cn.
    通讯作者: 张磊, leizhang@xidian.edu.cn
  • 基金项目:

    国家自然科学基金(61771372, 61771367),上海市自然科学基金(16ZR1434900)

  • 中图分类号: TN975

Forward-looking Multi-channel SAR Adaptive Identification to Suppress Deception Jamming

    Corresponding author: ZHANG Lei, leizhang@xidian.edu.cn
  • Fund Project: The National Natural Science Foundation of China(61771372, 61771367), Sponsored by Natural Science Foundation of Shanghai (16ZR1434900)

    CLC number: TN975

  • 摘要: 前视多通道SAR成像存在回波左右模糊的问题,在成像时需要利用空域资源进行解模糊处理,使得其成像过程较一般侧视SAR更加复杂。在复杂电磁环境下,若要获得无干扰、不模糊的前视SAR图像,难度较大。该文提出一种基于方位向自适应波束形成(AADBF)技术的欺骗干扰自适应鉴别抑制算法。该算法首先采用AADBF技术对多通道接收回波信号进行地物相消,保留欺骗干扰样本。然后利用门限检测法在高分辨SAR图像上,鉴别出欺骗干扰的像素位置。最后,对存在干扰的像素点进行自适应空域滤波,从而达到抗欺骗干扰的目的。仿真结果表明,该方法能够有效地鉴别并抑制欺骗干扰,与此同时能够实现前视SAR无干扰聚焦成像。
  • 图 1  FLMC-SAR几何模型

    Figure 1.  FLMC-SAR geometric model

    图 2  欺骗干扰与真实回波空时关系差异

    Figure 2.  Difference between deception jamming and real echo space time relationship

    图 3  欺骗干扰抑制算法流程

    Figure 3.  Deception interference suppression algorithm flow

    图 4  加入欺骗干扰前后对比

    Figure 4.  Comparison before and after adding spoofing jamming

    图 5  地物相消与干扰位置定位

    Figure 5.  Location cancellation and interference position locating

    图 6  干扰角度估计结果

    Figure 6.  Interference angle estimation results

    图 7  欺骗干扰抑制解模糊对比

    Figure 7.  Comparison of deception interference suppression ambiguity resolution

    表 1  系统参数

    Table 1.  System parameters

    参数数值
    载频35 GHz
    信号带宽50 MHz
    脉冲宽度0.5 μs
    采样带宽60 MHz
    平台高度4000 m
    脉冲重复频率9 kHz
    波束方位宽度20°
    波束俯仰宽度22°
    场景中心斜距${R_s}$8 km
    阵列实孔径长度0.0642 m
    下载: 导出CSV

    表 2  处理性能对比

    Table 2.  The comparison of processing performance

    处理结果信干比增益(dB)干扰抑制比(dB)信号抑制比(dB)
    含干扰图像–15.6000
    全域干扰抑制030.7530.75
    局部干扰抑制32.0631.75–0.31
    下载: 导出CSV
  • [1] 张英杰, 王彦平, 谭维贤, 等. 机载前视阵列SAR运动补偿研究[J]. 雷达学报, 2013, 2(2): 168–179. doi: 10.3724/SP.J.1300.2013.20067ZHANG Ying-jie, WANG Yan-ping, TAN Wei-xian, et al. Motion compensation for airborne forward looking synthetic aperture radar with linear array antennas[J]. Journal of Radars, 2013, 2(2): 168–179. doi: 10.3724/SP.J.1300.2013.20067
    [2] DAI S L, LIU M, SUN Y J, et al. The latest development of high resolution imaging for forward looking SAR with multiple receiving antennas[C]. Proceedings of Scanning the Present and Resolving the Future. Proceedings. IEEE 2001 International Geoscience and Remote Sensing Symposium, Sydney, NSW, Australia, 2001: 1433–1435.
    [3] 冯德军, 王俊杰, 王俊卿. 移频导前假目标的特性分析及其鉴别方法[J]. 雷达学报, 2017, 6(4): 325–331. doi: 10.12000/JR17026FENG Dejun, WANG Junjie, and WANG Junqing. Signature analysis and discrimination method of preceded frequency-shift false target[J]. Journal of Radars, 2017, 6(4): 325–331. doi: 10.12000/JR17026
    [4] 齐维孔, 禹卫东. 距离向DBF-SAR系统抗欺骗干扰研究[J]. 电子与信息学报, 2010, 32(12): 2830–2835. doi: 10.3724/SP.J.1146.2010.00074QI Wei-kong and YU Wei-dong. Research on countering deceptive jamming with range DBF-SAR system[J]. Journal of Electronics &Information Technology, 2010, 32(12): 2830–2835. doi: 10.3724/SP.J.1146.2010.00074
    [5] 张双喜, 孙光才, 刘艳阳, 等. 基于四通道SAR欺骗式干扰抑制算法[J]. 现代雷达, 2011, 33(2): 22–26, 30. doi: 10.3969/j.issn.1004-7859.2011.02.006ZHANG Shuang-xi, SUN Guang-cai, LIU Yan-yang, et al. Deception interference supression algorithm based on four-channel for SAR[J]. Modern Radar, 2011, 33(2): 22–26, 30. doi: 10.3969/j.issn.1004-7859.2011.02.006
    [6] 冯祥芝, 许小剑. 随机线性调频斜率SAR抗欺骗干扰方法研究[J]. 系统工程与电子技术, 2009, 31(1): 69–73. doi: 10.3321/j.issn:1001-506X.2009.01.018FENG Xiang-zhi and XU Xiao-jian. Study of countermeasures to deceptive jamming using random linear modulation frequency ratio SAR[J]. Systems Engineering and Electronics, 2009, 31(1): 69–73. doi: 10.3321/j.issn:1001-506X.2009.01.018
    [7] KRIEGER G, GEBERT N, and MOREIRA N. Multidimensional waveform encoding: A new digital beamforming technique for synthetic aperture radar remote sensing[J]. IEEE Transactions on Geoscience and Remote Sensing, 2008, 46(1): 31–46. doi: 10.1109/TGRS.2007.905974
    [8] 谢文冲, 段克清, 王永良. 机载雷达空时自适应处理技术研究综述[J]. 雷达学报, 2017, 6(6): 575–586. doi: 10.12000/JR17073XIE Wenchong, DUAN Keqing, and WANG Yongliang. Space time adaptive processing technique for airborne radar: An overview of its development and prospects[J]. Journal of Radars, 2017, 6(6): 575–586. doi: 10.12000/JR17073
    [9] 王永良, 丁前军, 李荣锋. 自适应阵列处理[M]. 北京: 清华大学出版社, 2009: 66–72.WANG Yong-liang, DING Qian-jun, and LI Rong-feng. Adaptive Array Processing[M]. Beijing: Tsinghua University Press, 2009: 66–72.
    [10] 龚耀寰. 自适应滤波-时域自适应滤波和智能天线[M]. 第2版, 北京: 电子工业出版社, 2003: 288–297.GONG Yao-huan. Adaptive Filtering-Time Domain Adaptive Filtering and Smart Antenna[M]. 2nd ed., Beijing: Public House of Electronic Industry, 2003: 288–297.
    [11] YEGULALP A F. Fast backprojection algorithm for synthetic aperture radar[C]. Proceedings of 1999 IEEE Radar Conference, Radar into the Next Millennium, Waltham, Massachusetts, USA, 1999: 60–65. doi: 10.1109/NRC.1999.767270.
    [12] 包敏, 周鹏, 史林. 双天线前视弹载SAR解模糊算法研究[J]. 电子与信息学报, 2013, 35(12): 2857–2862. doi: 10.3724/SP.J.1146.2013.00083BAO Min, ZHOU Peng, and SHI Lin. Study on deambiguity algorithm for double antenna forward looking missile borne SAR[J]. Journal of Electronics &Information Technology, 2013, 35(12): 2857–2862. doi: 10.3724/SP.J.1146.2013.00083
    [13] 卢景月, 张磊, 王冠勇. 前视多通道合成孔径雷达解模糊成像方法[J]. 电子与信息学报, 2018, 40(12): 2820–2825. doi: 10.11999/JEIT180177LU Jing-yue, ZHANG Lei, and WANG Guan-yong. Ambiguity resolving and imaging algorithm for multi-channel forward-looking synthetic aperture radar[J]. Journal of Electronics &Information Technology, 2018, 40(12): 2820–2825. doi: 10.11999/JEIT180177
  • [1] 周春晖李飞李宁郑慧芳王翔宇 . 改进的基于特征子空间的SAR图像射频干扰抑制算法. 雷达学报, 2018, 7(2): 235-243. doi: 10.12000/JR17025
    [2] 李伟王兴亮邹鲲徐艺盟张群 . 基于数据融合和陷波滤波的MIMO 雷达抗欺骗干扰算法 (英文). 雷达学报, 2012, 1(3): 246-252. doi: 10.3724/SP.J.1300.2012.20060
    [3] 周超刘泉华胡程 . 间歇采样转发式干扰的时频域辨识与抑制. 雷达学报, 2019, 8(1): 100-106. doi: 10.12000/JR18080
    [4] 陈希信 . 天波雷达后多普勒自适应波束形成. 雷达学报, 2016, 5(4): 373-377. doi: 10.12000/JR15124
    [5] 邹鲲吴德伟张斌李伟 . 一种针对失配信号的可调自适应检测器. 雷达学报, 2015, 4(4): 411-417. doi: 10.12000/JR14129
    [6] 李永祯胡万秋孙豆李重威 . 一种基于极化信息的机载拖曳式诱饵存在性检测与抑制方法研究. 雷达学报, 2016, 5(6): 666-672. doi: 10.12000/JR16115
    [7] 刘宁赵博黄磊 . 单通道SAR抗欺骗干扰方法. 雷达学报, 2019, 8(1): 73-81. doi: 10.12000/JR18072
    [8] 王文钦陈慧郑植张顺生 . 频控阵雷达技术及其应用研究进展. 雷达学报, 2018, 7(2): 153-166. doi: 10.12000/JR18029
    [9] 刘永才王伟潘小义代大海 . 基于延迟-移频的SAR 有源欺骗干扰有效区域研究. 雷达学报, 2013, 2(1): 46-53. doi: 10.3724/SP.J.1300.2012.13001
    [10] 曾铮张福博陈龙永卜祥玺周思言 . 基于多输入多输出合成孔径雷达的二维混合基线抗欺骗干扰方法. 雷达学报, 2019, 8(1): 90-99. doi: 10.12000/JR18118
    [11] 黄翀鹏王剑徐保国 . 旁瓣抑制对线性调频脉冲移频干扰的影响. 雷达学报, 2012, 1(3): 270-276. doi: 10.3724/SP.J.1300.2012.20028
    [12] 段晓菲刘志文徐友根 . 四元数域宽带鲁棒自适应波束形成. 雷达学报, 2019, 8(1): 117-124. doi: 10.12000/JR18083
    [13] 张驰李悦丽周智敏 . 基于独立分量分析法的稀疏阵列穿墙成像雷达直达波干扰抑制. 雷达学报, 2014, 3(5): 524-532. doi: 10.3724/SP.J.1300.2014.14066
    [14] 秦记东赖涛赵拥军黄洁白冰 . 基于通道误差校准的空域导向矢量多通道SAR-GMTI 杂波抑制方法. 雷达学报, 2014, 3(1): 70-77. doi: 10.3724/SP.J.1300.2014.130118
    [15] 刘天鹏魏玺章刘振丁锦灿 . 交叉眼干扰研究综述. 雷达学报, 2019, 8(1): 140-153. doi: 10.12000/JR19013
    [16] 周阳毕大平沈爱国房明星 . 基于运动调制的SAR-GMTI间歇采样遮蔽干扰方法. 雷达学报, 2017, 6(4): 359-367. doi: 10.12000/JR16075
    [17] 丁永超刘成城赵拥军刘亚奇 . 基于反对角单位阵的窄带多径信号接收波束形成算法. 雷达学报, 2015, 4(2): 192-198. doi: 10.12000/JR14012
    [18] 许京伟廖桂生 . 前视阵FDA-STAP雷达距离模糊杂波抑制方法. 雷达学报, 2015, 4(4): 386-392. doi: 10.12000/JR15101
    [19] 许述文石星宇水鹏朗 . 复合高斯杂波下抑制失配信号的自适应检测器. 雷达学报, 0, (): 1-9. doi: 10.12000/JR19030
    [20] 陈伟万显荣张勋饶云华程丰 . 外辐射源雷达多通道时域杂波抑制算法并行实现. 雷达学报, 2014, 3(6): 686-693. doi: 10.12000/JR14157
  • 加载中
图(7)表(2)
计量
  • 文章访问数:  212
  • HTML浏览量:  91
  • PDF下载量:  65
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-09-27
  • 录用日期:  2018-10-23
  • 网络出版日期:  2018-11-22
  • 刊出日期:  2019-02-28

前视多通道SAR自适应鉴别抑制欺骗干扰

    通讯作者: 张磊, leizhang@xidian.edu.cn
    作者简介: 孟智超(1997–),男,陕西宝鸡人,博士生。2018年6月在西安电子科技大学电子工程学院获得学士学位,现为西安电子科技大学电子工程学院博士研究生。主要研究方向为合成孔径雷达干扰与抗干扰;卢景月(1994–),男,河北沧州人,博士生。2016年在中国海洋大学信息科学与工程学院获得学士学位,现为西安电子科技大学电子工程学院博士研究生。主要研究方向为合成孔径雷达成像;张 磊(1984–),男,浙江金华人,博士,副教授。2012年在西安电子科技大学电子工程学院获得博士学位,现担任西安电子科技大学电子工程学院副教授。研究方向为雷达成像、目标识别。E-mail: leizhang@xidian.edu.cn
  • 西安电子科技大学雷达信号处理国家重点实验室   西安   710071
基金项目:  国家自然科学基金(61771372, 61771367),上海市自然科学基金(16ZR1434900)

摘要: 前视多通道SAR成像存在回波左右模糊的问题,在成像时需要利用空域资源进行解模糊处理,使得其成像过程较一般侧视SAR更加复杂。在复杂电磁环境下,若要获得无干扰、不模糊的前视SAR图像,难度较大。该文提出一种基于方位向自适应波束形成(AADBF)技术的欺骗干扰自适应鉴别抑制算法。该算法首先采用AADBF技术对多通道接收回波信号进行地物相消,保留欺骗干扰样本。然后利用门限检测法在高分辨SAR图像上,鉴别出欺骗干扰的像素位置。最后,对存在干扰的像素点进行自适应空域滤波,从而达到抗欺骗干扰的目的。仿真结果表明,该方法能够有效地鉴别并抑制欺骗干扰,与此同时能够实现前视SAR无干扰聚焦成像。

English Abstract

    • 前视合成孔径雷达能够对平台正前方区域进行大范围、高分辨成像[1,2],在军用以及民用中具有良好的应用前景。在军事侦查中,前视SAR系统往往要面对复杂电磁环境。此时,SAR系统将受到一定程度的干扰,当干扰严重时系统将无法有效成像。在众多的干扰方法中,有源欺骗干扰是较难鉴别与抑制的干扰之一[3]。有源欺骗干扰能够在SAR图像上模拟出偏离干扰机位置的大面积虚假目标来掩盖真实目标,使得对真实目标的识别变得更为困难。由于前视SAR多用于弹载制导领域,虚假目标的出现将对弹载信号与信息处理机能否正确识别真实目标产生致命影响,因而鉴别并抑制有源欺骗干扰是前视SAR成像必不可少的处理步骤。

      在战场侦察、精确制导等领域,亟需对欺骗干扰进行鉴别和抑制的有效方法。文献[4]提出利用距离向DBF-SAR系统进行欺骗干扰抑制,但实际工作中干扰回波与真实目标回波必然是面目标且有多个不同方向,文中所述的方法将无法确定干扰与真实目标各自的方向,因而无法进行干扰抑制。文献[5]提出基于4通道相干积累的抗欺骗干扰方法,但是实际中通道数有限且通道误差不可避免,使得相干积累后信干比提升有限。文献[6]提出SAR通过发射随机线性调频斜率信号,使得干扰机侦测并发射的信号落后于真实回波信号,从而在匹配滤波时得不到相应的处理增益来提高信干比。但是可以看到,当调频斜率变化系数增大时成像质量严重下降,干扰抑制能力有限。此外,在前视成像模式下抗干扰与解模糊所产生的空域资源竞争也是一个关键问题,只有合理的分配空域资源才能利用有限通道做到抗干扰、解模糊同时进行。

      距离向DBF-SAR[7]系统的提出,为波束形成技术在SAR领域的应用起到了举足轻重的作用。本文针对前视SAR模式下的欺骗干扰抑制问题,提出一种欺骗干扰自适应鉴别抑制算法。其主要思想是利用欺骗干扰回波与真实地物回波空时关系的差异进行干扰与真实目标的分离,利用门限检测法在SAR图像上进行干扰样本及干扰像素位置的鉴别,然后以空间谱估计法得到干扰方位角度,最后利用自适应空域滤波法进行干扰抑制,实现抗欺骗干扰的目的。该方法有效地解决了文献[4]不能对面目标欺骗干扰进行鉴别与抑制的问题,相对于文献[5]而言大大提高了处理后图像信干比,同时不会出现文献[6]中图像质量下降的问题。本文在信号处理时,利用在阵列通道上以步长为1的滑窗法依次取出多通道信号进行干扰抑制与解模糊成像,使得空域资源充分利用,在通道数较少的情况下实现最大阵列孔径干扰抑制与解模糊。针对前视SAR具有高分辨的特性以及欺骗干扰偏离真实位置的特性,提出局部干扰抑制法,在干扰抑制的同时完整的保留图像信息。

    • FLMC-SAR (Forward-Looking Multi-Channel SAR)几何模型如图1所示。以雷达航迹在地面的投影为$X$轴,以阵列排布方向为$Y\,$轴,平台离地高度向为$Z$轴。假设平台做匀速直线运动,航迹为图中红线所示。天线波束下视角为$\alpha $,目标在地面的方位角表示为$\theta $${R_{xi}}$为接收天线阵列,间距为$d$,均匀分布在${Y_1}{Y_N}$上,发射天线${T_x} $位于阵列中心。即阵列中心的天线为收发共用天线,其余天线为接收天线。

      图  1  FLMC-SAR几何模型

      Figure 1.  FLMC-SAR geometric model

      图1所示,场景中有一干扰机位于${P_ {\rm{J}}}\left( {{x_ {\rm{J}}}, {y_ {\rm{J}}}, 0} \right)$。当平台对场景成像时,干扰机欲在${P_ {\rm{D}}}\left( {{x_ {\rm{D}}}, {y_ {\rm{D}}}, 0} \right)$点形成欺骗干扰点目标。为了简单起见,首先考虑单天线情况,即只有中间一个收发共用天线。此时天线的位置坐标为$\left( {v{t_ {\rm{m}}}, 0, h} \right)$。则干扰机与欺骗点斜距方程如下:

      ${R_ {\rm{J}}}\left( {{t_ {\rm{m}}}} \right) = \sqrt {{{\left( {{x_ {\rm{J}}} - v{t_ {\rm{m}}}} \right)}^2} + {y_ {\rm{J}}}\!\!^{{{2}}} + {h^2}} $

      ${R_ {\rm{D}}}\left( {{t_ {\rm{m}}}} \right) = \sqrt {{{\left( {{x_ {\rm{D}}} - v{t_ {\rm{m}}}} \right)}^2} + y_ {\rm{D}}^2 + {h^2}} $

      其中,$v$为平台速度,${t_ {\rm{m}}}$为慢时间变量,$h$为平台高度。则对应于干扰机位置与欺骗点位置回波如下:

      $ \begin{aligned} {S_ {\rm{J}}}\left( {t, {t_ {\rm{m}}}} \right) = & {A_ {\rm{J}}}{S_0}\left( {t, {t_ {\rm{m}}}} \right) = {A_ {\rm{J}}}{w_ {\rm{r}}}\left( {t - \frac{{2{R_ {\rm{J}}}\left( {{t_ {\rm{m}}}} \right)}}{C}} \right)\\ {\rm{}}& \cdot \exp\left[ {{\rm{j}}{{π}} \gamma {{\left( {t - \frac{{2{R_ {\rm{J}}}\left( {{t_ {\rm{m}}}} \right)}}{C}} \right)}^2}} \right] \\ {\rm{}}&\cdot {w_ {\rm{a}}}\left( {{t_ {\rm{m}}} - \frac{{{x_ {\rm{J}}}}}{v}} \right)\exp\left[ { - {\rm{j}}\frac{{4{{π}} }}{\lambda }{R_ {\rm{J}}}\left( {{t_ {\rm{m}}}} \right)} \right] \end{aligned} $

      $ \begin{align} {S_ {\rm{D}}}\left( {t, {t_ {\rm{m}}}} \right) =& {A_ {\rm{D}}}{S_0}\left( {t, {t_ {\rm{m}}}} \right) = {A_ {\rm{D}}}{w_ {\rm{r}}}\left( {t - \frac{{2{R_ {\rm{D}}}\left( {{t_ {\rm{m}}}} \right)}}{C}} \right)\\ {\rm{}} & \cdot \exp\left[ {{\rm{j}}{{π}} \gamma {{\left( {t - \frac{{2{R_ {\rm{D}}}\left( {{t_ {\rm{m}}}} \right)}}{C}} \right)}^2}} \right] \\ {\rm{}} &\cdot{w_ {\rm{a}}}\left( {{t_ {\rm{m}}} - \frac{{{x_ {\rm{D}}}}}{v}} \right)\exp\left[ { - {\rm{j}}\frac{{4{{π}} }}{\lambda }{R_ {\rm{D}}}\left( {{t_ {\rm{m}}}} \right)} \right] \end{align} $

      其中,${A_ {\rm{J}}}$为干扰机位置处目标回波幅度,${A_ {\rm{D}}}$为欺骗干扰点位置处目标回波幅度。${w_ {\rm{r}}}$${w_ {\rm{a}}}$分别为距离向与方位向窗函数,本文将其简化为矩形窗。

      由式(3)与式(4)可以看出,欺骗干扰回波与干扰机处真实目标回波的差异在于斜距历史的不同,也即斜距历史引起的回波距离向滞后时间与方位向多普勒相位不同。由此可知,为了使干扰机能够发射欺骗干扰,只需对干扰机接收的目标SAR信号进行距离向与方位向调制即可,调制的大小由二者实时距离差决定。对此,干扰生成的方法具体如下:

      对于转发式干扰机形成欺骗干扰的方法为,在回波${S_ {\rm{J}}}\left( {t, {t_ {\rm{m}}}} \right)$的基础上进行调制得到${S_ {\rm{D}}}\left( {t, {t_ {\rm{m}}}} \right)$。对两回波做距离向傅里叶变换得:

      $ \begin{align} {S_ {\rm{J}}}\left( {f, {t_ {\rm{m}}}} \right) =& {A_ {\rm{J}}}{w_ {\rm{r}}}\left( f \right)\exp\left( { - {\rm{j}}{{π}} \frac{{{f^2}}}{\gamma }} \right){w_ {\rm{a}}}\left( {{t_ {\rm{m}}} - \frac{{{x_ {\rm{J}}}}}{v}} \right)\\ {\rm{}}& \cdot \exp\left[ { - {\rm{j}}4{{π}} \frac{{\left( {f + {f\!_{\rm c}}} \right)}}{C}{R_ {\rm{J}}}\left( {{t_ {\rm{m}}}} \right)} \right] \end{align} $

      $ \begin{align} {S_ {\rm{D}}}\left( {f, {t_ {\rm{m}}}} \right) =& {A_ {\rm{D}}}{w_ {\rm{r}}}\left( f \right)\exp\left( { - {\rm{j}}{{π}} \frac{{{f^2}}}{\gamma }} \right){w_ {\rm{a}}}\left( {{t_ {\rm{m}}} - \frac{{{x_ {\rm{D}}}}}{v}} \right)\\ {\rm{}}&\cdot \exp\left[ { - {\rm{j}}4{{π}} \frac{{\left( {f + {f\!_{\rm c}}} \right)}}{C}{R_ {\rm{D}}}\left( {{t_ {\rm{m}}}} \right)} \right] \end{align} $

      其中,$f$为快时间对应频率,${f\!_{\rm c}}$为载频频率。令:

      $\Delta R\left( {{t_ {\rm{m}}}} \right) = {R_ {\rm{D}}}\left( {{t_ {\rm{m}}}} \right) - {R_ {\rm{J}}}\left( {{t_ {\rm{m}}}} \right)$

      ${S_ {\rm{D}}}\left( {f, {t_ {\rm{m}}}} \right)$可由${S_ {\rm{J}}}\left( {f, {t_ {\rm{m}}}} \right)$表示为

      $ \begin{align} {S_ {\rm{D}}}\left( {f, {t_ {\rm{m}}}} \right) =& {S_ {\rm{J}}}\left( {f, {t_ {\rm{m}}}} \right)\\ {\rm{}}& \cdot \exp\left( { - {\rm{j}}4{{π}} \frac{{\left( {f + {f\!_{\rm c}}} \right)}}{C}\Delta R\left( {{t_ {\rm{m}}}} \right)} \right) \end{align} $

      由式(8)可知,干扰机只要对接收到的SAR信号在距离频域调制一个与欺骗点到干扰机之间斜距差有关的线性相位项,再进行逆傅里叶变换得到时域信号发射出去,即可形成欺骗干扰。下面,我们写出前视SAR系统接收到的多通道干扰回波信号模型。

      在满足窄带信号与平台飞行距离很小的假设下,干扰机方位角度为${\theta _ {\rm{J}}}$且在成像期间保持不变。则此时干扰机发射的欺骗信号被每个接收通道接收,各通道欺骗干扰回波信号如下:

      $ \begin{align} {S_{{\rm{r}} i}}\left( {t, {t_ {\rm{m}}}} \right) =& {S_ {\rm{D}}}\left( {t, {t_ {\rm{m}}}} \right)\exp\left[ {\frac{{ - {\rm{j}}2{{π}} \sin\left( {{\theta _ {\rm{J}}}} \right)\!{Y_i}}}{\lambda }} \right],\\ {\rm{}}& \quad i = 1, 2,·\!·\!· , N \end{align} $

      ${Y_i}$为第$i$个接收通道的$Y\;$轴坐标。分析式(5)、式(8)、式(9)可得如下结论:欺骗干扰回波对每一个单通道接收来说都是与真实回波无差异的,但是将多个接收通道回波联合来看,欺骗干扰回波成像位置在欺骗干扰点${P_ {\rm{D}}}$处,但是其回波到达角却与干扰机真实位置${P_ {\rm{J}}}$点一致。这就为我们鉴别欺骗干扰提供了思路,即查看每一像素点的成像位置是否与其波达角一致,不一致则是欺骗干扰。

    • 上节中,通过对比干扰机处真实SAR回波与欺骗干扰回波的差异,发现欺骗干扰是可以由干扰机对真实SAR信号进行调制得到的。由式(8)中与${f\!_{\rm c}}$相关的相位项可知,对于方位向信号来说,由于干扰信号是对真实目标信号调制了一个附加多普勒相位项。使得对于前视SAR接收到的多通道回波而言,干扰信号的空时关系不同于真实目标回波信号的空时关系[7,8]。如图2所示。

      图  2  欺骗干扰与真实回波空时关系差异

      Figure 2.  Difference between deception jamming and real echo space time relationship

      图2中,欺骗干扰空时谱由于附加了一个与实时斜距差有关的多普勒相位,导致欺骗干扰回波空时关系偏离真实空时谱并在多普勒轴上切出一个截距,这一截距正是附加的多普勒值。此时欺骗干扰的空时关系如下:

      $ {f\!_ {\rm{dD}}} = \frac{{2v{\cos}\left( \theta \right)\sin \left( \alpha \right)}}{\lambda } + \Delta {f\!_ {\rm{d}}} = {f\!_ {\rm{dJ}}} + \Delta {f\!_ {\rm{d}}} $

      其中,${f\!_ {\rm{dD}}}$${f\!_ {\rm{dJ}}}$分别为${P_ {\rm{D}}}$${P_ {\rm{J}}}$处目标回波的多普勒频率,$\Delta {f\!_ {\rm{d}}}$为附加多普勒值。

      同时,可以看到图中的欺骗干扰空时谱的斜率为真实回波空时谱的一半。这是由于,当干扰机已知SAR信号参数时,仅对目标SAR发射干扰而非调制转发SAR回波。因此干扰回波为单程,多普勒与速度为1倍关系:

      $ {f\!_ {\rm{dD}}} = \frac{{v{\cos}\left( \theta \right)\sin \left( \alpha \right)}}{\lambda } + \Delta {f\!_ {\rm{d}}} $

      从这一差异图可以明显看出,利用其中任何一种空时谱就可以轻而易举地提取该信号并抑制另一信号,或者抑制该信号提取另一信号。由于欺骗干扰空时关系未知,因此,可以选择利用真实地物回波空时谱进行真实地物回波提取,并抑制欺骗干扰。也可以先提取欺骗干扰,在获得干扰信息后进行干扰抑制。

      为了能够获得干扰更多的信息,本文选择先提取干扰,在获得干扰信息后进行干扰抑制。从图2中可以看到,欺骗干扰空时谱与真实地物回波空时谱无重叠,也就是在进行地物相消时不损失干扰信号,且干扰信号中不含有地物回波信号。

      前视情况下,由于方位采样的原因从回波信号中不能直接得到真实的空时关系。一种免于求解空时谱关系的方法是对回波进行BP (Back Projection)成像,对方位角度进行网格划分,此时每个方位角已知,因而不必计算。首先,构造可以抑制当前角度真实地物回波的正交投影权矢量[9,10]。由于前视回波的左右模糊特性,需要得到能够同时在对称角度上形成凹口的正交投影权。角度${\theta _i}$对应的空域导向矢量为

      $ \begin{align} {{{V}}'}\!\!_{1i} =& \left[ 1\ \exp\left( {\frac{{{\rm{j}}2{{π}} \sin\left( {{\theta _i}} \right)}}{\lambda }} \right)\; ·\!··\right. \\ {\rm{}}& \quad\left.\exp\left( {\frac{{{\rm{j}}2{{π}} \sin\left( {{\theta _i}} \right)\left( {N - 1} \right)}}{\lambda }} \right) \right]^{\rm{T}}, \\ {\rm{}}& \quad i = 1, 2, ·\!·\!·, M \end{align} $

      则,角度${\theta _i}$所对应的正交投影矩阵为

      ${{{P}}_{1i}} = {{I}} - {{{V}}\!\!_{1i}}{{V}}\!_{1i}^{\,\rm{H}}$

      ${{I}}$为一$N \times N$的单位矩阵,H为共轭转置,${{{V}}_{1i}}$${{{V}}'}\!\!_{1i}$的归一化导向矢量。与其对称角度$ - {\theta _i}$的导向矢量计算如下:

      ${{V}}{'\!\!_{2i}} = {{{P}}\!_{1i}}{\rm{conj}}\left( {{{{V}}\!\!_{1i}}} \right)$

      ${\rm{conj}}\left( \cdot \right)$为取共轭。则此角度对应的正交投影矩阵${{{P}}\!_{2i}}$仍利用式(13)所给形式计算。则同时正交于${\theta _i}$$ - {\theta _i}$的投影矩阵为

      ${{{P}}\!_i} = {{{P}}\!_{2i}}{{{P}}\!_{1i}}$

      这一投影矩阵的每一行均与导向矢量${{{V}}\!\!_{1i}}$${{{V}}\!\!_{2i}}$正交,所以仅选取其中任意一行即可得到权矢量。

      $ {{{W}}\!\!_i} = {{N}}{{{P}}\!_i} $

      ${{N}}$为一可选取${{{P}}\!_i}$任意一行的常数矩阵。

      利用上述步骤便可以得到抑制真实地物回波的空域正交权矢量。然后对BP成像后的每一方位角上所有距离单元进行空域滤波便可以抑制真实地物回波。假设第$i$列含有欺骗干扰,其对应的方位角为${\theta _i}$,则这一角度上距离单元为$R$的回波快拍为

      $ \begin{align} {{{S}}_i}\left( {{\theta _i}, R} \right) =& \sigma \left( {{\theta _{1i}}, R} \right) \cdot {{{V}}'}\!\!_{1i} + \sigma \left( {{\theta _{2i}}, R} \right) \cdot {{{V}}'}\!\!_{2i} \\ {\rm{}}&+ {\sigma _ {\rm{D}}} \cdot {{{V}}'}\!\!_ {\rm{D}} \end{align} $

      式中,$\sigma $为与距离$R$,角度$\theta $有关的真实地物散射系数。${\sigma _ {\rm{D}}}$为欺骗干扰散射系数,由干扰发射功率决定,${{{V}}'}\!\!_ {\rm{D}}$为干扰机处的导向矢量。

      对式(17)进行空域滤波,即右乘以权矢量${{{W}}\!\!_i}$

      $ \begin{align} {{{S}}_{i\ {\rm res}}} =& {{{W}}\!\!_i}{{{S}}_i}\left( {\theta , R} \right) = \sigma \left( {{\theta _{1i}}, R} \right) \cdot {{{NP}}\!\!_i} \cdot {{{V}}'}\!\!_{1i} \\ {\rm{}}&+ \sigma \left( {{\theta _{2i}}, R} \right) \cdot {{N}}{{{P}}\!\!_i} \cdot {{{V}}'}\!\!_{2i} + {\sigma _ {\rm{D}}} \cdot {{{NP}}\!_i} \cdot {{{V}}'}\!\!_ {\rm{D}} \\ = &{\sigma _ {\rm{D}}} \cdot {{N}}{{{P}}\!_i} \cdot {{{V}}'}\!\!_ {\rm{D}} \end{align} $

      从式(18)可以看出,利用权矢量抑制真实地物回波的同时由于欺骗干扰导向矢量与权矢量不正交,所以干扰不受影响,被保留下来。

      在利用上述地物相消法后,回波数据中的真实地物回波将被完全抑制,使得干信比大大提升。这时可利用门限检测法获得具有高干信比的干扰样本,利用高干信比的干扰样本将得到精确的干扰估计角度。最后,利用空域自适应波束形成对含干扰数据进行空域滤波便能够有效地抑制所有干扰,得到无干扰图像。

      在上述空域滤波时,一般做法是将得到的空域滤波器对全空域数据进行滤波。但是,此方法在抑制欺骗干扰的同时,不可避免地将与干扰机同一方向上的真实地物回波一同滤除,图像上将出现因信号损失而形成的黑色区域。由前面几节所述,欺骗干扰在图像上所覆盖的位置并非其真实位置。如式(17)所示,图像上含干扰的方位位置所对应的角度${\theta _{1i}}$${\theta _{2i}}$并不是欺骗干扰波达方向${\theta _ {\rm{D}}}$。此时利用估计得到的干扰方位角${\hat \theta _ {\rm{D}}}$获得干扰导向矢量${\hat {{V}}'}\!\!_ {\rm{D}}$,利用式(13)计算得到正交投影矩阵,同时利用式(16)获得一个正交投影权矢量${\hat {{W}}\!_ {\rm{D}}}$。此时的权矢量仅抑制欺骗干扰而对真实地物回波不加抑制,其表达式如下:

      $ \begin{align} {{{S}}_{i \ {\rm res}}} = & {{\hat {{W}}}\!_ {\rm{D}}}{{{S}}_i}\left( {\theta , R} \right) = \sigma \left( {{\theta _{1i}}, R} \right) \cdot {{N}}{{\hat {{P}}}\!_ {\rm{D}}} \cdot {{{V}}'}\!\!_{1i} \\ {\rm{}}&+ \sigma \left( {{\theta _{2i}}, R} \right) \cdot {{N}}{{{P}}\!_ {\rm{D}}} \cdot {{{V}}'}\!\!_{2i} + {\sigma _ {\rm{D}}} \cdot {{N}}{{\hat {{P}}}\!_ {\rm{D}}} \cdot {{\hat {{V}}}'}\!\!_ {\rm{D}} \\ = &\sigma \left( {{\theta _{1i}}, R} \right) \cdot {{N}}{{\hat {{P}}}\!_ {\rm{D}}} \cdot {{{V}}'}\!\!_{1i} \\ {\rm{}}& + \sigma \left( {{\theta _{2i}}, R} \right) \cdot {{N}}{{{P}}\!_ {\rm{D}}} \cdot {{{V}}'}\!\!_{2i} \end{align} $

      可以看到,在对干扰所在像素区域进行空域滤波时,由于真实地物回波的到达角不同于欺骗干扰,所以,真实地物回波信号不受损失而干扰被抑制。如果干扰抑制仅在干扰所在的像素位置上进行,则与欺骗干扰同一方向的地物回波不受损失,且被欺骗干扰掩盖的地物回波也不受损失,所以图像将完整保留。本文将此方法称之为自适应欺骗干扰局部抑制法。这一方法行之有效的基础是前视SAR系统可以得到高分辨的图像,使得在地物抑制后可以利用门限检测法来对欺骗干扰进行像素位置的定位。

      前视多通道合成孔径雷达成像的基本步骤为:回波距离脉压、BP成像[11]、解模糊成像[12,13]。由于干扰不能被解模糊,且解模糊之后将影响干扰样本的提取,所以干扰抑制需在解模糊之前进行。则前视多通道合成孔径雷达抗欺骗干扰成像的算法具体流程见图3

      图  3  欺骗干扰抑制算法流程

      Figure 3.  Deception interference suppression algorithm flow

    • 在真实环境中当弹载雷达锁定真实舰船时,干扰机发射欺骗信号形成欺骗目标。若弹载计算机无法判断欺骗干扰与真实目标,则将误导弹载信号与信息处理机做出错误判断。因此本文仿真设置见表1

      参数数值
      载频35 GHz
      信号带宽50 MHz
      脉冲宽度0.5 μs
      采样带宽60 MHz
      平台高度4000 m
      脉冲重复频率9 kHz
      波束方位宽度20°
      波束俯仰宽度22°
      场景中心斜距${R_s}$8 km
      阵列实孔径长度0.0642 m

      表 1  系统参数

      Table 1.  System parameters

      采用前视9通道线性阵列接收回波,信号处理时通道资源配置如下:干扰抑制用5通道进行,波束形成解模糊也用5通道进行,多通道数据获取方法为如前所述的步长为1的滑窗法。目标设置如下:利用真实SAR图像作为背景仿真得到真实目标回波,仿真得到真实方位角为5.17°虚假方位角为(–2°, –3°)之间的3艘舰船欺骗干扰目标。

      实验模拟直线航迹下的FLMC-SAR系统进行回波仿真,采用上述算法对回波进行欺骗干扰鉴别抑制成像如下:

      加入欺骗干扰后进行前视单通道BP模糊成像的结果如图4(b)所示,可以看到在泽碧岛下部小船的附近出现了3艘反射很强的舰船目标,即为欺骗舰船。

      图  4  加入欺骗干扰前后对比

      Figure 4.  Comparison before and after adding spoofing jamming

      图5(a)所示,地物抑制后图像域仅保留欺骗干扰。此时,干信比大大提升使得干扰样本的提取更加容易。由图5(b)可以看出,在利用门限检测法后几乎所有欺骗干扰都被正确定位,图5(a)中船只之间干扰旁瓣所在区域也被正确选出。

      图  5  地物相消与干扰位置定位

      Figure 5.  Location cancellation and interference position locating

      图6即为利用获得的高干信比的干扰样本估计的欺骗干扰方向,与设置值5.17°对比来看误差为0.006°,这一误差值在实际数据处理中远小于误差容忍值。

      图  6  干扰角度估计结果

      Figure 6.  Interference angle estimation results

      图7(a)图7(b)为干扰抑制后解模糊成像结果,图7(a)为全空域干扰抑制结果,图7(b)为局部干扰抑制结果。两图对比来看,局部干扰抑制解模糊完整地保存了原始图像,被欺骗干扰覆盖的区域以及与欺骗干扰同方向的区域,地物真实回波信号都未损失图像被完整保留。而对于全空域欺骗干扰抑制,被欺骗干扰覆盖的区域地物真实回波未损失。干扰机位置附近的地物回波损失严重,使得右下角的岛屿及道路完全被抑制,图像信息未能完整保存。

      图  7  欺骗干扰抑制解模糊对比

      Figure 7.  Comparison of deception interference suppression ambiguity resolution

      处理前后的性能指标测得如表2所示。其中信干比增益为干扰机方向上干扰抑制后信干比与抑制前信干比之比的分贝值,干扰抑制比为处理前后欺骗干扰功率比,信号抑制比为处理前后干扰机方向上期望信号功率比。

      处理结果信干比增益(dB)干扰抑制比(dB)信号抑制比(dB)
      含干扰图像–15.6000
      全域干扰抑制030.7530.75
      局部干扰抑制32.0631.75–0.31

      表 2  处理性能对比

      Table 2.  The comparison of processing performance

      从结果图可以看出,传统的全空域干扰抑制法在去除干扰的同时,使得场景右边干扰机附近的图像被抑制,形成一条黑带区域,图像信息损失。这一点从表2也可以看出,其中信号与干扰都被抑制了30.75 dB。而本文所提的局部干扰抑制法几乎对图像无损失,图像信息被完整保留。且从表2可以看出,此时的干扰抑制性能良好。

    • 本文针对前视SAR系统面临欺骗干扰难以有效鉴别并抑制的问题,提出了一种欺骗干扰自适应鉴别抑制算法。该算法根据欺骗干扰与真实地物回波空时关系的差异进行干扰的分离与鉴别,利用AADBF法进行局部空域干扰抑制,从而实现了前视SAR系统抗欺骗干扰的目的。该方法不仅能够准确地鉴别抑制面目标欺骗干扰,而且相对于传统的全空域干扰抑制法,不会造成与干扰同一方向上的真实地物回波信号的损失,使得图像信息完整保存。本文方法是基于欺骗干扰回波与真实地物回波空时关系差异较大的情况下提出的,需要进一步研究如何在这一差异较小时进行欺骗干扰的鉴别与抑制。

参考文献 (13)

目录

    /

    返回文章
    返回