专利名称：一种块稀疏度未知的稀疏信号压缩感知重构方法
技术领域：
本发明涉及压缩感知技术领域，具体涉及对块稀疏信号的重构方法。
背景技术：
传统的信号采样理论是基于奈奎斯特采样定理，即为保证不丢失源信号的信息， 实现无失真恢复源信号，采样率至少需要两倍信号带宽。这对于宽带模拟信号的数字化往 往需要很高的采样率，增加了物理器件的负担。且对于数据量很大的信号，存储容量和处理 速度都受到进一步限制。压缩感知(Compressed Sensing，CS)是2004年提出的一个全新的信号采样理论， 其思想是对稀疏信号能以远低于奈奎斯特采样率的速度对信号进行全局观测，然后通过适 当的重构算法从观测信号中重构出源信号。CS理论将传统的对信号的采样转化成了对信息 的采样，把采样和压缩结合成一步对信号进行编码，极大地降低信号的采样频率及数据存 储和传输代价，显著地降低信号处理时间和计算成本。由于其特殊的性质，CS理论在信源 编码、数据挖掘、雷达信号处理、医学信号处理以及图像处理等领域有广泛的应用前景。目前块稀疏信号的重构方法主要有混合厶/Λ优化(Mixed IJI1 Optimization Program, L-0PT)算法，±夬稀疏匹配追踪(Block-sparse Matching Pursuit, BMP)算法，土夬 稀疏正交匹配追踪(Block-sparse Orthogonal Matching Pursuit, Β0ΜΡ)算法。L-0PT 算 法是基于混合型范数的凸优化算法，其复杂度较高；BMP算法和BOMP算法一旦找到匹配原 子后便不再改变，所以容易造成过匹配现象，且这些算法都需要块稀疏度f作为先验知识， 而在实际中却很难得到块稀疏度
发明内容
本发明为了解决现有的块稀疏信号匹配追踪类重构方法存在过匹配现象，并且在 未知块稀疏度^时无法重构的问题，提出一种块稀疏度未知的稀疏信号压缩感知重构方 法。本发明方法的具体过程如下
采集获得源信号X的观测信号为y，观测信号y表示为/=LF1, _F2，^3……，/J，其中 为观测信号y的长度，
步骤一、稀疏信号压缩感知重构过程初始化设定初始化块稀疏度k,
l<k<K I为源信号X的块稀疏度，Xeiami ,初始化测量矩阵ΦEITwa『，设定迭代误 差err，分块向量Group E Jt M ,分块向量Group如下表示
其中，J为分块向量Group的子块长度，设定残差初始值r(1=y，恢复矩阵％ = 0 ,步长
对印=1，信号支撑集大小5·=1源信号χ的重构向量^ = 0 ；
步骤二、计算获得与残差I^1最匹配的子空间i7，其中I^1表示第次迭代产生的残 差，i7表示第1次迭代的最匹配的子空间J,具体过程为
Γ[Φμ ))
Ii = Hrgmp^an,
即测量矩阵Φ每一块的转置与上一次迭代的残差IV1进行乘积操作后，对每一块的d 个数据分别取绝对值后，再求获得的J个绝对值的平均值，最后从获得的#个平均值中选择 最大的S个值的标号赋值给所述最匹配的子空间i7，的值对应分块向量Group中的分组 号，所述分组号为1，2，···，#;
步骤三、将步骤二获得的最匹配的子空间i7中对应的测量矩阵Φ的列向量组成的矩
阵Oii与上一次迭代获得的恢复矩阵Tw的并集赋值给初始恢复矩阵 ；
Tj “ Tj^i U Φ恕；
步骤四、计算获得第ι次迭代与观测信号1最匹配的子空间t7， 体过程为
arg max (mean i
φ:[
τ；
其中矩阵中7为无的分块标号，Φ《=(Φ^Φ^的伪
逆矩阵，即矩阵φ &的每一块与观测信号y进行乘积操作后，对每一块的i/个数据分别取
绝对值后，再将获得的J个绝对值求平均值，最后从获得的多个平均值中选择最大的S个值 的标号赋值给所述最匹配的子空间t7，t7的值对应分块向量Group中的分组号，所述分组号 为 1，2，···,#；
步骤五、将步骤四获得的最匹配的子空间t7中多个分组号所对应的测量矩阵Φ中的 列向量组成的矩阵,赋值给恢复矩阵Ti
X -屯、；
步骤六、计算本次迭代的残差r7
Γ -；
其中Φ·^为的伪逆矩阵；步骤七、判断Ihll2 ^ lh-ι Il2成立与否，其中Ihii2表示2-范数；
如果成立，则令对印=step^l,S= StepX义然后执行步骤八；如果不成立，则直接执行 步骤八；
步骤八、令/4+1，返回执行步骤二，直到迭代次数1大于分组数见或者残差小于算法 迭代误差err时，迭代结束，执行步骤九；
步骤九、输出源信号χ的重构向量;^，重构向量&为
完成对源信号χ的压缩感知重构。本发明通过初始化块稀疏度t对每一个块稀疏信号的迭代，找到信号支撑集的一 个子集，随着不断迭代，块稀疏度随之增加，最后找到整个源信号X的支撑集，从而达到重 构源信号X的目的，本发明采用多次迭代并修正支撑集的思想来重构信号，重构的精度高， 与现有的块稀疏度匹配追踪方法以及正交匹配追踪方法相比，对块稀疏信号的重构概率 大，不存在过匹配现象。且本发明不需要以块稀疏度作为先验知识，尤其适用于块稀疏度未 知信号重构领域。
图1为压缩感知基本模型的示意图。图2为观测信号幅度为高斯分布信号的重构 概率对比图。图3为观测信号幅度为二值信号的重构概率对比图。
具体实施例方式
具体实施方式
一、结合图1说明本实施方式，一种块稀疏度未知的稀疏信 号压缩感知重构方法，具体过程如下
采集获得源信号X的观测信号为y，观测信号y表示为/=LF1, _F2，^3……，/J，其中 为观测信号y的长度，
步骤一、稀疏信号压缩感知重构过程初始化设定初始化块稀疏度H Κ^Γ , f为
源信号X的块稀疏度，X e Rm ,初始化测量矩阵Φ e mrN ,设定迭代误差err，分块向
量Group E麗1，分块向量Group如下表示
Group' = [1,...,1,2,...,2,...,M5....,Μ]
d d d
其中，J为分块向量Group的子块长度，设定残差初始值r(l=y，恢复矩阵Tq = 0 ,步长
对卬=1，信号支撑集大小源信号χ的重构向量;-0 ；
步骤二、计算获得与残差I^1最匹配的子空间i7，其中Iv1表示第次迭代产生的残 差，i7表示第1次迭代的最匹配的子空间，i e{1,2,...,AO，具体过程为
即测量矩阵Φ每一块的转置与上一次迭代的残差IV1进行乘积操作后，对每一块的d 个数据分别取绝对值后，再求获得的J个绝对值的平均值，最后从获得的#个平均值中选择 最大的S个值的标号赋值给所述最匹配的子空间i7，i7的值对应分块向量Group中的分组 号，所述分组号为1，2，···，#;
步骤三、将步骤二获得的最匹配的子空间i7中对应的测量矩阵Φ的列向量组成的矩
阵1与上一次迭代获得的恢复矩阵Tp1的并集赋值给初始恢复矩阵 ；
X 二 T1L.
步骤四、计算获得第1次迭代与观测信号y最匹配的子空间t7，/ e {13 2, 体过程为
,具
其中矩阵Φ中j·为χ的分块标号
%的伪
逆矩阵，即矩阵Φ^的每一块与观测信号y进行乘积操作后，对每一块的J个数据分别取
绝对值后，再将获得的J个绝对值求平均值，最后从获得的多个平均值中选择最大的S个值 的标号赋值给所述最匹配的子空间t7，t7的值对应分块向量Group中的分组号，所述分组号 为 1，2，···,#；
步骤五、将步骤四获得的最匹配的子空间t7中多个分组号所对应的测量矩阵Φ中的 列向量组成的矩阵*赋值给恢复矩阵 ;
步骤六、计算本次迭代的残差r7
的伪逆矩阵；
步骤七、判断Ih I2 > IKr-IIl2成立与否，其中IPII2表示2-范数；
如果成立，则令对印
义然后执行步骤八；如果不成立，则直接执行 步骤八；
步骤八、令/4+1，返回执行步骤二，直到迭代次数1大于分组数见或者残差小于算法 迭代误差err时，迭代结束，执行步骤九；步骤九、输出源信号X的重构向量χ，重构向量χ为
完成对源信号X的压缩感知重构。压缩感知处理信号的基本模型如下
y = Φχ
Il-1I1 xe R^xl, |supp(x)| 4K《N
表示信号的ο-范数，即信号值不为ο的个数。# e 为测量矩阵，且
Bid在拥有了观测信号y和测量矩阵Φ的情况下，可以用0-范数意义下的优化问题来 重构或逼近源信号X。考虑另外一种类型的稀疏信号-块稀疏信号(Block-sparse Signal)，如下所示
T γι
A—J, .·,, A^cI 3J ·. .J Λ》『1
χ[1]s[2]χ[Μ]
其中Ar=IXi/, χ[刀，σ=1，···，#)为一子块(Sub-block)。当 /=1时，块稀疏退化成公 式一所描述的一般意义下的稀疏。如果向量X称为块I稀疏信号，则χ[7]，σ=ι，···，#)至 多有f个不为O的欧几里德(Euclidean)范数，定义
M^ 二 >°)
=1
0
> 也
__ Au
1乐
IN
J 5 「 ◎
i-v----
11
_
Jti
N
Ic
中 其
{
所以，向量χ称为块ι稀疏信号，则满足1Mb ^。初始化块稀疏*的选取应该满足1< K J为源信号X的真实块稀疏度，一般 地，可以选A为较小的值，当左=1时是最安全的选择，k值的大小和算法运行的效率有关，k 越大本发明的效率越高。
具体实施方式
二、本实施方式是对具体实施方式
一的进一步说明，步骤一中初始 化块稀疏度K，迭代误差err=10-5，测量矩阵《||E Rmxji服从高斯分布。
具体实施方式
三、本实施方式是对具体实施方式
二的进一步说明，测量矩阵 Φ = [^5 φ2,伊财],其中炉广jE{l,2,...，iV}为测量矩阵φ的归一化列向量，
为原始测量矩阵的列向量，表示2-范数。
IHL
具体实施方式
四、结合图2和图3说明本实施方式，本实施方式是将本发明方法 和混合Λ/Λ优化(Mixed IJI1 Optimization Program, L-0PT)算法、块稀疏匹配追踪 (Block-sparse Matching Pursuit, BMP)算法、±夬稀疏正交匹配追踪法(Block-sparse Orthogonal Matching Pursuit, BOMP)分别应用于块稀疏信号的重构，并计算每一种方法 的重构概率进行对比。仿真实验按以下步骤进行
首先、随机产生一个高斯分布测量矩阵^e Mmxw ,给定分块向量GrouP的分组数目
#以及块稀疏度&随机选定^个分块，分别在这^个块上赋值得到所需的仿真观测信号y (幅度采用高斯分布或0-1信号)；
然后、观测信号 拖，利用上述4种重构方法得到重构信号~，若 y = Φχχ
广则
重构成功；
最后、对每种重构方法运行500次，并计算重构概率。本实验过程中，观测信号y幅值分别采用高斯分布的信号和0-1的二值信号进行 实验。测量矩阵Φ行数 =80，列数#=160，分组数1=20，源信号χ的块稀疏度#=1，2，…，12，
计算每种算法在不同f值下的重构概率，并绘制重构概率随块稀疏度的变化曲线。实验 结果如图2和图3所示所示图2为观测信号y幅值为高斯分布信号的实验结果，图2中 一为曲线1，为本发明的方法的重构概率曲线，—为曲线2，为BMP方法的重构概 率曲线，-O ‘曲线3，为BOMP方法的重构概率曲线，_χ-曲线4，为L-OPT方法的重构概 率曲线，图3为观测信号y幅值为0-1 二值信号的实验结果，图3中^^为曲线5，为本 发明的方法的重构概率曲线，-O-为曲线6，为BMP方法的重构概率曲线，· -曲线7， 为BOMP方法的重构概率曲线，-M -曲线8，为L-OPT方法的重构概率曲线。从图2和图3中可见，无论对于哪类块稀疏信号，本发明方法的重构概率比L-OPT 方法都有大幅提高；对于幅值为高斯分布的信号，本发明方法略高于BMP和BOMP方法，而对 于0-1的二值信号，本发明方法的重构概率有显著提高，因此本发明方法尤其适用于实际 中的脉冲、黑白图像类等二值信号。
权利要求
一种块稀疏度未知的稀疏信号压缩感知重构方法，其特征在于具体过程如下采集获得源信号x的观测信号为y，观测信号y表示为yT=[y1, y2, y3……, ym]，其中m为观测信号y的长度，步骤一、稀疏信号压缩感知重构过程初始化设定初始化块稀疏度k，，K为源信号x的真实块稀疏度，，初始化测量矩阵，设定迭代误差err，分块向量Group，分块向量Group如下表示其中，d为分块向量Group的子块长度，设定残差初始值r0=y，恢复矩阵，步长step=1，信号支撑集大小S=k，源信号x的重构向量；步骤二、计算获得与残差rl 1最匹配的子空间il，其中rl 1表示第l 1次迭代产生的残差，il表示第l次迭代的最匹配的子空间，，具体过程为即测量矩阵每一块的转置与上一次迭代的残差rl 1进行乘积操作后，对每一块的d个数据分别取绝对值后，再求获得的d个绝对值的平均值，最后从获得的M个平均值中选择最大的S个值的标号赋值给所述最匹配的子空间il，il的值对应分块向量Group中的分组号，所述分组号为1,2,…,M；步骤三、将步骤二获得的最匹配的子空间il中对应的测量矩阵的列向量组成的矩阵与上一次迭代获得的恢复矩阵的并集赋值给初始恢复矩阵；步骤四、计算获得第l次迭代与观测信号y最匹配的子空间tl，，具体过程为其中矩阵中j为的分块标号，为的伪逆矩阵，即矩阵的每一块与观测信号y进行乘积操作后，对每一块的d个数据分别取绝对值后，再将获得的d个绝对值求平均值，最后从获得的多个平均值中选择最大的S个值的标号赋值给所述最匹配的子空间tl，tl的值对应分块向量Group中的分组号，所述分组号为1,2,…,M；步骤五、将步骤四获得的最匹配的子空间tl中多个分组号所对应的测量矩阵中的列向量组成的矩阵赋值给恢复矩阵；步骤六、计算本次迭代的残差rl；其中为的伪逆矩阵；步骤七、判断成立与否，其中表示2 范数；如果成立，则令step= step+1，S= step×S，然后执行步骤八；如果不成立，则直接执行步骤八；步骤八、令l=l+1，返回执行步骤二，直到迭代次数l大于分组数M，或者残差小于算法迭代误差err时，迭代结束，执行步骤九；步骤九、输出源信号x的重构向量，重构向量为，完成对源信号x的压缩感知重构。691584dest_path_image001.jpg,828167dest_path_image002.jpg,683997dest_path_image003.jpg,556324dest_path_image005.jpg,231019dest_path_image006.jpg,944284dest_path_image007.jpg,311812dest_path_image008.jpg,740388dest_path_image009.jpg,953195dest_path_image010.jpg,783616dest_path_image010.jpg,56466dest_path_image011.jpg,910021dest_path_image012.jpg,172857dest_path_image013.jpg,608517dest_path_image014.jpg,567115dest_path_image015.jpg,455436dest_path_image016.jpg,524892dest_path_image017.jpg,815059dest_path_image018.jpg,944558dest_path_image019.jpg,320176dest_path_image020.jpg,930673dest_path_image021.jpg,340926dest_path_image010.jpg,454376dest_path_image022.jpg,566557dest_path_image023.jpg,728548dest_path_image024.jpg,976996dest_path_image025.jpg,199030dest_path_image026.jpg,798507dest_path_image027.jpg,764189dest_path_image028.jpg,129793dest_path_image029.jpg,257149dest_path_image030.jpg,343922dest_path_image030.jpg,847716dest_path_image031.jpg
2.根据权利要求1所述的一种块稀疏度未知的稀疏信号压缩感知重构方法，其特征在 于步骤一中初始化块稀疏度K，迭代误差err=10-5，测量矩阵<iie Rmx^服从高斯分布。
3.根据权利要求2所述的一种块稀疏度未知的稀疏信号压缩感知重构方法，其特征在 于测量矩阵Φ =[街，釣，-·, φΝ] ’ Iljl1 Ψ, ’ ye {1,2,...,Λ7}为测丨ι!:矩阵φ的D丨■化列向量，炉;j 为原始测量矩阵的列向量，表示2-范数.
全文摘要
一种块稀疏度未知的稀疏信号压缩感知重构方法，它涉及压缩感知技术领域，具体涉及对块稀疏信号的重构方法。本发明通过初始化块稀疏度k，对每一个块稀疏信号的迭代，找到信号支撑集的一个子集，随着不断迭代，块稀疏度随之增加，最后找到整个源信号x的支撑集，从而达到重构源信号x的目的，本发明采用多次迭代并修正支撑集的思想来重构信号，重构的精度高，与现有的块稀疏度匹配追踪以及正交匹配追踪方法相比，对块稀疏信号的重构概率大，不存在过匹配现象。且本发明不需要以块稀疏度作为先验知识，尤其适用于块稀疏度未知信号重构领域。
文档编号H03M7/30GK101908889SQ20101024032
公开日2010年12月8日 申请日期2010年7月30日 优先权日2010年7月30日
发明者乔立岩, 付宁, 彭喜元, 曹离然, 马云彤 申请人:哈尔滨工业大学