基于互质阵列的稀疏频率波形产生方法、装置及存储介质

技术领域

本发明涉及通信传输领域，特别涉及一种基于互质阵列的稀疏频率波形产生方法、装置及存储介质。

背景技术

互质阵列是一种结构稀疏阵列，由一对阵元数满足互质条件的稀疏均匀线性阵列构成。互质阵列相较于传统均匀阵列，能够实现入射信号的欠采样，从而突破奈奎斯特定理对天线阵元间距的限制；可以在相同阵元的条件下比常规的均匀线性阵列拥有更多的孔径，实现更高的分辨率。

现今，由雷达、通信系统等电磁波收发设备的增加而导致的频谱拥挤情况已经不容忽视。对于阵列而言，采用传统波形的雷达系统与同频工作的相邻阵元可能会产生互相干扰的情况。

发明内容

本发明旨在至少解决现有技术中存在的技术问题之一。

为此，本发明提出一种基于互质阵列的稀疏频率波形产生方法，能够很好地减少波形之间相互干扰的情况，提升检测效果。

本发明还提出一种应用上述基于互质阵列的稀疏频率波形产生方法的基于互质阵列的稀疏频率波形产生装置。

本发明还提出一种应用上述基于互质阵列的稀疏频率波形产生方法的计算机可读存储介质。

根据本发明第一方面实施例的基于互质阵列的稀疏频率波形产生方法，包括：基于互质阵列获取期望方向图；

通过所述期望方向图获取自相关矩阵；

通过所述自相关矩阵并且结合恒模约束和稀疏频率约束，获取结果波形。

根据本发明实施例的基于互质阵列的稀疏频率波形产生方法，至少具有如下有益效果：首先利用互质阵列结构的特点，增加虚拟孔径以获取期望方向图；然后针对期望方向图，获取自相关矩阵，使得发射方向图接近期望方向图；然后针对期望方向图，结合恒模约束和稀疏频率约束，最后获得结果波形，能够很好地减少波形之间相互干扰的情况，提升检测效果。

根据本发明的一些实施例，所述基于互质阵列获取期望方向图，包括：基于互质阵列结构，得到增广虚拟阵列和协方差矩阵；

通过所述增广虚拟阵列和所述协方差矩阵，得到空间平滑虚拟域协方差矩阵；

通过所述空间平滑虚拟域协方差矩阵，利用MUSIC算法计算得到期望方向图。

根据本发明的一些实施例，所述通过所述期望方向图获取自相关矩阵，包括：通过半正定约束由所述期望方向图得到自相关矩阵。

根据本发明的一些实施例，所述通过所述自相关矩阵并且结合恒模约束和稀疏频率约束，获取结果波形，包括：

通过自相关矩阵并且结合恒模约束，获得雷达波形矩阵最优解；

通过所述雷达波形矩阵最优解并且结合稀疏频率约束，构建恒模波形目标函数；通过所述恒模波形目标函数结合循环算法，获得结果波形。

根据本发明的一些实施例，所述空间平滑虚拟域协方差矩阵可表示为：

其中，MN+1为增广虚拟阵列分割为MN+1个分别包含MN+1个虚拟阵元的线性子阵列，R

根据本发明的一些实施例，所述雷达波形矩阵最优解可表示为：

其中，X表示N×M雷达波形矩阵，R表示自相关矩阵，V表示变换矩阵，N表示发射信号在一个脉冲重复周期的采样数，V

根据本发明的一些实施例，所述自相关矩阵的目标函数可表示为

s.t r

R≥0

其中，优化变量α是与方向图幅度有关的比例值，r

根据本发明第二方面实施例的基于互质阵列的稀疏频率波形产生装置，包括：处理单元，用于基于互质阵列获取期望方向图；

运算单元，用于通过所述期望方向图获取自相关矩阵；

计算单元，用于通过所述自相关矩阵并且结合恒模约束和稀疏频率约束，获取结果波形。

根据本发明的一些实施例，所述处理单元包括：

第一处理单元，用于基于互质阵列结构，得到增广虚拟阵列和协方差矩阵；

第二处理单元，用于通过所述增广虚拟阵列和所述协方差矩阵，得到空间平滑虚拟域协方差矩阵；

第三处理单元，用于通过所述空间平滑虚拟域协方差矩阵，利用MUSIC算法计算得到期望方向图。

根据本发明实施例的基于互质阵列的稀疏频率波形产生装置，至少具有如下有益效果：首先利用互质阵列结构的特点，增加虚拟孔径以获取期望方向图；然后针对期望方向图，获取自相关矩阵，使得发射方向图接近期望方向图；然后针对期望方向图，结合恒模约束和稀疏频率约束，最后获得结果波形，能够很好地减少波形之间相互干扰的情况，提升检测效果。

根据本发明第三方面实施例的计算机可读存储介质，能够应用根据本发明上述第一方面实施例的基于互质阵列的稀疏频率波形产生方法。

根据本发明实施例的计算机可读存储介质，至少具有如下有益效果：首先利用互质阵列结构的特点，增加虚拟孔径以获取期望方向图；然后针对期望方向图，获取自相关矩阵，使得发射方向图接近期望方向图；然后针对期望方向图，结合恒模约束和稀疏频率约束，最后获得结果波形，能够很好地减少波形之间相互干扰的情况，提升检测效果。

本发明的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。

附图说明

本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：

图1为本发明实施例一的基于互质阵列的稀疏频率波形产生方法的工作流程图；

图2为本发明实施例一的基于互质阵列的稀疏频率波形产生方法中的基于互质阵列获取期望方向图的工作流程图；

图3为本发明实施例一的基于互质阵列的稀疏频率波形产生方法中的通过自相关矩阵获取结果波形的工作流程图；

图4为本发明实施例二的基于互质阵列的稀疏频率波形产生装置的结构示意图。

具体实施方式

下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。

本发明的描述中，除非另有明确的限定，设置、连接等词语应做广义理解，所属技术领域技术人员可以结合技术方案的具体内容合理确定上述词语在本发明中的具体含义。

实施例一

参照图1，本发明实施例一提供了一种基于互质阵列的稀疏频率波形产生方法，其中的一种实施例包括但不限于以下步骤：

步骤S100，基于互质阵列获取期望方向图。

在本实施例中，本步骤首先基于互质阵列结构特点，利用互质阵列增加虚拟孔径，获取期望方向图，提高了分辨率，为了后续的波形设计做好前提准备。

步骤S200，通过期望方向图获取自相关矩阵。

在本实施例中，本步骤通过期望方向图获取自相关矩阵，使得发射方向图更加接近期望方向图。

步骤S300，通过自相关矩阵并且结合恒模约束和稀疏频率约束，获取结果波形。

在本实施例中，本步骤通过自相关矩阵计算结果波形，并且在计算结果波形的过程中结合了恒模约束和稀疏频率约束，最后得到结果波形，能够很好地减少波形之间相互干扰的情况，提升检测效果。

参照图2，在步骤S100中，可以包括但不限于以下步骤：

步骤S110，基于互质阵列结构，得到增广虚拟阵列和协方差矩阵。

在本实施例中，本步骤首先基于互质阵列的结构特点，处理得到增广虚拟阵列和协方差矩阵，为了后续获得期望方向图做好准备。

步骤S120，通过增广虚拟阵列和协方差矩阵，得到空间平滑虚拟域协方差矩阵。

在本实施例中，本步骤通过对增广虚拟阵列和协方差矩阵进行运算，获得空间平滑虚拟域协方差矩阵。

步骤S130，通过空间平滑虚拟域协方差矩阵，利用MUSIC算法计算得到期望方向图。

在本实施例中，本步骤通过空间平滑虚拟域协方差矩阵，并且利用MUSIC算法获得期望方向图，在虚拟域匹配奈奎斯特方法实现精确的波达方向估计；其中，MUSIC算法即为传统多重信号分类方法。

在步骤S200中，可以包括但不限于以下步骤：

通过半正定约束由期望方向图得到自相关矩阵。

在本实施例中，本步骤通过半正定约束由期望方向图计算得到自相关矩阵，为了后续的稀疏波形产生做好前提准备。

参照图3，在步骤S300中，可以包括但不限于以下步骤：

步骤S310，通过自相关矩阵并且结合恒模约束，获得雷达波形矩阵最优解。

在本实施例中，本步骤通过自相关矩阵并且结合恒模约束，获得雷达波形矩阵最优解。

步骤S320，通过雷达波形矩阵最优解并且结合稀疏频率约束，构建恒模波形目标函数。

在本实施例中，本步骤通过雷达波形矩阵最优解并且结合稀疏频率约束，构建恒模波形目标函数。

步骤S330，通过恒模波形目标函数结合循环算法，获得结果波形。

在本实施例中，本步骤通过恒模波形目标函数结合循环算法，并且在利用循环算法进行计算的过程中，进行迭代更新，多次利用恒模波形目标函数进行计算，直到满足终止条件，最后输出结果波形，结果波形具有稀疏频率特性，其对应的方向图跟期望方向图之间非常接近，并且可以利用目标函数范数的形式来进行衡量，使得雷达波形在频域上满足稀疏频率约束，而且具有良好的自相关特性。

在本发明的一些具体实施例中，空间平滑虚拟域协方差矩阵可表示为：

其中，MN+1为增广虚拟阵列分割为MN+1个分别包含MN+1个虚拟阵元的线性子阵列，R

在本发明的一些具体实施例中，雷达波形矩阵最优解可表示为：

其中，X表示N×M雷达波形矩阵，R表示自相关矩阵，V表示变换矩阵，N表示发射信号在一个脉冲重复周期的采样数，V

在本发明的一些具体实施例中，自相关矩阵的目标函数可表示为

s.t r

R≥0

其中，优化变量α是与方向图幅度有关的比例值，r

下面利用一个具体实施例对本发明进行进一步的说明：

其中，基于互质阵列结构，采用虚拟阵列信号处理方法。该方法通过计算互质阵列差集数组在虚拟域上形成一个增广虚拟阵列，并利用其2阶等价虚拟域信号的处理，实现基于互质阵列的有效波达方向估计。因为虚拟阵元个数大于物理阵元个数，基于虚拟阵列信号处理的方法在自由度性能方面有着显著提升。虚拟域奈奎斯特匹配方法，通过增广虚拟均匀阵列结构和对应的等价虚拟域信号模型，引入空间平滑技术和虚拟域多重信号分类处理，得到超分辨的波达方向估计结果。在该方法中，互质阵列

0≤m≤2M-1,0≤n≤N-1}

其中，该数组中包含了一个虚拟阵元位置由-MNd到MNd的连续子集，这就表明物理阵元个数为2M+N-1，而扩展互质阵列可以获得2MN+1的自由度。对于空间中K个来自θ＝[θ

推导实际接收信号x(t)与构造的虚拟阵列等价信号y之间的数学映射关系

其中

其中，R

其中，对稀疏频率波形研究的基本方法是对功率谱密度(PSD)和波形自相关函数(ACF)进行约束。针对雷达方向图设计的扩展应用和波形设计的频谱需求，提出了一种互质阵列雷达方向图与稀疏频率波形联合设计方法。首先通过半正定规划由期望方向获得自相关矩阵；然后由自相关矩阵可以求得相应的波形，在求解波形的过程中考虑恒模约束与稀疏频率约束；最后当满足终止条件时循环迭代结束，输出结果波形。得到的具有稀疏频率特性的恒模波形，其对应的方向图跟期望方向图之间的接近程度，用目标函数范数的形式来衡量。雷达波形在频域上满足稀疏频率约束，而且具有良好的自相关特性。

假设互质阵列的均匀部分载频相同为f

对于互质阵列均匀阵元部分，发射信号在一个脉冲重复周期的采样数为N,用x

可得一个脉冲重复周期内辐射到目标θ处的功率为

P(θ)＝w(θ)

其中，w(θ)为阵列导向矢量，w(θ)＝[1 e

对于互质阵列，如果获得期望方向图功率分布，通过设计自相关矩阵R,可以使得发射方向接近期望方向图。其中目标函数为

s.t r

R≥0

优化变量α是与方向图幅度有关的比例值，r

由

V

对应离散信号表示为

A为离散傅里叶变换矩阵所以加上稀疏频率条件后，互质阵列第m个阵元的发射波形x

s.t s

X∈Q V

然后利用循环算法求得波形矩阵X，即获得结果波形。

通过上述方案可知，首先利用互质阵列结构的特点，增加虚拟孔径以获取期望方向图；然后针对期望方向图，获取自相关矩阵，使得发射方向图接近期望方向图；然后针对期望方向图，结合恒模约束和稀疏频率约束，最后获得结果波形，能够很好地减少波形之间相互干扰的情况，提升检测效果。

实施例二

参照图4，本发明实施例二提供了一种基于互质阵列的稀疏频率波形产生装置1000，包括：

处理单元1100，用于基于互质阵列获取期望方向图；

运算单元1200，用于通过所述期望方向图获取自相关矩阵；

计算单元1300，用于通过所述自相关矩阵并且结合恒模约束和稀疏频率约束，获取结果波形。

在本发明的一些具体实施例中，处理单元1100包括：

第一处理单元1110，用于基于互质阵列结构，得到增广虚拟阵列和协方差矩阵；

第二处理单元1120，用于通过所述增广虚拟阵列和所述协方差矩阵，得到空间平滑虚拟域协方差矩阵；

第三处理单元1130，用于通过所述空间平滑虚拟域协方差矩阵，利用MUSIC算法计算得到期望方向图。

需要说明的是，由于本实施例中的互质阵列的稀疏频率波形产生装置与上述实施例一中的互质阵列的稀疏频率波形产生方法基于相同的发明构思，因此，方法实施例一中的相应内容同样适用于本系统实施例，此处不再详述。

通过上述方案可知，首先利用互质阵列结构的特点，增加虚拟孔径以获取期望方向图；然后针对期望方向图，获取自相关矩阵，使得发射方向图接近期望方向图；然后针对期望方向图，结合恒模约束和稀疏频率约束，最后获得结果波形，能够很好地减少波形之间相互干扰的情况，提升检测效果。

实施例三

本发明实施例三还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有基于互质阵列的稀疏频率波形产生装置可执行指令，基于互质阵列的稀疏频率波形产生装置可执行指令用于使基于互质阵列的稀疏频率波形产生装置执行上述的基于互质阵列的稀疏频率波形产生方法，例如，执行以上描述的图1中的方法步骤S100至S300，实现图4中的单元1000-1300的功能。

在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示意性实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。

尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同物限定。