基于双层嵌套阵列的波达方向角估计方法技术

本发明专利技术公开了一种基于双层嵌套阵列的波达方向角估计方法，主要解决现有技术运算量大，识别信源数目少和估计误差大的问题，其实现步骤是：1)分别构建一层嵌套阵列和二层嵌套阵列；2)获取一层嵌套阵列的输出信号和二层嵌套阵列的输出信号；3)根据一层嵌套阵列输出信号与二层嵌套阵列输出信号计算所有四阶累量，形成四阶累量向量；4)根据四阶累量向量构造超完备基，并定义空域稀疏向量；5)将探测空域波达方向角估计转化为求解稀疏约束方程；6)利用凸优化方法，得到空域稀疏向量的最稀疏解；7)根据最稀疏解绘制幅度谱图，得到波达方向角。本发明专利技术在阵元数量有限的情况下大大提高了阵列可识别的信源数目，可用于目标侦察和无源定位。

【技术实现步骤摘要】

本专利技术属于信号处理
，特别涉及一种电磁信号的阵列信号波达方向角估计方法，可用于对飞机、舰船运动目标的侦察与无源定位。
技术介绍
信号的波达方向角DOA估计是阵列信号处理领域的一个重要分支，它是指利用天线阵列对空间声学信号、电磁信号进行感应接收，再运用现代信号处理方法快速准确的估计出信号源的方向，在雷达、声纳、无线通信等领域具有重要应用价值。随着科技的不断进步，对阵列在进行信号波达方向估计时达到的自由度也有越来越高的要求。针对该问题的研究中，出现较早、应用较为广泛的是多重信号分类MUSIC子空间的模型，对于一个L阵元的典型线性均匀阵列，传统的MUSIC类计算方法可检测的信源数目是L-1个。这些算法采用典型的线性均匀阵列，造成估计的信号数目低于阵元数目，目标个数很多时甚至无法识别，导致目标捕获失败。为了能够在阵元数量较少情况下获得尽量大的角度自由度，检测更多的信源，出现了嵌套阵列以及互质阵列。嵌套阵列能够使用M+N个阵元，生成2MN+2N-1个虚拟阵元，可检测MN+N-1个信号。互质阵列能够使用2M+N-1个阵元，生成2MN+2M-1个虚拟阵元，可检测MN+M-1个信号。基于这种信号检测效率较高的嵌套阵列，出现了进一步提高阵列虚拟阵元数的基于四阶累量方法，可使信号估计数量得到较大提升。QingShen等人在其发表的论文“Extensionofnestedarrayswiththefourth-orderdifferenceco-arrayenhancement”(the41stIEEEinternationalconferenceonacoustics,speechandsignalprocessing,shanghai,China,March2016,2991-2995)公开了一种基于四阶累量的嵌套阵列DOA估计方法，该方法使用M+N+P个阵元，生成2[(3P+2)N(M+1)-2P-1]-1个虚拟阵元，可检测(3P+2)N(M+1)-2P-2个信号。但该方法相当于用一个均匀阵列对嵌套阵列进行阵列扩展，没有合理设置阵元位置而不能使阵元达到最大利用率。在实际应用中，如果给定了一定数量的阵元，却不能合理利用这些阵元获得足够多的虚拟阵元，就不能估计足够多的信号，造成侦察和定位资源的浪费和因角度估计误差引起的目标侦察失误。
技术实现思路
本专利技术的目的在于针对上述现有技术存在的不足，提出一种基于双层嵌套阵列的波达方向角估计方法，用一个嵌套阵列对嵌套阵列进行阵列扩展，在阵元数量一定的情况下，提高了可估计信号的数量，避免了不合理利用阵元造成的资源浪费和因角度估计误差引起的目标侦察失误。为实现上述目的，本专利技术的实现步骤包括如下：(1)利用天线接收机形成第一层嵌套阵列，并获取虚拟阵元数：(1a)定义每个天线接收机为一个阵元，用M个天线接收机形成第一均匀线性阵列a，其阵元间距为d；用N个天线接收机形成第二均匀线性阵列b，其阵元间距为(M+1)d；定义第一均匀线性阵列a的第一个阵元为起始阵元，其中，M≥1,N≥1，0<d≤λ/2，λ为入射到嵌套阵列的窄带信号波长；(1b)将第一均匀线性阵列a与第二均匀线性阵列b组合为第一层嵌套阵列：将第二均匀线性阵列b的第一个阵元放置于与起始阵元相距为(M+1)d的位置，再将第二均匀线性阵列b的所有阵元依次插于第一均匀线性阵列a中，得到由M+N个天线接收机形成的第一层嵌套阵列；(1c)计算第一层嵌套阵列的虚拟阵元数：A＝2MN+2N-1；(2)利用天线接收机形成第二层嵌套阵列：(2a)用P个天线接收机形成第三均匀线性阵列f，其阵元间距为Ad，用Q个天线接收机形成第四均匀线性阵列h，其阵元间距为(P+1)Ad，其中，Q≥1，P≥1；(2b)将第三均匀线性阵列f与第四均匀线性阵列h组合为第二层嵌套阵列，即将第三均匀线性阵列f的第一个阵元放置于与起始阵元相距为(A+D)d的位置，将第四均匀线性阵列h的第一个阵元放置于与起始阵元相距为(P+1)Ad的位置，并将第四均匀线性阵列h的所有阵元依次插于第三均匀线性阵列f中，得到由P+Q个天线接收机形成的第二层嵌套阵列，其中D＝(AB-1)/2+MN+N，B＝2PQ+2Q-1；(3)由第一层嵌套阵列和第二层嵌套阵列的天线接收机对空间目标信号进行采样，分别得到第一层嵌套阵列输出信号Y(t)和第二层嵌套阵列输出信号Z(t)，其中，Y(t)＝[y1(t),…,yi(t),…,yM+N(t)]，yi(t)表示第一层嵌套阵列第i个阵元的输出信号，i的取值范围是1≤i≤M+N；Z(t)＝[z1(t),…,zj(t),…,zP+Q(t)]，zj(t)表示第二层嵌套阵列的第j个阵元的输出信号，j的取值范围是1≤j≤P+Q；(4)计算第一层嵌套阵列输出信号Y(t)与第二层嵌套阵列输出信号Z(t)的所有四阶累量R，形成四阶累量向量G；(5)对探测空域波达方向角观测空间进行网格划分，构造出超完备基Φ(θ)，并定义一个Q×1维波达角空域稀疏向量s；(6)根据四阶累量向量G和超完备基Φ(θ)建立约束关系，将探测空域波达方向角估计转化为求解如下稀疏约束方程：min||s||1subjectto||G-Φ(θ)s||2≤ηs≥0其中，||·1表示表示求矩阵一阶范数操作，||·||2表示求矩阵二阶范数操作，subjectto表示约束关系，η为现实环境中的噪声标准差；(7)利用凸优化方法求解上述稀疏约束方程，得到空域稀疏向量s的最稀疏解(8)以波达方向角范围θ＝[θ1,θ2,...,θe,...,θE]的值为x轴坐标，以空域稀疏向量s的最稀疏解的幅度值为y轴坐标，绘制幅度谱图，从该幅度谱图中按照从高到低的顺序寻找幅度值较大的前K个谱峰，这些谱峰的峰值点所对应的x轴坐标即为所求的波达方向角度值，其中，K表示入射到第一层互质阵列与第二层嵌套阵列的空间目标信号个数，且假设空间目标信号在传播过程中加入了均值为零的复高斯白噪声，K≥1。本专利技术与现有技术相比具有以下优点：1)本专利技术将四阶累量应用到第一层嵌套阵列和第二层嵌套阵列形成的双层嵌套阵列模型的DOA估计中，通过使用四阶累量，合理设置阵元位置，可使用M+N+P+Q个阵元获得(2MN+2N-1)(2PQ+2Q-1)+2(P+1)(2MN+2N-1)+2(MN+N-1)个连续虚拟阵元，大大提高了阵列利用率，同时获得更多的虚拟阵元，增加了阵列可识别的信源数目。2)本专利技术采用稀疏表示技术将波达方向角估计转化为稀疏信号的重构，是新理论技术与传统问题的结合；本专利技术利用信号源的空域稀疏性，得到的稀疏解中大系数对应的角度，突破了阵列分辨率的瑞利限，提高了目标侦察和无源定位在低信噪比、低快拍数、低先验知识条件下的DOA估计精度，避免了因角度估计误差引起的目标侦察失误。3)本专利技术利用信号源的空域稀疏性得到的稀疏解中大系数对应的角度，为信号源的波达方向角，无目标的方向角对应的系数约等于0，因此无需预先知道目标的数目。附图说明图1是本专利技术的实现流程图；图2是对本专利技术与现有基于四阶累量嵌套阵列的虚拟阵元数对比图。具体实施方式以下参照附图，对本专利技术的技术方案和效果作进一步的详细说明。参附图1，本专利技术的具体步骤如下：步骤1：利用天线接收机形成第一层嵌本文档来自技高网...

【技术保护点】
一种基于双层嵌套阵列的波达方向角估计方法，包括：(1)利用天线接收机形成第一层嵌套阵列，并获取虚拟阵元数：(1a)定义每个天线接收机为一个阵元，用M个天线接收机形成第一均匀线性阵列a，其阵元间距为d；用N个天线接收机形成第二均匀线性阵列b，其阵元间距为(M+1)d；定义第一均匀线性阵列a的第一个阵元为起始阵元，其中，M≥1,N≥1，0<d≤λ/2，λ为入射到嵌套阵列的窄带信号波长；(1b)将第一均匀线性阵列a与第二均匀线性阵列b组合为第一层嵌套阵列：将第二均匀线性阵列b的第一个阵元放置于与起始阵元相距为(M+1)d的位置，再将第二均匀线性阵列b的所有阵元依次插于第一均匀线性阵列a中，得到由M+N个天线接收机形成的第一层嵌套阵列；(1c)计算第一层嵌套阵列的虚拟阵元数：A＝2MN+2N‑1；(2)利用天线接收机形成第二层嵌套阵列：(2a)用P个天线接收机形成第三均匀线性阵列f，其阵元间距为Ad，用Q个天线接收机形成第四均匀线性阵列h，其阵元间距为(P+1)Ad，其中，Q≥1，P≥1；(2b)将第三均匀线性阵列f与第四均匀线性阵列h组合为第二层嵌套阵列，即将第三均匀线性阵列f的第一个阵元放置于与起始阵元相距为(A+D)d的位置，将第四均匀线性阵列h的第一个阵元放置于与起始阵元相距为(P+1)Ad的位置，并将第四均匀线性阵列h的所有阵元依次插于第三均匀线性阵列f中，得到由P+Q个天线接收机形成的第二层嵌套阵列，其中D＝(AB‑1)/2+MN+N，B＝2PQ+2Q‑1；(3)由第一层嵌套阵列和第二层嵌套阵列的天线接收机对空间目标信号进行采样，分别得到第一层嵌套阵列输出信号Y(t)和第二层嵌套阵列输出信号Z(t)，其中，Y(t)＝[y1(t),…,yi(t),…,yM+N(t)]，yi(t)表示第一层嵌套阵列第i个阵元的输出信号，i的取值范围是1≤i≤M+N；Z(t)＝[z1(t),…,zj(t),…,zP+Q(t)]，zj(t)表示第二层嵌套阵列的第j个阵元的输出信号，j的取值范围是1≤j≤P+Q；(4)计算第一层嵌套阵列输出信号Y(t)与第二层嵌套阵列输出信号Z(t)的所有四阶累量R，形成四阶累量向量G；(5)对探测空域波达方向角观测空间进行网格划分，构造出超完备基Φ(θ)，并定义一个Q×1维波达角空域稀疏向量s；(6)根据四阶累量向量G和超完备基Φ(θ)建立约束关系，将探测空域波达方向角估计转化为求解如下稀疏约束方程：min||s||1subject to||G‑Φ(θ)s||2≤ηs≥0其中，||·||1表示表示求矩阵一阶范数操作，||·||2表示求矩阵二阶范数操作，subjectto表示约束关系，η为现实环境中的噪声标准差；(7)利用凸优化方法求解上述稀疏约束方程，得到空域稀疏向量s的最稀疏解(8)以波达方向角范围θ＝[θ1,θ2,...,θe,...,θE]的值为x轴坐标，以空域稀疏向量s的最稀疏解的幅度值为y轴坐标，绘制幅度谱图，从该幅度谱图中按照从高到低的顺序寻找幅度值较大的前K个谱峰，这些谱峰的峰值点所对应的x轴坐标即为所求的波达方向角度值，其中，K表示入射到第一层互质阵列与第二层嵌套阵列的空间目标信号个数，且假设空间目标信号在传播过程中加入了均值为零的复高斯白噪声，K≥1。...

【技术特征摘要】
1.一种基于双层嵌套阵列的波达方向角估计方法，包括：(1)利用天线接收机形成第一层嵌套阵列，并获取虚拟阵元数：(1a)定义每个天线接收机为一个阵元，用M个天线接收机形成第一均匀线性阵列a，其阵元间距为d；用N个天线接收机形成第二均匀线性阵列b，其阵元间距为(M+1)d；定义第一均匀线性阵列a的第一个阵元为起始阵元，其中，M≥1,N≥1，0<d≤λ/2，λ为入射到嵌套阵列的窄带信号波长；(1b)将第一均匀线性阵列a与第二均匀线性阵列b组合为第一层嵌套阵列：将第二均匀线性阵列b的第一个阵元放置于与起始阵元相距为(M+1)d的位置，再将第二均匀线性阵列b的所有阵元依次插于第一均匀线性阵列a中，得到由M+N个天线接收机形成的第一层嵌套阵列；(1c)计算第一层嵌套阵列的虚拟阵元数：A＝2MN+2N-1；(2)利用天线接收机形成第二层嵌套阵列：(2a)用P个天线接收机形成第三均匀线性阵列f，其阵元间距为Ad，用Q个天线接收机形成第四均匀线性阵列h，其阵元间距为(P+1)Ad，其中，Q≥1，P≥1；(2b)将第三均匀线性阵列f与第四均匀线性阵列h组合为第二层嵌套阵列，即将第三均匀线性阵列f的第一个阵元放置于与起始阵元相距为(A+D)d的位置，将第四均匀线性阵列h的第一个阵元放置于与起始阵元相距为(P+1)Ad的位置，并将第四均匀线性阵列h的所有阵元依次插于第三均匀线性阵列f中，得到由P+Q个天线接收机形成的第二层嵌套阵列，其中D＝(AB-1)/2+MN+N，B＝2PQ+2Q-1；(3)由第一层嵌套阵列和第二层嵌套阵列的天线接收机对空间目标信号进行采样，分别得到第一层嵌套阵列输出信号Y(t)和第二层嵌套阵列输出信号Z(t)，其中，Y(t)＝[y1(t),…,yi(t),…,yM+N(t)]，yi(t)表示第一层嵌套阵列第i个阵元的输出信号，i的取值范围是1≤i≤M+N；Z(t)＝[z1(t),…,zj(t),…,zP+Q(t)]，zj(t)表示第二层嵌套阵列的第j个阵元的输出信号，j的取值范围是1≤j≤P+Q；(4)计算第一层嵌套阵列输出信号Y(t)与第二层嵌套阵列输出信号Z(t)的所有四阶累量R，形成四阶累量向量G；(5)对探测空域波达方向角观测空间进行网格划分，构造出超完备基Φ(θ)，并定义一个Q×1维波达角空域稀疏向量s；(6)根据四阶累量向量G和超完备基Φ(θ)建立约束关系，将探测空域波达方向角估计转化为求解如下稀疏约束方程：min||s||1subjectto||G-Φ(θ)s||2≤ηs≥0其中，||·||1表示表示求矩阵一阶范数操作，||·||2表示求矩阵二阶范数操作，subjectto表示约束关系，η为现实环境中的噪声标准差；(7)利用凸优化方法求解上述稀疏约束方程，得到空域稀疏向量s的最稀疏解(8)以波达方向角范围θ＝[θ1,θ2,...,θe,...,θE]的值为x轴坐标，以空域稀疏向量s的最稀疏解的幅度值为y轴坐标，绘制幅度谱图，从该幅度谱图中按照从高到...

【专利技术属性】
技术研发人员：蔡晶晶，宗汝，王振华，李鹏，赵晗希，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西;61