Study on Efficient High Frequency Method of Electromagnetic Scattering from Ships above Sea Surface
-
摘要: 3维随机粗糙海面与其上方复杂目标复合电磁(EM)散射特性的建模与分析在微波遥感、目标识别、雷达成像、导弹制导等领域中有着重要的研究价值。该文主要研究了基于高频算法的随机粗糙海面及舰船的复合电磁散射特性,开发了PO-IPO混合方法,为3维随机粗糙海面与复杂目标一体化高效求解提供了新思路。文中分别使用了物理光学方法(PO)、迭代物理光学方法(IPO)、PO-PO以及PO-IPO混合方法对海面及舰船进行了建模与仿真,其中,引入锥形波来代替平面波作为发射源,锥形波可以更好地抑制粗糙面在边缘位置被突然截断而形成的电磁反射和边缘绕射等效应。从数值仿真结果中可以看出,PO-IPO混合方法可将复杂物体本身面元间以及粗糙海面与物体间的耦合作用考虑在内,因此PO-IPO可以作为一种有效的途径来快速获取随机粗糙海面及舰船的复合电磁散射特性。Abstract: The study on the composite ElectroMagnetic (EM) scattering characteristics of electrically large ships over a random rough surface is significant in the fields of microwave remote sensing, target recognition, and radar imaging. In this paper, an efficient PO-IPO hybrid method is proposed to analyze the scattering from complex targets over the sea surface. The PO method can be used to analyze the scattering from an electrically large sea surface. However, low accuracy will be achieved for targets with complex structures. Therefore, the IPO is applied for computation herein. In the numerical results, several high-frequency methods (PO, IPO, PO-PO, and PO-IPO) are used to analyze the scattering characteristics of targets over the sea surface. Notably, the tapered wave is used as the source in the numerical results, which can restrain the EM reflection and edge diffraction. The results also indicate that the proposed PO-IPO can be used as an efficient tool to obtain the scattering characteristics of complex targets over the sea surface.
-
Key words:
- Random rough surface /
- Ship /
- EM scattering /
- Tapered wave
-
1. 引言
作为现代通信、探测和测试系统中的核心部件,频率综合器能够在一定带宽范围内灵活输出高精度、高稳定的频率信号,广泛应用于无线通信、雷达探测、测试与计量、卫星等领域[1–4]。随着技术的发展,高频率、大带宽、低相噪及低杂散成为新一代高性能频率综合系统的必然演进方向。传统基于微波技术的频率综合技术采用低频晶体振荡器作为参考源,通过倍频锁相的方式产生高频信号,但是所产生的高频信号的相位噪声在倍频锁相过程中被快速恶化(20lgN, N为倍频因子),难以进一步提升。
相比之下,光子技术具有低噪声、低损耗和大带宽等关键优势,因此通过光子技术实现频率综合已获得人们的广泛重视[5–9]。光域实现微波频率综合技术主要包括双光频梳锁定技术[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[10–12]相比,本文利用相位调制技术展宽光谱并保持光功率恒定,降低振荡环路中光纤引入的受激布里渊散射[13]、自相位调制[14]以及四波混频[15]等非线性影响。采用注入信号和光电振荡信号分别独立调制的方式,隔离注入信号对振荡环路微波器件的反向干扰,从而消除微波链路中反射干扰对振荡信号的影响。双输出马赫曾德尔干涉器(Mach-Zehnder Interferometer,MZI)实现相位调制到强度调制的转化,之后通过平衡探测的方式消除系统共模强度噪声,大幅提升系统信噪比。此外,本文基于所提出的OEO,联合直接数字频率合成(Direct Digital frequency Synthesis, DDS)与锁相环(Phase-Locked Loop, PLL)的混合锁相技术,实现了宽带、高性能频率综合信号的产生。本文中,首先提出并分析基于级联相位调制器的注入锁定OEO的结构、机理及其性能,然后提出基于上述OEO的频率综合系统结构与原理,最后搭建基于级联相位调制器的注入锁定OEO的频率综合实验系统,并对其产生信号的性能开展研究与分析,证实所提OEO及其频综系统的可行性与技术优势。
2. 基于级联PM的注入锁定光电振荡器
图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输出频率不固定,且边模抑制较差。
为了使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)=2ℜcos[θ(t−τ0)−θ(t−τ0−τd)+ωcτd+φd] (3) 其中,
ℜ 为BPD中单个探测器的响应度。由于光电振荡器输出的信号为单频信号,所以可以令θ(t)=βcos(ω0t) ,其中ω0 为振荡信号的角频率,β 为小信号调制系数。调节MZI的延时使得ω0τd=(2n−1)π ,其中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,表明本系统实现了超低相位噪声信号的产生。
为了分析光纤非线性效应对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的频率综合系统
目前,实现频率综合器的主要方式包括直接模拟频率合成、DDS、混频法频率合成以及锁相环频率合成等几种途径[19]。结合第2节所提出的基于级联PM的注入锁定OEO,本文利用DDS与锁相环技术实现宽带、高性能频率综合系统,输出频率为5.9~12.9 GHz,其系统结构如图4所示。
基于级联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. 实验及结果分析
基于图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,并优化锁相环的有源噪声。
改变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,说明该系统具有高频谱纯度。
由于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,表明该系统具有快速跳频性能。
5. 结束语
本文提出一种基于级联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 PO与JMCFIE计算资源对比
Table 1. Comparison of computational resources between PO and JMCFIE
数值方法 时间(s) 内存(MB) SBR 27 176 PO 23 174 JMCFIE 1956 5912 表 2 SBR, IPO与MoM-MLFMA计算资源对比
Table 2. Comparison of computational resources among SBR, IPO and MoM-MLFMA
数值方法 时间(s) 内存(MB) SBR 15 58 IPO 2593 74 MoM-MLFMA 112570 1178 表 3 PO-IPO与JMCFIE计算资源对比
Table 3. Comparison of computational resources between PO-IPO and JMCFIE
数值方法 时间(s) 内存(MB) PO-IPO 471.8 42.5 JMCFIE 8994.0 21645.0 -
[1] THORSOS E I. The validity of the Kirchhoff approximation for rough surface scattering using a Gaussian roughness spectrum[J]. The Journal of the Acoustical Society of America, 1988, 83(1): 78–92. doi: 10.1121/1.396188 [2] LEADER J. The relationship between the Kirchhoff approach and small perturbation analysis in rough surface scattering theory[J]. IEEE Transactions on Antennas and Propagation, 1971, 19(6): 786–788. doi: 10.1109/TAP.1971.1140044 [3] BASS F G and FUKS I M. Wave Scattering from Statistically Rough Surfaces[M]. Oxford: Pergamon Press, 1979. [4] FUKS I M. Wave diffraction by a rough boundary of an arbitrary plane-layered medium[J]. IEEE Transactions on Antennas and Propagation, 2001, 49(4): 630–639. doi: 10.1109/8.923325 [5] 魏鹏博, 江旺强, 张晓梨, 等. 基于小斜率近似方法的沥青路面电磁散射[C]. 第十二届全国电波传播学术讨论会论文集, 青岛, 2013.WEI Pengbo, JIANG Wangqiang, ZHANG Xiaoli, et al. Electromagnetic L scattering of asphalt pavement based on small slope approximation method[C]. National Symposium on Radio Communication, Qingdao, 2013. [6] JOHNSON J T, SHIN R T, KONG Jinau, et al. A numerical study of the composite surface model for ocean backscattering[J]. IEEE Transactions on Geoscience and Remote Sensing, 1998, 36(1): 72–83. doi: 10.1109/36.655319 [7] CHEN M F and FUNG A K. A numerical study of the regions of validity of the Kirchhoff and small-perturbation rough surface scattering models[J]. Radio Science, 1988, 23(2): 163–170. doi: 10.1029/RS023i002p00163 [8] THORSOS E I and JACKSON D R. The validity of the perturbation approximation for rough surface scattering using a Gaussian roughness spectrum[J]. The Journal of the Acoustical Society of America, 1989, 86(1): 261–277. doi: 10.1121/1.398342 [9] JOHNSON J T. A numerical study of scattering from an object above a rough surface[J]. IEEE Transactions on Antennas and Propagation, 2002, 50(10): 1361–1367. doi: 10.1109/TAP.2002.802152 [10] VORONOVICH A. Small-slope approximation for electromagnetic wave scattering at a rough interface of two dielectric half-spaces[J]. Waves in Random Media, 1994, 4(3): 337–367. doi: 10.1088/0959-7174/4/3/008 [11] NIETO-VESPERINAS M. Depolarization of electromagnetic waves scattered from slightly rough random surfaces: A study by means of the extinction theorem[J]. Journal of the Optical Society of America, 1982, 72(5): 539–547. doi: 10.1364/JOSA.72.000539 [12] BAHAR E and LEE B S. Radar scatter cross sections for two-dimensional random rough surfaces - full wave solutions and comparisons with experiments[J]. Waves in Random Media, 1996, 6(1): 1–23. doi: 10.1080/13616679609409792 [13] WINEBRENNER D P and ISHIMARU A. Application of the phase-perturbation technique to randomly rough surfaces[J]. Journal of the Optical Society of America A, 1985, 2(12): 2285–2294. doi: 10.1364/JOSAA.2.002285 [14] LIANG Y, ZENG X H, GUO L X, et al. Review on the pre-study about specific target and rough soil surface composite electromagnetic scattering basing on an efficient numerical strategy[C]. 2016 Progress in Electromagnetic Research Symposium (PIERS), Shanghai, China, 2016: 3198. [15] 叶红霞. 随机粗糙面与目标复合电磁散射的数值计算方法[D]. [博士论文], 复旦大学, 2007.YE Hongxia. Numerical methods for the composite electromagnetic scattering of random rough surface and targets[D]. [Ph.D. dissertation], Fudan University, 2007. [16] 朱国强, 孙劲, 郑立志, 等. 平板目标与随机粗糙面对电磁波的复合散射[J]. 武汉大学学报(自然科学版), 2000, 46(1): 99–103.ZHU Guoqiang, SUN Jin, ZHENG Lizhi, et al. Composite scattering of electromagnetic waves by the plate target and random rough surface[J]. Journal of Wuhan University (Natural Science Edition) , 2000, 46(1): 99–103. [17] 杨伟. 三维复杂粗糙海面电磁散射建模研究与特性分析[D]. [博士论文], 电子科技大学, 2012.YANG Wei. Charecteristics analysis and modeling research on electromagnetic scattering from 3D complicated rough sea surface[D]. [Ph.D. dissertation], University of Electronic Science and Technology of China, 2012. [18] 王强, 郭立新, 刘忠玉. 粗糙海面与上方多目标瞬态电磁散射混合算法[J]. 西安电子科技大学学报, 2018, 45(5): 57–63. doi: 10.3969/j.issn.1001-2400.2018.05.010WANG Qiang, GUO Lixin, and LIU Zhongyu. Hybrid method for transient electromagnetic scattering from multiple targets above the rough sea surface[J]. Journal of Xidian University, 2018, 45(5): 57–63. doi: 10.3969/j.issn.1001-2400.2018.05.010 [19] 叶红霞, 金亚秋. 三维介质粗糙面上导体目标散射的解析-数值混合算法[J]. 电波科学学报, 2008, 23(6): 1144–1153, 1187. doi: 10.3969/j.issn.1005-0388.2008.06.025YE Hongxia and JIN Yaqiu. A hybrid analytical-numerical algorithm for computation of scattering from a 3D PEC target above a dielectric rough surface[J]. Chinese Journal of Radio Science, 2008, 23(6): 1144–1153, 1187. doi: 10.3969/j.issn.1005-0388.2008.06.025 [20] 张民, 郭立新, 聂丁, 等. 海面目标雷达散射特性与电磁成像[M]. 北京: 科学出版社, 2015.ZHANG Min, GUO Lixin, NIE Ding, et al. Radar Scattering Characteristics and Electromagnetic Imaging of Targets on the Sea Surface[M]. Beijing: Science Press, 2015. [21] 金亚秋, 刘鹏, 叶红霞. 随机粗糙面与目标复合散射数值模拟理论与方法[M]. 北京: 科学出版社, 2008.JIN Yaqiu, LIU Peng, and YE Hongxia. Theory and Method of Numerical Simulation of Composite Scattering from the Object and Randomly Rough Surface[M]. Beijing: Science Press, 2008. [22] 郭立新, 王蕊, 吴振森. 随机粗糙面散射的基本理论与方法[M]. 北京: 科学出版社, 2010.GUO Lixin, WANG Rui, and WU Zhensen. Basics Theory and Method for the Scattering from the Random Rough Surface[M]. Beijing: Science Press, 2010. [23] 康士峰, 王显德. 粗糙面与目标电磁散射统计特性分析[J]. 微波学报, 2004, 20(3): 43–46. doi: 10.3969/j.issn.1005-6122.2004.03.010KANG Shifeng and WANG Xiande. Analysis of statistical characteristics for EM scattering from rough surface and target[J]. Journal of Microwaves, 2004, 20(3): 43–46. doi: 10.3969/j.issn.1005-6122.2004.03.010 [24] ZHU G Q, CAO Q F, and YANG H L. Analysis of double bounce scattering between target and rough surface[J]. Chinese Journal of Radio Science, 1996, 11(1): 80–84. [25] 邹高翔, 童创明, 高飞, 等. 复杂陆地粗糙面及其上方坦克目标复合散射研究[J]. 电波科学学报, 2017, 32(3): 261–272. doi: 10.13443/j.cjors.2017042801ZOU Gaoxiang, TONG Chuangming, GAO Fei, et al. Composite electromagnetic scattering from tank target above complicated ground rough surface[J]. Chinese Journal of Radio Science, 2017, 32(3): 261–272. doi: 10.13443/j.cjors.2017042801 [26] 潘小敏, 盛新庆, 孔雷. 电特大舰船电磁环境的并行多层快速多极子研究[J]. 装备环境工程, 2008, 5(1): 48–52, 60. doi: 10.3969/j.issn.1672-9242.2008.01.012PAN Xiaomin, SHENG Xinqing, and KONG Lei. Analysis on electromagnetic environment of an extremely large ship by parallel MLFMA[J]. Equipment Environmental Engineering, 2008, 5(1): 48–52, 60. doi: 10.3969/j.issn.1672-9242.2008.01.012 [27] 刁桂杰, 许小剑, 倪虹, 等. 舰船目标宽带单脉冲雷达三维成像建模与仿真[J]. 系统仿真学报, 2018, 30(7): 2515–2524. doi: 10.16182/j.issn1004731/x.joss.201807011DIAO Guijie, XU Xiaojian, NI Hong, et al. Three-dimensional imaging modeling and simulation of ships based on wide-band monopulse radar[J]. Journal of System Simulation, 2018, 30(7): 2515–2524. doi: 10.16182/j.issn1004731/x.joss.201807011 [28] YE H X and JIN Y Q. Fast iterative approach to electromagnetic scattering from the target above a rough surface[J]. IEEE Transactions on Geoscience and Remote Sensing, 2006, 44(1): 108–115. [29] DING D Z, LI G M, AN Y Y, et al. Application of hierarchical two-level spectral preconditioning method for electromagnetic scattering from the rough surface[J]. International Journal of Antennas and Propagation, 2014, 2014: 1–10. doi: 10.1155/2014/752418 [30] AN Yuyuan, WANG Daoxiang, and CHEN Rushan. A fast numerical algorithm for calculating electromagnetic scattering from an object above a rough surface[J]. Electromagnetics, 2013, 33(1): 10–22. doi: 10.1080/02726343.2013.751006 -