均匀圆阵ROOT

第四章均匀圆阵ROOT-MUSIC算法ROOT-MUSIC算法是一种MUSIC算法的变形，它的性能明显优于谱峰搜索MUSIC算法，一般MUSIC算法需要谱峰搜索，而ROOT-MUSIC算法是基于多项式求根，它的运算量小，并且能够改善分辨门限和估计性能。对均匀线阵列(ULA),它的阵列流行阵列是范德蒙结构的，所以MUSIC算法的谱峰搜索过程可以用多项式求根代替，即ULA-ROOT-MUSIC算法；而对于均匀圆阵列(UCA),阵列流形矩阵不再是范德蒙结构，所以利用范德蒙结构开发的ROOT-MUSIC算法无法直接应用于圆阵，为了解决这一问题可以利用相位模式激励法使圆阵阵列流形具有类线阵的形式(结构)，这样谱峰搜索过程也可以用多项式求根来代替，形成UCA-ROOT-MUSIC算法。ROOT-MUSIC算法扩大了估计范围,提高了估计精度,无需谱峰搜索,就可得到来波方向的估计值，扩大了估计范围,提高了估计精度。本章主要研究ROOT-MUSIC算法，首先简单回顾一下直线阵ROOT-MUSIC算法，然后使用相位模式激励的技术将ROOT-MUSIC算法应用于圆阵。通过仿真结果可得出ROOT-MUSIC算法具有比MUSIC算法更优的性能［1］。4.1ROOT-MUSIC算法的基本原理如图4.1所示，均匀线阵有N个各向同性的全向天线阵元，有P(P〈N)个互不相关的远场窄带平面波信号入射到此阵列上，入射方向分别为｛0,0,,0｝，噪声为互不相关的高TOC\o"1-5"\h\z12P斯白噪声，与信号源不相关，方差为b2，则第i个阵元的观测信号为x(t)=£s(t)exp[j2冗(d/九)(i一l)sir0]+n(t)，其中s?(t)为信号的包络，d为阵元ikiikk二1间距，n.(t)为第i个阵元上的噪声。将上式写成向量形式：ix(t)=As(t)+n(t)4-1V图2.1灼匀线阵(ULA)模型第p牛窄带平面波＞：c■7：XLCt)其中x(t)=[x(t),x(t),•••,x(t))]T为观测数据矢量12Ns(t)=[s(t),s(t),•••,s(t))]T为信号矢量12Pn(t)=[n(t),n(t),…,n(t))]t为噪声矢量，通常假定n(t)为独立的零均值白噪声，12Ni且与信号源无关。对于窄带信号，矩阵A是9的函数，有A(9)=[a(9),a(9)，…,a(9)]12Pa(9)=[1,e-j入sin9i，…,e-j入(N9「]t(i=1,2,…,P)为第i个信源的导向矢量。角i度超分辨算法所要解决的问题即为基于(4-1)式所描述的阵列接收信号来估计信号源的到达方向{9,9,…,9}。TOC\o"1-5"\h\z12P用上面的数学模型，可得输入协方差矩阵R，可 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 示为：XXR=EXXx0x1L*x*01x*]M—1=E[XXH]=AE[ssH]AH+E[nnH](3-8)xM—1TOC\o"1-5"\h\zR=ARAh+b21(3-9)XXssn其中R为信号自相关矩阵，R=E[ssh]。设R的特征值为[入，入,…,入],则由式ssssXX12P(3-9)，易知ARAh的特征值v为v八—b2(i=l,2，…,P)。ssiiin由于矩阵A的列向量(即响应信源的导向矢量)线性独立，于是其列满秩。只要信源之间不是强相关，信号自相关矩阵R为非奇异矩阵。阵列流形矩阵A满秩和R为非奇异ssss矩阵保证当信源个数P小于阵元个数N时，NxN阶矩阵ARAh半正定且秩为P。根据矩ss阵理论，易得ARAh有N-P个特征值v=0；R有N-P个特征值为b2，若按降序排列ssiXXn特征值，则有九,……,九=b2PN—Pn由上面对R的分析，对R进行特征值分解可得：XXXXRV=VM(3-10)xx其中，M=diag[久,久，…,久],久>久>……>久为特征值，diag表示对角阵，01M—101M—1V=[q,q，…,q]为相对应的特征向量。利用最小特征值九的重数K估计信号数01M-1minP=N-K。则相关矩阵的特征向量空间分为两个正交子空间。信号子空间为P个大的特征值对应的P个特征向量，V={q,q,…,q}。噪声子空间为N-P个小的相等特征值对应的N-Ps12P个特征向量，V={q,q，…,q}。设对应于特征值九的特征向量为q，有TOC\o"1-5"\h\znP+1P+2N—Pii(R-九I)q=0;将K个最小特征重值b2代入，则有XXin(R—b21)q=0=ARAhq+b21-b21(3.11)XXnissinn得到ARAHq=0，又因为A满秩且R为非奇异矩阵，故有Ahq=0;说明对应于N-Pssissi个小特征值的特征向量与阵列流形A的P个导向矢量正交。由此则可通过在所有可能的阵列导引向量中搜索与噪声子空间正交的向量，就可确定DOA角度。定义MUSIC空间谱为：P(9)=1(3.12)musicaH(9)EEHa(9)nn当9为信号的波达方向角时，aH(9)EEHa(9)=0。这种正交性使分母达到最小，于nn是通过谱峰搜索，便可求出DOA，这便是著名的MUSIC算法。如果令z=ej2兀dxsin(9丿/人，则有a(9)=[1,z-1,…，zN-1]T，上式为z的一个多项式，可通i过多项式求根得到波达方向，这就是ROOT-MUSIC算法。它将MUSIC算法的谱峰搜索转化为多项式求根,从而不再需要谱峰搜索就可以直接得到准确的DOA估计，使得计算量减小，大大的节省了MUSIC算法的估计时间，成为有效的DOA估计算法。下面给出ULA-ROOT-MUSIC算法步骤：计算接收数据协方差矩阵R。R矩阵的(i,k)元素通过下面的方法求出x(t)=As(t)+n(t)R=E[xxh]=AE(ssH)Ah+b21=ARAH+b21ss对R进行特征值分解，在由R所张成的实空间中，可以使特征值分解得到信号子空间E和噪声子空间估计E。噪声子空间矩阵为E={e，…,e}。snnP+1N解多项式p(z)=aH(z)EEHa(z)的P个根,其中nna(z)=[1,z-1,…,zN-1]tz=ej2兀d”sin(9；)/入p(z)=p+pz+…+p*zP-1+p*zP4)0110根据求得的p(z)选择P个信号根zi,i=1,2,…,P。求出目标方位估计值为0=arcsin(i九arg(z))2兀di仿真1：八元均匀直线阵。阵元间距选半波长，两个非相关窄带信号，第一个信号的入射来自30°,第二个信号的入射来自34。。信噪比均为20dB,比较MUSIC与ROOT-MUSIC算法分辨率。-0.8-0.6估计出来的角度为(29.9714，34.0114)分析与讨论：从结果中可以看到，当信噪比(SNR=20)比较高的时候，两个角度相差4度，MUSIC与ROOT-MUSIC算法都可以把角度分辨出来，说明他们在大信噪比的条件下分辨率相当。仿真2：八元均匀直线阵。阵元间距选半波长，两个非相关窄带信号，第一个信号的入射来自30°，第二个信号的入射来自34。。信噪比均为5dB,比较MUSIC与ROOT-MUSIC算法分辨率。FH=WHFH=WHCVH5-28)结论：当信噪比为5dB时，MUSIC算法对30,34的角度分辨已失效，而ROOT-MUSIC对角度的分辨还相当的明显，并且很精确。ROOT-MUSIC算法的性能明显优于MUSIC算法。从上图可以清楚看到，ROOT-MUSIC算法对MUSIC算法无法分辨的十分靠近的信号仍然具有分辨能力。4.2UCA-ROOT-MUSIC算法本文把ROOT-MUSIC算法推广到圆阵，利用相位模式激励可以把圆阵输出转换到相位模式空间输出，圆阵在相位模式空间可以等效为一个虚拟线阵结构，使得圆阵的阵列流形向量具有范德蒙结构，这样，通过相位模式转换，ROOT-MUSIC算法就可以很自然地应用到圆阵了。计算机仿真结果表明，UCA-ROOT-MUSIC算法运算量低，对目标方位的估计精度高于UCA-RB-MUSIC算法，如图2.10所示，具有N个阵元的均匀圆阵，其N个各向同性的全向天线阵元均匀分布在半径为r的圆上，选择圆心是接收天线阵的参考点，阵面与目标共面，即仰角为90度。虚拟线阵结构：为了提高波束空间MUSIC算法性能，利用上一章节的实波束空间投影矩阵Fh，fh。er首先Fh波束空间阵列流形为：ea=FHa=CVHa=i：NJv®)(5-25)eee上式中v(0)=［e-jM©,…,e-,ej0,e•神,…，ejM©］即为虚拟线阵结构。波束形成器FH是FH乘以一个具有中心希尔伯特特性的矩阵WH形成的，FH使其波rer束空间具有实值阵列流形［10］。形成的实值阵列流形矢量为：a=FHa=\NWhJv®)(5-29)rrq我们可以通过引入正交波束形成器FH，在求解信号子空间矩阵S时，使求解过程简化。r同时也可以提高估计精度。步骤如下：通过FH将阵元接收信号响应x(t)=As(t)+n(t)转变为：ry(t)=FHAs(t)+FHn(t)=As(t)+FHn(t)rrbr此时DOA矩阵由原阵元空间阵列流形矩阵A转换为实值波束空间阵列流形矩阵A=FHA。br求解FH波束空间的相关矩阵：rR=E[y(t)yH(t)]=APAT+bI=CVhARVCh+a2i，bbvsv2M+1取相关矩阵实值R=real(R)做特征值分解，得到噪声子空间估计矩阵G。e实导引向量a(q,©)=F5(q,Q)= 公式 小学单位换算公式大全免费下载公式下载行测公式大全下载excel公式下载逻辑回归公式下载 可以仅仅借助前向协方差矩阵获得11C=(QHRQ+QHJR*JQ)=(QHRQ+(Q*)HJR*JQ*)=Re{QHRQ}22特征值分解：C=ErEh=ErEH+b2EEhSSSNNE={e,…,e}S1PE={e,…,e}NP+1Nrs=diag{Y1,…，丫p}所以有：PU-MUSIC(z)=aT(1/z)EEHa(z)其中a(z)q=QHa(z)实值的协方差矩阵为：AC=QHARAHQ+Q2QHQ=ssARss这里：A(z)=QHA(z)qUnitary-ROOT-MUSIC算法步骤如下：1)通过FH将阵元接收信号响应x(t)=As(t)+n(t)转变为：y(t)=FHx(t)=FHAs(t)+FHn(t)=As(t)+FHn(t)eeee2)求解FH波束空间的相关矩阵：eR=E[y(t)yH(t)]=APAT+QI=CVHARAHVCH+Q2Ieeese2M+13)实值协方差矩阵为：NqC=QHRQ3)取相关矩阵实值R=real(C)做特征值分解，得到噪声子空间估计矩阵G。e实导引向量a(g,©)=QHa(g,©)=Qhjv3)与噪声子空间解多项式qgp(z)=aH(g,©)GGHa(g,©)=vH(©)J(g)[QGGHQH]J(g)v(©)qqv(©)=[e-jM©，…,e-j©,ej0,ej©,…,ejM©]=[z-M,…,z-1,1,z1,…,zM]，z=ej©p(z)=p+pz+…+p*zP-1+p*zP0110根据求得的p(z)选择P个信号根，z,i=1,2,…,P。求出目标方位估计值为i©=arg(z)i=1,2,…,Pii我们可以得出Unitary-ROOT-MUSIC算法整个的计算量只是传统ROOT-MUSIC算法计算量的1/4我们也要指出这个结论不能推广至谱峰搜索MUSIC算法,因为谱峰搜索MUSIC的计算量主要是花费在谱峰搜索上，而不是在特征值分解上。在这里我们对比一下Unitary-ROOT-MUSIC算法和波束空间ROOT-MUSIC算法，后者依靠中心阵元的假设(而这种假设是不容易达到的)，而前者不需要这样的假设。再者，Unitary-ROOT-MUSIC算法忽略了波束空间协方差矩阵的复杂性，最后前后向平滑ROOT-MUSIC算法Unitary-ROOT-MUSIC算法在性能上是一致的，他们都可以应用到有相关源信号的环境中，而波束空间ROOT-MUSIC算法则不能应用。仿真1：由波束形成器Fh形成Unitary-ROOT-MUSIC算法估计DOAr仿真条件：我们选用12元均匀圆阵，r/九=0.5，最大模式数M=6；选用的两个入射信号DOA分别为30°，40°；快拍数为100。采用DOA估计定位误差的均方差(RMS)作为精度的度量。信噪比从5dB到40dB变化，每隔2dB选一点。每个选定信噪比经10次Mote-Carlo独立计算平均得出结果。UCA-Unitary-ROOT-MUSIC(mITIPNELijSZH结论5dB时，UCA-ROOT-MUSIC算法对30,40的角度分辨已失效，而Unitary-ROOT-MUSIC算法在信噪比为5dB时，对角度的分辨还相当的明显，并且很精确。在10dB时，UCA-ROOT-MUSIC算法对30，40度的均方误差为0.9，1.1；而Unitary-ROOT-MUSIC算法在信噪比为10dB时，对同样的角度的均方误差为0.5，0.7,Unitary-ROOT-MUSIC算法性能明显优于UCA-ROOT-MUSIC算法。仿真2：小角度