互质面阵下降维Capon的二维DOA估计方法

摘要

本发明公开了在互质面阵下对多个信号进行二维DOA估计的降维Capon算法。本发明算法先求得接收信号的协方差矩阵，然后把二维谱函数降维到一维，利用真实和模糊DOA估计之间的线性关系，在局部区域上进行一维部分谱搜索，就可以得到一个与理论DOA相关的模糊DOA估计，从而显著降低了计算复杂度，最后，根据互质面阵的性质可以得到真实的DOA估计。与传统的互质面阵二维Capon算法相比，该方法可以获得几乎相同的DOA估计性能，同时又避免了二维谱峰搜索带来的巨大运算量，大大降低了计算复杂度。

说明书

技术领域

本发明涉及波达方向估计领域，具体的为互质面阵下降维Capon的二维DOA估计方法。

背景技术

多个窄带信号源的波达方向估计是阵列信号处理中的一个基本课题，对无线通信、医学成像、雷达、声纳以及导航等各种工程应用也有着重大的意义。互质阵列因为可以减少互耦效应，增加自由度，提高分辨率，引起了人们的极大关注。人们也相继提出了各种互质阵列的DOA估计方法。为了消除相位模糊问题，有基于类投影的无搜索DOA算法，可以大大降低计算的复杂度，还可以用谱峰搜索方法，利用两个子阵的结果和互质阵性质得到波达方向。

在估计二维信号源的波达方向时，可以把互质对称稀疏交叉阵列用于DOA估计，并通过构造高阶矩阵给出了一种有效的算法。也可以将传统的二维多信号分类算法(MUSIC算法)应用于互质平面阵列，这样的全谱搜索方法具有很好的DOA估计性能，但是计算复杂度很高，部分谱搜索方法虽然降低了一定的复杂度，但是因其仍然需要二维谱峰搜索，计算量依然很大，在实际工程中耗费的成本很高。

发明内容

本发明所要解决的技术问题针对背景技术中所涉及的二维谱峰搜索计算复杂度过大的问题，采用对谱函数降维处理的方法来避免二维谱峰搜索带来的巨大计算复杂度，减少算法的时间和空间成本，同时能够保持其良好的估计性能。

本发明为解决上述技术问题采用以下技术方案：互质面阵下降维Capon的二维DOA估计方法：

步骤1,将互质面阵分为两个均匀面阵，接下来分别考虑这两个均匀面阵，计算接收信号的协方差矩阵

步骤2，在v∈[0,2/M

步骤3，根据估计的v

步骤4，根据真实波达方向与模糊波达方向之间的相位差计算出所有可行的(u,v)对，再计算差分函数d

步骤5，求出信源方向的仰角和方位角得到二维DOA估计。

优选的，所述步骤1中构造协方差矩阵以及Q(v)的具体步骤如下：

步骤1.1，将互质面阵分为两个均匀面阵，接下来分别考虑这两个均匀面阵，接收信号的协方差矩阵

步骤1.2，二维Capon的空间谱函数为：

优选的，所述步骤2中一维谱峰搜索求取v

步骤2.1，

构造下面的代价函数

根据

所以

此外，我们还可以得到v

步骤2.2，假设只有一个信源入射到M

2πd

为了简化，还可以表示为

可以知道，在变换域中每个真实角度参数与相应的模糊角度参数之间具有线性关系。对于第i个子阵而言，差值为

在不丧失一般性的前提下，对于第i个子阵，我们选择在v∈[0,2/M

优选的，所述步骤3中最小二乘法求解u

根据

g

现在用最小二乘法来估计u

优选的，所述步骤4求解

首先利用线性特征，计算出其他所有可行的模糊(u,v)估计。

虽然当阵元间距大于半波长时，会出现相位模糊的情况，但是利用互质阵的特性，可以得到真实的波达方向估计。根据

定义差分函数d

其中(u

优选的，所述步骤5中求解信源方向的仰角和方位角的具体步骤如下：

可得信源方向的仰角和方位角分别为

其中k＝1,2,…,K。

有益效果：本发明采用以上技术方案与现有技术相比，具有以下技术效果：

1.该算法利用真实和模糊DOA估计之间的线性关系，在局部区域内只需要进行一维部分搜索，显著降低了计算成本。

2.该算法的估计性能与互质阵中的全局搜索或局部搜索性能几乎相同，但是计算复杂度要低很多。

3.该算法可以实现自动配对的二维角度估计。

附图说明

图1为本发明提供的互质面阵下降维Capon的二维DOA估计方法的实现流程图。

图2为互质面阵的阵列结构模型图。

图3为本发明算法估计的波达方向散点图。

图4为各种不同算法下仰角的均方根误差比较图。

图5为各种不同算法下方位角的均方根误差比较图。

图6为不同互质面阵参数下所提算法估计仰角的均方根误差比较图。

图7为不同互质面阵参数下所提算法估计方位角的均方根误差比较图。

图8为不同快拍数下所提算法估计仰角的均方根误差比较图。

图9为不同快拍数下所提算法估计方位角的均方根误差比较图。

具体实施方式

下面结合附图对本发明做更进一步的解释。

如图1所示为互质面阵下降维Capon的二维DOA估计方法的实现流程图，由图2的阵列结构可知互质面阵可以分为两个均匀子阵来考虑，先求得接收信号的协方差矩阵，然后把二维谱函数降维到一维，利用真实和模糊DOA估计之间的线性关系，在局部区域上进行一维部分谱搜索，就可以得到一个与理论DOA相关的模糊DOA估计，从而显著降低了计算复杂度，最后，根据互质面阵的性质可以得到真实的DOA估计。

如图3所示是使用所提算法DOA估计的散点图结果。仿真中假设互质面阵由两个具5×5和4×4个阵元的均匀平面子阵列组成，其中d

如图4和图5所示为本发明所提的互质面阵中降维Capon算法与其他相关算法的误差比较情况。为了公平比较，均匀面阵设计按照与互质面阵阵元总数40相同的间距为半波长的5×8阵列排布。定义均方根误差为

如图6和图7所示为不同M

如图8和图9所示为不同的快拍数L下的DOA估计性能对比，其中M

以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。