专利摘要

本发明公开了一种基于范德蒙分解的互质阵列波达方向估计方法，主要解决现有技术中自由度性能受限与计算复杂度高的问题，其实现步骤是：接收端天线按互质阵列结构进行架构；利用互质阵列接收入射信号并建模；计算互质阵列接收信号所对应的等价虚拟信号；构造虚拟阵列协方差矩阵；对虚拟阵列协方差矩阵进行去噪处理；将去噪后的虚拟阵列协方差矩阵进行具有唯一性的范德蒙分解；将范德蒙分解所得参量与理论值一一对应获得波达方向和功率估计结果。本发明充分利用了互质阵列能够增加自由度的优势和范德蒙分解方法无需预先设置网格点的特点，在降低计算复杂度的同时实现波达方向估计方法自由度的增加，可用于无源定位和目标探测。

权利要求

1.一种基于范德蒙分解的互质阵列波达方向估计方法，其特征在于，包含以下步骤：

(1)接收端使用2M+N-1个天线，并按照互质阵列结构进行架构；其中M与N为互质整数，且M<N；

(2)利用互质阵列接收K个方向为θ₁,θ₂,…,θ_K的远场窄带非相干信号源入射信号，则(2M+N-1)×1维互质阵列接收信号y(t)可建模为：

y(t)=Σk=1Kd(θk)sk(t)+n(t),

其中，s_k(t)为信号波形，n(t)为与各信号源相互独立的噪声分量，d(θ_k)为θ_k方向对应的互质阵列导引向量，可表示为

d(θk)=[1,e-jπu2sin(θk),...,e-jπu2M+N-1sin(θk)]T,

其中，u_i,i＝1,2,…,2M+N-1表示互质阵列中第i个物理天线阵元的实际位置，且u₁＝0，[·]^T表示转置操作。共采集T个采样快拍，得到采样协方差矩阵

R^=1TΣt=1Ty(t)yH(t),

这里(·)^H表示共轭转置；

(3)计算互质阵列接收信号所对应的等价虚拟信号。向量化互质阵列接收信号的采样协方差矩阵 获得虚拟阵列等价接收信号z：

z=vec(R^)=Vσ2+σn2i,

其中， 为(2M+N-1)²×K维矩阵， 包含K个入射信号源的功率， 为噪声功率，i＝vec(I_2M+N-1)。这里，vec(·)表示向量化操作，即把矩阵中的各列依次堆叠以形成一个新的向量，(·)^*表示共轭操作， 表示克罗内克积，I_2M+N-1表示(2M+N-1)×(2M+N-1)维单位矩阵。向量z对应的虚拟阵列中各虚拟阵元的位置为

去除集合 中各位置上重复的虚拟阵元，得到一个非均匀的虚拟阵列 其对应的等价虚拟信号 可通过选取向量z中相对应位置上的元素获得；

(4)构造虚拟阵列协方差矩阵。选取非均匀虚拟阵列 中以0为中心连续均匀排列的一段虚拟阵元，形成一个包含L个虚拟阵元的均匀虚拟阵列 其相应的虚拟阵元位置为(-L_v+1)d到(L_v-1)d之间的连续位置，其中，d为入射窄带信号波长的一半，且

Lv=L+12;

相应地，该均匀虚拟阵列的等价信号 可通过截取 中与该L个虚拟阵元所对应的位置上的元素获得，维度为L×1。虚拟阵列协方差矩阵R_v可通过空间平滑技术获得：将向量 分割为L_v个相互重叠的子向量，每个子向量的维度为L_v×1，包含向量 中的第i个至第i+L_v-1个元素，表示为

则R_v可通过取四阶统计量的主平方根获得：

其中，R_v的维度为L_v×L_v，对应于均匀虚拟阵列 中位置为0到(L_v-1)d的虚拟阵元，称为虚拟均匀线性阵列；R_v具有Toeplitz和Hermitian对称结构，在理论上可表示为：

其中， 为L_v×K维的虚拟均匀线性阵列导引矩阵，具有范德蒙结构。∑为对角元素为各信号功率 的对角矩阵， 表示L_v×L_v维单位矩阵；

(5)对虚拟阵列协方差矩阵R_v进行去噪处理，将其包含的噪声分量去除，获得去噪后的等价虚拟阵列协方差矩阵R_s；

(6)对去噪后的虚拟阵列协方差矩阵R_s进行具有唯一性的范德蒙分解。定义 其中，Λ_s为K×K维对角矩阵，包含将R_v特征值分解后的从大到小排列的前K个特征值，Ω_s为包含该K个特征值所对应特征向量的L_v×K维矩阵；同时，存在一个K×K维酉矩阵U满足X_＝X^-U，其中，X^-和X_分别为移除了X的最后一行和第一行得到的(L_v-1)×K维矩阵；对U进行特征值分解得到U＝QΛ_uQ^H，其中Λ_u为对角元素为[λ₁,λ₂,…,λ_K]的对角矩阵，Q为K×K维酉矩阵。定义 为波达方向估计值，λ_k可以表示为 同时定义p_k＝|x₁q_k|²，其中，x₁为矩阵X第一行形成的1×K维行向量，q_k为矩阵Q的第k列形成的K×1维列向量，去噪后的虚拟阵列协方差矩阵R_s可范德蒙分解为：

Rs=Σk=1Kpkdv(θ^k)dvH(θ^k)=Dv(θ^)PDvH(θ^),

其中， 为L_v×K维范德蒙矩阵，P为对角元素为[p₁,p₂,…,p_K]的对角矩阵；上述去噪后的虚拟阵列协方差矩阵R_s范德蒙分解唯一性的前提是K＜L_v；

(7)将范德蒙分解所得参量与理论值一一对应获得波达方向估计结果。步骤(6)中得到的范德蒙矩阵 的理论值为步骤(4)中虚拟均匀线性阵列导引矩阵定义中的参量 对角矩阵P的理论值为步骤(4)中的参量 信号的波达方向估计值为：

θ^k=arcsin(-1πimag(lnλk)), 2

其中，imag(·)表示复数的虚部。同时，相应波达方向的信号功率估计值为：

σ^k2=Lv|x1qk|2,

其中，未知参量λ_k和q_k可由下式获得：

(X^-HX_-λ_kX^-HX^-)q_k＝0，

即X^-HX_和X^-HX^-进行广义特征分解后得到的特征值和对应特征向量。

2.根据权利要求1所述的基于范德蒙分解的互质阵列波达方向估计方法，其特征在于：步骤(1)所述的互质阵列由一对稀疏均匀线性子阵列组合而成，其中第一个子阵列包含2M个天线阵元，阵元间距为Nd；第二个子阵列包含N个天线阵元，阵元间距为Md。将两个子阵列以首个天线阵元重叠的方式进行组合，得到包含2M+N-1个物理天线阵元的互质阵列架构。

3.根据权利要求1所述的基于范德蒙分解的互质阵列波达方向估计方法，其特征在于：步骤(5)所述的去噪后的虚拟阵列协方差矩阵R_s可由如下方法得到：将虚拟阵列协方差矩阵R_v进行特征值分解：

R_v＝Ω_sΛ_sΩ_s^H+Ω_nΛ_nΩ_n^H，

其中，Λ_n为(L_v-K)×(L_v-K)维对角矩阵，包含将R_v的特征值从大到小排列后的后(L_v-K)个特征值，Ω_n为包含该(L_v-K)个特征值所对应特征向量的L_v×(L_v-K)维矩阵。故R_v的信号子空间部分R_s＝Ω_sΛ_sΩ_s^H可看作去噪后的虚拟阵列协方差矩阵。

4.根据权利要求1所述的基于范德蒙分解的互质阵列波达方向估计方法，其特征在于：步骤(6)所述的虚拟阵列协方差矩阵R_s范德蒙分解的唯一性原理在于，当K<L_v时，任何K个不同的 均线性无关。

说明书

技术领域

本发明属于信号处理技术领域，尤其涉及对雷达信号、声学信号及电磁信号的波达方向估计，具体是一种基于范德蒙分解的互质阵列波达方向估计方法，可用于无源定位和目标探测。

背景技术

波达方向(Direction-of-Arrival,DOA)估计是阵列信号处理领域的一个重要分支，它是指利用阵列天线接收空域信号，并通过统计信号处理技术和各类优化方法对接收信号进行处理，以恢复信号中包含的DOA信息，在雷达、声呐、语音、无线通信等领域有着广泛的应用。

DOA估计方法的自由度是指其能够分辨的入射信号源个数，作为实际系统应用中的一个重要性能指标，决定着系统的总体复杂度。现有的DOA估计方法通常采用均匀线性阵列进行信号的接收与建模，但是基于均匀线性阵列方法的自由度受限于实际天线阵元个数。具体而言，对于一个包含L个天线阵元的均匀线性阵列，其自由度为L-1，即最多只能分辨L-1个入射信号。因此，当某个空域范围内入射信号源的个数大于或等于阵列中天线阵元的个数时，现有采用均匀线性阵列的方法将无法进行有效的DOA估计。为了增加自由度，传统方法需要通过增加物理天线阵元及相应的射频模块来实现，这造成了系统计算复杂度和硬件复杂度的增加。因此，现有采用均匀线性阵列的DOA估计方法在自由度性能与计算复杂度之间存在着一定的利弊权衡问题。

与均匀阵列相比，互质阵列能够采用相同个数的天线阵元取得更多的自由度，因而受到了广泛关注。作为互质采样技术在空间域上的一个典型表现形式，互质阵列提供了一个系统化的稀疏阵列架构方案，并能够突破传统均匀线性阵列自由度受限的瓶颈，实现DOA估计方法自由度性能的提升。现有的基于互质阵列的DOA估计方法可以通过利用质数的性质将互质阵列推导到虚拟域，并形成等价虚拟均匀线性阵列接收信号以实现DOA估计。由于虚拟阵列中包含的虚拟阵元数大于实际的天线阵元数，自由度因此得到了有效的提升。但是，现有基于等价虚拟信号的方法通常需要通过设置预先定义的网格点来设计优化问题，从而进行DOA估计。由于实际中信号的来波方向不会完全位于这些预先定义的网格点上，导致了DOA估计的精度受限于网格密度。尽管可以通过减小网格间的采样间隔实现估计精度的提升，网格点的密集化将导致方法计算复杂度的指数型增长。此外，现有的DOA估计方法的频谱图通常为伪谱，其DOA估计方向的波峰响应无法用来表示信号功率；而事实上，信号功率也是描述信号源的一个重要参量，我们希望频谱图能够在估计波达方向的同时反映其功率信息。

发明内容

本发明的目的在于针对上述现有技术存在的不足，提出一种基于范德蒙分解的互质阵列波达方向估计方法，通过利用互质阵列的特性推导虚拟域等价接收信号，以提升估计方法的自由度，并对虚拟阵列协方差矩阵进行具有唯一性的范德蒙分解，在进行无网格化波达方向估计的同时实现了相应波达方向上的功率估计，从而在实际应用过程中降低了计算复杂度和硬件复杂度，提高了系统的整体效率和估计精确度。

本发明的目的是通过以下技术方案来实现的：一种基于范德蒙分解的互质阵列波达方向估计方法，包含以下步骤：

(1)接收端使用2M+N-1个天线，并按照互质阵列结构进行架构；其中M与N为互质整数，且M<N；

(2)利用互质阵列接收K个方向为θ₁,θ₂,…,θ_K的远场窄带非相干信号源入射信号，则(2M+N-1)×1维互质阵列接收信号y(t)可建模为：

其中，s_k(t)为信号波形，n(t)为与各信号源相互独立的噪声分量，d(θ_k)为θ_k方向对应的互质阵列导引向量，可表示为

其中，u_i,i＝1,2,,2M+N-1表示互质阵列中第i个物理天线阵元的实际位置，且u₁＝0，[·]^T表示转置操作。共采集T个采样快拍，得到采样协方差矩阵

这里(·)^H表示共轭转置；

(3)计算互质阵列接收信号所对应的等价虚拟信号。向量化互质阵列接收信号的采样协方差矩阵 获得虚拟阵列等价接收信号z：

其中， 为(2M+N-1)²×K维矩阵， 包含K个入射信号源的功率， 为噪声功率，i＝vec(I_2M+N-1)。这里，vec(·)表示向量化操作，即把矩阵中的各列依次堆叠以形成一个新的向量，(·)^*表示共轭操作， 表示克罗内克积，I_2M+N-1表示(2M+N-1)×(2M+N-1)维单位矩阵。向量z对应的虚拟阵列中各虚拟阵元的位置为

去除集合 中各位置上重复的虚拟阵元，得到一个非均匀的虚拟阵列 其对应的等价虚拟信号 可通过选取向量z中相对应位置上的元素获得；

(4)构造虚拟阵列协方差矩阵。选取非均匀虚拟阵列 中以0为中心连续均匀排列的一段虚拟阵元，形成一个包含L个虚拟阵元的均匀虚拟阵列 其相应的虚拟阵元位置为(-L_v+1)d到(L_v-1)d之间的连续位置，其中，d为入射窄带信号波长的一半，且

相应地，该均匀虚拟阵列的等价信号 可通过截取 中与该L个虚拟阵元所对应的位置上的元素获得，维度为L×1。虚拟阵列协方差矩阵R_v可通过空间平滑技术获得：将向量 分割为L_v个相互重叠的子向量，每个子向量的维度为L_v×1，包含向量 中的第i个至第i+L_v-1个元素，表示为

则R_v可通过取四阶统计量的主平方根获得：

其中，R_v的维度为L_v×L_v，对应于均匀虚拟阵列 中位置为0到(L_v-1)d的虚拟阵元，称为虚拟均匀线性阵列；R_v具有Toeplitz和Hermitian对称结构，在理论上可表示为：

其中， 为L_v×K维的虚拟均匀线性阵列导引矩阵，具有范德蒙结构。∑为对角元素为各信号功率 的对角矩阵， 表示L_v×L_v维单位矩阵；

(5)对虚拟阵列协方差矩阵R_v进行去噪处理，将其包含的噪声分量去除，获得去噪后的等价虚拟阵列协方差矩阵R_s；

(6)对去噪后的虚拟阵列协方差矩阵R_s进行具有唯一性的范德蒙分解。定义 其中，Λ_s为K×K维对角矩阵，包含将R_v特征值分解后的从大到小排列的前K个特征值，Ω_s为包含该K个特征值所对应特征向量的L_v×K维矩阵；同时，存在一个K×K维酉矩阵U满足X_-＝X^-U，其中，X^-和X_-分别为移除了X的最后一行和第一行得到的(L_v-1)×K维矩阵；对U进行特征值分解得到U＝QΛ_uQ^H，其中Λ_u为对角元素为[λ₁,λ₂,…,λ_K]的对角矩阵，Q为K×K维酉矩阵。定义 为波达方向估计值，λ_k可以表示为 同时定义p_k＝|x₁q_k|²，其中，x₁为矩阵X第一行形成的1×K维行向量，q_k为矩阵Q的第k列形成的K×1维列向量，去噪后的虚拟阵列协方差矩阵R_s可范德蒙分解为：

其中， 为L_v×K维范德蒙矩阵，P为对角元素为[p₁,p₂,…,p_K]的对角矩阵；上述去噪后的虚拟阵列协方差矩阵R_s范德蒙分解唯一性的前提是K＜L_v；

(7)将范德蒙分解所得参量与理论值一一对应获得波达方向估计结果。步骤(6)中得到的范德蒙矩阵 的理论值为步骤(4)中虚拟均匀线性阵列导引矩阵定义中的参量 对角矩阵P的理论值为步骤(4)中的参量 信号的波达方向估计值为：

其中，imag(·)表示复数的虚部。同时，相应波达方向的信号功率估计值为：

其中，未知参量λ_k和q_k可由下式获得：

(X^-HX_{_}-λ_kX^-HX^-)q_k＝0，

即X^-HX_{_}和X^-HX^-进行广义特征分解后得到的特征值和对应特征向量。

进一步地，步骤(1)所述的互质阵列由一对稀疏均匀线性子阵列组合而成，其中第一个子阵列包含2M个天线阵元，阵元间距为Nd；第二个子阵列包含N个天线阵元，阵元间距为Md。将两个子阵列以首个天线阵元重叠的方式进行组合，得到包含2M+N-1个物理天线阵元的互质阵列架构。

进一步地，步骤(5)所述的去噪后的虚拟阵列协方差矩阵R_s可由如下方法得到：将虚拟阵列协方差矩阵R_v进行特征值分解：

R_v＝Ω_sΛ_sΩ_s^H+Ω_nΛ_nΩ_n^H，

其中，Λ_n为(L_v-K)×(L_v-K)维对角矩阵，包含将R_v的特征值从大到小排列后的后(L_v-K)个特征值，Ω_n为包含该(L_v-K)个特征值所对应特征向量的L_v×(L_v-K)维矩阵。故R_v的信号子空间部分R_s＝Ω_sΛ_sΩ_s^H可看作去噪后的虚拟阵列协方差矩阵。

进一步地，步骤(6)所述的虚拟阵列协方差矩阵R_s范德蒙分解的唯一性原理在于，当K<L_v时，任何K个不同的 均线性无关。

本发明与现有技术相比具有以下优点：

(1)本发明充分利用了互质阵列能够增加DOA估计的自由度这一优势，将获得的信号模型推导至虚拟域，并通过虚拟阵列协方差矩阵实现DOA估计，可分辨的入射信号源个数大于物理天线阵元个数，使得自由度得到提升，并降低了计算复杂度和硬件复杂度；

(2)本发明利用具有唯一性的范德蒙分解的方法，无需预先设置网格点来进行波达方向估计，避免了传统方法预先定义网格点所导致的固有估计误差；

(3)本发明在有效估计入射信号DOA的同时还能实现相应DOA的功率估计，所得的空间谱能够同时反映入射信号的波达方向信息和功率信息。

附图说明

图1是本发明的总体流程框图。

图2是本发明中组成互质阵列的一对稀疏均匀子阵列结构示意图。

图3是本发明中互质阵列的结构示意图。

图4是本发明所提方法与现有各方法的空间谱对比示意图；其中图4(a)是基于信号稀疏重建的互质阵列波达方向估计方法；图4(b)是基于多重信号子空间分类的互质阵列波达方向估计方法；图4(c)是基于虚拟阵列内插的互质阵列波达方向估计方法；图4(d)是本发明所提方法。

图5是本发明所提方法与现有各方法均方根误差和信噪比之间的关系曲线图。

图6是本发明所提方法与现有各方法均方根误差和采样快拍数之间的关系曲线图。

具体实施方式

以下参照附图，对本发明的技术方案和效果作进一步的详细说明。

对于DOA估计方法在实际系统中的应用，估计精度、自由度和计算复杂度是重要的技术性能指标。现有方法在自由度性能上受限于物理天线阵元个数，而计算复杂度和估计精度也大大受限于预先定义的网格点的设置，且获得的空间谱响应无法反映功率信息。为了能在增加自由度的同时进行无网格化的DOA估计，本发明提供了一种基于范德蒙分解的互质阵列波达方向估计方法，参照图1，本发明的实现步骤如下：

步骤一：在接收端使用2M+N-1个物理天线阵元架构互质阵列。首先选取一对互质整数M、N，且M<N；然后，参照图2，构造一对稀疏均匀线性子阵列，其中第一个子阵列包含2M个间距为Nd的天线阵元，其位置为0,Nd,…,(2M-1)Nd，第二个子阵列包含N个间距为Md的天线阵元，其位置为0,Md,…,(N-1)Md；单位间隔d取为入射窄带信号波长的一半；接着，将两个子阵列按照首个阵元重叠的方式进行子阵列组合，参照图3，获得实际包含2M+N-1个天线阵元的非均匀互质阵列。

步骤二：采用互质阵列接收信号并建模。假设有K个来自θ₁,θ₂,…,θ_K方向的远场窄带非相干信号源，采用步骤一架构的非均匀互质阵列接收入射信号，得到(2M+N-1)×1维互质阵列接收信号y(t)，可建模为：

其中，s_k(t)为信号波形，n(t)为与各信号源相互独立的噪声分量，d(θ_k)为θ_k方向的导引向量，表示为

其中，u_i,i＝1,2,…,2M+N-1表示互质阵列中第i个物理天线阵元的实际位置，且u₁＝0，[·]^T表示转置操作。采集T个采样快拍，得到采样协方差矩阵

这里(·)^H表示共轭转置。

步骤三：计算互质阵列接收信号所对应的等价虚拟信号。向量化互质阵列接收信号的采样协方差矩阵 获得虚拟阵列等价接收信号z：

其中， 为(2M+N-1)²×K维虚拟阵列导引矩阵， 包含K个入射信号源的功率， 为噪声功率，i＝vec(I_2M+N-1)。这里，vec(·)表示向量化操作，即把矩阵中的各列依次堆叠以形成一个新的向量，(·)^*表示共轭操作， 表示克罗内克积，i_2M+N-1表示(2M+N-1)×(2M+N-1)维单位矩阵。向量z对应的虚拟阵列中各虚拟阵元的位置为

去除集合 中各位置上重复的虚拟阵元，得到一个非均匀的虚拟阵列 其对应的等价虚拟信号 可通过选取向量z中相对应位置上的元素获得。

步骤四：构造虚拟阵列协方差矩阵。首先，选取非均匀虚拟阵列 中以0为中心连续均匀排列的一段虚拟阵元，形成一个包含L个虚拟阵元的均匀虚拟阵列 (由于 中的虚拟阵元以零位对称分布，L始终为奇数)，其相应的虚拟阵元位置为(-L_v+1)d到(L_v-1)d之间的连续位置，其中

相应地，该均匀虚拟阵列的等价信号 可通过截取 中与该L个虚拟阵元所对应的位置上的元素获得，维度为L×1。接着，虚拟阵列协方差矩阵R_v可通过空间平滑技术获得：将向量 分割为L_v个维度为L_v×1的子向量，每个子向量包含向量 中的第i个至第i+L_v-1个元素，即：

则R_v可通过取四阶统计量的主平方根获得：

其中，R_v的维度为L_v×L_v，对应于均匀虚拟阵列 中位置为0到(L_v-1)d的虚拟阵元，称为虚拟均匀线性阵列，因此R_v具有Toeplitz和Hermitian对称结构，且R_v在理论上可以等价地表示为：

其中， 为L_v×K维的虚拟均匀线性阵列导引矩阵，具有范德蒙结构。∑为包含各信号功率 的对角矩阵， 表示L_v×L_v维单位矩阵。

步骤五：将虚拟阵列协方差矩阵R_v进行去噪处理，得到去噪后的虚拟协方差矩阵R_s。由于R_v中的噪声分量将会影响波达方向估计的准确性，我们首先将获得的虚拟阵列协方差矩阵R_v中的噪声分量去除。将R_v进行特征值分解为

R_v＝Ω_sΛ_sΩ_s^H+Ω_nΛ_nΩ_n^H，

其中，Λ_s为K×K维对角矩阵，包含将R_v的特征值从大到小排列后的前K个特征值，Ω_s为包含该K个特征值所对应特征向量的L_v×K维矩阵；相应地，Λ_n为(L_v-K)×(L_v-K)维对角矩阵，包含余下的(L_v-K)个R_v的特征值，Ω_n为包含该(L_v-K)个特征值所对应特征向量的L_v×(L_v-K)维矩阵。于是，R_v的信号子空间部分R_s＝Ω_sΛ_sΩ_s^H可看作是去噪后的虚拟阵列协方差矩阵。

步骤六：对去噪后的虚拟阵列协方差矩阵R_s进行具有唯一性的范德蒙分解。首先，由R_s的半正定性，根据步骤五定义 使得R_s＝XX^H。其次，R_s的Toeplitz结构使得移除了X的最后一行和第一行得到的矩阵X^-和X_{_}满足等式X^-X^-H＝X_{_}X_{_}^H，这说明了存在一个K×K维酉矩阵U满足X_{_}＝X^-U。进一步对U进行特征值分解后得到，U＝QΛ_uQ^H，其中Λ_u包含对角元素为[λ₁,λ₂,…,λ_K]的对角矩阵，Q为K×K维酉矩阵。利用酉矩阵U的各特征值的模均为1性质，在引入一个新参数 之后，λ_k可以表示为 同时定义p_k＝|x₁q_k|²，其中，x₁为矩阵X第一行形成的1×K维行向量，q_k为矩阵Q的第k列形成的K×1维列向量，将去噪后的虚拟阵列协方差矩阵R_s进行范德蒙分解，得到如下形式：

其中， 是一个典型的L_v×K维范德蒙矩阵，P为对角元素为[p₁,p₂,…,p_K]的对角矩阵。当K<L_v时，任何K个不同的 均线性无关，则该范德蒙分解是唯一的。

步骤七：将范德蒙分解所得参量与理论值一一对应获得波达方向估计结果。将步骤六中R_s的范德蒙分解结果与步骤四中R_v的信号子空间部分理论值相对比，由于范德蒙分解具有唯一性，步骤六中得到的范德蒙矩阵 等价于步骤四中的虚拟阵列导引矩阵理论值 而步骤六中得到的对角矩阵P等价于步骤四中的 因此，步骤六中的 和 可以作为信号波达方向和功率的估计值。故我们可以得到信号的波达方向估计：

其中，imag(·)表示复数的虚部。同时可以得到对应的信号功率估计：

这里，未知参量λ_k和q_k可由如下方法得到：将步骤六中的U＝QΛ_uQ^H代入X_-＝X^-U后得到，

(X^-HX_{_}-λ_kX^-HX^-)q_k＝0，

故对X^-HX_{_}和X^-HX^-进行广义特征分解后得到的特征值和对应的特征向量即为λ_k和q_k。

本发明一方面充分利用了互质阵列能够增加DOA估计方法自由度的优势，突破了均匀线性阵列自由度受限的瓶颈，实现了在天线阵元个数一定的条件下估计更多个数的入射信号源；另一方面利用了具有唯一性的范德蒙分解方法，使DOA估计能在无网格的条件下进行，大大降低了计算复杂度，并在进行波达方向估计的同时频谱图能够直观而准确地反映功率信息。

下面结合仿真实例对本发明所提方法的效果做进一步的描述。

仿真条件：互质阵列的参数选取为M＝3，N＝5，即架构的互质阵列共包含2M+N-1＝10个天线阵元。波达方向角的角度域范围为[-90°,90°]，空间域网格点均匀采样间距设置为0.1°；正则化参数ξ设置为0.25。

仿真实例1：本发明所提出的基于范德蒙分解的互质阵列波达方向估计方法将与现有的基于信号稀疏重建的互质阵列波达方向估计方法，基于多重信号子空间分类的互质阵列波达方向估计方法和基于虚拟阵列内插的互质阵列波达方向估计方法进行空间谱图的比较，以说明本发明方法在功率谱特征和功率估计方面的优势。由于仿真条件中入射信号源的个数大于天线阵元的个数，采用均匀线性阵列的方法自由度受限而无法实现有效的DOA估计，因而在此不作比较。假定入射窄带信号个数为15，且入射方向均匀分布于-60°至60°这一区间范围内，信噪比为10dB，采样快拍数为T＝500。各方法的空间谱如图4所示，其中，虚线表示真实的入射信号波达方向。可以看出，由于所有方法均基于互质阵列，故均能仅用10个物理阵元对15个入射信号进行有效的波达方向估计，体现了互质阵列在自由度上的优势。但是，如图4(a)所示，基于信号稀疏重建的互质阵列波达方向估计方法的信号波峰响应周围存在相当数量的虚峰，这将导致估计误差。而如图4(b)所示的基于多重信号子空间分类的互质阵列波达方向估计方法的信号波峰响应并不完全对应于真实的入射信号波达方向。同时，基于多重信号子空间分类的互质阵列波达方向估计方法和基于虚拟阵列内插的互质阵列波达方向估计方法的空间谱均属于伪谱，信号的波峰响应并不能反映功率信息。相比之下，本发明所提方法能实现所有信号源波达方向的有效估计，且空间谱的波峰响应为相对应波达方向上的信号功率。

仿真实例2：本发明所提出的基于范德蒙分解的互质阵列波达方向估计方法将与现有的基于信号稀疏重建的互质阵列波达方向估计方法，基于多重信号子空间分类的互质阵列波达方向估计方法和基于虚拟阵列内插的互质阵列波达方向估计方法进行均方根误差的对比。每一次蒙特卡罗试验中的入射信号方向均由随机标准正态分布产生，均方根误差(root mean square error,RMSE)的计算公式为

其中， 和θ_q表示每一次蒙特卡罗试验中的波达方向估计值和真实值。对于每一组参数取值，蒙特卡罗试验的次数为Q＝500。图6为均方根误差和信噪比之间的关系曲线图，采样快拍数设置为T＝50；可以看出，本发明所提方法的RMSE性能明显优于其他方法，尤其当信噪比大于10dB时，基于范德蒙分解的特点保证了本方法无网格估计的优势，使得RMSE随信噪比的增大而持续减小。相较而言，基于多重信号子空间分类的互质阵列波达方向估计方法和基于信号稀疏重建的互质阵列波达方向估计方法都是需要预先设置网格点来进行DOA估计，所以它们的RMSE性能受限于网格点的采样间距，可以明显地看到，在信噪比大于10dB之后以上两种方法的RMSE始终维持在0.1左右。而基于虚拟阵列内插的互质阵列波达方向估计方法在该采样快拍数相对较少的情况下，并没有体现无网格化波达方向估计方法的优势。图6为均方根误差和采样快拍数之间的关系曲线图，信噪比设置为0dB；与图5相似，本发明所提方法的RMSE随采样快拍数的增大而持续减小，在与其他方法的比较中表现出了最好的性能。

综上所述，本发明主要解决了现有技术在DOA估计自由度性能与计算复杂度方面存在的不足，一方面充分利用互质阵列的特性在虚拟域进行信号处理以实现自由度的增加；另一方面基于范德蒙分解的无网格化DOA估计有效地降低了计算复杂度，且能够同时估计各相应信号源的功率，在无源定位和目标探测等实际应用中表现出突出的优势。

基于范德蒙分解的互质阵列波达方向估计方法专利购买费用说明

Q：办理专利转让的流程及所需资料

A：专利权人变更需要办理著录项目变更手续，有代理机构的，变更手续应当由代理机构办理。

1：专利变更应当使用专利局统一制作的“著录项目变更申报书”提出。

2：按规定缴纳著录项目变更手续费。

3：同时提交相关证明文件原件。

4：专利权转移的，变更后的专利权人委托新专利代理机构的，应当提交变更后的全体专利申请人签字或者盖章的委托书。

Q:专利著录项目变更费用如何缴交

A：（1）直接到国家知识产权局受理大厅收费窗口缴纳,（2）通过代办处缴纳，（3）通过邮局或者银行汇款,更多缴纳方式

Q：专利转让变更，多久能出结果

A：著录项目变更请求书递交后，一般1-2个月左右就会收到通知，国家知识产权局会下达《转让手续合格通知书》。