基于极化时变调控表面的有源欺骗干扰辨识方法

陈焱 王占领 庞晨 李永祯 王壮

薛策文, 冯晅, 李晓天, 等. 全极化探地雷达多极化数据融合分析研究[J]. 雷达学报, 2021, 10(1): 74–85. doi: 10.12000/JR20104
引用本文: 陈焱, 王占领, 庞晨, 等. 基于极化时变调控表面的有源欺骗干扰辨识方法[J]. 雷达学报(中英文), 2024, 13(4): 929–940. doi: 10.12000/JR24028
XUE Cewen, FENG Xuan, LI Xiaotian, et al. Multi-polarization data fusion analysis of full-polarimetric ground penetrating radar[J]. Journal of Radars, 2021, 10(1): 74–85. doi: 10.12000/JR20104
Citation: CHEN Yan, WANG Zhanling, PANG Chen, et al. Radar active deception jamming recognition method based on the time-varying polarization-conversion metasurface[J]. Journal of Radars, 2024, 13(4): 929–940. doi: 10.12000/JR24028

基于极化时变调控表面的有源欺骗干扰辨识方法

DOI: 10.12000/JR24028
基金项目: 国家自然科学基金(61971429, 61921001, 62301580),博士后基金(2022M723917)
详细信息
    作者简介:

    陈 焱,博士生,主要研究方向为极化调控技术、雷达对抗技术

    王占领,博士,助理研究员,主要研究方向为极化阵列雷达及对抗技术

    庞 晨,博士,副研究员,主要研究方向为雷达极化信息处理、雷达抗干扰与识别技术

    李永祯,博士,研究员,主要研究方向为雷达极化信息处理、空间电子对抗、目标检测与识别

    王 壮,博士,教授,主要研究方向为雷达信息处理、空间目标监视、自动目标识别

    通讯作者:

    王占领 wangzhanling17@nudt.edu.cn

  • 责任主编:全英汇 Corresponding Editor: QUAN Yinghui
  • 中图分类号: TN957

Radar Active Deception Jamming Recognition Method Based on the Time-varying Polarization-conversion Metasurface

Funds: The National Natural Science Foundation of China (61971429, 61921001, 62301580), China Postdoctoral Science Foundation (2022M723917)
More Information
  • 摘要: 聚焦雷达对抗中极化信息获取与利用的应用需求,该文研究了基于极化时变调控表面的有源欺骗干扰辨识方法。首先,设计了一套在9.6~10.1 GHz频带内支持3 bit相位量化的各向异性相位调制表面,通过优化相位调制编码序列,实现了极化态按需调控。然后,将极化调控表面加装在单极化雷达天线上,使天线发射和接收电磁波的极化态沿特定极化轨道变化,通过提取目标与有源欺骗干扰的极化域差异,实现两者辨识。仿真分析表明,在3种不同的极化轨道约束下,干扰与目标均具有显著的聚类效应,可获得稳定的干扰辨识效果。相较于依赖双极化或全极化雷达体制的干扰辨识方法,该文所提方法兼具低成本与高效性,在雷达抗干扰中具有很大的应用潜力。

     

  • 海洋涡旋是一种旋转的、以封闭环流为主要特征的水体,是由于各种气象因素作用和海洋动力不稳定性形成的。作为一种重要的海洋现象,涡旋不仅能够影响海洋流场与化学物质的输送,从而对海洋的环流结构和海洋生态等产生重要作用,还能通过海气相互作用,对风场、云及降雨等大气现象产生影响[1,2]

    合成孔径雷达(Synthetic Aperture Radar, SAR)具有全天时、全天候、高分辨率、广覆盖面等优点,对海洋涡旋探测具有特殊意义,受到国际海洋遥感界的重视。然而,涡旋在SAR成像时会受到各种海洋环境因素的影响,通过真实SAR图像难以完全解译涡旋的特征。利用仿真SAR图像可以为涡旋的SAR图像特征解译提供指导,但是目前利用SAR图像对涡旋的研究主要集中在涡旋的统计性研究[35]、涡旋的形成机制和成像分类[68]以及涡旋的检测和特征提取方面的研究[912],极少有关于涡旋SAR图像仿真方法的研究。

    由于海面随机运动且电磁散射特性复杂,难以进行时间和空间上的SAR原始回波仿真。海浪谱能够描述随机海面不同波长海浪的能量分布情况,因此利用海浪谱可以很好地描述不同海况下随机海面的统计特征。海洋涡旋、内波、浅海地形、锋面等都可看作是通过波流交互作用,即利用自身流场改变海浪谱分布,并经过海面电磁散射模型,进而得以在SAR图像上体现。目前SAR海面图像仿真常用的电磁散射模型包括Kirchhoff散射模型、Bragg散射模型以及组合表面散射模型[13]。这些模型只考虑了1阶Bragg散射,仅适用于低频(小于L波段)SAR海面图像仿真。1997年,Romeiser和Alpers[14,15]提出了改进的组合表面模型,该模型考虑了2阶Bragg波散射的影响,从而使仿真的SAR海面图像更接近实际情况。2002年,Romeiser[16]利用该模型研究了浅海地形在SAR图像上的特征,并与声学多普勒流速剖面仪测量的浅海地形进行对比,验证了该模型用于SAR海面图像仿真的合理性;2011年,欧阳越等[17]利用该模型仿真了不同雷达参数下海洋内波图像,并同实际内波SAR图像进行对比,发现二者具有较高的一致性。但是目前,利用海面电磁散射模型对海洋涡旋SAR图像仿真的研究尚未见报道。

    为此,本文提出了一种海洋涡旋SAR图像仿真方法,利用流体力学中典型的Burgers-Rott涡旋模型,建立涡旋的2维流场。利用SAR海洋成像仿真模型,仿真涡旋SAR图像。基于此方法,本文进行了气旋式涡旋与反气旋式涡旋SAR图像仿真实验,并将仿真SAR图像与ERS-2 SAR图像和ENVISAT-1 ASAR图像进行对比,从而验证该方法的有效性。

    本文建立的涡旋SAR图像仿真方法,是在给定2维涡旋流场和风场条件下,利用SAR海洋成像模型生成随机海面的2维海浪谱,再根据2维海浪谱与SAR图像之间的调制传递函数,生成仿真涡旋SAR图像。

    涡旋SAR图像仿真方法分为两步,如图1所示。首先,输入涡旋流场参数,基于涡旋动力学模型建立涡旋2维流场(于2.1节介绍)。然后,将仿真的涡旋流场和海面风场输入到SAR海洋成像仿真模型,通过设置SAR参数获得仿真涡旋SAR图像(于2.2节介绍)。

    图  1  涡旋SAR图像仿真方法流程图
    Figure  1.  Flow chart of the simulation method of SAR eddy image

    涡旋一般遵循流体力学的纳维-斯托克斯(Navier-Stokes,简写N-S)方程,根据方程中黏性力项、惯性力项以及离心力项的平衡关系,可以建立不同的涡旋模型。常见的涡旋模型包括Rankine涡旋、Oseen涡旋、Sullivan涡旋以及Burgers-Rott涡旋[1820],其中,Rankine涡旋模型没有考虑N-S方程中的黏性力项,流体以常角速度ω旋转,没有径向速度,因而不能产生涡旋的辐散、辐聚和上升运动;Oseen涡旋模型仅考虑N-S方程中惯性力项的局地项及黏性力项,其轨道是一个圆形涡旋,不符合实际SAR图像中涡旋的形态;Sullivan涡旋模型和Burgers-Rott涡旋模型考虑了N-S方程中全部的黏性力项、惯性力项及离心力项,但由于Sullivan涡旋模型的轨道是一个双螺旋涡旋,Burgers-Rott涡旋模型的轨道是一个螺旋形涡旋,后者与真实SAR图像所呈现的涡旋形状更为接近,因此本文选用Burgers-Rott涡旋模型来建立海洋涡旋的流场。

    Burgers-Rott涡旋模型是从N-S方程求得的一个涡旋解[19,20],假定涡旋是定常和轴对称的,涡旋速度场在柱坐标系下表示为

    {Vr=drdt=α2rVθ=rdθdt=Γ02πr(1eαr24υ)Vz=dzdt=αz
    (1)

    其中,Vr,Vθ,Vz分别是r,θ,z方向的速度分量,α为吸入强度,υ为黏性系数,Γ0r时的速度环量,Γ=2πrvθ

    将式(1)转化为直角坐标系,涡旋速度场可表示为

    {Vx=α2xΓ0α8πυyVy=Γ0α8πυxα2y
    (2)

    其中,Vx为涡旋速度场在x方向上的速度分量,Vy是涡旋速度场在y方向上的速度分量。

    通过设置参数α Γ0/υ的值,根据式(2)可以得到涡旋2维流场。通过仿真发现,α的值会影响涡旋流场流速的大小,α的值越大,涡旋流场流速越大,反之则越小;α的正负影响涡旋流场的旋向,α为正,流场顺时针旋转,α为负,流场逆时针旋转;Γ0/υ的值则会影响涡旋臂的曲率,Γ0/υ的值越大,涡旋臂的曲率越大。

    获得了涡旋的流场之后,下一步将进行涡旋SAR图像的仿真。本文使用SAR海洋成像仿真模型来仿真涡旋SAR图像。SAR海洋成像仿真模型主要分为波流交互作用模型、雷达后向散射模型和SAR成像模型3个部分,如图2所示。

    图  2  SAR海洋成像仿真模型示意图[22]
    Figure  2.  Schematic diagram of oceanic SAR imagery simulation model[22]

    首先,将仿真的涡旋2维流场和海面风场输入到波流交互作用模型,通过求解作用量谱平衡方程,计算给定海面流场和海面风场下被调制的海浪谱。作用量谱平衡方程如式(3)所示[21]

    N(x,k,t)t+[cg(k)+U(x,t)]N(x,k,t)xkU(x,t)xN(x,k,t)k=S(x,k,t)
    (3)

    其中,N为微尺度波作用量谱密度,x=(x,y)为空间位置矢量,k=(kx,ky)为波数矢量,U为表面流场,cg为被调制波浪的群速度,S为源函数(风场输入、非线性波-波作用和弥散等作用之和),在本文模型中采用的源函数表达式为

    S(x,k,t)=μ(k)N(x,k,t)(1N(x,k,t)N0(k))
    (4)

    其中,N0为不存在海流时平衡状态下的作用量谱密度,μ为松弛率。

    作用量谱密度与海浪谱的关系为[23]

    N(x,k,t)=ρω0(k)kψ(x,k,t)
    (5)

    其中,ω0(k)=gk+(τ/ρ)k3, τ为表面张力,ρ为海水密度,ψ为海浪谱。

    Q(x,k,t)=1/N(x,k,t), Q0(k)=1/N0(k),则被流场调制后的海浪谱为

    ψ(x,k,t)ψ0(k)=Q0(k)Q0(k)+δQ(x,k,t)=11+δQ(x,k,t)Q0(k)
    (6)

    其中,δQ表示调制引起的作用量谱变化量。

    然后,将计算得到的海浪谱输入到雷达后向散射模型,在给定雷达频率、入射角、极化方式及雷达视向等雷达参数下,仿真涡旋SAR图像后向散射强度。本文采用的雷达后向散射模型为改进的组合表面模型,是Romeiser和Alpers等在Bragg共振散射模型基础上的改进[14,15]。该模型同时考虑了长波和中波对短波的倾斜调制和水动力调制,所以从理论和试验研究上更能表现海面微波散射的实际情况,是目前最为完善的海面微波散射模型之一。该模型是基于2维海面坡度,通过傅里叶变换对后向散射截面进行泰勒级数展开,并对后向散射截面进行时间和空间上的平均。由于1阶项平均后为0,因此得到2阶Bragg散射后的海面归一化后向散射系数为[14,22]

    σ=σ(0)+σ(2)=σ|s=0+(2σspsp|s=0+2σspsp|s=0)k2pψ(k)d2k+(2σsnsn|s=0+2σsnsn|s=0)k2nψ(k)d2k+(2σspsn|s=0+2σsnsp|s=0+2σspsn|s=0+2σsnsp|s=0)kpknψ(k)d2k
    (7)

    其中,σ(0)为平静海面的归一化后向散射系数;σ(2)表示表面坡度引起的2阶Bragg散射之和;符号表示统计平均;s=(sp,sn)为海面坡度;kp, kn分别为平行和垂直于雷达视向的Bragg波波数分量;ψ(k)为海浪波数谱;符号分别表示σ对波数k的傅里叶变换及其共轭;表示σ对组合波数k1+k2的傅里叶变换及其共轭;表示σ对组合波数k1k2的傅里叶变换及其共轭。

    上述过程中,利用海浪谱与雷达后向散射模型得到仿真的海面归一化后向散射系数,但这是一个实孔径雷达成像过程,SAR图像仿真还需考虑海面运动的影响。当目标存在沿雷达视线方向的径向速度时,将在方位向上产生偏移Δx

    Δx=RV vr=Rλ2V fD
    (8)

    其中,R是雷达至目标的距离,V是平台飞行速度,λ是雷达波长,fD=2vrλ是目标速度导致的Doppler谱中心偏移。对于海面而言,由于其各点速度不同,在方位向上偏移量不同,导致SAR海面图像产生压缩或拉伸的现象,即速度聚束效应[24]。此外,分辨单元内不同散射点速度的分布方差将造成回波Doppler谱展宽,并导致分辨率下降。SAR成像模型通过计算每个分辨单元的平均Doppler谱中心和方差引入海面运动对SAR成像造成的影响。

    这里采用Romeiser和Thompson[25]给出的双高斯形Doppler谱模型计算Doppler谱中心和方差,该模型将海面回波Doppler谱分成朝向雷达和远离雷达两个传播方向的Bragg波Doppler谱的叠加,每个Doppler谱分量为高斯形,其具体表达式为

    W(fD)=σ+2πγ2D+e(fDfD+σ)2/γ2D++σ2πγ2De(fDfDσ)2/γ2D
    (9)

    其中,±表示远离雷达方向和朝向雷达方向的两组Bragg波分量,fD±σ表示经过归一化后向散射系数σ加权的平均Doppler中心;γD±表示Doppler谱的方差。fD±σγD±的具体计算过程可以参考文献[25],这里不再赘述。

    另外,仿真的SAR图像还需考虑噪声的影响,本文涡旋SAR图像仿真过程中,仅考虑热噪声对仿真SAR图像信噪比的影响。信噪比由噪声等效后向散射系数以及海面归一化后向散射系数所决定:

    SNR(dB)=σNEσ0
    (10)

    其中,海面归一化后向散射系数σ由入射角、雷达频率、极化方式、海面风速等参数所决定,NEσ0为噪声等效后向散射系数,由系统硬件参数所决定。因此,SAR成像模型根据给定的仿真输入参数计算信噪比,从而得到具有统计特性的仿真涡旋SAR图像。

    根据涡旋旋转方向的不同,可将涡旋分为气旋式涡旋与反气旋式涡旋[26]。气旋式涡旋在北半球逆时针旋转,在南半球顺时针旋转;反气旋式涡旋在北半球顺时针旋转,在南半球逆时针旋转。不同旋转方向的涡旋将产生不同的涡旋流场,从而在SAR图像中呈现不同的涡旋特征。下面,本文分别针对气旋式涡旋与反气旋式涡旋进行仿真实验。

    图3是一幅ERS-2 SAR图像,图像获取时间为2009.08.19, 02:23:50 UTC,获取地点为中国东海海域。图中方框1处为一个气旋式涡旋,旋转方向为逆时针。为了便于对比仿真SAR图像与真实SAR图像,将方框1处的涡旋截取出来,截取图像尺寸为18 km×24 km,如图4所示。ERS-2 SAR图像的具体雷达参数如表1所示。

    图  3  中国东海海域获取的ERS-2 SAR图像,获取时间为2009.08.19, 02:23:50 UTC
    Figure  3.  ERS-2 SAR image of the East China Sea obtained on August 19, 2009 at 02:23:50 UTC
    图  4  从方框1处截取的涡旋SAR图像
    Figure  4.  Enlargement of the eddy in Frame 1
    表  1  ERS-2 SAR参数
    Table  1.  SAR parameters of ERS-2
    参数数值
    极化方式VV
    波段C
    入射角23.0°
    平台高度780 km
    平台速度7500 m/s
    下载: 导出CSV 
    | 显示表格

    从欧洲中期天气预报中心(Europe Centre for Medium-Range Weather Forecasts, ECMWF)获取2009.08.19, 03:00:00时刻的风场再分析资料,分辨率为0.125°×0.125°。根据数据显示,涡旋区域附近的风速为1.4 m/s,风向为257.9°。从全球海洋数据同化系统(Global Ocean Data Assimilation System, GODAS)获取相同位置的5日平均流场再分析资料,分辨率为(1/3)°×1°。根据数据显示,涡旋区域附近的流速为0.61 m/s。因此,设置参数α为–0.003486,流场大小设置为18 km×24 km,空间分辨率为100 m,雷达参数设置为表1中ERS-2SAR参数。

    图5(a)图5(b)分别是该涡旋的仿真SAR图像与获取的真实SAR图像,仿真时设定的雷达参数、海面风场条件与真实SAR图像获取条件完全一致。对比图5(a)图5(b)两图发现,仿真SAR图像与真实SAR图像中的涡旋臂形状几乎一致,涡旋臂的亮暗特征也基本吻合。从逆时针方向看,涡旋臂由外到内的亮暗特征均为亮-暗-亮,这种亮暗特征的变化是由雷达后向散射引起的布拉格波谱密度变化导致的[7]。这初步验证了仿真方法的正确性。

    图  5  相同参数下仿真SAR图像与ERS-2 SAR图像对比图
    Figure  5.  Comparison of simulated SAR image and ERS-2 SAR image under the same parameters

    为了进一步验证仿真方法的正确性,定量地描述仿真SAR图像与真实SAR图像的中涡旋的相似程度,采用文献[9]中基于对数螺旋线边缘拟合的SAR图像涡旋信息提取方法,提取仿真SAR图像和真实SAR图像中涡旋的中心位置、直径及边缘长度,并加以比较。拟合及提取结果如图6所示,红色加号表示涡旋中心位置,黄色箭头表示涡旋直径,蓝色曲线表示涡旋边缘,具体数值如表2所示。

    图  6  仿真SAR图像与ERS-2 SAR图像涡旋信息提取
    Figure  6.  Eddy information extraction of simulated SAR image and ERS-2 SAR image
    表  2  涡旋信息提取结果
    Table  2.  Results of eddy information extraction
    SAR图像涡旋中心位置涡旋直径涡旋边缘长度
    仿真SAR图像(116,75)18.9 km35.7 km
    真实SAR图像(113,71)18.7 km35.4 km
    绝对/相对误差(3,4)/—0.2 km/0.0110.3 km/0.008
    下载: 导出CSV 
    | 显示表格

    对比仿真SAR图像与真实SAR图像的涡旋信息提取结果,可以发现两幅图像中涡旋的中心位置较为一致,方位向和距离向上仅相差3~4个像素点,涡旋直径及边缘长度的相对误差均不超过0.011,证明本文提出的基于Burgers-Rott涡旋模型的涡旋SAR图像仿真方法能够实现气旋式涡旋的SAR图像仿真,并且仿真SAR图像与真实SAR图像能够较好地吻合。

    3.1节对气旋式涡旋进行了仿真实验,本节将针对反气旋式涡旋进行仿真实验。图7是一幅ENVISAT-1 ASAR图像,图像获取时间为2010.06.11, 01:51:48 UTC,获取地点在吕宋海峡。图中方框2处为一个反气旋式涡旋,旋转方向为顺时针。将方框2处的涡旋截取出来,截取图像尺寸为24 km×24 km,如图8所示。ENVISAT-1 ASAR图像的具体雷达参数如表3所示。

    图  7  吕宋海峡获取的ENVISAT-1 ASAR图像,获取时间为2010.06.11, 01:51:48 UTC
    Figure  7.  ENVISAT-1 ASAR image of the Luson Strait obtained on June 11, 2010 at 01:51:48 UTC
    图  8  方框2处截取的涡旋SAR图像
    Figure  8.  Enlargement of the eddy in Frame 2
    表  3  ENVISAT-1 ASAR参数
    Table  3.  ASAR parameters of ENVISAT-1
    参数数值
    极化方式HH
    波段C
    入射角26.7°
    平台高度800 km
    平台速度7455 m/s
    下载: 导出CSV 
    | 显示表格

    从ECMWF获取2010.06.11, 03:00:00时刻的风场再分析资料,分辨率为0.125°×0.125°。根据数据显示,涡旋区域附近的风速为2.1 m/s,风向为45°。从GODAS获取相同位置的5日平均流场再分析资料,分辨率为(1/3)°×1°。根据数据显示,涡旋区域附近的流速为0.23 m/s。因此,设置参数α为0.000657,流场大小设置为24 km×24 km,空间分辨率为100 m,雷达参数设置为表3中ENVISAT-1 ASAR参数。

    图9(a)图9(b)分别是该涡旋的仿真SAR图像与真实SAR图像,仿真时设定的雷达参数、海面风场条件与真实SAR图像完全一致。对比图9(a)图9(b)两图发现,仿真SAR图像与真实SAR图像中的涡旋形状基本一致,涡旋臂的亮暗特征也较为吻合。从顺时针方向看,涡旋臂从外到内均呈现为由暗到亮的特征。该结果与Lyzenga等人[6]研究结果一致,初步验证了仿真方法的正确性。

    图  9  相同参数下仿真SAR图像与ENVISAT-1 ASAR图像对比图
    Figure  9.  Comparison of simulated SAR image and ENVISAT-1 ASAR image under the same parameters

    为了定量地描述仿真SAR图像与真实SAR图像的中涡旋的相似性,同样采用3.1节中的分析方法,得到涡旋拟合结果如图10所示,提取的涡旋信息如表4所示。

    图  10  仿真SAR图像与ENVISAT-1 ASAR图像涡旋信息提取
    Figure  10.  Eddy information extraction of simulated SAR image and ENVISAT-1 ASAR image
    表  4  涡旋信息提取结果
    Table  4.  Results of eddy information extraction
    SAR图像涡旋中心位置涡旋直径涡旋边缘尺寸
    仿真SAR图像(144,78)24.0 km49.4 km
    真实SAR图像(147,81)23.9 km49.7 km
    绝对/相对误差(3,3)/—0.1 km/0.0040.3 km/0.006
    下载: 导出CSV 
    | 显示表格

    对比仿真SAR图像与真实SAR图像的涡旋信息提取结果,可以发现两幅图像中涡旋的中心位置较为接近,方位向和距离向上均相差3个像素点,涡旋直径及边缘尺寸相对误差均不超过0.006,这进一步验证了仿真方法的正确性,说明本文提出的基于Burgers-Rott涡旋模型的涡旋SAR图像仿真方法能够实现反气旋式涡旋的SAR图像仿真。

    本文基于Burgers-Rott涡旋模型,提出了一种涡旋SAR图像仿真方法,并分别针对气旋式涡旋与反气旋式涡旋进行了仿真实验。通过将仿真SAR图像与真实SAR图像对比验证发现,本文提出的涡旋SAR图像仿真方法能够实现气旋式涡旋和反气旋式涡旋的SAR图像仿真,且仿真SAR图像与真实SAR图像能够较好地吻合。

    通过涡旋SAR图像仿真实验发现,无论气旋式涡旋还是反气旋式涡旋,其涡旋臂在SAR图像中都会呈现亮暗交替变化的特征。其中,气旋式涡旋臂呈现两个亮暗交替周期,即亮-暗-亮;反气旋式涡旋臂呈现一个亮暗交替周期,即暗-亮。这是由于这两个涡旋臂的曲率不同,气旋式涡旋臂曲率较大,亮暗交替周期较多,反气旋式涡旋臂曲率较小,亮暗交替周期较少。

    由于涡旋在SAR成像时会受到各种海洋环境因素的影响,通过真实SAR图像难以完全解译涡旋的特征。本文提出的SAR图像仿真方法能够弥补这种不足,可以清晰地获取涡旋的尺度、亮暗等特征,这为海洋涡旋特征的解译和提取提供了便利。

  • 图  1  基于相位调制编码的极化时变调控超表面及极化调控原理示意图

    Figure  1.  Schematic diagram of time-varying polarization-converting metasurface based on phase-modulated coding and its mechanism of polarization modulation

    图  2  不同相位量化位数及编码序列长度下,等效反射系数幅/相分布图

    Figure  2.  The amplitude-ratio and phase-difference distribution of the effective reflection coefficients with different coding bits and coding lengths

    图  3  各向异性相位调制表面单元的三维拓扑结构图及详细几何参数

    Figure  3.  3D topology expanded view of the anisotropy phase-modulated metasurface unit cell with the detailed geometricparameters

    图  4  相位调制表面各极化通道反射特性仿真分析

    Figure  4.  Simultated reflection characteristics of each polarization channel of the phase-modulated metasurface

    图  5  不同角度入射波照射下相位调制表面反射相位稳定性分析

    Figure  5.  Phase stability analysis of the metasurface for different incidence angles

    图  6  两种特定极化态对应的相位调制编码序列及偏置电压编码序列(入射波为45°线极化)

    Figure  6.  Coding sequence for the specific polarization states (with 45° linear polarization incidence)

    图  7  基于极化时变调控方法调控得到的3条极化轨道在Poincaré球上的分布

    Figure  7.  The distribution of three polarization trajectories on Poincaré sphere generated by the method proposed in this paper

    图  8  “大圆”轨道约束下的有源欺骗干扰与目标回波的极化域特征分布

    Figure  8.  Polarization domain characteristics of the target signal and the active detection jamming under the constraint of “great circle” trajectory

    图  9  “8”字轨道(过南北极)约束下的有源欺骗干扰与目标回波的极化域特征分布

    Figure  9.  Polarization domain characteristics of the target signal and the active detection jamming under the constraint of “8” likely trajectory through north and south poles

    图  10  “8”字轨道(未过南北极)约束下的有源欺骗干扰与目标回波的极化域特征分布

    Figure  10.  Polarization domain characteristics of the target signal and the active detection jamming under the constraint of “8” likely trajectory without north and south poles

    表  1  3 bit量化相位及偏置电压对应关系

    Table  1.   Corresponding relationship between 3 bit quantization phase and bias voltage

    水平极化 垂直极化
    VDC (V) Δφ (°) VDC (V) Δφ (°)
    0 0 0 0
    2.6 45 2.7 45
    3.2 89 3.5 90
    3.5 135 4.5 135
    3.9 180 4.6 179
    4.5 226 5.1 225
    6.0 270 7.0 270
    14.0 315 15.0 314
    下载: 导出CSV
  • [1] ROOME S J. Digital radio frequency memory[J]. Electronics & Communication Engineering Journal, 1990, 2(4): 147.
    [2] 曹旭源. 基于DRFM的雷达干扰技术研究[D]. [硕士论文], 西安电子科技大学, 2013: 5–19.

    CAO Xuyuan. Research on jamming against radar based on digital radio frequency memory[D]. [Master dissertation], Xidian University, 2013: 5–19.
    [3] 高佳旭. 雷达信号有源干扰技术研究及实现[D]. [硕士论文], 哈尔滨工程大学, 2021: 13–18.

    GAO Jiaxu. Research and realization of radar signal active jamming technology[D]. [Master dissertation], Harbin Engineering University, 2021: 13–18.
    [4] 肖顺平, 徐振海, 代大海, 等. 雷达极化技术[M]. 北京: 清华大学出版社, 2022: 1–51.

    XIAO Shunping, XU Zhenhai, DAI Dahai, et al. Radar Polarization Techniques[M]. Beijing: Tsinghua University Press, 2022: 1–51.
    [5] 王雪松. 雷达极化技术研究现状与展望[J]. 雷达学报, 2016, 5(2): 119–131. doi: 10.12000/JR16039.

    WANG Xuesong. Status and prospects of radar polarimetry techniques[J]. Journal of Radars, 2016, 5(2): 119–131. doi: 10.12000/JR16039.
    [6] 施龙飞, 马佳智, 庞晨, 等. 极化雷达信号处理与抗干扰技术[M]. 北京: 国防工业出版社, 2019: 1–13.

    SHI Longfei, MA Jiazhi, PANG Chen, et al. Signal Processing and Anti-interference Techniques for Polarimetric Radar[M]. Beijing: National Defense Industry Press, 2019: 1–13.
    [7] 马佳智, 施龙飞, 徐振海, 等. 单脉冲雷达多点源参数估计与抗干扰技术进展[J]. 雷达学报, 2019, 8(1): 125–139. doi: 10.12000/JR18093.

    MA Jiazhi, SHI Longfei, XU Zhenhai, et al. Overview of multi-source parameter estimation and jamming mitigation for monopulse radars[J]. Journal of Radars, 2019, 8(1): 125–139. doi: 10.12000/JR18093.
    [8] 施龙飞, 任博, 马佳智, 等. 雷达极化抗干扰技术进展[J]. 现代雷达, 2016, 38(4): 1–7, 29. doi: 10.16592/j.cnki.1004-7859.2016.04.001.

    SHI Longfei, REN Bo, MA Jiazhi, et al. Recent developments of radar anti-interference techniques with polarimetry[J]. Modern Radar, 2016, 38(4): 1–7, 29. doi: 10.16592/j.cnki.1004-7859.2016.04.001.
    [9] 李永祯, 王雪松, 肖顺平, 等. 基于IPPV的真假目标极化鉴别算法[J]. 现代雷达, 2004, 26(9): 38–42. doi: 10.3969/j.issn.1004-7859.2004.09.011.

    LI Yongzhen, WANG Xuesong, XIAO Shunping, et al. A new polarization discrimination algorithm for active decoy and radar target based on IPPV[J]. Modern Radar, 2004, 26(9): 38–42. doi: 10.3969/j.issn.1004-7859.2004.09.011.
    [10] 施龙飞, 王雪松, 肖顺平. 转发式假目标干扰的极化鉴别[J]. 中国科学 F辑: 信息科学, 2009, 39(4): 468–475.

    SHI Longfei, WANG Xuesong, and XIAO Shunping. Polarization discrimination between repeater false-target and radar target[J]. Science in China Series F: Information Sciences, 2009, 39(4): 468–475.
    [11] ZONG Zhiwei, SHI Longfei, LI Yongzhen, et al. Detection-discrimination method for multiple repeater false targets based on radar polarization echoes[J]. Radioengineering, 2014, 23(1): 104–110.
    [12] 施龙飞, 帅鹏, 王雪松, 等. 极化调制假目标干扰的鉴别[J]. 信号处理, 2008, 24(6): 894–899. doi: 10.3969/j.issn.1003-0530.2008.06.002.

    SHI Longfei, SHUAI Peng, WANG Xuesong, et al. Polarization discrimination between modulation polarization decoy and radar target[J]. Signal Processing, 2008, 24(6): 894–899. doi: 10.3969/j.issn.1003-0530.2008.06.002.
    [13] 王雪松. 宽带极化信息处理的研究[D]. [博士论文], 国防科学技术大学, 1999: 21–26, 127–179.

    WANG Xuesong. Study on wide-band polarization information processing[D]. [Ph.D. dissertation], National University of Defense Technology, 1999: 21–26, 127–179.
    [14] 蒋仁培, 苏丽萍. 雷达极化问题和铁氧体变极化技术[J]. 现代雷达, 2001, 23(1): 65–69, 72. doi: 10.3969/j.issn.1004-7859.2001.01.014.

    JIANG Renpei and SU Liping. Problem of radar polarization and technique of ferrite variable polarization[J]. Modern Radar, 2001, 23(1): 65–69, 72. doi: 10.3969/j.issn.1004-7859.2001.01.014.
    [15] 蒋仁培, 苏丽萍, 魏克珠. 广义铁氧体变极化理论[J]. 微波学报, 2000, 16(4): 336–342. doi: 10.3969/j.issn.1005-6122.2000.04.002.

    JIANG Renpei, SU Liping, and WEI Kezhu. Generalized theory of ferrite variable polarization[J]. Journal of Microwaves, 2000, 16(4): 336–342. doi: 10.3969/j.issn.1005-6122.2000.04.002.
    [16] KARAMIRAD M, GHOBADI C, and NOURINIA J. Metasurfaces for wideband and efficient polarization rotation[J]. IEEE Transactions on Antennas and Propagation, 2021, 69(3): 1799–1804. doi: 10.1109/TAP.2020.3012828.
    [17] WANG Yidan, SHI Hongyu, CHEN Juan, et al. Digital polarization programmable metasurface for continuous polarization angle rotation and radar applications[J]. Frontiers in Materials, 2022, 9: 931868. doi: 10.3389/fmats.2022.931868.
    [18] TIAN Jianghao, CAO Xiangyu, GAO Jun, et al. A reconfigurable ultra-wideband polarization converter based on metasurface incorporated with PIN diodes[J]. Journal of Applied Physics, 2019, 125(13): 135105. doi: 10.1063/1.5067383.
    [19] YANG Zhengyi, KOU Na, YU Shixing, et al. Reconfigurable multifunction polarization converter integrated with PIN diode[J]. IEEE Microwave and Wireless Components Letters, 2021, 31(6): 557–560. doi: 10.1109/LMWC.2021.3064039.
    [20] YANG Heng, WANG Shicong, LI Peng, et al. A broadband multifunctional reconfigurable polarization conversion metasurface[J]. IEEE Transactions on Antennas and Propagation, 2023, 71(7): 5759–5767. doi: 10.1109/TAP.2023.3266498.
    [21] LI You, WANG Yi, and CAO Qunsheng. Design of a multifunctional reconfigurable metasurface for polarization and propagation manipulation[J]. IEEE Access, 2019, 7: 129183–129191. doi: 10.1109/ACCESS.2019.2939200.
    [22] CERVENY M, FORD K L, and TENNANT A. Reflective switchable polarization rotator based on metasurface with PIN diodes[J]. IEEE Transactions on Antennas and Propagation, 2021, 69(3): 1483–1492. doi: 10.1109/TAP.2020.3026883.
    [23] LIN Baoqin, GUO Jianxin, LV Lintao, et al. Ultra-wideband and high-efficiency reflective polarization converter for both linear and circular polarized waves[J]. Applied Physics A, 2019, 125(2): 76. doi: 10.1007/s00339-018-2368-9.
    [24] HAO Jiaming, YUAN Yu, RAN Lixin, et al. Manipulating electromagnetic wave polarizations by anisotropic metamaterials[J]. Physical Review Letters, 2007, 99(6): 063908. doi: 10.1103/PhysRevLett.99.063908.
    [25] WANG Shuai, DENG Zilan, WANG Yujie, et al. Arbitrary polarization conversion dichroism metasurfaces for all-in-one full Poincaré sphere polarizers[J]. Light: Science & Applications, 2021, 10(1): 24. doi: 10.1038/s41377-021-00468-y.
    [26] KE Junchen, DAI Junyan, CHEN Mingzheng, et al. Linear and nonlinear polarization syntheses and their programmable controls based on anisotropic time-domain digital coding metasurface[J]. Small Structures, 2021, 2(1): 2000060. doi: 10.1002/sstr.202000060.
    [27] HU Qi, CHEN Ke, ZHANG Na, et al. Arbitrary and dynamic Poincaré sphere polarization converter with a time-varying metasurface[J]. Advanced Optical Materials, 2022, 10(4): 2101915. doi: 10.1002/adom.202101915.
    [28] ZHANG Xinge, YU Qian, JIANG Weixiang, et al. Polarization-controlled dual-programmable metasurfaces[J]. Advanced Science, 2020, 7(11): 1903382. doi: 10.1002/advs.201903382.
    [29] 胡琪, 陈克, 郑依琳, 等. 时变极化编码表面及其在无线通信中的应用[J]. 雷达学报, 2021, 10(2): 304–312. doi: 10.12000/JR21042.

    HU Qi, CHEN Ke, ZHENG Yilin, et al. Time-varying polarization-converting programmable metasurface and its application in wireless communication system[J]. Journal of Radars, 2021, 10(2): 304–312. doi: 10.12000/JR21042.
    [30] DAI Junyan, ZHAO Jie, CHENG Qiang, et al. Independent control of harmonic amplitudes and phases via a time-domain digital coding metasurface[J]. Light: Science & Applications, 2018, 7: 90.
    [31] 戴俊彦. 时域超表面理论研究与应用[D]. [博士论文], 东南大学, 2019: 41.

    DAI Junyan. Research and application of time-domain metasurface[D]. [Ph.D. dissertation], Southeast University, 2019: 41.
    [32] 程强, 崔铁军. 电磁超材料[M]. 南京: 东南大学出版社, 2022: 259–262.

    CHENG Qiang and CUI Tiejun. Metamaterials[M]. Nanjing: Southeast University Press, 2022: 259–262.
    [33] SUN Guang, WANG Junjie, XING Shiqi, et al. A flexible conformal multifunctional time-modulated metasurface for radar characteristics manipulation[J]. IEEE Transactions on Microwave Theory and Techniques, 2023: 1–15. doi: 10.1109/TMTT.2023.3337646.
    [34] MACOM. MA46H120 varactor datasheet[EB/OL]. https://www.macom.com/products/product-detail/MA46H120.
  • 期刊类型引用(3)

    1. 阮航,崔家豪,毛秀华,任建迎,罗镔延,曹航,李海峰. SAR目标识别对抗攻击综述:从数字域迈向物理域. 雷达学报. 2024(06): 1298-1326 . 本站查看
    2. 翟多才,黄贤俊,李元鑫,徐延林,姚理想,林铭团,刘培国. 一种角度与极化稳定的可调控电磁散射结构的设计与验证. 电波科学学报. 2023(02): 205-210 . 百度学术
    3. 朱泳庚,艾夏,王伟东,李猛猛,丁大志. 新型多维电磁特性可调控雷达角反射器设计方法研究. 电波科学学报. 2023(06): 1015-1021+1056 . 百度学术

    其他类型引用(1)

  • 加载中
图(10) / 表(1)
计量
  • 文章访问数: 490
  • HTML全文浏览量: 123
  • PDF下载量: 231
  • 被引次数: 4
出版历程
  • 收稿日期:  2024-02-23
  • 修回日期:  2024-04-25
  • 网络出版日期:  2024-05-30
  • 刊出日期:  2024-08-28

目录

/

返回文章
返回