交错阵列甚低频信号产生方法研究

崔岸婧 李道京 周凯 王宇 洪峻

崔岸婧, 李道京, 周凯, 等. 交错阵列甚低频信号产生方法研究[J]. 雷达学报, 2020, 9(5): 925–938. doi:  10.12000/JR20082
引用本文: 崔岸婧, 李道京, 周凯, 等. 交错阵列甚低频信号产生方法研究[J]. 雷达学报, 2020, 9(5): 925–938. doi:  10.12000/JR20082
CUI Anjing, LI Daojing, ZHOU Kai, et al. Research on the method of composing very low frequency signals based on the staggered array[J]. Journal of Radars, 2020, 9(5): 925–938. doi:  10.12000/JR20082
Citation: CUI Anjing, LI Daojing, ZHOU Kai, et al. Research on the method of composing very low frequency signals based on the staggered array [J]. Journal of Radars, 2020, 9(5): 925–938. doi:  10.12000/JR20082

交错阵列甚低频信号产生方法研究

doi: 10.12000/JR20082
基金项目: 国家部委基金,中国科学院空天信息创新研究院课题(Y910340Z2F)
详细信息
    作者简介:

    崔岸婧(1997–),女,中国科学院大学博士生。研究方向为阵列天线与信号处理。E-mail: ajcui@qq.com

    李道京(1964–),男,中国科学院空天信息创新研究院,研究员,博士生导师,研究方向为雷达系统与雷达信号处理。E-mail: lidj@aircas.ac.cn

    周 凯(1995–),男,中国科学院大学博士生,研究方向为信号与信息处理。E-mail: zk_6810@163.com

    王 宇(1976–),男,中国科学院空天信息创新研究院,研究员,研究方向为机/星载SAR定标,定量遥感应用。E-mail: wangyu@mail.ie.ac.cn

    洪 峻(1960–),男,中国科学院空天信息创新研究院,研究员,博士生导师,研究方向为SAR定标技术、定量化SAR图像应用等。E-mail: jhong@mail.ie.ac.cn

    通讯作者:

    李道京 lidj@aircas.ac.cn

  • 责任主编:廖桂生 Corresponding Editor: LIAO Guisheng
  • 中图分类号: TN951

Research on the Method of Composing Very Low Frequency Signals Based on the Staggered Array

Funds: The National Ministries Foundation, The Aerospace Information Research Institute, Chinese Academy of Science (Y910340Z2F)
More Information
  • 摘要: 10 kHz量级甚低频电磁波信号具有较强的地物穿透能力,可用于地质勘探。由于其天线尺寸在10 km量级,其应用场合受到限制,研究基于适当尺寸高频雷达天线的甚低频电磁波信号产生方法具有重要意义。该文提出基于高频阵列天线产生甚低频信号的概念,利用阵列天线合成产生近光速远离运动雷达多普勒信号,实现信号频率的大幅降低。给出了发射波形、交错阵列设计和阵列参数选择方法。将周期脉冲串信号作为辐射单元信号,增大合成信号脉宽。利用阵列产生的脉宽展宽量填补脉冲信号的休止期,在目标区合成时间连续的甚低频信号。采用峰值旁瓣比(PSLR)、积分旁瓣比(ISLR)、阵列发射信号与合成信号的频谱对比评价合成低频信号的性能和发射信号的能量利用率。该文仿真了百米量级阵列100 MHz辐射单元信号在目标区合成10 kHz甚低频信号的情况:9行阵列构成交错阵列、辐射单元信号脉宽设置为0.115 μs时,合成信号频谱的峰值旁瓣比和积分旁瓣比分别为–13.34 dB和–9.44 dB, 10 kHz低频信号在合成信号中的能量占比为89.79%。该文分析了辐射单元间距误差、辐射单元信号时间、相位与幅度误差以及目标偏离预定位置的影响。仿真结果表明了该文方法的有效性。
  • 图  1  雷达运动示意图

    Figure  1.  Radar motion diagram

    图  2  多普勒效应中的发射/接收信号波形与频谱

    Figure  2.  The emission/received signal waveform and spectrum of Doppler effect

    图  3  运动雷达发射信号过程分解和天线阵列等效雷达运动示意图

    Figure  3.  Schematic diagram of the decomposed process of radar transmitting signals and the equivalent radar motion generated by the antenna array

    图  4  交错阵列结构

    Figure  4.  Staggered array structure

    图  5  辐射单元信号示意图

    Figure  5.  Schematic diagram of the radiation unit signal

    图  6  辐射单元信号周期为${{{T_0}} / 2}$时合成信号的波形、频谱与频谱对比图

    Figure  6.  Waveform, spectrum and spectral contrast diagram of the composite signal when the radiation unit signal period being ${{{T_0}} / 2}$

    图  7  辐射单元信号周期为${{{T_0}} / 6}$时合成信号的波形、频谱与频谱对比图

    Figure  7.  Waveform, spectrum and spectral contrast diagram of the composite signal when the radiation unit signal period being ${{{T_0}} / 6}$

    图  8  辐射单元信号周期为${{{T_0}} / 6}$时9行阵列构成的交错阵列在目标区合成信号的波形、频谱与频谱对比图

    Figure  8.  Waveform, spectrum and spectral contrast diagram of the composite signal generated by the 9-row staggered array when the radiation unit signal period being ${{{T_0}} / 6}$

    图  9  辐射单元信号包络示意图

    Figure  9.  Schematic diagram of the envelope of radiation unit signals

    图  10  辐射单元间距误差分布、目标区合成信号的波形与频谱图

    Figure  10.  The distribution of spacing error, waveform and spectrum of the composite signal in target area

    图  11  辐射单元信号时间误差分布、目标区合成信号的波形与频谱图

    Figure  11.  The distribution of radiation unit signal time error, waveform and spectrum of the composite signal in target area

    图  12  辐射单元信号相位误差分布、目标区合成信号的波形与频谱图

    Figure  12.  The distribution of radiation unit signal phase error, waveform and spectrum of the composite signal in target area

    图  13  辐射单元信号幅度误差分布、目标区合成信号的波形与频谱图

    Figure  13.  The distribution of radiation unit signal amplitude error, waveform and spectrum of the composite signal in target area

    图  14  目标在xoy平面内45°扫描角方向距离阵列近端20 km时目标区合成信号的波形与频谱

    Figure  14.  Waveform and spectrum of the composite signal in the target area when the target is 20 km away from the near end of the array at 45° scanning angle in xoy plane

    图  15  目标在xoy平面内45°扫描角方向距离阵列近端40 km时目标区合成信号的波形与频谱

    Figure  15.  Waveform and spectrum of the composite signal in the target area when the target is 40 km away from the near end of the array at 45° scanning angle in xoy plane

    图  16  目标沿$z$轴偏离1 km时目标区合成信号的波形与频谱

    Figure  16.  Waveform and spectrum of the composite signal in the target area when the target deviates 1 km along the z-axis

    图  17  受综合误差影响时目标区合成信号的波形与频谱

    Figure  17.  Waveform and spectrum of the composite signal in the target area influenced by synthetical error

    表  1  辐射单元信号周期取${{{{{T_0}}}} / {{2}}}$${{{{{T_0}}}} / {{6}}}$时的仿真参数与结果

    Table  1.   Simulation parameters and results when the radiation unit signal period being ${{{{{T_0}}}} / {{2}}}$ and ${{{{{T_0}} / 6}}}$ respectively

    参数辐射单元信号周期取${{{T_0}} / 2}$辐射单元信号周期取${{{T_0}} / 6}$
    最大辐射单元信号周期${T_0}$(μs)1.381.38
    辐射单元信号周期数(个)200600
    辐射单元信号脉宽(μs)0.3450.115
    辐射单元信号脉宽展宽(μs)0.3450.115
    辐射单元信号周期(μs)0.690.23
    辐射单元信号休止期(μs)0.3450.115
    合成信号脉宽(μs)138.36138.59
    峰值旁瓣比(dB)–11.45–13.34
    积分旁瓣比(dB)–4.33–8.77
    低频信号能量占比(%)73.0588.29
    频谱对比图中的10 kHz分量(dB)–14.610–5.081
    发射信号能量利用率(%)18.6055.71
    下载: 导出CSV

    表  2  不同误差影响下合成信号频谱参数(dB)

    Table  2.   Spectrum parameters of the composite signal influenced by different errors (dB)

    误差峰值旁瓣比积分旁瓣比
    辐射单元间距误差–13.31–8.95
    辐射单元信号时间误差–13.34–7.64
    辐射单元信号相位误差–13.34–9.07
    辐射单元信号幅度误差–13.34–9.33
    目标在xoy平面内距离阵列近端20 km–13.34–9.44
    目标在xoy平面内距离阵列近端40 km–13.34–9.44
    目标在z轴偏离1 km–13.34–8.97
    综合误差–13.28–6.63
    下载: 导出CSV
  • [1] 汪谋. 地质雷达探测效果影响因素研究[J]. 雷达科学与技术, 2007, 5(2): 86–90.WANG Mou. Research on influencing factors of GPR’s detection effectiveness[J]. Radar Science and Technology, 2007, 5(2): 86–90.
    [2] 卓贤军, 陆建勋, 赵国泽, 等. 极低频探地(WEM)工程[J]. 中国工程科学, 2011, 13(9): 42–50.ZHUO Xianjun, LU Jianxun, ZHAO Guoze, et al. The extremely low frequency engineering project using WEM for underground exploration[J]. Strategic Study of CAE, 2011, 13(9): 42–50.
    [3] 卓贤军, 陆建勋. “极低频探地工程”在资源探测和地震预测中的应用与展望[J]. 舰船科学技术, 2010, 32(6): 3–7, 30.ZHUO Xianjun and LU Jianxun. Application and prospect of WEM to resource exploration and earthquake predication[J]. Ship Science and Technology, 2010, 32(6): 3–7, 30.
    [4] 施伟, 周强, 刘斌. 基于旋转永磁体的超低频机械天线电磁特性分析[J]. 物理学报, 2019, 68(18): 188401.SHI Wei, ZHOU Qiang, and LIU Bin. Performance analysis of spinning magnet as mechanical antenna[J]. Acta Physica Sinica, 2019, 68(18): 188401.
    [5] 周强, 姚富强, 施伟, 等. 机械式低频天线机理及其关键技术研究[J]. 中国科学: 技术科学, 2020, 50(1): 69–84. doi:  10.1360/SST-2019-0118ZHOU Qiang, YAO Fuqiang, SHI Wei, et al. Research on mechanism and key technology of mechanical antenna for a low-frequency transmission[J]. Scientia Sinica Technologica, 2020, 50(1): 69–84. doi:  10.1360/SST-2019-0118
    [6] 崔岸婧, 李道京, 周凯, 等. 阵列结构下的低频信号合成方法研究[J]. 物理学报, 2020, 69(19): 194101. doi:  10.7498/aps.69.20200501CUI Anjing, LI Daojing, ZHOU Kai, et al. Research on the method of composing low frequency signals based on array structures[J]. Acta Physica Sinica, 2020, 69(19): 194101. doi:  10.7498/aps.69.20200501
    [7] 王劲东, 薛洪波, 张艺腾, 等. 高精度航空地磁矢量测量技术[C]. 2018年中国地球科学联合学术年会论文集(二十八)——专题56: 煤炭资源与矿山地球物理、专题58: 地球物理探测方法与仪器新技术, 北京, 2018: 28, 44.WANG Jindong, XUE Hongbo, ZHANG Yiteng, et al. High precision airborne geomagnetic vector measurement technique[C]. Annual Meeting of Chinese Geoscience Union, Beijing, China, 2018: 28, 44.
    [8] HALLIDAY D, RESNICK R, and WALKER J. Fundamentals of Physics[M]. Weinheim: Wiley, 2014: 1135–1137.
    [9] 别业广. 电磁波的多普勒效应[J]. 物理与工程, 2003, 13(4): 62, 32.BIE Yeguang. The Doppler effect of electromagnetic waves[J]. Physics and Engineering, 2003, 13(4): 62, 32.
    [10] 高炳坤, 王凤林. 相对论多普勒效应的简易推导[J]. 大学物理, 2003, 22(8): 15–16.GAO Bingkun and WANG Fenglin. A concise deduction of relativistic Doppler effect[J]. College Physics, 2003, 22(8): 15–16.
    [11] 严欣达, 程先卿. 由相位不变性讨论光波的多普勒效应[J]. 大学物理, 1987, (11): 25.YAN Xinda and CHENG Xianqing. The Doppler effect of light wave is discussed from phase invariance[J]. College Physics, 1987, (11): 25.
    [12] 王景雪, 汤正新, 陈庆东, 等. 基于同时的相对性对钟慢尺缩效应的再认识[J]. 大学物理, 2009, 28(10): 24–27.WANG Jingxue, TANG Zhengxin, CHEN Qingdong, et al. Recognition on time dilation and length contraction effect based on relativity of simultaneity[J]. College Physics, 2009, 28(10): 24–27.
    [13] 张三慧. 时间膨胀与多普勒效应[J]. 大学物理, 1990, (9): 31–32.ZHANG Sanhui. Time dilation and Doppler effect[J]. College Physics, 1990, (9): 31–32.
    [14] 王建, 郑一农, 何子远. 阵列天线理论与工程应用[M]. 北京: 电子工业出版社, 2015: 8.WANG Jian, ZHENG Yinong, and HE Ziyuan. Antenna Array Theory and Engineering Applications[M]. Beijing: Publishing House of Electronics Industry, 2015: 8.
    [15] 丁鹭飞, 耿富录, 陈建春. 雷达原理[M]. 5版. 北京: 电子工业出版社, 2014: 253–254.DING Lufei, GENG Fulu, and CHEN Jianchun. Radar Principles[M]. 5th ed. Beijing: Publishing House of Electronics Industry, 2014: 253–254.
    [16] 吴翊, 朱炬波, 易东云, 等. 关于多普勒频率转换成距离变化率公式的讨论[J]. 中国空间科学技术, 1997, (6): 45–49.WU Yi, ZHU Jubo, YI Dongyun, et al. Discuss on formula of Doppler frequency to velocity[J]. Chinese Space Science and Technology, 1997, (6): 45–49.
    [17] 魏钟铨. 合成孔径雷达卫星[M]. 北京: 科学出版社, 2001: 204–206.WEI Zhongquan. Synthetic Aperture Radar Satellite[M]. Beijing: Science Press, 2001: 204–206.
    [18] 保铮, 邢孟道, 王彤. 雷达成像技术[M]. 北京: 电子工业出版社, 2005: 125–132.BAO Zheng, XING Mengdao, and WANG Tong. Radar Imaging Technology[M]. Beijing: Publishing House of Electronics Industry, 2005: 125–132.
  • [1] 郭忠义, 汪彦哲, 郑群, 尹超逸, 杨阳, 宫玉彬.  涡旋电磁波天线技术研究进展 . 雷达学报, 2019, 8(5): 631-655. doi: 10.12000/JR19091
    [2] 周成伟, 郑航, 顾宇杰, 王勇, 史治国.  互质阵列信号处理研究进展:波达方向估计与自适应波束成形 . 雷达学报, 2019, 8(5): 558-577. doi: 10.12000/JR19068
    [3] 刘玉琪, 易建新, 万显荣, 程丰, 饶云华, 龚子平.  数字电视外辐射源雷达多旋翼无人机微多普勒效应实验研究 . 雷达学报, 2018, 7(5): 585-592. doi: 10.12000/JR18062
    [4] 李世超, 侯培培, 屈俭, 郝丛静, 贾渠, 李刚, 李超.  基于波导缝隙阵列的新型太赫兹频率扫描天线 . 雷达学报, 2018, 7(1): 119-126. doi: 10.12000/JR17098
    [5] 刘向阳, 杨君刚, 孟进, 张晓, 牛德智.  低信噪比下基于Hough变换的前视阵列SAR稀疏三维成像 . 雷达学报, 2017, 6(3): 316-323. doi: 10.12000/JR17011
    [6] 韩壮志, 宋春吉, 侯建强.  伪码族复合连续波信号的多分辨率特性分析 . 雷达学报, 2016, 5(3): 278-283. doi: 10.12000/JR15100
    [7] 路满, 宋红军, 罗运华.  基于调频连续波信号的圆弧式合成孔径雷达成像方法 . 雷达学报, 2016, 5(4): 425-433. doi: 10.12000/JR16007
    [8] 李烈辰, 李道京, 黄平平.  基于变换域稀疏压缩感知的艇载稀疏阵列天线雷达实孔径成像 . 雷达学报, 2016, 5(1): 109-117. doi: 10.12000/JR14159
    [9] 陈莹, 钟菲, 郭树旭.  非合作跳频信号参数的盲压缩感知估计 . 雷达学报, 2016, 5(5): 531-537. doi: 10.12000/JR15106
    [10] 覃尧, 黄春琳, 陆珉, 徐伟.  基于小波变换与主成分分析的探地雷达自适应杂波抑制方法研究 . 雷达学报, 2015, 4(4): 445-451. doi: 10.12000/JR15013
    [11] 邹宁, 徐争光, 冉建华, 李朝阳.  基于带通欠采样的脉冲超宽带信号重建算法性能研究 . 雷达学报, 2015, 4(2): 185-191. doi: 10.12000/JR14103
    [12] 李万春, 黄成峰.  基于角度和多普勒频率的外辐射源定位系统的接收器最优航迹分析 . 雷达学报, 2014, 3(6): 660-665. doi: 10.12000/JR14118
    [13] 折小强, 仇晓兰, 韩冰, 雷斌.  一种基于变换域的滑动聚束SAR 调频率估计方法 . 雷达学报, 2014, 3(4): 419-427. doi: 10.3724/SP.J.1300.2014.14008
    [14] 曾操, 梁思嘉, 王威, 徐青.  基于频率步进信号的旋转式合成孔径雷达成像方法 . 雷达学报, 2014, 3(4): 401-408. doi: 10.3724/SP.J.1300.2014.14043
    [15] 关欣, 仲利华, 胡东辉, 丁赤飚.  一种基于RSPWVD-Hough 变换的无源雷达多普勒展宽补偿方法 . 雷达学报, 2013, 2(4): 430-438. doi: 10.3724/SP.J.1300.2013.13073
    [16] 孙红梅, 陈广东, 张弓.  基于测控电磁波极化特征信息的飞行器姿态估计 . 雷达学报, 2013, 2(4): 466-475. doi: 10.3724/SP.J.1300.2013.13087
    [17] 邓冬虎, 张群, 罗迎, 李松, 朱仁飞.  双基地ISAR 系统中分辨率分析及微多普勒效应研究(英文) . 雷达学报, 2013, 2(2): 152-167. doi: 10.3724/SP.J.1300.2013.13039
    [18] 万显荣.  基于低频段数字广播电视信号的外辐射源雷达发展现状与趋势 . 雷达学报, 2012, 1(2): 109-123. doi: 10.3724/SP.J.1300.2012.20027
    [19] 沈绍祥, 叶盛波, 方广有.  基于无载频脉冲雷达信号的高速数字采样方法与实现 . 雷达学报, 2012, 1(2): 136-142. doi: 10.3724/SP.J.1300.2012.20022
    [20] 龚斌, 王壮, 程翥.  基于功率聚焦的宽带阵列信号检测算法 . 雷达学报, 2012, 1(3): 253-261. doi: 10.3724/SP.J.1300.2012.20049
  • 加载中
图(17) / 表 (2)
计量
  • 文章访问数:  65
  • HTML全文浏览量:  44
  • PDF下载量:  23
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-06-18
  • 修回日期:  2020-09-13
  • 网络出版日期:  2020-10-12
  • 刊出日期:  2020-10-28

交错阵列甚低频信号产生方法研究

doi: 10.12000/JR20082
    基金项目:  国家部委基金,中国科学院空天信息创新研究院课题(Y910340Z2F)
    作者简介:

    崔岸婧(1997–),女,中国科学院大学博士生。研究方向为阵列天线与信号处理。E-mail: ajcui@qq.com

    李道京(1964–),男,中国科学院空天信息创新研究院,研究员,博士生导师,研究方向为雷达系统与雷达信号处理。E-mail: lidj@aircas.ac.cn

    周 凯(1995–),男,中国科学院大学博士生,研究方向为信号与信息处理。E-mail: zk_6810@163.com

    王 宇(1976–),男,中国科学院空天信息创新研究院,研究员,研究方向为机/星载SAR定标,定量遥感应用。E-mail: wangyu@mail.ie.ac.cn

    洪 峻(1960–),男,中国科学院空天信息创新研究院,研究员,博士生导师,研究方向为SAR定标技术、定量化SAR图像应用等。E-mail: jhong@mail.ie.ac.cn

    通讯作者: 李道京 lidj@aircas.ac.cn
  • 责任主编:廖桂生 Corresponding Editor: LIAO Guisheng
  • 中图分类号: TN951

摘要: 10 kHz量级甚低频电磁波信号具有较强的地物穿透能力,可用于地质勘探。由于其天线尺寸在10 km量级,其应用场合受到限制,研究基于适当尺寸高频雷达天线的甚低频电磁波信号产生方法具有重要意义。该文提出基于高频阵列天线产生甚低频信号的概念,利用阵列天线合成产生近光速远离运动雷达多普勒信号,实现信号频率的大幅降低。给出了发射波形、交错阵列设计和阵列参数选择方法。将周期脉冲串信号作为辐射单元信号,增大合成信号脉宽。利用阵列产生的脉宽展宽量填补脉冲信号的休止期,在目标区合成时间连续的甚低频信号。采用峰值旁瓣比(PSLR)、积分旁瓣比(ISLR)、阵列发射信号与合成信号的频谱对比评价合成低频信号的性能和发射信号的能量利用率。该文仿真了百米量级阵列100 MHz辐射单元信号在目标区合成10 kHz甚低频信号的情况:9行阵列构成交错阵列、辐射单元信号脉宽设置为0.115 μs时,合成信号频谱的峰值旁瓣比和积分旁瓣比分别为–13.34 dB和–9.44 dB, 10 kHz低频信号在合成信号中的能量占比为89.79%。该文分析了辐射单元间距误差、辐射单元信号时间、相位与幅度误差以及目标偏离预定位置的影响。仿真结果表明了该文方法的有效性。

注释:
1)  责任主编:廖桂生 Corresponding Editor: LIAO Guisheng

English Abstract

崔岸婧, 李道京, 周凯, 等. 交错阵列甚低频信号产生方法研究[J]. 雷达学报, 2020, 9(5): 925–938. doi:  10.12000/JR20082
引用本文: 崔岸婧, 李道京, 周凯, 等. 交错阵列甚低频信号产生方法研究[J]. 雷达学报, 2020, 9(5): 925–938. doi:  10.12000/JR20082
CUI Anjing, LI Daojing, ZHOU Kai, et al. Research on the method of composing very low frequency signals based on the staggered array[J]. Journal of Radars, 2020, 9(5): 925–938. doi:  10.12000/JR20082
Citation: CUI Anjing, LI Daojing, ZHOU Kai, et al. Research on the method of composing very low frequency signals based on the staggered array [J]. Journal of Radars, 2020, 9(5): 925–938. doi:  10.12000/JR20082
    • 10 kHz量级甚低频电磁波信号具有较强的地物穿透能力,可用于地质勘探[1-3]。传统天线的辐射单元尺寸需达到1/4波长,否则不能有效辐射电磁波。频率为10 kHz电磁波信号的波长为30 km,其天线辐射单元尺寸在7~8 km,天线尺寸在10 km量级,这使其应用受到限制,研究基于适当尺寸高频雷达天线的甚低频电磁波信号产生方法具有重要意义。

      甚低频电磁波信号产生一直是个难题,近年国外提出了一种采用机械旋转永磁体[4,5]的方式实现小型甚低频发射天线方法。该方法将机械能转换为电能,利用永磁体机械旋转直接激励电磁波,可突破天线物理尺寸限制。由于该方法涉及多学科交叉,国内外对该方法的研究均处于起步阶段,其性能尚未验证。

      在文献[6]的基础上,本文将电磁波多普勒效应与交错阵列结合,通过对交错阵列中各辐射单元信号的波形、时序、相位和周期,以及阵列数等参数的控制,提出了一种在目标区合成甚低频信号的方法。该甚低频信号与地物相互作用后,可通过磁探仪[7]接收处理,用于地质分析。

    • 文献[8-11]推导了电磁波多普勒效应的原理。以运动雷达为信号源,接收装置位于目标区,二者均位于静止坐标系$K$中,初始坐标分别为$\left( {{R_0},0,0} \right)$$\left( {0,0,0} \right)$。以雷达为原点建立运动坐标系$K'$,且在运动过程中,雷达始终位于$K'$系原点,$K'$系与$K$系的时间关系符合钟慢效应[12]。记$K'$系中雷达发射信号的频率为${f_0}$,由雷达运动所产生的多普勒频率为${f_{\rm{d}}}$, $K$系中接收信号频率为${f_{0{\rm{r}}}}$。如图1所示,左侧为接收装置,右侧为雷达与运动坐标系$K'$。雷达在$K$系中以速度$v$沿$x$轴正方向远离接收装置运动,且在运动过程中,$K'$系的$x'$轴与$K$系的$x$轴始终重合。

      图  1  雷达运动示意图

      Figure 1.  Radar motion diagram

      以下推导电磁波多普勒效应中脉宽与频率的变化[13]。设雷达发射信号在$K'$系中的脉宽为${\tau _{{\rm{e}}K'}}$,在$K$系中的对应脉宽为${\tau _{{\rm{e}}K}}$, $K$系中接收信号的脉宽为${\tau _{{\rm{r}}K}}$,其中${\tau _{{\rm{e}}K'}}$${\tau _{{\rm{r}}K}}$的关系符合常见的电磁波多普勒效应,信号频率的变化也由${\tau _{{\rm{e}}K'}}$${\tau _{{\rm{r}}K}}$推导,${\tau _{{\rm{e}}K}}$${\tau _{{\rm{r}}K}}$的关系表示$K$系中信号脉宽的变化。

      由多普勒效应,当雷达远离目标运动时,$K$系中接收信号脉宽${\tau _{{\rm{r}}K}}$$K'$系中发射信号脉宽${\tau _{{\rm{e}}K'}}$的关系为

      $${\tau _{{\rm{r}}K}} = \sqrt {\frac{{{\rm{c}} + v}}{{{\rm{c}} - v}}} {\tau _{{\rm{e}}K'}}$$ (1)

      其中,$v$为雷达运动速度,${\rm{c}}$为光速。由钟慢效应,雷达发射信号在$K$系与$K'$系中脉宽的对应关系为

      $$ {\tau _{{\rm{e}}K}} = \frac{{{\tau _{{\rm{e}}K'}}}}{{\sqrt {1 - {{\left( {\dfrac{v}{{\rm{c}}}} \right)}^2}} }} $$ (2)

      将式(2)代入式(1)可得K系中雷达发射信号与目标处接收信号脉宽的关系

      $${\tau _{{\rm{r}}K}} = \sqrt {\frac{{{\rm{c}} + v}}{{{\rm{c}} - v}}} {\tau _{{\rm{e}}K'}} = \frac{{{\rm{c}} + v}}{{\rm{c}}}{\tau _{{\rm{e}}K}}$$ (3)

      因此雷达发射信号在$K'$系和$K$系中的脉宽展宽量为分别为

      $$ \Delta {\tau _{K'}} = {\tau _{{\rm{r}}K}} - {\tau _{{\rm{e}}K'}} = \left( {\sqrt {\frac{{{\rm{c}} + v}}{{{\rm{c}} - v}}} - 1} \right){\tau _{{\rm{e}}K'}} $$ (4)
      $$ \Delta {\tau _K} = {\tau _{{\rm{r}}K}} - {\tau _{{\rm{e}}K}} = \frac{v}{{\rm{c}}}{\tau _{{\rm{e}}K}} $$ (5)

      式(5)表示静止坐标系$K$中雷达发射信号脉宽变化。

      由于$K'$系中发射信号与$K$系中接收信号脉冲中信号的周期数不变,因此信号频率变化由$K'$系中雷达发射信号脉宽和$K$系中接收信号脉宽可得

      $$ \frac{{{f_{0{\rm{r}}}}}}{{{f_0}}} = \frac{{{\tau _{{\rm{e}}K'}}}}{{{\tau _{{\rm{r}}K}}}} = \sqrt {\frac{{{\rm{c}} - v}}{{{\rm{c}} + v}}} $$ (6)

      对应的多普勒频率为

      $$ {f_{\rm{d}}} = {f_{0{\rm{r}}}} - {f_0} = \left( {\sqrt {\frac{{{\rm{c}} - v}}{{{\rm{c}} + v}}} - 1} \right){f_0} $$ (7)

      当雷达远离目标运动的速度接近电磁波速度${\rm{c}}$时,接收信号频率将会明显降低。

      若雷达与目标的初始距离为30 km,发射载频100 MHz脉宽为0.05 μs的信号,同时以速度$v = \dfrac{{1{0^8} - 1}}{{1{0^8} + 1}}{\rm{c}}$远离目标运动,多普勒频率为–99.99 MHz,目标区接收信号频率为10 kHz,脉宽为500 μs,雷达发射信号与目标区接收信号的波形和频谱如图2所示。

      图  2  多普勒效应中的发射/接收信号波形与频谱

      Figure 2.  The emission/received signal waveform and spectrum of Doppler effect

    • 根据对电磁波多普勒效应的理解,如图3(a)所示,在图1的基础上,在静止坐标系中分解运动雷达发射信号的过程。图中的$x$轴为静止坐标系$K$$x$轴,$t$轴表示雷达运动的时间,设雷达在$t = 0$时刻位于$x = {R_0}$处,并从该时刻开始在$K$系中以速度$v$沿$x$轴正方向运动并发射信号。以时间间隔$\Delta t$将雷达的运动过程进行分解,则每个时间间隔内雷达的运动距离为$v\Delta t$,此时雷达运动过程变为步进过程。如图3(b)所示,用天线阵列近似步进的雷达运动过程,天线阵列和接收装置均位于$K$系中。

      图  3  运动雷达发射信号过程分解和天线阵列等效雷达运动示意图

      Figure 3.  Schematic diagram of the decomposed process of radar transmitting signals and the equivalent radar motion generated by the antenna array

      采用多行阵列构成交错阵列,可缩短辐射单元间距,减小离散化的时间间隔,使得合成信号更接近所需的低频信号。如图4所示,在直角坐标系中建立交错阵列结构。${N_{\rm{l}}}$行阵列在$xoz$平面内错位排布,单行阵列沿$x$轴负方向排布,相邻阵列在$x$轴以距离${d_{\rm{s}}}$错位,沿$z $轴正方向以间距$h = {\lambda / 2}$排布,$\lambda $为辐射单元信号的载波波长。单行阵列的阵长为${L_0}$,辐射单元间距${d_0} = {\lambda / 2}$,辐射单元个数${N_0}$。接收装置位于目标区,目标位于45°扫描角[14,15]方向,坐标为$\left( {{x_0},{y_0},0} \right)$,与阵列近端距离为${R_0}$。当单行阵列辐射单元间距${d_0}$为错位距离${d_{\rm{s}}}$${N_{\rm{l}}}$倍时,由于$h \ll {R_0}$,阵列在$z$轴上的距离影响极小,所以${N_{\rm{l}}}$行阵列等效在$x$轴构成辐射单元间距${d_{\rm{s}}}$、阵长为${L_{\rm{s}}} = {L_0} + \left( {{N_{\rm{l}}} - 1} \right){d_{\rm{s}}}$的交错阵列。

      图  4  交错阵列结构

      Figure 4.  Staggered array structure

      辐射单元记为${T_{{n_{\rm{l}}}{n_0}}}$,其中${n_{\rm{l}}}$${n_0}$分别表示该辐射单元所在的阵列与阵列中的位置(${n_{\rm{l}}} = 0,1, ··· ,{N_{\rm{l}}} - 1$; ${n_0} = 0,1, ··· ,{N_0} - 1$)。通过交错阵列中近端至远端的辐射单元依次发射信号,等效雷达在阵列上的高速运动。单行阵列与交错阵列中,辐射单元发射信号的时间间隔分别为${{{d_0}} / v}$${{{d_{\rm{s}}}} / v}$, $0 \le v < {\rm{c}}$。图中辐射单元发射信号的顺序为${T_{00}}$, ${T_{10}}$, ···, ${T_{\left( {{N_{\rm{l}}} - 2} \right)\left( {{N_0} - 1} \right)}}$, ${T_{\left( {{N_{\rm{l}}} - 1} \right)\left( {{N_0} - 1} \right)}}$

    • (1) 辐射单元信号波形

      阵长产生的辐射单元信号脉宽展宽量为

      $$ \Delta {\tau _{\rm{L}}} = \frac{{\sqrt {{{\left( {{x_0} + {L_{\rm{s}}}} \right)}^2} + y_0^2} - {R_0}}}{{\rm{c}}} + \frac{{{L_{\rm{s}}}}}{v} $$ (8)

      令辐射单元信号脉宽等于雷达发射信号在$K$系中的脉宽${\tau _{{\rm{e}}K}}$, $\Delta {\tau _{\rm{L}}}$等效多普勒效应中信号脉宽展宽量$\Delta {\tau _K}$。由式(5),当$v \to {\rm{c}}$时,$\Delta {\tau _K} \approx {\tau _{{\rm{e}}K}}$,即辐射单元信号脉宽与其展宽量近似相等,因此当辐射单元发射单脉冲信号时,目标区合成信号的脉宽为

      $$ \tau ' \!=\! {\tau _{{\rm{e}}K}} + \Delta {\tau _K} \!=\! 2\Delta {\tau _{\rm{L}}} \approx 2\frac{{\sqrt {{{\left( {{x_0} + {L_{\rm{s}}}} \right)}^2} \!+\! y_0^2} - {R_0} \!+\! {L_{\rm{s}}}}}{{\rm{c}}} $$ (9)

      为获得甚低频信号,合成信号需要有足够大的脉宽。在阵长固定的条件下,直接增大辐射单元信号脉宽将导致合成信号中载波频率分量的增大,影响甚低频信号对阵列发射信号的能量利用率。

      为增大合成信号的脉宽,令辐射单元发射周期脉冲串信号,通过增加辐射单元信号的周期数实现合成信号脉宽的增大。利用阵列产生的脉宽展宽量填补周期脉冲串信号的休止期,使得目标区合成信号在其时宽内连续。辐射单元信号波形如图5所示。辐射单元信号周期和脉宽分别记为${T_0}$${\tau _0}$, 3.2.2节中将推导${T_0}$${\tau _0}$的取值。辐射单元信号每个脉冲前沿的相位根据多普勒效应设置。每个辐射单元信号脉冲内包含多个辐射单元信号载波周期${1}/{{{f_0}}}$

      图  5  辐射单元信号示意图

      Figure 5.  Schematic diagram of the radiation unit signal

      用快时间和慢时间描述辐射单元信号与目标区的合成信号。记同一辐射单元信号的时间为快时间$\hat t$,不同辐射单元信号的时间为慢时间${t_{\rm{m}}}$,目标区信号的时间为${t_{\rm{r}}}$,三者之间的关系为

      $${t_{\rm{r}}} = \hat t + {t_{\rm{m}}}$$ (10)

      图5中描述单个辐射单元信号,因此其横坐标为快时间$\hat t$

      (2) 辐射单元信号脉宽与周期

      由式(8)可得阵列产生的脉宽展宽量$\Delta {\tau _{\rm{L}}}$。根据式(5),当$v = {\rm{c}}$时,静止坐标系$K$中雷达发射信号脉宽与其脉宽展宽量相等,即${\tau _{{\rm{e}}K}} = \Delta {\tau _K}$,此时${\tau _{{\rm{e}}K}}$占接收信号脉宽${\tau _{{\rm{r}}K}}$的50%。令辐射单元信号脉宽${\tau _0}$、周期${T_0}$和阵列产生的脉宽展宽量$\Delta {\tau _{\rm{L}}}$分别等效${\tau _{{\rm{e}}K}}$, ${\tau _{{\rm{r}}K}}$$\Delta {\tau _K}$。由于合成甚低频信号的仿真中$v \to {\rm{c}}$,则有${\tau _0} \approx \Delta {\tau _{\rm{L}}}$, ${\tau _0} \approx \dfrac{1}{2}{T_0}$,因此将辐射单元发射的周期脉冲串信号的占空比选定为50%,辐射单元信号脉宽为

      $$ {\tau _0} = \Delta {\tau _{\rm{L}}} = \frac{{\sqrt {{{\left( {{x_0} + {L_{\rm{s}}}} \right)}^2} + y_0^2} - {R_0}}}{{\rm{c}}} + \frac{{{L_{\rm{s}}}}}{v} $$ (11)

      辐射单元信号周期为

      $$ {T_0} = 2{\tau _0} = 2\Delta {\tau _{\rm{L}}} $$ (12)

      此时$\Delta {\tau _{\rm{L}}}$恰好填补辐射单元信号的休止期,目标区合成信号的时间连续。

      若将辐射单元信号周期减小为

      $$ {T_0}' = \frac{1}{u}{T_0} $$ (13)

      其中$u > 1$,此时辐射单元信号的休止期$\dfrac{1}{2}{T_0}'$<$\Delta {\tau _{\rm{L}}}$,因此阵列产生的脉宽展宽量能够实现周期脉冲串信号休止期的填补,使得目标区合成信号的时间连续。

      通过增大$u$,减小辐射单元信号周期,增加辐射单元信号周期数,等效增加合成信号相位调制次数,实现谐波的抑制,同时减少叠加所导致的信号抵消,提高合成信号对阵列发射信号的能量利用率。4.1节中仿真分析了辐射单元信号周期对目标区合成信号性能的影响。

      (3) 合成信号脉宽

      合成信号由交错阵列中各辐射单元发射的周期脉冲串信号构成。记各辐射单元发射周期脉冲串信号的周期数为${N_{\rm{p}}}$,合成信号脉宽为

      $$ \begin{split} {\tau _{{\rm{Lp}}}} =\,& \left[ {\left( {{N_{\rm{p}}} - 1} \right){T_0} + {\tau _0}} \right] \\ & + \left[ {\frac{{\sqrt {{{\left( {{x_0} + {L_{\rm{s}}}} \right)}^2} + y_0^2} - {R_0}}}{{\rm{c}}} + \frac{{{L_{\rm{s}}}}}{v}} \right] \end{split} $$ (14)

      其中第1项为各辐射单元发射信号的时长,第2项为交错阵列近端与远端辐射单元信号传播至目标区的时间差。

    • 由多普勒频率与辐射单元斜距可得辐射单元信号[16]和目标区合成信号表达式。辐射单元${T_{{n_{\rm{l}}}{n_0}}}$与目标之间的距离为

      $$ {r_{{n_{\rm{l}}}{n_0}}} = \sqrt {{{\left( {{x_0} + {n_{\rm{l}}}{d_{\rm{s}}} + {n_0}{d_0}} \right)}^2} + y_0^2 + {{\left( {{n_{\rm{l}}}h} \right)}^2}} $$ (15)

      当目标位于$xoy$平面时,由于${\left( {{n_{\rm{l}}}h} \right)^2} \ll x_0^2 + y_0^2 = R_0^2$,因此可忽略辐射单元在$z$轴上距离的影响,即式(15)可简化为

      $${r_{{n_{\rm{l}}}{n_0}}} = \sqrt {{{\left( {{x_0} + {n_{\rm{l}}}{d_{\rm{s}}} + {n_0}{d_0}} \right)}^2} + y_0^2} $$ (16)

      设交错阵列近端的辐射单元${T_{00}}$${t_{\rm{m}}} = 0$时刻开始发射信号,${T_{{n_{\rm{l}}}{n_0}}}$开始发射信号的时刻为${t_{\rm{m}}} = \dfrac{{{n_{\rm{l}}}{d_{\rm{s}}} + {n_0}{d_0}}}{v}$,该辐射单元的信号在${t_{\rm{m}}} = \dfrac{{{r_{{n_{\rm{l}}}{n_0}}}}}{{\rm{c}}} + \dfrac{{{n_{\rm{l}}}{d_{\rm{s}}} + {n_0}{d_0}}}{v}$时传播至目标区。

      由式(6)和式(15),可得辐射单元信号的表达式

      $$ \begin{split} {s_{{n_{\rm{l}}}{n_0}}}\left( {\hat t} \right) =\,& \sum\limits_{{n_{\rm{p}}} = 0}^{{N_{\rm{p}}} - 1} {\rm{rect}}\left( {\frac{{\hat t - {n_{\rm{p}}}{T_0} - 0.5{\tau _0}}}{{{\tau _0}}}} \right)\\ & \cdot\exp \Bigr[ {{\rm{j}}\left( {2{\rm{\pi}} {f_0}\hat t + \varphi _0^{{n_{\rm{l}}}{n_0}{n_{\rm{p}}}}} \right)} \Bigr] \end{split} $$ (17)
      $$ \varphi _0^{{n_{\rm{l}}}{n_0}{n_{\rm{p}}}} = 2{\rm{\pi}} \sqrt {\frac{{{\rm{c}} - v}}{{{\rm{c}} + v}}} {f_0}t_{\rm{m}}^{{n_{\rm{l}}}{n_0}} + 2{\rm{\pi}} {f_{\rm{d}}}{n_{\rm{p}}}{T_0} $$ (18)

      其中$t_{\rm{m}}^{{n_{\rm{l}}}{n_0}} = \dfrac{{{r_{{n_{\rm{l}}}{n_0}}}}}{{\rm{c}}} + \dfrac{{{n_{\rm{l}}}{d_{\rm{s}}} + {n_0}{d_0}}}{v}$, ${n_{\rm{p}}}$为辐射单元信号的周期编号,$\varphi _0^{{n_{\rm{l}}}{n_0}{n_{\rm{p}}}}$为辐射单元${T_{{n_{\rm{l}}}{n_0}}}$发射第${n_{\rm{p}}}$个脉冲的初始相位。

      将式(10)代入式(17)可得目标区辐射单元${T_{{n_{\rm{l}}}{n_0}}}$的信号表达式

      $$ \begin{split} {s_{{\rm{r}}{n_{\rm{l}}}{n_0}}}\left( {{t_{\rm{r}}}} \right) =\,& \sum\limits_{{n_{\rm{p}}} = 0}^{{N_{\rm{p}}} - 1} {\rm{rect}}\left( {\frac{{{t_{\rm{r}}} - t_{\rm{m}}^{{n_{\rm{l}}}{n_0}} - {n_{\rm{p}}}{T_0} - 0.5{\tau _0}}}{{{\tau _0}}}} \right)\\ & \cdot\exp \Bigr\{ {{\rm{j}}\left[ {2{\rm{\pi}} {f_0}\left( {{t_{\rm{r}}} - t_{\rm{m}}^{{n_{\rm{l}}}{n_0}}} \right) + \varphi _0^{{n_{\rm{l}}}{n_0}{n_{\rm{p}}}}} \right]} \Bigr\} \\ = \,&\sum\limits_{{n_{\rm{p}}} = 0}^{{N_{\rm{p}}} - 1} {\rm{rect}}\left( {\frac{{{t_{\rm{r}}} - t_{\rm{m}}^{{n_{\rm{l}}}{n_0}} - {n_{\rm{p}}}{T_0} - 0.5{\tau _0}}}{{{\tau _0}}}} \right)\\ & \cdot\exp \Bigr\{ {{\rm{j}}\left[ {2{\rm{\pi}} {f_0}{t_{\rm{r}}} + 2{\rm{\pi}} {f_{\rm{d}}}\left( {t_{\rm{m}}^{{n_{\rm{l}}}{n_0}} + {n_{\rm{p}}}{T_0}} \right)} \right]} \Bigr\} \end{split} $$ (19)

      目标区的合成信号等于交错阵列中各辐射单元信号的叠加,因此目标区合成信号的表达式为

      $$ {s_{\rm{r}}}\left( {{t_{\rm{r}}}} \right) = \sum\limits_{{n_{\rm{l}}} = 0}^{{N_{\rm{l}}} - 1} {\sum\limits_{{n_0} = 0}^{{N_0} - 1} {{s_{{\rm{r}}{n_{\rm{l}}}{n_0}}}\left( {{t_{\rm{r}}}} \right)} } $$ (20)

      式(19)中,$\varphi _{\rm{d}}^{{n_{\rm{l}}}{n_0}{n_{\rm{p}}}} = 2\pi {f_{\rm{d}}}\left( {t_{\rm{m}}^{{n_{\rm{l}}}{n_0}} + {n_{\rm{p}}}{T_0}} \right)$等效多普勒相位。$\varphi _{\rm{d}}^{{n_{\rm{l}}}{n_0}{n_{\rm{p}}}}$的引入导致空间中合成信号的相位变化,因此合成信号为宽谱信号。

      由于辐射单元信号仅设置各脉冲的初始相位,脉冲内没有相位调制,因此辐射单元信号属于窄带信号,可有效辐射。增大辐射单元的带宽有益于保证合成信号的性能。

      若假设辐射单元带宽为中心频率的10%,为保证辐射单元信号带宽满足要求,限制辐射单元信号的周期与脉宽

      $$ \qquad\qquad\quad{\tau _0} = \frac{1}{{2u}}{T_0} $$ (21)
      $$ \qquad\qquad\quad\frac{1}{{{\tau _0}}} = \frac{{2u}}{{{T_0}}} \le \frac{{{f_0}}}{{10}} $$ (22)

      即辐射单元信号周期与脉宽的最小值分别为${{20}}/{{{f_0}}}$${{10}}/{{{f_0}}}$

      结合实现空间加密的交错阵列结构和设置确定相位的辐射单元信号,使得合成信号相位在快时间上快速变化,带宽有限的辐射单元信号在空间构成宽谱的合成信号。

    • (1) 辐射单元信号周期

      以下利用峰值旁瓣比(Peak Side Lobe Ratio, PSLR)、积分旁瓣比(Integrated Side Lobe Ratio, ISLR)[17]和发射信号能量利用率评价合成信号中谐波分量的影响。

      因为各辐射单元信号的形式相同,所以可用辐射单元${T_{{\rm{00}}}}$的信号等效其他辐射单元信号,由此可得交错阵列的发射信号。将归一化处理的阵列发射信号频谱和合成信号频谱进行对比,可分析发射信号的能量利用率。

      当目标位于$xoy$平面内且距离阵列近端30 km,辐射单元信号载波频率为100 MHz,单行阵列长度为120 m,交错阵列由5行阵列构成,其长度为121.2 m,单行阵列与交错阵列中辐射单元间距分别为1.5 m和0.3 m,辐射单元数分别为81和405,合成信号频率10 kHz时,分别仿真辐射单元周期为${{{T_0}} / 2}$${{{T_0}} / 6}$的情况下,合成信号的波形、频率与频谱对比图,仿真结果如图6图7表1所示。根据交错阵列阵长,可得辐射单元信号最大周期${T_0}$=1.38 μs。仿真中描述信号的采样频率为500 MHz(对应时间分辨率为2 ns),通过多行阵列交错排布,减少辐射单元间隔,实现空间加密,降低系统对时间分辨率的要求,提高该采样频率有助于对合成信号的描述。

      表 1  辐射单元信号周期取${{{{{T_0}}}} / {{2}}}$${{{{{T_0}}}} / {{6}}}$时的仿真参数与结果

      Table 1.  Simulation parameters and results when the radiation unit signal period being ${{{{{T_0}}}} / {{2}}}$ and ${{{{{T_0}} / 6}}}$ respectively

      参数辐射单元信号周期取${{{T_0}} / 2}$辐射单元信号周期取${{{T_0}} / 6}$
      最大辐射单元信号周期${T_0}$(μs)1.381.38
      辐射单元信号周期数(个)200600
      辐射单元信号脉宽(μs)0.3450.115
      辐射单元信号脉宽展宽(μs)0.3450.115
      辐射单元信号周期(μs)0.690.23
      辐射单元信号休止期(μs)0.3450.115
      合成信号脉宽(μs)138.36138.59
      峰值旁瓣比(dB)–11.45–13.34
      积分旁瓣比(dB)–4.33–8.77
      低频信号能量占比(%)73.0588.29
      频谱对比图中的10 kHz分量(dB)–14.610–5.081
      发射信号能量利用率(%)18.6055.71

      图  6  辐射单元信号周期为${{{T_0}} / 2}$时合成信号的波形、频谱与频谱对比图

      Figure 6.  Waveform, spectrum and spectral contrast diagram of the composite signal when the radiation unit signal period being ${{{T_0}} / 2}$

      图  7  辐射单元信号周期为${{{T_0}} / 6}$时合成信号的波形、频谱与频谱对比图

      Figure 7.  Waveform, spectrum and spectral contrast diagram of the composite signal when the radiation unit signal period being ${{{T_0}} / 6}$

      对比图6图7及其表1中的仿真结果可见,通过减小辐射单元信号周期,可抑制谐波和载波分量,减少由叠加导致的信号抵消,提高交错阵列发射信号能量的利用率。由于通常辐射单元带宽为中心频率的10%,为保证辐射单元信号脉宽满足要求,不再增大$u$或减小辐射单元信号周期。

      (2) 交错阵列行数和等效辐射单元间距

      当辐射单元信号周期取${{{T_0}} / 6}$=0.23 μs时,根据表1所示参数,仿真9行阵列构成的交错阵列在目标区合成信号的波形、频谱与频谱对比图。仿真结果如图8所示。

      图  8  辐射单元信号周期为${{{T_0}} / 6}$时9行阵列构成的交错阵列在目标区合成信号的波形、频谱与频谱对比图

      Figure 8.  Waveform, spectrum and spectral contrast diagram of the composite signal generated by the 9-row staggered array when the radiation unit signal period being ${{{T_0}} / 6}$

      图8中合成信号频谱的峰值旁瓣比为–13.34 dB,积分旁瓣比为–9.44 dB,低频信号在合成信号中的能量占比为89.79%。频谱对比图中,随机相位辐射单元信号所合成信号的频率分量分布于辐射单元信号载频附近,而当辐射单元信号相位根据多普勒效应设置时,其合成信号的10 kHz分量幅度为–2.542 dB,对应的发射信号能量利用率为74.63%。

      对比图7图8,通过多行阵列交错排布的方式减小辐射单元间距,可在步进方式下更精确模拟高速运动多普勒信号,抑制合成信号的谐波,并提高发射信号的能量利用率。

      上述信号产生过程中,信号包络也在同步移动,如图9所示,这与雷达成像[18]中的距离徙动类似。由于多普勒效应涉及频率的变化,所以本文甚低频信号产生过程已不属于线性过程。

      图  9  辐射单元信号包络示意图

      Figure 9.  Schematic diagram of the envelope of radiation unit signals

    • 分析辐射单元间距误差、辐射单元信号的时间、相位、幅度误差以及目标偏离预定位置对合成信号的影响,对实际应用具有意义。以下对合成信号的误差分析基于9行阵列构成交错阵列、辐射单元周期为${{{T_0}} / 6}$的前提。

      (1) 辐射单元间距误差

      当辐射单元间距误差服从均值为0,标准差为0.03 m的正态分布时,误差的分布、目标区合成信号的波形与频谱如图10所示。

      图  10  辐射单元间距误差分布、目标区合成信号的波形与频谱图

      Figure 10.  The distribution of spacing error, waveform and spectrum of the composite signal in target area

      (2) 辐射单元信号时间误差

      当辐射单元信号时间误差服从均值为0,标准差为1 ns的正态分布时,误差的分布、目标区合成信号的波形与频谱如图11所示。

      图  11  辐射单元信号时间误差分布、目标区合成信号的波形与频谱图

      Figure 11.  The distribution of radiation unit signal time error, waveform and spectrum of the composite signal in target area

      (3) 辐射单元信号相位误差

      当辐射单元信号相位误差服从均值为0,标准差${\pi }/{{{2^6}}}$的正态分布时,相位误差分布、目标区合成信号的波形与频谱如图12所示。

      图  12  辐射单元信号相位误差分布、目标区合成信号的波形与频谱图

      Figure 12.  The distribution of radiation unit signal phase error, waveform and spectrum of the composite signal in target area

      (4) 辐射单元信号幅度误差

      辐射单元增益的不同将导致辐射单元信号的幅度误差。当归一化电压幅度误差服从均值为0,标准差为0.05 V的正态分布时,对误差分布与合成信号的仿真如图13所示。

      图  13  辐射单元信号幅度误差分布、目标区合成信号的波形与频谱图

      Figure 13.  The distribution of radiation unit signal amplitude error, waveform and spectrum of the composite signal in target area

      (5) 目标偏离预定位置

      以上仿真中,目标均位于$xoy$平面内45°扫描角方向,与交错阵列近端距离30 km。若目标偏离该预定位置,将对合成信号产生影响。

      当目标位于$xoy$平面内45°扫描角方向,与交错阵列近端距离20 km和40 km时,合成信号的波形与频谱如图14图15所示。

      图  14  目标在xoy平面内45°扫描角方向距离阵列近端20 km时目标区合成信号的波形与频谱

      Figure 14.  Waveform and spectrum of the composite signal in the target area when the target is 20 km away from the near end of the array at 45° scanning angle in xoy plane

      图  15  目标在xoy平面内45°扫描角方向距离阵列近端40 km时目标区合成信号的波形与频谱

      Figure 15.  Waveform and spectrum of the composite signal in the target area when the target is 40 km away from the near end of the array at 45° scanning angle in xoy plane

      若目标在$xoy$平面上的投影为预定位置,在$z$轴上偏离1 km,由式(14)可得各辐射单元信号传播至实际目标区的慢时间为

      $$ \begin{split} t_{\rm{m}}^{{n_{\rm{l}}}{n_0}} =\,& \frac{{\sqrt {{{\left( {{x_0} + {n_{\rm{l}}}{d_{\rm{s}}} + {n_0}{d_0}} \right)}^2} + y_0^2 + {{\left( {{z_0} - {n_{\rm{l}}}h} \right)}^2}} }}{{\rm{c}}} \\ & + \frac{{{n_{\rm{l}}}{d_{\rm{s}}} + {n_0}{d_0}}}{v}\\[-17pt] \end{split} $$ (23)

      其中${z_0}$=1 km。由于辐射单元在$z$轴上的位置$z = {n_{\rm{l}}}h$${\left( {{z_0} - {n_{\rm{l}}}h} \right)^2}$有较大影响,因此当目标不在$xoy$平面内时,辐射单元在$z$轴上的排布不可忽略,仿真结果如图16所示。

      图  16  目标沿$z$轴偏离1 km时目标区合成信号的波形与频谱

      Figure 16.  Waveform and spectrum of the composite signal in the target area when the target deviates 1 km along the z-axis

      (6) 综合误差

      当目标位于$xoy$平面内45°扫描方向,距离交错阵列近端40 km时,不考虑交错阵列在$z$轴上的排布影响,若交错阵列及其信号同时受到辐射单元间距误差、辐射单元信号相位和幅度误差,目标区合成信号的波形和频谱如图17所示。

      图  17  受综合误差影响时目标区合成信号的波形与频谱

      Figure 17.  Waveform and spectrum of the composite signal in the target area influenced by synthetical error

      (7) 仿真结果分析

      表2给出了各误差影响下目标区合成信号频谱的峰值旁瓣比和积分旁瓣比。由仿真结果可见,辐射单元间距误差,辐射单元信号时间、相位、幅度误差和目标偏离预定位置均导致合成信号中谐波分量的增大,其中辐射单元信号时间误差和目标在$z$轴上偏离预定位置对合成信号的影响较大。

      表 2  不同误差影响下合成信号频谱参数(dB)

      Table 2.  Spectrum parameters of the composite signal influenced by different errors (dB)

      误差峰值旁瓣比积分旁瓣比
      辐射单元间距误差–13.31–8.95
      辐射单元信号时间误差–13.34–7.64
      辐射单元信号相位误差–13.34–9.07
      辐射单元信号幅度误差–13.34–9.33
      目标在xoy平面内距离阵列近端20 km–13.34–9.44
      目标在xoy平面内距离阵列近端40 km–13.34–9.44
      目标在z轴偏离1 km–13.34–8.97
      综合误差–13.28–6.63
    • 结合电磁波多普勒效应和交错阵列,本文研究了高频天线在目标区合成甚低频信号的问题。将辐射单元信号设计为周期脉冲串信号,通过阵列产生的脉宽展宽量填补信号的休止期,使得合成信号时间连续。通过多行阵列交错排布减小辐射单元间隔,降低系统对时间分辨率的要求。从仿真结果看,本文提出的交错阵列甚低频信号产生方法具有一定的可行性。但合成信号的谐波控制水平和发射信号能量利用率仍待进一步提高。

      致谢 感谢西安电子科技大学的水鹏朗老师对本文工作的帮助。

参考文献 (18)

目录

    /

    返回文章
    返回