基于窄带雷达网的弹道目标三维进动特征提取

赵双 鲁卫红 冯存前 王义哲

引用本文:
Citation:

基于窄带雷达网的弹道目标三维进动特征提取

    作者简介: 赵双 (1992-), 男, 湖南湘谭人。现为空军工程大学防空反导学院硕士研究生。主要从事雷达信号处理研究。E-mail:zhaoshuanglw@163.com;鲁卫红 (1968-), 男, 河北安平人。现为空军工程大学副教授。主要从事雷达对抗、引战配合技术研究。E-mail:lwh_001@126.com;冯存前 (1976-), 男, 陕西富平人。现为空军工程大学教授, 主要从事自适应信号处理研究。E-mail:fengcunqian@sina.com;王义哲 (1992-), 男, 河南驻马店人。现为空军工程大学防空反导学院硕士研究生。主要从事雷达信号处理研究。E-mail:wangyizhe813@163.com.
    通讯作者: 赵双zhaoshuanglw@163.com
  • 基金项目:

    陕西省自然科学基础研究计划 2014JM8308

    国家自然科学基金 61372166

  • 中图分类号: TN957

Three-Dimensional Precession Feature Extraction of Ballistic Targets Based on Narrowband Radar Network

    Corresponding author:
  • Fund Project: The Project Supported by Natural Science Basic Research Plan in Shaanxi Province of China 2014JM8308The National Natural Science Foundation of China 61372166

    CLC number: TN957

  • 摘要: 微动特征是弹道目标识别的重要特征之一。针对单一视角观测难以获取目标真实微动参数的问题,该文提出了一种基于窄带雷达网的3维进动特征提取方法。首先建立了锥体目标进动模型,在考虑散射中心遮挡的前提下,详细分析了进动引发的微多普勒频率调制特性。然后基于锥顶微多普勒频率调制系数比,实现了不同视角下散射中心匹配关联,通过构建多视角联合方程组获取了目标的3维锥旋矢量,进而利用各散射中心微多普勒频率相关性,结合频率补偿的方法对锥体特征参数进行了提取,在此基础上解算出每一时刻锥顶坐标,实现了目标空间位置的3维重构。仿真结果证明了该方法的有效性与实用性。
  • 图 1  锥形弹头进动模型

    Figure 1.  The precession model of conical warhead

    图 2  组网雷达示意图

    Figure 2.  The sketch map of netted radar

    图 3  3部雷达时频分布图

    Figure 3.  The Time-frequency Distribution figure of three radars

    图 4  3部雷达IF提取结果

    Figure 4.  The IF result extracted by three radars

    图 5  补偿系数求解结果

    Figure 5.  The solving results of compensation coefficient

    图 6  锥顶散射中心的轨迹

    Figure 6.  The trajectory of the top scattering center

    图 7  鲁棒性分析

    Figure 7.  Robustness analysis

    表 1  锥体弹头进动及结构参数估计结果

    Table 1.  The estimation result of cone-shaped warhead's parameters

    参数 理论值 估计值 相对误差 (%)
    α1(°) 70.7288 70.7025 0.30
    α2(°) 36.8974 36.9851 0.037
    α3(°) 30.8829 30.7848 0.32
    θ(°) 13 13.5579 4.29
    r(m) 0.5 0.5237 4.74
    h1(m) 2.0 1.9132 4.43
    h2(m) 0.5 0.4811 3.20
    下载: 导出CSV
  • [1] Victoria S. American Missile Defense[M]. California:United States of America, 2010:44-78.
    [2] Chen V C. Advances in applications of radar micro-Doppler signatures[C]. 2014 IEEE Conference on Antenna Measurements & Application, Antibes, Juan-les-pins, 2014:1-4.
    [3] Wang Jun, Lei Peng, Sun Jin-ping, et al.. Spectral characteristics of mixed micro-Doppler time-frequency data sequences in micro-motion and inertial parameter estimation of radar targets[J]. IET Radar, Sonar & Navigation, 2014, 8(4):275-281.
    [4] Smith G E, Woodbridge K, Baker C J, et al.. Multistatic micro-Doppler radar signatures of personnel targets[J]. IET Signal Processing, 2010, 4(3):224-233. doi: 10.1049/iet-spr.2009.0058
    [5] 贺思三, 赵会宁, 冯存前.多视角距离像序列弹道目标的进动参数估计[J].信号处理, 2013, 29(8):1027-1035.He Si-san, Zhao Hui-ning, and Feng Cun-qian. Precession parameter estimation for ballistic targets based on multiaspect range-profile sequence[J]. Journal of Signal Processing, 2013, 29(8):1027-1035.
    [6] Liu L H, McLernon D, Ghogho M, et al.. Ballistic missile detection via micro-Doppler frequency estimation from radar return[J]. Digital Signal Processing, 2012, 22(1):87-95. doi: 10.1016/j.dsp.2011.10.009
    [7] 罗迎, 张群, 封同安, 等. OFD-LFM MIMO雷达中旋转目标微多普勒效应分析与三维微动特征提取[J].电子与信息学报, 2011, 33(1):8-13.Luo Ying, Zhang Qun, Feng Tong-an, et al.. Micro-Doppler effect analysis of rotating target and three-dimensional micro-motion feature extraction in OFD-LFM MIMO radar[J]. Journal of Electronics & Information Technology, 2011, 33(1):8-13.
    [8] 罗迎, 张群, 李松, 等.基于分布式组网雷达的弹道目标三维进动特征提取[J].电子学报, 2012, 40(6):1080-1085.Luo Ying, Zhang Qun, Li Song, et al.. Three-dimensional precession feature extraction of ballistic targets in distributed radar networks[J]. Acta Electronica Sinica, 2012, 40(6):1080-1085.
    [9] 雷腾, 刘进忙, 杨少春, 等.基于三站一维距离像融合的弹道目标特征提取方法研究[J].宇航学报, 2012, 33(2):228-234.Lei Teng, Liu Jin-mang, Yang Shao-chun, et al.. Study on feature extraction method of ballistic target based on threestation range profiles[J]. Journal of Astronautics, 2012, 33(2):228-234.
    [10] 韩勋, 杜兰, 刘宏伟.基于窄带雷达组网的空间锥体目标特征提取方法[J].电子与信息学报, 2014, 36(12):2956-2962.Han Xun, Du Lan, and Liu Hong-wei. Feature extraction of space cone-shaped target based on narrow-band radar networks[J]. Journal of Electronics & Information Technology, 2014, 36(12):2956-2962.
    [11] 宁超, 黄璟, 黄培康.基于HRRP的进动锥体目标特征参数求解方法[J].系统工程与电子技术, 2014, 36(4):650-655.Ning Chao, Huang Jing, and Huang Pei-kang. Solution for characteristic parameters of precession cone-shaped target using HRRP[J]. Systems Engineering and Electronics, 2014, 36(4):650-655.
    [12] Fang Hong-qing, Chen Long, and Raghavendiran Srinivasan. Influence of time and length size feature selections for human activity sequences recognition[J]. ISA Transactions, 2014, 53(1):134-140. doi: 10.1016/j.isatra.2013.09.001
    [13] Li Po, Wang De-chun, and Wang Lu. Separation of microDoppler signals based on time frequency filter and Viterbi algorithm[J]. Signal, Image and Video Processing, 2013, 7(3):593-605. doi: 10.1007/s11760-011-0263-3
    [14] 关永胜, 左群声, 刘宏伟.高噪声环境下微动多目标分辨[J].电子与信息学报, 2010, 32(11):2630-2635.Guan Yong-sheng, Zuo Qun-sheng, and Liu Hong-wei. Micro-motion targets resolution in a high noise environment[J]. Journal of Electronics & Information Technology, 2010, 32(11):2630-2635.
    [15] 王辉, 黄成军, 姚林朋, 等.重排Cohen类时频分布用于GIS局部放电声信号时频分析[J].高电压技术, 2010, 36(11):2724-2730.Wang Hui, Huang Chen-jun, Yao Lin-peng, et al.. Application of reassigned cohen class time-frequency distribution to the analysis of acoustic emission partial discharge signal for GIS[J]. High Voltage Engineering, 2010, 36(11):2724-2730.
    [16] 葛哲生, 陈仲生. Matlab时频分析技术及其应用[M].北京:人民邮电出版社, 2006:107-112.Ge Zhe-sheng and Chen Zhong-sheng. Time-Frequency Analysis Techniques of Matlab and Its Application[M]. Beijing:Posts & Telecom Press, 2006:107-112.
  • [1] 杨利民苏卫民顾红薄超 . 基于二维映射的LFM UWB 雷达动目标速度估计. 雷达学报, 2012, 1(3): 232-237. doi: 10.3724/SP.J.1300.2012.20048
    [2] 许京伟朱圣棋廖桂生张玉洪 . 频率分集阵雷达技术探讨. 雷达学报, 2018, 7(2): 167-182. doi: 10.12000/JR18023
    [3] 周超伟李真芳王跃锟解金卫 . 联合多方位角调频率估计的星载SAR三维成像方法. 雷达学报, 2018, 7(6): 696-704. doi: 10.12000/JR18094
    [4] 曾操梁思嘉王威徐青 . 基于频率步进信号的旋转式合成孔径雷达成像方法. 雷达学报, 2014, 3(4): 401-408. doi: 10.3724/SP.J.1300.2014.14043
    [5] 刘晓明文贡坚钟金荣 . 基于SAR 数据的三维散射中心模型位置重构方法. 雷达学报, 2013, 2(2): 187-194. doi: 10.3724/SP.J.1300.2013.20080
    [6] 陈曾平张炜承林钱强 . 宽带雷达ISAR 成像相位补偿新方法(英文). 雷达学报, 2013, 2(1): 23-29. doi: 10.3724/SP.J.1300.2013.13023
    [7] 丁振宇谭维贤王彦平洪文吴一戎 . 基于波数域子孔径的机载三维SAR偏航角运动误差补偿. 雷达学报, 2015, 4(4): 467-473. doi: 10.12000/JR15016
    [8] 田鹤李道京 . 稀疏重航过阵列SAR运动误差补偿和三维成像方法. 雷达学报, 2018, 7(6): 717-729. doi: 10.12000/JR18101
    [9] 赵婉婉王鹏波门志荣李春升 . 一种基于二维信号稀疏重构的互质采样星载SAR成像处理方法. 雷达学报, 2019, 8(): 1-12. doi: 10.12000/JR19086
    [10] 饶云华明燕珍林静朱逢园万显荣龚子平 . WiFi外辐射源雷达参考信号重构及其对探测性能影响研究. 雷达学报, 2016, 5(3): 284-292. doi: 10.12000/JR15108
    [11] 刘奇勇张群洪文苏令华梁佳 . 基于参数估计的下视稀疏阵列三维SAR运动误差补偿和成像处理方法. 雷达学报, 2018, 7(6): 730-739. doi: 10.12000/JR18107
    [12] 赵雨露张群英李超纪奕才方广有 . 视频合成孔径雷达振动误差分析及补偿方案研究. 雷达学报, 2015, 4(2): 230-239. doi: 10.12000/JR14153
    [13] 洪永彬张勇鲁振兴黄巍 . 一种高效的基于对比度的步进频雷达运动补偿算法. 雷达学报, 2016, 5(4): 378-388. doi: 10.12000/JR16068
    [14] 赵军田斌朱岱寅 . 基于PAST处理的机载双基雷达自适应角度-多普勒补偿算法. 雷达学报, 2017, 6(6): 594-601. doi: 10.12000/JR17053
    [15] 关欣仲利华胡东辉丁赤飚 . 一种基于RSPWVD-Hough 变换的无源雷达多普勒展宽补偿方法. 雷达学报, 2013, 2(4): 430-438. doi: 10.3724/SP.J.1300.2013.13073
    [16] 胡程邓云开田卫明曾涛 . 地基干涉合成孔径雷达图像非线性大气相位补偿方法. 雷达学报, 2019, 8(6): 831-840. doi: 10.12000/JR19073
    [17] 赵耀东吕晓德李纪传向茂生 . 无源雷达多普勒谱分析实现动目标检测的方法. 雷达学报, 2013, 2(2): 247-256. doi: 10.3724/SP.J.1300.2012.20081
    [18] 陈小龙关键何友于晓涵 . 高分辨稀疏表示及其在雷达动目标检测中的应用. 雷达学报, 2017, 6(3): 239-251. doi: 10.12000/JR16110
    [19] 丁永超刘成城赵拥军刘亚奇 . 基于反对角单位阵的窄带多径信号接收波束形成算法. 雷达学报, 2015, 4(2): 192-198. doi: 10.12000/JR14012
    [20] 郑明洁闫贺张冰尘赵凤军杨汝良 . 基于DBS 图像的双通道广域监视雷达动目标检测和参数估计方法. 雷达学报, 2012, 1(1): 36-42. doi: 10.3724/SP.J.1300.2013.20007
  • 加载中
图(7)表(1)
计量
  • 文章访问数:  844
  • HTML浏览量:  138
  • PDF下载量:  2006
  • 被引次数: 0
出版历程
  • 收稿日期:  2015-12-23
  • 录用日期:  2016-03-22
  • 网络出版日期:  2016-05-09
  • 刊出日期:  2017-02-28

基于窄带雷达网的弹道目标三维进动特征提取

    通讯作者: 赵双zhaoshuanglw@163.com
    作者简介: 赵双 (1992-), 男, 湖南湘谭人。现为空军工程大学防空反导学院硕士研究生。主要从事雷达信号处理研究。E-mail:zhaoshuanglw@163.com;鲁卫红 (1968-), 男, 河北安平人。现为空军工程大学副教授。主要从事雷达对抗、引战配合技术研究。E-mail:lwh_001@126.com;冯存前 (1976-), 男, 陕西富平人。现为空军工程大学教授, 主要从事自适应信号处理研究。E-mail:fengcunqian@sina.com;王义哲 (1992-), 男, 河南驻马店人。现为空军工程大学防空反导学院硕士研究生。主要从事雷达信号处理研究。E-mail:wangyizhe813@163.com
  • ①. 空军工程大学防空反导学院 西安 710051
  • ②. 信息感知技术协同创新中心 西安 710077
基金项目:  陕西省自然科学基础研究计划 2014JM8308国家自然科学基金 61372166

摘要: 微动特征是弹道目标识别的重要特征之一。针对单一视角观测难以获取目标真实微动参数的问题,该文提出了一种基于窄带雷达网的3维进动特征提取方法。首先建立了锥体目标进动模型,在考虑散射中心遮挡的前提下,详细分析了进动引发的微多普勒频率调制特性。然后基于锥顶微多普勒频率调制系数比,实现了不同视角下散射中心匹配关联,通过构建多视角联合方程组获取了目标的3维锥旋矢量,进而利用各散射中心微多普勒频率相关性,结合频率补偿的方法对锥体特征参数进行了提取,在此基础上解算出每一时刻锥顶坐标,实现了目标空间位置的3维重构。仿真结果证明了该方法的有效性与实用性。

English Abstract

    • 随着诱饵技术和电磁控制技术等的不断发展,基于传统特征量的弹道目标识别技术已难以满足未来高科技战争的需求,而微动特征作为弹道目标的固有属性,难以被模仿,且弹头和诱饵、碎片等运动形式存在明显差异,因而被用来区分识别真弹头,近年来受到国内外研究学者的广泛关注[1–3]。由于目标微动会对雷达回波产生微多普勒调制,无论是宽带雷达还是窄带雷达,均可通过对回波信号参数进行估计来提取目标的微动特征,进而进行目标的分类识别。然而研究表明[4–6],单基雷达获取的目标微多普勒信息均具有较强的姿态敏感性,不同视角获得的微动特征不同,因而难以准确反映目标的真实空间结构。考虑到组网雷达能获得目标在不同视角上的观测信息,于是有望利用多部雷达联合观测对目标进行3维微动特征提取,从而克服单一视角的局限性,提高目标识别精度。

      考虑宽带雷达能得到目标的高分辨距离像,获取目标更细微的形状结构特征,目前基于宽带雷达组网的目标微动特征提取研究较多。文献[7]对MIMO雷达中的旋转目标微多普勒效应展开分析,并基于时间-距离像对目标的3维微动特征进行了提取;文献[8]利用分布式组网雷达对有翼锥体目标进行联合观测,并基于不同视角微多普勒特征的相关性,借助几何分析的方法,实现了目标的3维进动特征提取;文献[9]通过三站1维距离像融合实现了锥体目标的3维重构。由此可见,基于宽带雷达组网的目标微动特征提取技术研究已较为成熟,然而,就目前现有雷达实际装备情况来看,由于宽带高分辨雷达价格昂贵,一时间仍难以实现对窄带雷达的全面升级和替换,雷达网仍以窄带体制为主,因此研究利用窄带雷达网对弹道目标进行特征提取具有更大的实际意义,能有效解决目前理论研究与实际脱节的矛盾;文献[10]在窄带组网体制下推导了不同视角锥体散射中心瞬时频率变化关系,利用频谱熵实现了散射中心的匹配关联,并进一步提出了基于散射中心瞬时频率相关性的目标参数提取算法,但没能获取目标的3维微动特征,也无法实现空间目标的3维重构。

      本文在文献[10]的基础上,进一步对窄带雷达网中的锥体目标3维进动特征提取展开了研究。在详细分析了锥体进动引发的微多普勒频率调制特性的基础上,利用锥顶微多普勒频率调制系数比,实现了不同视角下散射中心匹配关联,并获取了目标的3维锥旋矢量,进而利用锥顶和底面边缘散射中心微多普勒频率相关性,结合频率补偿的方法对锥体特征参数进行了提取,在此基础上解算出每一时刻锥顶坐标,从而实现了目标空间位置的3维重构。最后,仿真分析了本文方法的有效性。

    • 以无翼锥形弹头为例,建立进动模型如图 1所示,由于其存在旋转对称性,因此仅考虑其做锥旋运动。假定目标以角速度ωc绕锥旋轴旋转,且锥旋轴与目标对称轴相交于点o,两者夹角为θ,以交点为坐标原点建立参考坐标系oxyz如下,目标对称轴初始方位为φ0, LOS表示雷达视线方向,其在参考坐标系中的方位角和俯仰角分别为 (ε, χ),与锥旋轴的夹角为α,与对称轴夹角为β,定义雷达视线方向与对称轴构成的平面为底面圆环的电磁波入射平面,该平面与圆环交于pq两点。锥体高度为h,底面半径为r,锥顶与进动中心的距离为h1,底面中心与进动中心的距离为h2,且目标满足远场条件,雷达与进动中心的距离为R0

      图  1  锥形弹头进动模型

      Figure 1.  The precession model of conical warhead

      根据散射中心理论,对于旋转对称目标,其高频散射特性主要由锥顶D及底面边缘两个散射中心p, q确定[1, 5]。由几何关系分析可得雷达视线与对称轴的夹角β满足:

      式中,φ为初始相位角,且由文献[6]可知是$\varphi ={{\phi }_{0}}-\varepsilon $。进一步对各散射中心在雷达视线上的投影关系分析可知,锥体3个散射中心到雷达的距离分别为:

      考虑到锥体目标在实际运动中各部分之间存在相互遮挡,目标上各散射中心不能始终保持同时可见,使得式 (2) 的使用范围受到限制。但锥顶D和近散射点p在大部分情况下都能被观测到[3],并能够获得二者的稳定连续观测信息,因此本文主要利用D, p的微动信息展开后继研究。假设雷达波长为λ,由式 (1) 和式 (2) 可得,进动引发的D, p两点的微多普勒调制为[11]

      由式 (3) 可以看出,D点的微多普勒频率变化服从正弦规律,而p点的微多普勒频率由两部分之和组成,不再服从简单的正弦调制规律,且两点微多普勒频率均与目标的进动和结构特征有关,共包含wc, h1, h2, r, θ, α, φ 7个未知参数,其中ωc, φ可通过提取正弦曲线特征得到,而θ, α两者之间存在耦合,仅通过单部雷达,仍无法实现对目标进动角及尺寸大小的求解。考虑到多视角观测能获得更加丰富的目标信息,具有较好的解耦合性能,因此本文将采用雷达组网方式对目标特征进行提取,并进一步实现3维重构。

    • 首先建立窄带雷达网系统观测模型如图 2所示,图中OXYZ为全局坐标系,与参考坐标系oxyz平行,假定系统中共有N部窄带雷达同时进行观测,并都已满足时空同步要求,各雷达视线在OXYZ坐标系中的方位角和俯仰角为 (εi, χi), ni (i=1, 2, ···, N) 为雷达视线方向,满足:

      图  2  组网雷达示意图

      Figure 2.  The sketch map of netted radar

      由第2节分析可知,当采用多部雷达同时进行观测时,由于各雷达观测视角不同,同一时刻目标各散射中心在雷达视线上的投影位置排列顺序将存在差异,相对应地,同一时刻各散射中心的微多普勒频率也会不同。因此,在利用组网雷达进行特征提取之前,首先得实现不同视角散射中心的匹配关联。

      由式 (3) 可以看出,对于同一观测目标而言,锥顶D的微多普勒频率调制系数A仅与雷达观测视角有关,任取雷达网中两部雷达,其调制系数比满足:

      p点调制规律更为复杂,不具备上述比例关系,因此,通过比较观察不同雷达间的调制系数比即可实现散射中心的匹配关联。文献[10]采用频谱分析的方法,通过计算不同散射中心的频谱熵来实现散射中心的匹配关联,然而在两个散射中心回波信号无法分离的情况下,散射中心无法与各自频谱一一对应起来,因此该方法存在较大的局限性。相比而言,本文方法则更加简单实用。

    • 为更好地实现对锥体目标的3维重构,首先对锥旋轴方向进行估计。采用Viterbi算法提取锥顶微多普勒曲线振幅得到:

      Viterbi算法[12–14]作为信号隐状态估计的有效手段之一,能够依据各信号成分强度对信号进行逐次分离,因此常被用来对多目标信号瞬时频率进行估计。与逆Randon变换、Hough变换等曲线参数提取方法相比,Viterbi算法对各信号分量形式依赖性不高,即使是非正弦信号,同样能实现瞬时频率的准确估计。

      此外,雷达观测视角αi还满足:

      式 (7) 中,$\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{\omega }\text{ = }{{\left( {{\omega }_{x}}, {{\omega }_{y}}, {{\omega }_{z}} \right)}^{\text{T}}}$,联立式 (4)、式 (6)、式 (7),令B=h1sinθ,此时方程组中共包含B, ωx, ωy, ωz 4个未知参数,因此,至少需要3部雷达同时进行观测才能实现对上述参数的求解。进一步将求得的参数回代到方程组,还可以确定sinai的大小。文献[10]采用基于视线角方差最小准则的频率搜索补偿方法对雷达视线角进行估计,需要经过多次循环迭代才能得到准确估计值,计算复杂,且易受噪声影响,而本文所提方法用到的只是锥顶微多普勒频率曲线的振幅和周期,且这两个曲线参数均可由Viterbi算法准确提取得到,在求得锥旋矢量的同时也能估计出各雷达视线角的大小,计算更为简单,算法稳定性更好。

    • 在上述分析的基础上,若要提取锥体弹头参数,还需p对点的微多普勒频率进行充分利用,观察式 (3) 可知,fdp由正弦部分和非正弦部分之和组成,且正弦部分${{f}_{k}}=-2{{\omega }_{\text{c}}}{{h}_{2}}\sin \theta \sin \alpha \cos ({{\omega }_{\text{c}}}t+\varphi )\text{ }/\text{ }\!\!\lambda\!\!\text{ }$满足fk=(-h2/h1)fd-D,而非正弦部分此时仅包含 (r, θ) 两个未知参数。若能将正弦部分完全补偿,便可利用多视角观测对 (r, θ) 联立求解。考虑到锥体目标尺寸信息仍然未知,先假设补偿系数为η,且η0=h2/h1,当η=η0时,便可实现完全补偿,于是p点补偿后的微多普勒频率满足:

      η进行遍历,利用补偿后的频率两两联立方程可求得:

      对每个η取值所对应求得的所有结果$\left( {{{\hat{r}}}_{\xi }}\left| \eta, {{{\hat{\theta }}}_{\xi }}\left| \eta \right. \right. \right)$做进一步处理,并定义归一化标准差σ为:

      式中,$\Delta r=\left[{{{\hat{r}}}_{1}}-\bar{r}\ \ {{{\hat{r}}}_{2}}-\bar{r}\cdots {{{\hat{r}}}_{C_{N}^{2}}}-\bar{r} \right], \Delta \theta =\left[{{{\hat{\theta }}}_{1}}-\bar{\theta }\ \ {{{\hat{\theta }}}_{2}}-\bar{\theta }\cdots {{{\hat{\theta }}}_{C_{N}^{2}}}-\bar{\theta } \right], \bar{r}, \bar{\theta }$为平均值,按照上述归一化标准差定义,对于每一个η取值均能得到对应的σ。若η=η1时,σ取得最小值,则说明此时fk被补偿得最完全,补偿系数η1也越接近η0,由此可求得:

      结合4.1节分析,将${\hat{\theta }} $代入B=h1sinθ中,于是求得${{h}_{1}}=B\sin \hat{\theta }, {{h}_{2}}={{h}_{1}}\cdot {{\eta }_{1}}$。

    • 在求得锥体目标结构参数及旋转轴方向的基础上,为实现对目标空间位置的3维重构,还需确定各散射中心的相对位置,由于底面边缘两个散射中心会随雷达视线方向改变产生滑动,位置坐标不易确定,因而本文从锥顶散射中心入手,在锥体目标结构参数已知的条件下,只要能够求得每一时刻锥顶坐标,同样能实现对目标空间位置的3维重构。由于窄带雷达距离分辨力较低,难以直接从目标回波中获得各散射中心的径向距离变化规律,因此本文考虑在已知各散射中心运动形式和参数基础上,通过微多普勒频率反推每一时刻各散射中心相对应的径向距离变化。

      d=RDR0,由式 (2) 可得:

      式中,${{\varphi }_{i}}={{\phi }_{0}}-{{\chi }_{i}}$,由于每一时刻D点的微多普勒频率fdD均已获得,且fdDRD满足导数关系,因而di在每一时刻的值也能求解得到。此外,结合图 1可知,dioD在第i部雷达视线上的投影,同时还应满足:

      若令oD=(Dx, Dy, Dz)T,通过3部雷达同时进行观测可以解算出oD为:

      $$ \rm{\pmb{OD}} = {\left[\begin{gathered} \mathit{\boldsymbol{n}}_1^{\rm{T}} \hfill \\ \mathit{\boldsymbol{n}}_2^{\rm{T}} \hfill \\ \mathit{\boldsymbol{n}}_3^{\rm{T}} \hfill \\ \end{gathered} \right]^{ - 1}} \cdot \left[\begin{gathered} {d_1} \hfill \\ {d_2} \hfill \\ {d_3} \hfill \\ \end{gathered} \right] $$

      综上所述,基于窄带雷达组网的弹道目标3维微动特征提取及重构步骤为:

      步骤1建立弹道目标进动模型,分析各散射中心微多普勒调制规律;

      步骤2对目标回波进行时频分析,采用Viterbi算法提取各散射中心微多普勒曲线;

      步骤3基于锥顶微多普勒频率调制系数比,实现不同视角下散射中心匹配关联;

      步骤4提取3维锥旋矢量和目标结构参数,在此基础上解算出每一时刻锥顶坐标,从而实现锥体目标空间3维重构。

    • 在下述仿真中设定目标为锥体,目标参数设置为:h1=2.0 m, h2=0.5 m, r=0.5 m, h=2.5 m, θ=13°,目标对称轴初始方位角φ0=60°,目标的锥旋频率为fc=4 Hz,锥旋矢量为$\left( 2\sqrt{3}\text{ }\!\!\pi\!\!\text{, 4}\sqrt{3}\text{ }\!\!\pi\!\!\text{ }, 2\text{ }\!\!\pi\!\!\text{ } \right) $。雷达参数设置为:载频f=8×109 Hz,信号带宽为5 MHz,雷达脉冲重复频率为2000 Hz,积累时间1 s,信噪比为10 dB。在全局坐标系中3部雷达M1, M2, M3测得的目标方位角和俯仰角 (εi, χi) 分别为 (40°, 84°), (45°, 48°), (30°, 17°)。图 3分别为该3部雷达获得的同一时间段内目标回波Cohen类时频分布重排结果,可以看出重排后的谱图不仅具有更好的时频聚集性,同时还有效抑制了各分量之间的交叉项[15, 16],有利于提高各散射中心瞬时频率的估计精度。

      图  3  3部雷达时频分布图

      Figure 3.  The Time-frequency Distribution figure of three radars

      在上述谱图重排的基础上,进一步采用Viterbi算法对各散射中心瞬时频率进行提取,并通过拟合更好地削减了交叉项所带来的不利影响,得到各雷达散射中心瞬时频率估计结果如图 4所示。然后依据第3节匹配关联准则,对各散射中心曲线幅度做进一步处理,区分出锥顶和锥底边缘散射中心。此时提取到各雷达锥顶正弦曲线的频率为$\left\| {\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{\omega }} \right\|$ =25.14 rad/s,振幅分别为567.6 Hz, 361.8 Hz, 307.8 Hz,而由式 (6) 计算得到的振幅理论值分别为569.2636 Hz, 362.0643 Hz, 309.5390 Hz,两者相当接近,代入方程组式 (7) 可以求解得到3维锥旋矢量$\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{\omega }={{\left( 10.8598, 21.7937, 6.2545 \right)}^{\text{T}}}$,与理论值基本吻合,同时可以求得各雷达视线角分别为${{{\hat{\alpha }}}_{1}}=70.7025, {{{\hat{\alpha }}}_{2}}=36.9851, {{{\hat{\alpha }}}_{3}}=30.7848$。

      图  4  3部雷达IF提取结果

      Figure 4.  The IF result extracted by three radars

      依据4.2节中锥体弹头参数提取算法,由求得的${{{\hat{\alpha }}}_{i}}$和各散射中心瞬时频率对 (r, θ) 进行估计,可以得到遍历时η归一化标准差σ随之变化的结果如图 5所示,当η=0.253时,σ取得最小值,接近于理论分析值η0=0.25,进一步将该η值代到方程式 (11),从而求得参数r, θ, h1, h2估计值。定义相对误差=|理论值–估计值|/理论值,各参数估计结果如表 1所示,其中α1, α2, α3估计精度相对较高,r, θ, h1, h2由于受到各散射中心瞬时频率提取误差的影响,估计精度则相对较低,通过进一步提高时频分辨率或者增加雷达网中雷达观测数量可实现上述参数估计精度的提升。但总的来说,各参数估计相对误差均小于5%,满足目标识别的精度要求,可用于下一步对锥顶坐标的求解。而在相同仿真条件下,文献[10]中的目标参数平均估计相对误差则接近于15%,明显高于本文误差,这也从侧面反映出本文参数提取算法的稳定性。

      图  5  补偿系数求解结果

      Figure 5.  The solving results of compensation coefficient

      参数 理论值 估计值 相对误差 (%)
      α1(°) 70.7288 70.7025 0.30
      α2(°) 36.8974 36.9851 0.037
      α3(°) 30.8829 30.7848 0.32
      θ(°) 13 13.5579 4.29
      r(m) 0.5 0.5237 4.74
      h1(m) 2.0 1.9132 4.43
      h2(m) 0.5 0.4811 3.20

      表 1  锥体弹头进动及结构参数估计结果

      Table 1.  The estimation result of cone-shaped warhead's parameters

      在上述目标特征参数提取的基础上,进一步按照4.3节所提算法对锥顶坐标进行求解。当时t=0.25 s,求得锥顶坐标$\left( {{{\hat{D}}}_{x}}, {{{\hat{D}}}_{y}}, {{{\hat{D}}}_{z}} \right)=\left( 1.1739, 1.1888, 0.5898 \right) $,与理论值基本相符,并最终得到在观测时间内0~0.25 s锥顶散射中心的实际轨迹如图 6所示,与其理论轨迹近乎重合,从而更加充分地说明了本文重构方法的准确性和有效性。

      图  6  锥顶散射中心的轨迹

      Figure 6.  The trajectory of the top scattering center

    • 为了充分验证本文所提算法的鲁棒性,仿真分析了曲线参数估计误差对目标参数提取及重构精度的影响。由于在4.1节中3维锥旋矢量的求解精度主要受到锥顶微多普勒幅度和微动周期的影响,而微动周期通常都能被准确估计,因而在此主要分析微多普勒幅度的影响。在4.2节中,目标参数的准确提取关键在于式 (9) 的求解,且求解精度主要受瞬时频率提取误差制约,因此也有必要对瞬时频率提取误差带来的影响进行分析。为便于分析,定义归一化误差如下:

      上式中为估计值,X为真实值,进一步定义|η|为归一化绝对误差。假设锥顶微多普勒曲线幅度AD提取值以及各散射中心瞬时频率的提取值归一化误差服从[–a a]上的均匀分布,采用蒙特卡洛方法进行分析,仿真100次,可以得到目标特征参数归一化绝对误差平均值的变化如图 7所示。从图 7(a)中可以看出,当a在区间[0, 0.1]变化时,锥旋矢量各方向分量估计误差均呈线性增加趋势,但总的来说,各方向分量估计精度仍然较高。而在图 7(b)中,当a在区间[0, 0.03]变化时,进动角θ变化较为平和,h2/h1则几乎不受影响,只有底面半径r值的估计精度变化最为敏感,随着各种高性能时频分析工具的发展,完全可以将瞬时频率提取精度进一步提高,从而满足目标参数的高精度提取。

      图  7  鲁棒性分析

      Figure 7.  Robustness analysis

      综上所述可知,本文所提锥体目标参数提取及重构算法在一定程度上受到微多普勒曲线参数提取精度的影响,但由于文中所采用的基于Cohen类时频重排的Viterbi算法较好地实现瞬时频率曲线的提取,因此本文算法能保证目标参数提取及重构的可靠性,可用于目标识别。

    • 本文对基于窄带雷达网的锥体目标3维进动特征提取问题展开了研究。依据目标的多视角微多普勒频率调制特性,利用3部雷达获取了目标的3维锥旋矢量及特征参数,并进一步解算出每一时刻锥顶坐标,实现了目标空间位置的3维重构。仿真结果表明,本文所提方法目标参数估计精度高,重构性能好,能够有效克服目标散射中心遮挡和姿态敏感性的不利影响,为基于窄带雷达的空间目标准确识别提供了解决方案。考虑到在实际的导弹防御体系当中,可能会出现窄带雷达和宽带雷达同时对目标进行观测识别的情况,后继工作将就不同体制雷达对目标3维特征参数的融合提取问题展开研究。

参考文献 (16)

目录

    /

    返回文章
    返回