一种基于阵列信号处理的波达方向估计方法

摘要

本发明公开了一种基于阵列信号处理的波达方向估计方法，所述方法针对室内定位场景，提出基于旋转不变子空间信息的循环平稳信号正交匹配追踪算法。引入循环平稳算法作为初始相关函数，将时间刻度上的循环平稳性与空间刻度上的阵列信号结合起来，从而有效解决正交匹配追踪算法对噪声敏感，抗噪性能不好的问题；采用旋转不变子空间信息，从而突破正交匹配追踪算法存在的瑞利极限，有效提高多物体定位分辨力并降低计算复杂度；采用基于压缩感知理论的正交匹配追踪算法，通过稀疏信号重构进一步提高算法的抗干扰性。满足了室内定位场景高精度，低时延的要求。

说明书

技术领域

本发明属于波达方向定位技术领域，尤其涉及一种基于阵列信号处理的波达方向估计方法。

背景技术

阵列信号处理作为信号处理领域一个重要的研究方向，近年来得到了广泛的关注与迅速的发展。该项技术通过在三维空间中按一定流型安放的信号接收传感器阵列来获得信号的离散观测数据，并利用不同传感器接收数据的差异，设计相应的算法以获得信号的各项参数。近年来，阵列信号处理已被广泛应用于声呐通信、射电天文、生物医学、定位导航等众多民用或军用领域。

阵列信号处理的研究主要包括自适应空域滤波和空间谱估计两个方面。其中，空间谱估计也称为波达方向(Direction Of Arrival,DOA)估计。DOA估计问题是根据阵列参考天线接收到的各个信号信息确定同时位于空间某片区域内多个信号的空间位置，也即是各个信号相对于阵列参考天线的方向角。DOA估计的分辨率和天线阵列的长度有关，当天线阵列的长度确定之后，其对应的分辨率也就被确定，称其为瑞利限。对于瑞利限而言，存在能够突破瑞利限限制的方法，这类方法称为超瑞利限方法或超分辨方法。

早期的DOA估计研究只研究空间来波信号的一维到达角度，接收阵列通常采用均匀线阵。最早提出的是常规波束形成法。随后，许多能够突破“瑞利限”且可以实现更高分辨率的估计方法先后被许多人提出来。其中以多重信号分类(MUSIC)算法最为出名。在1989年,旋转不变技术(ESPRIT)算法被Roy R等人提出来。20世纪80年代，出现了诸如最大似然(ML)算法、子空间拟合(SF)算法等子空间类拟合算法。这些算法与MUSIC算法和 ESPRIT等算法相比，估计的角度值更加准确，但是运算复杂度也更大。

近年来，压缩感知(CS)理论的发展也为阵列信号处理中DOA估计提供了新思路。CS理论跳出奈奎斯特定理的限制，当以低于奈奎斯特频率的采样频率对信号采样时，其能从维数不多的采样信号中完整恢复出原有信号。CS理论被应用在在空间谱估计中时，接收信号在整个空域来说可以看成是以低于奈奎斯特频率进行采样出的信号，然后通过求解最优解的方式即能找出信源到达角度。与传统的DOA估计方法相比基于压缩感知理论的DOA(CS-DOA)估计具有众多优点，包括采样率低使得采样时间更短，抗干扰能力强等优势。目前，CS-DOA估计仍处于不断探索的阶段，但其理论研究价值和实际工程意义不可估量。

发明内容

本发明目的在于提供一种基于阵列信号处理的波达方向估计方法,通过天线阵列获得来波采样数据，由循环平稳算法计算循环相关矩阵，接着利用奇异值分解获得对应的信号子空间，构建稀疏重构模型，然后基于旋转不变子空间信息获得每次迭代残差对应的字典向量，将字典向量匹配后加入支撑集中计算更新残差，不断迭代直至计算出所有来波方向，以解决不能满足室内定位对精度和时延所提出要求的技术问题。

为解决上述技术问题，本发明的具体技术方案如下：

一种基于阵列信号处理的波达方向估计方法，包括以下步骤：

步骤1、通过天线阵列获得来波采样数据；

步骤2、由循环平稳算法计算循环互相关矩阵；

步骤3、利用奇异值分解获得对应的信号子空间，构建稀疏重构模型；

步骤4、基于旋转不变子空间信息获得每次迭代残差对应的字典向量；

步骤5、将字典向量匹配后加入支撑集中计算更新残差；

步骤6、不断迭代直至计算出所有来波方向，重复步骤4和步骤5，不断迭代直至计算出所有来波方向。

进一步的，步骤1具体包括以下步骤：采用均匀线阵对信号进行接收，对于相隔d共有N个阵元的均匀线阵，f

进一步的，步骤2的计算方法包括以下步骤：

步骤2.1、计算第k个阵元和第i个信号源的时延τ

步骤2.2、由循环互相关的定义及阵元接收数据可得阵元s

步骤2.3、由

进一步的，所述步骤2.1中，

步骤2.2中，

其中

由于噪声循环频率不为α，其在α频率上的相关系数为0，得到循环互相关函数为：

其中s

步骤2.3中阵列接收信号的循环互相关矩阵计算为：

进一步的，步骤3包括以下步骤：

对步骤2所得循环互相关矩阵

进一步的，所述步骤3中，循环互相关矩阵

根据

其中

其中U

其中

进一步的，步骤4计算方法包括如下步骤：

步骤4.1、初始化状态支撑集S

步骤4.2、设此时为第k次迭代：支撑集为S

其中A(θ)表示接收信号角度矩阵，A

设过渡矩阵Ψ，其中A

进一步的，步骤5包括以下步骤：设第k次迭代的所选择的向量角标为i

进一步的，所述步骤5中，更新支撑集和稀疏表示矩阵：S

正交匹配追踪的残差迭代方程表示为：

其中下标k表示第k次迭代，B表示残差迭代方程的未知数，B

本发明的一种基于阵列信号处理的波达方向估计方法，具有以下优点：

本发明基于循环平稳算法，旋转不变子空间信息以及正交匹配追踪算法，得到一种基于阵列信号处理的波达方向估计方法。利用循环平稳算法，将时间刻度上的循环平稳性与空间刻度上的阵列信号相结合，使得算法的抗噪性能增强。所述方法通过旋转不变子空间信息与基于压缩感知理论的正交匹配追踪算法的结合，突破正交匹配追踪算法存在的瑞利极限，有效提高多物体定位分辨力并降低计算复杂度，并进一步提高算法的抗干扰性。达到了室内定位场景高精度，低时延的要求，改善用户体验。

附图说明

图1为本发明的基于阵列信号处理的波达方向估计方法的场景示意图；

图2为本发明的基于阵列信号处理的波达方向估计方法流程图；

图3为本发明的基于阵列信号处理的波达方向的Cyclic-正交匹配追踪算法流程图；

图4为本发明的基于阵列信号处理的波达方向的基于旋转不变子空间信息的精简字典正交匹配追踪算法流程图。

具体实施方式

为了更好地了解本发明的目的、结构及功能，下面结合附图，对本发明一种基于阵列信号处理的波达方向估计方法做进一步详细的描述。

如图1所示，本发明的研究基于室内定位场景。其中绑定在所需定位物体上的发信器用于发送信号，多个不同位置的天线阵列用于采样接收信号，处理器用于根据采样信号进行操作以执行所述的方法，将所得角度输入PC处理端后，根据多个角度数据获得高精度定位信息。

如图2所示，本发明设计了一种基于阵列信号处理的蓝牙波达方向的方法，通过该方法对信号进行处理，获得对应的来波方向，用于实现高精度，低时延的室内定位。实际应用当中，具体包括如下步骤。

步骤1、通过天线阵列获得来波采样数据。我们采用均匀线阵对信号进行接收，对于相隔d共有N个阵元的均匀线阵，当仅有一个信号源入射，忽略噪声时，其第一个阵元接收的信号为：

其中，f

当有M个循环频率为α的信号源入射时(M≤N)，第k个阵元在第t个快拍处接收的信号为：

其中

即：S(t)＝A(θ)X(t)+N(t)，t＝1,2,......,U。其中，U是快拍的数量， A(θ)＝[a(θ

在步骤2和步骤3中，阐述了本发明设计的基于阵列信号处理的蓝牙波达方向Cyclic-正交匹配追踪算法，如图3所示。

步骤2、由循环平稳算法计算循环相关矩阵。本算法采用了Cyclic算法进行计算。设时延τ

其中p＝1,2,…,N,q＝1,2,…,N,

由于噪声循环频率不为α，其在α频率上的相关系数为0，得到循环互相关函数为：

其中s

由于

步骤3、利用奇异值分解获得对应的信号子空间，构建稀疏重构模型，循环相关矩阵

根据

其中

其中U

其中

在步骤4，步骤5和步骤6中，阐述了本发明设计的基于阵列信号处理的蓝牙波达方向的基于旋转不变子空间信息的精简字典正交匹配追踪算法，如图4所示。

步骤4、基于旋转不变子空间信息获得每次迭代残差对应的字典向量。

设初始化状态支撑集S

设此时为第k次迭代：支撑集为S

令Span(A(θ))、Span(U

由残差迭代方程可知Span(C

由上述两式可得Span(C

由上式可得A(θ)·T＝C

其中A

同理，设

由上述三式可得

由上式可得

由上式可得

由于T

可得字典向量：

步骤5、将字典向量匹配后加入支撑集中计算更新残差。设第k次迭代的所选择的向量角标为i

S

OMP(正交匹配追踪)的残差迭代方程表示为：

其中下标k表示第k次迭代，B表示残差迭代方程的未知数，B

步骤6、不断迭代直至计算出所有来波方向。重复步骤4和步骤5，不断迭代直至计算出所有来波方向。

可以理解，本发明是通过一些实施例进行描述的，本领域技术人员知悉的，在不脱离本发明的精神和范围的情况下，可以对这些特征和实施例进行各种改变或等效替换。另外，在本发明的教导下，可以对这些特征和实施例进行修改以适应具体的情况及材料而不会脱离本发明的精神和范围。因此，本发明不受此处所公开的具体实施例的限制，所有落入本申请的权利要求范围内的实施例都属于本发明所保护的范围内。