doa估计小论文解读.docxVIP

基于不同阵列模型的目标（DOA）估计 黎建春 西安电子科技大学电子信息攻防对抗与仿真技术教育部重点实验室陕西西安710071 摘要：最近三十年是阵列信号发展迅猛的时期，出现了大量优秀的阵列信号处理 算法，为在此基础上的目标到达角（DOA的估计提供了理论依据，本文主要针 对不同的阵列模型，提出不同的 DOA古计算法，并进行仿真分析。 关键字：阵列，到达角估计 1引言 2阵列模型及其算法仿真 2.1 直线型线性阵列 DOA估计 2.1.1原理介绍 均匀线阵（ULA : Uniform Lin ear Array ）是一最简单常用的阵列形式， 如图1所示，将M个 阵元等距离排列成一直线， 阵元间距为d。假定一信源位于远场， 即其信号到达各阵元的波 前为平面波，其波达方向（DOA ）定义为与阵列法线的夹角 321 3 2 1 图1 ULA示意图 以第一个阵元为参考阵元，则各阵元相对参考阵元的时延为: (5) 由此可得等距线阵的方向向量为: p /尊旳叩,/欝2dsine,L 訐矽(M/ysineF 1訂奢呷,严哪扎，严心丿呻）［（6） 当波长和阵列的几何结构确定时， 该方向向量只与空间角 二有关，因此等距线阵的方向 向量记为ar，它与基准点的位置无关。若有 D个信号源，其波达方向分别为 可， i =1,2丄,D，则阵列流形矩阵为： A = a -1 i，a R ,|l|,a %」 - 1 1 L 1 1 e』萨ine） 「2严 L __j3jld sir e扣 e） M M M M e e-^M~*dsin (日) L sin? (7) 以上给出了等距线阵的方向向量的表示形式。实际使用的阵列结构要求方向向量 a二 与空间角二一一对应，不能出现模糊现象。 这里需要说明的是： 阵元间距d是不能任意选定 的，甚至有时需要非常精确的校准。 假设d很大，相邻阵元的相位延迟就会超过 2二，此时， 阵列方向向量无法在数值上分辨出具体的相位延迟， 就会出现相位模糊。可见，对于等距线 阵来说，为了避免方向向量和空间叫之间的模糊，其阵元间距不能大于半波长 0，以保证 2 阵列流形矩阵的各个列向量线性独立。天线阵列的输出为： 』 * —I ^fm 4 d sin 0 y t s t wme 0 mF ( 8) 其向量形式为： (9)y k =wHX k (9) 式中，w -〔\ ”1,W2丄,Wm 1为权重向量。  对于上述的直线型均匀纯属阵列的 DOA估计，可采用 MUSIC算法进行估计 Music算法利用了信号子空间和噪声子空间的正交性，构造空间谱函数，通过谱峰搜索，检 测信号的DOA。它是建立在以下假设基础上的： （1） 阵列形式为线性均匀阵，阵元间距不大于处理最高频率信号波长的二分之一； （2） 处理器的噪声为加性高斯分布，不同阵元间距噪声均为平稳随机过程，且相互独立，空 间平稳（各阵元噪声方差相等）； （3） 空间信号为零均值平稳随机过程，它与阵元噪声相互独立，且信号间相互独立； （4） 信号源少于阵元数，信号取样数大于阵元数。 在此假设基础上， Music算法对DOA的估计从理论上可以有任意高的分辨率。 Music算法原理如下： 由式（4）可得接收信号的协方差矩阵为： R x 二E X t XH t TOC \o 1-5 \h \z =ae sS〕￡ + a「s卜e ~ Ns〕hae ~ nTN （22） 由于假设信号与噪声是不相关的，且噪声为平稳的加性高斯白噪声，因此式（ 22）中 的二，三项为零，且有 E || NN H wN 1。则式（22）简化为式（23）： H 2 RX = AR sA +ctn I （ 23） 式（23）中的Rs是有用信号的协方差矩阵。 「J由于假设信号源之间互不相关，因此 Rs为满秩矩阵，其秩为 D。而A为M D维的 矩阵，其秩也是D，并且AR sA H是Hermite半正定矩阵，其秩也是 D。因此，令AR sA H 的特征值为 ％ _叫_L _』? 0，那么R X的M 「J ^0 , L ,D - 1 k=D,D 1_, M,— 1 它们对应的特征向量分别为 q°,qi丄,qDj,qd,L ,qm 4，其中前D个对应大特征值, 后M - D个对应小特征值。 由此可以看出，协方差矩阵 Rx经过特征值分解后可以产生 D个较大的特征值和 M -D个较小的特征值，并且这M -D个小特征值非常接近。 所以当这些小特征值的重数 K确定了，那么信号的个数就可以由式（ 24）估计出来： (24)3 =M -K (24) 对于与M _D个最小特征值对应的特征向量，有: Rx -」qi = 0 , i = D,D 1 丄，M -1 即： 2 h 2 2 R x - N 1 q i = AR sA ■二 N 1 - ；N 1