CN109946642B - 涡旋电磁波的到达角估计方法 - Google Patents
涡旋电磁波的到达角估计方法 Download PDFInfo
- Publication number
- CN109946642B CN109946642B CN201910173063.XA CN201910173063A CN109946642B CN 109946642 B CN109946642 B CN 109946642B CN 201910173063 A CN201910173063 A CN 201910173063A CN 109946642 B CN109946642 B CN 109946642B
- Authority
- CN
- China
- Prior art keywords
- transmitting end
- angle
- coordinate system
- transmitting
- theta
- 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
- 238000000034 method Methods 0.000 title claims abstract description 31
- 239000011159 matrix material Substances 0.000 claims abstract description 40
- 238000004891 communication Methods 0.000 claims abstract description 31
- 238000003491 array Methods 0.000 claims abstract description 5
- 239000013598 vector Substances 0.000 claims description 22
- 230000005684 electric field Effects 0.000 claims description 3
- 238000012549 training Methods 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims 1
- 239000000969 carrier Substances 0.000 abstract 1
- 238000001228 spectrum Methods 0.000 description 8
- 238000004088 simulation Methods 0.000 description 4
- 238000000354 decomposition reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000002441 reversible effect Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 241000764238 Isis Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种涡旋电磁波到达角估计方法,主要解决现有技术无法直接估计多载波‑多模态OAM通信***发送端到达角的问题。其实现方案是:利用两个具有相同参数的均匀圆形天线阵构成通信***的收发端,发送端利用多个载波及多个模态发送已知的信息信号;在收发端分别建立直角坐标系,确定收发端坐标系间的角度关系,得到接收信号在发送端坐标系下的信号形式并构成接收信号矩阵;基于2‑D ESPRIT算法及二分法利用接收信号矩阵估计发送端的方位角及俯仰角,即得到发送端所发送涡旋电磁波的到达角。本发明能精确估计出多载波‑多模态OAM通信***发送端的到达角,可用于涡旋电磁波无线通信。
Description
技术领域
本发明属于无线通信技术领域,具体涉及一种到达角估计方法,可用于涡旋电磁波无线通信。
背景技术
现阶段,通信技术的发展导致了无线频谱资源需求的急速上升,然而,传统对频谱资源的控制导致了频谱利用率的不足,频谱资源匮乏的问题越来越严重。电磁波属性中的振幅、频率、相位和偏振态等维度均已用于信号表征来提升传输容量。在现有基础上无法继续采用增加电磁波表征维度的方式来扩充信道容量,只能通过频谱压缩、提高调制速率或者调制阶数等方法来进一步提高频谱效率。轨道角动量OAM作为一个新的传输维度,能够在同一频带同时传输多路信息,可以有效地解决频谱资源短缺的问题,目前已引起了广泛关注。
实现OAM通信的关键技术之一就是如何获取涡旋电磁波的到达角。近年来,人们对涡旋电磁波的方位角估算方法进行了大量的研究。一些用来估计OAM波束方位角的方法已经提出,比如傅里叶变换算法、逆投影算法、MUSIC算法及ESPRIT算法。然而,傅里叶变换算法及逆投影算法需要利用大量OAM模态才能较为准确的估计涡旋电磁波的方位角。使用MUSIC算法及ESPRIT算法可利用相对少的OAM模态实现涡旋电磁波方位角的超分辨估计,但上述方法均无法估计涡旋电磁波的俯仰角,无法实现涡旋电磁波到达角的估计,使得涡旋电磁波在无线通信的实际应用受到局限,影响OAM通信***的通信质量。
发明内容
本发明的目的在于针对上述现有技术的不足,提出一种涡旋电磁波的到达角估计方法,改善OAM通信***的通信质量。
本发明的技术方案是:利用2-D ESPRIT算法及二分法来估计涡旋电磁波到达角的方法,以进行OAM通信***收发天线阵列的对准,提高通信***容量,其实现步骤包括如下:
(2)在发送端与接收端分别建立坐标系,将接收端圆心在发送端坐标中的坐标表示为(r,α,π-θ),其中r为发送端圆心与接收端圆心间的距离,α为发送端坐标系与接收端坐标系X轴间的夹角,θ为发送端的俯仰角,将接收端接收到的信号在发送端坐标系中表示为其中kp为波数,lq为所用模态;
(4)提取接收信号矩阵中各元素的辐角,构成辐角矩阵X,将辐角矩阵X的各列向量依次排列成向量X1,将辐角矩阵X的各行向量依次排列成向量X2,基于2-D ESPRIT算法利用这两个向量X1及X2信号子空间的旋转不变性,分别估计发送端圆心与接收端圆心间的距离r及发送端坐标系与接收端坐标系X轴间的夹角α;
本发明具有以下优点:
1.本发明通过观察OAM通信***接收信号的信号形式,确定了接收信号辐角与幅值的估计为两个独立过程,即利用2-D ESPRIT算法与二分法分别估计接收信号的辐角与幅值,解决了现有技术无法直接估计OAM通信***发送端方位角与俯仰角的问题。
2.本发明利用超分辨2-D ESPRIT算法估计接收信号的辐角,并在利用二分法估计接收信号幅值的过程中设置非常高的估算精度,因而能够实现发送端到达角的高精度估计。
3.本发明在估计发送端到达角的过程中运用了频率分集及OAM模态分集,故与传统多输入多输出MIMO通信***相比,本发明能够利用极少的时间采样、较少的子载波实现发送端到达角的估计,能够有效的节省时域及频谱资源。
附图说明
图1是本发明的实现流程图;
图2是本发明中发送端与接收端坐标系间的角度转化关系;
图3是本发明使用的OAM通信***模型示意图;
图4是在不同信噪比条件下利用本发明所提出的方法得到的发送端到达角的估计图。
具体实施方式
以下参照附图对本发明实施例和效果作进一步详细描述。
参照图1,本实施例的实现步骤如下:
步骤1,确定OAM通信***接收端圆心坐标。
由于结构简单且易于复用多种模态的涡旋电磁波,本步骤中的OAM通信***用均匀圆形天线阵列UCA作为***的发送端与接收端。在实际通信中,接收端与发送端的正对准很难实现,为便于分析,本步骤考虑了非平行错位情况。
参照图2,以发送端圆心为原点O,与接收端UCA平行的平面为XOY平面建立坐标系Z-XOY,以发送端圆心为原点O,发送端UCA所在平面为平面建立直角坐标系以接收端圆心为原点O’,接收端UCA所在平面为X’O’Y’平面建立直角坐标系Z’-X’O’Y’,由图2可知,轴与Z轴间的夹角为根据三余弦定理可知,X轴和轴间的夹角可表示为:
步骤2,构建接收信号矩阵。
参照图3,本实例中OAM通信***使用基于UCA的P路子载波-Q维模态多路复用通信***,其收发天线阵均采用阵元为N、半径为a的均匀圆形天线阵,根据常用OAM通信***模型,收发天线阵元的个数N取9;
其中,j是每个各阵元处的电流密度,d是电偶极子的长度,μ0为真空中的磁导率,ω为角频率,N为收发端UCA阵元数目,k=2πf/c为波数,f为频率,c为真空中的光速,l为模态,i为虚数单位,a为收发端阵列半径,为点P的位置矢量,r0为点P到点O的距离,为点P在Z-XOY坐标系下的方位角,θ0点P在Z-XOY坐标系下的俯仰角,Jl(·)为第一类贝塞尔函数;
当发送端在训练阶段发送已知的信息信号s(kp,lq)时,接收端所接收到的信号在Z-XOY坐标系下的表达示为:
由于发送端利用P路子载波及Q个模态产生P×Q个已知信息信号,因此,接收端将接收到P×Q个相互正交的接收信号,其构成的接收信号矩阵:
其中,kp+1-kp=1,lq+1-lq=1,p=1,2,…,P,q=1,2,…,Q,即发送信号所使用的波数kp及模态lq需为连续整数。
步骤3,利用接收信号矩阵构建辐角矩阵。
其中,(·)w表示共轭复数,sign(·)为符号函数;
步骤4,估计发送端圆心与接收端圆心间的距离r。
4.1)将式<3>中辐角矩阵X的各列向量依次排列构成向量X1:
X1=[x(k1,l1),x(k1,l2),…,x(kp,lq),…,x(kP,lQ)]T,
由于X1可分解为X1=A1+N1,则X1的协方差矩阵可表示为:
其中,Ψ=T-1ΦT为关系矩阵,Φ为Ψ的特征值;
X2=[x(k1,l1),x(k2,l1),…,x(kp,lq),…,x(kP,lQ)]T,
由于X2可分解为X2=A2+N2,则X2的协方差矩阵可表示为:
其中,∑=M-1δM为转化矩阵,δ为∑的特征值;
步骤6,估计发送端方位角及俯仰角。
根据实际情况,仅当接收端位于发送端主瓣波束范围[θa,θb]内时,才能够接收到清晰的信号,因此,发送端俯仰角需满足关系θ∈[θa,θb],此外,由于本方案所用场景为一对一通信,即当θ∈[θa,θb]时,零点函数有且仅有一个零点,且因此,可利用二分法估算零点函数的零点,将某一子载波kp及模态lq下零点函数的零点作为发送端俯仰角的估计值其具体步骤如下:
6.2a)求发送波束主瓣范围[θa,θb]的中点θc,给定估计精度ξ;
本发明的效果可以通过以下仿真结果进一步说明:
1.仿真条件:
取收发端UCA阵元个数N为9,发送信号所用波数为k=189,190,…,196,所用模态为l=-4,-3,…,0,…,3,二分法估计精度ξ=0.0001,发送端圆心与接收端圆心距离r为40m,发送端的方位角为7度,发送端的俯仰角θ为7度。
2.仿真内容:
仿真1,在信噪比为5dB,10dB,15dB及20dB的条件下,利用2-D ESPRIT算法及二分法估算发送端的方位角及俯仰角,结果如图4所示。
从图4表明,随着信噪比的增加,利用本发明方法估计得到的方位角及俯仰角逐渐趋近于真实值,当信噪比为20dB时,发送端方位角的估计值为7.000度,发送端俯仰角的估计值为6.993度,与发送端方位角与俯仰角的真实值非常接近。可见,利用本发明进行发送端方位角与俯仰角θ的估计能够很好的解决OAM通信***涡旋电磁波到达角估计的问题。
综上,本实例能够精确估计OAM通信***发送端的到达角,为OAM通信***的实际商用奠定了基础。
Claims (7)
1.一种涡旋电磁波的到达角估计方法,其特征在于,包括如下:
(1)利用两个具有相同参数的均匀圆阵构成轨道角动量通信***的发送端与接收端,将发送端在远场空间中任意一点所产生的电场强度表示为发送端在训练阶段发送已知的信息信号;其中,为点P的位置矢量,为点P到发送端圆心的距离,为点P的方位角,θ0点P的俯仰角;
(2)在发送端与接收端分别建立坐标系,将接收端圆心在发送端坐标中的坐标表示为(r,α,π-θ),其中r为发送端圆心与接收端圆心间的距离,α为发送端坐标系与接收端坐标系X轴间的夹角,θ为发送端的俯仰角,将接收端接收到的信号在发送端坐标系中表示为其中kp为波数,lq为所用模态;
(4)提取接收信号矩阵中各元素的辐角,构成辐角矩阵X,将辐角矩阵X的各列向量依次排列成向量X1,将辐角矩阵X的各行向量依次排列成向量X2,基于2-D ESPRIT算法利用这两个向量X1及X2信号子空间的旋转不变性,分别估计发送端圆心与接收端圆心间的距离r及发送端坐标系与接收端坐标系X轴间的夹角α;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910173063.XA CN109946642B (zh) | 2019-03-07 | 2019-03-07 | 涡旋电磁波的到达角估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910173063.XA CN109946642B (zh) | 2019-03-07 | 2019-03-07 | 涡旋电磁波的到达角估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109946642A CN109946642A (zh) | 2019-06-28 |
CN109946642B true CN109946642B (zh) | 2023-04-07 |
Family
ID=67008466
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910173063.XA Active CN109946642B (zh) | 2019-03-07 | 2019-03-07 | 涡旋电磁波的到达角估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109946642B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2022183469A1 (en) * | 2021-03-05 | 2022-09-09 | Qualcomm Incorporated | Multi-mode precoding matrix information report for orbital angular momentum based communication system |
CN113660012B (zh) * | 2021-09-28 | 2022-09-23 | 南京航空航天大学 | 基于随机频率分集阵列的距离敏感局部涡旋波产生方法 |
CN115015912B (zh) * | 2022-02-24 | 2024-07-19 | 西安电子科技大学 | 一种基于涡旋电磁波的旋转目标空间角估计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106526589A (zh) * | 2016-12-26 | 2017-03-22 | 浙江大学 | 一种基于涡旋电磁波的雷达目标二维成像方法 |
WO2017062902A1 (en) * | 2015-10-08 | 2017-04-13 | Polte Corporation | Angle of arrival positioning system for tracking objects |
JP2018037744A (ja) * | 2016-08-29 | 2018-03-08 | 日本電信電話株式会社 | アンテナ調整方法及び通信システム |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9998187B2 (en) * | 2014-10-13 | 2018-06-12 | Nxgen Partners Ip, Llc | System and method for combining MIMO and mode-division multiplexing |
-
2019
- 2019-03-07 CN CN201910173063.XA patent/CN109946642B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017062902A1 (en) * | 2015-10-08 | 2017-04-13 | Polte Corporation | Angle of arrival positioning system for tracking objects |
JP2018037744A (ja) * | 2016-08-29 | 2018-03-08 | 日本電信電話株式会社 | アンテナ調整方法及び通信システム |
CN106526589A (zh) * | 2016-12-26 | 2017-03-22 | 浙江大学 | 一种基于涡旋电磁波的雷达目标二维成像方法 |
Non-Patent Citations (3)
Title |
---|
ORBITAL ANGULAR MOMENTUM-BASED TWO-DIMENSIONAL SUPER-RESOLUTION TARGETS IMAGING;Rui Chen;<2018 IEEE Global Conference on Signal and Information Processing (GlobalSIP)>;20190221;全文 * |
涡旋电磁波及其在雷达中应用研究进展;刘康等;《电子学报》;20180915(第09期);全文 * |
涡旋电磁波的理论与应用研究进展;李龙等;《微波学报》;20180419(第02期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109946642A (zh) | 2019-06-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109946642B (zh) | 涡旋电磁波的到达角估计方法 | |
CN110266354B (zh) | 轨道角动量无线通信***的信号接收方法 | |
CN113114311B (zh) | 基于智能反射面和发射端的联合波束赋形与空间调制方法 | |
CN105319545B (zh) | 提高stap检测性能的mimo-ofdm雷达波形设计方法 | |
WO2021227482A1 (zh) | 一种大规模天线***中的安全传输方法 | |
CN111181671B (zh) | 一种基于深度学习的下行信道快速重建方法 | |
CN104833959A (zh) | 基于目标先验信息的mimo雷达波形优化方法 | |
CN108667493B (zh) | 一种面向大规模mimo中nlos场景下的波束赋形方法 | |
CN106209183B (zh) | 面向同心圆或同轴圆台均匀圆阵列天线的涡旋电磁波mimo***的涡旋信号分离方法与装置 | |
CN106998307B (zh) | 一种用于大规模天线***的盲信号检测及信道估计方法 | |
Wu et al. | Reconfigurable hybrid beamforming for dual-polarized mmWave MIMO channels: Stochastic channel modeling and architectural adaptation methods | |
CN105487054A (zh) | 提高基于mimo-ofdm雷达stap最差检测性能的稳健波形设计方法 | |
Jiang et al. | Hybrid beamforming for massive MIMO: A unified solution for both phase shifter and switch networks | |
CN104506256B (zh) | 一种mimo多天线***性能评估方法及多天线*** | |
CN108572347A (zh) | 基于通信信号信道状态响应信息的面阵二维测角方法及*** | |
CN104994041B (zh) | 三维大规模天线***基于导频的波达波离角联合估计方法 | |
CN110113281B (zh) | 一种mimo通信中多进制fsk非相干检测实现空分复用的方法 | |
CN114244658B (zh) | 大规模mimo***中的基于倍角估计的信道估计方法 | |
CN114154347B (zh) | 基于admm的双功能mimo雷达通信***波形设计方法 | |
CN114338303B (zh) | 大规模MIMO***中基于多维Hankel矩阵的信道估计方法及*** | |
CN110213186A (zh) | 基于混合结构的esprit的到达角与信噪比估计方法 | |
Feng et al. | Hybrid precoding for massive MIMO systems using partially-connected phase shifter network | |
CN112188530B (zh) | 一种基于透镜天线阵列提高无线通信***安全能量效率的传输方法 | |
KR102111855B1 (ko) | 그룹배열안테나 기반의 mimo 송수신 방법, mimo 송신기 및 mimo 수신기 | |
Yang et al. | Efficient multi-beam training for terahertz wireless communications |
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 |