本发明属于信号处理技术领域，尤其涉及对雷达信号、声学信号及电磁信号的波达方向估计，具体是一种基于多采样快拍互质阵列接收信号快速傅里叶逆变换的波达方向估计方法，可用于无源定位和目标探测。
背景技术：
波达方向(direction-of-arrival,doa)估计是阵列信号处理领域的一个基本问题，它是指利用传感器阵列接收信号，并通过一系列信号处理方法对接收信号进行统计处理，从而获得信号中所包含的波达方向信息，在雷达、声呐、语音、无线通信等领域均有着广泛的应用。
波达方向估计方法的分辨率是指其对入射角度相近的不同信号的分辨能力，通常由阵列孔径决定，阵列孔径越大则对应的分辨率越高。通常地，分辨率越高对应了波达方向估计空间谱中峰值宽度越窄。目前波达方向估计方法大多基于均匀线性阵列结构，其分辨率性能的提升主要通过增加更多的物理阵元以获得更大的阵列孔径实现的，而更多的物理阵元也引入了较高的系统软硬件复杂度。
互质阵列作为一种具有系统化结构的稀疏阵列，近年来受到了学术界的广泛关注，现有基于互质阵列的doa估计方法大多在互质阵列二阶等价虚拟信号的基础上进行较为复杂的运算，包括求逆、特征值分解等复杂的矩阵运算，以及凸优化问题的设计与求解等过程，这些运算过程导致了较高的计算复杂度，在实时性要求较高的应用场景下面临一定的挑战，且在实际系统中的硬件实现较为困难。
技术实现要素：
本发明的目的在于针对上述现有技术存在的不足，提出一种基于多采样快拍互质阵列接收信号快速傅里叶逆变换的波达方向估计方法，在互质阵列一阶多采样快拍接收信号的基础上进行快速傅里叶逆变换，保证波达方向估计结果正确性的同时可以获得比传统均匀线性阵列更高的分辨率，且有效地降低了系统软硬件复杂度。
本发明的目的是通过以下技术方案来实现的：一种基于多采样快拍互质阵列接收信号快速傅里叶逆变换的波达方向估计方法，包含以下步骤：
(1)接收端使用2m+n-1个阵元，并按照互质阵列结构进行架构；其中m与n为互质整数；
(2)假设有l个来自θ＝[θ1,θ2,…,θl]t方向的远场窄带非相干信号源，[·]t表示转置操作，利用步骤(1)中构建的互质阵列对入射信号进行接收，则在t时刻的互质阵列接收信号x(t)可建模为：
其中，x(t)为(2m+n-1)×1维向量，sl(t)为第l个入射信号的波形，n(t)为与各信号源相互独立的噪声分量，a(θl)为对应于θl方向信号源的互质阵列导引向量，可表示为
其中，μi,i＝1,2,3,…,2m+n-1表示互质阵列中第i个物理阵元的实际位置，且首个物理阵元的位置为μ1＝0，λ为入射窄带信号的波长，j为虚数单位；
(3)构造互质阵列一阶多采样快拍接收信号：采用连续t个单采样快拍接收信号x(t)的一阶统计量作为互质阵列一阶多采样快拍接收信号
(4)对互质阵列一阶多采样快拍接收信号进行补零操作：保持互质阵列中物理阵元的位置不变，向互质阵列一阶多采样快拍接收信号
(5)对补零后的互质阵列一阶多采样快拍接收信号
其中，k＝0,1,…,k-1，arccos(·)为反余弦函数，h为保证
(6)根据空间谱进行波达方向估计：对步骤(5)构建的空间谱进行谱峰搜索操作，将幅度最大的前l个峰值对应的角度，作为l个入射信号的波达方向估计。
进一步地，步骤(1)所述的互质阵列结构可具体描述为：首先，选取一组互质的整数m、n，构造一对稀疏均匀线性子阵列。第一个子阵列包含2m个间距为nd的阵元，其位置为0,nd,…,(2m-1)nd；第二个子阵列包含n个间距为md的阵元，其位置为0,md,…,(n-1)md，其中单位间距d为入射窄带信号的半波长即d＝λ/2。之后以两个子阵列的第一个阵元为参考阵元，将两个参考阵元重合使得两个子阵列组合为一个包含2m+n-1个实际阵元的非均匀互质阵列。
进一步地，所述步骤(3)中，一阶统计量采用连续t个单采样快拍接收信号x(t)的平均值或者加和形式。
进一步地，所述步骤(5)中，通过快速傅里叶逆变换得到接收信号
其中，
所得的空间响应
进一步地，所述步骤(5)中，所构建的空间谱反映了空间中各角度上的响应幅度，其中存在对应于l个入射信号的l个峰值。
本发明与现有技术相比具有以下优点：
(1)本发明所提方法在互质阵列一阶多采样快拍接收信号的基础上进行快速傅里叶逆变换，从而获得空间响应并基于此构建一个空间谱，通过对所构建空间谱的谱峰搜索过程获得波达方向估计结果，避免了传统波达方向估计方法中常用的凸优化问题设计求解以及矩阵求逆、矩阵特征值分解等复杂计算过程，有效地降低了计算复杂度，更好地满足了对实时性有较高要求的实际应用场景，同时易于在实际系统中硬件实现。
(2)利用互质阵列架构，在相同物理阵元数量的情况下，获得比传统均匀线性阵列更大的阵列孔径，在保证波达方向估计正确性的基础上，实现了更高的分辨率。
附图说明
图1是本发明的方法总体流程框图；
图2是本发明中组成互质阵列的一对稀疏均匀子阵列结构示意图；
图3是本发明中互质阵列的结构示意图；
图4是用于体现本发明所提方法分辨率性能的空间谱示意图。
具体实施方式
以下参照附图，对本发明的技术方案和效果作进一步的详细说明。
现有波达方向估计方法大多基于均匀线性阵列结构，主要通过增加更多的物理阵元增大阵列孔径，从而获得更高的分辨率，这同时增加了系统软硬件复杂度。针对以上问题，本发明提供了一种基于多采样快拍互质阵列接收信号快速傅里叶逆变换的波达方向估计方法，参照图1，本发明的实现步骤如下：
步骤一：在接收端使用2m+n-1个实际阵元构建互质阵列。首先，选取一组互质的整数m、n，构造一对稀疏均匀线性子阵列。第一个子阵列包含2m个间距为nd的阵元，其位置为0,nd,…,(2m-1)nd；第二个子阵列包含n个间距为md的阵元，其位置为0,md,…,(n-1)md，其中单位间距d为入射窄带信号的半波长即d＝λ/2。之后将两个子阵列按照首个阵元重合的方式进行组合，得到包含2m+n-1个实际阵元的非均匀互质阵列。
步骤二：采用互质阵列接收信号并对信号建模。假设有l个来自θ＝[θ1,θ2,…,θl]t方向的远场窄带非相干信号源，[·]t表示转置操作，利用步骤一中构建的互质阵列对入射信号进行接收，则在t时刻的互质阵列接收信号x(t)可建模为：
其中，x(t)为(2m+n-1)×1维向量，sl(t)为第l个入射信号的波形，n(t)为与各信号源相互独立的噪声分量，a(θl)为对应于θl方向信号源的互质阵列导引向量，可表示为
其中，μi,i＝1,2,3,…,2m＋n-1表示互质阵列中第i个物理阵元的实际位置，且首个物理阵元的位置为μ1＝0，λ为入射窄带信号的波长，j为虚数单位。
步骤三：构造互质阵列一阶多采样快拍接收信号。采用连续t个单采样快拍接收信号x(t)的平均值作为互质阵列一阶多采样快拍接收信号
步骤四：对互质阵列一阶多采样快拍接收信号进行补零操作。保持互质阵列中物理阵元的位置不变，向互质阵列一阶多采样快拍接收信号
其中，0表示满足要求的适当长度的零向量。
步骤五：对补零后的互质阵列一阶多采样快拍接收信号
其中，
所得的空间响应
其中，k＝0,1,…,k-1，arccos(·)为反余弦函数，h为保证
其中，[·]k表示向量中第k个元素，|·|表示复数的模。所构建的空间谱反映了空间中各角度上的响应幅度，其中存在对应于l个入射信号的l个峰值。
步骤六：根据所构建的空间谱进行波达方向估计。对步骤(5)中空间谱进行谱峰搜索操作，将其峰值按照幅度由高到低排列，则最大的前l个峰值对应的角度为l个入射信号的波达方向估计。
本发明所提方法在互质阵列一阶多采样快拍接收信号的基础上进行快速傅里叶逆变换，从而获得空间响应并基于此构建一个空间谱，通过对所构建空间谱的谱峰搜索过程获得波达方向估计结果，在保证波达方向估计结果正确性的同时，获得比传统均匀线性阵列更高的分辨率，且快速傅里叶逆变换的计算复杂度仅为o(klogk)，更好地满足了对估计实时性有较高要求的应用场景。进一步地，由于快速傅里叶逆变换的计算仅由复数加法和乘法组成，本发明所提波达方向估计方法相比于传统基于互质阵列的波达方向估计方法而言更容易在实际硬件系统中实现。
下面结合仿真实例对本发明的效果做进一步的描述。
采用互质阵列接收入射信号，其参数选取为m＝9，n＝10，即架构的互质阵列共包含2m+n-1＝27个物理阵元。固定采样快拍数t＝500，假定空间中有2个入射窄带信号，入射方向分别为45°和50°，补零之后的快速傅里叶逆变换点数k＝2048。本发明所提方法得到的空间谱如图4所示，其中垂直虚线代表入射信号源的真实方向。
综上所述，本发明所提方法通过在互质阵列一阶多采样快拍接收信号的基础上进行快速傅里叶逆变换得到其空间响应并基于此构建一个空间谱，通过对所构建空间谱的谱峰搜索过程得到doa估计，在保证波达方向估计结果正确性的同时，获得比传统均匀线性阵列更高的分辨率，同时有效地降低了计算复杂度。
上述实施例用来解释说明本发明，而不是对本发明进行限制，在本发明的精神和权利要求的保护范围内，对本发明作出的任何修改和改变，都落入本发明的保护范围。