《武汉工程大学学报》 2013年03期
75-78
出版日期：2013-03-31
ISSN:1674-2869
CN:42-1779/TQ
均匀圆阵下单个近场源信号四维参数估计快速算法
0引言利用空间位置分散的阵元组成的阵列来进行空间信号源定位在雷达、声纳和无线通信等方面有着重要的应用.而用均匀圆阵（Uniform Circle Array，以下简称：UCA）估计信号源位置相对于其它阵列结构有着360°的全方位角覆盖范围、几乎不变的方向图和额外的俯仰角信息等优点\[1-2\].近年来，文献\[1-4\]研究了基于圆阵的近场源无源定位问题.文献\[2\]中首先利用三维多重信号分类算法（3Dimensions Multiple Signal Classification Method，以下简称：3D MUSIC）来联合估计出近场源的方位角、俯仰角和距离\[2\]，但基于MUSIC的方法由于需要3D搜索计算因而计算复杂度高而难以工程实时实现.通过在UCA阵的中心增加一个阵元，Bea\[3\]等最近提出了一种计算量比较小的基于最小二乘（Least Square，LS）的近场源定位算法，但是这种方法只适用于非圆信号源.笔者提出了一种基于最小二乘的快速信源定位算法，该算法首先给出方位角、俯仰角和频率的估计，然后再估计出距离参数.和文献\[3\]不同的是，该方法无需利用UCA中心的参考阵元，且信号源无需非圆特性假设.该算法的核心思想是利用UCA的中心对称结构构造出两个相关序列，且该相关序列的相位上包含信号的位置信息，然后利用构造的相关序列的相位信息估计出信号源的4D 参数.值得注意的是所提算法可以认为是远场定位方法的一个推广.仿真结果表明提出的算法和单个信号源的4D MUSIC算法在性能上相当，但在计算复杂度上所提算法因不需要三维搜索和特征分解而具有显著优势.1提出的方法假定中心对称均匀圆阵的半径为R，且放置M（2的倍数）个阵元，所有的阵元分布在xy平面上，且UCA的中心为相位参考点.假设窄带近场信号源在球面坐标系统上的位置参数为(θ,,r)：θ∈\[0,2π)是信号的方位角，表示从x轴的正方向按逆时针方向所旋转的角度；∈\[0,π/2)是信号的俯仰角，表示从z轴的正方向按从上往下所旋转的角度；距离r是信号离UCA中心的距离，设w0是信号载波的频率，则第k个阵元上的接收数据可以写成\[2\]：xk(n)=s(n)ejw0nej2πλ(r-rk(θ,,r))+wk(n),k=1,2,3,…,M,n=1,2,3,…, N(1)式（1）中：s(n)为零均值的复窄带信号，其功率为E\[|s(n)|2\]=σ2s，E\[·\]为期望运算符；噪声{wk(n)}是各态历经、零均值、空域和时域上的复白高斯过程，且和信号源相互独立；N为快拍数；λ是信号的波长；rk(θ,,r)是信号源和第k个阵元间的距离，它可以写成如下的形式：rk(θ,,r)=r2+R2-2rRζk(θ,),k=1,2,3,…,M(2)这里，ζk(θ,)=cos(γk-θ)sin (3)γk=2π(k-1)/M为第k个阵元的方位角.我们的任务是在获得{xk(n)}后估计θ，，r,w0.当r足够大时，根据泰勒展开，式（2）可以近似为：rk(θ,,r)≈r-Rζk(θ,)+R22r(1-ζ2k(θ,))(4)第3期吴云韬，等：均匀圆阵下单个近场源信号四维参数估计快速算法武汉工程大学学报第35卷把式（4）带入到式（1）中，有xk(n)=s(n)ejw0nej2πRλ\[ζk(θ,)-R2r(1-ζ2k(θ,))\]+wk(n)(5)在阵列配置的中心对称性和阵元数为偶数的假设下，可以得到γM/2+k=γk+π，并由此可以得到：ak(θ,)=a*k+M/2(θ,), k=1,2,3,…, M/2(6)式（6）中：ak(θ,)为第k个阵元对信号源的复响应.(·)*为复共轭算子.同时有关系式γM/4+k=γk+π/2，由此可以得到：ξk+M/4(θ,,r)=-sin(γk-θ)sin ,k=1,2,3,…, 3M/4(7)根据式（6）和式（7），利用{xk(n)}构造两个相关序列且：p(k)=1N∑Nn=1xk(n)x*k+M/2(n),k=1,2,3,…,M/2(8)和q(k)=1N∑Nn=1xk(n)x*k+M/4(n),k=1,2,3,…,3M/4(9)当N和信噪比都足够大时，且当各态历经条件成立时，xk(n)≈xk(n+1)，p(k)可以近似表示为\[1\]，p(k)≈E{xk(n+1)x*k+M/2(n)}=σ2sak(θ,)a*k+M/2=ak(θ,)σ2sejw0ej4πRλcos(θ-γk)sin ,k=1,2,3,…,M/2(10)同样的，式（9）可以如下近似：q(k)≈E{xk(n)x*k+M/4(n)}=σ2sak(θ,)a*k+M/4=σ2sej2πRλ\[2cos(rk-θ-π/4)sin +R2rcos(2θ-2γk)sin2 \],k=1,2,3,…,3M/4(11)1.1信号源角度和载波频率的估计首先从式（10）中估计出θ、和w0，然后从式（11）中估计出r.令wk为p(k)的相角,从式（10）中可以得到ωk≈w0+4πRλcos(θ-γk)sin +2πmk,k=1,2,3,…, M/2(12)式（12）中mk为整数.假设式(12)中的相位模糊问题已经通过选择适当的mk得到解决\[5-6\]，这样就可以得到唯一的角度估计.不失一般性，假设mk=0，于是式（12）可以写成ωk≈4πRλ\[cos γk sin γk w0\]sin cos θsin sin θ1,k=1,2,3，…,M/2(13)写成矩阵形式：w≈Ab(14)这里w=\[ω1ω2…ωM/2\]Tb=4πRλ\[sin cos θ sin sin θ 1\]TA=cos γ1sin γ1w0cos γ2sin γ2w0………cos γM/2sin γM/2w0这里，(·)T为转置算子.于是b的最小二乘估计，写成b＾，为b＾=\[b＾1b＾2b＾3\]T=(ATA)-1ATw(15)对θ,,w0的估计，写成θ＾，＾，w＾0，可以由b＾给出：θ＾=arctan(b＾2/b＾1)(16)＾=arcsin(λb＾21+b＾22/4πR)(17)w＾0=b＾3(18)需要指出的是，上述估计步骤和文献\[4,67\]是一样的.1.2信号源距离的估计利用前面估计的角度结果，笔者将给出距离的估计，设uk为qk的相位，可以从式（11）中得到如下关系：uk≈2πRλ\[2cos(γk-θ＾-π/4)sin ＾+R2rcos(2θ＾-2γk)sin2 ＾\],k=1,2,3,…,3M/4(19)将式（19）表示成矩阵形式有，u≈d+R2rv(20)这里，u=\[u1u2…u3M/4\]v=sin2 ＾\[cos(2θ＾-2γ1)cos(2θ＾-2γ2)…cos(2θ＾-2γ3M/4)\](21)d=sin ＾\[2cos(θ＾-γ1-π/4)2cos(θ＾-γ2-π/4)…2cos(θ＾-γ3M/4-π/4)\]利用式（20）得到r的最小二乘估计r＾，即r＾=2Rv(u-d)T((u-d)(u-d)T)-1(22)2仿真结果笔者将所提算法和3D MUSIC\[2\]的性能在MATLAB环境下进行比较，入射信号为单个窄带近场信号，噪声为加性高斯白噪声.设定仿真参数M=8,R=λ=1，信号参数为(θ,,r,w0)=(70,10,1.5R,0.2)，快拍数N=200.用均方差来评价算法的估计性能，性能结果为300次独立运行的均值.图1~图4给出了所提算法proposed和传统的基于MUSIC方法的参数估计性能比较，图中均方误差越小越好.由图1~图4可以看出，两种方法的估计性能大致相当，但在计算复杂度方面，比较300次独立运算的总时间即单个信噪比（Signal to Noise Ratio,以下简称：SNR），所提算法只要0.224 6 s，而3D MUSIC方法却要349825 4 s, 表明本文算法具有显著优势. 图1θ的均方误差随着SNR的变化曲线Fig.1Mean square error of θ versus SNR图2的均方误差随着SNR的变化曲线Fig.2Mean square error of versus SNR图3r的均方误差随着SNR的变化曲线Fig.3Mean square error of r versus SNR图4w0的均方误差随着SNR的变化曲线Fig.4Mean square error of w0 versus SNR3结语本文针对UCA阵列结构，提出了一种对单个近场源信号进行定位的快速算法，并给出了信号的方位角、俯仰角、距离和载波频率参数的闭式解.仿真结果表明，所提算法在性能上和3D MUSIC相近，但具有更低的计算复杂度而利于工程实时应用.该方法只能处理单个近场信号源，而无法解决存在多个信号源情况下的定位问题，推广到多个信号源是笔者需要进一步探讨的问题.致谢本工作得到国家自然科学基金委员会和湖北省高等学校优秀中青年科技创新团队计划项目(T201206)提供的联合资助，武汉工程大学计算机学院、智能机器人湖北省重点实验室提供了有利的实验条件，在此致以衷心的感谢！