一种用于波达方向估计的双边二次拓展嵌套阵设计方法技术

本发明专利技术公开了一种用于波达方向估计的双边二次拓展嵌套阵设计方法，该策略在现有嵌套阵相关结构的基础上进一步增加了两侧子阵列的元素间距，产生更长的无孔虚拟连续段，使自由度得到了极大的提升。为了保证差分共阵连续性的最大化，提出了四种条件，完善了实现物理阵列结构的步骤，并在此策略的基础上，提出了一种全新的具有闭合表达式的阵列结构，在给定阵元数不变的情况下，具有更高的自由度和更好的目标探测和方位估计性能。的目标探测和方位估计性能。的目标探测和方位估计性能。

【技术实现步骤摘要】
一种用于波达方向估计的双边二次拓展嵌套阵设计方法


[0001]本专利技术涉及一种用于波达方向估计的雷达相控稀疏阵列设计方案，主要用于被动稀疏阵列目标探测与方位估计，在传感器总数有限且固定的情况下，通过对阵元布放位置结构进行设计和优化，能够提高所产生的虚拟差分共阵孔径大小，最终实现被动探测情况下的较优估计方案。

技术介绍

[0002]稀疏阵列在雷达相控阵波达方向估计上有着重要作用，是相关稀疏算法估计的关键结构模型。对稀疏矩阵的每个阵元的接收信号进行协方差处理并矢量化，得到等效的虚拟差分阵列，其虚拟阵元位于物理阵元的差分值处，选取差分阵列中最长的对称连续段，就可以用相关稀疏算法进行处理，从而实现相控阵波达方向估计的目的。相较于传统均匀线性阵列，稀疏阵列传感器数量更少，阵列物理成本大幅降低，硬件复杂度和数据处理复杂度降低；同时在性能方面，方位估计性能和空间探测能力远超传统线性阵。由于稀疏阵列的良好性能，近年来稀疏阵列结构的研究受到广泛关注。最小冗余阵及最小孔洞阵是两种早期被提出的非均匀线列阵列。它们的结构能够检测比同数量均匀线性阵列更多的空间目标并实现结构上的冗余度最小化，但其缺点是没有解析表达式，在阵元数较多的情况下只能通过计算机暴力求解，不符合实际的需求。最近几年提出的具有解析表达式的两种新型阵列为嵌套阵和互质阵。其中，文献[1](P.Pal and P.P.Vaidyanathan,“Multiple level nested array:An efficient geometry for 2qth order cumulant based array processing,”IEEE Trans.Signal Process.,vol.60,no.3,pp.1253
–
1269,Mar.2012.)中的嵌套阵能获得无孔差分虚拟阵，阵列孔径大小比传统均匀线阵提高了一倍以上，使嵌套阵的进一步优化布局成为趋势。目前最新的一些成果包括文献[2](J.Liu,Y.Zhang,Y.Lu,S.Ren,and S.Cao,“Augmented nested arrays with enhanced DOF and reduced mutual coupling,”IEEE Trans.Signal Process.,vol.65,no.21,pp.5549
–
5563,Nov.2017.)中的ANA、文献[3](Z.Zheng,W.
‑
Q.Wang,Y.Kong,and Y.D.Zhang,“MISC array:A new sparse array design achieving increased degrees of freedom and reduced mutual coupling effect,”IEEE Trans.Signal Process.,vol.67,no.7,pp.1728
–
1741,Apr.2019.)中的MISC、文献[4](S.Ren,W.Dong,X.Li,W.Wang and X.Li,“Extended Nested Arrays for Consecutive Virtual Aperture Enhancement,”IEEE Signal Processing Letters,vol.27,pp.575
‑
579,2020.)中的OS
‑
ENA和TS
‑
ENA等多种基于嵌套阵的结构被先后发现，阵列整体自由度被进一步提升。TS
‑
ENA相较于最小冗余阵，在结构两侧自由度提升、拓展方面仍有较大的改进空间，并且其中稀疏阵元的提升规律存在且目前未被完整地总结出来。
[0003]在这种自由度提升规律的指导下，如何发现更多具有更高自由度的一维稀疏嵌套系结构并且设计出来，从而令虚拟差分阵元的最大对称连续段更长，使阵列的总体自由度进一步增加，成为了当前阵列结构设计中的主要问题。为此，稀疏阵列的设计问题变得十分
重要，是阵列信号处理其他领域突破的重要支撑。

技术实现思路

[0004]本专利技术的目的在于针对现有技术的不足，通过研究适用于雷达相控阵波达方向估计的阵列优化准则和基于此准则的具体闭合结构，在双边稀疏嵌套阵结构的基础上，进一步拓展两侧子阵的间距并进一步细化两侧子阵传感器数量的排布，从而在传感器数目一定的情况下，获得比现有稀疏嵌套阵相关优化阵型更大的虚拟差分共阵和更大的阵列自由度。
[0005]为实现上述目的，本专利技术的设计思路如下：延续TS
‑
ENA中嵌套三段式结构，将右侧子阵的部分阵元与中间子阵最右侧阵元位置的距离由l增大到2N+l，并且将左侧子阵的部分阵元与Nd位置处阵元的距离由l'增大到2N+l'，使得虚拟差分阵元的最大对称连续段更长，从而使得阵列的总体自由度进一步增加。
[0006]本专利技术是通过以下技术方案实现的，该方案包括以下步骤：
[0007](1)用d表示阵元间距的基本单元，将其取为半波长，即d＝λ/2，λ表示阵列入射信号的波长；
[0008](2)根据阵列给定的总阵元数R，确定双边二次拓展嵌套阵中大间隔均匀线列子阵的阵元间距N和阵元数M，满足R＝N+M，其中N≥31,M≥1；
[0009](3)根据N计算阵元间隔参数l＝(N
‑
3)/4，其中l＞6；
[0010](4)将第一个阵元放置在1位置处；根据阵元间隔参数l和M，确定双边二次拓展嵌套阵的所有阵元位置集合其中
[0011][0012](5)K个远场窄带信号以角度入射到阵元位置集合为的稀疏阵上，t时刻接收阵列的信号表示为x(t)＝As(t)+n(t)，其中A为阵列流型矩阵，s(t)为入射信号矢量，n(t)为均匀白噪声；对x(t)的协方差矩阵进行矢量化，可以得到矢量化后的协方差矩阵其中R
x
为x(t)的协方差矩阵，[
·
]*
表示矩阵的共轭，
⊙
表示KR积，为第k个信号的功率，1
R
＝vec(I
R
)，为噪声功率；找出稀疏阵的虚拟差分共阵中的最大连续段，然后可将SS
‑
MUSIC等算法应用于该虚拟差分共阵进行波达方向估计，求解待测方位角度。
[0013]进一步地，步骤(4)中所述稀疏阵列包含三个组成部分：第一部分表示双边二次拓展嵌套阵位于[1,12l+6]d范围内的阵元位置；第二部分是中间部分间隔稀疏均匀线阵，包含M个间距为Nd的均匀分布阵元，其阵元位置分布在(12l+6+Nr)d,r∈[1,M]；第三部分是位于[14l+8+(4l+3)M,24l+12+(4l+3)M]d范围内的阵元位置；和分别位于的左侧和右侧。
[0014]进一步地，步骤(5)中所述稀疏阵列的差分虚拟阵位置可以通过计算获得；从中找出最大连续段[
‑
L
u
,L
本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种用于波达方向估计的双边二次拓展嵌套阵设计方法，其特征在于，包含以下步骤：(1)用d表示阵元间距的基本单元，将其取为半波长，即d＝λ/2，λ表示阵列入射信号的波长；(2)根据阵列给定的总阵元数R，确定双边二次拓展嵌套阵中大间隔均匀线列子阵的阵元间距N和阵元数M，满足R＝N+M，其中N≥31,M≥1；(3)根据N计算阵元间隔参数l＝(N
‑
3)/4，其中l＞6；(4)将第一个阵元放置在1位置处；根据阵元间隔参数l和M，确定双边二次拓展嵌套阵的所有阵元位置集合其中(5)K个远场窄带信号以角度入射到阵元位置集合为的稀疏阵上，t时刻接收阵列的信号表示为x(t)＝As(t)+n(t)，其中A为阵列流型矩阵，s(t)为入射信号矢量，n(t)为均匀白噪声；对x(t)的协方差矩阵进行矢量化，可以得到矢量化后的协方差矩阵其中R
x
为x(t)的协方差矩阵，[
·
]
*
表示矩阵的共轭，
⊙
表示KR积，KR积，为第k个信号的功率，1
R
＝vec(I
R
)，为噪声功率；找出稀疏阵的虚拟差分共阵中的最大连续段，然后可将SS
‑
MUSIC等算法应用于该虚拟差分共阵进行波达方向估计，求解待测方位角度。2.根据权...

【专利技术属性】
技术研发人员：任仕伟，王树建，王贵愚，李翔南，王卫江，
申请(专利权)人：北京理工大学，
类型：发明
国别省市：