涡旋电磁波的到达角估计方法与流程

本发明属于无线通信技术领域，具体涉及一种到达角估计方法，可用于涡旋电磁波无线通信。

背景技术：

现阶段，通信技术的发展导致了无线频谱资源需求的急速上升，然而，传统对频谱资源的控制导致了频谱利用率的不足，频谱资源匮乏的问题越来越严重。电磁波属性中的振幅、频率、相位和偏振态等维度均已用于信号表征来提升传输容量。在现有基础上无法继续采用增加电磁波表征维度的方式来扩充信道容量，只能通过频谱压缩、提高调制速率或者调制阶数等方法来进一步提高频谱效率。轨道角动量oam作为一个新的传输维度，能够在同一频带同时传输多路信息，可以有效地解决频谱资源短缺的问题，目前已引起了广泛关注。

实现oam通信的关键技术之一就是如何获取涡旋电磁波的到达角。近年来，人们对涡旋电磁波的方位角估算方法进行了大量的研究。一些用来估计oam波束方位角的方法已经提出，比如傅里叶变换算法、逆投影算法、music算法及esprit算法。然而，傅里叶变换算法及逆投影算法需要利用大量oam模态才能较为准确的估计涡旋电磁波的方位角。使用music算法及esprit算法可利用相对少的oam模态实现涡旋电磁波方位角的超分辨估计，但上述方法均无法估计涡旋电磁波的俯仰角，无法实现涡旋电磁波到达角的估计，使得涡旋电磁波在无线通信的实际应用受到局限，影响oam通信系统的通信质量。

技术实现要素：

本发明的目的在于针对上述现有技术的不足，提出一种涡旋电磁波的到达角估计方法，改善oam通信系统的通信质量。

本发明的技术方案是：利用2-desprit算法及二分法来估计涡旋电磁波到达角的方法，以进行oam通信系统收发天线阵列的对准，提高通信系统容量，其实现步骤包括如下：

(1)利用两个具有相同参数的均匀圆阵构成轨道角动量通信系统的发送端与接收端，将发送端在远场空间中任意一点所产生的电场强度表示为发送端在训练阶段发送已知的信息信号；

(2)在发送端与接收端分别建立坐标系，将接收端圆心在发送端坐标中的坐标表示为(r，α，π-θ)，其中r为发送端圆心与接收端圆心间的距离，α为发送端坐标系与接收端坐标系x轴间的夹角，θ为发送端的俯仰角，将接收端接收到的信号在发送端坐标系中表示为其中kp为波数，lq为所用模态；

(3)利用构成接收信号矩阵

(4)提取接收信号矩阵中各元素的辐角，构成辐角矩阵x，将辐角矩阵x的各列向量依次排列成向量x1，将辐角矩阵x的各行向量依次排列成向量x2，基于2-desprit算法利用这两个向量x1及x2信号子空间的旋转不变性，分别估计发送端圆心与接收端圆心间的距离r及发送端坐标系与接收端坐标系x轴间的夹角α；

(5)提取接收信号的幅值

其中，是与发送端方位角及俯仰角θ估算过程无关的常数项，kp＝2πfp/c为波数，fp为所用频率，c为光速，a为收发端圆阵半径；

(6)将式<1>右边的各项移至式<1>的左边，将这一函数表示为利用二分法在发送端主瓣波束的俯仰角范围内估计的零点，估计得到发送端的俯仰角θ；

(7)利用发送端坐标系与接收端坐标系x轴间夹角α的估计值和发送端俯仰角θ的估计值，估计得到发送端的方位角

本发明具有以下优点：

1.本发明通过观察oam通信系统接收信号的信号形式，确定了接收信号辐角与幅值的估计为两个独立过程，即利用2-desprit算法与二分法分别估计接收信号的辐角与幅值，解决了现有技术无法直接估计oam通信系统发送端方位角与俯仰角的问题。

2.本发明利用超分辨2-desprit算法估计接收信号的辐角，并在利用二分法估计接收信号幅值的过程中设置非常高的估算精度，因而能够实现发送端到达角的高精度估计。

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轴和轴间的夹角可表示为：

其中，θ为发送端的俯仰角，为发送端的方位角。因此，o’在z-xoy坐标系下的坐标为(r，α，π-θ)，其中，r为发送端圆心与接收端圆心间的距离。

步骤2，构建接收信号矩阵。

参照图3，本实例中oam通信系统使用基于uca的p路子载波-q维模态多路复用通信系统，其收发天线阵均采用阵元为n、半径为a的均匀圆形天线阵，根据常用oam通信系统模型，收发天线阵元的个数n取9；

将发送端在远场空间中任意一点处所产生的电场强度在z-xoy坐标系下表示为：

其中，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坐标系下的表达示为：

其中，是与发送端方位角及俯仰角θ估算过程无关的常数项，为附加噪声，kp为所用波数，lq为所用模态。

由于发送端利用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，利用接收信号矩阵构建辐角矩阵。

3.1)提取接收信号矩阵中各元素的辐角：

其中，(·)w表示共轭复数，sign(·)为符号函数；

3.2)利用的辐角x(kp，lq)构成辐角矩阵：

步骤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的协方差矩阵可表示为：

其中，n1为附加噪声，为噪声方差；

4.2)将特征分解为：其中，q是pq×pq维的酉矩阵，为特征值构成的对角阵；

4.3)定义与λ1max所对应的特征向量为其中为x1的信号子空间，利用的第1至q×(p-1)个元素构成向量的第(q+1)至q×p个元素构成向量根据信号子空间的旋转不变性，得到：

其中，aa为a1的前(p-1)×q个元素构成的向量，ab为a1的后(p-1)×q个元素构成的向量，aa与ab满足ab＝aaφ，φ＝e^ir，t为可逆矩阵，且满足

4.4)根据式<4>，得到与间的关系式：

其中，ψ＝t^-1φt为关系矩阵，φ为ψ的特征值；

4.5)根据上述推导，利用及得到关系矩阵ψ：

其中(·)⁺表示矩阵的广义逆；

4.6)对关系矩阵ψ进行特征分解得到特征值φ，提取特征值φ的辐角作为收发端圆心距离r的估计值

步骤5，估计x轴和轴间的夹角α。

5.1)在估算x轴和轴间的夹角α时，将式<3>中辐角矩阵x的各行元素依次排列构成向量x2：

x2＝[x(k1，l1)，x(k2，l1)，…，x(kp，lq)，…，x(kp，lq)]^t，

由于x2可分解为x2＝a2+n2，则x2的协方差矩阵可表示为：

其中，n2为附加噪声，为噪声方差；

5.2)将特征分解为其中，p是pq×pq维的酉矩阵，为特征值构成的对角阵；

5.3)定义与λ2max所对应的特征向量为其中为x2的信号子空间，利用的第1至p×(q-1)个元素构成向量利用的第(p+1)至q×p个元素构成向量根据信号子空间的旋转不变性，得到：

其中，ac为a2的前(p-1)×q个元素构成的向量，ad为a2的后(p-1)×q个元素构成的向量，ac与ad满足ad＝acδ，δ＝e^iα，m为可逆矩阵，且满足

5.4)根据式<5>，得到与间的关系式：

其中，∑＝m^-1δm为转化矩阵，δ为∑的特征值；

5.5)根据上述推导，利用及得到转化矩阵∑：

其中(·)⁺表示矩阵的广义逆；

5.6)对转化矩阵∑进行特征分解得到特征值δ，提取特征值δ的辐角作为x轴和轴夹角α的估计值

步骤6，估计发送端方位角及俯仰角。

6.1)当接收端信噪比远大于1时，将接收信号的幅值近似表示为：

将式<6>中右边的各项移到式<6>的左边，整理得到零点函数

其中，的下标kp，lq表示接收信号的子载波为kp、模态为lq；

6.2)利用二分法估算的零点：

根据实际情况，仅当接收端位于发送端主瓣波束范围[θa，θb]内时，才能够接收到清晰的信号，因此，发送端俯仰角需满足关系θ∈[θa，θb]，此外，由于本方案所用场景为一对一通信，即当θ∈[θa，θb]时，零点函数有且仅有一个零点，且因此，可利用二分法估算零点函数的零点，将某一子载波kp及模态lq下零点函数的零点作为发送端俯仰角的估计值其具体步骤如下：

6.2a)求发送波束主瓣范围[θa，θb]的中点θc，给定估计精度ξ；

6.2b)计算θ＝θc时零点函数的函数值判断是否为0：若则以θc作为俯仰角的估计值否则，执行6.2c)；

6.2c)分别计算θ＝θa和θ＝θb时零点函数的函数值和

6.2d)判断是否满足若满足，则令θb＝θc，否则，令θa＝θc；

6.2e)判断发送端俯仰角取值范围[θa，θb]是否满足|θa-θb|＜ξ：若满足，则以θa作为发送端俯仰角的估计值否则，返回6.2c)；

6.3)对于不同波数kp及模态lq下的函数均执行6.2)得到对应的发送端俯仰角估计值

其中，θ为发送端俯仰角的真实值，为真实值θ与估算值间的误差p＝1，2，…，p，q＝1，2，…，q，由于各个的估算过程相同，则不同的具有相同的方差；

6.4)计算步骤6.3)中得到的全部的均值将的方差表示为：

其中，d(ε)为的方差，根据式<7>可知则本步骤以均值作为发送端俯仰角的估计值；

6.5)根据式<2>计算发送端方位角的估计值其中，为由步骤5)计算得到的x轴和轴夹角α的估计值。

本发明的效果可以通过以下仿真结果进一步说明：

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-desprit算法及二分法估算发送端的方位角及俯仰角，结果如图4所示。

从图4表明，随着信噪比的增加，利用本发明方法估计得到的方位角及俯仰角逐渐趋近于真实值，当信噪比为20db时，发送端方位角的估计值为7.000度，发送端俯仰角的估计值为6.993度，与发送端方位角与俯仰角的真实值非常接近。可见，利用本发明进行发送端方位角与俯仰角θ的估计能够很好的解决oam通信系统涡旋电磁波到达角估计的问题。

综上，本实例能够精确估计oam通信系统发送端的到达角，为oam通信系统的实际商用奠定了基础。