基于精炼和分块搜索的重构稀疏信号方法
【专利摘要】本发明公开了一种基于精炼和分块搜索的重构稀疏信号方法，用以重构压缩感知(CS)技术框架下的稀疏信号。本发明方法利用MP算法简单的更新准则进行分块搜索，减少了运用最小二乘方法的次数，从而降低了OMP算法的复杂度；采用了精炼的方法提高重构性能，先进行扩展的索引集搜索，再通过最小二乘方法精炼出已知个数的索引；通过调节分块大小和扩展索引集维数实现重构性能与复杂度的权衡。本发明方法能够有效提高OMP算法的重构性能，并大大提高了重构速度。
【专利说明】基于精炼和分块搜索的重构稀疏信号方法
【技术领域】
[0001]本发明属于信号处理技术，涉及压缩感知技术及数据压缩【技术领域】。
【背景技术】
[0002]传统的信号采样以奈奎斯特采样定理为基础。在获取信号时，为了不丢失信号的信息，采样频率必须大于信号中最高频率的两倍，才能精确重构信号。但是随着科技的迅速发展，高分辨率的数码装置的采样产生了庞大的数据，如何更高效地处理这些数据并最大限度地节省存储传输的成本是一大难题。实际上采样得到的大部分数据是不重要的，在信号或图像的处理过程中，只保留了某些重要的数据，舍弃了大量的剩余数据，重构后的信号或图像并不会引起视觉上的差异。于是科学家们提出一个构想，既然采集到的数据大部分都是不重要的，可以被丢弃，能否直接地采集那部分重要的、最后没有被丢弃的数据，并且能够精确地重构原始信号或图像。
[0003]在2004年，由Donoho等人提出了压缩感知(CS)技术。压缩感知技术表明:如果信号通过某种变换(如傅立叶变换，小波变换等)后，是可稀疏表示或可压缩的，则可设计一个与变换基不相关的测量矩阵测量信号，得到的测量值通过求解优化问题，可实现信号的精确或近似重构。测量后，信号X由N维减少到信号y的M维(M〈〈N)，这M个测量值只包含了信号的重要信息。信号的观测过程是非自适应的，测量矩阵的设计不依赖于信号的结构。压缩感知的应用很大程 度地减少测量时间、采样速率及测量设备的数量。
[0004]信号重构是压缩感知领域的核心技术之一，从测量值中求得稀疏解是一个多项式复杂程度的非确定性(NP)问题。在文献[D.L.Donoho, “Compressed sensing, ^IEEE Trans.1nform.Theory, vol.52, n0.4, pp.1289-1306，2006.]中，方程组 y = Φχ 的最小范数解，在测量矩阵满足限制等距准则(RIP)条件时，可等同于求解X的最小A范数。最小A范数解可通过专门的凸优化迭代技术求得，但是其涉及的复杂度较大，并不适用于实际应用。近年来，凸优化技术的一种替代方法贪婪算法得到了广泛的关注。
[0005]贪婪算法的优点是较低的计算复杂度和简单的几何解释。匹配追踪(MP)算法[S.G.MalIat and Z.Zhang, “Matching pursuit with time-frequencydictionaries,，，IEEE Transactions on Signal Processing, vol.41, n0.12, pp.3397-3415,Dec.1993.]是一种较早的贪婪算法，该算法在迭代过程中逐一挑选出候选索引，并直接把最大的投影值作为非零元素的值。与此对比的是，正交匹配追踪(OMP)算法[J.A.Tropp andA.C.Gilbert, “Signal recovery from random measurements via orthogonal matchingpursuit, ” IEEE Trans.1nform.Theory, vol.53, n0.12, pp.4655-4666，Dec.2007.]通过最小二乘方法修正这些非零值。另外还有一些OMP算法的改进算法，如ROMP算法[D.Needelland R.Vershynin, “Signal recovery from incomplete and inaccurate measurementsvia regularized orthogonal matching pursuit, ” IEEE Journal of Selected TopicsinmSignal Processing, vol.4, n0.2, pp.310-316，Apr.2010.], StOMP算法[D.L.Donoho, Y.Tsaig, 1.Drori, and J.-L.Starck, “Sparse solution of underdetermined systems oflinear equations by stagewise orthogonal matching pursuit,，’IEEE Trans.1nform.Theory, vol.58, n0.2, pp.1094-1121，2012.], SP算法[ff.Dai and 0.Milenkovic, “Subspacepursuit for compressive sensing signal reconstruction, ” IEEE Trans.1nform.Theory, vol.55，n0.5，pp.2230-2249，May.2009.]等。ROMP 算法和 StOMP 算法在每次迭代中选择多个索引，因此能降低复杂度，但并不能带来明显的性能提升。SP算法引进了回溯追踪和精炼策略，性能得到提升的同时带来了复杂度的增加，特别是在稀疏度大的情况。
[0006]无噪情况的压缩感知问题就是从线性测量值y e RM(y = Φχ)中重构出原始稀疏信号X e Rn，其中X中只包含K (K << N)个非零值，这些非零值的位置定义为索引集Λ =U e [1，N]: Xi# 0}，稀疏度定义为 I A I = K,
[0007]Φ e Rmxn为测量矩阵。贪婪算法以一种次优化迭代的方法解决上诉问题，在迭代过程中依次寻找一个或多个索引(基于某种选取准则)，再将这些索引对测量值y的影响消除(基于某种更新准则)。最简单的选取准则是选取Φτι.中幅值最大的索引，其中r表示残差。
[0008]MP算法是一种比OMP算法更早的贪婪算法，该算法直接利用最大的投影值作为非零元素的值，用来更新残差。
[0009]作为对比的是，OMP算法利用最小二乘方法优化选中索引对应的非零值。这样做的好处就是误差很小，但计算量会大大增加。MP算法比OMP算法快得多，因为它不需要最小二乘计算，但是重构性能会变差。
[0010]对于OMP算法，一旦选中的原子加入索引集，将永远不会被移除。SP算法采用了一种简单的方法重新判断已加入索引集原子的可靠性。在每步迭代中，SP算法测试2K个索弓丨，并通过最小二乘方法选出K个最优的索引。某个在上一次迭代中认为是正确的索引，有可能在下次迭代中被认为是错误的，从而从索引集中移除。这种精炼索引集的方法，一方面可以提高重构性能；另一方面，最小二乘方法涉及的索引维数增长到了 2K，将导致计算量大大增加。
【发明内容】
[0011]为了解决上述问题，从用压缩感知(CS)技术压缩后的信号中重构出原始稀疏信号，本发明提供了一种基于精炼和分块搜索的重构稀疏信号方法(Block-refined 0ΜΡ，简称BR0MP)，包括如下步骤:
[0012](a)接收压缩后的测量信号，设定分块搜索的索引数量和扩展的索引集维数；
[0013](b)利用匹配追踪(MP)算法的更新准则找出若干索引，其数目等于分块搜索的索引数量；
[0014](C)将得到的索引加入索引集，利用最小二乘方法更新残差；
[0015](d)当索引集的维数超过扩展的索引集维数或残差满足精度要求时，进行下面步骤；否则，继续执行步骤(b)?步骤(c)；
[0016](e)利用最小二乘方法求得索引集对应的稀疏信号的数值，并选取幅值最大的索弓K该索引可以是若干个)作为最终的索引集，索引集的维数等于原信号的稀疏度；
[0017](f)再次利用最小二乘方法求得精炼后的索引集对应位置的数值；
[0018](g)输出索引集及其对应位置的数值。[0019]其中，步骤(b)利用了 MP算法的更新准则进行分块搜索，从而加快了正交匹配追踪(OMP)算法的重构速度；步骤(d)和步骤(e)进行了扩展的索引集搜索，并对此索引集进行精炼，从而提高了 OMP算法的重构性能；步骤(a)可通过调节分块大小和扩展索引集维数实现重构性能与复杂度的权衡。
[0020]本发明要用到的符号含义说明:
[0021]K:信号稀疏度；
[0022]Φ:MXN维测量矩阵，其每列是归一化的单位向量；
[0023]y:测量值；
[0024]B:分块大小；
[0025]σ:精度要求；
[0026]1:外循环迭代次数；
[0027]s:扩展的额外索引集维数；
[0028]OMP 算法:文献[J.A.Tropp and A.C.Gilbert, “Signal recovery from randommeasurements via orthogonal matching pursuit, ” IEEE Trans.1nform.Theory, vol.53，n0.12，pp.4655-4666，Dec.2007.]中给出的一种重构稀疏信号方法。
[0029]SP 算法:文献[W.Dai and 0.Milenkovic, “Subspace pursuit for compressivesensing signal reconstruction, ” IEEE Trans.1nform.Theory, vol.55, n0.5, pp.2230-2249，May.2009.]中给出的一种重构稀疏信号方法。
[0030]本发明方法利用了 MP算法的更新准则进行分块搜索，重构速度快；并通过对索引集进行精炼，提高了 OMP算法的重构性能；通过调节分块大小和扩展索引集维数，实现重构性能与复杂度的权衡；本发明方法应用于典型的压缩感知重构，其性能优于OMP算法，接近SP算法，复杂度低于OMP算法和SP算法。
【专利附图】
【附图说明】
[0031]图1是本发明BROMP方法的总流程图。
[0032]图2是本发明BROMP方法与OMP算法和SP算法重构高斯稀疏信号时的性能比较图。
[0033]图3是本发明BROMP方法与OMP算法和SP算法重构O -1稀疏信号时的性能比较图。
[0034]图4是本发明BROMP方法与OMP算法和SP算法重构高斯稀疏信号时的耗时比较图。
[0035]图5是本发明BROMP方法与OMP算法和SP算法重构O -1稀疏信号时的耗时比较图。
【具体实施方式】
[0036]下面结合附图和实施例对本发明作进一步详细说明。
[0037]如图1所示，本发明方法流程如下所示:
[0038](I)输入利用压缩感知技术压缩后的信号y，设定分块大小B、扩展的索引集维数s以及重构精度σ (稀疏度K和归一化的测量矩阵Φ已知)；[0039](2)外循环初始化:残差r° = y,索引集Λ° = 0，外循环迭代次数￡ = I ；
[0040](3)内循环初始化:残差< =rM，索引集Λ°ρ =0，内循环迭代次数q = I ;
[0041](4)内循环残差< 与测量矩阵Φ进行相关运算，得到相关值V = ;
[0042](5)将相关结果hq中的最大值对应的索引加入内循环的索引集，K = N'丨,lIJ ja = argmax,.| h"( j) | j ；
[0043](6)利用MP算法的更新准则更新内循环残差< =r；—1 -沪(α)Φ(α)，内循环迭代次数q+Ι ;如果q = B，执行下面步骤；否则，继续循环执行步骤(4)~(6)；
[0044](7)将内循环得到的索引集Λ=加入外循环索引集Λ._ n|i, Λ__ =Λ" 'IJa；:；
[0045](8)利用最小二乘方法更新外循环残差，y =y-Φ、(Φ； V),外循环迭代次数￡ + 1 ；其中+表示最小二乘求违逆，φΛ?表示索引集Ai在测量矩阵Φ中对应的列；
[0046](9)如果残差满足精度要求Il r Il >^7或外循环索引集的维数|Λ< |>尤+ (S，执行下面步骤；否则，继续循环执行步骤(3)~(9)；
[0047](?ο)执行精炼操作:如果外循环索引集的维数|λ \>κ最终索引Λ =.!115〗y中最大的K个元素对应的索引}，并计算？ iA =Φ> ；
[0048](11)输出最终索引集Λ以及其对应位置原信号的值iA。
[0049]本发明方法执行「(K+S)/B]次外循环(其中O < s < K,例如「&/16],「尤/8],「尤/4]，卜]表示不小于a的最小整数)，与OMP算法不同的是，本发明方
法选取B个索引，OMP算法只选取一个。在每次外循环中，本发明方法利用MP算法的更新准则选取出B个索引，因为MP算法在特别稀疏信号的重构上快速有效。一旦选出的B个索引加入到索引集中，本发明方法利用最小二乘方法求得最小误差的稀疏信号元素值。当索引集的维数超过K+s时，外循环迭代停止。扩展的索引集搜索允许索引集中的有少许错误的索引，如果某个错误的索引在某次迭代中加入了索引集，本发明方法利用精炼策略剔除这些错误索引，所以该方法能够提升重构性能。本发明方法计算最小二乘的次数为「(K+s)/B>l,与OMP算法相比，减少了尤-「(K+S)/Bl-1次，所以该方法能够加快重构速度。
[0050]本发明中，仿真参数M = 256，N = 128，K ( 64，仿真次数500，精确重构概率指输出SNR≥IOOdB占总次数的比例，重构时间是500次重构耗时总和。SP算法中的停止迭代
准则:||/Il2 <口或|<1//8。采用参数，B = 2，4，8, σ =10'
[0051]从图2和图4可以看出，本发明方法在重构高斯稀疏信号时的性能和时间上都较OMP算法有明显优势，重构时间随着分块大小B的增大而减少。与SP算法相比，本发明方法在重构时间较为理想，特别是稀疏度K较大时。当B=2，4时，本发明方法与SP算法性能很接近。
[0052]从图3和图5可以看出，本发明方法在重构O -1稀疏信号时的性能和时间上都较OMP算法较有优势，B的改变对重构性能影响很小。
【权利要求】
1.一种基于精炼和分块搜索的重构稀疏信号方法，其特征在于，包括如下步骤: Ca)接收压缩后的测量信号，设定分块搜索的索引数量和扩展的索引集维数； (b)利用MP算法的更新准则找出若干索引，其数目等于分块搜索的索引数量； (c)将得到的索引加入索引集，利用最小二乘方法更新残差； Cd)当索引集的维数超过扩展的索引集维数或残差满足精度要求时，进行步骤(e);否贝U，继续执行步骤(b)?步骤(C); (e)利用最小二乘方法求得索引集对应的稀疏信号的数值，并选取幅值最大的索引作为最终的索引集，索引集的维数等于原信号的稀疏度； Cf)再次利用最小二乘方法求得精炼后的索引集对应位置的数值； (g)输出索引集及其对应位置的数值。
【文档编号】H03M7/30GK103944581SQ201410151554
【公开日】2014年7月23日 申请日期:2014年4月15日 优先权日:2014年4月15日
【发明者】吴晓富, 嵇赢, 颜俊, 朱卫平 申请人:南京邮电大学