均匀圆阵下的DOA估计算法研究.pdf

1、均匀圆阵下的 DOA 估计算法研究*杨华 邱天爽 （大连理工大学 电子与信息工程学院，辽宁 大连 116023）E-mail： 摘要：波达方向估计(DOA)在移动通信中一直具有举足轻重的作用。本文基于第三代移动通信系统 TD-SCDMA 采用的 8 阵元均匀圆阵智能天线，首先给出一种利用内插阵列变换的圆阵平移一维 DOA 估计方法，该算法精度较高，并且能够分辨相干信号；接着介绍了UCA-RB-MUSIC 算法，该算法基于空间模式变换对信号源可进行一维和二维 DOA 估计。通过仿真实验说明内插阵列变换圆阵平移一维 DOA 估计方法对相干信号源可以有效估计，而UCA-RB-MUSIC 算法则不行；

2、但是 UCA-RB-MUSIC 算法可以实现对信号源的二维 DOA 估计。关键词：关键词：均匀圆阵；波达方向估计；MUSIC 中图分类号：中图分类号：TN91 Researches on DOA Estimation Algorithms for Uniform Circular Arrays Yang Hua Qiu Tianshuang (School of Electronic and Information Engineering Dalian University of Technology,Liaoning Dalian 116023)E-mail： Abstract:Direct

3、ion of arrival estimation always plays an important part in mobile communications.On the basis of uniform circular arrays(UCA)with 8 sensors adopted by the third-generation telecommunication system TD-SCDMA,this paper introduces a 1-D algorithm which takes advantage of interpolated arrays to realize

4、 the shifts of the UCA and get DOA estimations.This algorithm can tell two coherent signals.Whats more,the paper also introduces a 2-D DOA estimation algorithm for UCA,which is based on mode space transformation to get 1-D and 2-D DOA estimation.Simulations show that the algorithm based on interpola

5、ted arrays works better than UCA-RB-MUSIC algorithm when the signals are coherent in 1-D DOA estimation;while UCA-RB-MUSIC algorithm can realize effective 2-D DOA estimation instead of the algorithm based on interpolated arrays.Key Words:Uniform Circular Arrays（UCA）;DOA;MUSIC 1.引言1.引言 TD-SCDMA1是由我国自

6、主研发的移动通信系统标准，因此研究基于TD-SCDMA标准下的无线定位技术具有广泛的意义和价值。而在TD-SCDMA无线定位的智能天线圆-角定位技术中，DOA算法为其核心。TD-SCDMA采用的是8阵元均匀圆阵智能天线，与均匀线阵相比，均匀圆阵具有可同时估计方位角和俯仰角，具有360的方位角观察范围等优点。但是均 *本课题得到国家自然科学基金（60372081），高等学校博士学科点专项科研基金（20050141025）资助 http:/ 1匀圆阵的阵列流形为形式复杂的指数形式，传统的适用于均匀线阵的空间谱估计方法如MUSIC2、ESPIRIT3等不能直接应用到均匀圆阵中来。实现均匀圆阵DOA估

7、计比较常用的方法有基于内插阵列4和基于空间模式变换的DOA算法。本文首先给出了利用内插阵列变换的圆阵平移一维DOA估计算法，仿真试验表明该算法精度高并且能够分辨相干信号。接着介绍了基于空间模式变换的二维DOA估计算法：UCA-RB-MUSIC5算法，该算法可以有效的实现二维DOA估计。2.2.UCA的信号模型 的信号模型 如图1所示，信号入射到8N=阵元的均匀圆阵中，方位角和俯仰角分别为和。假设有D个非相干窄带信源入射到均匀圆阵中，分别记为()()tstsD,1?。它们的方位角和俯仰角分别记为()()DD,11?。其中，()2,0，2,0。那么圆阵的输出可以记为()()()()tntstxkk

8、kDk+=,1a （1）xyz14786532iin 图1 均匀圆阵示意图 Fig.1 Sketch map for UCA 其中，导向矢量()()()TjjNee,coscos0=?a，sin2 r=，Nnn21=+，1,0=Nn?。写成向量形式可得()()()tttnAsx+=（2）()()TNtxtx,1?=x，(),),(11DDaaA?=()()TNtntn,1?=n，()()TDtsts,1?=s （3）式中，噪声n是均值为零，方差为2的加性高斯白噪声。3.均匀圆阵下的3.均匀圆阵下的DOA估计算法 估计算法 3.1 基于内插阵列的圆阵平移3.1 基于内插阵列的圆阵平移DOA估计算

9、法 估计算法 内插阵列变换的基本思想6是将整个天线扫描区域划分为若干个子区域，再将某个子区域细分，假设信号位于区域，将区域划分为 111 2 rr=+?（4）1和r为的左右边界，为步长。则实际阵列的阵列流形矩阵为 http:/ 2111()()(2)()raaaa=+A?（5）而在同一区域内，内插阵列的阵列流形矩阵为 )()2()()(111raaaa?+=A （6）则内插阵列与实际阵列之间存在着一个固定的变换关系B，使得 1=AABABA （7）该算法的核心思想是构造多个内插阵列，使这些阵列为要进行 DOA 估计均匀圆阵的平移形式，从而可以实现对均匀圆阵的 DOA 估计。假设实际阵列的数据协

10、方差阵为xR，噪声协方差阵为nR，则内插阵列的数据协方差阵为()21121 KKHHHkkkkkkKHHkkkkk=+=+ssRB RBBAR AI BA R AB B （8）其中，kA是第k个内插阵列的导向矢量。那么我们可以得到这样的关系：kk=AAD （9）其中，()()11sincossincosdiag,kkkNkNjxyjxykee+=D?，kx，ky分别表示第k个内插阵列在x和y轴方向上的平移。将式（9）代入式（8）可得 2112 KKHHHkkkkkkHn=+=+ssRAD R DAB BAR AR （10）接下来再进行一下白化处理 1 22Hnn=RRRR?（11）然后对R?进

11、行特征值分解 H=RUU?（12）用nE?表示R?中包含了MN个小特征值的矩阵，那么根据MUSIC算法只需要对下式进行谱峰搜索就可以得到DOA估计值，()()1 21Hnn=SaRE?（13）在实际情况中我们利用样本协方差R?代替R，所以上式变为()()1 21Hnn=SaRE?。由于式（10）实际上相当于对阵列进行了空间平滑处理，因此该方法可以有效地分辨两个相干信号。在实际应用环境中，上述内插阵列变换方法面临的首要问题就是区域的划分。这里可先计算)(BAA的Frobenius范数和A的模的比值。如果这个比值足够小，就可以接受这个变换矩阵B，否则可对观察区域进一步细分,重新计算这个变换矩阵。虽

12、然为了得到精度较高的变换矩阵所需的运算量较大，但对于实际系统来说，这往往是一个离线计算的过程（即可以将观察区域划分为不同的观察区间TI，针对不同的观察区间计算出对应的变换矩阵B，并存储到系统中）。http:/ 3但是该算法只能进行一维DOA估计，在近场信号的情况下我们需要同时知道信号的方位角和俯仰角，因此下面给出一种二维DOA估计算法。3.2 3.2 UCA-RB-MUSIC算法 算法 UCA-RB-MUSIC算法由Cherian P.Mathews于1994年提出5，是著名的MUSIC算法的波束空间形式。该算法使用波束形成器HrF将UCA的导向矢量(),a映射到波束空间阵列流形(),b。()

13、,b在整个到达角的范围内都是实值的。首先引入基于相位模式激励5理论的波束形成器HeF，()(),eeaaF=H，(),ea是中心厄米特矩阵。将HeF之前再乘以一个具有中心厄米特列的矩阵就可以得到实值波束空间阵列流形，因此定义HHHerFWF=。HW具有中心厄米特列。HHVCFve=（14）其中,MMjjjjj=,diag101vC，MMNwwwV,0=（15）110,1=NjmjmjmHmeeeNw （16）式（5）中M为激励的最大模式，当智能天线为8阵元均匀圆阵时，3=M。这样，()()()()vJVCaFaNHH=,vee （17）()TjMjjjjMeeeee,0=v （18）()()(

14、)()()MMJJJJJ,diag101=J （19）设Q为置换矩阵，除副对角线元素为1外，其它元素均为0。这样，(),ea就满足()(),ee=aQa。同时易知IFF=eHe。令()MMiMii,122+=，则()()()MMMvvvW,10=（20）由 以 上 可 得HrF。若AFBrH=，()()()()ttttHHnFBsxFyrr+=，则 其 协 方 差 矩 阵 为()()IBPByyRy+=THttE。那么 IBBPRRy+=TRRe （21）其中，PPRe=R。对R进行实值特征值分解可以得到相互正交的信号子空间S和噪声子空间G。1211,+=MDDggGssS （22）这样，UC

15、A-RB-MUSIC谱就为()()(),1,bGGbTTS=（23）http:/ 4之后进行谱峰搜索可得DOA估计值。该算法的优点是可以既可以进行二维DOA估计，又可以进行一维DOA估计，缺点是不能处理相干信号。4.仿真实验 4.仿真实验 本实验为了验证了算法的有效性，分别将两种算法进行了仿真。实验中使用TD-SCDMA所采用的8阵元均匀圆阵智能天线，为保证相邻阵元间距为2，避免模糊产生，圆阵的半径设为()Nrsin4=，信号为窄带QPSK信号，快拍数p为1000；噪声为均值为0，方差为1的高斯白噪声；独立蒙特卡洛仿真100次。信噪比定义为2210log10SNRns=。4.1 仿真实验1 4

16、.1 仿真实验1 本实验进行一维DOA估计。图2（a）和图2（b）分别是内插阵列圆阵平移算法和UCA-RB-MUSIC算法对两个入射到均匀圆阵中的独立窄带信号源的一维DOA估计。dB10SNR=，两路信号源的入射方位角为（10，30）。从图中可以看出，两种方法都能较好的分辨两窄带独立信号，峰值尖锐，估计较为准确。下面我们考察当信号源相干的情况下两种算法的性能。其它条件同图二中的实验，将独立信号源换成相干信号源。图中实线为内插阵列圆阵平移算法的估计结果，虚线为UCA-RB-MUSIC算法估计结果，从图中我们可以得出结论，对于相干信号源内插阵列圆阵平移算法能够良好的分辨而UCA-RB-MUSIC算法已经失效。010203040506005101520253035MUSIC功 率 谱方 位 角（0-60）一维空间谱函数P/dB 0102030405060-20-100102030405060MUSIC功率谱方位角（0-60）一维空间谱函数P/dB（a）（b）图2 两独立窄带信号源，入射方位角分别为（10，30），dB10SNR=，峰值为估计值。（a）内插阵列圆阵平移算法估计结果（b）UCA-