CN110988868A - 一种等效分数阶模式涡旋电磁波产生及成像方法 - Google Patents

一种等效分数阶模式涡旋电磁波产生及成像方法 Download PDF

Info

Publication number
CN110988868A
CN110988868A CN201911134517.9A CN201911134517A CN110988868A CN 110988868 A CN110988868 A CN 110988868A CN 201911134517 A CN201911134517 A CN 201911134517A CN 110988868 A CN110988868 A CN 110988868A
Authority
CN
China
Prior art keywords
mode
imaging
vortex electromagnetic
fractional order
order mode
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201911134517.9A
Other languages
English (en)
Other versions
CN110988868B (zh
Inventor
陈如山
何姿
丁大志
国少卿
印必还
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201911134517.9A priority Critical patent/CN110988868B/zh
Publication of CN110988868A publication Critical patent/CN110988868A/zh
Application granted granted Critical
Publication of CN110988868B publication Critical patent/CN110988868B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging

Abstract

本发明公开了一种等效分数阶模式涡旋电磁波产生及成像方法。该方法首次采用同心圆环阵列得到等效分数阶模式并用于雷达成像。其实现方案为:首先利用同心均匀圆环阵列产生整数阶涡旋波模式,通过设定收、发模式为相邻整数从而等效得到分数阶模式,针对单圈圆环阵列产生的不同模式涡旋波最大辐射方向不同的问题,本发明采用多圈同心圆环阵列,通过优化各圈馈电幅度相位,使得各模式最大辐射方向相同,得到更好的成像效果。本发明可以解决现有涡旋电磁波成像在多发多收模式下存在的方位角混叠问题。

Description

一种等效分数阶模式涡旋电磁波产生及成像方法
技术领域
本发明属于雷达成像技术领域,具体是一种产生等效分数阶涡旋波模式的方法。
背景技术
涡旋电磁波由于具有独特的螺旋型相位波前分布,近年来引起了研究者的广泛关注。当涡旋电磁波携带的轨道角动量模式为α时,在与传播方向垂直的平面内,其相位绕传播轴旋转一周将变化2πα。由于涡旋电磁波的不同模式之间相互正交,可以用于通信以提高频谱效率。另外,由于不同模式的涡旋电磁波的相位波前分布不同,有望可以用于雷达探测领域以获取更丰富的目标信息。
目前文献公开报道的产生涡旋电磁波的方式有很多种,根据产生原理大致可分为三类:透射型、反射型以及相位调制型。其中利用相位调制原理产生涡旋波如均匀圆环阵列,已被用于雷达成像中。与前面两种产生方式相比,其优势在于,可以用同一个阵列产生多个涡旋波模式,且容易实现。基于涡旋电磁波的雷达成像与传统成像方法相比,雷达与目标之间不需要相对运动,而是利用涡旋波的模式数与目标方位角的对偶关系便可实现成像。利用涡旋电磁波成像有两种方式:多发多收及多发单收。前者的方位角分辨率优于后者,但存在方位角混叠现象。为实现方位角的高分辨并消除方位角混叠的问题,可在多发多收成像方式中引入分数阶模式。但直接产生的分数阶模式并不稳定。
发明内容
本发明的目的在于提出一种等效分数阶模式涡旋电磁波产生及成像方法。
实现本发明目的的技术解决方案为:一种产生等效分数阶涡旋波模式的方法,步骤如下:
第一步,采用多圈同心均匀圆环阵列产生若干辐射方向相同且旁瓣较低的整数阶模式的涡旋电磁波;
第二步,发送某一整数阶模式的涡旋电磁波,并采用相邻的另一整数阶模式接收,得到等效分数阶模式的涡旋电磁波;
第三步,利用整数及分数阶模式的涡旋电磁波并采用多发多收的方式进行成像。
本发明与现有的多发多收成像相比,其显著优点是:能够消除方位混叠现象;同时,与多发单收成像相比,其方位分辨率更高。
附图说明
图1是本发明所用的同心圆环阵列示意图。
图2是本发明所用的同心圆环阵列产生的不同模式的方向图。
图3是本发明方法所仿真的理想散射点构成的两个“飞机”模型。
图4是采用整数阶模式的涡旋波多发多收成像与采用本发明方法进行成像的结果对比图。
具体实施方式
下面结合说明书附图对本发明作进一步描述。
本发明一种等效分数阶模式涡旋电磁波产生及成像方法,步骤如下:
第一步,采用多圈同心均匀圆环阵列产生若干辐射方向相同且旁瓣较低的整数阶模式的涡旋电磁波。具体方法如下:
同心圆环阵列结构如图1所示。记构成同心圆环阵列中的圆环数为N,半径分别为a1~aN,各圈圆环上的单元均匀分布,记第j圈圆环上的单元数位Nj。当要产生模式为α的涡旋电磁波时,第j圈圆环上天线第t个单元激励相位应为为
Figure BDA0002279224600000021
t=1,2,…, Nj,激励信号幅度记为Ij,当圆环上单元数Nj足够多时,同心圆环阵列的辐射方向图表达式为
Figure BDA0002279224600000022
其中,f0(θ,φ)为阵列单元方向图,
Figure BDA0002279224600000023
为第j圈圆环阵归一化系数,Jα(kajsinθ)为α阶第一类贝塞尔函数,aj为第j圈圆环半径,k为波数。本发明中所用的单元方向图为
Figure BDA0002279224600000024
若要考虑用实际的天线单元(如微带、喇叭天线等)进行组阵,只需将上式换为所用到的实际天线单元方向图即可。单圈圆环阵列所产生的涡旋电磁波的最大辐射方向与圆环半径及模式的关系如下:
θmax≈arcsin[(1.0509α+1.1562)/(ka)] (3)
可以看出,对相同半径的单圈圆环阵列,高模式的涡旋电磁波的主瓣指向将变大。采用多圈同心圆环阵列可以解决这一问题。为使得同心圆环阵列所产生的所有模式的涡旋电磁波辐射方向相同,并且旁瓣电平较低。需要对各圆环馈电幅度进行优化。定义目标函数为
fit=w1objmax|+w2|SLLobj-SLLmax| (4)
其中θobj和θmax分别表示设计的主瓣目标指向角度值和优化过程中计算得到的方向图主瓣指向角度,SLLobj和SLLmax分别为设计的归一化最大副瓣电平值与优化过程中计算得到的归一化最大副瓣电平值。w1和w2为权重系数,可以根据需要产生的不同模式进行调整,以得到最优结果。
本发明采用自适应协方差矩阵进化策略(CMA-ES)对各圆环激励幅度进行优化。CMA-ES是一种可以从复杂参数的相关性中进行学习的自适应算法,这种自适应算法通过学习获得较优的进化参数,可以显着影响优化性能,与常用的基因算法相比,其优化速度更快并且鲁棒性更好,适合解决具有较大维度变量的优化问题。
本发明所设计的同心圆环阵列中心频率为fc=6GHz,共10圈,所生成的涡旋电磁波最大整数阶模式为10,由于正负模式的馈电幅值相同,馈电相位相反,因此可以得到的整数阶模式为-10~10。设定各模式的主波束指向为θobj=11°。经过优化后,各圈圆环激励幅度如表1所示。
Figure BDA0002279224600000031
表1本发明所优化得到的同心圆环阵列产生不同模式涡旋波时各圈馈电幅值
图为优化后的不同模式的相位分布特性,图2为优化后的同心圆环阵列所产生的不同模式数涡旋电磁波的辐射方向图。可以看出,不同模式的涡旋电磁波,其最大辐射方向基本一致,旁瓣电平基本小于-20dB。
第二步,发送某一整数阶模式的涡旋电磁波并采用相邻的另一整数阶模式接收,得到等效分数阶模式的涡旋电磁波。同心均匀圆环阵列共产生-10~10的整数阶模式。具体实施方法如下:
对空间中的M个点目标,当采用第一步中的同心圆环阵列发射某一模式为α的涡旋波时,在接收端单元的接收信号可表示为:
Figure BDA0002279224600000041
其中,M为空间中的散射点的个数,
Figure BDA0002279224600000042
为第m个散射点的散射系数,(rmmm)为第 m个散射点在空间中的位置。由于采用多发多收的成像方式,当以α模式发射,并以β模式接收,回波信号可表示为:
Figure BDA0002279224600000043
从上述回波信号表达式可以看出,其中含有相位因子
Figure BDA0002279224600000044
即等效于采用了(α+β)/2模式进行收发。在本发明中,取α=β+1,可以等效得到模式为α+0.5的涡旋电磁波。
第三步,利用一定范围内的整数及分数阶模式采用多发多收模式进行成像。本发明采用多发多收方式,利用快速傅立叶变换进行二维成像。本发明共产生的整数阶模式为-10~10。当发射、接收模式分别为-10和-9时,可以等效获得模式为-9.5的涡旋电磁波。将所有(-L,-L+0.5,-L+1,……L-1,L-0.5,L)用于成像。为验证本发明所提出的方法的有效性,设置了方位角相距180度的散射点组成的两组“飞机”模型,如图3所示。图4为本发明实现二维成像结果,其中中心频率f=6e9,带宽500Mhz。当采用传统的整数阶模式的多发多收进行成像时,成像结果如图4(a)所示,由于方位混叠现象,只能观察到一个“飞机”模型。当采用本发明的方法,成像结果如图4(b)所示,能正确得到两个“飞机”的位置。从而证明了本发明的有效性。

Claims (4)

1.一种等效分数阶模式涡旋电磁波产生及成像方法,其特征在于,步骤如下:
第一步,采用多圈同心均匀圆环阵列产生若干辐射方向相同且旁瓣较低的整数阶模式的涡旋电磁波;
第二步,发送某一整数阶模式的涡旋电磁波,并采用相邻的另一整数阶模式接收,得到等效分数阶模式的涡旋电磁波;
第三步,利用整数及分数阶模式的涡旋电磁波并采用多发多收的方式进行成像。
2.根据权利要求1所述的等效分数阶模式涡旋电磁波产生及成像方法,其特征在于,所述第一步中:通过采用自适应协方差矩阵进化策略CMA-ES优化同心均匀圆环阵列各圈的单元馈电幅值,使产生的所有整数阶模式的涡旋电磁波最大辐射方向相同,并且最大旁瓣电平低于主瓣电平20dB。
3.根据权利要求1所述的等效分数阶模式涡旋电磁波产生及成像方法,其特征在于,所述第二步中:当发送整数阶模式为α时,设置接收模式为α+1或α-1,得到等效分数阶模式α+0.5或α-0.5。
4.根据权利要求1所述的等效分数阶模式涡旋电磁波产生及成像方法,其特征在于,所述第三步中:成像模式采用多发多收方式,发射信号的中心频率为f c ,带宽为B,当产生的整数阶模式为(-L~L)时,结合分数阶模式,利用得到的(-L,-L+0.5,-L+1,……L-1,L-0.5,L)模式,通过二维快速傅立叶变换得到二维图像。
CN201911134517.9A 2019-11-19 2019-11-19 一种等效分数阶模式涡旋电磁波产生及成像方法 Active CN110988868B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911134517.9A CN110988868B (zh) 2019-11-19 2019-11-19 一种等效分数阶模式涡旋电磁波产生及成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911134517.9A CN110988868B (zh) 2019-11-19 2019-11-19 一种等效分数阶模式涡旋电磁波产生及成像方法

Publications (2)

Publication Number Publication Date
CN110988868A true CN110988868A (zh) 2020-04-10
CN110988868B CN110988868B (zh) 2023-09-26

Family

ID=70084907

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911134517.9A Active CN110988868B (zh) 2019-11-19 2019-11-19 一种等效分数阶模式涡旋电磁波产生及成像方法

Country Status (1)

Country Link
CN (1) CN110988868B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112068129A (zh) * 2020-08-21 2020-12-11 西安空间无线电技术研究所 一种同心圆环阵列电磁涡旋雷达成像方法
CN112083430A (zh) * 2020-09-10 2020-12-15 天津水聿方舟海洋工程技术有限公司 一种适用于轨道角动量三维成像声呐的旁瓣抑制方法
CN112505689A (zh) * 2020-08-21 2021-03-16 西安空间无线电技术研究所 一种电磁涡旋雷达同心圆环阵列设计方法
CN112764041A (zh) * 2021-01-19 2021-05-07 鹏城实验室 成像系统和方法
CN117706490A (zh) * 2024-02-06 2024-03-15 南京理工大学 一种基于单站雷达的金属目标间的耦合散射中心建模方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103954367A (zh) * 2014-04-18 2014-07-30 河南科技大学 一种测量分数阶光学涡旋拓扑荷值的装置及其测量方法
CN105445943A (zh) * 2015-12-24 2016-03-30 河南科技大学 一种分数阶完美涡旋光束的产生装置及产生方法
CN106059675A (zh) * 2016-05-24 2016-10-26 陕西理工学院 圆阵列旋转分形嵌套高阶涡旋电磁波生成分离方法与装置
CN106058490A (zh) * 2016-06-01 2016-10-26 电子科技大学 一种产生涡旋电磁波的方法
CN107238933A (zh) * 2017-07-07 2017-10-10 苏州大学 产生局部分数阶贝塞尔涡旋光束的方法和系统
CN109616777A (zh) * 2018-12-28 2019-04-12 西安电子科技大学 基于时间调制阵列的任意模态涡旋电磁波产生方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103954367A (zh) * 2014-04-18 2014-07-30 河南科技大学 一种测量分数阶光学涡旋拓扑荷值的装置及其测量方法
CN105445943A (zh) * 2015-12-24 2016-03-30 河南科技大学 一种分数阶完美涡旋光束的产生装置及产生方法
CN106059675A (zh) * 2016-05-24 2016-10-26 陕西理工学院 圆阵列旋转分形嵌套高阶涡旋电磁波生成分离方法与装置
CN106058490A (zh) * 2016-06-01 2016-10-26 电子科技大学 一种产生涡旋电磁波的方法
CN107238933A (zh) * 2017-07-07 2017-10-10 苏州大学 产生局部分数阶贝塞尔涡旋光束的方法和系统
CN109616777A (zh) * 2018-12-28 2019-04-12 西安电子科技大学 基于时间调制阵列的任意模态涡旋电磁波产生方法

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112068129A (zh) * 2020-08-21 2020-12-11 西安空间无线电技术研究所 一种同心圆环阵列电磁涡旋雷达成像方法
CN112505689A (zh) * 2020-08-21 2021-03-16 西安空间无线电技术研究所 一种电磁涡旋雷达同心圆环阵列设计方法
CN112083430A (zh) * 2020-09-10 2020-12-15 天津水聿方舟海洋工程技术有限公司 一种适用于轨道角动量三维成像声呐的旁瓣抑制方法
CN112083430B (zh) * 2020-09-10 2023-11-07 天津水聿方舟海洋工程技术有限公司 一种适用于轨道角动量三维成像声呐的旁瓣抑制方法
CN112764041A (zh) * 2021-01-19 2021-05-07 鹏城实验室 成像系统和方法
CN112764041B (zh) * 2021-01-19 2023-05-26 鹏城实验室 成像系统和方法
CN117706490A (zh) * 2024-02-06 2024-03-15 南京理工大学 一种基于单站雷达的金属目标间的耦合散射中心建模方法

Also Published As

Publication number Publication date
CN110988868B (zh) 2023-09-26

Similar Documents

Publication Publication Date Title
CN110988868B (zh) 一种等效分数阶模式涡旋电磁波产生及成像方法
CN106886020B (zh) 一种单天线接收条件下的电磁涡旋成像方法
CN103592642B (zh) Mimo雷达波形的设计方法
CN106486784A (zh) 反射式天线阵列及波束扫描方法
CN104615854B (zh) 一种基于稀疏约束的波束展宽和旁瓣抑制方法
CN105929384B (zh) 基于十字阵发射圆环阵接收的联合低旁瓣多波束形成方法
EP0807992B1 (en) Logarithmic spiral array
CN107645068A (zh) 一种旋转圆极化阵元产生oam波束的圆形天线阵设计方法
CN112083430B (zh) 一种适用于轨道角动量三维成像声呐的旁瓣抑制方法
CN114357767B (zh) 一种宽带宽角波束扫描的和差波束稀疏阵列综合方法
Guo et al. CUCA based equivalent fractional order OAM mode for electromagnetic vortex imaging
CN111352080A (zh) Papr和相似性约束下低截获频控阵mimo雷达系统的设计方法
Khalid et al. Synthesis of linear antenna array using genetic algorithm to reduce peak sidelobe level
CN110850383B (zh) 一种基于共形阵的mimo雷达信号处理方法
CN104808178B (zh) 一种机载雷达发射方向图设计方法
CN116224240A (zh) 基于最大输出信噪比的共形极化阵发射方向图赋形算法
CN113540791B (zh) 一种孔径级收发同时阵列优化方法
CN114386271A (zh) 一种考虑互耦效应的随机阵列天线方向图综合方法
CN112505689A (zh) 一种电磁涡旋雷达同心圆环阵列设计方法
Gu et al. An Efficient Rotationally Symmetric Approach for the Design of Sparse Conformal Arrays in Wide Angle Scanning.
Sun et al. Shared aperture multibeam forming of time-modulated linear array
Mahato et al. Null synthesis of scanned linear array antenna with minimum side lobe level and fixed dynamic range ratio using iterative fast Fourier transform
Ren et al. A Novel Optimal Subarray Partition Method for Large Uniform Linear Array based on Hierarchical Clustering Algorithm
Sun et al. Circular truncated cone conformal array beamforming based on directional array element
Xuejun et al. Design of sparse concentric multi-circular array using micro-genetic optimum algorithm

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant