一种基于离格稀疏贝叶斯学习的nested阵列波达方向估计方法技术

本发明专利技术公开了一种基于离格稀疏贝叶斯学习的nested阵列波达方向估计方法，包括以下步骤：1：nested阵列接收到的窄带高斯信号经过匹配滤波后，得到在t时刻包含DOA信息的数据向量x(t)。2：利用x(t)求得在T快拍数下的接收数据协方差矩阵

【技术实现步骤摘要】
一种基于离格稀疏贝叶斯学习的nested阵列波达方向估计方法
本专利技术属于阵列信号处理领域，涉及阵列信号的波达方向估计，具体地说是一种适用于非均匀的nested阵列信号的波达方向估计的方法。
技术介绍
采用天线阵列进行多窄带信号的波达方向(DirectionofArrival,DOA)估计的方法已经被广泛应用于雷达、声纳和通信等领域。在过去的几十年里，人们提出了大量行之有效的DOA估计方法。由于非均匀阵列自身在增加信号处理自由度方面具有很大优势，目前，基于nested阵列的目标角度估计算法也成为研究热点。例如在文献：P.Pal,P.Vaidyanathan,Nestedarrays:Anovelapproachtoarrayprocessingwithenhanceddegreesoffreedom,IEEETransactionsonSignalProcessing,58(8)(2010)4167-4181中，提出了一种基于空间平滑的多重信号分类(SpatialSmoothingbasedMultipleSignalClassification，SS-MUSIC)方法；在文献：J.Yang,G.Liao,J.Li,Anefficientoff-griddoaestimationapproachfornestedarraysignalprocessingbyusingsparsebayesianlearningstrategies,SignalProcessing,128(2016)110-122中，提出了一种基于二阶泰勒展开的稀疏贝叶斯学习(SparseBayesianLearning,SBL)方法。然而，传统的基于子空间处理的算法易受信噪比和快拍数的影响，现有的基于稀疏贝叶斯学习的DOA估计方法主要缺点在于：有效阵列孔径损失较严重，计算复杂度较高。
技术实现思路
针对现有方法的不足，本专利技术提出了一种新型的基于离格(off-grid)稀疏贝叶斯学习的nested阵列DOA估计方法，该方法可以有效避免阵列孔径的损失，简化稀疏贝叶斯学习估计的步骤，并采用一种新的网格更新方法来减小由off-grid模型引起的误差，可明显改善DOA的估计性能。用于实现本专利技术的技术解决方案包括如下步骤：步骤1：nested阵列接收到的远场窄带高斯信号经过匹配滤波后，得到在t时刻包含DOA信息的数据向量x(t)。步骤2：利用步骤1中接收到的数据向量x(t)，求得在T快拍数下的接收数据协方差矩阵对向量化，得到一个一维的数据向量步骤3：在的范围均匀划分出个网格点设置迭代次数计数变量l＝1，初始化方差向量δ和角度偏移向量β，构造测量矩阵步骤4：利用期望最大化(ExpectationMaximization，EM)准则，更新维的方差向量δ和角度偏离值β。步骤5：利用步骤4中求出的β值更新网格如果在的范围中，则更新网格点否则不更新。步骤6：判断迭代计数变量l是否达到上限L或δ是否收敛，如果都不满足，则迭代计数变量l＝l+1，并令β＝0，利用更新的网格更新并返回步骤4。步骤7：对方差向量δ进行谱峰搜索，得到K个极大值点对应的角度，即为目标角度的最终估计值。本专利技术的有益效果：本专利技术提出了一种基于off-grid稀疏贝叶斯学习的nested阵列系统的DOA估计方法，有效的避免阵列孔径减小和复杂的矩阵变换，能够自动求出噪声方差，新的网格更新方法可以获得更加精确的角度估计值。附图说明图1是本专利技术实施流程图。图2是200次蒙特卡洛实验条件下，nested阵列内外层阵元数均为3，快拍数为100，网格间距为5°，信噪比由-10dB到10dB变化，检测2个目标时本专利技术与二阶SBL方法估计DOA的均方根误差(RMSE)比较。图3是在信噪比为0dB，快拍数由50到800变化，检测2个目标时本专利技术与二阶SBL方法估计DOA的均方根误差比较。具体实施方式下面结合附图和具体实施例对本专利技术作进一步说明。如图1所示，本专利技术的方法包括如下步骤：(1)nested阵列接收到的远场窄带高斯信号经过匹配滤波后，得到在t时刻包含DOA信息的数据向量x(t)＝As(t)+n(t)，t＝1,2,...,T，式中：T表示快拍数；s(t)＝[s1(t),s2(t),...,sK(t)]T表示在t时刻发射的K个不相关窄带信号，其中sk(t)满足均值为0，方差为k＝1,2,...,K的复高斯分布，(·)T表示转置；A＝[a(θ1),a(θ2),...,a(θK)]表示M×K维的阵列流型矩阵，其中M＝M1+M2为nested阵列阵元个数，M1和M2分别表示nested阵列内外层阵元个数，内外层阵元间距分别为d和(M1+1)d，令[r1,r2,...,rM]＝[0,1,...,(M1-1),M1,2(M1+1)-1,M2(M1+1)-1]，则第m个阵元的位置可以表示为d·rm，m＝1,2,...,M。阵列流型向量a(θk)＝[1,exp(j2πdr2sinθk/λ),...,exp(j2πdrMsinθk/λ)]T，k＝1,2,...,K，θk为第k个真实的DOA，λ为电磁波的波长；n(t)表示t时刻一个M维的均值为0，方差为的高斯白噪声。(2)在T快拍数下，求数据向量x(t)的协方差矩阵(·)H表示共轭转置，将向量化得到vec(·)表示向量化操作。(3)在的范围均匀划分出个网格点在此网格上的虚拟阵列流形矩阵为(·)*表示共轭运算，表示Kronecker积。一阶泰勒展开的Off-grid虚拟阵列流型矩阵为：其中(·)′表示一阶导数运算，diag(·)表示取对角运算，βi表示网格点上的角度偏移值。进一步，将步骤(2)中的数据模型表示为：其中表示T次快拍数下的协方差矩阵与真实值Rx之间的近似误差，满足均值为零，方差为的复高斯分布，向量em表示除第m个元素为1，其余元素均为零，维向量的非零元素对应于发射信号方差设置迭代次数计数变量l＝1，初始化方差向量中的各元素为1，同时初始化β为全零向量(β＝0)，定义测量矩阵为(4)利用期望最大化(ExpectationMaximization，EM)准则，更新δ和β：β＝P-1v其中Ξ＝μμH+Σ，[·]ii表示矩阵的第(i,i)个元素，ρ为一个较小正数(例如ρ＝0.01)，μ和Σ分别表示d后验分布的均值和方差，μ-表示向量μ的前个元素，μ0表示μ的最后一个元素，BW＝W-1/2B，1W＝W-1/21n，Δ＝diag(δ)。其中Σ(l1:l2,c1:c2)表示Σ的第l1到l2行和c1到c2列组成的子矩阵，⊙表示Khatri-Rao积，Re{·}表示取实部操作。(5)将网格看作可变参数，利用步骤(4)中求出的β值更新网格如果在的范围中，则更新网格点否则不更新。(6)判断迭代计数变量l是否达到上限L或方差δ是否收敛，如果都不满足，则迭代计数变量l＝l+1，并令β＝0，然后利用更新的网格更新测量矩阵并返回(4)。(7)对方差向量δ进行谱峰搜索，得到K个极大值点对应的角度，即为目标角度的最终估计值。下面结合仿真实验对本专利技术的效果做进一步说明。为了评估本方法的性能，考虑一nested阵列，阵元个数M＝6，其中内外层阵元数M1＝M2＝3，网格间距为5°，假设远场有本文档来自技高网...

【技术保护点】
1.一种基于离格稀疏贝叶斯学习的nested阵列波达方向估计方法，其特征在于，包括如下步骤：步骤1：nested阵列接收到的远场窄带高斯信号经过匹配滤波后，得到在t时刻包含DOA信息的数据向量x(t)；步骤2：利用步骤1中接收到的数据向量x(t)，求得在T快拍数下的接收数据协方差矩阵

【技术特征摘要】
1.一种基于离格稀疏贝叶斯学习的nested阵列波达方向估计方法，其特征在于，包括如下步骤：步骤1：nested阵列接收到的远场窄带高斯信号经过匹配滤波后，得到在t时刻包含DOA信息的数据向量x(t)；步骤2：利用步骤1中接收到的数据向量x(t)，求得在T快拍数下的接收数据协方差矩阵对向量化，得到一个一维的数据向量步骤3：在的范围均匀划分出个网格点设置迭代次数计数变量l＝1，初始化方差向量δ和角度偏移向量β，构造测量矩阵步骤4：利用期望最大化(ExpectationMaximization，EM)准则，更新维的方差向量δ和角度偏离值β；步骤5：利用步骤4中求出的β值更新网格步骤6：判断迭代计数变量l是否达到上限L或δ是否收敛，如果都不满足，则迭代计数变量l＝l+1，并令β＝0，利用更新的网格更新并返回步骤4；步骤7：对方差向量δ进行谱峰搜索，得到K个极大值点对应的角度，即为目标角度的最终估计值。2.根据权利要求1所述的一种基于离格稀疏贝叶斯学习的nested阵列波达方向估计方法，其特征在于，所述步骤1中数据向量x(t)的表达式为：x(t)＝As(t)+n(t)，t＝1,2,...,T，式中：T表示快拍数；s(t)＝[s1(t),s2(t),...,sK(t)]T表示在t时刻发射的K个不相关窄带信号，其中sk(t)满足均值为0，方差为k＝1,2,...,K的复高斯分布，(·)T表示转置；A＝[a(θ1),a(θ2),...,a(θK)]表示M×K维的阵列流型矩阵，其中M＝M1+M2为nested阵列阵元个数，M1和M2分别表示nested阵列内外层阵元个数，内外层阵元间距分别为d和(M1+1)d，令[r1,r2,...,rM]＝[0,1,...,(M1-1),M1,2(M1+1)-1,M2(M1+1)-1]，则第m个阵元的位置可以表示为d·rm，m＝1,2,...,M。阵列流型向量a(θk)＝[1,exp(j2πdr2sinθk/λ),...,exp(j2πdrMsinθk/λ)]T，k＝1,2,...,K，θk为第k个真实的DOA，λ为电磁波的波长；n(t)表示t时刻一个M维的均值为0，方差为的高斯白噪声。3.根据权利要求1所述的一种基于离格稀疏贝叶斯学习的nested阵列波达方向估计方法，其特征在于，所述步骤2中接收数据的协方差矩阵的表...

【专利技术属性】
技术研发人员：戴继生，陈方方，鲍煦，张文策，
申请(专利权)人：江苏大学，
类型：发明
国别省市：江苏,32