MIMO双基地机载雷达距离模糊杂波抑制方法

王宇卓 朱圣棋 许京伟

引用本文:
Citation:

MIMO双基地机载雷达距离模糊杂波抑制方法

    作者简介: 王宇卓(1992–),男,陕西西安人,西安电子科技大学在读博士。主要研究方向为机载雷达空时信号处理。E-mail: 1145416259@qq.com;朱圣棋(1984–),男,江西赣州人,西安电子科技大学教授,博士生导师,西安电子科技大学科研院副院长,科技部综合电子信息系统国家级国际合作基地副主任,雷达信号处理国防科技重点实验室教授,中国电子学会高级会员,中国电子学会雷达信号与系统专委会副秘书长,国家自然科学基金创新群体骨干成员,陕西省质量监督局国家质量技术基础科技委专家,发表学术论文90余篇,其中40余篇发表在SCI期刊,申请专利40余项,其中授权20项。主要研究方向为阵列信号处理、空时2维信号处理、机载/星载雷达地面运动目标检测、参数估计以及成像。E-mail: zhushengqi8@163.com;许京伟(1987–),男,山东日照人,西安电子科技大学讲师,博士,主要研究方向为空时自适应处理、频率分集阵列信号处理等。E-mail: xujingwei1987@163.com.
    通讯作者: 朱圣棋, zhushengqi8@163.com
  • 基金项目:

    国家自然科学基金委创新研究群体项目(61621005),国家重点研发计划战略性国际科技创新合作重点专项(2016YFE0200400),陕西省国际科技合作重点研发计划(2017KW-ZD-12)

A Range-ambiguous Clutter Suppression Method for MIMO Bistatic Airborne Radar

    Corresponding author: Zhu Shengqi, zhushengqi8@163.com ;
  • Fund Project: The National Natural Science Foundation of China Innovation Research Group Project (61621005), The National Key RD Program of China (2016YFE0200400), The Key RD Program of Shaanxi Province (2017KW-ZD-12)

  • 摘要: 正侧视阵MIMO(Multiple-Input Multiple-Output)双基地雷达系统具有严重的距离依赖问题和杂波谱展宽现象,该文分析其杂波特性,提出了一种距离模糊杂波抑制方法。该方法在传统的MIMO双基雷达中引入距离相关的相位项,从而提供额外的自由度,通过该相位项的距离信息实现不同距离的杂波分离,并进一步对其进行有效抑制。从而解决了双基地预警雷达大场景下地面杂波距离杂波模糊的问题,仿真结果显示该方法的正确性,并且和已有的DW(Doppler Warping)距离模糊杂波抑制方法进行比较,使用IF(Improve Factor)谱线对抑制结果进行对比,显示该方法有更优的性能。
  • 图 1  FDA-MIMO双基地雷达几何结构示意图

    Figure 1.  Geometric schematic diagram of FDA-MIMO bistatic radar

    图 2  接收信号处理框图

    Figure 2.  Block diagram of receiving signal processing

    图 3  双基地雷达杂波谱

    Figure 3.  Bistatic radar clutter spectrum

    图 4  雷达数据处理流程

    Figure 4.  Flow chart of radar data processing

    图 6  $\Delta f$ 不同值时对应的杂波谱

    Figure 6.  The corresponding clutter spectrum of $\Delta f$ at different values

    图 5  双基地雷达杂波谱图

    Figure 5.  Bistatic radar clutter spectrum

    图 7  双基地雷达杂波谱

    Figure 7.  Bistatic radar clutter spectrum

    图 8  雷达杂波IF特性

    Figure 8.  IF characteristics of radar clutter

    图 9  双基地MIMO雷达杂波脊

    Figure 9.  Bistatic MIMO radar clutter ridge

    图 10  FDA-MIMO双基地雷达DW补偿杂波谱

    Figure 10.  FDA-MIMO bistatic radar clutter spectrum of DW compensation method

    图 11  FDA-MIMO双基地雷达IF谱线对比图

    Figure 11.  Comparison of FDA-MIMO bistatic radar IF lines

    表 1  双基地雷达仿真参数

    Table 1.  Simulation parameters of bistatic radar

    参数 参数值 参数 参数值
    发射阵元数 6 发射机高度 12 km
    接收阵元数 6 接收机高度 12 km
    参考工作频率 1.104 GHz 双基基线长度 50 km
    阵列阵元间距 0.0543 m PRF 2434.8 Hz
    发射机速度 140 m/s 相干脉冲数 8
    目标距离 64.964 km 频率增量 1217.4 Hz
    注:在传统MIMO双基地雷达中没有 $\Delta $ f 项或者 $\Delta $ f 为0
    下载: 导出CSV
  • [1] Willis N J and Griffiths H D. Advances in Bistatic Radar[M]. Raleigh: Scitech Pubishing, 2007: 15–31.
    [2] Behner F, Reuter S, Nies H, et al. Synchronization and processing in the HITCHHIKER bistatic SAR experiment[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2016, 9(3): 1028–1035. DOI: 10.1109/JSTARS.2015.2471082
    [3] Hajj G A and Zuffada C. Theoretical description of a bistatic system for ocean altimetry using the GPS signal[J]. Radio Science, 2003, 38(5): 1089. DOI: 10.1029/2002RS002787
    [4] Xie W C and Wang Y L. Range-dependence compensation method for bistatic STAP radar[C]. Proceedings of 2008 Congress on Image and Signal Processing, Sanya, Hainan, 2008: 503–506.
    [5] Samuele Gelli, Alessio Bacci, Marco Martorella, et al.. Clutter suppression and high-resolution imaging of noncooperative ground targets for bistatic airborne radar[J]. IEEE Transactions on Aerospace and Electronic Systems, 2018, 54(2): 932–949.
    [6] Melvin W L, Callahan M J, and Wicks M C. Adaptive clutter cancellation in bistatic radar[C]. Proceedings of the Conference Record of the 34th Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, 2000: 1125–1130.
    [7] Borsari G K. Mitigating effects on STAP processing caused by an inclined array[C]. Proceedings of 1998 IEEE Radar Conference, Dallas, TX, 1998: 135–140.
    [8] Himed B, Zhang Y, and Hajjari A. STAP with angle-Doppler compensation for bistatic airborne radars[C]. Proceedings of 2002 IEEE Radar Conference, Long Beach, CA, 2002: 311–317.
    [9] Antonik P, Wicks M C, Griffiths H D, et al.. Range-dependent beamforming using element level waveform diversity[C]. Proceedings of 2006 International Waveform Diversity & Design Conference, Lihue, HI, 2006: 1–6.
    [10] Antonik P, Wicks M C, Griffiths H D, et al.. Multi-mission multi-mode waveform diversity[C]. Proceedings of 2006 IEEE Conference on Radar, Verona, NY, 2006: 580–582. DOI: 10.1109/RADAR.2006.1631858.
    [11] Antonik P, Wicks M C, Griffiths H D, et al.. Frequency diverse array radars[C]. Proceedings of 2006 IEEE Conference on Radar, Verona, NY, 2006: 215–217. DOI: 10.1109/RADAR.2006.1631800.
    [12] Wu J X, Wang T, Zhang L F, et al. Range-dependent clutter suppression for airborne sidelooking radar using MIMO technique[J]. IEEE Transactions on Aerospace and Electronic Systems, 2012, 48(4): 3647–3654. DOI: 10.1109/TAES.2012.6324751
    [13] Xie W C, Wang Y L, Zhang B H, et al.. Clutter suppression for bistatic airborne radar with range ambiguity[C]. Proceedings of 2011 IEEE CIE International Conference on Radar, Chengdu, 2011, 2: 1893–1897. DOI: 10.1109/CIE-Radar.2011.6159944.
    [14] Sammartino P F, Baker C J, and Griffiths H D. Frequency diverse MIMO techniques for radar[J]. IEEE Transactions on Aerospace and Electronic Systems, 2013, 49(1): 201–222. DOI: 10.1109/TAES.2013.6404099
    [15] Li J, Liao G S, and Griffiths H. Bistatic MIMO radar space-time adaptive processing[C]. Proceedings of 2011 IEEE Radar Conference, Kansas City, MO, 2011: 498–502. DOI: 10.1109/RADAR.2011.5960587.
    [16] Xu J W, Liao G S, Huang L, et al. Robust adaptive beamforming for fast-moving target detection with FDA-STAP radar[J]. IEEE Transactions on Signal Processing, 2017, 65(4): 973–984. DOI: 10.1109/TSP.2016.2628340
    [17] Klemm R. Comparison between monostatic and bistatic antenna configurations for STAP[J]. IEEE Transactions on Aerospace and Electronic Systems, 2000, 36(2): 596–608. DOI: 10.1109/7.845248
    [18] Xu J W, Zhu S Q, and Liao G S. Range ambiguous clutter suppression for airborne FDA-STAP radar[J]. IEEE Journal of Selected Topics in Signal Processing, 2015, 9(8): 1620–1631. DOI: 10.1109/JSTSP.2015.2465353
  • [1] 韩金旺张子敬刘军赵永波 . 基于贝叶斯的高斯杂波背景下MIMO雷达自适应检测算法. 雷达学报, 2019, (): 1-9. doi: 10.12000/JR18090
    [2] 李伟王兴亮邹鲲徐艺盟张群 . 基于数据融合和陷波滤波的MIMO 雷达抗欺骗干扰算法 (英文). 雷达学报, 2012, 1(3): 246-252. doi: 10.3724/SP.J.1300.2012.20060
    [3] 谷文堃王党卫马晓岩 . 分布式MIMO-ISAR子图像融合方法. 雷达学报, 2017, 6(1): 90-97. doi: 10.12000/JR16042
    [4] 曾铮张福博陈龙永卜祥玺周思言 . 基于多输入多输出合成孔径雷达的二维混合基线抗欺骗干扰方法. 雷达学报, 2019, 8(1): 90-99. doi: 10.12000/JR18118
    [5] 魏民李小波王理 . 减小距离模糊影响的机载双基地雷达配置方法. 雷达学报, 2017, 6(1): 106-113. doi: 10.12000/JR16082
    [6] 许京伟廖桂生 . 前视阵FDA-STAP雷达距离模糊杂波抑制方法. 雷达学报, 2015, 4(4): 386-392. doi: 10.12000/JR15101
    [7] 李军党博刘长赞廖桂生 . 利用发射角度的双基地MIMO 雷达杂波抑制方法. 雷达学报, 2014, 3(2): 208-216. doi: 10.3724/SP.J.1300.2014.13148
    [8] 许京伟朱圣棋廖桂生张玉洪 . 频率分集阵雷达技术探讨. 雷达学报, 2018, 7(2): 167-182. doi: 10.12000/JR18023
    [9] 胡程刘长江曾涛 . 双基地前向散射雷达探测与成像. 雷达学报, 2016, 5(3): 229-243. doi: 10.12000/JR16058
    [10] 曾涛 . 双基地合成孔径雷达发展现状与趋势分析. 雷达学报, 2012, 1(4): 329-341. doi: 10.3724/SP.J.1300.2012.20093
    [11] 林春风黄春琳粟毅 . 双基地雷达Radon-Fourier变换弱目标积累检测. 雷达学报, 2016, 5(5): 526-530. doi: 10.12000/JR16049
    [12] 钱李昌许稼胡国旭 . 非合作无源双基地雷达弱目标长时间积累技术. 雷达学报, 2017, 6(3): 259-266. doi: 10.12000/JR16137
    [13] 李学仕孙光才邵鹏吴玉峰邢孟道 . 基于数字阵列雷达的同时多模式SAR 成像体制研究. 雷达学报, 2014, 3(4): 480-489. doi: 10.3724/SP.J.1300.2014.13113
    [14] 祝晓静李飞王宇王伟孙翔 . 基于改进方位相位编码的全极化SAR距离模糊抑制方法. 雷达学报, 2017, 6(4): 420-431. doi: 10.12000/JR17015
    [15] 饶云华朱逢园张修志万显荣龚子平 . WiFi 外辐射源雷达信号模糊函数及副峰抑制分析. 雷达学报, 2012, 1(3): 225-231. doi: 10.3724/SP.J.1300.2012.20061
    [16] 王本静易建新万显荣但阳鹏 . LTE外辐射源雷达帧间模糊带分析与抑制. 雷达学报, 2018, 7(4): 514-522. doi: 10.12000/JR18025
    [17] 王龙刚李廉林 . 基于互质阵列雷达技术的近距离目标探测方法(英文). 雷达学报, 2016, 5(3): 244-253. doi: 10.12000/JR16022
    [18] 艾小锋曾勇虎高磊王肖洋汪连栋 . 飞机目标全极化双基地散射特性研究. 雷达学报, 2016, 5(6): 639-646. doi: 10.12000/JR16070
    [19] 陈伟万显荣张勋饶云华程丰 . 外辐射源雷达多通道时域杂波抑制算法并行实现. 雷达学报, 2014, 3(6): 686-693. doi: 10.12000/JR14157
    [20] 韩雁飞吴仁彪李海 . 基于双门限控制的机载气象雷达地杂波抑制方法. 雷达学报, 2013, 2(1): 97-104. doi: 10.3724/SP.J.1300.2012.20057
  • 加载中
图(11)表(1)
计量
  • 文章访问数:  348
  • HTML浏览量:  139
  • PDF下载量:  297
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-02-10
  • 录用日期:  2018-04-20
  • 刊出日期:  2018-04-28

MIMO双基地机载雷达距离模糊杂波抑制方法

    通讯作者: 朱圣棋, zhushengqi8@163.com
    作者简介: 王宇卓(1992–),男,陕西西安人,西安电子科技大学在读博士。主要研究方向为机载雷达空时信号处理。E-mail: 1145416259@qq.com;朱圣棋(1984–),男,江西赣州人,西安电子科技大学教授,博士生导师,西安电子科技大学科研院副院长,科技部综合电子信息系统国家级国际合作基地副主任,雷达信号处理国防科技重点实验室教授,中国电子学会高级会员,中国电子学会雷达信号与系统专委会副秘书长,国家自然科学基金创新群体骨干成员,陕西省质量监督局国家质量技术基础科技委专家,发表学术论文90余篇,其中40余篇发表在SCI期刊,申请专利40余项,其中授权20项。主要研究方向为阵列信号处理、空时2维信号处理、机载/星载雷达地面运动目标检测、参数估计以及成像。E-mail: zhushengqi8@163.com;许京伟(1987–),男,山东日照人,西安电子科技大学讲师,博士,主要研究方向为空时自适应处理、频率分集阵列信号处理等。E-mail: xujingwei1987@163.com
  • 西安电子科技大学雷达信号处理国家重点实验室   西安   710071
基金项目:  国家自然科学基金委创新研究群体项目(61621005),国家重点研发计划战略性国际科技创新合作重点专项(2016YFE0200400),陕西省国际科技合作重点研发计划(2017KW-ZD-12)

摘要: 正侧视阵MIMO(Multiple-Input Multiple-Output)双基地雷达系统具有严重的距离依赖问题和杂波谱展宽现象,该文分析其杂波特性,提出了一种距离模糊杂波抑制方法。该方法在传统的MIMO双基雷达中引入距离相关的相位项,从而提供额外的自由度,通过该相位项的距离信息实现不同距离的杂波分离,并进一步对其进行有效抑制。从而解决了双基地预警雷达大场景下地面杂波距离杂波模糊的问题,仿真结果显示该方法的正确性,并且和已有的DW(Doppler Warping)距离模糊杂波抑制方法进行比较,使用IF(Improve Factor)谱线对抑制结果进行对比,显示该方法有更优的性能。

English Abstract

    • 双基地雷达具有良好的特性以应对传统单基地雷达所面临的“四大威胁”:电子干扰、超低空突防、反辐射导弹、隐身武器,因此在现代军事领域双基地雷达被美、俄、英等越来越多的国家应用在机载、车载、陆基等雷达系统中,用于完善现有雷达系统的侦测、监视和跟踪功能[1],在成像[2],海洋高度测量[3]方向也有大量的工作成果,而MIMO双基地雷达作为一种新型雷达,优化了传统双基地雷达的性能。相比传统双基地雷达,MIMO双基地雷达提高了雷达的目标参数估计精度并且由于MIMO体系的正交波形在空间形成较宽的波束,简化了双基地雷达空间同步问题,同时因为信号的相互独立性,MIMO双基地雷达信号有更低的截获概率。但是由于机载双基地雷达杂波具有距离依赖性,MIMO双基地雷达依然无法改善距离模糊杂波引起的多普勒频率展宽的问题,其与传统相控阵双基地雷达都会因距离模糊杂波而性能下降,因此距离模糊非均匀杂波[4,5]仍旧是MIMO双基地雷达面临的一个重要问题。

      对于由于杂波依赖性而导致的性能下降,已存在的方法有多普勒频移补偿法(Doppler Warping, DW)[6,7]和角度多普勒补偿法(Angle Doppler Compensation, ADC)[8]。DW将非均匀杂波在主瓣方向上进行杂波的多普勒维平移,ADC则是在角度维和多普勒维2维进行平移,但因为平移后多重距离模糊杂波依旧重叠在一起,所以该方法改善效果有限,且当存在多重距离模糊时由于算法稳健性难以保证,所以导致杂波对齐后杂波抑制效果不明显。2006年,Antonik等人[911]提出了频率分集阵列(Frequency Diverse Array, FDA)的概念,FDA雷达通过添加阵元载频间的微小增量,空间频率相位增加了距离耦合的相位分量,使得角度距离维的波束指向随距离的变化而发生周期性的改变,从而在距离域增加了额外的自由度,在杂波抑制方面,FDA的特性使得多重距离模糊信号可以在距离的维度被加以区分。

      本文将讨论一种新的MIMO双基地雷达体系FDA-MIMO双基地雷达,通过相邻阵元发射频率设置微小的频率增量,并且在发射时加入相互正交的波形,将距离耦合相位分量引入双基地雷达中,该方法使用距离耦合的相位实现不同距离的距离模糊杂波在角度维实现分离,从而起到距离模糊杂波抑制效果,用以解决MIMO双基地雷达中距离模糊杂波的多普勒展宽问题[12,13],仿真结果显示该方法能够改善双基地雷达因距离模糊杂波引起的性能下降现象,且性能优于DW技术。

    • 在传统MIMO双基地雷达系统中,由于脉冲周期长度的限制,位于相应位置不同距离环的模糊杂波在回波信号中会重叠在一起,导致杂波展宽现象严重,雷达性能下降,这是原系统的不足之处。针对上述问题,本文将在传统MIMO的基础上进行改进。本节将给出FDA-MIMO双基地雷达的基本模型[14,15],分析其信号模型的变化。

      图1为双基地雷达的几何结构示意图,在(x, y, z)坐标系中,发射阵列在(BL, 0, ht)处,接收阵列在(0, 0, hr)处,hthr分别为发射接收机载平台的高度,BL为基线长度,且阵列均为正侧视阵,目标ci与发射端距离设为Rt,目标距接收端距离设为Rr,目标和发射阵列夹角设为 ${\psi _{\rm{t}}}$ ,和接收阵列夹角设为 ${\psi _{\rm{r}}}$ ,发射阵列阵元数为M,接收阵列阵元数为N。在原双基地MIMO雷达的基础上,FDA-MIMO双基地雷达将每一个发射阵列阵元的发射频率f0依次累加相同的频率增量 $\Delta $ f,即:

      图  1  FDA-MIMO双基地雷达几何结构示意图

      Figure 1.  Geometric schematic diagram of FDA-MIMO bistatic radar

      ${f_m}{\rm{ = }}{f_{\rm{0}}}{\rm{ + }}m\Delta f,\ {\rm{ }}m = 0,1, ·\!·\!· ,M - 1$

      其中,m表示发射阵元的编号,f0是参考频率, $\Delta $ f是频率增量,且要求 $\Delta $ f<<f0。与传统MIMO双基地雷达信号形式相似,由式(1)得到FDA-MIMO双基地雷达第m通道发射信号为 $ {\hat x_m}(t) = {s_m}(t) $ . $\exp ({\rm j}2{{π}} {f_m}t){\rm{ }},m = 0,1,2, ·\!·\!· M - 1$ ,可以看出与传统MIMO双基地雷达不同的是不同通道的频率由f 变为fm。其中,sm(t)为正交波形,其满足正交条件:

      $\int\nolimits_{{T_{\rm{p}}}} {{s_i}(t)s_j^{\rm{H}}(t){\rm d}t = } \left\{ {\begin{array}{*{20}{c}} 1,&{i = j} \\ 0,&{i \ne j} \end{array}} \right.$

      式中,Tp为发射脉冲宽度,H表示共轭转置运算符。

      到达目标时窄带信号为 $ {\tilde x_{}}(t,{\psi _{\rm{t}}},{R_{\rm{t}}}) =$ $ \displaystyle\sum\limits_{m \!=\! 0}^{M \!-\! 1} \!{{s\!_m}(t)\exp\! \left({\rm \!j}2{{π}} {f\!_m}\!\left(\!\!{t \!-\! \frac{{{R_{\rm t}} \!-\! m{d_{\rm t}}\cos {\psi _{\rm t}}{\rm{ }}}}{c}}\!\!\right) \!\!+\!\! {f\!_{m{\rm{d}}}}t\!\!\right) \!+\! {{\tilde n}_{\rm c}}} $ ,其中m=0, 1, ···, M–1, Rt为目标与发射端的距离,dt为发射阵列阵元间距, ${\psi _{\rm t}}$ 为目标和发射阵列夹角, ${\tilde n_{\rm{c}}}$ 表示传播过程中的高斯白噪声,fmd为对应于发射平台频率fm引起的多普勒频率偏移。信号达到目标后经过目标反射,被接收端阵列进行接收,通过MIMO正交波形滤波处理,实现不同通道信号分离。

      图2xn表示第n个接收通道接收的回波信号, $s_m^{\rm H} $ 表示第m个发射通道所对应的正交波形的共轭转置,则经过式(2)的正交波形积分和匹配滤波后,在FDA-MIMO双基地雷达第k个脉冲第m个通道发射第n个通道接收的回波信号为:

      图  2  接收信号处理框图

      Figure 2.  Block diagram of receiving signal processing

      $\begin{align} {x_{mnk}}({R_{\rm{t}}} + {R_{\rm{r}}},{\psi _{\rm{t}}},{\psi _{\rm{r}}},v,k) =& {\rho _{mn}}\exp \left\{ {{\rm j}2{{π}} {f_m}\left( { - \frac{{{R_{\rm{t}}} - m{d_{\rm{t}}}\cos {\psi _{\rm t}}}}{c} - \frac{{{R_{\rm{r}}} - n{d_{\rm{r}}}\cos {\psi _{\rm{r}}}}}{c}} \right) + {\rm j}2{{π}} {f_m}_{\rm{d}}/{f_{{\rm{prf}}}}k} \right\} + {{\tilde n}_{\rm{c}}} \\ =& {\rho _{mn}}\exp \left\{ {{\rm j}2{{π}} {f_m}\left( { - \frac{{{R_{\rm{t}}} + {R_{\rm{r}}}}}{c} + \frac{{m{d_{\rm{t}}}\cos {\psi _{\rm{t}}}}}{c} + \frac{{n{d_{\rm{r}}}\cos {\psi _{\rm{r}}}}}{c}} \right) + {\rm j}2{{π}} {f_{m{\rm{d}}}}k} \right\} + {{\tilde n}_{\rm{c}}} \\ =& {\rho _{mn}}\exp \left\{ { - {\rm j}2{{π}} {f_m}\frac{{{R_\sum }}}{c} + {\rm j}2{{π}} {f_m}\left( {\frac{{m{d_{\rm{t}}}\cos {\psi _{\rm{t}}}}}{c} + \frac{{n{d_{\rm{r}}}\cos {\psi _{\rm{r}}}}}{c}} \right) + {\rm j}2{{π}} {f_{m{\rm{d}}}}k} \right\} + {{\tilde n}_{\rm{c}}} \end{align} $

      其中,ρmn为照射目标对应的散射系数,Rr表示目标与接收端的距离,dr为接收阵列阵元间距, ${\psi _{\rm{r}}}$ 是目标和接收阵列夹角, ${R_\sum } = {R_{\rm{t}}} + {R_{\rm{r}}}$ 表示发射阵元和接收阵元到目标的距离和,fmd表示频率fmd对应的归一化多普勒频率,即f m d= $ \displaystyle\frac{{{v_{\rm{t}}}}}{{c{f_{\rm{r}}}}}{f\!_m}\cos {\psi _{\rm{t}}} + \displaystyle\frac{{{v_{\rm{r}}}}}{{c{f_{\rm{r}}}}}{f\!_m}\cos {\psi _{\rm{r}}} = \displaystyle\frac{{{v_{\rm{t}}}}}{{{\lambda _m}c{f_{\rm{r}}}}}\cos {\psi _{\rm{t}}} + \displaystyle\frac{{{v_{\rm{r}}}}}{{{\lambda _m}{f_{\rm{r}}}}}$ $ \cos {\psi _{\rm{r}}}$ ,由于 $\Delta $ f<<f0,所以这里 ${\lambda _m} \approx {\lambda _0}$ , ${f\!_{m{\rm{d}}}} \approx {f_{\rm{d}}} = \displaystyle\frac{{{v_{\rm{t}}}}}{{{\lambda _0}c{f_{\rm{r}}}}}\cos {\psi _{\rm{t}}} + \displaystyle\frac{{{v_{\rm{r}}}}}{{{\lambda _0}{f_{\rm{r}}}}}\cos {\psi _{\rm{r}}}$ [14,16]。式(3)中设 ${\phi _n} = - 2{{π}} {f\!_m}\displaystyle\frac{{{R_\sum }}}{c}$ ,则有:

      $\begin{align} \Delta {\phi _m} = &{\phi _m} - {\phi _0} = - 2{{π}} m\Delta f\frac{{{R_\sum }}}{c} \\ \end{align} $

      所以式(4)可进一步表示为:

      $\begin{align} {x_{mnk}}({R_{\rm{t}}} + {R_{\rm{r}}},{\psi _{\rm{t}}},{\psi _{\rm{r}}},v,k) =& {\rho _{mn}}\exp \left\{ { - {\rm j}2{{π}} {f\!_m}\frac{{{R_\sum }}}{c} + {\rm j}2{{π}} {f\!_m}\left( {\frac{{m{d_{\rm{t}}}\cos {\psi _{\rm{t}}}}}{c} + \frac{{n{d_{\rm{r}}}\cos {\psi _{\rm{r}}}}}{c}} \right) + {\rm j}2{{π}} {f\!_{\rm{d}}}k} \right\} + {{\tilde n}_{\rm{c}}}\\ =& {{Γ} _0} \cdot {\rho _{mn}}\exp \left\{ { - {\rm j}2{{π}} m\Delta f\frac{{{R_\sum }}}{c} \!+\!\! {\rm j}2{{π}} {f\!_m}\left( {\frac{{m{d_{\rm{t}}}\cos {\psi _{\rm{t}}}}}{c} \!+\! \frac{{n{d_{\rm{r}}}\cos {\psi _{\rm{r}}}}}{c}} \right) \!+\! {\rm j}2{{π}} {f\!_{\rm{d}}}k} \right\}\!+\! {{\tilde n}_{\rm{c}}} \end{align}$

      其中 ${{Γ}_0}$ 为各项包含f0的共有项:

      $\begin{align} {{Γ}_0} =& \exp \left[ {{\rm j}{\phi _{\rm{0}}}} \right] = \exp \left[ { - {\rm j}2{{π}} {f_{\rm{0}}}\frac{{{R_\sum }}}{c}} \right] \\ \end{align} $

      相干处理脉冲数为K,接收阵元发射阵元分别为MN,接收信号的空时快拍表示为:

      $\begin{align} {s} =& {{{Γ}}_0}{ρ} {{s}_{\rm{d}}}({{f}_{\rm{d}}}) \otimes \left\{ {{{s}_{\rm{r}}}({{f}_{\rm{r}}}) \odot {{s}_{{\rm{at}}}}({{f}_{{\rm{at}}}})} \right\} + {{\tilde n}_{\rm{c}}}\\ = &{{{Γ}}_0}{ρ} {{s}_{\rm{d}}}({{f}_{\rm{d}}}) \!\otimes \!\left\{ \!{{{s}_{\rm{r}}}({{f}_{\rm{r}}}) \!\odot \!\left(\!{{{P}}\!_{{{N}}1}} \!\otimes\! {{[{{s}_{\rm{R}}}({{f}_{\rm{R}}}\!) \!\odot\!\! {{s}_{\rm{t}}}({{f}_{\rm{t}}})]}^{\rm T}}\!\right)}\!\! \right\} \\ &+ {{\tilde n}_{\rm{c}}} \end{align}$

      其中,T表示转置符号, ${ρ} \in {\mathbb{C}^{NMK \times 1}}$ 为照射目标对应的散射系数。式(7)中

      $ \begin{align} {{{s}}_{\rm{d}}}({f_{\rm{d}}}) =& \left[ 1,\exp \left\{ {{\rm j}2{{π}} {{f}_{\rm{d}}}} \right\},\exp \left\{ {{\rm j}2{{π}} 2{{f}_{\rm{d}}}} \right\}, ·\!·\!· ,\right.\\ & \quad \left.\exp \left\{ {{\rm j}2{{π}} {{f}_{\rm{d}}}(K - 1)} \right\} \right]^\rm{T} \end{align} $

      式(8)为多普勒频移对应的相位增量。

      $ \begin{align} {{{s}}_{\rm{R}}}({f_{\rm{R}}}) \!=\! &\left[\! 1,\exp \left\{\!\! {\! - {\rm j}2{{π}} \frac{{{R_\sum }}}{c }\Delta f} \! \right\},\right.\\ &\exp \left\{\!\! { \!- {\rm j}2{{π}} \frac{{{R_\sum }}}{c }2\Delta f} \! \right\}, \cdots ,\\ &\left. \exp \left\{ { - {\rm j}2{{π}} \frac{{{R_\sum }}}{c}(M - 1)\Delta f} \right\} \right]^{\rm T} \end{align} $

      式(9)为发射导向矢量中距离导向矢量,其因频率步进量而产生的与距离耦合的相位分量。

      $ \begin{align} {{{s}}_{\rm{t}}}({f_{\rm{t}}}) =& \left[ 1,\exp \left\{ {{\rm j}2{{π}} {{f}_{\rm{t}}}} \right\},\exp \left\{ {{\rm j}2{{π}} 2{{f}_{{\rm{2t}}}}} \right\}, ·\!·\!· ,\right.\\ &\quad \left. \exp \left\{ {{\rm j}2{{π}} (M - 1){f_{(M - 1){\rm{t}}}}} \right\} \right]^{\rm T} \end{align} $

      式(10)为发射导向矢量中发射角度空间角频率导向矢量。

      ${{{s}}_{\rm{at}}}({{f}_{\rm{at}}}) = {{{P}}_{{{N}}{1}}} \otimes {[{{{s}}_{\rm{R}}}({{f}_{\rm{R}}}) \odot {{{s}}_{\rm{t}}}({{f}_{\rm{t}}})]^{\rm T}}$

      式(11)为发射信号的导向矢量。式中 ${{{P}}_{N{1}}} \in {\mathbb{C}^{N \times 1}}$ ,为 $N \times 1$ 维全1矩阵。

      $ \begin{align} {{{s}}_{\rm{r}}}({{f}_{\rm{r}}}) = & \exp \left(\! {{\rm j}2{{π}} \frac{{{d_{\rm{r}}}}}{c}\cos {\psi _{\rm r}}{{[0,1, ·\!·\!· ,N \!- \!1]}^{\rm T}}} \!\right. \\ & \quad \otimes {{[{f_0},{f_1}, ·\!·\!· ,{f_{M - 1}}]}^\rm{T}}\biggr)\\ = & \exp \left\{ {{\rm j}2{{π}} \frac{{{d_{\rm r}}}}{c}\cos {\psi _{\rm r}}{{h}}} \right\} \\ \end{align} $

      其中, ${{h} = \left[ {{{0\ \;\; 0\ ·\!·\!·\ \; \;0}\ \ \; {{f_0}}\ \;{{f_1}}\ ·\!· ·}\ \; {{f_M}}\ \;{2{f_0}}\;\ {2f ·\!· · }}\right.$ $\left. { { {2{f_{M - 1}}}·\!·\!· \, {\left( {N \!-\! 1} \right){f_0}}\ {\left( {N\! - \!1} \right)f}\, { ·\!·\!· {\left( {N \!-\! 1} \right){f_{M \!-\! 1}}}}}} \right]^{\rm T}\, $ , $ {{{s}}_{\rm{r}}}({{f}_{\rm{r}}}) $ $\in {\mathbb{C}^{NM \times 1}}$ , ${{{s}}_{\rm{r}}}({{f}_{\rm{r}}})$ $NM \times 1$ 维接收角度空间角频率导向矢量。

    • 在普通MIMO双基地雷达中,距离模糊杂波因回波信号在时域重叠在一起,故在杂波谱中无法分离,造成杂波谱严重展宽,这是传统的双基地雷达面临的问题。本节将在上一节FDA-MIMO双基地雷达基本的信号模型的基础上,具体阐述其如何根据距离耦合的相位项来实现不同距离的距离模糊杂波在角度和多普勒域进行分离,从而减少距离模糊杂波重叠带来性能下降的影响,并且进一步分析了FDA-MIMO双基地雷达系统参数对分离效果的影响。

      注意到,式(11)中发射信号的导向矢量有距离导向矢量成分和发射角度导向矢量成分,这里去掉PN1项,即:

      ${{{s}}_{{\rm{at}}}}({{f}_{{\rm{at}}}}) = {{{s}}_{\rm{R}}}({{f}_{\rm{R}}}) \odot {{{s}}_{\rm{t}}}({{f}_{\rm{t}}})$

      其中,第1项为FDA的距离导向矢量,第2项为发射角度导向矢量。由式(13)可得:

      $\begin{align} {f_{{\rm{at}}}} =& {f_{\rm{R}}} + {f_{\rm{t}}} = - \frac{{{R_\sum }}}{c}\Delta f + {f_{m{\rm{t}}}} \\ =& - \frac{{{R_\sum }}}{c}\Delta f + \frac{{{d_{\rm{t}}}}}{c}{f_m}\cos {\theta _{\rm{t}}}\cos {\phi _{\rm{t}}} \\ \end{align} $

      观察式(14)可以看出信号快拍中距离相位和发射角度相位是相互非耦合的关系。我们将需要观测的距离范围称为无距离模糊区域,其次是距离模糊区域。模糊区域又可以分为二重模糊区域,三重模糊区域等等,将同一雷达观察角度中,无距离模糊区域和距离模糊区域相应的距离导向可以分别写为:

      $\left.\begin{array}{*{20}{c}} {{f_{{\rm{at}}0}} = - \frac{{{R_{\sum 0}}}}{c}\Delta f + \frac{{{d_{\rm{t}}}}}{c}{f_m}\cos {\theta _{\rm{t}}}\cos {\phi _{\rm{t}}}} \\ {{f_{{\rm{at1}}}} = - \frac{{{R_{\sum 1}}}}{c}\Delta f + \frac{{{d_{\rm{t}}}}}{c}{f_m}\cos {\theta _{\rm{t}}}\cos {\phi _{\rm{t}}}} \\ { ·\!·\!· } \\ {{f_{{\rm{at}}U}} = - \frac{{{R_{\sum {{U}}}}}}{c}\Delta f + \frac{{{d_{\rm{t}}}}}{c}{f_m}\cos {\theta _{\rm{t}}}\cos {\phi _{\rm{t}}}} \end{array}\right\}$

      其中,fat0表示无距离模糊区域对应的导向矢量,fat1表示第1模糊区域即二重杂波模糊区域的导向矢量,R $\sum \! q$ 表示第q模糊区域的发射与接收阵列到目标的距离和,U 表示双基地雷达威力范围内最大距离模糊区域数。下面我们以无模糊和第1模糊区域为例进行分析,已知无模糊区域和第1模糊区域之间相差一个最大不模糊距离Ru,即:

      ${R_{\sum 1}} = {R_{\sum 0}} + {R_u}$

      式中,Ru=c/fr,其中fr为脉冲重复频率,所以无模糊和第1模糊区域导向矢量的相位差可以表示为:

      $\begin{align} \Delta {f_{{\rm{at}}}} =& {f_{{\rm{at}}0}} - {f_{{\rm{at}}1}} \\ =& - \frac{{{R_{\sum 0}}}}{c}\Delta f + \frac{{{d_{\rm{t}}}}}{c}{f\!_m}\cos {\theta _{\rm{t}}}\cos {\phi _{\rm{t}}} \\ & - \left( { - \frac{{{R_u} + {R_{\sum 0}}}}{c}\Delta f + \frac{{{d_{\rm{t}}}}}{c}{f\!_m}\cos {\theta _{\rm{t}}}\cos {\phi _{\rm{t}}}} \right) \\ =& \frac{{{R_u}}}{c}\Delta f = \frac{1}{{{f_{\rm{r}}}}}\Delta f \\ \end{align} $

      由式(17)可以看出尽管无距离模糊区域的杂波和模糊区域的杂波在时间域内混叠在一起,但在空域由于有距离耦合相位项的不同,存在有相位差,由此可以实现在角度域的分离,而这一点是FDA-MIMO双基地雷达区别于传统相控阵和MIMO相控阵双基地雷达的地方,利用该特性,可以实现距离模糊杂波的分离,提高改善因子,完善雷达性能。在发射阵列接收阵列同速,同一高度的情况下,假设接收角度已被扫描估计并在快拍中补偿消除,其杂波谱示意图表示如图3。其中红色,黄色,褐色分别表示无模糊距离,第1模糊距离和第2模糊距离的杂波谱。由图3可以看出FDA-MIMO双基地雷达中因距离导向矢量中不同距离模糊区域存在相位差,所以不同距离模糊区域杂波在杂波谱中相互分离,而MIMO双基地雷达并不能区分,因此距离模糊杂波相互混叠在一起。相位差表示为:

      图  3  双基地雷达杂波谱

      Figure 3.  Bistatic radar clutter spectrum

      $\Delta {f_{{\rm{at}}}} = \frac{{{R_u}}}{c}\Delta f = \frac{{\Delta f}}{{{f_{\rm{r}}}}}$

      观察式(18)可以看出距离模糊杂波分离的程度和 $\Delta $ f 成正比和fr成反比。对于给定距离的无模糊杂波区域,我们可以构建一个FDA-MIMO双基地杂波距离相位补偿导向矢量:

      $ \begin{align} {{{q}}_{\rm{R}}}({{f}_{{\rm{R0}}}}) =& \left[ 1,\exp \left\{ { - {\rm j}2{{π}} \frac{{{R_{\sum 0}}}}{c}\Delta f} \right\},\right.\\ &\quad \exp \left\{ { - {\rm j}2{{π}} \frac{{{R_{\sum 0}}}}{c}2\Delta f} \right\}, ·\!·\!· ,\\ & \quad\left.\exp \left\{ { - {\rm j}2{{π}} \frac{{{R_{\sum 0}}}}{c}(M - 1)\Delta f} \right\} \right]^\rm{T}\quad\quad \quad\quad \end{align} $

      $ \begin{align} {{{e}}_{\rm{R}}}({{f}_{{\rm{R0}}}}) =& {{{P}}_{{{K 1}}}} \otimes {{{P}}_{{{N 1}}}} \otimes \left[ 1,\exp \left\{ { - {\rm j}2{{π}} \frac{{{R_{\sum 0}}}}{c}\Delta f} \right\},\right.\\ & \;\;\;\exp \left\{ { - {\rm j}2{{π}} \frac{{{R_{\sum 0}}}}{c}2\Delta f} \right\}, ·\!·\!· ,\\ & \left. \;\;\;\exp \left\{ { - {\rm j}2{{π}} \frac{{{R_{\sum 0}}}}{c}(M - 1)\Delta f} \right\} \right]^\rm{T} \end{align} $

      式中, ${{{P}}_{{{K 1}}}} \in {\mathbb{C}^{K \times 1}}$ ${{{P}}_{{{N 1}}}} \in {\mathbb{C}^{N \times 1}}$ 为全1矩阵。将式(20)代入空时快拍中进行2次距离杂波补偿,则有:

      $ \begin{align} Y =& \sum\limits_{q = 0}^U {{{{s}}_q} \odot {{{e}}_{\rm{R}}}({{f}_{{\rm{R0}}}})} = \sum\limits_{q = 0}^U {\sum\limits_{i = 1}^{{N_{\rm c}}} {{\varUpsilon _0}{{{ρ }}_{qi}}{{{e}}_{\rm{R}}}({{f}_{{\rm{R0}}}}) \odot } \left\{ {{{{s}}_{\rm{d}}}({f_{{\rm{d}}q}}) \otimes {{\bigr[ {{{{s}}_{\rm{r}}}({{f}_{{\rm{r}}q}}) \odot {{{s}}_{{\rm{at}}}}({f_{{\rm{at}}q}})} \bigr]}^\rm{T}}} \right\}} \\ =& \sum\limits_{q = 0}^U {\sum\limits_{i = 1}^{{N_{\rm c}}} {{\varUpsilon _0}{{{ρ }}_{qi}}} {{{s}}_{\rm{d}}}({f_{{\rm{d}}q}}) \otimes \left\{ {{{{s}}_{\rm{r}}}({{f}_{{\rm{r}}q}}) \odot \left[ {{{{P}}_{{{N 1}}}} \otimes {{\bigr[ {{{{s}}_{\rm{R}}}({{f}_{{\rm{r}}q}} + {{f}_{{\rm{R0}}}}) \odot {{{s}}_{\rm{t}}}({{f}_{{\rm{t}}q}})} \bigr]}^\rm{T}}} \right]} \right\}} \\ =& \sum\limits_{q = 0}^U {\sum\limits_{i = 1}^{{N_{\rm c}}} {{\varUpsilon _0}{{{ρ }}_{qi}}} {{{s}}_{\rm{d}}}({f_{{\rm{d}}q}}) \otimes \left\{ {{{{s}}_{\rm{r}}}({{f}_{{\rm{r}}q}}) \odot \left[ {{{{P}}_{{{N 1}}}} \otimes {{\left[ {{{{s}}_{\rm{R}}}\left( - \frac{{{R_{\sum q}} - {R_{\sum 0}}}}{c}\Delta f\right) \odot {{{s}}_{\rm{t}}}({{f}_{{\rm{t}}q}})} \right]}^\rm{T}}} \right]} \right\}} \\ =& \sum\limits_{q = 0}^u {\sum\limits_{i = 1}^{{N_{\rm c}}} {{\varUpsilon _0}{{{ρ }}_{qi}}} {{{s}}_{\rm{d}}}({f_{{\rm{d}}q}}) \otimes \left\{ {{{{s}}_{\rm{r}}}({{f}_{{\rm{r}}q}}) \odot \left[ {{{{P}}_{{{N 1}}}} \otimes {{\left[ {{{{s}}_{\rm{R}}}\left( - q\frac{{{R_u}}}{c}\Delta f \right) \odot {{{s}}_{\rm{t}}}\left({{f}_{{\rm{t}}q}}\right)} \right]}^\rm{T}}} \right]} \right\}} \end{align} $

      2次距离杂波补偿后,我们所关心的给定距离的无模糊杂波区域q=0,因此被放置在杂波谱的中心,而其他多重距离模糊的杂波被偏置于杂波谱的两旁,从而不影响对所关心区域的观察。整个FDA-MIMO双基地雷达的系统处理流程框图如图4,其中距离模糊区域杂波分离、2次距离杂波补偿和杂波依赖性补偿是本方法的核心步骤。补偿的效果可以从下一节仿真结果中直接观察到。

      图  4  雷达数据处理流程

      Figure 4.  Flow chart of radar data processing

    • 本节将对文中提到的方法进行仿真,首先分析距离模糊杂波和步进频率 $\Delta $ f与脉冲重复频率fr的关系,其次使用IF谱线比较2次距离杂波补偿后的FDA-MIMO双基地雷达与传统的MIMO双基地雷达的性能,最后将本文提出的方法和已有的DW距离模糊杂波抑制方法进行对比。在仿真的讨论中发射与接收阵列均视为等距线阵。由于MIMO双基地雷达杂波谱在发射角度域、接收角度域、多普勒域3维分布,不利于展示说明,因此在本文的仿真中,杂波已经进行过接收角度域的搜索,并进行补偿,图中杂波仅在发射角度域和多普勒域2维分布。双基地雷达接收机与发射机在同一高度和飞行速度,接收与发射阵元均为线性阵列,具有相同阵元数和阵元间距,且飞行速度与基线平行[17],即Vt=Vr, Ht=Hr图5展示了MIMO双基地雷达和未补偿的FDA-MIMO双基地雷达距离模糊下杂波的发射角度域和多普勒频率域2维扫描杂波谱图。仿真参数如表1所示。

      参数 参数值 参数 参数值
      发射阵元数 6 发射机高度 12 km
      接收阵元数 6 接收机高度 12 km
      参考工作频率 1.104 GHz 双基基线长度 50 km
      阵列阵元间距 0.0543 m PRF 2434.8 Hz
      发射机速度 140 m/s 相干脉冲数 8
      目标距离 64.964 km 频率增量 1217.4 Hz
      注:在传统MIMO双基地雷达中没有 $\Delta $f 项或者 $\Delta $f 为0

      表 1  双基地雷达仿真参数

      Table 1.  Simulation parameters of bistatic radar

      图5可以看出FDA-MIMO双基地雷达中因距离导向矢量中不同距离模糊区域存在相位差,所以不同距离模糊区域杂波在杂波谱中相互分离,而MIMO双基地雷达并不能区分因此相互混叠在一起。相位差可表示为:

      $\Delta {f_{{\rm{at}}}} = \frac{{{R_u}}}{c}\Delta f = \frac{{\Delta f}}{{{f_{\rm r}}}}$

      因此距离模糊杂波分离的程度和 $\Delta $ f成正比和fr成反比。

      图5展示了传统MIMO双基地雷达的杂波谱和经过FDA杂波分离后杂波谱的对比情况。图6中绘制了目标区域距离64.964 km的距离杂波空时2维杂波谱和2次距离模糊区域188.174 km的杂波谱,图6(a)图6(c)对应不同 $\Delta $ f 时杂波谱空域上角频率的变化,随着 $\Delta $ f 的逐步增大,距离模糊杂波谱在空域维度向角频率减小的方向平移。具体来看,图6(a) $\Delta $ f=0 Hz,杂波和普通双基地杂波一样,两距离杂波谱分布在空域角频率cos $\psi $ R为0.4~–0.4之间。图6(b) $\Delta $ f=1217.40 Hz,目标距离杂波空域角频率变小,根据式(15)可知,在角度域向下平移R $\sum $ 0/c $\Delta $ f, 2次距离模糊区域杂波向下平移R $\sum $ 1/c $\cdot \Delta $ f,由于空间角频率的周期性,2次距离模糊区域杂波出现在目标距离杂波的上头,根据式(22)两杂波在空域相位之间相位差为 $\Delta $ w= ${\rm{{{π}} }}$ ,从图6(b)也可直观地看出,同一多普勒频率下的两杂波在角度域的间隔长度为整个2 ${\rm{{{π}} }}$ 空间角度频率域的一半。图6(c) $\Delta $ f=2434.80 Hz,此时的 $\Delta $ ffr相同,目标距离杂波和2次距离模糊杂波在角度域向下分别平移R $\sum $ 0/c $\Delta $ fR $\sum $ 1/c $\Delta $ f,而两杂波在空域频率之间相位差为 $\Delta $ w=2 ${\rm{{{π}} }}$ ,这时虽然平移的相位不同但由于相差为2 ${\rm{{{π}} }}$ ,所以在周期的空间角频率中两杂波又重叠在一起,如图6(c)所示,只不过和图6(a) $\Delta $ f =0 Hz相比,两杂波有一定的偏移。由上面的分析可知,杂波在空域角频率维度的间距 $\Delta $ w的大小和 $\Delta $ ffr的比值成正比,由于空域角频率的周期性,在 $w \in [ - {{π}} ,{{π}} ]$ 的范围内,在 $\Delta $ w = ${\rm{{{π}} }}$ 时,两杂波分离间距最大,当 $\Delta $ w< ${\rm{{{π}} }}$ 时,杂波距离靠近,当 $\Delta $ w> ${\rm{{{π}} }}$ 时,角频率周期性使得2次杂波移动至目标上方,同样相位间距小于 ${\rm{{{π}} }}$ ,所以当 $\Delta $ w = ${\rm{{{π}} }}$ , $\Delta $ f 取 1/2 fr时,实现了两杂波的最大分离,利用此特性将距离模糊杂波进行两侧分离,从而减小了杂波展宽,实现了杂波抑制。文献[18]也对此作出了分析。

      图  6  $\Delta f$ 不同值时对应的杂波谱

      Figure 6.  The corresponding clutter spectrum of $\Delta f$ at different values

      图  5  双基地雷达杂波谱图

      Figure 5.  Bistatic radar clutter spectrum

      下面将展示说明FDA-MIMO双基地雷达和传统MIMO双基地雷达的区别,并将其效果与DW杂波抑制方法作对比。DW杂波模糊抑制的基本思路是以雷达波束主瓣照射方向的待测单元样本为基准,将相邻的非均匀训练样本在主瓣照射方向上进行多普勒维的平移。其使在主瓣方向的多普勒杂波重叠在一点,从而抑制距离耦合效应。

      图7(a)表示了FDA-MIMO双基地雷达由于距离相位的引入,实现了不同距离的距离模糊杂波的分离,图7(b)则是加入了构造的杂波距离相位补偿导向矢量 ${{{e}}_{\rm R}}({f_{{\rm R}0}})$ ,补偿后的杂波谱将我们关心的无模糊距离的杂波谱放置在观察的中心。将图7(b)补偿后的FDA-MIMO双基地雷达杂波谱和图7(c)距离模糊杂波未分离前的MIMO双基地雷达杂波谱作对比,可以看出图7(b)中不同重数的距离模糊杂波在谱中实现分离,因其杂波强度和距离成反比所以所需的无模糊区域距离杂波强度最大在中心零度周围,第1模糊区域的二重距离模糊杂波和第2模糊区域的三重距离模糊杂波被分离在两处。而图7(c)中并未实现分离。

      图7(b)图7(c)可以得到0°角频率处多普勒IF曲线图。IF谱线是用于评价雷达效能的一项指标,其定义为输出端和输入端信噪功率比的比值:

      $ \begin{align} {\rm{IF}} =& \frac{{p_{\rm{s}}^{{\rm{out}}}}}{{p_{\rm{n}}^{{\rm{out}}}}}\bigg/\frac{{p_{\rm{s}}^{{\rm{in}}}}}{{p_{\rm{n}}^{{\rm{in}}}}} = \frac{{{{{w}}^*}{{s}}{{{s}}^*}{{w}}}}{{{{{w}}^*}{{Qw}}}}\bigg/\frac{{{{{s}}^*}{{s}}}}{{{\rm{tr}}({{Q}})}} \\ =& \frac{{{{{w}}^*}{{s}}{{{s}}^*}{{w}} \;{\bf \cdot} \;{\rm{tr}}({{Q}})}}{{{{{w}}^*}{{Qw}} \;{\bf \cdot}\; {{{s}}^*}{{s}}}} \end{align}$

      其中,s表示信号矢量,w表示接收端导向矢量,Q为噪声协方差矩阵。

      图8图7(b)图7(c)的0°角频率处多普勒IF曲线图进行对比,图中红色为FDA-MIMO双基地雷达,蓝色为MIMO双基地雷达,由于MIMO双基地雷达杂波具有距离依赖性所以其在空间角度和多普勒频率域杂波谱展开严重,从图8看出传统MIMO双基地雷达有明显的两处凹陷,分别在归一化多普勒频率的0.17和–0.22处,–0.22处凹陷源于无模糊杂波区域的杂波凹陷,而0.17处源于第1模糊杂波区域的杂波凹陷,由于杂波的距离依赖性,两处杂波凹陷不在同一位置,相互叠加形成了蓝色IF谱线的性能衰减,而FDA-MIMO双基地雷达通过2次距离杂波补偿将不同模糊区域的杂波分离,在观察处仅有一处源于无模糊区域杂波的–0.22频率的凹陷,因此性能好于MIMO双基地雷达。

      图  7  双基地雷达杂波谱

      Figure 7.  Bistatic radar clutter spectrum

      图  8  雷达杂波IF特性

      Figure 8.  IF characteristics of radar clutter

      对比完,FDA-MIMO和传统MIMO双基地雷达后,下面将FDA-MIMO双基地的杂波抑制效果和使用DW的MIMO双基地杂波抑制效果进行对比。仿真中频率增量取 $\Delta $ f = 811.6 Hz。图9表示了MIMO双基地雷达杂波脊在接收角频域,发射角频域和归一化多普勒频率域的3维关系图,图中可以看出,双基地杂波具有距离依赖性,随着距离的变化杂波关系分布为非均匀。在下面的讨论中,假设接收角度已被扫描估计且补偿过,杂波的讨论在发射角度和多普勒2维进行。

      图  9  双基地MIMO雷达杂波脊

      Figure 9.  Bistatic MIMO radar clutter ridge

      图10(a)图10(b)展示了FDA-MIMO双基地雷达的杂波在DW补偿前后杂波谱的改变情况,DW补偿使得杂波谱收窄,尤其在wt=–0.8处到wt=0处能明显观察到杂波的变化。

      图  10  FDA-MIMO双基地雷达DW补偿杂波谱

      Figure 10.  FDA-MIMO bistatic radar clutter spectrum of DW compensation method

      图11分别展示了杂波谱MIMO双基地雷达IF图线(紫色)、MIMO双基地雷达经过DW补偿后IF图线(绿色)、FDA-MIMO双基地雷达IF图线(蓝色)及FDA-MIMO双基地雷达经过DW补偿后IF图线(红色)。可以看出和FDA-MIMO双基地雷达比较,首先,由于距离模糊杂波的距离依赖性,传统MIMO双基地雷达杂波谱在多普勒和空间角度2维域中杂波展宽严重,其使得目标检测性能下降。通过图11的IF谱线可以看出,FDA-MIMO双基地雷达检测性能优于传统MIMO双基地雷达。此外,由于没有实现多重距离模糊区域杂波的分离,MIMO双基地IF谱线在–0.48和–0.31两处存在凹陷,而FDA-MIMO双基地雷达仅出现一处凹陷。对比雷达杂波DW补偿前后的IF谱线,和补偿前比较,无论是FDA-MIMO还是MIMO双基地雷达,其IF谱线均有明显的收窄,杂波谱IF性能得到改善。

      图  11  FDA-MIMO双基地雷达IF谱线对比图

      Figure 11.  Comparison of FDA-MIMO bistatic radar IF lines

    • 本文针对MIMO双基地雷达存在的距离杂波模糊问题提出了新的解决方法,将距离耦合相位分量引入双基地雷达中,使用FDA技术实现了双基地MIMO雷达距离模糊杂波的抑制,明显地改善了MIMO双基地雷达因距离模糊而产生的非均匀杂波谱展宽严重的现象。并且文中将FDA-MIMO双基地雷达技术和已有的DW技术进行比较,使用IF指标,对比了两种技术在杂波抑制上性能的表现。FDA-MIMO技术改善效果优于DW的性能,最后将FDA-MIMO双基地雷达和DW结合,使用FDA-MIMO双基地雷达再进行DW处理,实现了一种对MIMO双基地雷达杂波更优的抑制方法。

参考文献 (18)

目录

    /

    返回文章
    返回