基于变分模态分解与优选的超高分辨ISAR成像微多普勒抑制方法

李中余 桂亮 海宇 武俊杰 王党卫 王安乐 杨建宇

刘世锋, 徐晓瑞, 张方正, 等. 超低相噪光电振荡器及其频率综合技术研究[J]. 雷达学报, 2019, 8(2): 243–250. doi: 10.12000/JR19029
引用本文: 李中余, 桂亮, 海宇, 等. 基于变分模态分解与优选的超高分辨ISAR成像微多普勒抑制方法[J]. 雷达学报(中英文), 2024, 13(4): 852–865. doi: 10.12000/JR24043
LIU Shifeng, XU Xiaorui, ZHANG Fangzheng, et al. Ultralow phase noise optoelectronic oscillator and its application to a frequency synthesizer[J]. Journal of Radars, 2019, 8(2): 243–250. doi: 10.12000/JR19029
Citation: LI Zhongyu, GUI Liang, HAI Yu, et al. Ultrahigh-resolution ISAR micro-Doppler suppression methodology based on variational mode decomposition and mode optimization[J]. Journal of Radars, 2024, 13(4): 852–865. doi: 10.12000/JR24043

基于变分模态分解与优选的超高分辨ISAR成像微多普勒抑制方法

DOI: 10.12000/JR24043
基金项目: 国家自然科学基金(62171084),衢州市财政资助科研项目(2022D014)
详细信息
    作者简介:

    李中余,教授,博士生导师,研究方向为新体制雷达成像技术等

    桂 亮,硕士生,研究方向为逆合成孔径雷达成像及雷达微多普勒效应等

    海 宇,博士生,研究方向为超高分辨率雷达成像、微波光子雷达成像、稀疏信号恢复等

    武俊杰,教授,博士生导师,研究方向为合成孔径雷达成像、双/多基合成孔径雷达、雷达信号处理等

    王党卫,教授,研究方向为新体制雷达与信号处理等

    王安乐,副教授,研究方向为微波光子雷达等

    杨建宇,教授,博士生导师,研究方向为雷达信号处理、合成孔径雷达成像等

    通讯作者:

    李中余 zhongyu_li@uestc.edu.cn

    桂亮 lianggui46@qq.com

  • 责任主编:罗迎 Corresponding Editor: LUO Ying
  • 中图分类号: TN957

Ultrahigh-resolution ISAR Micro-Doppler Suppression Methodology Based on Variational Mode Decomposition and Mode Optimization

Funds: The National Natural Science Foundation of China (62171084), The Municipal Government of Quzhou (2022D014)
More Information
  • 摘要: 逆合成孔径雷达(ISAR)在对空中目标成像时,目标自身的转动、振动等局部微动将产生微多普勒效应,回波将附加额外的多普勒调制,造成频谱展宽。在超高分辨条件下,这一微动特性将会影响主体散射点的聚焦,导致目标图像局部散焦模糊,严重影响成像质量。并且,微多普勒相位还具有时变非平稳特性,难以从ISAR目标回波中准确估计或分离出微多普勒。为了解决上述问题,该文利用目标主体回波和微多普勒分量的时频分布差异,提出一种基于变分模态分解(VMD)与优选的非参数化方法抑制了回波中的微多普勒分量,消除了微多普勒对成像的影响,获得超高分辨率的无人机ISAR成像结果。该文首先引入VMD算法并将其扩展到复数域,将ISAR目标回波数据沿方位向分解为若干个中心频率均匀分布于多普勒采样带宽中的模函数,在此基础上利用图像熵指标优化分解参数和筛选成像模态,以保证微多普勒的良好抑制和主体回波的较完整保留。与现有基于经验模态分解(EMD)和局部均值分解(LMD)的方法相比,所提方法在超大带宽条件下对旋翼微动引起的微多普勒干扰有着更为出色的抑制效果,而且对机身部分的保留更为完整。最后,通过仿真对比和超宽带微波光子ISAR无人机实测数据处理,证明了该文所提方法的有效性和优势。

     

  • 作为现代通信、探测和测试系统中的核心部件,频率综合器能够在一定带宽范围内灵活输出高精度、高稳定的频率信号,广泛应用于无线通信、雷达探测、测试与计量、卫星等领域[14]。随着技术的发展,高频率、大带宽、低相噪及低杂散成为新一代高性能频率综合系统的必然演进方向。传统基于微波技术的频率综合技术采用低频晶体振荡器作为参考源,通过倍频锁相的方式产生高频信号,但是所产生的高频信号的相位噪声在倍频锁相过程中被快速恶化(20lgN, N为倍频因子),难以进一步提升。

    相比之下,光子技术具有低噪声、低损耗和大带宽等关键优势,因此通过光子技术实现频率综合已获得人们的广泛重视[59]。光域实现微波频率综合技术主要包括双光频梳锁定技术[5],飞秒锁模激光器锁相技术[6],半导体激光器注入锁定技术[7]和基于光电振荡器(Opto-Electronic Oscillator, OEO)的锁相技术[8,9]等。双光频梳锁定技术和飞秒锁模激光器锁相技术本质上都需要选出两个光波长进行拍频,将拍频产生的微波信号作为频率综合器的参考信号。这两种技术都能产生高稳定的参考信号,但为了保证选出的两个波长的频率稳定性,需要对光频梳对或者锁模激光器进行复杂的控制,很难具有实用性。而基于半导体激光器的注入锁定技术[7]能够实现宽带频综信号产生,但是由于主光波长和从光波长之间无法严格相位锁定,因此产生的频综信号相位噪声性能较差。与前3种技术不同,基于OEO的频率综合技术采用OEO输出的高频、低相位噪声的微波频率作为参考信号,利用锁相环技术将频率综合器的输出信号和光电振荡频率进行相位锁定,进而实现宽带高性能的频率综合信号产生。基于该机理,Peng等人[8]利用锁相OEO实现了高频率分辨率的X波段频率综合器,输出频率范围覆盖8~10 GHz。该技术避免了传统倍频锁相技术引入的相噪恶化,但其OEO输出频率的相位噪声仅为–100 dBc/Hz@10 kHz,极大地限制了频率综合器性能的提升。此外,Beltchicov等人[9]基于振荡频率为10.2 GHz的OEO及锁相技术也实现了9~18 GHz的频率综合器。该系统通过频率合成输出的10 GHz信号相位噪声约为–120 dBc/Hz@1 kHz以及–140 dBc/Hz@10 kHz,具有很好的相噪性能。值得注意的是,该方案中采用的OEO相位噪声约为–122 dBc/Hz@1 kHz和–151 dBc/Hz@10 kHz,从根本上保证了频率综合器输出信号的频谱纯度及噪声性能。综上所述,基于超低相位噪声OEO的频率综合技术为实现高性能先进频率综合系统提供了有效新思路。

    本文提出了一种基于级联相位调制器(Phase Modulator,PM)的注入锁定OEO及其频率综合系统。与传统的OEO[1012]相比,本文利用相位调制技术展宽光谱并保持光功率恒定,降低振荡环路中光纤引入的受激布里渊散射[13]、自相位调制[14]以及四波混频[15]等非线性影响。采用注入信号和光电振荡信号分别独立调制的方式,隔离注入信号对振荡环路微波器件的反向干扰,从而消除微波链路中反射干扰对振荡信号的影响。双输出马赫曾德尔干涉器(Mach-Zehnder Interferometer,MZI)实现相位调制到强度调制的转化,之后通过平衡探测的方式消除系统共模强度噪声,大幅提升系统信噪比。此外,本文基于所提出的OEO,联合直接数字频率合成(Direct Digital frequency Synthesis, DDS)与锁相环(Phase-Locked Loop, PLL)的混合锁相技术,实现了宽带、高性能频率综合信号的产生。本文中,首先提出并分析基于级联相位调制器的注入锁定OEO的结构、机理及其性能,然后提出基于上述OEO的频率综合系统结构与原理,最后搭建基于级联相位调制器的注入锁定OEO的频率综合实验系统,并对其产生信号的性能开展研究与分析,证实所提OEO及其频综系统的可行性与技术优势。

    图1所示为基于级联相位调制器的注入锁定OEO结构示意图,主要由光源,相位调制器(PM1, PM2),光纤,双输出MZI,平衡探测器(Balanced PhotoDetector, BPD),放大器,窄带带通滤波器(BandPass Filter, BPF),压控移相器(Voltage-Controlled Phase Shifter, VCPS),注入源以及锁相环等各部分组成。PM1将注入源输出的信号相位调制到光源输出的光载波上,输出的相位调制光信号在级联的PM2中相位调制上光电振荡信号,并输出级联相位调制后的光信号。PM2输出的光信号经过长光纤延时后,在双输出MZI中实现相位调制到强度调制的转化,输出两路差分的强度调制光信号,并在BPD中进行平衡探测,实现光信号到电信号的转化。平衡探测后的微波信号经过放大器、带通滤波器及压控移相器后连接至PM2的射频端口形成振荡环路。若注入源信号关闭,当反馈环路满足自由振荡条件后[16],OEO将输出频率为带通滤波器中心频率附近的振荡信号,但是自由振荡的OEO输出频率不固定,且边模抑制较差。

    图  1  基于级联PM的注入锁定OEO结构图
    Figure  1.  Configuration of the injection-locked OEO based on the cascaded phase modulators

    为了使OEO输出固定频率且具有高杂散抑制比,本文将注入源信号打开,调制到PM1上,并通过PM2、光纤、双输出MZI和BPD进入至光电振荡环路,且注入信号频率与带通滤波器的中心频率保持一致。得益于注入信号的牵引与对其它振荡模式的抑制作用,OEO将在注入信号频率处形成振荡,并在带通滤波器输出端口输出高杂散抑制比的振荡信号[10]。此外,为了提高环路信号的信噪比以及降低振荡信号的相位噪声,光源输出功率需要尽可能提高。但是基于传统强度调制的OEO方案,随着光源功率提高,进入到光纤中的强度调制光信号功率也将提升,容易激发光纤中的各类非线性效应[14,17],比如受激布里渊散射、自相位调制及四波混频等。本文在进入长光纤前采用相位调制方式,将光载波的能量分散到各个调制边带上,降低了光纤中受激布里渊散射的影响[17]。同时相位调制不影响输出光信号强度,保持光能量在光纤传播中不随时间变化,从而降低光纤克尔非线性引入的噪声影响。此外,为了提升OEO的频率稳定性,本文将OEO输出频率相位锁定到注入源信号上,即采用PLL技术将OEO的输出频率与注入源信号鉴相并将获得的误差信号控制VCPS的相位变化,从而反馈控制光电振荡环路的相位,最终实现OEO输出频率的稳定。值得注意的是,本文引入的双输出MZI级联平衡探测器的方法拥有两倍于基于传统强度调制OEO的光功率,光电转化后产生的微波功率提升6 dB[18],可大幅提升振荡信号的信噪比。接下来对其进行理论分析,令通过PM1, PM2以及光纤后,输入MZI的光场为

    E(t)=E0exp[jωc(tτ0)+jθ(tτ0)]
    (1)

    其中,E0为进入MZI的光场幅度,ωc为光载波角频率,θ(t)为经过PM1和PM2后相位调制总和,τ0为光纤引入的时延。假设双输出MZI的两臂之间引入的时延差为τd, φd为MZI中两臂引入的相位差,那么MZI两臂输出光场E1(t), E2(t)分别为(忽略MZI本征长度)

    E1(t)=22[ejφdE(tτd)+E(t)]E2(t)=22[ejφdE(tτd)E(t)]}
    (2)

    双输出MZI的两路输出光信号分别注入到BPD中,通过BPD转化成电信号,并输出差分后的光电流i(t),可表示为

    i(t)=2cos[θ(tτ0)θ(tτ0τd)+ωcτd+φd]
    (3)

    其中,为BPD中单个探测器的响应度。由于光电振荡器输出的信号为单频信号,所以可以令θ(t)=βcos(ω0t),其中ω0为振荡信号的角频率,β为小信号调制系数。调节MZI的延时使得ω0τd=(2n1)π,其中n为整数,同时调节MZI的相对相位差使得ωcτd+φd=π/2,最终由式(3)可得(忽略高阶分量)

    i(t)=4βsin(ωtωτ0)
    (4)

    由式(2)、式(3)、式(4)式可得,BPD输出的信号比单路PD输出微波信号功率提升了4倍,从而大幅度提升了系统的信噪比。

    根据图1的原理结构,搭建了基于级联PM的注入锁定OEO,输出信号频率为9.9999914 GHz,其结果如图2所示。图2(a)为OEO输出信号在频率范围为1 MHz、RBW为1 kHz的频谱图。由图2(a)可知,OEO输出信号杂散抑制比大于85 dB,表明该系统具有很低的杂散分量。这些杂散分量主要来自于光电振荡环路中残余的竞争模式。接下来对本文所提的OEO输出频率的相位噪声进行分析,如图2(b)中蓝色曲线所示,系统产生的光电振荡信号在1 kHz频偏处的相位噪声为–127.2 dBc/Hz, 10 kHz频偏处的相位噪声为–153.1 dBc/Hz,振荡模式模式间隔为46 kHz左右。作为对比,红色曲线为商用基于微波技术实现的高性能微波源(E8257D, Keysight)在同频点处的相位噪声曲线,其相位噪声为–105 dBc/Hz@1 kHz以及–114.4 dBc/Hz@10 kHz。对比分析可知,本文所提出的基于级联PM的注入锁定OEO在1 kHz和10 kHz频偏处的相位噪声比E8257D分别低22.2 dB和38.7 dB,表明本系统实现了超低相位噪声信号的产生。

    图  2  所提出OEO输出信号结果图
    Figure  2.  Output results of the proposed OEO

    为了分析光纤非线性效应对OEO相位噪声的影响,本文在相位调制器为小信号调制状态且输出光功率为15.5 dBm情况下,依次在环路光纤前端和后端增加光衰减器,测试不同光衰减量下的相位噪声情况。如图3所示,蓝色曲线和红色曲线分别为在环路光纤的前端和后端增加光衰减器、对应OEO在10 kHz频偏处的相位噪声随光衰减量的变化情况。对比两条曲线可知,当光衰减量从0~5 dB变化时,红色曲线对应的相噪变化较小,而蓝色曲线对应的相位噪声值比红色曲线有了极大的提升,且逐步接近–150 dBc/Hz。表明高功率的调制光信号激发了光纤中的非线性效应,严重地恶化了OEO的相位噪声。通过在光纤前端增加光衰减,降低进入光纤的光功率,可达到降低光纤非线性效应影响、提升OEO性能的效果。当光衰减量从5.0~8.8 dB变化时,蓝色曲线对应的OEO逐渐从最优的相位噪声恶化至–130 dBc/Hz左右,恶化的主要来源为环路衰减引起的信噪比损失。而红色曲线对应的相位噪声仍然比蓝色曲线高近20 dB,表明光纤非线性效应对相位噪声的影响在红色曲线中仍然表现显著。所以,为了实现最优化的相位噪声性能,需要提高系统信噪比同时尽可能降低光纤中非线性效应的影响。

    图  3  光纤非线性效应对OEO相位噪声影响
    Figure  3.  Influence of the nonlinear effects in the optical fiber on the phase noise of the OEO

    目前,实现频率综合器的主要方式包括直接模拟频率合成、DDS、混频法频率合成以及锁相环频率合成等几种途径[19]。结合第2节所提出的基于级联PM的注入锁定OEO,本文利用DDS与锁相环技术实现宽带、高性能频率综合系统,输出频率为5.9~12.9 GHz,其系统结构如图4所示。

    图  4  基于超低相噪OEO的频率综合器系统结构示意图
    Figure  4.  Configuration of the frequency synthesizer based on the ultra-low phase noise OEO

    基于级联PM的注入锁定OEO为本文所提出频率综合系统提供高性能参考信号。如图4所示,OEO输出的9.9999914 GHz信号通过功分器分为A, B, C 3路,其中A路输出信号与压控振荡器(Voltage-Controlled Oscillator, VCO)在混频器1中混频产生0.1~1.1 GHz的中频信号。B路输出信号经过10分频器分频后,产生多个次谐波分频分量,其中1 GHz, 2 GHz, 3 GHz的分频分量分别由对应的带通滤波器选出。频率为3 GHz的分频分量被选出用作DDS的外部参考时钟。DDS的输出频率通过编程控制,可实现0.2~1.2 GHz的频率输出。DDS输出信号与来自混频器1中的中频信号再次混频,输出频率为100 MHz的中频信号,并与C路的光电振荡频率100次分频后的信号一起输入至鉴相鉴频器中(Phase-Frequency Detector, PFD)。PFD输出的误差信号经由带宽为1 MHz的环路滤波器(Loop Filter, LF)提取,并反馈至VCO,进行频率稳定控制,从而实现VCO到OEO的相位锁定。当锁相环路稳定后,VCO可输出频率为8.9~9.9 GHz的信号。然后,通过将VCO输出信号与OEO的10分频信号中1 GHz, 2 GHz或3 GHz的分频分量在混频器3处混频,实现频率为5.9~12.9 GHz的频率输出。

    此外,为了实现宽频谱范围内频率的快速切换,系统引入射频开关(SW1, SW2, SW3及SW4),并通过FPGA对其快速控制。由于频综系统中引入了多个混频器,混频器输出的交调分量将会恶化系统的杂散抑制度,因此在信号输出端引入带宽为500 MHz的带通滤波器组(BPF Group)对杂散信号进行滤除,抑制残余的杂散分量。本系统通过引入高频、超低相位噪声参考源,可实现频率综合器锁相环路中分频系数的大幅降低,避免相位噪声的倍频恶化,从而保证了系统信号的高性能输出。

    基于图4所示结构图,搭建了基于超低相位噪声OEO的频率综合器系统,其中OEO采用基于级联PM的注入锁定技术,输出频率与注入源频率同频(9.9999914 GHz)。OEO中激光器(Emcore 1782)输出波长为1550.2 nm,输出光功率为20 dBm。环路中光纤为单模光纤,类型为G652D,模场直径为9.3 μm。光纤长度约为4.4 km且注入至光纤中的光功率约为 13 dBm。为保证注入信号不对光电振荡器相位噪声造成影响,注入功率低于–15 dBm,对应信号的相位噪声为–118 dBc/Hz @1 kHz以及–122 dBc/Hz @10 kHz。此外,该频率综合系统主要还包括VCO(Hittite-HMC511), DDS(ADI-AD9914), PFD(Hittile-HMC439),分频器,混频器,低噪声放大器,带通滤波器,功分器和射频开关等。实验过程中,采用相噪分析仪(FSWP-50, R&S)对频率综合信号的频谱和相位噪声进行分析,使用实时示波器(Keysight DSOX93304,采样率80 GSa/s)对跳频信号进行采样与分析。

    为了验证以基于级联PM的注入锁定OEO为参考源,锁相后VCO的输出信号性能,调节DDS参数,使其输出599.9914 MHz信号用于VCO的鉴相鉴频,从而使VCO输出9.4 GHz的微波信号,其频谱图如图5(a)所示。由图5(a)可知,VCO成功产生了频率为9.4 GHz的微波信号,且杂散抑制比在10.0 MHz的频率范围内达到74.3 dB。需要注意的是,在9.4 GHz频率两侧频偏大约1 MHz处有两个较高的噪声包络,该噪声包络主要由带宽为1 MHz的锁相环引入。此外,本文还对VCO输出9.4 GHz信号的相位噪声进行了分析,图5(b)给出了锁相前VCO输出相噪曲线(黑线),锁相后VCO输出相噪曲线(红线)以及OEO信号的相噪曲线(蓝线)。由图可知,锁相后输出信号的相位噪声在1 kHz和10 kHz频偏处分别为–115.3 dBc/Hz和–130.7 dBc/Hz,比锁定前信号的相位噪声分别低65.2 dB和47.3 dB,表明VCO输出信号锁相至OEO后相位噪声比锁相前有了极大的提升。值得注意的是,锁相后VCO信号的相位噪声相对于OEO相噪仍然有较大的差距,在1 kHz和10 kHz频偏处比OEO的相位噪声分别高11.9 dB和22.4 dB,其主要原因在于VCO本征噪声、DDS本征噪声和锁相环带内噪声的总噪底比OEO的噪底要高。为进一步降低VCO输出信号的相位噪声,可通过引入更低本征噪声的VCO和DDS,并优化锁相环的有源噪声。

    图  5  VCO输出频率为9.4 GHz信号的结果图
    Figure  5.  Output results of the 9.4 GHz signal generated by the VCO

    改变DDS输出频率,VCO输出频率也会随之变化。图6为VCO输出信号频率为8.9~9.9 GHz、步进为0.1 GHz的相位噪声和杂散抑制情况。由图6(a)可知,VCO输出各频点在10 kHz频偏处的相位噪声值为–130 dBc/Hz左右,表明VCO锁相至超低相噪OEO后实现了宽带、低相噪信号的产生。图6(b)为VCO输出各频点的杂散抑制情况,从图中可知,各频点杂散抑制比都大于70 dB,说明该系统具有高频谱纯度。

    图  6  频率合成输出为8.9~9.9 GHz的情况
    Figure  6.  Frequency synthesyzing performance within 8.9~9.9 GHz

    由于VCO输出的频率范围有限,为进一步拓展频率综合系统的输出频率范围,将VCO的输出信号与OEO的10分频后的次谐波分量(分别为1 GHz, 2 GHz, 3 GHz)进行混频,从而实现了5.9~12.9 GHz的频率合成,其结果如图7所示。图7(a)为本文所提出频率综合系统输出频率为5.9~12.9 GHz、步进为1 GHz情况下的相位噪声曲线。由图7(a)可知,频率综合系统输出信号的相位噪声在5.9~12.9 GHz范围内具有高一致性,且在近频偏处相对于VCO有极大提升。图7(b)分析了该频率综合系统的在不同频率下、10 kHz频偏处的相位噪声值和杂散抑制情况。由图可得,其10 kHz频偏处相位噪声在5.9~12.9 GHz范围内均低于–125 dBc/Hz,最低能到–135 dBc/Hz,表明该系统在更大工作频率范围下仍然能保持低相位噪声的信号输出。值得注意的是,输出信号的杂散抑制比随频率在65~70 dB范围内变化,比图6(b)恶化了5 dB左右,主要原因在于混频器是非线性器件,混频后的信号中包含众多交调分量,恶化了系统杂散性能。此外,还对频率综合系统的跳频时间进行了分析,如图7(c)所示,从9.04 GHz到9.02 GHz的跳频时间低于1.48 μs,表明该系统具有快速跳频性能。

    图  7  频率合成输出为5.9~12.9 GHz的情况
    Figure  7.  Frequency synthesyzing performance within 5.9~12.9 GHz

    本文提出一种基于级联PM的注入锁定OEO及其频率综合系统。该OEO利用相位调制输出光信号具有光谱宽、功率恒定的特点,降低了光纤中的非线性效应引入的强度噪声。还利用双输出MZI级联BPD的结构实现相位调制到强度调制的转化并有效提高了系统的信噪比。为了降低光电振荡信号的边模抑制比,采用在PM1中注入信号的方式,实现了注入锁定OEO,其输出频率为9.9999914 GHz,其相位噪声分别为–127.2 dBc/Hz@1 kHz和–153.1 dBc/Hz@10 kHz,比Keysight E8257D在相同频点、相同频偏处分别低22.2 dB和38.7 dB,具有超低相位噪声性能。此外,本文还基于所提出的超低相位噪声OEO构建了一个宽带频率综合系统。该频率综合系统采用DDS和PLL混合技术,实现了VCO锁相至OEO中,并输出8.9~9.9 GHz信号,且相位噪声保持在–130 dBc/Hz@10 kHz附近,杂散抑制比优于70 dB。最后,本文通过引入电混频器、带通滤波器组及射频开关,将频率综合系统的工作频率进行了拓展,实现了5.9~12.9 GHz的频率输出,对应的10 kHz频偏处相位噪声仍然为–130 dBc/Hz附近,跳频时间小于1.48 μs,且杂散抑制比优于65 dB。值得注意的是,通过改善或者引入更高性能VCO、DDS和锁相环,可进一步提升该频率综合系统的相位噪声性能。综上所述,本文所提出的基于级联PM的注入锁定OEO及其频率综合系统具有宽带、超低相位噪声、低跳频时间等特性,可应用于现代卫星通信、电子战系统、雷达探测以及高端仪器仪表等领域。

  • 图  1  几种典型目标所占距离门数目随信号带宽的变化

    Figure  1.  Variation of the number of range cells for several typical targets with bandwidth

    图  2  ISAR转台模型

    Figure  2.  Turntable model of ISAR

    图  3  含微动部件目标的ISAR几何模型

    Figure  3.  ISAR geometric model of targets with micro-motion components

    图  4  微多普勒信号的STFT谱和频谱图

    Figure  4.  STFT spectrogram and spectrum of the micro-Doppler signal

    图  5  扩展VMD原理示意图

    Figure  5.  Schematic diagram of the extended VMD

    图  6  所提微多普勒抑制方法流程图

    Figure  6.  Flow diagram of the proposed micro-Doppler suppression method

    图  7  多点目标布局及直接成像图

    Figure  7.  Multi-point target layout and direct imaging diagram

    图  8  各种算法处理后的成像结果对比

    Figure  8.  Comparison of the imaging results processed by various algorithms

    图  9  各种算法处理后的STFT谱对比

    Figure  9.  Comparison of the STFT spectrogram processed by various algorithms

    图  10  不同信噪比下的仿真成像结果对比

    Figure  10.  Comparison of simulated imaging results at different SNR

    图  11  成像场景及目标示意图

    Figure  11.  Radar system and target diagram

    图  12  各种算法处理无人机数据的成像结果对比

    Figure  12.  Comparison of the effects of different algorithms on the processing of measured data

    图  13  处理后旋翼中心点对应方位STFT谱

    Figure  13.  STFT spectrogram of the center point on the rotor after processing

    图  14  旋翼中心点残余相位与方位剖面对比图

    Figure  14.  Azimuth residual phase and profile comparison of the center point on the rotor

    表  1  仿真参数表

    Table  1.   Simulation parameter table

    参数 数值
    载波频率fc 35 GHz
    信号带宽Br 10 GHz
    脉冲重复频率PRF 833 Hz
    成像积累时间 2.46 s
    无人机运动速度 5 m/s
    目标俯仰角 15°
    旋翼个数 4
    单旋翼叶片个数 2
    叶片长度l 0.1 m
    单叶片微动散射点个数 5
    旋翼转速ωr 20π rad/s
    下载: 导出CSV

    表  2  各种算法处理后的图像质量比较

    Table  2.   Comparison of image quality processed by various algorithms

    算法 能量相似比(区域1) 能量相似比(区域2) 图像熵 对比度 锐度
    直接成像 8.64 3.87 8.84 18.65 1.03E+10
    PGA 9.60 4.11 8.75 20.19 9.54E+09
    基于EMD 0.22 1.23 8.31 21.44 8.72E+09
    基于LMD 0.16 0.84 8.28 20.45 8.58E+09
    本文方法 0.15 0.16 8.16 22.71 9.15E+09
    下载: 导出CSV

    表  3  六旋翼无人机实测实验参数表

    Table  3.   Experimental parameter table for the hexacopter UAV

    参数 数值
    载波频率fc 35 GHz
    信号带宽Br 10 GHz
    脉冲重复频率PRF 833 Hz
    成像积累时间 0.72 s
    目标俯仰角 15°
    下载: 导出CSV

    表  4  六旋翼无人机成像质量比较

    Table  4.   Comparison of imaging quality of the hexacopter

    算法 图像熵 对比度 锐度
    直接成像 8.19 22.12 5.44E+11
    PGA 7.98 23.59 2.90E+11
    基于EMD 7.63 23.28 4.17E+11
    基于LMD 7.59 23.14 4.11E+11
    基于CVMD 7.67 23.66 4.05E+11
    本文方法 7.05 25.79 4.33E+11
    下载: 导出CSV

    表  5  不同算法的运算复杂度

    Table  5.   Computational complexity of different algorithms

    算法 运算复杂度
    本文算法 O(MNKL1L2log2M)L1为VMD循环次数,L2为差分进化迭代次数
    基于LMD[23] O(MNL3L4)L3为调频信号迭代次数,L4为子信号分离次数
    基于EMD[20] O(MNL5L6)L5为EMD循环次数,L6为复信号投影维度数
    PGA[29] O(ML7log2M+MNL7)L7为PGA迭代次数
    下载: 导出CSV
  • [1] 保铮, 邢孟道, 王彤. 雷达成像技术[M]. 北京: 电子工业出版社, 2005: 239–241.

    BAO Zheng, XING Mengdao, and WANG Tong. Radar Imaging Techniques[M]. Beijing: Publishing House of Electronics Industry, 2005: 239–241.
    [2] 李源. 逆合成孔径雷达理论与对抗[M]. 北京: 国防工业出版社, 2013: 48–55.

    LI Yuan. Inverse Synthetic Aperture Radar Theory and Confrontation[M]. Beijing: National Defense Industry Press, 2013: 48–55.
    [3] 杨建宇. 雷达技术发展规律和宏观趋势分析[J]. 雷达学报, 2012, 1(1): 19–27. doi: 10.3724/SP.J.1300.2013.20010.

    YANG Jianyu. Development laws and macro trends analysis of radar technology[J]. Journal of Radars, 2012, 1(1): 19–27. doi: 10.3724/SP.J.1300.2013.20010.
    [4] 张群, 胡健, 罗迎, 等. 微动目标雷达特征提取、成像与识别研究进展[J]. 雷达学报, 2018, 7(5): 531–547. doi: 10.12000/JR18049.

    ZHANG Qun, HU Jian, LUO Ying, et al. Research progresses in radar feature extraction, imaging, and recognition of target with micro-motions[J]. Journal of Radars, 2018, 7(5): 531–547. doi: 10.12000/JR18049.
    [5] LIU Zheng and SUN Huixia. Micro-Doppler analysis and application of radar targets[C]. IEEE International Conference on Information and Automation, Changsha, China, 2008: 1343–1347. doi: 10.1109/ICINFA.2008.4608210.
    [6] 张群, 罗迎, 何劲. 雷达目标微多普勒效应研究概述[J]. 空军工程大学学报: 自然科学版, 2011, 12(2): 22–26. doi: 10.3969/j.issn.1009-3516.2011.02.005.

    ZHANG Qun, LUO Ying, and HE Jin. Overview of research on micro-Doppler effect of radar targets[J]. Journal of Air Force Engineering University: Natural Science Edition, 2011, 12(2): 22–26. doi: 10.3969/j.issn.1009-3516.2011.02.005.
    [7] CHEN V C. Analysis of radar micro-Doppler with time-frequency transform[C]. The 10th IEEE Workshop on Statistical Signal and Array Processing, Pocono Manor, USA, 2000: 463–466. doi: 10.1109/SSAP.2000.870167.
    [8] WANG Anle, ZHENG Daikun, DU Shirui, et al. Microwave photonic radar system with ultra-flexible frequency-domain tunability[J]. Optics Express, 2021, 29(9): 13887–13898. doi: 10.1364/OE.423952.
    [9] LUO Xiong, WANG Anle, WO Jianghai, et al. Microwave photonic video imaging radar with widely tunable bandwidth for monitoring diverse airspace targets[J]. Optics Communications, 2019, 451: 296–300. doi: 10.1016/j.optcom.2019.06.073.
    [10] CHEN V C, TAHMOUSH D, and MICELI W J. Radar Micro-Doppler Signatures: Processing and Applications[M]. Stevenage: The Institution of Engineering and Technology, 2014: 187–225. doi: 10.1049/pbra034e.
    [11] CHEN V C, LI Fayin, 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. doi: 10.1109/TAES.2006.1603402.
    [12] TUSZYNSKI M, WOJTKIEWICZ A, and KLEMBOWSKI W. Bimodal clutter MTI filter for staggered PRF radars[C]. IEEE International Conference on Radar, Arlington, USA, 1990: 176–180. doi: 10.1109/RADAR.1990.201158.
    [13] 万显荣, 谢德强, 易建新, 等. 基于STFT谱图滑窗相消的微动杂波去除方法[J]. 雷达学报, 2022, 11(5): 794–804. doi: 10.12000/JR22157.

    WAN Xianrong, XIE Deqiang, YI Jianxin, et al. Micro-Doppler clutter removal method based on the cancelation of sliding STFT spectrogram[J]. Journal of Radars, 2022, 11(5): 794–804. doi: 10.12000/JR22157.
    [14] WANG Yong, ZHOU Xingyu, LU Xiaofei, et al. An approach of motion compensation and ISAR imaging for micro-motion targets[J]. Journal of Systems Engineering and Electronics, 2021, 32(1): 68–80. doi: 10.23919/JSEE.2021.000008.
    [15] 何其芳, 张群, 罗迎, 等. 正弦调频Fourier-Bessel变换及其在微动目标特征提取中的应用[J]. 雷达学报, 2018, 7(5): 593–601. doi: 10.12000/JR17069.

    HE Qifang, 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, 2018, 7(5): 593–601. doi: 10.12000/JR17069.
    [16] 符吉祥, 邢孟道, 徐丹, 等. 一种基于微波光子超高分辨雷达机翼振动参数估计方法[J]. 雷达学报, 2019, 8(2): 232–242. doi: 10.12000/JR19001.

    FU Jixiang, XING Mengdao, XU Dan, et al. Vibration-parameters estimation method for airplane wings based on microwave-photonics ultrahigh-resolution radar[J]. Journal of Radars, 2019, 8(2): 232–242. doi: 10.12000/JR19001.
    [17] STANKOVIC L, DJUROVIC I, and THAYAPARAN T. Separation of target rigid body and micro-Doppler effects in ISAR imaging[J]. IEEE Transactions on Aerospace and Electronic Systems, 2006, 42(4): 1496–1506. doi: 10.1109/TAES.2006.314590.
    [18] LI Kaiming, LIANG Xianjiao, ZHANG Qun, et al. Micro-Doppler signature extraction and ISAR imaging for target with micromotion dynamics[J]. IEEE Geoscience and Remote Sensing Letters, 2011, 8(3): 411–415. doi: 10.1109/LGRS.2010.2081660.
    [19] CHOI I, KANG K, KIM K, et al. Use of ICA to separate micro-Doppler signatures in ISAR images of aircraft that has fast-rotating parts[J]. IEEE Transactions on Aerospace and Electronic Systems, 2022, 58(1): 234–246. doi: 10.1109/TAES.2021.3098110.
    [20] BAI Xueru, XING Mengdao, ZHOU Feng, et al. Imaging of micromotion targets with rotating parts based on empirical-mode decomposition[J]. IEEE Transactions on Geoscience and Remote Sensing, 2008, 46(11): 3514–3523. doi: 10.1109/TGRS.2008.2002322.
    [21] FLANDRIN P, RILLING G, and GONCALVES P. Empirical mode decomposition as a filter bank[J]. IEEE Signal Processing Letters, 2004, 11(2): 112–114. doi: 10.1109/LSP.2003.821662.
    [22] GAO Yunchao, GE Guangtao, SHENG Zhengyan, et al. Analysis and solution to the mode mixing phenomenon in EMD[C]. International Congress on Image and Signal Processing, Sanya, China, 2008: 223–227. doi: 10.1109/CISP.2008.193.
    [23] YUAN Bin, CHEN Zengping, and XU Shiyou. Micro-Doppler analysis and separation based on complex local mean decomposition for aircraft with fast-rotating parts in ISAR imaging[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(2): 1285–1298. doi: 10.1109/TGRS.2013.2249588.
    [24] DRAGOMIRETSKIY K and ZOSSO D. Variational mode decomposition[J]. IEEE Transactions on Signal Processing, 2014, 62(3): 531–544. doi: 10.1109/TSP.2013.2288675.
    [25] 杨利超. 超高分辨ISAR成像技术研究[D]. [博士论文], 西安电子科技大学, 2021. doi: 10.27389/d.cnki.gxadu.2021.000089.

    YANG Lichao. Study on ISAR ultrahigh-resolution imaging techniques[D]. [Ph.D. dissertation], Xidian University, 2021. doi: 10.27389/d.cnki.gxadu.2021.000089.
    [26] 邵帅. 高分辨ISAR成像与精细化运动补偿技术研究[D]. [博士论文], 西安电子科技大学, 2020. doi: 10.27389/d.cnki.gxadu.2020.003431.

    SHAO Shuai. Study on high resolution ISAR imaging and fine motion compensation techniques[D]. [Ph.D. dissertation], Xidian University, 2020. doi: 10.27389/d.cnki.gxadu.2020.003431.
    [27] YANG Degui, LI Jin, LIANG Buge, et al. A multi-rotor drone micro-motion parameter estimation method based on CVMD and SVD[J]. Remote Sensing, 2022, 14(14): 3326. doi: 10.3390/rs14143326.
    [28] DAS S and SUGANTHAN P N. Differential evolution: A survey of the state-of-the-art[J]. IEEE Transactions on Evolutionary Computation, 2011, 15(1): 4–31. doi: 10.1109/TEVC.2010.2059031.
    [29] EICHEL P H and JAKOWATZ C V. Phase-gradient algorithm as an optimal estimator of the phase derivative[J]. Optics Letters, 1989, 14(20): 1101–1103. doi: 10.1364/OL.14.001101.
  • 加载中
图(14) / 表(5)
计量
  • 文章访问数: 708
  • HTML全文浏览量: 94
  • PDF下载量: 205
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-03-20
  • 修回日期:  2024-05-24
  • 网络出版日期:  2024-06-19
  • 刊出日期:  2024-08-28

目录

/

返回文章
返回