一种基于多层互质阵列的相干源波达方向估计方法技术

本发明专利技术提供一种基于多层互质阵列的相干源波达方向估计方法，包含步骤：S1、接收端利用M个阵元构建包含Q层稀疏均匀直线子阵列的多层互质阵列；S2、对每层稀疏均匀直线子阵列的接收信号建模，根据所述建模采集T个采样快拍，得到每层稀疏均匀直线子阵列接收信号的采样协方差矩阵；S3、重构所述采样协方差矩阵的秩，获得对应的平滑协方差矩阵；S4、对所述平滑协方差矩阵使用求根多重信号分类方法，获得模糊波达方向集合；S5、构建代价函数，从模糊波达方向集合中选出每层稀疏均匀直线子阵列对目标测得的模糊波达方向估计值；S6、对S5的结果求均值，获得所述多层互质阵列对目标测得的波达方向估计值。

A method of DOA estimation of coherent sources based on multi-layer coprime array

【技术实现步骤摘要】
一种基于多层互质阵列的相干源波达方向估计方法
本专利技术涉及信号处理
，具体是一种基于多层互质阵列的相干源波达方向估计方法。
技术介绍
目标的波达方向(Directionofarrival，DOA)估计为阵列信号处理的一个重要分支，其在雷达、声呐以及无线电通信领域有着广泛的应用，其主要工作原理为通过空间分布的传感器对目标的回波信号进行空间采样，然后利用回波在不同传感器之间的时延关系获得目标的DOA估计。阵列的角度分辨力是指阵列对两个临近目标的分辨能力，该能力由阵列的孔径决定，阵列的孔径越大，对应的角度分辨力越高。在满足空间采样的情况下，可以通过增加阵列中的阵元数目提高阵列的孔径，从而提高阵列的角度分辨力，而更多的阵元数增加了系统的成本和复杂度。对于成本敏感型系统，在其设计过程中，可以使用的阵元数是固定的，如何在给定阵元数目下提高阵列系统的角度分辨力具有重要意义。一方面，采用稀疏均匀直线阵列配置可以获得大的阵列孔径，然而，稀疏均匀直线阵列不满足空间采样定理，存在空间模糊问题，不利于目标的探测；另一方面，可以利用非均匀阵列配置获得大的物理孔径，同时解决空间模糊问题。互质阵列作为一种具有系统化结构的稀疏非均匀阵列，近年来备受学术界关注，现有的基于互质阵列的DOA估计方法主要是在互质阵列二阶等价虚拟信号的基础上进行的。相比于均匀直线阵列，互质阵列具有较高的阵列自由度，更高的分辨力，然而，获取二阶等价虚拟信号前提条件是目标的回波信号是非相干的。当多径效应出现时，目标的回波之间往往是相干的，此时，不能对互质阵列接收信号进行二阶等价变换。传统的多重信号分类方法MUSIC不能直接对相干源信号进行处理，需要先进行空间平滑预处理。然而，空间平滑预处理需要在均匀直线阵列上进行，由于互质阵列为稀疏非均匀阵列，因此只能将互质阵列分解成为两个稀疏均匀直线子阵列分别进行空间平滑处理，然后进行DOA估计。由于两个子阵列都是空间欠采样的，因此每个子阵列都存在空间模糊问题。专利号为201710041947.0[P]的一种基于压缩感知的互质阵列高精度波达方向估计方法，该专利使用压缩感知技术解决现有技术中估计精度与计算复杂度之间的问题。该方法可以处理相干源目标，但是存在网格失配问题，且在低信噪比情况下该方法存在较大的估计误差。本专利中提出的方法不存在网格失配问题，在低信噪比情况下依旧具有良好的估计精度。专利号为201710302914.7[P]的基于内插虚拟阵列信号原子范数最小化的互质阵列波达方向估计方法，该专利使用内插技术解决现有技术中虚拟阵列的非均匀性所导致的信息损失问题。该方法受限于非相干信号源。此外，该方法使用了范数优化方法，存在计算复杂度高、不利于工程实现的问题。专利号为201810464602.0[P]的基于角度一空间频率域快速傅里叶变换的互质阵列波达方向估计方法，该专利使用频率域快速傅里叶变换实现互质阵列波达方向估计,解决现有技术计算复杂度高的问题，但是该方法使用的是非相干信号模型，不能处理相干源目标。因此该方法不适用于本专利技术。论文名为“SufficientConditionsforAliasing-freeDirectionofarrivalEstimationinPeriodicSpatialSpectra”(《IEEETransactiononAntennasandPropagation》，VOL41，NO.4，1993)指出，利用互质阵列两个子阵列阵元间距的互质特性，通过子阵列的MUSIC谱峰匹配处理，可以实现解模糊处理。该论文进一步指出，对于由M+1稀疏子阵列构成的稀疏非均匀阵列，最多能处理M个目标。由于互质阵列只有2个子稀疏均匀子阵列，因此，其最多只能处理1个信号，不能对相干信号进行解模糊处理。此外，MUSIC方法对全空间进行谱搜索需要较高的计算量，在实时性要求较高的应用场景下面临一定的挑战，且在实际系统中的硬件实现较为困难。论文名为“PartialSpectralSearch-BasedDOAestimationMethodforCo-primeLinearArrays”(《ElectronicsLetters》，VOL.51，NO.24，2015).利用对子阵列空间谱的周期特性，利用局部谱峰搜索代替全局谱峰搜索，同时提出了模糊谱峰匹配方案。虽然该方法在一定程度上降低了计算量，但是仍然需要进行谱峰搜索，计算量依旧很大。此外，该文献中的模糊谱峰匹配仅适用于双子阵列情况，难以推广到多子阵列情形。
技术实现思路
本专利技术的目的在于提供一种基于多层互质阵列的相干源波达方向估计方法，通过建立多层互质阵列，并结合求根MUSIC(multiplesignalclassification多重信号分类)方法进行相干源波达方向的估计，降低了相干源波达方向估计的计算复杂度。为了达到上述目的，本专利技术提供一种基于多层互质阵列的相干源波达方向估计方法，用于估计空间内来自K个方向的远场窄带相干信号源的波达方向，包含步骤：S1、接收端利用M个阵元构建包含Q层稀疏均匀直线子阵列的多层互质阵列；所述多层互质阵列中的第q层稀疏均匀直线子阵列包含N＝M/Q个阵元，该N个阵元位于同一条直线上，相邻阵元的间距为Mqd；其中q∈[1,Q]，d＝λ/2，λ为所述窄带相干信号源的波长；M1，…，Mq均为正整数且两两互质；Q＞K且S2、对第q层稀疏均匀直线子阵列的接收信号建模，根据所述建模采集T个采样快拍，得到第q层稀疏均匀直线子阵列接收信号的采样协方差矩阵q∈[1,Q]；S3、重构所述采样协方差矩阵的秩，获得维度为L×L的平滑协方差矩阵q∈[1,Q]；S4、对所述平滑协方差矩阵使用求根多重信号分类方法，获得模糊波达方向集合其中为K个信号源各自的角度模糊值，k∈[1,K]；为与对应的角度真实值，m∈[1,Mq]；q∈[1,Q]；S5、构建代价函数F(θ)；选取满足为中的K个最小值；将作为第q层稀疏均匀直线子阵列对K个信号源的模糊波达方向估计值；S6、将作为所述多层互质阵列对K个信号源的波达方向估计值。所述步骤S1中，不同层的稀疏均匀直线子阵列互相平行；每层稀疏均匀直线子阵列的首个阵元重叠。所述步骤S2具体包含：S21、对第q层稀疏均匀直线子阵列的接收信号建模为：其中，s0(t)为接收信号波形，ρk、φk分别为第k个相干信号波形的幅度、附加相位，nq(t)为与信号源相互独立的高斯白噪声分量，θk为信号源的波达方向角，aq(θk)为第q层稀疏均匀直线子阵列在θk方向上的导向矢量，其表达式为：其中[·]T表示转置；S22、共采集T个采样快拍，得到第q层稀疏均匀直线子阵列接收信号的采样协方差矩阵其中[·]H表示共轭转置，q∈[1,Q]。所述步骤S3中具体通过前后向空间平滑方法实现重构所述采样协方差矩阵的秩。步骤S4具体包含：S41、对所述平滑协方差矩阵进行特征分解，获得与对应的求根多项式本文档来自技高网...

【技术保护点】
1.一种基于多层互质阵列的相干源波达方向估计方法，用于估计空间内来自K个方向的远场窄带相干信号源的波达方向，其特征在于，包含步骤：/nS1、接收端利用M个阵元构建包含Q层稀疏均匀直线子阵列的多层互质阵列；所述多层互质阵列中的第q层稀疏均匀直线子阵列包含N＝M/Q个阵元，该N个阵元位于同一条直线上，相邻阵元的间距为M

【技术特征摘要】
1.一种基于多层互质阵列的相干源波达方向估计方法，用于估计空间内来自K个方向的远场窄带相干信号源的波达方向，其特征在于，包含步骤：
S1、接收端利用M个阵元构建包含Q层稀疏均匀直线子阵列的多层互质阵列；所述多层互质阵列中的第q层稀疏均匀直线子阵列包含N＝M/Q个阵元，该N个阵元位于同一条直线上，相邻阵元的间距为Mqd；其中q∈[1,Q]，d＝λ/2，λ为所述窄带相干信号源的波长；M1，…，Mq均为正整数且两两互质；Q＞K且
S2、对第q层稀疏均匀直线子阵列的接收信号建模，根据所述建模采集T个采样快拍，得到第q层稀疏均匀直线子阵列接收信号的采样协方差矩阵q∈[1,Q]；
S3、重构所述采样协方差矩阵的秩，获得维度为L×L的平滑协方差矩阵q∈[1,Q]；
S4、对所述平滑协方差矩阵使用求根多重信号分类方法，获得模糊波达方向集合其中为K个信号源各自的角度模糊值，k∈[1,K]；为与对应的角度的模糊值，m∈[1,Mq-1]；q∈[1,Q]；
S5、构建代价函数F(θ)；选取满足为中的K个最小值；将作为第q层稀疏均匀直线子阵列对K个信号源的模糊波达方向估计值；
S6、将作为所述多层互质阵列对K个信号源的波达方向估计值。


2.如权利要求1所述的基于多层互质阵列的相干源波达方向估计方法，其特征在于，所述步骤S1中，不同层的稀疏均匀直线子阵列互相平行；每层稀疏均匀直线子阵列的首个阵元重叠。


3.如权利要求1所述的基于多层互质阵列的相干源波达方向估计方法，其特征在于，所述步骤S2具体包含：
S21、对第q层稀疏均匀直线子阵列的接收信号建模为：<...

【专利技术属性】
技术研发人员：刘爱华，李亚军，王鹏飞，朱剑，王阳阳，成思文，杨娜，
申请(专利权)人：上海无线电设备研究所，
类型：发明
国别省市：上海;31