基于跟踪脉冲的MIMO雷达多目标微动特征提取

罗迎 龚逸帅 陈怡君 张群

引用本文:
Citation:

基于跟踪脉冲的MIMO雷达多目标微动特征提取

    作者简介: 陈怡君(1989–),女,陕西西安人,现为武警工程大学信息工程学院讲师,研究方向为认知成像与目标识别。E-mail: chenyijun519@126.com.
    通讯作者: 罗迎  luoying2002521@163.com
  • 基金项目:

    国家自然科学基金(61571457,61631019)

Multi-target Micro-motion Feature Extraction Based on Tracking Pulses in MIMO Radar

    Corresponding author:
  • Fund Project: The National Natural Science Foundation of China (61571457, 61631019)

  • 摘要: 微动特征是空间目标识别的重要特征信息之一。然而,现有的多功能多输入多输出(Multi-Input Multi-Output, MIMO)雷达通常需要在完成目标搜索和跟踪任务之后为目标微动特征提取分配大量连续的时间资源,导致目标识别实时性能和雷达系统整体工作性能均不高。针对该问题,该文提出了一种基于跟踪脉冲的MIMO雷达多目标微动特征提取方法。首先依据各目标的方位信息对MIMO雷达发射波形进行设计,为不同方向目标同时发射跟踪脉冲;在此基础上,综合考虑目标微动特征提取性能以及目标跟踪性能的需求,对跟踪脉冲的发射时间序列进行优化设计;最后,直接利用窄带跟踪脉冲实现对不同方向目标微动特征的同时提取,无需再为目标微动特征提取分配额外的时间资源,有效提升目标识别实时性和雷达工作效率。仿真实验表明,在信噪比大于–10 dB时,所提方法能够实现多目标微动特征的准确提取,具有良好的有效性和鲁棒性。
  • 图 1  MIMO雷达发射阵列示意图

    Figure 1.  Transmit array of MIMO radar

    图 2  旋转目标的几何模型

    Figure 2.  Geometry of a rotating target

    图 3  算法流程

    Figure 3.  Flow chart of the algorithm

    图 4  发射方向图

    Figure 4.  The transmit beampattern

    图 5  无噪声条件下文献[18]方法重构的时频分布

    Figure 5.  The reconstructed time-frequency distribution using method in Ref. [18] without noise

    图 7  0°方向所提方法重构的时频分布

    Figure 7.  The reconstructed time-frequency distribution in 0° direction using the proposed method

    图 6  –30°方向所提方法重构的时频分布

    Figure 6.  The reconstructed time-frequency distribution in –30° direction using the proposed method

    图 8  算法性能随信噪比变化曲线

    Figure 8.  The performance curve versus SNR

    表 1  目标参数

    Table 1.  The parameters of the targets

    目标方位(°) 参考点坐标(km) 旋转角速度 $\left( {\rm{rad}} \cdot {{\rm{s}}^{{\rm{ - 1}}}}\right) $ 目标散射点坐标(m)
    –30 (15, 26, 0) ${\rm{(}}{π} {\rm{,}}{π} {\rm{/2,}}{π} {\rm{)}}$ $(0,0,0),({\rm{1}},{\rm{1}},{\rm{0}}),( - {\rm{1}}, - {\rm{1}},{\rm{0}})$
    0 (0, 40, 0) ${\rm{(0,0,}}{π} {\rm{)}}$ $(2,0,0),( - 2,0,0),(0, - 2,0),(0,2,0)$
    下载: 导出CSV

    表 2  文献[18]方法目标微动特征提取结果

    Table 2.  The micro-motion features of targets obtained by method in Ref. [18]

    目标方向(°) 特征参数序号 目标散射点坐标 ${{r}_0}$ (m) 旋转角速度 ${{w}}$ $\left( {\rm{rad}} \cdot {{\rm{s}}^{{\rm{ - 1}}}}\right) $ 系数A
    –30 1 $(0.1 \ 0 \ - 0.1)$ $(3.14 \ 1.57 \ 3.14)$ 30.40
    2 $({\rm{1 \ 1 \ 0}})$ $(3.14 \ 1.57 \ 3.14)$ 27.07
    3 $( - {\rm{1.0 }}\ - {\rm{1.1}} \ {\rm{ 0}})$ $(3.14 \ 1.57 \ 3.14)$ 25.11
    4 $({\rm{0}}{\rm{.3 }}\ 0\ {\rm{ 0}})$ $(3.14 \ 1.57 \ 3.14)$ 5.96
    5 $(1.5 \ - {\rm{1.2}}\ {\rm{0}})$ $(2.83 \ 1.57 \ 3.14)$ 4.29
    0 1 $( - 1.8 \ 0.1 \ 0)$ $(0.31 \ 0 \ 2.51)$ 29.15
    2 $(0.1 \ - 1.9 \ 0)$ $(0.31 \ 0 \ 2.51)$ 28.10
    3 $(0.5 \ 2.4 \ 0)$ $(0.31 \ 0 \ 2.51)$ 24.27
    4 $(2.6 \ - 0.2 \ 0)$ $(0.31 \ 0 \ 2.51)$ 23.31
    5 $( - 1.9 \ - {\rm{0.4 }}\ {\rm{0}})$ $(0 \ 0.31 \ 2.83)$ 3.15
    下载: 导出CSV

    表 3  无噪声条件下所提方法目标微动特征提取结果

    Table 3.  The micro-motion features of targets obtained by the proposed method without noise

    目标方向(°) 特征参数序号 目标散射点坐标 ${{r}_0}$ (m) 旋转角速度 ${{w}}$ ( ${\rm{rad}} \cdot {{\rm{s}}^{{\rm{ - 1}}}}$ ) 系数A
    –30 1 $(0.1 \ 0 \ - 0.1)$ $(3.14 \ 1.57 \ 3.14)$ 27.23
    2 $({\rm{1 \ 1 \ 0}})$ $(3.14 \ 1.57 \ 3.14)$ 26.47
    3 $( - {\rm{1.0}} \ - {\rm{1.1}} \ {\rm{ 0}})$ $(3.14 \ 1.57 \ 3.14)$ 26.01
    4 $( - 1.3 \ - {\rm{1.2}}\ {\rm{0}})$ $(2.83 \ 1.57 \ 3.14)$ 6.40
    5 $({\rm{0}}{\rm{.3 }} \ 0 \ {\rm{ 0}})$ $(3.14 \ 1.88 \ 3.14)$ 6.06
    0 1 $(2.0 \ 0.1 \ 0)$ $(0 \ 0 \ 3.14)$ 28.63
    2 $( - 1.9 \ 0 \ 0)$ $(0 \ 0 \ 3.14)$ 28.39
    3 $(0.1 \ - 2.0 \ 0)$ $(0 \ 0 \ 3.14)$ 27.71
    4 $(0 \ 1.9 \ 0)$ $(0 \ 0 \ 3.14)$ 27.19
    5 $(1.8 \ 0.1 \ 0.3)$ $(0 \ 0.31 \ 2.83)$ 6.72
    下载: 导出CSV

    表 4  SNR=–10 dB条件下所提方法目标微动特征提取结果

    Table 4.  The micro-motion features of targets obtained by the proposed method when SNR=–10 dB

    目标方向(°) 特征参数序号 目标散射点坐标 ${{r}_0}$ (m) 旋转角速度 ${{w}}$ $\left( {\rm{rad}} \cdot {{\rm{s}}^{{\rm{ - 1}}}}\right) $ 系数A
    –30 1 $(0.1 \ 0 \ - 0.1)$ $(3.14 \ 1.57 \ 3.14)$ 26.83
    2 $({\rm{1 \ 1 \ 0}})$ $(3.14 \ 1.57 \ 3.14)$ 25.48
    3 $( - {\rm{1.0}} - {\rm{1.1}}\ {\rm{ 0}})$ $(3.14 \ 1.57 \ 3.14)$ 25.09
    4 $( - 1.3 \ - {\rm{1.2}}\ {\rm{ 0}})$ $(2.83 \ 1.57 \ 3.14)$ 7.78
    5 $({\rm{0}}{\rm{.3 }}\ 0 \ {\rm{ 0}})$ $(3.14 \ 1.88 \ 3.14)$ 7.31
    0 1 $(2.0 \ 0.1 \ 0)$ $(0 \ 0 \ 3.14)$ 27.52
    2 $( - 1.9 \ 0 \ 0)$ $(0 \ 0 \ 3.14)$ 27.28
    3 $(0.1 \ - 2.0 \ 0)$ $(0 \ 0 \ 3.14)$ 26.59
    4 $(0 \ 1.9 \ 0)$ $(0 \ 0 \ 3.14)$ 26.08
    5 $(1.8 \ 0.1 \ 0.3)$ $(0 \ 0.31 \ 2.83)$ 9.15
    下载: 导出CSV

    表 5  SNR=–15 dB条件下所提方法目标微动特征提取结果

    Table 5.  The micro-motion features of targets obtained by the proposed method when SNR=–15 dB

    目标方向(°) 特征参数序号 目标散射点坐标 ${{r}_0}$ (m) 旋转角速度 ${{w}}$ $\left( {\rm{rad}} \cdot {{\rm{s}}^{{\rm{ - 1}}}}\right) $ 系数A
    –30 1 $(0.2 \ 0.1 \ - 0.1)$ $(3.14 \ 1.57 \ 3.14)$ 25.73
    2 $({\rm{0.9}}{\rm{ \ 1.0 \ 0}})$ $(3.14 \ 1.57 \ 3.14)$ 23.27
    3 $( - {\rm{1}}{\rm{.2 }}\ - {\rm{1.0}}\ {\rm{0.1}})$ $(3.14 \ 1.57 \ 3.14)$ 22.16
    4 $( - 1.5 \ - {\rm{1.3}}\ {\rm{ 0}})$ $(2.83 \ 1.57 \ 3.14)$ 15.04
    5 $({\rm{0}}{\rm{.8 }} \ 0 \ {\rm{ 0}}{\rm{.5}})$ $(2.83 \ 2.19 \ 3.14)$ 10.58
    0 1 $(2.2 \ - 0.1 \ 0)$ $(0 \ 0 \ 3.14)$ 26.81
    2 $( - 2.0 \ 0.1 \ 0)$ $(0 \ 0 \ 3.14)$ 25.35
    3 $(0.2 \ - 1.8 \ 0)$ $(0 \ 0 \ 3.14)$ 25.19
    4 $(0.1 \ 1.8 \ 0.2)$ $(0 \ 0 \ 3.14)$ 24.07
    5 $(1.7 \ 0.3 \ 0.4)$ $(0 \ 0.31 \ 2.83)$ 13.15
    下载: 导出CSV
  • [1] Zhang Q, Luo Y, and Chen Y A. Micro-Doppler Characteristics of Radar Targets[M]. Jonathan Simpson, United Kingdom: Elsevier Science, 2016: 4–6.
    [2] 杨琪, 邓彬, 王宏强, 等. 太赫兹雷达目标微动特征提取研究进展[J]. 雷达学报, 2018, 7(1): 22–45. DOI: 10.12000/JR17087Yang Qi, Deng Bin, Wang Hong-qiang, et al. Advancements in research on micro-motion feature extraction in the terahertz region[J]. Journal of Radars, 2018, 7(1): 22–45. DOI: 10.12000/JR17087
    [3] 陈小龙, 董云龙, 李秀友, 等. 海面刚体目标微动特征建模及特性分析[J]. 雷达学报, 2015, 4(6): 630–638. DOI: 10.12000/JR15079Chen Xiao-long, Dong Yun-long, Li Xiu-you, et al. Modeling of micromotion and analysis of properties of rigid marine targets[J]. Journal of Radars, 2015, 4(6): 630–638. DOI: 10.12000/JR15079
    [4] 冯存前, 李靖卿, 贺思三, 等. 组网雷达中弹道目标微动特征提取与识别综述[J]. 雷达学报, 2015, 4(6): 609–620. DOI: 10.12000/JR15084Feng Cun-qian, Li Jing-qing, He Si-san, et al. Micro-Doppler feature extraction and recognition based on netted radar for ballistic targets[J]. Journal of Radars, 2015, 4(6): 609–620. DOI: 10.12000/JR15084
    [5] Chen V C, Li F, Ho S S, et al. Micro-Doppler effect in radar: Phenomenon, model, and simulation study[J]. IEEE Transactions on Aerospace and Electronic Systems, 2006, 42(1): 2–21.
    [6] Bai X R, Zhou F, and Bao Z. High-resolution three-dimensional imaging of space targets in micromotion[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2015, 8(7): 3428–3440. DOI: 10.1109/JSTARS.2015.2431119
    [7] Li G and Varshney P K. Micro-Doppler parameter estimation via parametric sparse representation and pruned orthogonal matching pursuit[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2014, 7(12): 4937–4948. DOI: 10.1109/JSTARS.2014.2318596
    [8] Du L, Ma Y Y, Wang B S, et al. Noise-robust classification of ground moving targets based on time-frequency feature from micro-Doppler signature[J]. IEEE Sensors Journal, 2014, 14(8): 2672–2682. DOI: 10.1109/JSEN.2014.2314219
    [9] 牛杰, 刘永祥, 秦玉亮, 等. 一种基于经验模态分解的锥体目标雷达微动特征提取新方法[J]. 电子学报, 2011, 39(7): 1712–1715Niu Jie, Liu Yong-xiang, Qin Yu-liang, et al. A new method of radar micro-motion feature extraction of cone target based on empirical mode decomposition[J]. Acta Electronica Sinica, 2011, 39(7): 1712–1715
    [10] 何其芳, 张群, 罗迎, 等. 正弦调频Fourier-Bessel变换及其在微动目标特征提取中的应用[J]. 雷达学报, 待出版. DOI: 10.12000/JR17069He Qi-fang, Zhang Qun, Luo Ying, et al. A sinusoidal frequency modulation Fourier-Bessel transform and its application to micro-Doppler feature extraction[J]. Journal of Radars, in press. DOI: 10.12000/JR17069
    [11] Li J and Stoica P. MIMO radar with colocated antennas[J]. IEEE Signal Processing Magazine, 2007, 24(5): 106–114. DOI: 10.1109/MSP.2007.904812
    [12] Ma C Z, Yeo T S, Liu Z F, et al. Target imaging based on l1l0 norms homotopy sparse signal recovery and distributed MIMO antennas[J]. IEEE Transactions on Aerospace and Electronic Systems, 2015, 51(4): 3399–3414. DOI: 10.1109/TAES.2015.140939
    [13] Li J, Stoica P, and Zheng X Y. Signal synthesis and receiver design for MIMO radar imaging[J]. IEEE Transactions on Signal Processing, 2008, 56(8): 3959–3968. DOI: 10.1109/TSP.2008.923197
    [14] 李慧, 赵永波, 程增飞. 基于线性调频时宽的MIMO雷达正交波形设计[J]. 电子与信息学报, 2018, 40(5): 1151–1158. DOI: 10.11999/JEIT170426Li Hui, Zhao Yong-bo, and Cheng Zeng-fei. MIMO radar orthogonal waveform set design based on chirp durations[J]. Journal of Electronics&Information Technology, 2018, 40(5): 1151–1158. DOI: 10.11999/JEIT170426
    [15] 李玉翔, 任修坤, 孙扬, 等. 一种循环迭代的宽带MIMO雷达正交稀疏频谱波形设计方法[J]. 电子与信息学报, 2017, 39(4): 953–959. DOI: 10.11999/JEIT160597Li Yu-xiang, Ren Xiu-kun, Sun Yang, et al. Cyclic iterative method for wideband MIMO radar orthogonal sparse frequency waveform design[J]. Journal of Electronics&Information Technology, 2017, 39(4): 953–959. DOI: 10.11999/JEIT160597
    [16] Mallat S G and Zhang Z F. Matching pursuits with time-frequency dictionaries[J]. IEEE Transactions on Signal Processing, 1993, 41(12): 3397–3415. DOI: 10.1109/78.258082
    [17] Chen Y J, Zhang Q, Ma C Z, et al. Micromotion feature extraction of radar target using tracking pulses with adaptive pulse repetition frequency adjustment[J]. Journal of Applied Remote Sensing, 2014, 8(1): 083569. DOI: 10.1117/1.JRS.8.083569
    [18] 叶淋美. 基于压缩感知的雷达信号处理应用研究[D]. [硕士论文], 厦门大学, 2014: 17–34.Ye Lin-mei. Study on application of compressive sensing in radar signal processing[D]. [Master dissertation], Xiamen University, 2014: 17–34.
  • [1] 冯存前李靖卿贺思三张豪 . 组网雷达中弹道目标微动特征提取与识别综述. 雷达学报, 2015, 4(6): 609-620. doi: 10.12000/JR15084
    [2] 许小剑刘永泽 . MIMO雷达三维干涉诊断成像方法. 雷达学报, 2018, 7(6): 655-663. doi: 10.12000/JR18088
    [3] 巩朋成刘刚黄禾王文钦 . 频控阵MIMO雷达中基于稀疏迭代的多维信息联合估计方法. 雷达学报, 2018, 7(2): 194-201. doi: 10.12000/JR16121
    [4] 梁浩崔琛余剑 . 十字型阵列MIMO雷达高精度二维DOA估计. 雷达学报, 2016, 5(3): 254-264. doi: 10.12000/JR16016
    [5] 王珽赵拥军胡涛 . 机载MIMO雷达空时自适应处理技术研究进展. 雷达学报, 2015, 4(2): 136-148. doi: 10.12000/JR14091
    [6] 许京伟朱圣棋廖桂生张玉洪 . 频率分集阵雷达技术探讨. 雷达学报, 2018, 7(2): 167-182. doi: 10.12000/JR18023
    [7] 鲁彦希何子述程子扬刘爽利 . 多目标跟踪分布式MIMO雷达收发站联合选择优化算法. 雷达学报, 2017, 6(1): 73-80. doi: 10.12000/JR16106
    [8] 杨军张群罗迎邓冬虎 . 基于压缩感知的认知雷达多目标跟踪方法. 雷达学报, 2016, 5(1): 90-98. doi: 10.12000/JR14107
    [9] 苏剑宋志勇付强龙照飞 . 单脉冲雷达不可分辨目标与诱饵联合跟踪方法. 雷达学报, 2015, 4(2): 160-171. doi: 10.12000/JR14094
    [10] 尉强刘忠 . 多普勒盲区下基于GM-PHD的雷达多目标跟踪算法. 雷达学报, 2017, 6(1): 34-42. doi: 10.12000/JR16125
    [11] 王祥丽易伟孔令讲 . 基于多目标跟踪的相控阵雷达波束和驻留时间联合分配方法. 雷达学报, 2017, 6(6): 602-610. doi: 10.12000/JR17045
    [12] 张群胡健罗迎陈怡君 . 微动目标雷达特征提取、成像与识别研究进展. 雷达学报, 2018, 7(5): 531-547. doi: 10.12000/JR18049
    [13] 杨琪邓彬王宏强秦玉亮 . 太赫兹雷达目标微动特征提取研究进展. 雷达学报, 2018, 7(1): 22-45. doi: 10.12000/JR17087
    [14] 陈小龙董云龙李秀友关键 . 海面刚体目标微动特征建模及特性分析. 雷达学报, 2015, 4(6): 630-638. doi: 10.12000/JR15079
    [15] 方琳琳周超王锐胡程 . 昆虫目标雷达散射截面积特征辅助跟踪算法. 雷达学报, 2019, 8(5): 598-605. doi: 10.12000/JR19067
    [16] 吴伟尹成友 . 一种用于多目标跟踪的增强型SMC-PHD 滤波算法. 雷达学报, 2012, 1(4): 406-413. doi: 10.3724/SP.J.1300.2012.20094
    [17] 金镇谢良贵文树梁 . 分布式MIMO 雷达单脉冲测角. 雷达学报, 2014, 3(4): 474-479. doi: 10.3724/SP.J.1300.2014.13077
    [18] 何其芳张群罗迎李开明 . 正弦调频Fourier-Bessel变换及其在微动目标特征提取中的应用. 雷达学报, 2018, 7(5): 593-601. doi: 10.12000/JR17069
    [19] 章鹏飞李刚霍超颖殷红成 . 基于双雷达微动特征融合的无人机分类识别. 雷达学报, 2018, 7(5): 557-564. doi: 10.12000/JR18061
    [20] 罗迎张群王国正管桦柏又青 . 基于复图像OMP分解的宽带雷达微动特征提取方法. 雷达学报, 2012, 1(4): 361-369. doi: 10.3724/SP.J.1300.2012.20065
  • 加载中
图(8)表(5)
计量
  • 文章访问数:  638
  • HTML浏览量:  418
  • PDF下载量:  361
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-04-25
  • 录用日期:  2018-06-02
  • 刊出日期:  2018-10-28

基于跟踪脉冲的MIMO雷达多目标微动特征提取

    通讯作者: 罗迎  luoying2002521@163.com
    作者简介: 陈怡君(1989–),女,陕西西安人,现为武警工程大学信息工程学院讲师,研究方向为认知成像与目标识别。E-mail: chenyijun519@126.com
  • ①. 空军工程大学信息与导航学院   西安   710077
  • ②. 中国人民解放军95561部队   拉萨   850000
  • ③. 武警工程大学信息工程学院   西安   710086
基金项目:  国家自然科学基金(61571457,61631019)

摘要: 微动特征是空间目标识别的重要特征信息之一。然而,现有的多功能多输入多输出(Multi-Input Multi-Output, MIMO)雷达通常需要在完成目标搜索和跟踪任务之后为目标微动特征提取分配大量连续的时间资源,导致目标识别实时性能和雷达系统整体工作性能均不高。针对该问题,该文提出了一种基于跟踪脉冲的MIMO雷达多目标微动特征提取方法。首先依据各目标的方位信息对MIMO雷达发射波形进行设计,为不同方向目标同时发射跟踪脉冲;在此基础上,综合考虑目标微动特征提取性能以及目标跟踪性能的需求,对跟踪脉冲的发射时间序列进行优化设计;最后,直接利用窄带跟踪脉冲实现对不同方向目标微动特征的同时提取,无需再为目标微动特征提取分配额外的时间资源,有效提升目标识别实时性和雷达工作效率。仿真实验表明,在信噪比大于–10 dB时,所提方法能够实现多目标微动特征的准确提取,具有良好的有效性和鲁棒性。

English Abstract

    • 微动是指目标或目标部件除主体平动之外的振动、旋转等微小运动[1]。微动特征是目标独一无二的特征,对各类航天器、卫星、弹道导弹、太空碎片等空间目标的分类与识别具有重要意义。近年来,雷达目标的微动特征提取已成为雷达信号处理领域的一个研究热点[27],如文献[5]详细分析了时频分析方法在微动特征提取中的应用;文献[6]利用空间目标的微动特征实现了对目标的高分辨3维成像;文献[7]将参数化稀疏表征的思想引入微动参数的估计中,准确获得了目标微动特征;文献[8]利用微动特征实现了地面车辆目标和行人目标的分类识别;文献[9]利用经验模态分解算法成功提取了锥体目标的进动特征以及目标纵横比;文献[10]提出了一种正弦调频Fourier-Bessel变换用于精确提取目标的微动频率;等等。然而,现有目标微动特征提取方法需要为实现该功能单独发射大量连续的脉冲,而且在一个时刻只能针对一个方向的目标进行微动特征提取,对雷达资源消耗大,目标识别实时性不强,工作效率不高。尤其是在多任务条件下,现有微动特征提取方法难以满足空天目标探测与识别的需求。

      通常,雷达发射窄带脉冲进行目标跟踪,若能利用这些窄带跟踪脉冲同时实现目标跟踪和微动特征提取,这样就无需再为目标微动特征提取分配额外的雷达资源,从而显著提高目标识别实时性和雷达整体工作效率。此外,多输入多输出(Multi-Input Multi-Output, MIMO)雷达具有多个发射和接收天线[11,12],与传统相控阵雷达相比,其各发射阵元可独立发射信号,因而可以根据实际应用灵活地设计发射波形[1315]。通过对发射波形进行优化设计,MIMO雷达能够发射具有特定方向图的波形以实现对不同方向目标的同时观测。因此,本文希望利用MIMO雷达发射端的高自由度优势,同时为不同方向的目标发射跟踪脉冲,并进一步利用跟踪脉冲实现目标微动特征提取,从而节约雷达资源,提高目标识别实时性和雷达工作效率。

      然而,雷达发射窄带脉冲对目标进行跟踪时,跟踪脉冲重复周期通常较大,若基于所发射的跟踪脉冲直接采用时频分析技术进行微动特征提取会产生微多普勒信号频域混叠现象,无法获得理想的结果。因此,现有窄带信号微动特征提取方法大多不再适用,需要研究低重频条件下的微动特征提取方法。此外,由于跟踪脉冲既用于目标跟踪,也用于微动特征提取,基于对目标微动特征提取性能以及目标跟踪性能的综合考虑,对脉冲发射的时间序列进行优化设计可进一步提升雷达系统的整体工作性能。需要说明的是,为了节约雷达资源和保持现有目标跟踪模式基本不变,在对跟踪脉冲发射时间序列进行优化设计时,通常需要保证与现有目标跟踪方式相比,跟踪脉冲重复周期不会发生较大改变。

      综上所述,本文提出一种基于跟踪脉冲的MIMO雷达多目标微动特征提取方法。该方法首先根据各目标方位信息对MIMO雷达发射波形进行设计,从而获得期望的发射方向图,为不同方向的目标同时发射跟踪脉冲;在此基础上,综合考虑目标微动特征提取性能和目标跟踪性能,将微动特征提取性能作为优化目标函数,将目标跟踪性能作为约束条件,建立跟踪脉冲发射时间序列优化模型,实现跟踪脉冲发射时刻的优化;最后,直接利用窄带跟踪脉冲,采用正交匹配追踪(Orthogonal Matching Pursuit, OMP)算法实现目标微动特征提取。

    • MIMO雷达系统的发射阵列为包含 $M$ 个阵元的均匀线阵,阵元间距为 $d$ ,如图1所示。

      图  1  MIMO雷达发射阵列示意图

      Figure 1.  Transmit array of MIMO radar

      $m$ 个阵元发射的信号为:

      ${s_m}(t) = {x_m}(t){{\mathop{\rm e}\nolimits} ^{{\rm{j2}}{{π}} {f\!_{\rm{c}}}t}}\;,\quad 0 \le t \le {T_{\rm{p}}}$

      其中, ${x_m}(t)$ 为第 $m$ 个阵元发射的基带信号, ${f\!_{\rm{c}}}$ 为信号载频, ${T_{\rm{p}}}$ 为脉冲宽度。则远场 $\theta $ 方向处的合成信号可表示为:

      $ \begin{aligned} s(\theta ,t) =& \sum\limits_{m = 1}^M {{s_m}} \left(t + \frac{{(m - 1)d\sin \theta }}{c}\right)\\ = & \sum\limits_{m \!=\! 1}^M {{x_m}} \left(t \!+\! \frac{{(m - 1)d\sin \theta }}{c}\right)\\ & \cdot {{\mathop{\rm e}\nolimits} ^{{\rm{j}}2{{π}} {f\!_{\rm{c}}} \left(t + \scriptsize\displaystyle\frac{{(m - 1)d\sin \theta }}{c} \right)}} \end{aligned} $

      其中, $c$ 为光速。

      在实际应用中,考虑离散基带发射信号,即 ${x_m}(l) = {x_m}(t){|_{t = (l - 1){T_{\rm{s}}}}}, \ l = 1,2, ·\!·\!· ,L$ ,其中 $L$ 为脉冲内的采样点数, ${T_{\rm{s}}}$ 为采样间隔,则脉冲宽度为 ${T_{\rm{p}}} = {T_{\rm{s}}} \cdot L$ 。则远场 $\theta $ 方向处的合成信号可表示为:

      ${{S}}(\theta ) = {{{a}}^{\rm H}}(\theta ){{X}}$

      其中, ${{a}}(\theta ) \!=\! {[1\; {\operatorname{e} ^{{\rm{j}}2{{π}} d\sin \theta /\lambda }}\, ·\!·\!· \, {\operatorname{e} ^{{\rm{j}}2{π} (M - 1)d\sin \theta /\lambda }}]^{\rm H}}$ 为发射导向矢量, $\lambda $ 为发射信号波长。 ${{X}} = {[{{x}}_1^{\rm H}\,{{x}}_2^{\rm H}\, ·\!·\!·\,{{x}}_M^{\rm H}]^{\rm H}}$ 为阵列发射的信号矩阵, ${{{x}}_m}$ 为第 $m$ 个阵元发射的基带信号且 ${{{x}}_m} \in {\mathbb{C}^{1 \times L}}$

      远场 $\theta $ 方向处的合成信号的功率为:

      $P(\theta ) = \frac{1}{L}{{S}}(\theta ){{{S}}^{\rm{H}}}(\theta ) = {{{a}}^{\rm{H}}}(\theta ){{Ra}}(\theta )$

      其中, ${{R}} = {{X}}{{{X}}^{\rm{H}}}/L$ 为阵列发射信号的协方差矩阵。 $\{ P({\theta _k})\} _{k = 1}^K$ 反映了发射信号功率在空域的分布,即为MIMO雷达的发射方向图, $\{ {\theta _k}\} _{k = 1}^K$ 为离散化方位角, ${\theta _k}$ 处的半波束宽度记为 $B({\theta _k})$ 。因此,可以通过对阵列发射信号的协方差矩阵 ${{R}}$ 进行设计以获取实际所需的发射方向图。

      根据需要同时观测的目标数量和方位来确定期望发射方向图 $\{ {P_{\rm{d}}}({\theta _k})\} _{k = 1}^K$ 。假设需要同时对H个目标发射跟踪脉冲,第 $h$ 个目标的方位角为 ${\theta _h}$ ,第 ${h_1}$ 个目标和第 ${h_2}$ 个目标之间的方位角之差记为 $\Delta \theta ({h_1},{h_2})$ ,则期望发射方向图 $\{ {P_{\rm{d}}}({\theta _k})\} _{k = 1}^K$ 应满足

      $\left. {\begin{array}{*{20}{l}} {\forall h \in 1, ·\!·\!· ,H,}\quad {{P_{\rm{d}}}({\theta _h}) \!=\! 1}\\ {\forall {h_1},{h_2} \in 1, ·\!·\!· ,H,\;{h_1} \!\ne\! {h_2},} \quad \!\!\!{B({\theta _h}) \!<\! \Delta \theta ({h_1},{h_2})} \end{array}}\!\!\! \right\} $

      采用方向图匹配方法对MIMO雷达发射波形进行设计。利用所设计方向图与期望方向图逼近的思想,建立如下优化模型:

      $ \left. {\begin{array}{*{20}{l}} {\mathop {\min }\limits_{\alpha ,{{R}}} \;\;\sum\limits_{k = 1}^K {{w_k}} {{\Bigr[\alpha {P_{\rm{d}}}({\theta _k}) - {{{a}}^{\rm H}}({\theta _k}){{Ra}}({\theta _k})\Bigr]}^2}}\\ {{\mathop{\rm s}\nolimits} .{\rm t}.\;\;\;{R_{mm}} = 1,\;\;\;m = 1,2, ·\!·\!· ,M}\\ \quad\quad\ {{{R}} \ge 0} \end{array}} \!\!\!\right\} $

      其中, ${w_k} \ge 0$ 为方位角 ${\theta _k}$ 处的逼近加权, $\alpha $ 为尺度因子, ${R_{mm}}$ 为协方差矩阵 ${{R}}$ 的第 $(m,m)$ 个元素。约束条件保证了各阵元的发射功率相同, ${{R}} \ge 0$ 保证了雷达发射信号协方差矩阵的半正定性。该优化模型可以用凸优化工具包进行求解。

      通过对上述优化模型进行求解可以得到雷达发射信号的协方差矩阵 ${{R}}$ ,但阵列发射的信号矩阵 ${{X}}$ 仍然未知,无法获取每个阵元具体的发射信号。为了获得恒模发射波形,基于循环算法对协方差矩阵匹配下的发射波形进行优化

      $ \left. {\begin{array}{*{20}{l}} {\min \;\;\left\|{{X}} - {{{R}}^{1/2}}{{U}}\right\|}\\ {{\mathop{\rm s}\nolimits} .{\rm t}.\;\;\;|X(m,l)| = 1/M}\\ \quad\quad\ {m = 1,2, ·\!·\!· ,M,l = 1,2, ·\!·\!· ,L} \end{array}}\right\} $

      其中, ${{{R}}^{1/2}}$ 为协方差矩阵的Hermite均方根, ${{U}}$ 为半正交矩阵,维数为 $M \times L$

      通过对发射信号矩阵 ${{X}}$ 和半正交矩阵 ${{U}}$ 的交替求解,可以最终得到雷达的发射信号矩阵,从而在空域合成所需的发射方向图以实现对不同方向目标同时发射跟踪脉冲。

    • 以旋转目标为例,基于第2节中设计的发射波形对目标的窄带雷达回波信号进行分析。设在第 $\{ h\} _{h = 1}^H$ 个目标方向,即在远场 ${\theta _h}$ 方向处的旋转目标为散射点模型且已完成了平动补偿,如图2所示。经过时间 $ \tau $ 的运动,目标上第 $i$ 个散射点 ${E_i}$ 以角速度 ${{w}}$ 转动到 ${E'\!\!_i}$ ,目标的旋转矢量为 ${{w}} = {({\omega _x}\ {\omega _y}\ {\omega _z})^{\rm T}}$ ,在起始时刻目标散射点 ${E_i}$ 的坐标向量为 ${{{r}}_0} = $ $ {({r_{x0}}\ {r_{y\,0}}\ {r_{z0}})^{\rm T}}$ , ${E'\!\!_i}$ 的坐标向量为 ${{r}} = {({r_x}\ {r_y}\ {r_z})^{\rm T}}$ ,本地坐标系 $(x,y,z)$ 中的向量乘以欧拉旋转矩阵 ${{{R}}_{{\rm{init}}}}$ 就可以变换为参考坐标系 $(X,Y,Z)$ 中的向量,即 $ {{\overset{\frown}{w}}}=$ ${({\omega _X}\ {\omega _Y}\ {\omega _Z})^{\rm T}} \!=\!\! {{{R}}_{{\rm{init}}}}{{w}}$ , ${{{\overset{\frown}{r} }}_0} =\!\! {({r_{X0}}\ {r_{Y\,0}}\ {r_{Z0}})^{\rm T}} \!\!=\!\! {{{R}}_{{\rm{init}}}}{{{r}}_0}$ 。旋转运动可以用一个参考坐标系中的旋转矩阵 ${{{R}}_{{\rm{rotating}}}}$ 来描述,即有 ${{\overset{\frown}{r} }} = {{{R}}_{{\rm{rotating}}}}{{{\overset{\frown}{r} }}_0}$

      图  2  旋转目标的几何模型

      Figure 2.  Geometry of a rotating target

      由式(2)可知,在 ${\theta _h}$ 方向处,目标上第 $i$ 个散射点 ${E_i}$ 的慢时间基带回波信号为:

      $\begin{align} {s_{ri}}({\theta _h},\tau )=& {\sigma _i}\sum\limits_{m = 1}^M {{x_m}} \left(t \!+\! \frac{{(m \!-\! 1)d\sin {\theta _h}}}{c} \!-\! \frac{{2{R_i}(\tau )}}{c}\right) \\ & \cdot{{\mathop{\rm e}\nolimits} ^{{\rm{j}}2{{π}} {f\!_{\rm{c}}}\left(\scriptsize\displaystyle\frac{{(m - 1)d\sin {\theta _h}}}{c} - \frac{{2{R_i}(\tau )}}{c}\right)}}, \tau \in [0,{T\,_{\rm{c}}}] \end{align}$

      其中, ${\sigma _i}$ 为该散射点的反射系数, ${R_i}(\tau )$ 为瞬时斜距, ${T\,_{\rm{c}}}$ 为方位向积累时间。在远场条件下,基于平面波近似, ${R_i}(\tau )$ 可写为:

      ${R_i}(\tau ) = \left\| {{{{R}}_{{\rm{rotating}}}}{{{{\overset{\frown}{r} }}}_0}} \right\|$

      其中, ${{{R}}_{{\rm{rotating}}}} = \exp \left(\left\| {{{\overset{\frown}{w} }}} \right\|{{\overset{\frown}{w}}' }\tau \right)$ ,旋转角速度单位向量为 ${{\overset{\frown}{{w}}' }} = {{\overset{\frown}{w} }}/\left\| {{{\overset{\frown}{w} }}} \right\| = {({\omega '\!\!_X}\ {\omega '\!\!_Y}\ \ {\omega '\!\!_Z})^{\rm T}}$ ,则

      ${\widehat{ w'}} = \left[ {\begin{array}{*{20}{c}} 0&{ - {{\omega '}\!\!_Z}}&{{{\omega '}\!\!_Y}} \\ {{{\omega '}\!\!_Z}}&0&{ - {{\omega '}\!\!_X}} \\ { - {{\omega '}\!\!_Y}}&{{{\omega '}\!\!_X}}&0 \end{array}} \right]$

      ${\theta _h}$ 方向处目标由J个散射点组成,则可以得到目标的基带回波信号为:

      $ \begin{align} {s_{\rm{r}}}({\theta _h},\tau ) =& \sum\limits_{i = 1}^J {\sigma _i}\\ &\cdot \sum\limits_{m = 1}^M {{x_m}} \left(t + \frac{{(m - 1)d\sin {\theta _h}}}{c} - \frac{{2{R_i}(\tau )}}{c}\right) \\ & \cdot {\mathop{\rm e}\nolimits} ^{{\rm{j}} 2{{π}}{f\!_{\rm{c}}}\left(\scriptsize\displaystyle\frac{{(m - 1)d\sin {\theta _h} - 2{R_i}(\tau )}}{c}\right)} \end{align} $

      近年来,OMP算法被广泛应用在稀疏信号的分解与重构中。因此,本节中利用OMP算法来实现雷达目标微动特征提取。由OMP算法原理可知,根据信号本身的固有特性构造字典中的原子,可以实现信号的分解与重构[16,17]。从式(11)中可以看出,目标各微动点回波由参数 $\left( {{{w,}}{{{r}}_0}} \right)$ 确定,因此,通过设定 ${{w}} = {({\omega _x}\ {\omega _y}\ {\omega _z})^{\rm T}}$ , ${{{r}}_0} = {({r_{x0}}\ {r_{y0}}\ {r_{z0}})^{\rm T}}\,$ 的不同取值来构造原子集D中的原子

      $ \begin{align} d \left({{{w}}_{{n_x}{n_y}{n_z}}},{{{r}}_0}_{{n_{rx}}{n_{ry}}{n_{rz}}}\right) =& \sum\limits_{m = 1}^M {{x_m}} \left(t + \frac{{(m - 1)d\sin {\theta _h}}}{c} - \frac{{2\left\| {\exp \left({{{R}}_{{\rm{init}}}}{{w} _{{n_x}{n_y}{n_z}}}n\Delta \tau \right){{{R}}_{{\rm{init}}}}{{{r}}_0}_{{n_{rx}}{n_{ry}}{n_{rz}}}} \right\|}}{c} \right. \\ &\left. \cdot \exp \left( {{\rm{j}}2{{π}} {f\!_{\rm{c}}}\left( {\frac{{(m - 1)d\sin {\theta _h} - 2\left\| {\exp ({{{R}}_{{\rm{init}}}}{{{w}}_{{n_x}{n_y}{n_z}}}n\Delta \tau ){{{R}}_{{\rm{init}}}}{{{r}}_0}_{{n_{rx}}{n_{ry}}{n_{rz}}}} \right\|}}{c}} \right)} \right)\right), \\ & \quad\quad\quad\quad\quad\quad n = 1, ·\!·\!· ,N \end{align} $

      其中

      $ \begin{align} {{{w}}_{{n_x}{n_y}{n_z}}} =& {\left[{\omega_{x{n_x}}}\quad {\omega _{y{n_y}}}\quad {\omega _{z{n_z}}}\right]^{\rm T}} = \left[\!{\varOmega _{x\min }} \!+\! ({n_x} - 1) \cdot \frac{{{\varOmega _{x\max }} \!-\! {\varOmega _{x\min }}}}{{{N_x}}} \right.\\ & {\varOmega _{y\min}}+ ({n_y}- 1)\cdot \left. \frac{{{\varOmega _{y\max }} \!-\! {\varOmega _{y\min }}}}{{{N_y}}} \quad\quad {\varOmega _{z\min }} \!+\! ({n_z} \!-\! 1) \cdot \frac{{{\varOmega _{z\max }} \!-\! {\varOmega _{z\min }}}}{{{N_z}}}\! \right]^{\rm T}, \\ & \quad\quad\quad {n_x} =1,2, ·\!·\!· ,{N_x};\;{n_y} = 1,2, ·\!·\!· ,{N_y};\;{n_z} = 1,2, ·\!·\!· ,{N_z} \end{align} $

      $ \begin{aligned} {{{r}}_0}_{{n_{rx}} {n_{ry}} {n_{rz}}} =& {\left[{r_{x0{n_{rx}}}}\quad {r_{y0{n_{ry}}}}\quad {r_{z0{n_{rz}}}}\right]^{\rm T}}= \left[{R_{x0\min }} + ({n_{rx}} - 1) \cdot \frac{{{R_{x0\max }} - {R_{x0\min }}}}{{{N_{rx}}}} \right. \\ & {R_{y0\min}}+ ({n_{ry}} - 1) \cdot \left. \frac{{{R_{y\,0\max }}- {R_{y\,0\min }}}}{{{N_{ry}}}}\quad\quad {R_{z0\min }} + ({n_{rz}} - 1) \cdot \frac{{{R_{z0\max }} - {R_{z0\min }}}}{{{N_{rz}}}} \right]^{\rm T} , \\ & \quad\quad\quad{n_{rx}} = 1,2, ·\!·\!· ,{N_{rx}};\;{n_{ry}} = 1,2, ·\!·\!· ,{N_{ry}};\;{n_{rz}} = 1,2, ·\!·\!· ,{N_{rz}} \end{aligned} $

      在式(12)到式(14)中, $([{\varOmega _{x\min }},{\varOmega _{x\max }}] ,$ $ [{\varOmega _{y\min }},{\varOmega _{y\max }}],[{\varOmega _{z\min }},{\varOmega _{z\max }}])$ , $([{R_{x0\min }},{R_{x0\max }}] ,$ $ [{R_{y\,0\min }},{R_{y\,0\max }}], [{R_{z0\min }},{R_{z0\max }}])$ 分别表示 ${{w}}$ ${{{r}}_0}$ 的取值范围; $\left(\displaystyle\frac{{{\varOmega _{x\max }} - {\varOmega _{x\min }}}}{{{N_x}}},\;\displaystyle\frac{{{\varOmega _{y\max }} - {\varOmega _{y\min }}}}{{{N_y}}},\right.$ $ \left.\displaystyle\frac{{{\varOmega _{z\max }} \!-\! {\varOmega _{z\min }}}}{{{N_z}}}\!\!\right)$ $\left(\!\!\displaystyle\frac{{{R_{x0\max }} \!-\! {R_{x0\min }}}}{{{N_{rx}}}}\!,\!\frac{{{R_{y\,0\max }} \!- \!{R_{y\,0\min }}}}{{{N_{ry}}}}\!,\right.$ $ \left.\displaystyle\frac{{{R_{z0\max }} \!-\! {R_{z0\min }}}}{{{N_{rz}}}}\right)$ 分别为 ${{w}}$ ${{{r}}_0}$ 的搜索步进增量。脉冲发射时间间隔 $\Delta \tau {\rm{ = }}{{1/F}_2}$ , ${F_2}$ 为微多普勒信号不存在频域混叠现象时的脉冲重复频率,总脉冲数 $N = {F_2} \cdot {T_{\rm{c}}}$ 。利用跟踪脉冲进行微动特征提取时,设跟踪脉冲重复频率为 ${F_1}$ , $Q \!=\!\! {F_1} \cdot {T_{\rm{c}}} + 1$ 为实际发射跟踪脉冲数。记 $\varepsilon \!=\! \left\lfloor {{{{F_2}}/{{F_1}}}} \right\rfloor $ ,则目标跟踪脉冲回波信号 ${s_{\rm{r}}}({\theta _h},\tau )$ 可表示为:

      ${{s}_{\rm{r}}} = {Φ} {D}\beta $

      其中, $\beta $ 表示脉冲重复频率取 ${F_2}$ 时回波信号在原子集D上的投影系数。观测矩阵 ${Φ} $ 为一个0-1矩阵,具体形式为:

      $ {Φ} \!=\! \left\{ {{\phi _{i,i'}}} \right\} \!=\! \left\{ {\begin{array}{*{20}{l}} \!\!\!\!{1,}&{\left\{ {\left(i,i'\right)|i' \!=\! \varepsilon \cdot i} \right\}}\\ \!\!\!\!{0,}&{{\rm{others}}} \end{array}} \right., \,i \!=\! 1, ·\!·\!·, Q $

      通常,目标散射分布具有空域稀疏特性,因此 $\beta $ 中只有少量非零元素。显然,式(15)是一个典型的压缩感知模型,当满足有限等距性质(Restricted Isometry Property, RIP)时,通过求解一个最优化问题即可实现对 $\beta $ 的重构

      $\min {\left\| \beta \right\|_1},\;\;\;\;{\rm{s}}{\rm{.t}}{\rm{.}}\;\;\;{Φ} {D}\beta = {{s}_{\rm{r}}}\;$

      在压缩感知重构过程中,感知矩阵 ${Θ} = {Φ} {{D}}$ 的正交性越好,算法的重构概率越高,微动特征的提取精度也就越高。实际上,观测矩阵 ${{Φ}}$ 的具体形式与脉冲信号的发射时间序列一一对应。基于对微动特征提取性能以及目标跟踪性能的考虑建立观测矩阵优化模型(本质上即跟踪脉冲发射时间序列优化模型):

      $ \left. \begin{array}{l} \mathop {\min }\limits_{Φ} {\left\| {{({{Φ} {{D}})}^{\rm{H}}}({Φ} {{D}}) - {{I}}} \right\|_2}\\ {\rm{s}}.{\rm{t}}.\;\;\;{Φ} (i,j) = 0,1,\;\;\;i = 1,2, ·\!·\!· ,Q;\\ \quad\quad\quad\quad\quad\quad\quad\quad\quad j = 1,2, ·\!·\!· ,N\\ \;\;\;\;\;\;\; \sum\limits_{i = 1}^Q {\sum\limits_{j = 1}^N {{Φ} (i,j) = } } Q\\ \;\;\;\;\;\;\; \sum\limits_{i = 1}^Q {{Φ} (i,j)} \le 1,\;\;j = 1,2, ·\!·\!· ,N\\ \;\;\;\;\;\;\; {ρ} (1) = 1\\ \;\;\;\;\;\;\; {ρ} (\varepsilon \cdot Q) = 1\\ \;\;\;\;\;\;\; \sum\limits_{i = - \kappa }^{\kappa + 1} {{ρ} (\varepsilon \cdot \mu + i)} \ge 1,\mu = 1,2, ·\!·\!· ,Q - 1 \end{array}\!\!\!\! \right\}$

      其中,I为维数为 ${N_{{{r}}0}}{N_{{w}}}$ 的单位阵, ${N_{{{r}}0}} = {N_{rx}} \cdot {N_{ry}} $ $ \cdot{N_{rz}} $ , ${N_{{w}}} = {N_x} \cdot {N_y} \cdot {N_z}$ 。第1个约束条件限制了观测矩阵 ${{Φ}}$ 为0-1矩阵;第2个约束条件限制了实际发射总脉冲数为Q,即保证发射脉冲数与目标跟踪所需的脉冲数保持一致,不额外发射脉冲;第3个约束条件保证了在每个时刻最多只发射一个脉冲;第4个约束条件和第5个约束条件保证了观测时间长度为 ${T_{\rm{c}}}$ ,即在第1个时刻和最后一个时刻必须发射脉冲,其中 ${ρ} = \displaystyle\sum\nolimits_{i = 1}^Q {{{Φ} _i}} $ , ${{Φ} _i}$ 为观测矩阵 ${{Φ}}$ 中的第i行;第6个约束条件限制了实际脉冲发射时刻仅在时间窗w允许的范围内进行调整,即保证与现有目标跟踪方式相比,跟踪脉冲重复周期不发生较大改变,从而避免目标失跟,其中, $\kappa \!=\! \left\lfloor {{w/{{T_2}}}} \right\rfloor $ , ${T_2} =$ $ 1/{F_2} $

      针对式(18)所示的优化模型,本文采用遗传算法对观测矩阵 ${{Φ}}$ 进行求解。在得到最优观测矩阵 ${{Φ}}$ 后,即可得到所设计的跟踪脉冲最优发射时间序列 ${{{t}}_{\rm{m}}} = {γ} \cdot {T_2}$ ,其中 ${γ} \!=\! \{ j|{ρ} (j) \!=\! 1,\;j \!=\! 1,2, ·\!·\!· ,N\} $ 。此时,目标基带回波信号和观测矩阵可分别表示为:

      $ {s_{\rm{s}}}({\theta _h},\tau ) =\!\! \sum\limits_{i = 1}^J {{\sigma _i}\sum\limits_{m = 1}^M {{{{x}}_m}} } \!\cdot\! {\operatorname{e} ^{{\rm{j}}2{π} {f\!_{\rm{c}}}\left(\scriptsize\displaystyle\frac{{(m \!-\! 1)d\sin {\theta _h}}}{c} \!\!-\!\! \frac{{2{R_i}({{{t}}_m})}}{c}\right)}} $

      ${{Φ} _{\rm{s}}} \!=\! \left\{ {{\phi _{i,i'}}} \right\} \!=\! \left\{ {\begin{array}{*{20}{l}} \!\!\!\!{1,}&{\left\{ {\left(i,i'\right)|i' \!=\! {γ} (i)} \right\}}\\ \!\!\!\!{0,}&{{\rm{others}}} \end{array}}\!\!\! \right., i = 1, ·\!·\!·, Q $

      利用OMP算法求解优化模型 $\min {\left\| \beta \right\|_1},$ $ \;{\rm{s}}{\rm{.t}}{\rm{.}}$ ${{Φ} _{\rm{s}}}{D}\beta = {{s}_{\rm{s}}}\; $ ,根据 $\beta $ 中非零元素所对应的原子 ${{w}}$ ${{{r}}_0}$ 的取值,即可实现对目标微动特征参数 ${{w}}$ ${{{r}}_0}$ 的有效提取。

      综上所述,图3给出了所提微动特征提取方法的完整流程。

      图  3  算法流程

      Figure 3.  Flow chart of the algorithm

    • 设MIMO雷达系统包含 $M = 20$ 个发射阵元,阵元间距 $d = c/\left( {2{f\!_{\rm{c}}}} \right)$ ,发射信号载频 ${f\!_{\rm{c}}} = 1\;{\rm{GHz}}$ ,观测时间 ${T_{\rm{c}}} = 3\;{\rm{s}}$ ,雷达发射信号脉宽 ${T_{\rm{p}}} = $ 1 μs,每个脉冲内的采样点数 $L = 400$ ,离散化方位角总数 $K = 181$ ,时间窗 $w = 20\;{\rm{ms}}$ ,跟踪脉冲发射重复频率 ${F_1} = 8\;{\rm{Hz}}$ ,微多普勒信号不存在频域混叠现象时的脉冲重复频率 ${F_2} = 800\;{\rm{Hz}}$

      假设在 $ - {30^ \circ }$ ${0^ \circ }$ 两个方向分别存在旋转目标,其参数如表1所示,其中目标参考点坐标以雷达所在位置为坐标原点,目标散射点坐标以目标参考点为原点。

      目标方位(°) 参考点坐标(km) 旋转角速度 $\left( {\rm{rad}} \cdot {{\rm{s}}^{{\rm{ - 1}}}}\right) $ 目标散射点坐标(m)
      –30 (15, 26, 0) ${\rm{(}}{π} {\rm{,}}{π} {\rm{/2,}}{π} {\rm{)}}$ $(0,0,0),({\rm{1}},{\rm{1}},{\rm{0}}),( - {\rm{1}}, - {\rm{1}},{\rm{0}})$
      0 (0, 40, 0) ${\rm{(0,0,}}{π} {\rm{)}}$ $(2,0,0),( - 2,0,0),(0, - 2,0),(0,2,0)$

      表 1  目标参数

      Table 1.  The parameters of the targets

      所设计的发射方向图和期望方向图如图4所示。期望方向图由采用30 dB切比雪夫权的传统波束的主瓣构成。根据目标的方位信息,期望发射方向图的主瓣区域设置为 $[ - {40^ \circ }, - {20^ \circ }] \cup [ - {10^ \circ },{10^ \circ }]$ 。由图可以看出所设计波形在 $ - {30^ \circ }$ ${0^ \circ }$ 两个目标方向形成波束,与期望发射方向图逼近,利用该发射波形可实现对 $ - {30^ \circ }$ 方向和 ${0^ \circ }$ 方向目标的同时观测。

      图  4  发射方向图

      Figure 4.  The transmit beampattern

      图5给出了在满足时间窗约束的条件下,采用文献[18]中所提随机稀疏采样多普勒解模糊算法重构出的信号时频分布图。由表2所示的微动特征提取结果可以看出该方法成功提取出了 $ - {30^ \circ }$ 方向目标的微动特征,而对于 ${0^ \circ }$ 方向的目标,该方法提取出了错误的特征。这是因为文献[18]中的脉冲发射时间序列是随机设定的,而本文所提方法考虑到脉冲发射时间序列是与观测矩阵 ${{Φ}}$ 的具体形式一一对应的,因此,将脉冲发射时间序列优化问题转化为观测矩阵优化问题,以最大化感知矩阵正交性为目标,对观测矩阵 ${{Φ}}$ 进行优化,本质上实现了脉冲发射时间序列的优化。因此,与文献[18]方法相比,本文所提方法的信号重构性能显著提高。对于0°方向目标,在脉冲数较少时,文献[18]方法得到的感知矩阵的列相关性较大,从而无法正确提取目标微动特征,而本文所提方法能够最大化感知矩阵列之间的正交性,因此能够有效提取目标微动特征。

      图  5  无噪声条件下文献[18]方法重构的时频分布

      Figure 5.  The reconstructed time-frequency distribution using method in Ref. [18] without noise

      目标方向(°) 特征参数序号 目标散射点坐标 ${{r}_0}$ (m) 旋转角速度 ${{w}}$ $\left( {\rm{rad}} \cdot {{\rm{s}}^{{\rm{ - 1}}}}\right) $ 系数A
      –30 1 $(0.1 \ 0 \ - 0.1)$ $(3.14 \ 1.57 \ 3.14)$ 30.40
      2 $({\rm{1 \ 1 \ 0}})$ $(3.14 \ 1.57 \ 3.14)$ 27.07
      3 $( - {\rm{1.0 }}\ - {\rm{1.1}} \ {\rm{ 0}})$ $(3.14 \ 1.57 \ 3.14)$ 25.11
      4 $({\rm{0}}{\rm{.3 }}\ 0\ {\rm{ 0}})$ $(3.14 \ 1.57 \ 3.14)$ 5.96
      5 $(1.5 \ - {\rm{1.2}}\ {\rm{0}})$ $(2.83 \ 1.57 \ 3.14)$ 4.29
      0 1 $( - 1.8 \ 0.1 \ 0)$ $(0.31 \ 0 \ 2.51)$ 29.15
      2 $(0.1 \ - 1.9 \ 0)$ $(0.31 \ 0 \ 2.51)$ 28.10
      3 $(0.5 \ 2.4 \ 0)$ $(0.31 \ 0 \ 2.51)$ 24.27
      4 $(2.6 \ - 0.2 \ 0)$ $(0.31 \ 0 \ 2.51)$ 23.31
      5 $( - 1.9 \ - {\rm{0.4 }}\ {\rm{0}})$ $(0 \ 0.31 \ 2.83)$ 3.15

      表 2  文献[18]方法目标微动特征提取结果

      Table 2.  The micro-motion features of targets obtained by method in Ref. [18]

      图6(a)图7(a)给出了在无噪声环境下,利用本文所设计的脉冲发射时刻 ${{{t}}_m}$ 对目标回波信号进行采样后,采用OMP算法重构出的信号时频分布图,由表3所示的微动特征提取结果可以看出算法成功提取出了 $ - {30^ \circ }$ 方向和 ${0^ \circ }$ 方向目标中旋转散射点的微动特征,微动特征参数的取值均与真实值接近,证明了本文所提算法的有效性。

      目标方向(°) 特征参数序号 目标散射点坐标 ${{r}_0}$ (m) 旋转角速度 ${{w}}$ ( ${\rm{rad}} \cdot {{\rm{s}}^{{\rm{ - 1}}}}$) 系数A
      –30 1 $(0.1 \ 0 \ - 0.1)$ $(3.14 \ 1.57 \ 3.14)$ 27.23
      2 $({\rm{1 \ 1 \ 0}})$ $(3.14 \ 1.57 \ 3.14)$ 26.47
      3 $( - {\rm{1.0}} \ - {\rm{1.1}} \ {\rm{ 0}})$ $(3.14 \ 1.57 \ 3.14)$ 26.01
      4 $( - 1.3 \ - {\rm{1.2}}\ {\rm{0}})$ $(2.83 \ 1.57 \ 3.14)$ 6.40
      5 $({\rm{0}}{\rm{.3 }} \ 0 \ {\rm{ 0}})$ $(3.14 \ 1.88 \ 3.14)$ 6.06
      0 1 $(2.0 \ 0.1 \ 0)$ $(0 \ 0 \ 3.14)$ 28.63
      2 $( - 1.9 \ 0 \ 0)$ $(0 \ 0 \ 3.14)$ 28.39
      3 $(0.1 \ - 2.0 \ 0)$ $(0 \ 0 \ 3.14)$ 27.71
      4 $(0 \ 1.9 \ 0)$ $(0 \ 0 \ 3.14)$ 27.19
      5 $(1.8 \ 0.1 \ 0.3)$ $(0 \ 0.31 \ 2.83)$ 6.72

      表 3  无噪声条件下所提方法目标微动特征提取结果

      Table 3.  The micro-motion features of targets obtained by the proposed method without noise

      图6(b)图7(b)给出了在信噪比(Signal-to-Noise Ratio, SNR)为–10 dB的环境下,采用所提方法重构出的信号时频分布图。由表4所示的微动特征提取结果可以看出算法在低信噪比环境中仍成功提取出了 $ - {30^ \circ }$ 方向和 ${0^ \circ }$ 方向目标中旋转散射点的微动特征,证明了本文算法在低信噪比环境下的鲁棒性。当信噪比进一步降为–15 dB时,由表5可以看出,由于噪声的影响,图6(c)图7(c)所示的时频分布图与图6(a)图7(a)相比能量聚集性也有所下降,提取出4个散射点的微动特征之外还提取出了虚假的微动参数。

      目标方向(°) 特征参数序号 目标散射点坐标 ${{r}_0}$ (m) 旋转角速度 ${{w}}$ $\left( {\rm{rad}} \cdot {{\rm{s}}^{{\rm{ - 1}}}}\right) $ 系数A
      –30 1 $(0.1 \ 0 \ - 0.1)$ $(3.14 \ 1.57 \ 3.14)$ 26.83
      2 $({\rm{1 \ 1 \ 0}})$ $(3.14 \ 1.57 \ 3.14)$ 25.48
      3 $( - {\rm{1.0}} - {\rm{1.1}}\ {\rm{ 0}})$ $(3.14 \ 1.57 \ 3.14)$ 25.09
      4 $( - 1.3 \ - {\rm{1.2}}\ {\rm{ 0}})$ $(2.83 \ 1.57 \ 3.14)$ 7.78
      5 $({\rm{0}}{\rm{.3 }}\ 0 \ {\rm{ 0}})$ $(3.14 \ 1.88 \ 3.14)$ 7.31
      0 1 $(2.0 \ 0.1 \ 0)$ $(0 \ 0 \ 3.14)$ 27.52
      2 $( - 1.9 \ 0 \ 0)$ $(0 \ 0 \ 3.14)$ 27.28
      3 $(0.1 \ - 2.0 \ 0)$ $(0 \ 0 \ 3.14)$ 26.59
      4 $(0 \ 1.9 \ 0)$ $(0 \ 0 \ 3.14)$ 26.08
      5 $(1.8 \ 0.1 \ 0.3)$ $(0 \ 0.31 \ 2.83)$ 9.15

      表 4  SNR=–10 dB条件下所提方法目标微动特征提取结果

      Table 4.  The micro-motion features of targets obtained by the proposed method when SNR=–10 dB

      目标方向(°) 特征参数序号 目标散射点坐标 ${{r}_0}$ (m) 旋转角速度 ${{w}}$ $\left( {\rm{rad}} \cdot {{\rm{s}}^{{\rm{ - 1}}}}\right) $ 系数A
      –30 1 $(0.2 \ 0.1 \ - 0.1)$ $(3.14 \ 1.57 \ 3.14)$ 25.73
      2 $({\rm{0.9}}{\rm{ \ 1.0 \ 0}})$ $(3.14 \ 1.57 \ 3.14)$ 23.27
      3 $( - {\rm{1}}{\rm{.2 }}\ - {\rm{1.0}}\ {\rm{0.1}})$ $(3.14 \ 1.57 \ 3.14)$ 22.16
      4 $( - 1.5 \ - {\rm{1.3}}\ {\rm{ 0}})$ $(2.83 \ 1.57 \ 3.14)$ 15.04
      5 $({\rm{0}}{\rm{.8 }} \ 0 \ {\rm{ 0}}{\rm{.5}})$ $(2.83 \ 2.19 \ 3.14)$ 10.58
      0 1 $(2.2 \ - 0.1 \ 0)$ $(0 \ 0 \ 3.14)$ 26.81
      2 $( - 2.0 \ 0.1 \ 0)$ $(0 \ 0 \ 3.14)$ 25.35
      3 $(0.2 \ - 1.8 \ 0)$ $(0 \ 0 \ 3.14)$ 25.19
      4 $(0.1 \ 1.8 \ 0.2)$ $(0 \ 0 \ 3.14)$ 24.07
      5 $(1.7 \ 0.3 \ 0.4)$ $(0 \ 0.31 \ 2.83)$ 13.15

      表 5  SNR=–15 dB条件下所提方法目标微动特征提取结果

      Table 5.  The micro-motion features of targets obtained by the proposed method when SNR=–15 dB

      图  7  0°方向所提方法重构的时频分布

      Figure 7.  The reconstructed time-frequency distribution in 0° direction using the proposed method

      图  6  –30°方向所提方法重构的时频分布

      Figure 6.  The reconstructed time-frequency distribution in –30° direction using the proposed method

      为了对算法性能进行定量分析,定义两个性能指标:真实散射点的特征参数提取精度E;虚假散射点与真实散射点的能量比 $\xi $ 。假设第h个目标上共有 ${g_h}$ 个散射点,各散射点的真实坐标和旋转角速度分别为 ${{{r}}_{0,{g_h}}}$ ${{{w}}_{{g_h}}}$ ,所提方法获得的相应散射点的特征提取结果分别为 ${{{r}'}\!_{0,{g_h}}}$ ${{{w}'}\!\!_{{g_h}}}$ ,则特征参数提取精度E定义如下:

      $E{\rm{ =\! }}\frac{1}{H}\sum\limits_{h = 1}^H {\sum\limits_{g \!=\! 1}^{{g_h}} {\frac{1}{{{g_h}}}\left( {\frac{{{{\left\| {{{{r}}_{0,g}}\! -\! {{{{r}'}}\!_{0,g}}} \right\|}_2}}}{{{{\left\| {{{{r}}_{0,g}}} \right\|}_2}}} \!+\! \frac{{{{\left\| {{{{w}}_g} \!-\! {{{{w}'}}\!\!_g}} \right\|}_2}}}{{{{\left\| {{{{w}}_g}} \right\|}_2}}}} \right)} } $

      假设对于第h个目标,所提方法共提取出 ${B_h}$ 组微动特征参数,相应的强度系数记为 ${A_{{B_h}}}$ 。显然, ${B_h}$ 组微动特征参数中有 ${g_h}$ 组对应目标真实散射点,其余 ${B_h} - {g_h}$ 组对应虚假散射点,则虚假散射点与真实散射点的能量比 $\xi $ 定义如下:

      $ \xi {\rm{ = }}\frac{1}{H}\sum\limits_{h = 1}^H {\left( {\sum\limits_{b = {g_h}{\rm{ + }}1}^{{B_h}} {{A_b}}\Biggr/\sum\limits_{b = 1}^{{g_h}} {{A_b}} } \right)} $

      显然,E $\xi $ 越小,算法性能越好。图8给出了E $\xi $ 随信噪比的变化曲线,可以看出,当信噪比大于–10 dB时,虽然虚假散射点与真实散射点的能量比 $\xi $ 随着信噪比的降低而缓慢增加,但 $\xi $ 的取值始终较小,而特征参数提取精度E则始终保持不变,说明算法能够有效提取出目标真实散射的微动特征参数,具有良好的鲁棒性。当信噪比小于–10 dB时,E $\xi $ 快速增大,算法性能显著退化,不仅真实散射点的特征提取误差较大,而且由于虚假散射点与真实散射点的能量比 $\xi $ 显著增大,导致虚假散射点的出现,难以满足雷达目标微动特征提取需求。

      图  8  算法性能随信噪比变化曲线

      Figure 8.  The performance curve versus SNR

    • 本文针对现有多功能MIMO雷达的多目标微动特征提取效率不高的问题,提出了一种基于跟踪脉冲的MIMO雷达多目标微动特征提取方法。首先基于设计方向图与期望方向图逼近的思想,对MIMO雷达发射波形进行优化设计;在此基础上,建立跟踪脉冲发射时间序列优化模型并求解;最后直接利用窄带跟踪脉冲,采用OMP算法实现目标微动特征提取。该方法能够同时对不同方向的目标发射跟踪脉冲,从而同时提取多个目标的微动特征参数,相较于传统雷达对不同方向目标进行逐个微动特征提取的工作模式,提高了雷达系统的工作效率。本文可为提升MIMO雷达的空间目标探测能力提供参考。在后续的研究工作中,需进一步考虑空间目标微动形式的多样性,研究能够适用于各类微动目标的多目标特征提取方法。

参考文献 (18)

目录

    /

    返回文章
    返回