基于非圆信号的单基地展开互质阵列MIMO雷达DOA估计方法与流程

本发明涉及雷达技术领域，尤其涉及基于非圆信号的单基地展开互质阵列mimo雷达doa估计方法。

背景技术：

近年来，将通信领域的多输入多输出(multipleinputmultipleoutput，mimo)技术思想引入到雷达领域，提出的mimo雷达系统引起了广泛的关注。mimo雷达充分利用信号分集、空域分集增益，获得较传统雷达更高的自由度，其目标检测能力和参数估计能力等均获得广泛的认可。由于mimo雷达潜在的诸多优点，关于mimo雷达的研究快速发展，其中mimo雷达的波达方向(directionofarrival，简称doa)估计问题是其中的重点研究课题。

当前大多数针对mimo雷达doa估计的研究都是建立在均匀线阵的基础上，利用旋转不变性(esprit)或多重信号分类的方式获得相较于传统体制雷达更优异的性能(空间分辨率、抗噪声性能等)。但是由于mimo雷达在发射阵列和接收阵列的阵元间距上一般为半波长，限制了阵列孔径，因此性能受到一定的限制。互质阵列是一种非均匀线阵，相较于均匀线阵，拥有更大的阵列孔径，互质阵空间谱估计算法利用子阵元数目互质特性辨识目标源，因此有更高的侧向精度与分辨率。但是现有方法均是将互质阵列当成两个稀疏的均匀阵列分别进行处理，因此存在以下问题：一是自由度被子阵阵元数目所限制，相较于均匀阵列降低了至少一半；二是需要进一步的运算来消除模糊性问题。展开互质阵列区别于传统互质阵列之处则是将两个子阵按照相反的方向展开排列，阵列孔径进一步扩展，使用所有的阵元接收数据构造空间协方差矩阵，随后使用music算法得到doa估计，阵元数目的互质抑制了相位模糊问题。

基于传统相控阵雷达的doa估计方法往往是利用圆信号构造接收数据矩阵，没有考虑到非圆信号的存在，然而在实际应用中，bpsk、am、mask等非圆信号的使用越来越广泛，因此，充分利用非圆信号的特性进行doa估计是空间谱估计理论的一个重要课题。对于基于mimo雷达的doa估计，一般性方法是将传统方法进行移植，如将子空间类的music方法应用于mimo雷达的二维music方法、降维music方法等，但此类方法往往计算量较大，并且不适用于采用非圆信号测向的方法，当存在相干信源时性能急剧恶化，doa估计性能低或者难以实现doa估计。因此，如何利用非圆信号的特性，进行基于单基地mimo雷达的相干信源doa估计是亟待解决的技术问题。

技术实现要素：

有鉴于此，本发明的目的是提供基于非圆信号的单基地展开互质阵列mimo雷达doa估计方法，将展开互质阵列作为mimo雷达的发射阵列和接收阵列，使得收发阵列的阵列空间大大扩展，并同时将非圆信号作为mimo雷达的发射阵列的发射信号，结合非圆信号的非圆特性、mimo雷达优异的空间分辨率和参数识别性、展开互质阵列扩展的阵列孔径的三重优点，极大提高了doa估计性能。

本发明通过以下技术手段解决上述技术问题：

基于非圆信号的单基地展开互质阵列mimo雷达doa估计方法，所述单基地展开互质阵列mimo雷达的doa估计方法包括以下步骤：

s1.对接收阵元匹配滤波后的数据进行l次采样，得到等价虚拟阵列的接收数据x(t)；

s2.由接收数据x(t)得到其共轭x^*(t)，进而将x(t)与x^*(t)组合为新的扩展接收数据y(t)；

s3.对新扩展接收数据y(t)进行空间协方差矩阵的估计；

s4.对空间协方差矩阵进行特征值分解，并获取噪声子空间；

s5.根据噪声子空间与导向矢量的正交性构造空间谱的计算表达式，按照空间谱的计算表达式进行谱峰搜索，k个谱峰的位置即为k个来波的波达方向估计。

进一步，所述单基地展开互质阵列mimo雷达包含发射阵列和接收阵列，所述发射阵列与接收阵列各包含两个稀疏均匀线阵的子阵1和子阵2，子阵1和子阵2按照相反方向完全展开排列，所述子阵1的最后一个阵元与子阵2的第一个阵元相重合。

进一步，所述单基地展开互质阵列mimo雷达中，设子阵1和子阵2的阵元数目分别为m、n，且m、n互质，m<n，子阵1阵元间距为nλ/2，子阵2阵元间距为mλ/2，λ表示入射到展开互质阵列的电磁信号波长，所述发射阵列与接收阵列各自的阵元数均为m+n-1。

进一步，所述接收数据x(t)如下：

x(t)＝as(t)+n(t)

式中，s(t)是回波信号矢量，a表示单基地展开互质阵列mimo雷达的虚拟阵列流形，n(t)是均值为0，方差为σ²的加性高斯白噪声矢量。

进一步，所述空间协方差矩阵如下：

式中，∑表示求和操作，h表示共轭转置操作。

进一步，所述空间协方差矩阵进行特征值分解获得

其中es是ryy的k个主特征值所对应的特征向量，即信号子空间，en是2(m+n-1)²-k个其他特征值所对应的特征向量，即噪声子空间，us和un分别是k个主特征值以及2(m+n-1)²-k个其他特征值所组成的对角矩阵。

进一步，所述空间谱的计算表达式如下

式中，

本发明相较于现有技术具有以下优点：

(1)通过将传统互质阵列展开成两个相反的方向，得到展开的互质阵列，从而，阵列孔径得到提高，进而将展开互质阵列分别作为发射阵和接收阵，本发明展开互质阵列与mimo雷达相结合，通过利用展开互质阵列的大阵列孔径以及mimo雷达优异的空间分辨率，自由度，以及参数的可识别性等性能，结合经典music算法实现优异的doa估计性能。

(2)本发明不同于传统互质阵列doa估计算法将两个子阵的接收数据分别进行处理，使得自由度下降了至少一半，而是将整个接收阵列作为一个整体，使用所有接收阵元的接收数据来进行doa估计，由于同时利用了阵列的自信息以及互信息，因此自由度以及doa估计精度相较于传统互质阵列doa方法均有大幅的提高。

(3)将非圆信号作为展开互质阵列mimo雷达的发射阵列的发射信号，进而将接收信号与其共轭组合为新的扩展接收信号，使得原接收信号维度倍增，充分利用了非圆信号的非圆特性，可在展开互质阵列mimo雷达优异的doa估计性能下进一步提高性能，在理论上，将自由度提高了一倍。

附图说明

图1是本发明展开互质mimo雷达阵列几何结构图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

一.建立系统数学模型

用2m+2n-1个天线阵元形成展开互质阵列mimo雷达，其中m，n表示两个互质的整数，mimo雷达的接收阵列接收来自k个目标反射的发射阵列发射信号的回波信号，构造展开互质阵列mimo雷达的具体步骤如下：

本发明实施例中的单基地展开互质阵列mimo雷达包含发射阵列和接收阵列，其几何结构如图1所示，发射阵列与接收阵列各包含两个稀疏均匀线阵的子阵1和子阵2，子阵1和子阵2按照相反方向完全展开排列，子阵1的最后一个阵元与子阵2的第一个阵元相重合。设子阵1和子阵2的阵元数目分别为m、n，且m、n互质，不失一般性设m<n。子阵1阵元间距为nλ/2，子阵2阵元间距为mλ/2，其中λ表示入射到展开互质阵列的电磁信号波长，发射阵列与接收阵列各自的阵元数均为m+n-1。

那么，发射阵列和接收阵列的阵元位置可表示为

pt＝pr＝{mnd0|0≤n≤(n-1)}∪{m(n-1)d0+nmd0|0≤m≤(m-1)}(1)

其中d0为半波长，由此可知展开互质阵列的孔径为m(n-1)d0+n(m-1)d0，相较于均匀阵列孔径nd0，传统互质阵列孔径max(m(n-1)d0，n(m-1)d0)可知展开互质阵列的孔径得到了较大的提升。

现各个发射阵元同时发射同频正交的周期相位编码信号，发射信号满足条件公式(2)

其中，si、sj分别是第i个和第j个发射阵元的信号，l为每个重复周期的相位编码个数，∑表示求和操作。

获取发射子阵与接收子阵的导向矢量：

假设存在k个互不相关的远场目标，且满足k<(m+n-1)²，目标的回波波达角度分别为θ1,θ2…θk，因此发射阵以及接收阵关于第k个目标的导向矢量可分别由公式(3)，(4)表示

其中，at1(θk)、at2(θk)分别是发射阵子阵1、子阵2的导向矢量，ar1(θk),ar2(θk)则分别是接收子阵1、子阵2的导向矢量，k＝1,2…k。发射子阵与接收子阵的导向矢量表达式分别为

因此得到发射阵与接收阵的阵列流形at、ar如下

at＝[at(θ1)at(θ2)…at(θk)](7)

ar＝[ar(θ1)ar(θ2)…ar(θk)](8)

进而得到单基地展开互质阵列mimo雷达的虚拟阵列流形

其中

表示khatri-rao积，表示kronecker积。由此可得接收阵元经过匹配滤波之后的接收数据

x(t)＝as(t)+n(t)(11)

其中s(t)是回波信号矢量

βk为第k个点目标的雷达截面系数(radarcrosssection，rcs)，fdk为第k个点目标的多普勒频率，fs为发射波形的脉冲重复频率。n(t)是均值为0，方差为σ²的加性高斯白噪声矢量。

二.基于非圆信号的单基地展开互质阵列mimo雷达doa估计方法

取匹配滤波后的接收信号x(t)的共轭从而得到x^*(t)，进而将x(t)与x^*(t)进行联合，得到扩展的接收信号y(t)

对扩展的接收信号y(t)计算空间协方差矩阵ryy

ryy＝e[y(t)y^h(t)](15)

实际工程中的ryy可由l个有限快拍估计得到。

对空间协方差矩阵ryy进行特征值分解

其中es是ryy的k个主特征值所对应的特征向量，即信号子空间，en是2(m+n-1)²-k个其他特征值所对应的特征向量，即噪声子空间，us和un分别是k个主特征值以及2(m+n-1)²-k个其他特征值所组成的对角矩阵。根据噪声子空间与导向矢量的正交性，可构造空间谱的计算表达式

k个谱峰所对应的位置即为估计出的doa，其中

以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的宗旨和范围，其均应涵盖在本发明的权利要求范围当中。本发明未详细描述的技术、形状、构造部分均为公知技术。