均匀圆阵模式空间下的宽带信号DOA估计.docVIP

均匀圆阵模式空间下的宽带信号DOA估计 徐青 电子科技大学电子工程学院，四川成都（610054） E-mail：alex357@ 摘 要：基于均匀圆阵（UCA）提出了一种宽带信号测向新方法。该方法首先将阵元空间数据变换为圆阵模式空间数据，使宽带信号有效带宽内各频点上的数据满足虚拟均匀线阵（VULA）模型，然后利用双边聚集变换（TCT）消除VULA阵列流形中的频率依赖性并进行能量积累，最终实现宽带信号测向。计算机仿真表明了新方法的有效性。 关键词：宽带信号；双边聚焦变换；到达角估计；均匀圆阵；模式空间 中图分类号：TP271+.5 1. 引言 宽带信号测向一直是阵列信号处理的主要课题之一。但是目前针对宽带信号的DOA算法绝大多数都是利用均匀线阵（ULA）进行估计，用均匀圆阵列天线（UCA）估计的并不多。而均匀圆阵（UCA）与均匀线阵（ULA）相比在DOA估计中具有很大的优越性：它不仅可以提供360度全方位、无模糊的方位角信息，还可以同时提供俯仰角信息，同时无论在任何方向上都具有近似相同的估计精度和分辨力。因此UCA在通信，雷达探测，导航制导等许多领域都得到了广泛的应用。 模式空间算法是针对均匀圆阵(UCA)所提出的一种空间谱估计算法。因为均匀圆阵阵列流型的表达方式比较复杂，且不具备均匀线阵的Vandermonde形式，所以经过模式空间变换之后，一些只适用于等距均匀线阵的算法可以适用于等距均匀圆阵。在宽带信号测向算法中，子空间类算法的双边相关变换（TCT）算法相比于其他相干信号的处理方法（CSM）[2决宽带信号的DOA估计问题。 在宽带信号TCT算法的聚焦过程中，由于均匀圆阵模式空间转换导致聚焦变换矩阵的维数随频率变化而变化，这是本算法的难点，为克服这个问题我们选取信号有效带宽内的最低频率作为聚集频率。 ～3]来说均方误差最小、稳定性高。因此本文在模式空间转换的基础上，结合经典的TCT算法来解 2. 宽带信号UCA模型 假设N个远场窄带信号源（θ，…，θ）入射到M个阵元的均匀圆阵上，r为圆阵半1N 径，圆阵阵元各向同性，则无噪声时，阵元l上的快拍数据为 N2lπs l=0,1,…,M-1 （1） x=∑exp(-jβcos(θ-liMii=1 其中β=2πr/λ，λ为入射信号的波长。该圆阵的阵列导向矢量为 a(θ)=[exp(?jβcos(θ))......exp(?jβcos(θ?2π(M?1))]T （2） M 如果入射信号为宽带信号源，则均匀圆阵的接收信号经过离散傅氏变换可分解为J个互不重叠的窄带部分： X(f)=A(f)S(f)+N(f) j=1,2,…,J （3） θjjjj -1- 式中，A(f)是均匀圆阵在频率为fj时的M×N方向矩阵，N×1源向量S(fj)为源信号fjθj 频率段的输出，M×1噪声向量N(fj)代表M个阵元上附加噪声的fj频段数据。 由于X(fj)不具备均匀线阵的Vandermonde矩阵形式，所以我们采用以模式空间为桥梁 将阵元空间中的均匀圆阵（UCA）换成模式空间的虚拟均匀线阵(ULA)的思路。根据文献[4]，对圆阵的输出X(fj)作如下变换： Y(fj)=T?X(fj) （4） 其中T为模式转换矩阵，得到的Y(fj)即为虚拟线阵（ULA）的输出，其模式空间中的数据 协方差矩阵为 R(f)=Y(f)?Y(f)=T?X(f)?X(f)?T=T?R(f)?T （5） 3. 基于UCA的宽带信号测向 3.1 模式空间下虚拟阵列模型 对于窄带信号而言，在均匀圆阵（UCA）接收到数据X后，首先要通过模式转换（4） 式将其变成虚拟线阵的输出形式Y，根据文献[4]可知，其模式转换矩阵为 T=1?1HJF （6） M 转换后Y(fj)的矩阵形式可表示为： S(f)Y(fj)=Ajj ?exp(?jKθ1)…exp(?jKθN)??s1(fj)? ? （7）???=? ????exp(jKθ) exp(jKθ)??s(f)?1N???Nj? 是模式转换后虚拟线阵（VULA）的方向矩阵，可以看出其具有Vandermonde结构，其中Aj (θ)=[exp(?jKθ)...exp(jKθ)]T。式（6）中矩阵J和F可分别由虚拟的导向矢量为 a 下两式表出: ?j?KJ?K(?β)0?J=?0% ?00? F=[ω?K?0?0 ? （8）jKJK(?β)??ω?K+1...ωK] ??2π?exp(jK) ? （9）M? ?2π?K(M?1))? exp(jM?111??2π2π?exp(j(?K))exp(j(?K+1))?MM=? ?2π2π?(?K)(M?1))exp(j(?K+1)(M?1))?exp(jMM? -2- 其中Jk(?)为k阶第一类bessel函数，β=2πr/λ，K=?