专利名称：二维到达角的获取方法和装置的制作方法
技术领域：
本发明涉及无线通讯领域，尤其涉及一种三维波束赋形(Three DimensionBeamforming, 3D BF)的到达角获取方法和装置。
背景技术：
来自空间的无线电波,其到达方向(Direction of arrival, D0A)可以用水平方位角Θ和垂直俯仰角φ表示，如图1所示。传统的处理方法认为无线电波来自远场，其俯仰角很小，如图2所示，可以忽略，从而只考虑水平方位角。这种不考虑俯仰角的方法，对移动站离通信站比较远时可以提供非常好的近似。但对于移动站离通信站比较近时，垂直方向的俯仰角一般比较大，不能忽略，如密集城区的微基站下的用户。我们称这种既考虑水平方位角Θ和垂直俯仰角Φ的到达方向为二维D0A,或者二维到达角(Angle of Arrival,A0A),如果无线电波是从天线阵发送出去的，也叫二维离开角(Angle of Departure, A0D)。本文统一称二维到达角。到达角估计广泛应用于波束赋形技术中。波束赋形是一种信号处理技术，它基于自适应天线原理，利用天线阵列通过先进的信号处理算法分别对各天线单元加权处理，使阵列实时对准有用信号方向，而在干扰方向形成零点以抑制干扰信号。从而提高信噪比，提升系统性能，增加系统的覆盖范围。传统的波束赋形或者到达角估计，都只进行水平方位角的估计。从而能区分水平方位角不同的移动站，如图3(a)所示。而对于水平方位角一样，垂直俯仰角不同的移动站是不能区分的，如图3(b)所示。这限制了波束赋形的自由度。而估计二维到达角，并将它应用到波束赋形中，就形成了既能自适应地调整水平方位角度，又能自适应地调整垂直俯仰角度的三维波束赋形。三维波束赋形在垂直覆盖(如图4所示)、垂直方向的干扰抑制(如图5所示)等方面有着二维波束赋形无法比拟的优势，它能更好地满足第四代无线通信系统对覆盖和干扰抑制的要求。二维到达角的另外一个应用是移动站的定位。在知道通信站的位置、高度以及二维到达角后，也可以很容易地确定移动站的位置，而移动站的定位，是非常有意义和价值的。总之，二维到达角的估计在三维波束赋形，移动站定位等领域有着广泛的应用，是一种非常值得研究的技术。
发明内容
本发明要解决的技术问题是提供一种二维到达角的获取方法和装置，可简化二维到达角获取的计算复杂度。为了解决上述问题，本发明提供了一种二维到达角的获取方法，包括:获取移动站和通信站之间的第一整体信道系数矩阵，根据所述第一整体信道系数矩阵提取垂直方向的信道系数矩阵，根据所述垂直方向的信道系数矩阵计算俯仰角；以及，获取移动站和通信站之间的第二整体信道系数矩阵，根据所述第二整体信道系数矩阵提取水平方向的信道系数矩阵，根据所述水平方向的信道系数矩阵计算方位角；其中，所述第一整体信道系数矩阵和第二整体信道系数矩阵为所述通信站和移动站之间相同时刻或不同时刻的整体信道系数矩阵。进一步的，上述方法还可具有以下特点，所述第一整体信道系数矩阵或第二整体信道系数矩阵是用来计算二维到达角的资源块上的信道系数矩阵H(I)，所述H(I)中的元
表示所述通信站的第m根天线到所述移动站的第n根天线在所述资源块上第I个载波上的信道系数，n = 1，…，Nk，当所述通信站的天线阵列索引以水平方向为主序时，m = (i_l)NH+j,当所述通信站的天线阵列索引以垂直方向为主序时，m = (j_l)Nv+i, i =I，…，Nv, j = I,…，Nh, Nk表示所述移动站的天线数目，Nh表示所述通信站的天线阵列在水平方向上的天线数目，Nv表示所述通信站的天线阵列在垂直方向的天线数目；1 = I,-,Ns表示所述资源块上的载波个数。进一步的，上述方法还可具有以下特点，所述根据第一整体信道系数矩阵提取垂直方向的信道系数矩阵包括:当所述通信站的天线阵列索引以水平方向为主序时，所述垂直方向的信道系数矩阵Hv(l, j)为所述H(I)中所有列索引满足Hii = (i_l)NH+j,i = I,...,Nv的列构成的矩阵；或者，当所述通信站的天线阵列索引以垂直方向为主序时，所述Hv(l，j)为H(I)中所有列索引满足Ini = (j-l)Nv+i, i = I,…，Nv的列构成的矩阵；其中，j= I，…，Nh, 1 = 1, NSo进一步的，上述方法还可具有以下特点，所述根据垂直方向的信道系数矩阵计算俯仰角包括:`根据所述垂直方向的信道系数矩阵计算垂直方向的相关矩阵Rv，计算垂直方向的激励向量ev与所述垂直方向的相关矩阵Rv的二次型的绝对值化)，将使得
取最大值时的角度0 v,_作为垂直方向的角频率，根据所述垂直方向的角频率
、■得到最终的俯仰角，所述 =[I ^…进一步的，上述方法还可具有以下特点，所述根据第二整体信道系数矩阵提取水平方向的信道系数矩阵包括:当所述通信站的天线阵列索引以水平方向为主序时，所述水平方向的信道系数矩阵^^(1,；0为所述H(I)中所有列索引满足nij = (1-l)NH+j, j = I,…,乂的列构成的矩阵；或者，当所述通信站的天线阵列索引以垂直方向为主序时，所述水平方向的信道系数矩阵4(1，i)为所述H(I)中所有列索引满足nij = (j-l)Nv+i, j = I,…，Nh的列构成的矩阵；其中，i= 1，…，Nv, 1 = 1,...，NS。进一步的，上述方法还可具有以下特点，所述根据水平方向的信道系数矩阵计算方位角包括:根据所述水平方向的信道系数矩阵计算水平方向的相关矩阵Rh，计算水平方向的激励向量eH与所述水平方向的相关矩阵Rh的二次型的绝对值，将使得取最大值时的角度θ Η,_作为水平方向的角频率，根据所述水平方向的角频率9H,max求得最终的方位角，所述水平方向的激励向量4=[1 一...e][N^ J 0
进一步的，上述方法还可具有以下特点，所述方法还包括:以第一周期周期性获取所述第一整体信道系数矩阵并计算俯仰角，以第二周期周期性获取所述第二整体信道系数矩阵并计算方位角，且所述第一周期大于第二周期。本发明还提供一种二维到达角的获取装置，包括:整体信道系数矩阵获取模块，用于获取移动站和通信站之间的第一整体信道系数矩阵和第二整体信道系数矩阵，其中，所述第一整体信道系数矩阵和第二整体信道系数矩阵为所述通信站和移动站之间相同时刻或不同时刻的整体信道系数矩阵；垂直方向信道系数矩阵提取模块，用于根据所述第一整体信道系数矩阵提取垂直方向的信道系数矩阵；俯仰角计算模块，用于根据所述垂直方向的信道系数矩阵计算俯仰角；水平方向信道系数矩阵提取模块，用于根据所述第二整体信道系数矩阵提取水平方向的信道系数矩阵；方位角计算模块，用于根据所述水平方向的信道系数矩阵计算方位角。进一步的，上述装置还可具有以下特点，所述整体信道系数矩阵获取模块获取的所述第一整体信道系数矩阵或第二整体信道系数矩阵是用来计算二维到达角的资源块上的信道系数矩阵H(I)，所述H(1)中的元素hn,m(l)表示所述通信站的第m根天线到所述移动站的第η根天线在所述资源块上第I个载波上的信道系数，η = 1，-,Ne,当所述通信站的天线阵列索引以水平方向为主序时，m= (i_l)NH+j，当所述通信站的天线阵列索引以垂直方向为主序时，m = (j-l)Nv+i, i = I,..., Nv, j = I,..., Nh, Nk表示所述移动站的天线数目，Nh表示所述通信站的天线阵列在水平方向上的天线数目，Nv表示所述通信站的天线阵列在垂直方向的天线数目；1 = 1，…，Ns表示所述资源块上的载波个数。进一步的，上述装置还可具有以下特点，所述垂直方向信道系数矩阵提取模块根据第一整体信道系数矩阵提取垂直方向的信道系数矩阵包括:当所述通信站的天线阵列索引以水平方向为主序时，所述垂直方向的信道系数矩阵Hv(l, j)为所述H(I)中所有列索引满足Hii = (i_l)NH+j,i = I,...,Nv的列构成的矩阵；或者，当所述通信站的天线阵列索引以垂直方向为主序时，所述Hv(l，j)为H(I)中所有列索引满足Ini = (j-l)Nv+i, i = l,..., Nv的列构成的矩阵；其中，j= I，…，Nh, 1 = 1,..., NSo进一步的，上述装置还可具有以下特点，所述俯仰角计算模块根据垂直方向的信道系数矩阵计算俯仰角包括:根据所述垂直方向的信道系数矩阵计算垂直方向的相关矩阵Rv，计算垂直方向的激励向量ev与所述垂直方向的相关矩阵Rv的二次型的绝对值，将使得
取最大值时的角度θ ν,_作为垂直方向的角频率，根据所述垂直方向的角频率
Θ v,max得到最终的俯仰角，所述 =[I ^…进一步的，上述装置还可具有以下特点，所述水平方向信道系数矩阵提取模块根据第二整体信道系数矩阵提取水平方向的信道系数矩阵包括:当所述通信站的天线阵列索引以水平方向为主序时，所述水平方向的信道系数矩阵^^(1,；0为所述H(I)中所有列索引满足nij = (1-l)NH+j, j = I,…,乂的列构成的矩阵；或者，当所述通信站的天线阵列索引以垂直方向为主序时，所述水平方向的信道系数矩阵4(1，i)为所述H(I)中所有列索引满足nij = (j-l)Nv+i, j = I,…，Nh的列构成的矩阵；其中，i= 1，…，Nv, 1 = 1, NSo进一步的，上述装置还可具有以下特点，所述方位角计算模块根据水平方向的信道系数矩阵计算方位角包括:根据所述水平方向的信道系数矩阵计算水平方向的相关矩阵Rh，计算水平方向的激励向量eH与所述水平方向的相关矩阵Rh的二次型的绝对值，将使得
取最大值时的角度0 H，_作为水平方向的角频率，根据所述水平方向的角频率
0H,max求得最终的方位角，所述水平方向的激励向量一 - e][N^ J 0进一步的，上述装置还可具有以下特点，所述整体信道系数矩阵获取模块以第一周期周期性计算所述第一整体信道系数矩阵，以第二周期周期性计算所述第二整体信道系数矩阵，且所述第一周期大于第二周期。
采用本发明所述的方法和装置，可以简化二维到达角的获取。
图1是无线电波到达天线阵列的示意图；图2是无线电波到达天线阵列俯仰角很小的场景；图3是传统的波束赋形示意图；图4是用二维波束赋形实现垂直方向覆盖不意图；图5是用三维波束赋形实现垂直方向干扰抑制示意图；图6是本发明所述计算二维到达角方法的流程图；图7是本发明所述的二维到达角获取装置一实施方式；图8是本发明所述的天线阵列索引以水平方向为主序的天线阵列；图9是本发明所述的天线阵列索引以垂直方向为主序的天线阵列；图10是本发明所述的二维到达角获取装置另一实施方式。
具体实施例方式为使本发明的目的、技术方案和优点更加清楚明白，下文中将结合附图对本发明的实施例进行详细说明。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互任意组合。本发明用到的一些术语定义如下。通信站:宏基站、微基站、直放站、中继、拉远设备、无线接入点等各种无线通信设备。
移动站:数据卡、手机、笔记本电脑、个人电脑、平板电脑、个人数字助理、蓝牙等各种终端。本发明实施例提供了一种二维到达角的获取方法，应用于包括通信站和移动站的无线通信系统，如图6所示，包括:获取通信站和移动站之间的第一整体信道系数矩阵，根据所述第一整体信道系数矩阵提取垂直方向的信道系数矩阵，根据垂直方向的信道系数矩阵计算俯仰角；获取通信站和移动站之间的第二整体信道系数矩阵，根据第二整体信道系数矩阵提取水平方向的信道系数矩阵，根据水平方向的信道系数矩阵计算方位角；所述第一整体信道系数矩阵和第二整体信道系数矩阵为所述通信站和移动站之间相同时刻或不同时刻的整体信道系数矩阵。如果是相同时刻的，则表示基于同一整体信道系数矩阵计算俯仰角和方位角。如果不同，则表示基于不同的整体信道系数矩阵计算俯仰角和方位角。
·
进一步的，可以周期性计算俯仰角和方位角。一种实施方式中，计算俯仰角的周期为第一周期，计算方位角的周期为第二周期，且所述第一周期大于第二周期。统计发现，垂直方向的角度改变比水平方向慢，可以不用过快的更新垂直方向的俯仰角，设置第一周期大于第二周期可以减少计算量。当然，在本发明其它实施方式中，也可以第一周期和第二周期一样大，或者，第一周期小于第二周期。进一步地，所述第一整体信道系数矩阵是指用来计算二维到达角的资源块上的信道系数矩阵H (I)。其中，H(I)中的第η行第m列元素hn,m(l)表示通信站第m根天线到移动站第η根天线在所述资源块上第I个载波上的信道系数，η = 1，-,Ne,当所述通信站的天线阵列索引以水平方向为主序时，m= (i_l)NH+j，当所述通信站的天线阵列索引以垂直方向为主序时，m = (j-l)Nv+i, i = l,…，Nv, j = I,..., Nh, Nk表示移动站的天线数目，Nh表示通信站天线阵列在水平方向上的天线数目，Nv表示通信站天线阵列在垂直方向的天线数目；1 =1，…，Ns表示所述资源块上的载波个数，它可以仅是导频载波或者数据载波，也可以是导频载波和数据载波。天线阵列在水平方向和垂直方向的天线数目在基站摆放天线时就能够确定。其中，所述资源块是一个时频二维结构，在时域上占据一个或多个连续的符号，在频域上占据一个或多个连续的子载波。所述符号包括正交频分复用(OrthogonalFrequency Division Multiplexing, OFDM)或者正交频分多址(Orthogonal FrequencyDivision Multiple Access, 0FDMA)。所述数据载波或者导频载波是指资源块上I个OFDM (或者0FDMA)上的I个子载波。进一步地，根据第一整体信道系数矩阵提取垂直方向的信道系数矩阵包括:提取垂直方向的信道系数矩阵Hv(l，j)，当所述通信站的天线阵列索引以水平方向为主序时，Hv(l, j)为H(I)中所有列索引满足Hii = (1-l)NH+j, i = l,..., Nv的列构成的矩阵；当以垂直方向为主序时，Hv(l, j)为H(I)中所有列索引满足Hii = (j-l)Nv+i, i =1，…，Nv的列构成的矩阵。其中，j = l,…，Nh, 1 = 1,…，队。进一步地，由垂直方向的信道系数矩阵计算俯仰角包括:根据垂直方向的信道系数矩阵计算垂直方向的相关矩阵Rv，计算垂直方向的激励向量ev与垂直方向的相关矩阵Rv的二次型的绝对值，将使得取最大值时的角度9 v,_作为垂直方向的角频率，根据垂直方向的角频率9V,_求得最终的俯仰角。其中，垂直方向的相关矩阵为
权利要求
1.一种二维到达角的获取方法,包括: 获取移动站和通信站之间的第一整体信道系数矩阵，根据所述第一整体信道系数矩阵提取垂直方向的信道系数矩阵，根据所述垂直方向的信道系数矩阵计算俯仰角；以及， 获取移动站和通信站之间的第二整体信道系数矩阵，根据所述第二整体信道系数矩阵提取水平方向的信道系数矩阵，根据所述水平方向的信道系数矩阵计算方位角； 其中，所述第一整 体信道系数矩阵和第二整体信道系数矩阵为所述通信站和移动站之间相同时刻或不同时刻的整体信道系数矩阵。
2.如权利要求1所述的方法，其特征在于，所述第一整体信道系数矩阵或第二整体信道系数矩阵是用来计算二维到达角的资源块上的信道系数矩阵H(I)，所述H(I)中的元素hn,ffl(l)表示所述通信站的第m根天线到所述移动站的第n根天线在所述资源块上第I个载波上的信道系数，n = 1，…，Nk，当所述通信站的天线阵列索引以水平方向为主序时，m =(1-l)NH+j,当所述通信站的天线阵列索引以垂直方向为主序时,m = (j-l)Nv+i, i = I,…，Nv, j = I,…，Nh, Ne表示所述移动站的天线数目，Nh表示所述通信站的天线阵列在水平方向上的天线数目，Nv表示所述通信站的天线阵列在垂直方向的天线数目；1 = 1，…，队表示所述资源块上的载波个数。
3.如权利要求2所述的方法，其特征在于，所述根据第一整体信道系数矩阵提取垂直方向的信道系数矩阵包括: 当所述通信站的天线阵列索引以水平方向为主序时，所述垂直方向的信道系数矩阵Hv(l, j)为所述H(I)中所有列索引满足Hii = (1-1)Nh+j, i = I,…，Nv的列构成的矩阵；或者， 当所述通信站的天线阵列索引以垂直方向为主序时，所述Hv(l，j)为H(I)中所有列索引满足Hii = (j-l)Nv+i, i = I,…，Nv的列构成的矩阵； 其中，j = l，…，Nh, 1 = 1,…，队。
4.如权利要求3所述的方法，其特征在于，所述根据垂直方向的信道系数矩阵计算俯仰角包括: 根据所述垂直方向的信道系数矩阵计算垂直方向的相关矩阵Rv，计算垂直方向的激励向量ev与所述垂直方向的相关矩阵Rv的二次型的绝对值化)，将使得取最大值时的角度0V,_作为垂直方向的角频率，根据所述垂直方向的角频率得到最终的俯仰角，所述 =[I ^…
5.如权利要求2所述的方法，其特征在于，所述根据第二整体信道系数矩阵提取水平方向的信道系数矩阵包括: 当所述通信站的天线阵列索引以水平方向为主序时，所述水平方向的信道系数矩阵Hh(I, i)为所述H(I)中所有列索引满足nij = (1-1)Nh+j, j = I,，Nh的列构成的矩阵；或者， 当所述通信站的天线阵列索引以垂直方向为主序时，所述水平方向的信道系数矩阵Hh(I, i)为所述H(I)中所有列索引满足nij = (j-l)Nv+i, j = I,…，Nh的列构成的矩阵； 其中，i = I,…，Nv, 1 = 1,…，Ns。
6.如权利要求5所述的方法，其特征在于，所述根据水平方向的信道系数矩阵计算方位角包括: 根据所述水平方向的信道系数矩阵计算水平方向的相关矩阵rh，计算水平方向的激励向量eH与所述水平方向的相关矩阵Rh的二次型的绝对值，将使得取最大值时的角度ΘΗ,_作为水平方向的角频率，根据所述水平方向的角频率ΘΗ,_求得最终的方位角，所述水平方向的激励向量eH=[l e他…了。
7.如权利要求1至6任一所述的方法，其特征在于，所述方法还包括:以第一周期周期性获取所述第一整体信道系数矩阵并计算俯仰角，以第二周期周期性获取所述第二整体信道系数矩阵并计算方位角，且所述第一周期大于第二周期。
8.—种二维到达角的获取装置，包括: 整体信道系数矩阵获取模块，用于获取移动站和通信站之间的第一整体信道系数矩阵和第二整体信道系数矩阵，其中，所述第一整体信道系数矩阵和第二整体信道系数矩阵为所述通信站和移动站之间相同时刻或不同时刻的整体信道系数矩阵； 垂直方向信道系数矩阵提取模块，用于根据所述第一整体信道系数矩阵提取垂直方向的信道系数矩阵； 俯仰角计算模块，用于根据所述垂直方向的信道系数矩阵计算俯仰角； 水平方向信道系数矩阵提取模块，用于根据所述第二整体信道系数矩阵提取水平方向的信道系数矩阵； 方位角计算模块，用于根据所述水平方向的信道系数矩阵计算方位角。
9.如权利要求8所述的装置，其特征在于，所述整体信道系数矩阵获取模块获取的所述第一整体信道系数矩阵 或第二整体信道系数矩阵是用来计算二维到达角的资源块上的信道系数矩阵H(I)，所述H(I)中的元素hn,m(l)表示所述通信站的第m根天线到所述移动站的第η根天线在所述资源块上第I个载波上的信道系数，η = 1，-,Ne,当所述通信站的天线阵列索引以水平方向为主序时，m= (i_l)NH+j，当所述通信站的天线阵列索引以垂直方向为主序时，m = (j-l)Nv+i, i = I,..., Nv, j = I,..., Nh, Nk表示所述移动站的天线数目，Nh表示所述通信站的天线阵列在水平方向上的天线数目，Nv表示所述通信站的天线阵列在垂直方向的天线数目；1 = 1，…，Ns表示所述资源块上的载波个数。
10.如权利要求9所述的装置，其特征在于，所述垂直方向信道系数矩阵提取模块根据第一整体信道系数矩阵提取垂直方向的信道系数矩阵包括: 当所述通信站的天线阵列索引以水平方向为主序时，所述垂直方向的信道系数矩阵Hv(l, j)为所述H(I)中所有列索引满足Hii = (1-l)NH+j, i = l,..., Nv的列构成的矩阵；或者， 当所述通信站的天线阵列索引以垂直方向为主序时，所述Hv(l，j)为H(I)中所有列索引满足Hii = (j-l)Nv+i, i = l,..., Nv的列构成的矩阵； 其中，j = l，…，Nh, 1 = 1,…，队。
11.如权利要求10所述的装置，其特征在于，所述俯仰角计算模块根据垂直方向的信道系数矩阵计算俯仰角包括: 根据所述垂直方向的信道系数矩阵计算垂直方向的相关矩阵Rv，计算垂直方向的激励向量ev与所述垂直方向的相关矩阵Rv的二次型的绝对值，将使得取最大值时的角度0V,_作为垂直方向的角频率，根据所述垂直方向的角频率得到最终的俯仰角，所述4= [I W…
12.如权利要求9所述的装置，其特征在于，所述水平方向信道系数矩阵提取模块根据第二整体信道系数矩阵提取水平方向的信道系数矩阵包括: 当所述通信站的天线阵列索引以水平方向为主序时，所述水平方向的信道系数矩阵Hh(I, i)为所述H(I)中所有列索引满足nij = (1-1)Nh+j, j = I,，Nh的列构成的矩阵；或者， 当所述通信站的天线阵列索引以垂直方向为主序时，所述水平方向的信道系数矩阵Hh(I, i)为所述H(I)中所有列索引满足nij = (j-l)Nv+i, j = I,…，Nh的列构成的矩阵；其中，i = I,…，Nv, 1 = 1,…，Ns。
13.如权利要求12所述的装置，其特征在于，所述方位角计算模块根据水平方向的信道系数矩阵计算方位角包括: 根据所述水平方向的信道系数矩阵计算水平方向的相关矩阵Rh，计算水平方向的激励向量eH与所述水平方向的相关矩阵Rh的二次型的绝对值，将使得取最大值时的角度0H,_作为水平方向的角频率，根据所述水平方向的角频率求得最终的方位角，所述水平方向的激励向量eH=[l e他…了。
14.如权利要求7至13任一所述的装置，其特征在于，所述整体信道系数矩阵获取模块以第一周期周期性计算所述第一整体信道系数矩阵，以第二周期周期性计算所述第二整体信道系数矩阵，且所 述第一周期大于第二周期。
全文摘要
本发明提供一种二维到达角的获取方法，包括获取移动站和通信站之间的第一整体信道系数矩阵，根据所述第一整体信道系数矩阵提取垂直方向的信道系数矩阵，根据所述垂直方向的信道系数矩阵计算俯仰角；以及，获取移动站和通信站之间的第二整体信道系数矩阵，根据所述第二整体信道系数矩阵提取水平方向的信道系数矩阵，根据所述水平方向的信道系数矩阵计算方位角；第一整体信道系数矩阵和第二整体信道系数矩阵为通信站和移动站之间相同时刻或不同时刻的整体信道系数矩阵。本发明还提供一种二维到达角获取装置。本发明简化了二维到达角的获取。
文档编号H04B7/08GK103152093SQ201210037119
公开日2013年6月12日 申请日期2012年2月17日 优先权日2011年12月6日
发明者肖华华, 朱登魁, 宁迪浩, 鲁照华 申请人:中兴通讯股份有限公司