CN110988868B - 一种等效分数阶模式涡旋电磁波产生及成像方法 - Google Patents
一种等效分数阶模式涡旋电磁波产生及成像方法 Download PDFInfo
- Publication number
- CN110988868B CN110988868B CN201911134517.9A CN201911134517A CN110988868B CN 110988868 B CN110988868 B CN 110988868B CN 201911134517 A CN201911134517 A CN 201911134517A CN 110988868 B CN110988868 B CN 110988868B
- Authority
- CN
- China
- Prior art keywords
- mode
- vortex electromagnetic
- imaging
- order mode
- alpha
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种等效分数阶模式涡旋电磁波产生及成像方法。该方法首次采用同心圆环阵列得到等效分数阶模式并用于雷达成像。其实现方案为:首先利用同心均匀圆环阵列产生整数阶涡旋波模式,通过设定收、发模式为相邻整数从而等效得到分数阶模式,针对单圈圆环阵列产生的不同模式涡旋波最大辐射方向不同的问题,本发明采用多圈同心圆环阵列,通过优化各圈馈电幅度相位,使得各模式最大辐射方向相同,得到更好的成像效果。本发明可以解决现有涡旋电磁波成像在多发多收模式下存在的方位角混叠问题。
Description
技术领域
本发明属于雷达成像技术领域,具体是一种产生等效分数阶涡旋波模式的方法。
背景技术
涡旋电磁波由于具有独特的螺旋型相位波前分布,近年来引起了研究者的广泛关注。当涡旋电磁波携带的轨道角动量模式为α时,在与传播方向垂直的平面内,其相位绕传播轴旋转一周将变化2πα。由于涡旋电磁波的不同模式之间相互正交,可以用于通信以提高频谱效率。另外,由于不同模式的涡旋电磁波的相位波前分布不同,有望可以用于雷达探测领域以获取更丰富的目标信息。
目前文献公开报道的产生涡旋电磁波的方式有很多种,根据产生原理大致可分为三类:透射型、反射型以及相位调制型。其中利用相位调制原理产生涡旋波如均匀圆环阵列,已被用于雷达成像中。与前面两种产生方式相比,其优势在于,可以用同一个阵列产生多个涡旋波模式,且容易实现。基于涡旋电磁波的雷达成像与传统成像方法相比,雷达与目标之间不需要相对运动,而是利用涡旋波的模式数与目标方位角的对偶关系便可实现成像。利用涡旋电磁波成像有两种方式:多发多收及多发单收。前者的方位角分辨率优于后者,但存在方位角混叠现象。为实现方位角的高分辨并消除方位角混叠的问题,可在多发多收成像方式中引入分数阶模式。但直接产生的分数阶模式并不稳定。
发明内容
本发明的目的在于提出一种等效分数阶模式涡旋电磁波产生及成像方法。
实现本发明目的的技术解决方案为:一种产生等效分数阶涡旋波模式的方法,步骤如下:
第一步,采用多圈同心均匀圆环阵列产生若干辐射方向相同且旁瓣较低的整数阶模式的涡旋电磁波;
第二步,发送某一整数阶模式的涡旋电磁波,并采用相邻的另一整数阶模式接收,得到等效分数阶模式的涡旋电磁波;
第三步,利用整数及分数阶模式的涡旋电磁波并采用多发多收的方式进行成像。
本发明与现有的多发多收成像相比,其显著优点是:能够消除方位混叠现象;同时,与多发单收成像相比,其方位分辨率更高。
附图说明
图1是本发明所用的同心圆环阵列示意图。
图2是本发明所用的同心圆环阵列产生的不同模式的方向图。
图3是本发明方法所仿真的理想散射点构成的两个“飞机”模型。
图4是采用整数阶模式的涡旋波多发多收成像与采用本发明方法进行成像的结果对比图。
具体实施方式
下面结合说明书附图对本发明作进一步描述。
本发明一种等效分数阶模式涡旋电磁波产生及成像方法,步骤如下:
第一步,采用多圈同心均匀圆环阵列产生若干辐射方向相同且旁瓣较低的整数阶模式的涡旋电磁波。具体方法如下:
同心圆环阵列结构如图1所示。记构成同心圆环阵列中的圆环数为N,半径分别为a1~aN,各圈圆环上的单元均匀分布,记第j圈圆环上的单元数位Nj。当要产生模式为α的涡旋电磁波时,第j圈圆环上天线第t个单元激励相位应为为t=1,2,…,Nj,激励信号幅度记为Ij,当圆环上单元数Nj足够多时,同心圆环阵列的辐射方向图表达式为
其中,f0(θ,φ)为阵列单元方向图,为第j圈圆环阵归一化系数,Jα(kajsinθ)为α阶第一类贝塞尔函数,aj为第j圈圆环半径,k为波数。本发明中所用的单元方向图为
若要考虑用实际的天线单元(如微带、喇叭天线等)进行组阵,只需将上式换为所用到的实际天线单元方向图即可。单圈圆环阵列所产生的涡旋电磁波的最大辐射方向与圆环半径及模式的关系如下:
θmax≈arcsin[(1.0509α+1.1562)/(ka)] (3)
可以看出,对相同半径的单圈圆环阵列,高模式的涡旋电磁波的主瓣指向将变大。采用多圈同心圆环阵列可以解决这一问题。为使得同心圆环阵列所产生的所有模式的涡旋电磁波辐射方向相同,并且旁瓣电平较低。需要对各圆环馈电幅度进行优化。定义目标函数为
fit=w1|θobj-θmax|+w2|SLLobj-SLLmax| (4)
其中θobj和θmax分别表示设计的主瓣目标指向角度值和优化过程中计算得到的方向图主瓣指向角度,SLLobj和SLLmax分别为设计的归一化最大副瓣电平值与优化过程中计算得到的归一化最大副瓣电平值。w1和w2为权重系数,可以根据需要产生的不同模式进行调整,以得到最优结果。
本发明采用自适应协方差矩阵进化策略(CMA-ES)对各圆环激励幅度进行优化。CMA-ES是一种可以从复杂参数的相关性中进行学习的自适应算法,这种自适应算法通过学习获得较优的进化参数,可以显着影响优化性能,与常用的基因算法相比,其优化速度更快并且鲁棒性更好,适合解决具有较大维度变量的优化问题。
本发明所设计的同心圆环阵列中心频率为fc=6GHz,共10圈,所生成的涡旋电磁波最大整数阶模式为10,由于正负模式的馈电幅值相同,馈电相位相反,因此可以得到的整数阶模式为-10~10。设定各模式的主波束指向为θobj=11°。经过优化后,各圈圆环激励幅度如表1所示。
表1本发明所优化得到的同心圆环阵列产生不同模式涡旋波时各圈馈电幅值
图为优化后的不同模式的相位分布特性,图2为优化后的同心圆环阵列所产生的不同模式数涡旋电磁波的辐射方向图。可以看出,不同模式的涡旋电磁波,其最大辐射方向基本一致,旁瓣电平基本小于-20dB。
第二步,发送某一整数阶模式的涡旋电磁波并采用相邻的另一整数阶模式接收,得到等效分数阶模式的涡旋电磁波。同心均匀圆环阵列共产生-10~10的整数阶模式。具体实施方法如下:
对空间中的M个点目标,当采用第一步中的同心圆环阵列发射某一模式为α的涡旋波时,在接收端单元的接收信号可表示为:
其中,M为空间中的散射点的个数,为第m个散射点的散射系数,(rm,θm,φm)为第m个散射点在空间中的位置。由于采用多发多收的成像方式,当以α模式发射,并以β模式接收,回波信号可表示为:
从上述回波信号表达式可以看出,其中含有相位因子即等效于采用了(α+β)/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 (3)
1.一种等效分数阶模式涡旋电磁波产生及成像方法,其特征在于,步骤如下:
第一步,采用多圈同心均匀圆环阵列产生若干辐射方向相同且旁瓣低的整数阶模式的涡旋电磁波;通过采用自适应协方差矩阵进化策略CMA-ES优化同心均匀圆环阵列各圈的单元馈电幅值,使产生的所有整数阶模式的涡旋电磁波最大辐射方向相同,并且最大旁瓣电平低于主瓣电平20dB;
第二步,发送某一整数阶模式的涡旋电磁波,并采用相邻的另一整数阶模式接收,得到等效分数阶模式的涡旋电磁波;
具体实施方法如下:
对空间中的M个点目标,当采用第一步中的同心圆环阵列发射某一模式为α的涡旋波时,在接收端单元的接收信号可表示为:
其中,构成同心圆环阵列中的圆环数为N,M为空间中的散射点的个数,为第m个散射点的散射系数,(rm,θm,φm)为第m个散射点在空间中的位置,激励信号幅度记为Ij,Jα(kajsinθ)为α阶第一类贝塞尔函数,aj为第j圈圆环半径,k为波数;由于采用多发多收的成像方式,当以α模式发射,并以β模式接收,回波信号可表示为:
从上述回波信号表达式看出,其中含有相位因子即等效于采用了(α+β)/2模式进行收发;
第三步,利用整数及分数阶模式的涡旋电磁波并采用多发多收的方式进行成像。
2.根据权利要求1所述的等效分数阶模式涡旋电磁波产生及成像方法,其特征在于,所述第二步中:当发送整数阶模式为α时,设置接收模式为α+1或α-1,得到等效分数阶模式α+0.5或α-0.5。
3.根据权利要求1所述的等效分数阶模式涡旋电磁波产生及成像方法,其特征在于,所述第三步中:成像模式采用多发多收方式,发射信号的中心频率为fc,带宽为B,当产生的整数阶模式为-L~L时,结合分数阶模式,利用得到的-L,-L+0.5,-L+1,……L-1,L-0.5,L模式,通过二维快速傅立叶变换得到二维图像。
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 CN110988868A (zh) | 2020-04-10 |
CN110988868B true 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) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112505689A (zh) * | 2020-08-21 | 2021-03-16 | 西安空间无线电技术研究所 | 一种电磁涡旋雷达同心圆环阵列设计方法 |
CN112068129A (zh) * | 2020-08-21 | 2020-12-11 | 西安空间无线电技术研究所 | 一种同心圆环阵列电磁涡旋雷达成像方法 |
CN112083430B (zh) * | 2020-09-10 | 2023-11-07 | 天津水聿方舟海洋工程技术有限公司 | 一种适用于轨道角动量三维成像声呐的旁瓣抑制方法 |
CN112764041B (zh) * | 2021-01-19 | 2023-05-26 | 鹏城实验室 | 成像系统和方法 |
CN117706490B (zh) * | 2024-02-06 | 2024-05-17 | 南京理工大学 | 一种基于单站雷达的金属目标间的耦合散射中心建模方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103954367A (zh) * | 2014-04-18 | 2014-07-30 | 河南科技大学 | 一种测量分数阶光学涡旋拓扑荷值的装置及其测量方法 |
CN105445943A (zh) * | 2015-12-24 | 2016-03-30 | 河南科技大学 | 一种分数阶完美涡旋光束的产生装置及产生方法 |
CN106058490A (zh) * | 2016-06-01 | 2016-10-26 | 电子科技大学 | 一种产生涡旋电磁波的方法 |
CN106059675A (zh) * | 2016-05-24 | 2016-10-26 | 陕西理工学院 | 圆阵列旋转分形嵌套高阶涡旋电磁波生成分离方法与装置 |
CN107238933A (zh) * | 2017-07-07 | 2017-10-10 | 苏州大学 | 产生局部分数阶贝塞尔涡旋光束的方法和系统 |
CN109616777A (zh) * | 2018-12-28 | 2019-04-12 | 西安电子科技大学 | 基于时间调制阵列的任意模态涡旋电磁波产生方法 |
-
2019
- 2019-11-19 CN CN201911134517.9A patent/CN110988868B/zh active Active
Patent Citations (6)
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 | 西安电子科技大学 | 基于时间调制阵列的任意模态涡旋电磁波产生方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110988868A (zh) | 2020-04-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110988868B (zh) | 一种等效分数阶模式涡旋电磁波产生及成像方法 | |
CN106886020B (zh) | 一种单天线接收条件下的电磁涡旋成像方法 | |
CN105785328B (zh) | 基于子阵划分的fda距离-角度解耦合波束形成方法 | |
CN106486784A (zh) | 反射式天线阵列及波束扫描方法 | |
CN107290728B (zh) | 一种等效电磁波轨道角动量脉冲雷达探测方法及系统 | |
CN103630887B (zh) | 一种多波束形成方法及使用该方法的多波束声纳 | |
EP0807992B1 (en) | Logarithmic spiral array | |
CN105929384B (zh) | 基于十字阵发射圆环阵接收的联合低旁瓣多波束形成方法 | |
CN110058218B (zh) | 一种基于四维天线阵的射频隐身发射波束形成方法及系统 | |
US3877033A (en) | Nonuniformly optimally spaced array with uniform amplitudes | |
CN106324602B (zh) | 一种mimo声呐系统 | |
CN110850383A (zh) | 一种基于共形阵的mimo雷达信号处理方法 | |
CN111007504A (zh) | 基于最小冗余的mimo三维成像雷达稀疏阵列设计方法 | |
CN104808178B (zh) | 一种机载雷达发射方向图设计方法 | |
Deng et al. | Frequency diverse arc array beampattern synthesis analysis with nonlinear frequency offset | |
CN113281743A (zh) | 一种快速高分辨的基于oam雷达成像系统的成像方法 | |
Li et al. | Fast sparse synthesis of large-scale arrays in wideband and wide angle scanning | |
Rahaman et al. | Performane analysis of linearly-arranged concentric circular antenna array using robust ODL technique | |
Xuejun et al. | Design of sparse concentric multi-circular array using micro-genetic optimum algorithm | |
CN112926261B (zh) | 一种相控阵列天线副瓣快速计算方法及降副瓣方法 | |
Ren et al. | A Novel Optimal Subarray Partition Method for Large Uniform Linear Array based on Hierarchical Clustering Algorithm | |
Graham et al. | Radar architecture using MIMO transmit subarrays | |
He et al. | Analysis of over-the-horizon target detection based on the small aperture array | |
He et al. | Design of Sparse MIMO Array for Imaging Fuze | |
Ren et al. | Investigation on grating lobe suppression of uniform antenna array |
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 |