[发明专利]杂波背景下天线阵列的信号波达方向估计方法

摘要

本发明公开了一种杂波背景下天线阵列的信号波达方向估计方法。其特征在于，包括：接收数个快拍向量，构造空间抽样阵列流型矩阵及其求导矩阵；初始化稀疏支撑向量、正则化参数、杂波‑噪声联合功率向量与波达方向偏差向量；依次迭代更新信号后验均值向量、协方差矩阵以及上述各变量；迭代停止条件满足时，波达方向估计由稀疏支撑向量与波达方向偏差向量给出。本发明的有益效果：(1)本发明利用杂波与噪声污染的各快拍中稀疏表示信号之间的联合稀疏性，对问题进行建模并求解，从而利用较少快拍即可得到较精确的估计；(2)本发明的正则化参数在迭代过程中自动更新，不需要经验预设；(3)本发明估计了各快拍中时变的杂波‑噪声联合功率。

主权项

1.一种杂波背景下天线阵列的信号波达方向估计方法，其特征在于，包括：步骤1.接收天线阵列的M个天线单元采集的T个快拍，构成T个快拍数据向量x_t,t＝1,2,...,T；设置一个对空间角度全覆盖的均匀抽样角度网格向量并由此构造阵列流型矩阵与阵列流型求导矩阵步骤2.迭代程序参数a与b的值设置为趋近于0的某个正实数，对迭代程序中的各变量初始化：稀疏支撑向量α初始化为正则化参数μ初始化为杂波‑噪声联合功率向量γ初始化为波达方向偏差向量β初始化为步骤3.利用T个快拍数据向量x_t,t＝1,2,...,T、阵列流型矩阵阵列流型求导矩阵旧的稀疏支撑向量旧的杂波‑噪声联合功率向量和旧的波达方向偏差向量计算各快拍中稀疏表示信号向量的后验均值向量ζ_t,t＝1,2,...,T与后验协方差矩阵Λ_t,t＝1,2,...,T；步骤4.利用旧的正则化参数各快拍中稀疏表示信号向量的后验均值向量ζ_t,t＝1,2,...,T与后验协方差矩阵Λ_t,t＝1,2,...,T，更新稀疏支撑向量α；步骤5.利用旧的稀疏支撑向量更新正则化参数μ；步骤6.利用T个快拍数据向量x_t,t＝1,2,...,T、阵列流型矩阵阵列流型求导矩阵旧的波达方向偏差向量各快拍中稀疏表示信号向量的后验均值向量ζ_t,t＝1,2,...,T与后验协方差矩阵Λ_t,t＝1,2,...,T，更新杂波‑噪声联合功率向量γ；步骤7.在稀疏支撑向量α上搜索出K个峰值，其峰值位置索引构成η；利用α的峰值位置索引η、T个快拍数据向量x_t,t＝1,2,...,T、阵列流型矩阵阵列流型求导矩阵旧的波达方向偏差向量各快拍中稀疏表示信号向量的后验均值向量ζ_t,t＝1,2,...,T与后验协方差矩阵Λ_t,t＝1,2,...,T，更新波达方向偏差向量β；步骤8.计算是否满足迭代终止条件：若未满足，令回到步骤3继续执行迭代；否则结束迭代并基于当前的稀疏支撑向量α与波达方向偏差向量β，给出波达方向估计结果。