一种基于离格稀疏贝叶斯学习的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°，假设远场有2个相互独立的目标，分别随机取自范围[-30°,-20°]和[0°,10°]。实验为检测2个目标时，本发明与二阶sbl方法估计doa的均方根误差比较。在所有的实验中，背景噪声假设为高斯白噪声，蒙特卡洛实验200次。

实验条件

实验1，采用本发明和现有的二阶sbl方法，在快拍数t＝100，信噪比从-10db到10db变化时，目标角度估计均方根误差(rmse)的比较实验，仿真结果如图2所示。

实验2，是在信噪比为0db，快拍数由50到800变化时，本发明与二阶sbl方法的目标角度估计rmse比较，仿真结果如图3所示。

实验分析

从图2可以看出，本发明估计目标角度的rmse随信噪比的增加而降低，与二阶sbl方法相比，本发明能够精确地估计出目标doa，即使在低信噪比下也能有很好的性能。

从图3可以看出，采用本发明估计目标角度的rmse随快拍数的增加而降低，与二阶sbl方法相比，本发明能以更高的精度对目标doa的估计。

上文所列出的一系列的详细说明仅仅是针对本发明的可行性实施方式的具体说明，它们并非用以限制本发明的保护范围，凡未脱离本发明技艺精神所作的等效实施方式或变更均应包含在本发明的保护范围之内。