技术领域
[0001] 本
发明属于
信号处理技术领域，尤其涉及一种波达方向估计方法，用于对目标方位的估计。
背景技术
[0002] 波达方向DOA估计是雷达、声纳
信号处理中的一个重要研究方向。众所周知，对于一个具有N个阵元的均匀线性阵列，采用传统的波达方向估计方法，如MUSIC，ESPRIT等所能分辨的最大目标个数为N-1个。而欠定波达方向估计问题，即目标个数大于阵元数的DOA估计经常出现并引起广泛的研究兴趣。解决该问题的其中一种有效方法是利用一个等效的虚拟阵列来提高波达方向估计的
自由度，该虚拟阵列是通过对一个特殊设计的非均匀线性阵列接收信号的协方差矩阵向量化来构造的。最小冗余阵列MRA、嵌套式阵列NA、互质阵列CA及嵌套式最小冗余阵列NMRA等都是针对欠定波达方向估计而设计的非均匀阵列。但是在利用阵列接收数据的协方差矩阵构造虚拟阵列时，虚拟阵列中等效信源被实际信源的功率所替代，因此这些等效信源就表现为完全相干的信号，这样就无法直接利用传统的DOA估计方法来进行波达方向估计。
[0003] 近来，许多利用非均匀线阵的虚拟阵列来解决欠定波达方向估计问题的新
算法被提出。其中一种代表性方法是P.Pal和P.P.Vaidyanathan在《Nested arrays:Anovel approach to array processing with enhanced degrees of freedom》提出的空间平滑(SS)-MUSIC算法，但该方法在进行DOA估计时没能用到虚拟阵列的所有阵元，从而导致一定的
信噪比损失。另一种具有代表性的方法是Y.D.Zhang、M.G.Amin和B.Himed在《Sparsity-based DOAestimationusing co-prime arrays》提出的
稀疏信号重构法，该方法利用了信号谱的稀疏性来进行DOA估计，然而实现这种方法的算法通常需要巨大的计算量。
发明内容
[0004] 本发明的目的在于克服上述已有DOA估计算法的不足，提出一种新的非均匀阵列欠定波达方向估计方法，以充分利用虚拟阵列的数据，减少计算量。
[0005] 为实现上述目的，本发明的技术思路是:利用接收阵列数据的协方差矩阵来构造虚拟阵列；利用全部虚拟阵元的数据来获得差分合成阵列的观测数据，而不是删除重复的阵元数据；通过构造一个观测数据的Toeplitz矩阵代替传统的空间平滑操作来解相干，并估计信号和噪声的子空间；通过构造一个关于MUSIC谱的多项式并利用其根来估计DOA，从而避免了大运算量的
角度栅格搜索。其实现步骤包括如下：
[0006] 1)根据阵列接收数据x(t)估计非均匀阵列的协方差矩阵 并对该该协方差矩2
阵 向量化，得到K×1的向量
[0007]
[0008] 其中 t表示
采样时刻，t＝1,2,…,N，N表示快拍数，(·)H表示共轭转置，K表示阵元数；
[0010] R＝(ETE)-1ET
[0011] 其中 是一个在第pl处为1，其余为0的2
fV×1向量，l＝1,2,…,K，i＝1,2,…,K，fV表示非均匀阵列的自由度，fV＝2lV+1，lV是T
非均匀阵列的孔径长度，符号 表示复数域，(·)表示转置；
[0012] 3)根据降维矩阵R和向量 计算非均匀阵列的差分合成阵列的观测数据z：
[0013]
[0014] 4)由观测数据z，构造Toeplitz矩阵Y：
[0015]
[0016] 其中zm是观测数据z的第m个元素，m＝1,2,…,fV；
[0017] 5)对Toeplitz矩阵Y进行特征值分解，即：
[0018]
[0019] 其中US为信号子空间，ΛS表示信号的特征值，σ2表示噪声功率，UN表示噪声子空间；
[0020] 6)根据5)所获得的噪声子空间UN，构造root-MUSIC多项式fMUSIC(r)：
[0021]
[0022] 其中r是关于目标到达角θ的未知量，L(r)表示r的lV次多项式所构成的向量；
[0023] 7)求解6)所构造的多项式fMUSIC(r)，得到多项式的解rn，n＝1,2,…,Q，Q为目标数，根据rn，进而获得目标到达角θn：
[0024]
[0025] 其中arg(·)表示取相角运算，d为阵元间隔。
[0026] 本发明与现有DOA估计算法相比具有以下优点：
[0027] 1.本发明综合利用了所有虚拟阵元接收到的数据，避免了接收信号的信噪比SNR的损失；
[0028] 2.本发明通过构造观测数据的Toeplitz矩阵代替传统的空间平滑操作来解相干，节省空间平滑解相干的运算量；
[0029] 3.本发明通过构造一个关于MUSIC谱的多项式并利用其根来估计DOA，从而避免了普通谱MUSIC算法复杂的角度栅格搜索，在减少计算复杂度的同时还能获得更好地DOA估计性能。
附图说明
[0031] 图2是本发明仿真使用的12阵元嵌套式最小冗余阵列NMRA的几何结构示意图；
[0032] 图3是用本发明和SS-MUSIC算法对波达方向DOA的估计均方根误差与信噪比关系的比较图。
具体实施方式
[0033] 下面结合附图详细说明本发明的内容和效果。
[0034] 参照图1，本发明的实现步骤如下：
[0035] 步骤1：根据阵列接收数据x(t)估计非均匀阵列的协方差矩阵
[0036] 1a)对于给定K个阵元的非均匀线性天线阵列，设阵元
位置为：
[0037] v＝[v1,v2,…,vi,…,vK]d,
[0038] 其中vi为第i个阵元位置系数，i＝1,2,…,K，d为阵元间隔，取值为半个信号
波长；
[0039] 1b)将目标n从角度θn入射到阵列的导向矢量表示为：
[0040]
[0041] 其中 表示第n个信号在第i个阵元处的阵因子，n＝1,2,…,Q，Q为信号数T目，κ为半
波数， (·)表示转置；
[0042] 根据入射信号的导向矢量a(θn)，构造导向矢量矩阵A：
[0043] A＝[a(θ1),a(θ2),…,a(θn),…,a(θQ)]，
[0044] 1d)设Q个不相关的窄带信号从角度θn入射到阵列，将阵列的接收数据表示为：
[0045] x(t)＝As(t)+n(t)，
[0046] 其中s(t)表示Q个不相关的窄带信号，s(t)＝[s1(t),s2(t),…,sn(t),…,sQ(t)]T，sn(t)表示第n个入射信号，各信号在时间上相互独立且非相关，且服从复高斯分布表示第n个信号sn(t)的功率；n(t)表示均值为0，方差为σ2高斯白噪声，其满足独立同分布，并且与入射信号不相关；t表示采样时刻，t＝1,2,…,N，N表示快拍数；
[0047] 1e)根据阵列接收数据x(t)，估计非均匀阵列的协方差矩阵
[0048]
[0049] 其中，IK是K阶单位方阵，(·)H表示共轭转置。
[0050] 步骤2，向量化阵列协方差矩阵 计算虚拟阵列的接收数据
[0051]
[0052] 其中B表示虚拟的差分合成阵列的方位矢量矩阵，
[0053] c表示等效信源矢量， 1n为单位矩阵， 表示第i
个元素为1，其余元素为0的单位列向量；符号vec(·)表示矩阵的向量化操作，符号 代表Kronecker积，(·)*表示共轭。
[0054] 步骤3：构造降维矩阵R。
[0055] 3a)规定单位向量
[0056]
[0057] 其中 在第pl处为1，其余为0的单位向量，pl表示虚拟阵元位置，pl＝2
1,2,…,fV，l＝1,2,…,K，fV表示非均匀阵列的自由度，fV＝2lV+1，lV是非均匀阵列的孔径长度；
[0058] 3b)根据3a)所规定的向量 生成单位矩阵E：
[0059]
[0060] 其中E是K2×fV的单位矩阵；
[0061] 3c)由3b)所生成的单位矩阵E，得到降维矩阵R：
[0062] R＝(ETE)-1ET；
[0063] 其中R是fV×K2的矩阵。
[0064] 步骤4：根据接收数据 和降维矩阵R，计算非均匀阵列的差分合成阵列观测的数据z：
[0065]
[0066] 这里根据观测数据z得到了所有虚拟阵元的数据，但并不删除其中重复的量。
[0067] 步骤5：由观测数据，构造Toeplitz矩阵Y。
[0068] 5a)展开观测数据z的各个元素zm，m＝1,2,…,2lV+1；
[0069] 5b)根据5a)中的zm，从m＝1开始，选取连续的lV+1个元素，生成第一个向量m逐次后移，生成第m个向量 直至m＝lV+1，最后生成第lV+1个向量
[0070] 5c)根据5b)中生成的lV+1个列向量ym，得到Toeplitz矩阵Y：
[0071]
[0072] 这里矩阵Y是个秩为Q的矩阵，如果Q＜lV+1，则等效信源的相干性就可以被解掉，矩阵Y就能被用来正确的估计信号和噪声子空间。由于矩阵Y，能够代替传统的空间平滑操作来解相干，而且由于矩阵Y是对观测数据重新排列来构造的，所以不需要乘法，因而能够节省空间平滑解相干的运算量。
[0073] 步骤6：对Toeplitz矩阵Y进行特征值分解，获得噪声子空间UN。
[0074] 定义矩阵Y的特征值分解为：
[0075]
[0076] 其中 分别为信号和噪声子空间，ΛS表示信号的特征值， λn是第n个信号的特征值，符号Diag{·}表示对角
矩阵。
[0077] 步骤7：根据步骤6所获得的噪声子空间UN，构造root-MUSIC多项式fMUSIC(r)。
[0078] 7a)设关于目标波达方向θ的未知量r：
[0079] r＝ej(2π/λ)dsinθ，
[0080] 其中r＝[r1,r2,…,rn,…,rQ]T，rn表示关于目标n的未知量，θ＝[θ1,θ2,…,θn,…,θQ]T，θn表示目标n的波达方向，λ为信号波长；
[0081] 7b)生成r的lV次多项式L(r)：
[0082]
[0083] 7c)结合L(r)和噪声子空间UN，得到root-MUSIC多项式fMUSIC(r)：
[0084]
[0085] 步骤8：求解步骤7构造的关于MUSIC谱的多项式fMUSIC(r)，利用其根来估计DOA。
[0086] 8a)求解步骤7所构造的多项式fMUSIC(r)，得到其根rn；
[0087] 8b)选取离单位圆最近的Q个根rn，估计目标n的到达角θn：
[0088]
[0089] 其中arg(·)表示取相角运算。
[0090] 通过求解关于MUSIC谱的多项式并利用其根来估计DOA，避免了普通谱MUSIC算法复杂的角度栅格搜索，在减少计算复杂度的同时还能获得更好地DOA估计性能。
[0091] 本发明的效果通过以下计算仿真进一步说明：
[0092] 仿真1：对本发明和SS-MUSIC算法的波达方向DOA估计的均方根误差RMSE与信噪比之间的关系进行仿真。
[0093] 1.1)仿真条件：图2是本发明仿真使用的12阵元嵌套式最小冗余阵列NMRA的几何结构示意图，阵元位置为[0,1,4,6,13,14,17,19,39,40,43,45]d，阵元个数为12，阵元间隔d＝λ/2，目标个数Q＝16，快拍数N＝500，500次Monte-Carlo实验。
[0094] 1.2)仿真内容与结果：
[0095] 用现有的SS-MUSIC算法和本发明的波达方向DOA估计性能进行比较，在上述1.1)仿真条件下，对本发明和SS-MUSIC算法的波达方向DOA估计均方根误差与信噪比关系做仿真，仿真结果如图3所示，其中横轴表示信噪比从-20dB至10dB变化，纵轴表示均方根误差。
[0096] 从图3可以看出，本发明DOA估计性能随着信噪比的提高而提高，而且本发明角度估计的均方根误差小于传统SS-MUSIC，说明本发明具有更好的DOA估计
精度。
[0097] 仿真2：计算复杂度比较。
[0098] 列举非均匀阵列的阵元数分别为12和24时，将本发明和现有的SS-MUSIC算法的计算复杂度进行比较，结果如表1。
[0099] 表1本发明和SS-MUSIC算法的计算复杂度
[0100]
[0101] 注：1)J表示非均匀阵列
波速宽度中角度栅格的个数；
[0102] 2)在12阵元嵌套式最小冗余阵列NMRA中，lV＝45,J＝1588。
[0103] 3)在24阵元嵌套式最小冗余阵列NMRA中，lV＝157,J＝6176。
[0104] 从表1中可以看出，在阵元数为12的情况下，本发明的计算复杂度不足SS-MUSIC算法的一半。而且阵元数越多，复杂度降低越多。在阵元数为24的情况下，本发明的计算复杂度仅为SS-MUSIC算法的计算复杂度的1/3，可见本发明能明显降
低信号处理的运算量，提高目标波达方向估计的计算速度。
[0105] 综上，本发明能够以更低的计算量获得更好的DOA估计精度。
