[發明專利]一種基于迭代重賦權的塊稀疏信號重構方法

在審

• 申請號: 202010583891.3
• 申請日: 2020-06-24
• 公開（公告）號: CN111900994A
• 公開（公告）日: 2020-11-06
• 發明（設計）人: 呂明久;陳文峰;趙麗;趙欣;楊軍
• 申請（專利權）人: 中國人民解放軍空軍預警學院雷達士官學校
• 主分類號: H03M7/30
• 分類號: H03M7/30
• 代理公司: 北京翔石知識產權代理事務所(普通合伙) 11816
• 代理人: 李勇

1.一種基于迭代重賦權的塊稀疏信號重構方法，其特征在于，包括：

步驟S1：進入外循環并初始化；

步驟S2：利用塊稀疏解的相鄰系數計算權值；

步驟S3：進入內循環并初始化；

步驟S4：利用復數近似消息傳遞算法求解1范數最小化問題的稀疏解；

步驟S5：判斷是否滿足內循環停止條件，若滿足，執行步驟S6；若不滿足，令t＝t+1，執行步驟S4；

步驟S6：利用權值和1范數最小化問題的稀疏解更新塊稀疏解；

步驟S7：判斷是否滿足外循環停止條件，若滿足，執行步驟S8；若不滿足，令k＝k+1，執行步驟S2；

步驟S8：輸出塊稀疏解。

2.根據權利要求1所述的基于迭代重賦權的塊稀疏信號重構方法，其特征在于，所述塊稀疏解問題表示公式為：

y＝As+w

其中，為量測值，為量測矩陣，為塊稀疏信號，w為零均值高斯白噪聲。

3.根據權利要求1所述的基于迭代重賦權的塊稀疏信號重構方法，其特征在于，所述步驟S1進入外循環并初始化，其中，輸入量測值y，感知矩陣A，令θ⁰＝1，β＝1，k＝1，k_max＝10，ε_k＝10^-5，ε＝10^-16。

4.根據權利要求1所述的基于迭代重賦權的塊稀疏信號重構方法，其特征在于，所述步驟S2，利用所述塊稀疏解的相鄰系數計算所述權值，計算所述權值公式如下：

5.根據權利要求1所述的基于迭代重賦權的塊稀疏信號重構方法，其特征在于，所述步驟S3進入內循環并初始化，令z⁰＝y，t＝1，t_max＝200，ε_t＝10^-5。

6.根據權利要求1所述的基于迭代重賦權的塊稀疏信號重構方法，其特征在于，步驟S4利用復數近似消息傳遞算法求解1范數最小化問題的稀疏解，利用如下公式計算中間變量v^t

利用如下公式計算殘差z^t

其中，是η^Re對實部的偏導數，是η^Im對虛部的偏導數，·表示對一個向量求平均；

利用式如下公式計算所述稀疏解

其中，所述復數軟閾值函數為：

η^Re和η^Im分別表示η(u+iv；μ)的實部和虛部。

7.根據權利要求1所述的基于迭代重賦權的塊稀疏信號重構方法，其特征在于，所述步驟S5判斷是否滿足內循環停止條件，即t＞t_max或若滿足，執行所述步驟S6；若不滿足，令t＝t+1，執行所述步驟S4。

8.根據權利要求1所述的基于迭代重賦權的塊稀疏信號重構方法，其特征在于，所述步驟S6利用權值和1范數最小化問題的稀疏解更新所述塊稀疏解，令利用如下公式計算塊稀疏解θ^k

9.根據權利要求1所述的基于迭代重賦權的塊稀疏信號重構方法，其特征在于，所述步驟S7判斷是否滿足外循環停止條件，即k＞k_max或||θ^k-θ^k-1||₂＜ε_k||θ^k-1||₂，若滿足，執行所述步驟S8；若不滿足，令k＝k+1，執行所述步驟S2。

10.根據權利要求1所述的基于迭代重賦權的塊稀疏信號重構方法，其特征在于，所述步驟S8輸出塊稀疏解θ^k。