基于贝叶斯的高斯杂波背景下MIMO雷达自适应检测算法

韩金旺 张子敬 刘军 赵永波

引用本文:
Citation:

基于贝叶斯的高斯杂波背景下MIMO雷达自适应检测算法

    作者简介: 韩金旺(1993–),男,河北石家庄人,西安电子科技大学硕士生,研究方向为MIMO雷达信号处理。E-mail: jwhan0828@163.com;张子敬(1967–),男,西安电子科技大学教授,博士生导师,IEEE会员。研究方向为雷达目标检测、雷达极化信息处理、多速率滤波器组设计等。E-mail: zjzhang@xidian.edu.cn;刘 军(1983–),男,现为中国科学技术大学信息技术学院副教授,中国电子学会高级会员。研究方向为雷达信号处理和机器学习。E-mail: junliu@ustc.edu.cn;赵永波(1972–),男,现为西安电子科技大学教授,中国电子学会高级会员。研究方向为新体制雷达系统和MIMO雷达。E-mail: ybzhao@xidian.edu.cn.
    通讯作者: 刘军, junliu@ustc.edu.cn
  • 基金项目:

    国家自然科学基金(61871469, 61571349),陕西省自然科学基金(2018JM6051)

  • 中图分类号: TN957.51

Adaptive Bayesian Detection for MIMO Radar in Gaussian Clutter

    Corresponding author: LIU Jun, junliu@ustc.edu.cn ;
  • Fund Project: The National Natural Science Foundation of China (61871469, 61571349), The Natural Science Foundation of Shaanxi Province (2018JM6051)

    CLC number: TN957.51

  • 摘要: 对于集中式多输入多输出(Multiple-Input Multiple-Output, MIMO)雷达,该文研究了高斯杂波背景下的目标检测问题。该文假设杂波的协方差矩阵是未知随机的,且服从逆复Wishart分布,基于贝叶斯方法和广义似然比检验准则设计了两种新型自适应检测器。该文提出的贝叶斯检测器具有两个显著的优点:(1)不需要训练数据;(2)杂波的先验知识体现在设计方案中,从而提高了检测性能。仿真结果显示该文提出的贝叶斯检测器的检测性能优于目前常用的非贝叶斯检测器,特别是在发射波形采样数较少时。另外,该贝叶斯检测器在参数失配条件下的性能会有一定程度下降。
  • 图 1  K 取不同值时的ROC曲线

    Figure 1.  ROC curves for different K

    图 2  K 取不同值时${{{P}}_{\rm{d}}}$关于SCR的曲线

    Figure 2.  ${{{P}}_{\rm{d}}}$ versus SCR for different K

    图 3  L 取不同值时检测器的ROC曲线

    Figure 3.  ROC curves for different L

    图 4  1阶迟滞相关系数失配时检测器的ROC曲线

    Figure 4.  ROC curves when one-lag correlation coefficient is mismatched

    图 5  自由度L失配时检测器的ROC曲线

    Figure 5.  ROC curves when L is mismatched

  • [1] LI Jian and STOICA P. MIMO Radar Signal Processing[M]. Hoboken, NJ: John Wiley & Sons, 2009.
    [2] TANG Bo, NAGHSH M M, and TANG Jun. Relative entropy-based waveform design for MIMO radar detection in the presence of clutter and interference[J]. IEEE Transactions on Signal Processing, 2015, 63(14): 3783–3796 doi: 10.1109/TSP.2015.2423257
    [3] 周伟, 刘永祥, 黎湘, 等. MIMO-SAR技术发展概况及应用浅析[J]. 雷达学报, 2014, 3(1): 10–18 doi: 10.3724/SP.J.1300.2013.13074ZHOU Wei, LIU Yongxiang, LI Xiang, et al. Brief analysis on the development and application of multi-input multi-output synthetic aperture radar[J]. Journal of Radars, 2014, 3(1): 10–18 doi: 10.3724/SP.J.1300.2013.13074
    [4] 王珽, 赵拥军, 胡涛. 机载MIMO雷达空时自适应处理技术研究进展[J]. 雷达学报, 2015, 4(2): 136–148 doi: 10.12000/JR14091WANG Ting, ZHAO Yongjun, and HU Tao. Overview of space-time adaptive processing for airborne MIMO radar[J]. Journal of Radars, 2015, 4(2): 136–148 doi: 10.12000/JR14091
    [5] HAIMOVICH A M, BLUM R S, and CIMINI L J. MIMO radar with widely separated antennas[J]. IEEE Signal Processing Magazine, 2008, 25(1): 116–129 doi: 10.1109/MSP.2008.4408448
    [6] ZHOU Shenghua, LIU Hongwei, ZHAO Yongbo, et al. Target spatial and frequency scattering diversity property for diversity MIMO radar[J]. Signal Processing, 2011, 91(2): 269–276 doi: 10.1016/j.sigpro.2010.07.004
    [7] TANG Bo and TANG Jun. Joint design of transmit waveforms and receive filters for MIMO radar space-time adaptive processing[J]. IEEE Transactions on Signal Processing, 2016, 64(18): 4707–4722 doi: 10.1109/TSP.2016.2569431
    [8] TANG Bo, ZHANG Yu, and TANG Jun. An efficient minorization maximization approach for MIMO radar waveform optimization via relative entropy[J]. IEEE Transactions on Signal Processing, 2018, 66(2): 400–411 doi: 10.1109/TSP.2017.2771726
    [9] 金镇, 谢良贵, 文树梁. 分布式MIMO雷达单脉冲测角[J]. 雷达学报, 2014, 3(4): 474–479 doi: 10.3724/SP.J.1300.2014.13077JIN Zhen, XIE Lianggui, and WEN Shuliang. Distributed MIMO radar monopulse angular estimation[J]. Journal of Radars, 2014, 3(4): 474–479 doi: 10.3724/SP.J.1300.2014.13077
    [10] LI Jian and STOICA P. MIMO radar with colocated antennas[J]. IEEE Signal Processing Magazine, 2007, 24(5): 106–114 doi: 10.1109/MSP.2007.904812
    [11] AUBRY A, DE MAIO A, and HUANG Yongwei. MIMO radar beampattern design via PSL/ISL optimization[J]. IEEE Transactions on Signal Processing, 2016, 64(15): 3955–3967 doi: 10.1109/TSP.2016.2543207
    [12] CUI Guolong, YU Xianxiang, CAROTENUTO V, et al. Space-time transmit code and receive filter design for colocated MIMO radar[J]. IEEE Transactions on Signal Processing, 2017, 65(5): 1116–1129 doi: 10.1109/TSP.2016.2633242
    [13] CUI Guolong, LI Hongli, and RANGASWAMY M. MIMO radar waveform design with constant modulus and similarity constraints[J]. IEEE Transactions on Signal Processing, 2014, 62(2): 343–353 doi: 10.1109/TSP.2013.2288086
    [14] 梁浩, 崔琛, 余剑. 十字型阵列MIMO雷达高精度二维DOA估计[J]. 雷达学报, 2016, 5(3): 254–264 doi: 10.12000/JR16016LIANG Hao, CUI Chen, and YU Jian. Two-dimensional DOA estimation with high accuracy for MIMO radar using cross array[J]. Journal of Radars, 2016, 5(3): 254–264 doi: 10.12000/JR16016
    [15] HE Qian, LEHMANN N H, BLUM R S, et al. MIMO radar moving target detection in homogeneous clutter[J]. IEEE Transactions on Aerospace and Electronic Systems, 2010, 46(3): 1290–1301 doi: 10.1109/TAES.2010.5545189
    [16] 程子扬, 何子述, 王智磊, 等. 分布式MIMO雷达目标检测性能分析[J]. 雷达学报, 2017, 6(1): 81–89 doi: 10.12000/JR16147CHENG Ziyang, HE Zishu, WANG Zhilei, et al. Detection Performance analysis for distributed MIMO radar[J]. Journal of Radars, 2017, 6(1): 81–89 doi: 10.12000/JR16147
    [17] HE Qian, BLUM R S, and HAIMOVICH A M. Noncoherent MIMO radar for location and velocity estimation: More antennas means better performance[J]. IEEE Transactions on Signal Processing, 2010, 58(7): 3661–3680 doi: 10.1109/TSP.2010.2044613
    [18] FUHRMANN D R and SAN ANTONIO G. Transmit beamforming for MIMO radar systems using signal cross-correlation[J]. IEEE Transactions on Aerospace and Electronic Systems, 2008, 44(1): 171–186 doi: 10.1109/TAES.2008.4516997
    [19] BEKKERMAN I and TABRIKIAN J. Target detection and localization using MIMO radars and sonars[J]. IEEE Transactions on Signal Processing, 2006, 54(10): 3873–3883 doi: 10.1109/TSP.2006.879267
    [20] XU Luzhou, LI Jian, and STOICA P. Target detection and parameter estimation for MIMO radar systems[J]. IEEE Transactions on Aerospace and Electronic Systems, 2008, 44(3): 927–939 doi: 10.1109/TAES.2008.4655353
    [21] DE MAIO A, LOPS M, and VENTURINO L. Diversity-integration tradeoffs in MIMO detection[J]. IEEE Transactions on Signal Processing, 2008, 56(10): 5051–5061 doi: 10.1109/TSP.2008.928693
    [22] CUI Guolong, KONG Lingjiang, YANG Xiaobo, et al. The Rao and Wald tests designed for distributed targets with polarization MIMO radar in compound-Gaussian clutter[J]. Circuits, Systems, and Signal Processing, 2012, 31(1): 237–254 doi: 10.1007/s00034-010-9250-0
    [23] ZHANG Tianxian, CUI Guolong, KONG Lingjiang, et al. Adaptive Bayesian detection using MIMO radar in spatially heterogeneous clutter[J]. IEEE Signal Processing Letters, 2013, 20(6): 547–550 doi: 10.1109/LSP.2013.2255272
    [24] LIU Jun, LI Hongbin, and HIMED B. Persymmetric adaptive target detection with distributed MIMO radar[J]. IEEE Transactions on Aerospace and Electronic Systems, 2015, 51(1): 372–382 doi: 10.1109/TAES.2014.130652
    [25] LI Jian, XU Luzhou, STOICA P, et al. Range compression and waveform optimization for MIMO radar: A CramÉr-Rao bound based study[J]. IEEE Transactions on Signal Processing, 2008, 56(1): 218–232 doi: 10.1109/TSP.2007.901653
    [26] LIU Weijian, WANG Yongliang, LIU Jun, et al. Adaptive detection without training data in colocated MIMO radar[J]. IEEE Transactions on Aerospace and Electronic Systems, 2015, 51(3): 2469–2479 doi: 10.1109/TAES.2015.130754
    [27] LIU Jun, ZHOU Shenghua, LIU Weijian, et al. Tunable adaptive detection in colocated MIMO radar[J]. IEEE Transactions on Signal Processing, 2018, 66(4): 1080–1092 doi: 10.1109/TSP.2017.2778693
    [28] ZHU Xumin, LI Jian, and STOICA P. Knowledge-aided space-time adaptive processing[J]. IEEE Transactions on Aerospace and Electronic Systems, 2011, 47(2): 1325–1336 doi: 10.1109/TAES.2011.5751261
    [29] WANG Yikai, XIA Wei, and HE Zishu. CFAR knowledge-aided radar detection with heterogeneous samples[J]. IEEE Signal Processing Letters, 2017, 24(5): 693–697 doi: 10.1109/LSP.2017.2688386
    [30] GAO Yongchan, LI Hongbin, and HIMED B. Knowledge-aided range-spread target detection for distributed MIMO radar in nonhomogeneous environments[J]. IEEE Transactions on Signal Processing, 2017, 65(3): 617–627 doi: 10.1109/TSP.2016.2625266
    [31] KONG Lingjiang, LI Na, CUI Guolong, et al. Adaptive Bayesian detection for multiple-input multiple-output radar in compound-gaussian clutter with random texture[J]. IET Radar, Sonar & Navigation, 2016, 10(4): 689–698 doi: 10.1049/iet-rsn.2015.0241
    [32] DE MAIO A, FARINA A, and FOGLIA G. Knowledge-aided Bayesian radar detectors & their application to live data[J]. IEEE Transactions on Aerospace and Electronic Systems, 2010, 46(1): 170–183 doi: 10.1109/TAES.2010.5417154
    [33] BILLINGSLEY J B, FARINA A, GINI F, et al. Statistical analyses of measured radar ground clutter data[J]. IEEE Transactions on Aerospace and Electronic Systems, 1999, 35(2): 579–593 doi: 10.1109/7.766939
    [34] MELVIN W L and SHOWMAN G A. An approach to knowledge-aided covariance estimation[J]. IEEE Transactions on Aerospace and Electronic Systems, 2006, 42(3): 1021–1042 doi: 10.1109/TAES.2006.248216
  • [1] 丁昊王国庆刘宁波关键 . 逆Gamma纹理背景下两类子空间目标的自适应检测方法. 雷达学报, doi: 10.12000/JR16088
    [2] 侯育星徐刚 . InSAR通道联合稀疏贝叶斯特征化成像. 雷达学报, doi: 10.12000/JR18100
    [3] 谷文堃王党卫马晓岩 . 分布式MIMO-ISAR子图像融合方法. 雷达学报, doi: 10.12000/JR16042
    [4] 丁昊薛永华黄勇关键 . 均匀和部分均匀杂波中子空间目标的斜对称自适应检测方法. 雷达学报, doi: 10.12000/JR14133
    [5] 邹鲲吴德伟张斌李伟 . 一种针对失配信号的可调自适应检测器. 雷达学报, doi: 10.12000/JR14129
    [6] 李伟王兴亮邹鲲徐艺盟张群 . 基于数据融合和陷波滤波的MIMO 雷达抗欺骗干扰算法 (英文). 雷达学报, doi: 10.3724/SP.J.1300.2012.20060
    [7] 王宇卓朱圣棋许京伟 . MIMO双基地机载雷达距离模糊杂波抑制方法. 雷达学报, doi: 10.12000/JR18016
    [8] 周明马亮王宁杨予昊 . 面向海面目标检测的陆海分离和海面分区算法研究. 雷达学报, doi: 10.12000/JR19036
    [9] 钟金荣文贡坚 , . 基于块稀疏贝叶斯学习的雷达目标压缩感知(英文). 雷达学报, doi: 10.12000/JR15056
    [10] 曾铮张福博陈龙永卜祥玺周思言 . 基于多输入多输出合成孔径雷达的二维混合基线抗欺骗干扰方法. 雷达学报, doi: 10.12000/JR18118
    [11] 赵军香梁兴东李焱磊 . 一种基于似然比统计量的SAR相干变化检测. 雷达学报, doi: 10.12000/JR16065
    [12] 黄剑胡卫东 . 基于贝叶斯框架的空间群目标跟踪技术. 雷达学报, doi: 10.3724/SP.J.1300.2012.20079
    [13] 王永良刘维建谢文冲段克清高飞王泽涛 . 机载雷达空时自适应检测方法研究进展. 雷达学报, doi: 10.3724/SP.J.1300.2014.13081
    [14] 闫敏韦顺军田博坤张晓玲师君 . 基于稀疏贝叶斯正则化的阵列SAR高分辨三维成像算法. 雷达学报, doi: 10.12000/JR18067
    [15] 黄瑞杜小勇胡卫东 . OFDM雷达多目标运动参数的近似最大似然估计. 雷达学报, doi: 10.12000/JR17116
    [16] 陈希信 . 天波雷达后多普勒自适应波束形成. 雷达学报, doi: 10.12000/JR15124
    [17] 刘维建王利才狄源水简涛谢谠王永良 . 自适应能量检测器及在失配信号检测中的应用(英文). 雷达学报, doi: 10.12000/JR14132
    [18] 王珽赵拥军胡涛 . 机载MIMO雷达空时自适应处理技术研究进展. 雷达学报, doi: 10.12000/JR14091
    [19] 胡锡坤金添 . 基于自适应小波尺度选择的生物雷达呼吸与心跳分离方法. 雷达学报, doi: 10.12000/JR16103
    [20] 孙铭才张秦陈光陆 . 动态时间窗下的相控阵雷达自适应调度算法. 雷达学报, doi: 10.12000/JR17104
  • 加载中
图(5)
计量
  • 文章访问数:  358
  • HTML浏览量:  113
  • PDF下载量:  49
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-10-25
  • 录用日期:  2019-01-03

基于贝叶斯的高斯杂波背景下MIMO雷达自适应检测算法

    通讯作者: 刘军, junliu@ustc.edu.cn
    作者简介: 韩金旺(1993–),男,河北石家庄人,西安电子科技大学硕士生,研究方向为MIMO雷达信号处理。E-mail: jwhan0828@163.com;张子敬(1967–),男,西安电子科技大学教授,博士生导师,IEEE会员。研究方向为雷达目标检测、雷达极化信息处理、多速率滤波器组设计等。E-mail: zjzhang@xidian.edu.cn;刘 军(1983–),男,现为中国科学技术大学信息技术学院副教授,中国电子学会高级会员。研究方向为雷达信号处理和机器学习。E-mail: junliu@ustc.edu.cn;赵永波(1972–),男,现为西安电子科技大学教授,中国电子学会高级会员。研究方向为新体制雷达系统和MIMO雷达。E-mail: ybzhao@xidian.edu.cn
  • ①. 西安电子科技大学雷达信号处理国家重点实验室   西安   710071
  • ②. 中国科学技术大学信息科学技术学院   合肥   230027
基金项目:  国家自然科学基金(61871469, 61571349),陕西省自然科学基金(2018JM6051)

摘要: 对于集中式多输入多输出(Multiple-Input Multiple-Output, MIMO)雷达,该文研究了高斯杂波背景下的目标检测问题。该文假设杂波的协方差矩阵是未知随机的,且服从逆复Wishart分布,基于贝叶斯方法和广义似然比检验准则设计了两种新型自适应检测器。该文提出的贝叶斯检测器具有两个显著的优点:(1)不需要训练数据;(2)杂波的先验知识体现在设计方案中,从而提高了检测性能。仿真结果显示该文提出的贝叶斯检测器的检测性能优于目前常用的非贝叶斯检测器,特别是在发射波形采样数较少时。另外,该贝叶斯检测器在参数失配条件下的性能会有一定程度下降。

English Abstract

    • 近年来多输入多输出(Multiple-Input Multiple-Output, MIMO)雷达是雷达领域的一个研究热点[14]。根据天线配置,MIMO雷达分为两类:一类是分布式MIMO雷达,其特点是发射和接收天线被分置于不同的空间位置[59],另一类是集中式MIMO雷达[1014],其将发射和接收天线集中放置。分布式MIMO雷达的空间分集性可以有效减轻目标的雷达横截面积(Radar Cross Sections, RCSs)闪烁,从而可以提高检测性能[15,16]和估计精度[17]。集中式MIMO雷达发射相互独立或正交的波形,从而获得了波形分集性能,提高了系统自由度,因此可以获得更优异的性能,如更加灵活的波束形成性能[18]、更优的检测性能[19]和更好的干扰抑制性能[20]

      雷达的主要功能之一是目标检测。许多文献对分布式MIMO雷达在协方差相关性未知情况下的检测问题进行了研究[2123]。为了估计未知的杂波协方差矩阵,一般要遵循如下假设:有充足的均匀的训练数据可用。实际中,由于杂波谱的快变性,往往无法得到充足的均匀性训练数据,由此会带来检测器的性能下降。可以采用一些方法来减轻性能损失,例如,利用杂波协方差矩阵的斜对称结构[24]。然而,由均匀性训练数据不足引起的问题不能被完全解决,因为这些方法仍然需要训练数据。

      为了彻底解决均匀性训练数据不足的问题,一个可能的方法是不使用训练数据进行目标检测。集中式MIMO雷达可以直接使用自适应处理技术而不需要训练数据[25]。本文研究的检测问题主要针对集中式MIMO雷达。除非特别说明,以下将集中式MIMO雷达简称为MIMO雷达。

      文献[20]从MIMO雷达参数估计和干扰抑制的角度,给出了一种基于广义似然比检验 (Generalized Likelihood Ratio Test, GLRT)的检测器,文献[26]推导出这种检测器虚警概率和检测概率的闭集表达式并分析其性能,注意到这种检测器不需要训练数据,从而有效解决了均匀训练数据匮乏的问题。根据Rao和 Wald准则,文献[26]设计了两个不需训练数据的自适应检测器。几种检测器各有优势,在信号匹配情况下,文献[20]中的GLRT检测器的检测性能最好。在信号失配情况下,与GLRT检测器相比,Rao检测器抑制失配信号的性能更好,而Wald检测器的鲁棒性更好。进一步地,文献[27]推导了一种参数可调的MIMO雷达检测器,从而可以灵活地调节选择性和鲁棒性。

      需要指出的是,许多检测器将杂波协方差矩阵建模为未知确定的矩阵[20,26,27],此方法的一个缺点是无法将协方差矩阵的先验信息应用到自适应检测算法中。贝叶斯方法将杂波协方差矩阵建模为未知随机的,并为其指定一个合适的先验分布,从而可以有效地利用协方差矩阵的先验知识。

      贝叶斯方法已经应用到了一些场景,如多通道雷达的空时自适应处理[28,29]和高斯杂波背景[30]或复合高斯杂波背景[31]下的分布式MIMO雷达的自适应检测。对于分布式MIMO雷达,文献[30]针对距离扩展目标推导了一种基于贝叶斯理论的GLRT检测器,雷达的检测性能得到了较大提高。文献[31]研究了随机纹理的复合高斯杂波背景下的分布式MIMO雷达检测问题。然而,目前还没有基于贝叶斯理论的集中式MIMO雷达检测器。

      本文借鉴在分布式MIMO雷达中基于贝叶斯框架的GLRT检测器的理念,结合集中式MIMO雷达的特点,研究与设计集中式MIMO雷达的贝叶斯检测器。本文假设杂波协方差矩阵服从逆复Wishart分布,因为这种分布可以满足数学上的易处理性而且可以反映雷达工作的环境。根据GLRT准则,本文利用杂波协方差矩阵的先验分布设计了两种贝叶斯检测器。仿真结果显示相比传统的非贝叶斯检测器,本文的贝叶斯检测器具有更好的检测性能。

      文中(·)T,(·)*和(·)H分别表示转置、复共轭和复共轭转置,符号“~”表示“服从分布于”,${\mathcal{W}^{\ - 1}}$表示逆复Wishart分布,$\Gamma \left( \cdot \right)$表示Gamma函数。det(·)和tr(·)分别表示矩阵的行列式和矩阵的迹。|·|表示复数的模。j表示虚数单位,${\mathbb{C}^{m \times n}}$表示m×n维的复矩阵集合,I 表示单位阵。

    • 假设MIMO雷达发射和接收天线数分别为MN,令atar分别表示发射和接收导向矢量,第m个发射天线发射的波形记为${s_m}(k),k = 1,2, ·\!·\!· ,K$,其中K是样本数。这里定义

      $ {{s}_k} = {[{s_1}(k)\ {s_2}(k) \ ·\!·\!·\ {s_M}(k)]^{\rm{T}}} $

      在窄带模型下,目标位置处的基带信号可以表示为

      $ \sum\limits_{m = 1}^M {{{\rm{e}}^{ - {\rm{j}}2{{π}}{f_0}{\tau _m}(\theta )}}} {s_m}(k) = {a}_{\rm{t}}^{\rm{H}}{{s}_k},\;\;k = 1,2, ·\!·\!· ,K $

      其中,${f_0}$是载频,${\tau _m}(\theta )$是第$m$个发射天线发射的信号到达目标所需的时间。发射导向矢量为

      $ {{a}_{\rm{t}}} = {\left[ {\begin{array}{*{20}{c}} {{{\rm{e}}^{{\rm{j}}2{{π}}{f_0}{\tau _1}(\theta )}}}&{{{\rm{e}}^{{\rm{j}}2{{π}}{f_0}{\tau _2}(\theta )}}}& ·\!·\!· &{{{\rm{e}}^{{\rm{j}}2{{π}}{f_0}{\tau _M}(\theta )}}} \end{array}} \right]^{\rm{T}}} $

      ${x_n}(k)$表示第$n$个接收天线收到的信号,并且定义

      $ \begin{align} {{x}_k} =& {[\begin{array}{*{20}{c}} {{x_1}(k)}&{{x_2}(k)}& ·\!·\!· &{{x_N}(k)} \end{array}]^{\rm{T}}},\\ {\rm{}}&\quad k = 1,2, ·\!·\!· ,K \end{align} $

      接收导向矢量为

      $ {{a}_{\rm{r}}} = {\left[ {\begin{array}{*{20}{c}} {{{\rm{e}}^{{\rm{j}}2{{π}}{f_0}{{\tilde \tau }_1}(\theta )}}}&{{{\rm{e}}^{{\rm{j}}2{{π}}{f_0}{{\tilde \tau }_2}(\theta )}}}& ·\!·\!· &{{{\rm{e}}^{{\rm{j}}2{{π}}{f_0}{{\tilde \tau }_N}(\theta )}}} \end{array}} \right]^{\rm{T}}} $

      其中,${\tilde \tau _n}(\theta )$是被目标反射的信号从目标位置到达第$n$个接收阵元所需的时间。基于点目标假设,接收的数据矢量可以表示为

      $ {{x}_k} = \alpha {{a}_{\rm{r}}}{a}_{\rm{t}}^{\rm{H}}{{s}_k} + {{v}_k},\;\;k = 1,2, ·\!·\!· ,K $

      其中,$\alpha $是未知确定的复幅度,与目标散射系数和传播损耗有关,${{{v}}_k}$代表干扰和噪声。进一步地,矩阵形式的接收信号为

      ${X} = \alpha {{a}_{\rm{r}}}{a}_{\rm{t}}^{\rm{H}}{S} + {V}$

      其中,

      $ {X} = [\!\!\begin{array}{*{20}{c}} {{x}(1)}&{{x}(2)}& ·\!·\!· &{{x}(K)} \end{array}\!\!] \in {\mathbb{C}^{N \times K}} $

      $ {S} = [\!\!\begin{array}{*{20}{c}} {{s}(1)}&{{s}(2)}& ·\!·\!· &{{s}(K)} \end{array}\!\!] \in {\mathbb{C}^{M \times K}}\ \ $

      $ {V} = [\!\!\begin{array}{*{20}{c}} {{v}(1)}&{{v}(2)}& ·\!·\!· &{{v}(K)} \end{array}\!\!] \in {\mathbb{C}^{N \times K}} $

      与文献[19,20]相同,这里假设${V}$的各列相互独立且服从均值为零、协方差矩阵为${{R}}$的高斯分布。

      目标检测问题可以转化为式(11)的二元假设检验问题:

      $\left. {\begin{array}{*{20}{l}} {{\rm{H}_0}:{{X}} = {{V}}} \\ {{\rm{H}_1}:{{X}} = \alpha {{{a}}_{\rm{r}}}{{a}}_{\rm{t}}^{\rm{H}}{{S}} + {{V}}} \end{array}} \right\}$

      其中,${\rm{H}_1}$表示目标存在的假设,${\rm{H}_0}$表示目标不存在的假设。

      为解决上述假设检验问题,所有已存在的方法是将协方差矩阵${{R}}$建模为未知确定的参数[20,26,27],这种模型没有采用任何杂波的先验信息,并且要求样本数满足$K \ge N + 1$。本文采用贝叶斯方法,即认为协方差矩阵${{R}}$是未知随机的矩阵,并为其指定一个先验分布。

      一个合适的先验分布要能反映雷达的工作环境,并且在数学上易处理,考虑以上两点,本文认为杂波协方差矩阵${{R}}$服从逆复Wishart分布[32],即

      $ {{R}} \sim {\mathcal{W}^{ \ - 1}}(L,L{{{R}}_0}) $

      其中,$L$是自由度,${{{R}}_0}$$N \times N$维的正定矩阵,并且$L \ge N$。协方差矩阵${{R}}$的概率密度函数为

      $ {f_{{R}}}({{R}}) = \frac{{{{\det }^{L + N}}\left({{{R}}^{ - 1}}\right)}}{{G(L,L{{{R}}_0})}}\exp \left[ - {\rm tr}\left({{{R}}^{ - 1}}L{{{R}}_0}\right)\right] $

      其中,det表示矩阵的行列式,

      $ \begin{align} G(L,L{{{R}}_0}) =& {{{π}}^{N(N - 1)/2}}{\det{^{\!- L}}}(L{{{R}}_0})\\ {\rm{}}& \cdot \mathop \prod \limits_{n = 1}^N \Gamma (L - n + 1) \end{align} $

      在实际应用中,矩阵${{{R}}_0}$可以通过地理信息系统、数字地形高程模型和综合脉冲孔径雷达SIAR成像等手段获取[33,34]。自由度$L$决定了协方差矩阵结构${{R}}$${{{R}}_0}$的距离,$L$越大表明先验信息越有价值。

    • 这一部分将推导出MAP-GLRT(Maximum A Posteriori GLRT)检测器,接收信号在假设${\rm{H}_1}$${\rm{H}_0}$下的概率密度函数分别表示为

      $ f({{X}}|{{R}},{\rm{H}_0}) = \frac{1}{{{{{π}}^{NK}}{{\det }^K}({{R}})}}\exp \left[ - \rm{tr} ({{{R}}^{ - 1}}{{{T}}_0})\right]\quad\quad $

      $ f({{X}}|\alpha ,{{R}},{\rm{H}_1}) =\frac{1}{{{{{π}}^{NK}}{{\det }^K}({{R}})}} \exp \left[ - \rm{tr} ({{{R}}^{ - 1}}{{{T}}_1})\right] $

      (16)

      其中,${{{T}}_0}$${{{T}}_1}$分别定义为

      $ {{{T}}_0} = {{X}}{{{X}}^{\rm{H}}} $

      $ {{{T}}_1} = \left({{X}} - \alpha {{{a}}_{\rm{r}}}{{a}}_{\rm{t}}^{\rm{H}}{{S}}\right){\left({{X}} - \alpha {{{a}}_{\rm{r}}}{{a}}_{\rm{t}}^{\rm{H}}{{S}}\right)^{\rm{H}}} $

      实际上,逆复Wishart分布是式(15)和式(16)的共轭先验分布。也就是说,如果先验分布是逆复Wishart分布,则式(15)和式(16)是后验概率密度函数(Probability Density Function, PDF),这是易处理性的基础。

      根据奈曼皮尔逊(Neyman-Pearson, NP)准则,最优的检测器可以通过似然比检验获得。然而,由于存在未知参数,使得无法采用似然比检验推导出最优检测器。因此,这里采用广义似然比检验设计检测器。广义似然比检验可以描述问题为

      $ \varLambda = \mathop {\max }\limits_\alpha \frac{{{{\max }_{{{R}}}}[f({{X}}|\alpha ,{{R}},{\rm{H}_1}){f_{{{R}}}}({{R}})]}}{{{{\max }_{{{R}}}}[f({{X}}|{{R}},{\rm{H}_0}){f_{{{R}}}}({{R}})]}}\mathop {\mathop \gtrless \limits^{{\rm{H}_1}} }\limits_{{\rm{H}_0}} \eta $

      其中,$\eta $表示检测门限。首先,需要分别计算${{R}}$${\rm{H}_0}$${\rm{H}_1}$假设下的最大后验概率(Maximum A Posteriori, MAP)估计,分别记为${\hat {{R}}_0}$${\hat {{R}}_1}$,即

      $ {\hat {{{R}}}_0} = \mathop {\arg \max }\limits_{{{R}}} \left[f({{X}}|{{R}},{\rm{H}_0}){f_{{{R}}}}({{R}})\right] $

      ${\hat {{{R}}}_1} = \mathop {\arg \max }\limits_{{{R}}} [f({{X}}|\alpha ,{{R}},{\rm{H}_1}){f_{{{R}}}}({{R}})]$

      将式(15)和式(13)代入到式(20)中,并取对数得

      $ \begin{align} {\hat {{R}}_0} = &\mathop {\arg \max }\limits_{{{R}}} \left\{ \left(K + L + N\right)\ln \left[\det \left({{{R}}^{ - 1}}\right)\right] \right.\\ {\rm{}}& \left.- {\rm{tr}} \left[({{{T}}_0} + L{{{R}}_0}){{{R}}^{ - 1}}\right]\right\} \end{align} $

      利用式(23),式(24)代数不等式

      $ {\frac{1}{n}}\sum\limits_{i = 1}^n {{x_i} \ge {{\left(\prod\limits_{i = 1}^n {{x_i}} \right)}^{{\displaystyle\scriptsize\frac{1}{n}}}}} $

      $ a\ln (x) - N{x^{{1/N}}} \le a(N\ln (a) - N) $

      从而得到

      $ \begin{aligned} {\rm{}} & (K + L + N)\ln \bigr[\det ({{{R}}^{ - 1}})\bigr] - {\rm{tr}} \bigr[({{{T}}_0} + L{{{R}}_0}){{{R}}^{ - 1}}\bigr] \\ {\rm{}} & \quad \le (K + L + N)\ln \left\{ {\det [({{{T}}_0} + L{{{R}}_0}){{{R}}^{ - 1}}]} \right\} \\ {\rm{}}&\quad \quad- (K + L + N)\ln \bigr[\det ({{{T}}_0} + L{{{R}}_0})\bigr] \\ {\rm{}} & \quad\quad - N{\det {\rm{}}^{\!\!1/N}}\bigr[({{{T}}_0} + L{{{R}}_0}){{{R}}^{ - 1}}\bigr] \\ {\rm{}} & \quad \le (K + L + N)\bigr\{ N\ln (K + L + N) - N \\ {\rm{}}& \quad\quad - \ln \bigr[\det ({{{T}}_0} + L{{{R}}_0})\bigr] \bigr\} \end{aligned} $

      进一步地,等式在式(26)情况下成立

      ${{R}} = \frac{1}{{K + L + N}}({{{T}}_0} + L{{{R}}_0})$

      ${{R}}$${\rm{H}_0}$假设下的MAP估计值为

      ${\hat {{R}}_0} = \frac{1}{{K + L + N}}({{{T}}_0} + L{{{R}}_0})$

      采用与上述类似的计算方法,可以得到${\rm{H}_1}$假设下${{R}}$的MAP估计值为

      $ {{\hat{ R}}_1} = \frac{1}{{K + L + N}}({{{T}}_1} + L{{{R}}_0}) $

      将式(27)和式(28)代入式(19)中并化简,则GLRT检测器可表示为

      $ {\varLambda _1} = \mathop {\max }\limits_\alpha \frac{{\det ({{{T}}_0} + L{{{R}}_0})}}{{\det ({{{T}}_1} + L{{{R}}_0})}}\mathop {\mathop \gtrless \limits^{{\rm{H}_1}} }\limits_{{\rm{H}_0}} {\eta _1} $

      注意到式(29)的分子与$\alpha $无关,因此式(29)检验不等式可以表示为

      $ {\varLambda _1} = \frac{{\det ({{{T}}_0} + L{{{R}}_0})}}{{\mathop {\min }\limits_\alpha \det ({{{T}}_1} + L{{{R}}_0})}}\mathop {\mathop \gtrless \limits^{{\rm{H}_1}} }\limits_{{\rm{H}_0}} {\eta _1} $

      显然,

      $ \begin{aligned} {\rm{}} & {{{T}}_1} + L{{{R}}_0} \\ {\rm{}}&\quad = {{X}}{{{X}}^{\rm{H}}} - {\alpha ^*}{{X}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}{{a}}_{\rm{r}}^{\rm{H}} - \alpha {{{a}}_{\rm{r}}}{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{X}}^{\rm{H}}} \\ {\rm{}} &\quad\quad+ |\alpha {|^2}{{{a}}_{\rm{r}}}{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}{{a}}_{\rm{r}}^{\rm{H}} + L{{{R}}_0} \\ {\rm{}} &\quad= {{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}} + L{{{R}}_0} + \left({{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}\right)\\ {\rm{}}& \quad\quad\cdot\!\left( {\alpha {{{a}}_{\rm{r}}} - \frac{{{{X}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}}{{{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}}} \right)\!{\left( {\alpha {{{a}}_{\rm{r}}} - \frac{{{{X}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}}{{{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}}} \right)^{\rm{H}}} \quad\quad \end{aligned} $

      (31)

      其中,

      $ {{P}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} = {{I}} - {{{P}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}} $

      ${{{P}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}}$是向${{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}$的列空间进行投影的正交投影矩阵,即

      $ {{{P}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}} = \frac{{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}{{a}}_{\rm{t}}^{\rm{H}}{{S}}}}{{{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}} $

      使用等式$\det ({{I}} + {{AB}}) = \det ({{I}} + {{BA}})$,则

      $ \begin{aligned} {\rm{}}&\det ({{{T}}_1} + L{{{R}}_0}) \\ {\rm{}}&\quad = \det \left( {{{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}} + L{{{R}}_0}} \right)\\ {\rm{}}& \quad\quad \cdot\det \left[ {{{I}} + \Bigr({{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}\Bigr){{\left( {{{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}} + L{{{R}}_0}} \right)}^{ - 1}}} \right. \\ {\rm{}}&\quad\quad \left. \cdot \left( {\alpha {{{a}}_{\rm{r}}} - \frac{{{{X}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}}{{{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}}} \right){{\left( {\alpha {{{a}}_{\rm{r}}} - \frac{{{{X}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}}{{{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}}} \right)}^{\rm{H}}} \right] \\ {\rm{}}&\quad = \det \left( {{{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}} + L{{{R}}_0}} \right)(1 + \gamma ) \end{aligned} $

      其中,

      $ \begin{align} \gamma =& \Bigr({{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}\Bigr){\left( {\alpha {{{a}}_{\rm{r}}} - \frac{{{{X}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}}{{{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}}} \right)^{\rm{H}}}\\ {\rm{}}& \cdot{\left( {{{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}} + L{{{R}}_0}} \right)^{ - 1}}\\ {\rm{}}& \cdot\left( {\alpha {{{a}}_{\rm{r}}} - \frac{{{{X}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}}{{{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}}} \right) \end{align} $

      现在需要计算$1 + \gamma $的值,其对$\alpha $求导数并令结果为0,有

      $ \hat \alpha = \frac{{{{a}}_{\rm{r}}^{\rm{H}}{{\left( {{{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}} + L{{{R}}_0}} \right)}^{ - 1}}{{X}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}}{{\Bigr({{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}\Bigr){{a}}_{\rm{r}}^{\rm{H}}{{\left( {{{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}} + L{{{R}}_0}} \right)}^{ - 1}}{{{a}}_{\rm{r}}}}} $

      因此,可以得到

      $ \begin{aligned} {\rm{}}& \mathop {\min }\limits_\alpha (1 + \gamma ) \\ {\rm{}}& \quad = 1 + \frac{{{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{X}}^{\rm{H}}}}}{{{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}}{\left( {{{XP}}_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^ \bot {{{X}}^{\rm{H}}} + L{{{R}}_0}} \right)^{ - 1}}\\ {\rm{}}& \quad\quad\cdot({{I}} - {{Z}}){{X}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}} \end{aligned} $

      其中,

      $ {{Z}} = \frac{{{{{a}}_{\rm{r}}}{{a}}_{\rm{r}}^{\rm{H}}{{\left( {{{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}} + L{{{R}}_0}} \right)}^{ - 1}}}}{{{{a}}_{\rm{r}}^{\rm{H}}{{\left( {{{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}} + L{{{R}}_0}} \right)}^{ - 1}}{{{a}}_{\rm{r}}}}} $

      将式(37)代入式(34),得

      $ \begin{aligned} {\rm{}}& \mathop {\min }\limits_\alpha \det ({{{T}}_1} + L{{{R}}_0}) \\ {\rm{}}& \quad = \det \left( {{{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}} + L{{{R}}_0}} \right)\mathop {\min }\limits_\alpha (1 + \gamma ) \\ {\rm{}}& \quad= \det \left( {{{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}} + L{{{R}}_0}} \right)\\ {\rm{}}& \quad\quad\cdot\det \Biggr[ {{{I}} + {{\left( {{{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}} + L{{{R}}_0}} \right)}^{ - 1}}} \\ {\rm{}}& \quad \quad\cdot {\frac{{({{I}} - {{Z}}){{X}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{X}}^{\rm{H}}}}}{{{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}}} \Biggr] \\ {\rm{}}& \quad= \det \Biggr[\!\! {{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}} \!+\! L{{{R}}_0} \\ {\rm{}}& \quad\quad+\! \frac{{({{I}} \!-\! {{Z}}){{X}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}{{a}}_{\rm{t}}^{\rm{H}}{{S}}\!{{{X}}^{\rm{H}}}}}{{{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}} \Biggr] \\ {\rm{}}&\quad= \det \left( {{{X}}\!{{{X}}^{\rm{H}}} + L{{{R}}_0} - \frac{{{{ZX}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}{{a}}_{\rm{t}}^{\rm{H}}{{S}}\!{{{X}}^{\rm{H}}}}}{{{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}}} \right) \\ {\rm{}}&\quad = \det \left({{X}}\!{{{X}}^{\rm{H}}} + L{{{R}}_0}\right)\det ({{I}} - {{U}}) \end{aligned} $

      其中,

      $ {{U}} = \frac{{{{\Bigr({{X}}{{{X}}^{\rm{H}}} + L{{{R}}_0}\Bigr)}^{ - 1}}{{ZX}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{X}}^{\rm{H}}}}}{{{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}} $

      现在,需要计算$\det ({{I}} - {{U}})$的值,注意到它可以写为

      $ \det ({{I}} \!-\! {{U}}) = 1 \!-\! \frac{\omega }{{{{a}}_{\rm{r}}^{\rm{H}}{{\left( {{{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}} \!+\! L{{{R}}_0}} \right)}^{ - 1}}{{{a}}_{\rm{r}}}}} $

      其中,

      $ \begin{align} \omega =& {{a}}_{\rm{r}}^{\rm{H}}{\left( {{{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}} + L{{{R}}_0}} \right)^{ - 1}}\\ {\rm{}}& \cdot\left( {{{X}}{{{X}}^{\rm{H}}} - {{XP}}_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot }{{{X}}^{\rm{H}}}} \right){\Bigr( {{{X}}{{{X}}^{\rm{H}}} + L{{{R}}_0}} \Bigr)^{ - 1}}{{{a}}_{\rm{r}}} \end{align} $

      (42)

      注意到

      $ \begin{align} {{X}}{{{X}}^{\rm{H}}} - {{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}} = &\Bigr({{X}}{{{X}}^{\rm{H}}} + L{{{R}}_0}\Bigr) \\ {\rm{}}& - \Bigr({{XP}}_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^ \bot {{{X}}^{\rm{H}}} + L{{{R}}_0}\Bigr) \end{align} $

      将式(43)和式(42)代入式(41),得到

      $ \det ({{I}} - {{U}}) = \frac{{{{a}}_{\rm{r}}^{\rm{H}}{{\Bigr({{X}}{{{X}}^{\rm{H}}} + L{{{R}}_0}\Bigr)}^{ - 1}}{{{a}}_{\rm{r}}}}}{{{{a}}_{\rm{r}}^{\rm{H}}{{\left( {{{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}} + L{{{R}}_0}} \right)}^{ - 1}}{{{a}}_{\rm{r}}}}} $

      将式(44)代入式(39),可以得到

      $ \begin{align} {\rm{}}&\mathop {\min }\limits_\alpha \det ({{{T}}_1} + L{{{R}}_0}) \\ {\rm{}}&\quad= \det \Bigr({{X}}{{{X}}^{\rm{H}}} + L{{{R}}_0}\Bigr)\\ {\rm{}}& \quad \quad\cdot\frac{{{{a}}_{\rm{r}}^{\rm{H}}{{\Bigr({{X}}{{{X}}^{\rm{H}}} + L{{{R}}_0}\Bigr)}^{ - 1}}{{{a}}_{\rm{r}}}}}{{{{a}}_{\rm{r}}^{\rm{H}}{{\left( {{{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}} + L{{{R}}_0}} \right)}^{ - 1}}{{{a}}_{\rm{r}}}}} \end{align} $

      将式(45)代入式(30),最终得到GLRT检测器,记为MAP-GLRT:

      $ {\varLambda _2} = \frac{{{{a}}_{\rm{r}}^{\rm{H}}{{\left( {{{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}} + L{{{R}}_0}} \right)}^{ - 1}}{{{a}}_{\rm{r}}}}}{{{{a}}_{\rm{r}}^{\rm{H}}{{({{X}}{{{X}}^{\rm{H}}} + L{{{R}}_0})}^{ - 1}}{{{a}}_{\rm{r}}}}}\mathop {\mathop \gtrless \limits^{{\rm{H}_1}} }\limits_{{\rm{H}_0}} {\eta _2} $

      其中,${\eta _2}$表示式(19)中门限的修正值。

    • 这一部分采用两步法来设计基于贝叶斯方法的BAMF(Bayesian AMF)检测器。首先,假设协方差矩阵${{R}}$是已知的,则GLRT可以表示为

      $ \varDelta = \frac{{\mathop {\max }\limits_\alpha f({{X}}|\alpha ,{{R}},{\rm{H}_1})}}{{f({{X}}|{{R}},{\rm{H}_0})}}\mathop {\mathop \gtrless \limits^{{\rm{H}_1}} }\limits_{{\rm{H}_0}} \delta $

      其中,$\delta $表示检测门限。

      式(16)取对数,并对$\alpha $求导得到

      $ {{\partial \ln f({{X}}|\alpha ,{{R}},{\rm{H}_1})}/{\partial \alpha }} \!=\! {{a}}_{\rm{t}}^{\rm{H}}{{S}}\!\!\left({{X}} - \alpha {{{a}}_{\rm{r}}}{{a}}_{\rm{t}}^{\rm{H}}{{S}}\right){{{R}}^{ - 1}}{{{a}}_{\rm{r}}} $

      (48)

      令式(48)为$0$,则得到${\rm{H}_1}$假设下$\alpha $的最大似然估计为

      ${\hat \alpha _1} = \frac{{{{a}}_{\rm{r}}^{\rm{H}}{{{R}}^{ - 1}}{{X}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}}{{{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}{{a}}_{\rm{r}}^{\rm{H}}{{{R}}^{ - 1}}{{{a}}_{\rm{r}}}}}$

      将式(49)代入式(47),经过化简,得到${{R}}$已知情形下的GLRT检测器为

      $ {\varDelta _1} = \frac{{\Bigr|{{a}}_{\rm{r}}^{\rm{H}}{{{R}}^{ - 1}}{{X}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}{\Bigr|^2}}}{{{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}{{a}}_{\rm{r}}^{\rm{H}}{{{R}}^{ - 1}}{{{a}}_{\rm{r}}}}}\mathop {\mathop \gtrless \limits^{{\rm{H}_1}} }\limits_{{\rm{H}_0}} {\delta _1} $

      为了得到一个完全自适应的AMF检测器,需要求出${\rm{H}_1}$假设下${{R}}$的MAP估计,即

      $ {\hat {{R}}_1} = \mathop {\arg \max }\limits_{{{R}}} [f({{X}}|{\hat \alpha _1},{{R}},{\rm{H}_1}){f_{{{R}}}}({{R}})] $

      其中,$\arg {\max _{{{R}}}}[ \cdot ]$表示使得表达式最大化的${{R}}$的值。将式(13),式(16)和式(49)代入式(51)并取对数,得到

      $ \begin{aligned} {\rm{}}& \ln [f({{X}}|{{\hat \alpha }_1},{{R}},{\rm{H}_1}){f_{{{R}}}}({{R}})] \\ {\rm{}}&\quad = - NK\ln {{π}} - (K + L + N)\ln [\det ({{R}})] \\ {\rm{}}&\quad\quad- \ln [G(L,L{{{R}}_0})] - {\rm{tr}} \Bigr({{{X}}^{\rm{H}}}{{{R}}^{ - 1}}{{X}} + L{{{R}}^{ - 1}}{{{R}}_0}\Bigr) \\ {\rm{}}&\quad\quad - \frac{{\Bigr|{{a}}_{\rm{r}}^{\rm{H}}{{{R}}^{ - 1}}{{X}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}{\Bigr|^2}}}{{{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}{{a}}_{\rm{r}}^{\rm{H}}{{{R}}^{ - 1}}{{{a}}_{\rm{r}}}}} \end{aligned} $

      式(52)对${{R}}$求导并令结果等于0,经过化简,得到

      $ \begin{aligned} {\rm{}} & (K + L + N){{R}} \\ {\rm{}} &\quad= {{X}}{{{X}}^{\rm{H}}} + L{{{R}}_0} + \frac{{|{{a}}_{\rm{r}}^{\rm{H}}{{{R}}^{ - 1}}{{X}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}{|^2}{{{a}}_{\rm{r}}}{{a}}_{\rm{r}}^{\rm{H}}}}{{{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}{{({{a}}_{\rm{r}}^{\rm{H}}{{{R}}^{ - 1}}{{{a}}_{\rm{r}}})}^2}}} \\ {\rm{}} & \quad\quad -\left[ {{{a}}_{\rm{r}}}{{a}}_{\rm{r}}^{\rm{H}}{{{R}}^{ - 1}}\!{{X}}{{{S}}^{\rm{H}}}\!{{{a}}_{\rm{t}}}{{a}}_{\rm{t}}^{\rm{H}}\!{{S}}{{{X}}^{\rm{H}}}\right. \\ {\rm{}}& \quad\quad\left.+ {{X}} {{{S}}^{\rm{H}}}\!{{{a}}_{\rm{t}}}{{a}}_{\rm{t}}^{\rm{H}}\!{{S}}{{{X}}^{\rm{H}}}\!{{{R}}^{ - 1}}\!{{{a}}_{\rm{r}}}{{a}}_{\rm{r}}^{\rm{H}}\right]\\ {\rm{}}&\quad\quad /\left[{{{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}} {{{a}}_{\rm{t}}}\!{{a}}_{\rm{r}}^{\rm{H}}\!{{{R}}^{ - 1}}\!{{{a}}_{\rm{r}}}}} \right] \end{aligned} $

      式(53)等号两边都右乘${{{R}}^{ - 1}}{{{a}}_{\rm{r}}}$

      $ \begin{align} (K + L + N){{{a}}_{\rm{r}}} =& {{X}}{{{X}}^{\rm{H}}}{{{R}}^{ - 1}}{{{a}}_{\rm{r}}} + L{{{R}}_0}{{{R}}^{ - 1}}{{{a}}_{\rm{r}}} \\ {\rm{}} & - \frac{{{{X}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{X}}^{\rm{H}}}{{{R}}^{ - 1}}{{{a}}_{\rm{r}}}}}{{{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}} \\ {\rm{}} =& \left( {{{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}} + L{{{R}}_0}} \right){{{R}}^{ - 1}}{{{a}}_{\rm{r}}} \quad\quad \end{align} $

      ${{R}}$${\rm{H}_1}$假设下的MAP估计为

      $ {\hat {{{R}}}_1} = \frac{{{{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}} + L{{{R}}_0}}}{{K + L + N}} $

      将式(55)代入式(50)并化简,最终得到BAMF检测器为

      $ {\varDelta _2} = \frac{{{{\left| {{{a}}_{\rm{r}}^{\rm{H}}{{\left( {{{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}} + L{{{R}}_0}} \right)}^{ - 1}}{{X}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}} \right|}^2}}}{{{{a}}_{\rm{r}}^{\rm{H}}{{\left( {{{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}} + L{{{R}}_0}} \right)}^{ - 1}}{{{a}}_{\rm{r}}}{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}}\mathop {\mathop \gtrless \limits^{{\rm{H}_1}} }\limits_{{\rm{H}_0}} {\delta _2} $

      (56)

      其中,${\delta _2}$是式(47)中门限的修正值。

    • 这一节,将提出的检测器与传统的非贝叶斯检测器进行检测性能比较。因为MAP-GLRT检测器与BAMF检测器的闭集形式的虚警概率(Probability of False Alarm, ${{P}_{\rm{fa}}}$) 和检测概率(Probability of Detection, ${{P}_{\rm{d}}}$)表达式无法得到,所以采用蒙特卡罗方法评估检测性能。检测门限和检测概率分别通过$100/{{P} _{\rm fa}}$${10^4}$次独立实验得到。在以下仿真中,发射和接收阵元数设置为$M = N = 8$,阵元间距为半波长。假设${{R}} \sim {\mathcal{W}^{\ - 1}}(L,L{{{R}}_0})$,其中$L = 10$, ${{{R}}_0}$采用指数相关的协方差矩阵,即第$(m,n)$个元素是

      $ {{{R}}_0}(m,n) = \sigma _{\rm c}^2{{\rho }^{|m - n|}} $

      这里$\rm{\rho } = 0.9$表示1阶迟滞相关系数,$\sigma _{\rm c}^2$为杂波功率,信杂比(Signal-to-Clutter Ratio, SCR)定义为

      $\rm{SCR} = |\alpha {|^2}{{a}}_{\rm{r}}^{\rm{H}}{{R}}_0^{ - 1}{{{a}}_{\rm{r}}}{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}$

      为了对比,这里给出传统的非贝叶斯GLRT检测器[20]

      $ {t_{\rm GLRT}} = 1 - \frac{{{{a}}_{\rm{r}}^{\rm{H}}{{\Bigr({{X}}{{{X}}^{\rm{H}}}\Bigr)}^{ - 1}}{{{a}}_{\rm{r}}}}}{{{{a}}_{\rm{r}}^{\rm{H}}{{\left( {{{XP}}_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^ \bot {{{X}}^{\rm{H}}}} \right)}^{ - 1}}{{{a}}_{\rm{r}}}}} $

      Rao检测器[26]

      $ {t_{\rm Rao}} = \frac{{\Bigr|{{a}}_{\rm{r}}^{\rm{H}}{{\Bigr({{X}}{{{X}}^{\rm{H}}}\Bigr)}^{ - 1}}{{X}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}{\Bigr|^2}}}{{{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}{{a}}_{\rm{r}}^{\rm{H}}{{\Bigr({{X}}{{{X}}^{\rm{H}}}\Bigr)}^{ - 1}}{{{a}}_{\rm{r}}}}} $

      Wald检测器[26]

      $ {t_{\rm Wald}} = \frac{{{{\left| {{{a}}_{\rm{r}}^{\rm{H}}{{\Bigr({{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}}\Bigr)}^{ - 1}}{{X}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}} \right|}^2}}}{{{{a}}_{\rm{t}}^{\rm{H}}{{S}}{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}{{a}}_{\rm{r}}^{\rm{H}}{{\Bigr({{XP}}\!\!_{{{{S}}^{\rm{H}}}{{{a}}_{\rm{t}}}}^{\ \ \bot} {{{X}}^{\rm{H}}}\Bigr)}^{ - 1}}{{{a}}_{\rm{r}}}}} $

      图1给出了K取不同值时的接收器特性(Receiver Operating Characteristic, ROC)曲线,其中,$L \!=\! 10$, $\rm{SCR} \!=\! 10\;dB$。仿真参数设置为(a): $K \!=\! 1.5N$; (b): $K = 2N$; (c): $K = 6N$。本文提出的贝叶斯检测器的性能明显优于传统的非贝叶斯检测器,尤其是在样本数$K$较小时,这是因为贝叶斯检测器使用了杂波协方差矩阵的先验知识。随着$K$的增大,所有检测器的性能都有所提高,与此同时,本文的贝叶斯检测器相对其他检测器的性能增益在逐渐降低。当发射波形采样数足够大(例如,$K = 6N$) 时,几种检测器的性能是相当的。这是预料之中的,因为在此情形下即使不采用先验信息估计的协方差矩阵也足够可靠了。

      图  1  K 取不同值时的ROC曲线

      Figure 1.  ROC curves for different K

      图2给出了${{{P}}_{\rm{d}}}$关于SCR的曲线图,其中$L = 10$, ${{{P}}_{{\rm{fa}}}} = 1{0^{ - 3}}$。仿真参数设置为(a): $K = 1.5N$; (b): $K = 2N$; (c): $K = 6N$。可以看到,与图1类似,当发射波形样本数较小时本文的检测器具有较大的性能优势,当样本数足够大时$5$种检测器的性能几乎是相同的。上述仿真结果说明,当样本数较小时利用协方差矩阵${{R}}$的先验知识是十分有效的策略。

      图  2  K 取不同值时${{{P}}_{\rm{d}}}$关于SCR的曲线

      Figure 2.  ${{{P}}_{\rm{d}}}$ versus SCR for different K

      图3显示了自由度$L$取不同值时MAP-GLRT检测器与BAMF检测器的ROC曲线,仿真参数设置为$K = 1.5\;N$$\rm{SCR} = 10\;dB$。可以看到,检测性能随着$L$的增大而增强,这是因为$L$越大,协方差矩阵的离差越小,表明获得的先验信息越有价值。

      图  3  L 取不同值时检测器的ROC曲线

      Figure 3.  ROC curves for different L

      注意到上述仿真均假设逆复 Wishart 分布的参数 (即$L$${{{R}}_0}$) 是完全已知的。实际上,这些参数是未知的,可以使用文献[32]中的方法进行估计,这就不可避免地产生估计误差,接下来研究贝叶斯检测器在参数$L$${{{R}}_0}$失配情况下的检测性能。首先,假设$L$是确定已知的而${{{R}}_0}$存在失配,图4仿真了这种情况下检测器的ROC曲线。这里采用两个$1$阶迟滞相关系数${\rho _{\rm a}}$${\rho _{\rm n}}$去仿真这种失配,其中${\rho _{\rm a}}$${\rho _{\rm n}}$分别用在实际的${{{R}}_0}$和名义上的${{{R}}_0}$中。仿真参数设置为:真实的$1$阶迟滞相关系数${\rho _{\rm a}} = 0.90$, 设定的$1$阶迟滞相关系数分别为${\rho _{\rm n}} = 0.60$, ${\rho _{\rm n}} = 0.90$, ${\rho _{\rm n}} = 0.98$,其他参数设置为$M = N = 8$, $K = 1.5N$$L = 10$图4显示$1$阶迟滞相关系数的失配导致了检测器性能的下降。

      图  4  1阶迟滞相关系数失配时检测器的ROC曲线

      Figure 4.  ROC curves when one-lag correlation coefficient is mismatched

      现在,假设${{{R}}_0}$是精确估计的而参数$L$存在失配。图5仿真了这种情形下检测器的ROC曲线。仿真参数设置为:真实自由度${L_{\rm a}} = 10$,名义上的自由度分别为${L_{\rm n}} = 2$, ${L_{\rm n}} = 10$${L_{\rm n}} = 100$,其他参数设置为$M = N = 8$, $K = 1.5N$, $\rm{SCR} = 10\;dB$$\rho = 0.9$。仿真结果表明自由度的失配会导致检测性能的下降。图4图5的仿真结果表明本文的贝叶测器对参数$L$${{{R}}_0}$比较敏感,所以实际中需要尽可能精确地估计这两个参数。

      图  5  自由度L失配时检测器的ROC曲线

      Figure 5.  ROC curves when L is mismatched

    • 本文研究了协方差矩阵未知的高斯杂波背景下 MIMO 雷达的自适应检测问题,杂波的协方差矩阵被建模为服从逆复 Wishart 分布的随机矩阵,然后基于贝叶斯方法和广义似然比检验准则设计了MAP-GLRT 检测器与BAMF检测器,这两个检测器不需要训练数据。研究结果表明利用协方差矩阵的先验信息可以有效提高检测器的检测性能,尤其是在发射波形采样数较少时。在1阶迟滞相关系数或自由度失配时,贝叶斯检测器的检测性能会有所下降。

参考文献 (34)

目录

    /

    返回文章
    返回