一种新型稀疏阵列的DOA估计方法技术

本发明专利技术公开了一种新型稀疏阵列设计方法，属于阵列信号处理领域。包括步骤：对原始嵌套阵列的两组子阵列进行重新排列设计，以构建新型嵌套阵列；计算新型嵌套阵列的差分共阵和求和共阵及其连续自由度；对目标信号源进行空间采样，计算新型阵列接收信号的协方差矩阵并将其向量化处理；构造选择矩阵，提取向量化矩阵中的非零元素；利用MUSIC算法实现目标信号源的高精度DOA估计。的高精度DOA估计。的高精度DOA估计。

【技术实现步骤摘要】
一种新型稀疏阵列的DOA估计方法


[0001]本专利技术涉及目标检测和信号处理领域，具体涉及DOA(Directionofarrival)估计的新型稀疏阵列设计方法。

技术介绍

[0002]新技术革命浪潮下，高超声速目标、声呐雷达信号检测、复杂电子对抗等现代化技术手段的快速发展给现有空间和安全防控领域带来严峻考验。因此，利用新型阵列设计理论开展多目标DOA估计算法研究，实现对敌方目标的快速有效定位，对于空间和安全防控领域具有重要的现实意义。
[0003]相比于均匀线阵，稀疏阵列具有更大的虚拟延迟、更高的自由度和更弱的互耦效应，带来了估计性能的极大提升，从而成为近年来DOA估计问题的研究热点。在所有的稀疏阵列设计中，嵌套阵列和互质阵列是最具代表性的阵列设计方法，已有方法大多是基于上述方法的改进。嵌套阵列由单位间距的均匀子阵和较大间距的稀疏子阵排列构成，可以使用N个传感器估计O(N2)个信源，但嵌套阵列中的均匀子阵存在耦合效应，从而影响其估计性能。

技术实现思路

[0004]本专利技术解决的技术问题是：针对现有模型DOA估计问题，提出一种新型稀疏阵列设计方法，实现高效的DOA估计。
[0005]本专利技术的技术解决方案是：一种新型稀疏阵列设计方法，步骤如下：
[0006]步骤1：对原始嵌套阵列的两组子阵列进行重新排列设计，以构建新型嵌套阵列；
[0007]步骤2：计算新型嵌套阵列的差分共阵和求和共阵及其连续自由度；
[0008]步骤3：对目标信号源进行空间采样，计算新型阵列接收信号的协方差矩阵并将其向量化处理；
[0009]步骤4：构造选择矩阵，提取向量化矩阵中的非零元素；
[0010]步骤5：利用MUSIC算法实现目标信号源的DOA估计。
[0011]进一步地，步骤1所述的新型阵列由两部分构成，在一条直线上排列的非均匀嵌套阵列，分别分为子阵一和子阵二。子阵一包含N个阵元，阵元间距为Md；子阵二包含M个阵元，阵元间距为其中M＝2,d代表半波长。新型阵列的构型表达式为S＝S1∪S2,
[0012][0013]其中,S1和S2分别代表子阵1和子阵2的物理阵元排列位置。
[0014]进一步地，由于虚拟阵元的位置关于0位置对称，为了便于说明，仅取虚拟阵元的
非负部分进行分析。步骤2中的稀疏阵列总延迟由差分共阵与求和共阵组成，而差分共阵包括自差集和互差集而求和共阵包括自和集与互和集相应的集合可表示为：
[0015][0016][0017][0018]进一步地，步骤3中的接收信号表示为：其中，x(t)为接收信号矢量，s
k
(t)为发送信源，为相应于θ
k
的导向矢量，d
m
代表第m个阵元与参考阵元之间的距离，A是阵列流形，为相位偏移矩阵,Φ
k
是第k个信源的相位，S(t)为信源矢量，n(t)是功率为的加性高斯白噪声，N(t)为相应的噪声矩阵；
[0019]协方差矩阵可表示为：
[0020][0021]矢量化后的协方差矩阵为：
[0022][0023]进一步地，步骤4中采用选择矩阵提取矢量化矩阵中的非零元素：
[0024]Z
c
＝J
·
z，
[0025]其中J代表选择矩阵。
[0026]本专利技术的有益效果：针对现有设计方法存在孔洞，连续滞后少的问题，通过改变阵元排布方式，设计一种虚拟阵元无孔的嵌套阵列。给出了联合利用差分共阵与求和共阵条件下连续延迟的数学表达式，大大提高了自由度，增加了DOA估计精度。
附图说明
[0027]为使得本专利技术中的技术方案，下面将会选取一些附图对本专利技术进行阐述。下述附图仅为本专利技术的一些实例；使用者可根据该附图获得其他类似的附图。其中：
[0028]图1为本专利技术的流程示意图；
[0029]图2为同条件下本专利技术与原始互质阵列、原始嵌套阵列的空间谱图；
[0030]图3为同条件下本专利技术与原始互质阵列、原始嵌套阵列的均方根误差随信噪比变
化图；
[0031]图4为同条件下本专利技术与原始互质阵列、原始嵌套阵列的均方根误差随快拍变化图；
具体实施方式
[0032]下面将更加深刻，完整的阐述本专利技术实施案例中的技术方案。所述实施案例为部分实施案例。使用者可根据本专利技术，不需付出创造性劳动的获得其他实施案例。该“其他实施案例”均属于本专利技术保护的范围。
[0033]如图1所示，一种无孔阵列DOA估计流程图，包括以下步骤：
[0034]对原始嵌套阵列进行重新排列设计构建新型嵌套阵列，新阵列包括子阵一和子阵二。子阵一包含N个阵元，阵元间距为Md；子阵二包含M个阵元，阵元间距为其中M＝2,d代表半波长。新型阵列的构型表达式为S＝S1∪S2,
[0035][0036]步骤2：计算新模型的差分共阵、求和共阵、以及连续自由度，包括：
[0037]步骤2.1：计算差分共阵与求和共阵：稀疏阵列总延迟由差分共阵与求和共阵组成，而差分共阵包括自差集和互差集而求和共阵包括自和集与互和集由于虚拟阵元的位置关于0位置对称，为了便于说明，仅取虚拟阵元的非负部分进行分析。
[0038]相应的集合可表示为：
[0039][0040][0041][0042]步骤2.2：计算连续自由度：
[0043]对于阵元数为Q＝M+N
‑
1的新型嵌套阵列，其连续自由度为：4MN
‑
3。
[0044]步骤3：对目标信号源进行采样，并对接收信号的协方差矩阵进行向量化处理，包括：
[0045]步骤3.1：对目标信号源进行采样，得到接收信号向量表示为：
[0046][0047]其中，x(t)为接收信号，s
k
(t)为发送信号，为相应于θ
k
的导向矢量，d
m
代表第m个阵元与参考阵元之间的距离，A是阵列流形，为相位偏移矩阵,Φ
k
是第k个信源的相位，S(t)为信源矢量，n(t)是功率为的加性高斯白噪声，N(t)为相应的噪声矩阵；
[0048]步骤3.2：计算接收到的协方差矩阵，可表示为：
[0049][0050]步骤3.3：对协方差矩阵矢量化处理：
[0051][0052]步骤4：采用选择矩阵提取矢量化矩阵中的非零元素：
[0053]Z
c
＝J
·
z，
[0054]其中J代表选择矩阵。
[0055]步骤5：利用MUSIC算法实现目标信号源的DOA估计。
[0056]具体实施例
[0057]现在结合附图2——4对本专利技术作进一步详细的说明。
[0058]采用阵元数目Q＝M+N
‑
1＝6的物理阵列，其中M＝2,N＝5，子阵1包含N个间距为Md的物理阵元，子阵2包含M个间距为的物理阵元，物理阵元排列位置本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种新型稀疏阵列的DOA估计方法，其特征在于，包括以下步骤：步骤1：对原始嵌套阵列的两组子阵列进行重新排列设计，以构建新型嵌套阵列；步骤2：计算新型嵌套阵列的差分共阵和求和共阵及其连续自由度；步骤3：对目标信号源进行空间采样，计算新型阵列接收信号的协方差矩阵并将其向量化处理；步骤4：构造选择矩阵，提取向量化矩阵中的非零元素；步骤5：利用MUSIC算法实现目标信号源的DOA估计。2.根据权利要求1所述的一种新型稀疏阵列的DOA估计方法，其特征在于，步骤1所述新型阵列由两部分构成，在一条直线上排列的非均匀嵌套阵列，分别分为子阵一和子阵二，子阵一包含N个阵元，阵元间距为Md；子阵二包含M个阵元，阵元间距为其中M＝2,d代表半波长。新型阵列的构型表达式为S＝S1∪S2,其中,S1和S2分别代表子阵1和子阵2的物理阵元排列位置。3.根据权利要求1所述的一种新型稀疏阵列的DOA估计方法，其特征在于，步骤2中所述的稀疏阵列总延迟由差分共阵与求和共阵组成，而差分共阵包括自差集和互差集而求和共阵包括自和集与...

【专利技术属性】
技术研发人员：周彬，刘可，
申请(专利权)人：哈尔滨理工大学，
类型：发明
国别省市：