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

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


背景技术：

2.稀疏阵列在雷达相控阵波达方向估计上有着重要作用，是相关稀疏算法估计的关键结构模型。对稀疏矩阵的每个阵元的接收信号进行协方差处理并矢量化，得到等效的虚拟差分阵列，其虚拟阵元位于物理阵元的差分值处，选取差分阵列中最长的对称连续段，就可以用相关稀疏算法进行处理，从而实现相控阵波达方向估计的目的。相较于传统均匀线性阵列，稀疏阵列传感器数量更少，阵列物理成本大幅降低，硬件复杂度和数据处理复杂度降低；同时在性能方面，方位估计性能和空间探测能力远超传统线性阵。由于稀疏阵列的良好性能，近年来稀疏阵列结构的研究受到广泛关注。最小冗余阵及最小孔洞阵是两种早期被提出的非均匀线列阵列。它们的结构能够检测比同数量均匀线性阵列更多的空间目标并实现结构上的冗余度最小化，但其缺点是没有解析表达式，在阵元数较多的情况下只能通过计算机暴力求解，不符合实际的需求。最近几年提出的具有解析表达式的两种新型阵列为嵌套阵和互质阵。其中，文献[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个信号的功率，1r＝vec(ir)，为噪声功率；找出稀疏阵的虚拟差分共阵中的最大连续段，然后可将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)中所述稀疏阵列的差分虚拟阵位置可以通过计算获得；从中找出最大连续段[-lu,lu]d，即可得到阵列自由度lu，且稀
疏阵列的自由度为lu＝mn 6n-12。
[0015]
进一步地，步骤(4)中所述稀疏阵列，在给定总阵元数r＝n m的情况下，当n和m分别取时，虚拟阵列的自由度lu最大，此时该结构能达到的最大自由度为
[0016]
本发明与现有技术相比具有以下优点：
[0017]
(1)本发明基于v＝v1n v2＝(v
1-2)n (2n v2)的数学事实，部分虚拟位置的构建采用少量的中间稀疏均匀阵元和更长间距的左右侧子阵阵元拉伸两侧密集阵。该结构保证了阵元总数不变，但阵元排布更加稀疏，虚拟差分阵的自由度更大。本发明在双边稀疏嵌套阵结构的基础上，进一步稀疏化设计两侧子阵间距，使得两侧子阵的扩展间距由n变为2n，同时保证其虚拟阵最长连续段增加。其虚拟孔径长度由双边稀疏嵌套阵的nm 3.5n-3.5提高到nm 6n-12。n越大，其自由度提高地越多；
[0018]
(2)本发明还给出了双边拓展两次嵌套阵的自由度理论上限，为估计该类结构的探测能力提出理论指导并对双边边拓展多次间距增强提出修改方式。
附图说明
[0019]
图1是本发明提出的双边二次拓展嵌套阵结构示意图，间隔为归一化间隔，省略基本单元d。
[0020]
图2是阵元数为r＝32(n＝31,m＝1)时的双边二次拓展嵌套阵结构图，间隔为归一化间隔。
[0021]
图3是本发明设计结构与其他五种稀疏结构的波达方向估计均方根误差随信噪比变化对比图。
[0022]
图4是本发明设计结构与其他五种稀疏结构的波达方向估计均方根误差随快拍数变化对比图。
[0023]
图5是本发明设计结构与其他五种稀疏结构的波达方向估计均方根误差随待测信源数变化对比图。
具体实施方式
[0024]
下面结合附图及实施例对本发明的内容做进一步的说明。
[0025]
以入射信号频率20ghz，阵元总数r＝32为例，对本发明的双边二次拓展嵌套阵(如图1所示)构建进行详细说明。
[0026]
步骤一，确定间距基本单元长度：
[0027]
计算入射信号波长λ＝c/f＝3.0
×
108/20
×
109＝1.5cm；用d表示阵元间距的基本单元，将其取为半波长，即d＝λ/2＝7.5mm；
[0028]
步骤二，计算相关阵列参数：
[0029]
根据阵列给定的总阵元数r＝32和r＝n m,选择一组n和m，满足n≥31,m≥1；本实施例选择n＝31,m＝1，结构如图2所示，再根据n计算阵元间隔参数l＝(n-3)/4＝7；
[0030]
步骤三，计算阵元位置：
array，ts2-sna)选取n＝25,m＝25，本发明的双边稀疏嵌套阵(two-side sparse nested array,ts-sna)选取n＝28,m＝22。
[0044]
假设200个不相关信号入射到上述六种稀疏结构上，入射角度均匀分布在-60
°
到60
°
之间。实施300次蒙特卡洛实验，将信号快拍数为1000次情况下的估计均方根误差(root-mean-square error，rmse)随入射信号信噪比(signal to noise ratio，snr)变化的曲线绘制在图3中。rmse随快拍数变化的曲线绘制在图4中，此时信噪比设为5db。rmse随待检测信源数数变化的曲线绘制在图5中，此时信噪比设为5db，信号快拍数为1000次。由三张图可见，本发明提出的ts2-sna比其他五种结构的在信噪比、快拍数、待测信源数方面的实验rmse结果都是最好的。这是因为，na、anai-2、os-ena、misc、ts-sna和ts2-sna的阵列自由度分别为650，662,696,696,710和769。本发明所提的ts2-sna具有相比于前五种结构更高的自由度。因此对目标方位估计的精度更高。
[0045]
综上所述，本发明所提设计方法针对嵌套阵系列，在双边稀疏嵌套阵结构的基础上，进一步拓展扩大了两侧子阵的间距并进一步细化两侧子阵传感器数量的排布，给出了一种关于双边拓展两次嵌套阵结构的完整构建准则。同时，在只给定阵元总数的情况下，通过解析闭合表达式的计算，获得一种具有比现有嵌套阵类结构更高自由度的稀疏阵列结构。在参数选定方面，给出最优选型的配置方案。该结构在自由度方面较目前其它嵌套阵系列有显著提升，在阵列波达方向估计方面和空间目标探测能力方面具有更优越的性能。