基于广义正交匹配追踪的波达方向估计方法技术

本发明专利技术提供的是一种基于广义正交匹配追踪的波达方向估计方法。将若干信号入射到均匀阵列接收模型上，得到均匀阵列接收信号，通过等角度空间网格划分方式将阵列流型矩阵转换为稀疏矩阵；生成一个高斯随机测量矩阵，进行SVD矩阵分解，并对得到的奇异值在合适的范围重新赋值，获得测量矩阵；使用测量矩阵对稀疏矩阵进行采样得到观测矩阵，使用测量矩阵对阵列接收信号进行采样，得到采样后的阵列接收信号；将采样后的阵列接收信号及观测矩阵作为输入，采用广义正交匹配追踪算法进行信号的重构，通过多次迭代使估计信号逐渐逼近原信号，至方差到达阈值时停止迭代，估计出来波信号的入射角度。本发明专利技术能够缩短信号重构时间，提高DOA估计的成功率。

【技术实现步骤摘要】
基于广义正交匹配追踪的波达方向估计方法
本专利技术涉及的是一种阵列信号处理方法，具体涉及一种波达方向估计方法。
技术介绍
波达方向估计(directionofarrival，DOA)是阵列信号处理中的重要研究内容之一，被广泛应用于声呐、雷达、通信、射电天文、地震勘探、医学诊疗等众多经济领域中。自20世纪60年代以来，研究者们提出了大量有效的DOA估计算法，主要有Capon提出的最小方差谱估计法(minimumvariancedistortionlessresponse，MVDR)和以Schimidt提出的多重信号分类法(multiplesignalclassification，MUSIC)为代表的子空间算法。但是这些算法受限于奈奎斯特采样定理，在处理快拍数少、信号相干等问题时性能不太好并且存在处理时间长等问题。近年来，信号处理领域中提出压缩感知(CompressingSensing，CS)理论，已广泛应用于图像处理和无线通信等诸多领域。压缩感知在DOA估计的应用中，学者们也是致力于信号稀疏化、测量矩阵以及稀疏重构算法的研究中。在重构算法一定时，测量矩阵的选择影响着重构的成功率，测量矩阵的列独立性、测量矩阵与稀疏基的非相干性等均会对重构信号的质量产生影响。目前，学者们也提出了众多性能优越的测量矩阵，但是这些测量矩阵存在相关性比较高等缺陷，导致重构时误差较大。常见的稀疏重构算法有正交匹配追踪算法(OrthogonalMatchingPursuit，OMP)，压缩采样匹配追踪算法(CompressiveSamplingMP，CoSaMP)以及子空间追踪算法(Subspacepursuit，SP)等。这些方法虽然能够在一定程度上将数据恢复出来，但是存在恢复时间长、成功率比较低等缺陷。
技术实现思路
本专利技术的目的在于提供一种能够缩短信号重构时间，提高DOA估计的成功率的基于广义正交匹配追踪的波达方向估计方法。本专利技术的目的是这样实现的：步骤一：将若干信号入射到均匀阵列接收模型上，得到均匀阵列接收信号，通过等角度空间网格划分方式将阵列流型矩阵转换为稀疏矩阵；步骤二：生成一个高斯随机测量矩阵，对高斯随机测量矩阵进行SVD矩阵分解，并对得到的奇异值在合适的范围重新赋值，降低高斯随机测量矩阵的相关性，获得测量矩阵；步骤三：使用测量矩阵对稀疏矩阵进行采样得到观测矩阵，使用测量矩阵对阵列接收信号进行采样，得到采样后的阵列接收信号；步骤四：将采样后的阵列接收信号及观测矩阵作为输入，采用广义正交匹配追踪算法进行信号的重构，通过多次迭代使估计信号逐渐逼近原信号，至方差到达阈值时停止迭代，估计出来波信号的入射角度。针对已有技术中存在的恢复时间长、成功率比较低等缺陷，本专利技术提出了一种基于广义正交匹配追踪的DOA估计方法，首先利用矩阵的SVD分解来降低测量矩阵的相关性，增强测量矩阵的独立性，然后采用广义正交匹配追踪算法进行信号的重构，最终估计出来波信号的角度。本专利技术所提方法缩短了信号重构时间，提高了DOA估计的成功率。与现有技术相比，本专利技术的优点在于：(1)在稀疏化处理过程中，对测量矩阵进行SVD分解，对奇异值在合适范围内重新赋值，且将其奇异值平均化，可以增大测量矩阵的最小奇异值，降低测量矩阵的相关性，在进行重构时能够获得更好的效果；(2)采用广义正交匹配追踪算法进行信号重构，从而能够实现准确的DOA估计，相较于现有技术而言，每次迭代时选出多个原子进行数据重构，在一定程度上避免了单个原子重构时造成的误差，加快了算法的收敛，从而降低了算法的复杂度，最终提高了重构的成功率，缩短了重构所需时间。附图说明图1为均匀线形阵列接收模型；图2为等角度空间网格划分图；图3为现有技术与本专利技术在不同信噪比下重构成功率的对比图；图4为现有技术与本专利技术在不同信噪比下均方误差的对比图。具体实施方式本专利技术的基于广义正交匹配追踪的波达方向估计方法主要包括如下步骤：步骤一：构造均匀线性阵列流型矩阵，假设若干信号入射到均匀阵列接收模型上，得到均匀阵列接收信号，通过等角度空间网格划分方式将阵列流型矩阵转换为稀疏矩阵，得到稀疏条件下DOA估计数学模型；步骤二：生成一个高斯随机测量矩阵，对测量矩阵进行SVD矩阵分解，并对得到的奇异值在合适的范围重新赋值，降低测量矩阵的相关性，获得性能更好的测量矩阵；步骤三：使用测量矩阵对稀疏矩阵进行采样得到观测矩阵，使用测量矩阵对阵列接收信号进行采样，得到采样后的阵列接收信号；步骤四：将采样后的阵列接收信号及观测矩阵作为输入，采用广义正交匹配追踪算法进行信号的重构，通过多次迭代使估计信号逐渐逼近原信号，至方差到达阈值时停止迭代，估计出来波信号的入射角度。下面举例对本专利技术做更详细的描述。步骤一：构造均匀线性阵列流型矩阵，假设若干信号入射到均匀阵列接收模型上，得到均匀阵列接收信号，通过等角度空间网格划分方式将得到稀疏化后的阵列流型矩阵，即稀疏变换矩阵，得到稀疏条件下DOA估计数学模型。图1为均匀线形阵列接收模型，假设K个远场窄带信号入射到M个各向同性的均匀线阵上，阵元间距为d，各信号来波方向为θi(i＝1,2,…,K)。则阵列接收信号为：式中，为M×1维的阵列接收信号矢量，s(t)＝[s1(t),s2(t),…,sK(t)]T为K×1维的空间信号矢量，为M×K维的阵列流型矩阵，为M×1维的方向矢量，θk为第k个信号的来波方向，n(t)＝[n1(t),n2(t),…,nM(t)]T为M×1维的噪声矢量。图2为等角度空间网格划分示意图，按照等角度划分方式进行空间网格划分，空间信号是N＞＞K可稀疏的，通过等角度网格划分可将空间信号稀疏化。将空间均匀划分为N份{θ1,θ2,…,θN}，N＞＞K。假设每一个θn(n＝1,2,…,N)都对应一个空间信号xn(n＝1,2,…,N)，这样便构造出了稀疏度为K的N×1维稀疏空间信号矢量x(t)＝[x1(t),x2(t),…,xN(t)]T，x(t)中只有K个位置的元素是非零的，对应着空间中实际存在的信号s(t)，其余N-K个位置的元素为零。则稀疏化后的空间信号矢量对应的M×N维阵列流型矩阵为：则稀疏表示模型下的DOA估计数学模型为：式中，为远场窄带信号入射到阵列上的接收信号，为稀疏变换矩阵，x(t)为待求稀疏系数矢量，n(t)为阵列噪声矢量。步骤二：生成一个高斯随机测量矩阵Φ，对测量矩阵Φ进行SVD矩阵分解，并对得到的奇异值在合适的范围重新赋值，降低测量矩阵的相关性，获得性能更好的测量矩阵。对测量矩阵Φ进行SVD分解，可得：Φ＝UΣVT(4)式中，U∈RM×N和V∈RN×N为正交矩阵，Σ∈RM×N为对角阵，Σ中的对角线上元素为矩阵Φ的奇异值。Σ可表示为：式中，Δ＝diag(δ1,δ2,…,δr)，且δ1＞δ2＞…＞δr。对得到的奇异值在合适的范围内重新赋值，所以构造一个新的对角矩阵Δ*＝diag(a1,a2,…,aM)，其中δ1＞1＞a1＞a2＞…＞aM＞0.5＞δr，可得：构造出新的矩阵对矩阵进行SVD分解且平均化其奇异值后，可降低矩阵的相关性，因此对矩阵的特征值取平均，即：Φnew＝UΣnewVT(10)Φnew即为经过矩阵分解优化后的新测量矩阵。步骤三：使用测量矩阵对稀疏矩阵进行采样得到观测矩阵，使用测本文档来自技高网...

【技术保护点】
1.一种基于广义正交匹配追踪的波达方向估计方法，其特征是：步骤一：将若干信号入射到均匀阵列接收模型上，得到均匀阵列接收信号，通过等角度空间网格划分方式将阵列流型矩阵转换为稀疏矩阵；步骤二：生成一个高斯随机测量矩阵，对高斯随机测量矩阵进行SVD矩阵分解，并对得到的奇异值在合适的范围重新赋值，获得测量矩阵；步骤三：使用测量矩阵对稀疏矩阵进行采样得到观测矩阵，使用测量矩阵对阵列接收信号进行采样，得到采样后的阵列接收信号；步骤四：将采样后的阵列接收信号及观测矩阵作为输入，采用广义正交匹配追踪算法进行信号的重构，通过多次迭代使估计信号逐渐逼近原信号，至方差到达阈值时停止迭代，估计出来波信号的入射角度。

【技术特征摘要】
1.一种基于广义正交匹配追踪的波达方向估计方法，其特征是：步骤一：将若干信号入射到均匀阵列接收模型上，得到均匀阵列接收信号，通过等角度空间网格划分方式将阵列流型矩阵转换为稀疏矩阵；步骤二：生成一个高斯随机测量矩阵，对高斯随机测量矩阵进行SVD矩阵分解，并对得到的奇异值在合适的范围重新赋值，获得测量矩阵；步骤三：使用测量矩阵对稀疏矩阵进行采样得到观测矩阵，使用测量矩阵对阵列接收信号进行采样，得到采样后的阵列接收信号；步骤四：将采样后的阵列接收信号及观测矩阵作为输入，采用广义正交匹配追踪算法进行信号的重构，通过多次迭代使估计信号逐渐逼近原信号，至方差到达阈值时停止迭代，估计出来波信号的入射角度。2.根据权利要求1所述的基于广义正交匹配追踪的波达方向估计方法，其特征是所述对高斯随机测量矩阵进行SVD矩阵分解具体包括：步骤二：对高斯随机测量矩阵Φ进行SVD分解，得：Φ＝UΣVT其中，U∈RM×N和V∈RN×N为正交矩阵，Σ∈RM×N为对角阵，M为各向同性的均匀线阵数，N为将空间均匀划分的份数，Σ中的对角线上元素为矩阵Φ的奇异值，Σ表示为：式中，Δ＝diag(δ1,δ2,…,δr)，且δ1＞δ2＞…＞δr。3.根据权利要求2所述的基于广义正交匹配追踪的波达方向估计方法，其特征是所述对得到的奇异值在合适的范围内重新赋值具体包括：构造一个新的对角矩...

【专利技术属性】
技术研发人员：郜丽鹏，刁弘扬，高敬鹏，赵忠凯，蒋伊琳，
申请(专利权)人：哈尔滨工程大学，
类型：发明
国别省市：黑龙江,23