CN110221293B - 一种基于涡旋电磁波的俯仰角成像方法 - Google Patents
一种基于涡旋电磁波的俯仰角成像方法 Download PDFInfo
- Publication number
- CN110221293B CN110221293B CN201910387874.XA CN201910387874A CN110221293B CN 110221293 B CN110221293 B CN 110221293B CN 201910387874 A CN201910387874 A CN 201910387874A CN 110221293 B CN110221293 B CN 110221293B
- Authority
- CN
- China
- Prior art keywords
- target
- pitch angle
- antenna
- theta
- imaging
- 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.)
- Expired - Fee Related
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 39
- 238000000034 method Methods 0.000 claims abstract description 28
- 239000013598 vector Substances 0.000 claims description 10
- 238000001228 spectrum Methods 0.000 claims description 5
- 230000010363 phase shift Effects 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 2
- 238000005259 measurement Methods 0.000 abstract description 4
- 230000005284 excitation Effects 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 2
- 230000007123 defense Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
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)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供一种基于涡旋电磁波的俯仰角成像方法,属于雷达目标成像领域。传统方法是基于电磁波幅度扫描,其幅度扫描范围受到天线工作带宽的限制,与传统基于涡旋电磁波的俯仰角成像不同,本发明基于回波信号的相位信息,在单个频点即可工作,因此本发明的俯仰角成像范围广,并且对天线带宽要求较低。传统方法幅度扫描的精度取决于电磁波束的宽度,电磁波束越窄,成像精度越高,然而较窄的电磁波束难以实现。本方法的成像精度取决于对俯仰角离散化时的网格粗细,采用较细的网格可以有效提高成像精度。本发明能实现5°~85°的俯仰角测量,测量范围广,并且在俯仰角5°~85°的范围内成像误差低于1.2°。
Description
技术领域
本发明属于雷达目标成像领域,具体涉及一种基于涡旋电磁波的俯仰角成像方法。
背景技术
目前,雷达成像主要是基于电磁波的幅度、频率和相位三个特征进行的调幅、调频和调相信息调制,其主要利用的是远场平面波近似。而在1992年,荷兰物理学家L.Allen发现Laguerre-Gaussian激光束携带轨道角动量(OAM),即Laguerre-Gaussian激光束的波前不再是平面,而是绕着传播方向旋转,形成涡旋状。因此,涡旋电磁波是一种携带轨道角动量的电磁波,其相位波前为螺旋形状。由于涡旋电磁波为信息调制提供了额外的自由度,即在传统的幅度、频率和相位的基础上,增加了轨道角动量,有望应用在雷达和通信领域,增强雷达信息多样性,提高信道容量。
2013年,国防科技大学的郭桂蓉教授(郭桂蓉,胡卫东,杜小勇.基于电磁涡旋的雷达目标成像[J].国防科技大学学报,2013,35(6):71-76.)利用环形天线阵列产生涡旋电磁波,建立了涡旋电磁波照射下的理想散射点回波模型,分析了轨道角动量拓扑荷α与目标方位角之间的近似对偶关系,通过在轨道角动量域对目标回波信号进行采样,求解目标所处的方位角。该方法为新体系雷达设计和目标识别技术提供了参考。该理论基础为:
如图1所示,N个天线阵元在半径为a的圆周上排列,以圆心o为原点、圆周法向为z轴建立坐标系。将目标当做理想散射点,假设目标所处位置在极坐标系下表示为其中r为目标到原点的距离,θ表示目标所处的俯仰角,表示目标所处的方位角。第n个天线阵元的位置矢量为其中和分别为x方向和y方向的单位矢量,相应的,为z方向的单位矢量,φn=2nπ/N为第n个天线阵元所处的方位角。给每个天线阵元附加的相移,对于空间中某一点其场强可以写为每个天线阵元在该点的辐射场强之和。
其中,μ0为真空磁导率,ω为电磁波角频率,k=ω/c为波束,c为真空中的光速,j为每个天线阵元出的电流密度,d为天线阵元振子长度,∫(·)dVn表示对阵元上的电流元积分。式(1)中利用了当阵元数充分大时的近似,α为轨道角动量拓扑荷,为Bessel函数。式(1)还利用了幅度近似|r-rn|≈r和相位近似其中为场点矢径的单位向量,r为场点的位置矢量,考虑理想的点散射行为,记P点的散射强度为第n个天线阵元接收的回波可以写为:
式(4)中出现了项,而其他项均与无关,此时轨道角动量拓扑荷α与第m(m=1,2,...,M)号目标方位角构成了傅里叶变换对,即可以在轨交角动量域采样(即收集不同α时的信号回波),通过谱估计方法(如傅里叶变换、压缩感知等)估计出目标的方位角。
2017年,浙江大学的武华阳在他的硕士论文中给出了一种基于涡旋电磁波的俯仰角成像方法(无线轨道角动量通信与雷达目标成像技术研究[D].2017.),拓展了涡旋电磁波雷达成像的成像范围。该方法同样采用如图1的成像系统,且推导方法相同,得出环形天线阵列产生的涡旋电磁波空间分布公式(1)后,从中提取信息:涡旋电磁波在θ面内的方向图与贝塞尔函数有关,即Jα(kasinθ),此时涡旋电磁波的最大方向指向对应的俯仰角为
其中f为激励信号频率,当f变化时,θmax也会变化,如图2所示,因此可以通过改变涡旋电磁波的激励信号频率,从而改变涡旋电磁波的最大增益方向,继而实现涡旋电磁波在俯仰角范围内的扫描。该方法的缺点在于:(1)如图2所示,当激励信号频率高于2.4GHz或低于1GHz时,最大增益角度θmax的变化趋势已经非常缓慢,且现实中的天线频带宽度有限,无法实现较大跨度的激励频率发生,因此无法实现较宽的扫描范围。(2)如图3所示,当激励频率较小时,虽然能够使涡旋电磁波的最大增益方向朝向较大的俯仰角,但是波束的主瓣宽度也相应增加,导致波束的方向性降低,不利于精确成像;
发明内容
针对背景技术所存在的问题,本发明的目的在于提供一种基于涡旋电磁波的俯仰角成像方法,该方法在涡旋电磁波回波信号模型中构造与目标俯仰角相关的傅里叶变换对,通过谱估计方法计算目标的俯仰角,该方法精确度更高,俯仰角成像范围更广。
为实现上述目的,本发明的技术方案如下:
一种基于涡旋电磁波的俯仰角成像方法,包括以下步骤:
步骤2:步骤1所述涡旋电磁波照射雷达目标,目标对来波产生散射,散射波由均匀环形天线阵列接收;
步骤4:对天线阵列接收到的回波信号进行谱估计:
步骤5:通过Θ(θ)值计算出目标俯仰角θ:目标俯仰角θ的取值范围为0:π/2,而sinθ在θ=0:π/2范围内是单值函数,即Θ(θ)唯一对应一个θ值。
进一步地,步骤1所述雷达天线具有收发功能。
进一步地,步骤1所述天线阵元为喇叭天线、贴片天线或偶极子天线等。
进一步地,步骤6所述谱估计方法为傅里叶变换法或压缩感知法等。
综上所述,由于采用了上述技术方案,本发明的有益效果是:传统方法是基于电磁波幅度扫描,其幅度扫描范围受到天线工作带宽的限制,与传统基于涡旋电磁波的俯仰角成像不同,本发明基于回波信号的相位信息,在单个频点即可工作,因此本发明的俯仰角成像范围广,并且对天线带宽要求较低。传统方法幅度扫描的精度取决于电磁波束的宽度,电磁波束越窄,成像精度越高,然而较窄的电磁波束难以实现。本方法的成像精度取决于对俯仰角离散化时的网格粗细,采用较细的网格可以有效提高成像精度。本发明能实现5°~85°的俯仰角测量,测量范围广,并且在俯仰角5°~85°的范围内成像误差低于1.2°。
附图说明
图1为涡旋电磁波照射下的理想散射点的回波信号模型。
图2为涡旋电磁波最大增益角度θmax随激励频率f的变化曲线图。
图3为不同激励信号频率下涡旋电磁波在俯仰角θ平面内的方向图。
图4为本发明的俯仰角成像结果与目标实际俯仰角对比图。
图5为本发明的俯仰角成像结果与目标实际俯仰角之间的误差图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面结合实施方式和附图,对本发明作进一步地详细描述。
本发明使用仿真软件FEKO,并在软件中建立涡旋电磁波成像模型。本发明中的天线阵元可以是喇叭天线、贴片天线等,仿真中用理想偶极子天线作为天线阵元,假设该目标所处的方位角已知。
一种基于涡旋电磁波的俯仰角成像方法,包括以下步骤:
步骤1:设置雷达天线:雷达由N=30个天线阵元在半径为a=50mm的圆周上组成的均匀环形天线阵列构成,以圆心为坐标原点、圆周法向为z轴建立坐标系,给每个天线阵元附加的相移,使其产生轨道角动量拓扑荷为α的涡旋电磁波,α可以选择除零以外的任意值,此处选择α=1;
步骤2:步骤1所述涡旋电磁波照射目标时产生散射,得到的散射波由均匀环形天线阵列接收;
步骤3:计算天线阵列接收的散射回波信号:
步骤4:通过压缩感知谱估计方法计算得出俯仰角Θ(θ):
将俯仰角域(0,π/2)离散成Q个点,两点之间的距离为π/(2Q-1),此时我们得到了Θ(θ)的离散值,记为Θq(θq)(q=1,2,...,Q)。可以将式(11)写为矩阵形式
其中Φn(n=0,1,...,N-1)为的简写,xq(q=1,2,...,Q)为待求系数,式(12)可以通过正交匹配追踪算法求解。式(11)是针对只有一个理想散射点的情况推导的,因此式(12)求解出的xq(q=1,2,...,Q)中只有一个值为非零值,其余全部为零,该非零值唯一对应一个Θq(θq)。目标俯仰角θ的取值范围为0:π/2,而sinθ在θ=0:π/2范围内是单值函数,即Θq(θq)也唯一对应一个θq,该θq即为目标所处的俯仰角。
仿真过程中使目标从俯仰角θ=5°到θ=85°范围内移动,每间隔为5°,分别对目标进行俯仰角成像。成像结果如图4所示,可见在俯仰角θ=5°到θ=85°范围内,成像结果与实际结果吻合良好。成像结果和目标实际俯仰角之间的差值如图5所示,可见在俯仰角θ=5°到θ=85°范围内,成像误差低于1.2°。
以上所述,仅为本发明的具体实施方式,本说明书中所公开的任一特征,除非特别叙述,均可被其他等效或具有类似目的的替代特征加以替换;所公开的所有特征、或所有方法或过程中的步骤,除了互相排斥的特征和/或步骤以外,均可以任何方式组合。
Claims (4)
1.一种基于涡旋电磁波的俯仰角成像方法,其特征在于,包括以下步骤:
步骤2:获取回波信息:步骤1所述涡旋电磁波照射目标产生散射,散射波由均匀环形天线阵列接收;
步骤3:根据回波信息求解目标的俯仰角,具体步骤如下:
其中,i为虚数单位,r为目标到原点的距离,θ表示目标所处的俯仰角,表示目标所处的方位角,rn为第n个天线阵元的位置矢量,φn=2nπ/N为第n个天线阵元所处的方位角,μ0为真空磁导率,ω为电磁波角频率,k=ω/c为波束,c为真空中的光速,j为每个天线阵元出的电流密度,d为天线阵元振子长度,∫(·)dVn表示对阵元上的电流元积分,为Bessel函数,r为场点的位置矢量, 为场点矢径的单位向量;
步骤3.2:对天线阵列接收到的回波信号进行谱估计计算得到Θ(θ)的值:
步骤3.3:根据Θ(θ)值计算出目标俯仰角θ。
2.一种如权利要求1所述的基于涡旋电磁波的俯仰角成像方法,其特征在于,步骤1所述雷达天线具有收发功能。
3.一种如权利要求1所述的基于涡旋电磁波的俯仰角成像方法,其特征在于,步骤1所述天线阵元为喇叭天线、贴片天线或偶极子天线。
4.一种如权利要求1所述的基于涡旋电磁波的俯仰角成像方法,其特征在于,步骤3.2所述谱估计方法为傅里叶变换法或压缩感知法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910387874.XA CN110221293B (zh) | 2019-05-10 | 2019-05-10 | 一种基于涡旋电磁波的俯仰角成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910387874.XA CN110221293B (zh) | 2019-05-10 | 2019-05-10 | 一种基于涡旋电磁波的俯仰角成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110221293A CN110221293A (zh) | 2019-09-10 |
CN110221293B true CN110221293B (zh) | 2022-11-08 |
Family
ID=67820659
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910387874.XA Expired - Fee Related CN110221293B (zh) | 2019-05-10 | 2019-05-10 | 一种基于涡旋电磁波的俯仰角成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110221293B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110824433B (zh) * | 2019-10-09 | 2021-11-23 | 清华大学 | 电磁波量子态轨道角动量雷达探测和方法 |
CN110967692B (zh) * | 2019-12-24 | 2023-04-11 | 上海无线电设备研究所 | 一种成像方法 |
CN111474543B (zh) * | 2020-04-13 | 2022-08-09 | 中国人民解放军国防科技大学 | 一种涡旋电磁波干涉的目标三维成像方法及装置 |
CN112751197B (zh) * | 2020-12-25 | 2021-11-09 | 电子科技大学 | 一种相控涡旋电磁波产生系统及方法 |
CN112327299B (zh) * | 2021-01-04 | 2021-03-16 | 中国人民解放军国防科技大学 | 电磁涡旋前视目标三维成像方法 |
CN112904347B (zh) * | 2021-01-19 | 2023-05-26 | 鹏城实验室 | 成像系统和方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109490881A (zh) * | 2018-11-21 | 2019-03-19 | 中国科学院电子学研究所 | 基于涡旋电磁波的干涉sar高程测量系统及测量方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5115243A (en) * | 1991-04-16 | 1992-05-19 | General Electric Co. | Radar system with active array antenna, beam multiplex control and pulse integration control responsive to azimuth angle |
US10698101B2 (en) * | 2015-04-14 | 2020-06-30 | Northeastern University | Compressive coded antenna/meta-antenna |
DE102016005910B4 (de) * | 2016-05-17 | 2024-02-15 | Rheinmetall Air Defence Ag | Antennenanordnung eines Lenkflugkörpers mit einer Radarantenne |
CN106526589B (zh) * | 2016-12-26 | 2018-11-30 | 浙江大学 | 一种基于涡旋电磁波的雷达目标二维成像方法 |
CN107394408A (zh) * | 2017-07-24 | 2017-11-24 | 电子科技大学 | 基于全息张量阻抗表面的涡旋电磁波束天线及其实现方法 |
CN108594221A (zh) * | 2018-04-27 | 2018-09-28 | 中国人民解放军国防科技大学 | 基于同心圆环阵列的涡旋电磁波产生与优化方法 |
-
2019
- 2019-05-10 CN CN201910387874.XA patent/CN110221293B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109490881A (zh) * | 2018-11-21 | 2019-03-19 | 中国科学院电子学研究所 | 基于涡旋电磁波的干涉sar高程测量系统及测量方法 |
Non-Patent Citations (1)
Title |
---|
基于轨道角动量的雷达关联成像技术研究;陈亚南等;《电子设计工程》;20180320(第06期);第115-119+125页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110221293A (zh) | 2019-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110221293B (zh) | 一种基于涡旋电磁波的俯仰角成像方法 | |
CN105044684B (zh) | 基于射频隐身的mimo跟踪雷达发射波束的形成方法 | |
JP2006270847A (ja) | アンテナ装置 | |
Lin et al. | Improved OAM-based radar targets detection using uniform concentric circular arrays | |
CN109143230B (zh) | 一种应用于s波段穿墙雷达的二维波达方向估计方法 | |
KR101807499B1 (ko) | 상호간섭효과를 보상한 원형 배열 안테나 기반의 도래각 추정 장치 및 추정 방법 | |
CN108089161B (zh) | 一种基于功率反馈的天线阵列合成波束空域选择方法 | |
Tsai et al. | Beam AoD-based indoor positioning for 60 GHz MmWave system | |
Alawsh et al. | Sparse direction of arrival estimation using sparse arrays based on software-defined-radio platform | |
RU2713503C1 (ru) | Способ углового сверхразрешения в приемных цифровых антенных решетках | |
Deng et al. | Frequency diverse arc array beampattern synthesis analysis with nonlinear frequency offset | |
Wang et al. | A practical calibration method of linear UHF Yagi arrays for ship target detection application | |
Shi et al. | Wideband DOA estimation based on A-shaped array | |
Liu et al. | Beamforming of joint polarization-space matched filtering for conformal array | |
Dmitriev et al. | Synthesis of Directivity Pattern for Various Antenna Arrays Configurations of Vertical Atmospheric Sensing Radar Station | |
CN111736111A (zh) | 一种基于同心均匀分布双圆阵的相干源doa的估计方法 | |
Yang et al. | Direction of arrival estimation on cylindrical conformal array using RARE | |
Yunhe et al. | Direction finding for bistatic MIMO radar with uniform circular array | |
Khedekar et al. | Digital beamforming to reduce antenna side lobes and minimize DOA error | |
Zhang et al. | Study on DOA estimation algorithm for smart antenna of uniform circular | |
Zhang et al. | Intermediate frequency digital receiver based on multi-FPGA system | |
Peshkov et al. | Choice of Antenna Array Geometry for Digital Spatial Filtering in Azimuth and Elevation | |
Alsalti et al. | Direction of arrival for uniform circular array using directional antenna elements | |
RU2746063C1 (ru) | Способ углового сверхразрешения в приемной цифровой антенной решётке | |
Mei et al. | Joint Estimation of DOA and Polarization Based on Cylindrical Conformal Array Antenna |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20221108 |