[发明专利]适用于强干扰环境下的波束域变分贝叶斯方位估计方法

权利要求书

1.适用于强干扰环境下的波束域变分贝叶斯方位估计方法，其特征在于，包括以下步骤：

步骤1：假设K_S个目标信号和K_D个干扰信号分别从和方向上入射至M元的均匀线阵列，其中表示第个目标信号方位，表示第个干扰信号方位，目标信号和干扰信号之间互不相关；当阵列接收到信号后，将接收信号划分为N段，对每一段进行傅里叶变换后，宽带信号划分为L个子带；第l个子带上第n段阵列接收信号对应的傅里叶变换系数记为采样协方差矩阵计算为

上标“H”为共轭转置运算；

步骤2：建立波束域模型，包括以下步骤：

步骤2.1：通过常规波束形成算法确定目标信号所在方位区域Θ_S＝[Θ_SL,Θ_SR]，其中Θ_SL和Θ_SR分别为区域Θ_S的左右边界；将该区域均匀划分为K_B个网格对于第l个子带，该区域上K_B个MVDR-DL波束形成器组成的波束形成矩阵W_l为

其中为指向φ_k的MVDR-DL波束形成器的加权量，a_l(φ_k)为第l个子带上指向φ_k的阵列流形，为通过求解的(M-K_S-K_D)个小特征值的平均值得到的噪声功率估计值，I_M为M维的单位矩阵，上标“-1”表示矩阵求逆；

通过该矩阵对步骤1得到的采样协方差矩阵进行滤波，得到波束域上协方差矩阵为

其中P_l^S和P_l^D表示第l个子带上的目标信号和干扰信号的协方差矩阵，σ_l表示第l个子带上的噪声功率，和表示第l个子带上目标信号和干扰信号的阵列流形矩阵，E_l为第l个子带上的扰动误差；

步骤2.2：对步骤2.1得到的协方差矩阵进行按列向量化运算，得到

其中和分别表示目标信号和干扰信号的功率向量，和分别表示W_l^HW和W_l^HE_lW矩阵按列向量化的向量，⊙表示Khatri-Rao乘积，上标“*”为共轭运算；

步骤2.3：定义矩阵对于第m行第n列的元素[J]_mn，若n＝K_B(m-1)+m，[J]_mn＝1，否则[J]_mn＝0；将步骤2.2得到结果与该矩阵相乘，因MVDR-DL在区域Θ_S外对干扰的波束响应较低，因此波束功率输出的线性关系式表示为：

其中

步骤2.4：将区域Θ_S均匀划分为K_G个网格，网格点组成的向量记为基于该离散网格，步骤2.3中的公式重新表示为

式中为该网格上的阵列流形矩阵，h＝1,...,K_G表示的第h列；p_l是一个稀疏向量，当p_l第m个元素等于第n个元素，否则为0；

步骤3：在贝叶斯框架下进行迭代计算，最终输出DOA估计值，包括以下子步骤：

步骤3.1：构建后验概率密度函数模型，表示为：

式中为联合概率密度函数，N(·)表示实高斯分布，表示Hadamard积；p(p_l；γ)＝N(0,Γ^-1)，为信号稀疏参数，上标“T”为转置运算，Γdiag(γ)为以γ中元素为对角元素的对角矩阵；为噪声功率的方差；为边缘概率密度函数；

构建好后，给定参数迭代初始值：

信号稀疏参数初始值其中表示元素为1的K_G×1维向量；

噪声功率后验均值初始值

噪声功率方差

步骤3.2：进行迭代更新，分别完成对信号功率后验协方差矩阵Σ_l，噪声功率后验方差信号功率后验均值向量μ_l，信号稀疏参数γ，噪声功率后验均值以及噪声功率方差更新；

在第i次迭代中，第l个子带上信号功率后验协方差矩阵更新：

其中上标(i)表示第i次迭代，Γ^(i-1)＝diag(γ^(i-1))；

第l个子带上噪声功率后验方差更新：

第l个子带上信号功率后验均值向量更新：信号稀疏参数更新：

其中表示γ⁽ⁱ⁾第p个元素，为第p个元素，为第p行第p列的元素；第l个子带上噪声功率后验均值更新：

第l个子带上噪声功率方差更新：

更新后的参数值若满足迭代终止条件，所估计的功率谱为

该功率谱中峰值对应的方位即为目标信号的DOA估计值；

若不满足迭代终止条件，则继续进行更新，直至满足迭代终止条件。

2.如权利要求1所述的适用于强干扰环境下的波束域变分贝叶斯方位估计方法，其特征在于，所述步骤3.2中的迭代终止条件为：当迭代满足||γ⁽ⁱ⁾-γ^(i-1)||₂/||γ^(i-1)||₂≤10^-3，其中||·||₂表示l₂范数，或者迭代次数大于Iter_max＝1000时，迭代终止。