《一种压缩感知信号重构方法.pdf》由会员分享，可在线阅读，更多相关《一种压缩感知信号重构方法.pdf（10页珍藏版）》请在专利查询网上搜索。
1、(10)申请公布号 CN 102882530 A (43)申请公布日 2013.01.16 C N 1 0 2 8 8 2 5 3 0 A *CN102882530A* (21)申请号 201210343893.0 (22)申请日 2012.09.17 H03M 7/30(2006.01) (71)申请人南京邮电大学 地址 210003 江苏省南京市新模范马路66 号 (72)发明人杨真真 杨震 (74)专利代理机构南京经纬专利商标代理有限 公司 32200 代理人杨楠 (54) 发明名称 一种压缩感知信号重构方法 (57) 摘要 本发明公开了一种压缩感知信号重构方法， 属于信号处理技术领域。。
2、本发明针对-正则项 化问题求解比较困难，将压缩感知信号的稀疏域 的-正则化问题通过变量分裂技术转化为与之 等价的，比-正则项化问题更能体现信号稀疏 特性的约束-正则项化问题，使得重构信号的 精度更高。本发明进一步利用快速交替方向乘子 法对约束-正则项化问题进行求解，对交替方 向乘子法算法的变量进行了二次更新，并更新了 乘子，加快了优化求解的收敛速度。相比现有技 术，本发明方法具有更高的重构精度及更快的收 敛速度。 (51)Int.Cl. 权利要求书1页 说明书7页 附图1页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书 1 页 说明书 7 页 附图 1 页 1/1页 。
3、2 1.一种压缩感知信号重构方法，其特征在于，首先通过求解以下优化问题得到原始信 号的稀疏系数 ： 其中，为压缩感知采样的观测矩阵，为辅助变量，为对原始信号进行压缩感知采 样获得的观测向量，为正则化参数； 然后根据得到的稀疏系数重构出原始信号。 2.如权利要求1所述压缩感知信号重构方法，其特征在于，采用交替方向乘子法进行 所述优化问题的求解。 3.如权利要求1所述压缩感知信号重构方法，其特征在于，采用快速交替方向乘子法 进行所述优化问题的求解，具体包括以下步骤： 步骤1、初始化：设置初始惩罚参数，初始辅助变量，初始乘子，初始变量 ，初始迭代次数； 步骤2、更新稀疏系数变量： ，其中为单位矩阵；。
4、 步骤3、更新变量：； 步骤4、二次更新稀疏系数变量：； 步骤5、更新辅助变量：，其中为阈值为的硬阈 值函数； 步骤6、二次更新辅助变量：； 步骤7、更新乘子：； 步骤8、判断是否满足终止条件，满足则停止；否则，令，转步骤 4。 权 利 要 求 书CN 102882530 A 1/7页 3 一种压缩感知信号重构方法 技术领域 0001 本发明涉及一种压缩感知(Compressed Sensing，简称CS)信号重构方法，属于信 号处理技术领域。 背景技术 0002 针对具有稀疏特性信号的压缩感知(Compressed Sensing，CS)理论由Donoho等 人于2004年提出。它在保证信号。
5、不受损失的情况下，用远低于Nyquist采样定理要求的速 率采集信号，同时又不损失信息，能够完全恢复信号，是信号处理领域的一大革命性成果。 在CS理论框架下，采样速率不再取决于信号的带宽，而取决于信息在信号中的结构和内 容。 0003 CS理论中的核心问题是压缩信号的重构问题，如何设计复杂度低、收敛速度快、鲁 棒性强的重构算法一直是CS理论重构算法研究的目标。另一方面，实际的信号处理环境 中，采集（或者压缩感知得到）的信号，经常是含有噪声的。在使用优化方法重构含噪信号 时，不同之处在于重构过程所使用的优化目标函数的形式不同，参数的设置不同。应用不 同的优化目标函数，导致信号的重建效果也不尽相同。
6、。CS中重构含噪信号的原始模型为 l 0 -正则化问题，但该问题是一个非凸问题，求解十分困难，迭代硬阈值(Iterative Hard Thresholding,IHT)算法虽然可以通过迭代的方法求解该问题，但IHT算法求解的是该问 题局部极小值，重构精度低且算法收敛速度慢。所以，通常情况下，将其转化为l 1 -正则化问 题进行求解。虽然l 1 -正则化问题是一个严格凸问题，可以求得其全局唯一解，但是l 1 -正 则化问题只有在一定条件下才与l 0 -正则化问题等价，所求的解也只有在一定条件下才是 原问题的解，此时才能精确的重构原信号。 0004 目前关于含噪信号重构的领域已经有一些研究成果。。
7、交替方向乘子法 (Alternating Direction Method of Multipliers,ADMM)以其较好的性能和严格的理论 保证，成为众多学者研究的热点。例如用ADMM解决l 1 -正则化问题，提出了SALSA(Split augmented Lagrangian Shrinkage Algorithm)和C-SALSA(Constrained-SALSA)算法， 并把其应用于图像的去模糊、解卷积等问题，取得了较好的效果，但当观测矩阵为随机矩阵 时这两种算法的计算复杂度较高；又如用最速下降法更新迭代变量代替计算量大的求逆运 算，但其算法的收敛性没法保证；也有用ADMM来解决。
8、TV问题,并把其应用于图像的去模 糊、去噪问题；再如用ADMM来解决TVl 1 -l 2 问题，使用部分Fourier矩阵作为其观测矩阵，但 当观测矩阵为随机矩阵时收敛速度较慢；再如用IADMM(Inexact ADMM)来解决TV问题，加 快了算法的收敛速度，但其过程中求得的解只是原问题的近似解，导致重构的精度不高。大 部分都是将ADMM应用于图像的去噪、去模糊、解卷积等问题，很少有将ADMM应用于投影矩 阵为随机矩阵的信号重构问题；另外，目前尚无用ADMM来解决约束l 0 -正则化问题。 发明内容 0005 本发明所要解决的技术问题在于克服现有技术的不足，提供一种压缩感知信号重 说 明 书。
9、CN 102882530 A 2/7页 4 构方法，使得重构的信号的精度更高、算法的收敛速度更快。 0006 本发明具体采用以下技术方案解决上述技术问题： 0007 一种压缩感知信号重构方法，首先通过求解以下优化问题得到原始信号的稀疏系 数： 0008 0009 s.t.= 0010 其中，A为压缩感知采样的观测矩阵，为辅助变量，y为对原始信号进行压缩感 知采样获得的观测向量，0,)为正则化参数； 0011 然后根据得到的稀疏系数重构出原始信号。 0012 针对l 0 -正则项化问题求解比较困难，本发明将压缩感知信号的稀疏域的l 0 -正则 化问题通过变量分裂(Variable Splitti。
10、ng,VS)技术转化为与其等价的约束l 0 -正则项化 问题。由于l 0 -正则化问题比l 1 -正则化问题更能体现信号的稀疏特性，因此使得重构的信 号的精度更高。对于转化后的约束l 0 -正则项问题，可采用交替方向乘子法、对偶交替方向 乘子法、乘子法、非精确交替方向乘子法、交替线性化法等现有的方法进行求解。本发明为 了加快优化求解算法的迭代速度，提高信号重构的效率和实时性，进一步对交替方向乘子 法进行改进，具体包括以下步骤： 0013 步骤1、初始化：设置初始惩罚参数0，初始辅助变量 0 ，初始乘子d 0 ，初始变量 t 0 =1，初始迭代次数k=0； 0014 步骤2、更新稀疏系数变量： 。
11、0015 其中I为单位矩阵； 0016 步骤3、更新变量t： 0017 步骤4、二次更新稀疏系数变量： 0018 步骤5、更新辅助变量：其中hard(,)为阈值为 的硬阈值函数； 0019 步骤6、二次更新辅助变量： 0020 步骤7、更新乘子d：d k+1 =d k -( k+1 - k+1 )； 0021 步骤8、判断是否满足终止条件，满足则停止；否则，令k=k+1，转步骤4。 0022 相比现有技术，本发明具有以下有益效果： 0023 (1)通过变量分裂技术将压缩感知信号重构的l 0 -正则化问题转化为约束l 0 -正 则化问题，该正则化问题比l 1 -正则化问题更体现信号的稀疏特性，使。
12、得重构的信号的精度 更高； 0024 (2)在对约束l 0 -正则化问题求解时，用加速策略对ADMM的变量进行了二次更新， 加快了算法的收敛速度。 说 明 书CN 102882530 A 3/7页 5 附图说明 0025 图1为本发明的压缩感知信号重构系统框图； 0026 图2为采用各种重构算法得到的含噪图像压缩感知重构结果。 具体实施方式 0027 下面结合附图对本发明的技术方案进行详细说明： 0028 为了便于公众理解本发明的技术方案，下面先对压缩感知的基本知识做一简要介 绍。 0029 现代信号处理的一个关键基础是Shannon采样理论：一个信号可以无失真重建所 要求的离散样本数由其带宽。
13、决定。但是Shannon采样定理是一个信号重建的充分非必要条 件。而压缩感知作为一个新的采样理论，它可以在远小于Nyquist采样率的条件下获取信 号的离散样本，保证信号的无失真重建。压缩感知理论的核心思想主要包括两点。第一个 是信号的稀疏结构。传统的Shannon信号表示方法只开发利用了很少的被采样信号的先验 信息，即信号的带宽。但是，现实生活中很多广受关注的信号本身具有一些结构特点。相对 于带宽信息的自由度，这些结构特点是由信号的更小的一部分自由度所决定。换句话说，在 很少的信息损失情况下，这种信号可以用很少的数字编码表示。所以，在这种意义上，这种 信号是稀疏信号（或者可压缩信号）。另外一。
14、点是不相关特性。稀疏信号的有用信息的获取 可以通过一个非自适应的采样方法将信号压缩成较小的样本数据来完成。理论证明压缩感 知的采样方法只是一个简单的将信号与一组确定的波形进行相关的操作。这些波形要求是 与信号所在的稀疏空间不相关的。压缩感知可广泛应用于信息论、图像处理、地球科学、光 学/微波成像、模式识别、无线通信、大气、地质等领域。 0030 具体地讲，对于任意信号xR N ，可以用N1维正交基向量的线性组合来表示 0031 0032 其中，为投影系数，= T x为投影系数矢量。若是K-稀疏的，可以用一 个与不相关的矩阵AR MN (其中MN)对进行线性测量，得到观测向量yR M 为 003。
15、3 y=A+e (2) 0034 其中，e为噪声(当不含有噪声时，令e=0即可)，AR MN 为与不相关的观测 矩阵。重构的方法是通过求解如下最优化问题 0035 0036 其中0,)为正则化参数(当=0即为不含噪声时的重构模型)。通过 求解该优化问题即可得到稀疏系数，进而可得到原始信号x。 0037 正如背景技术部分所述，式（3）的l 0 -正则化问题是一个非凸问题，求解十分困难， 因此通常将其转化为l 1 -正则化问题进行求解。虽然l 1 -正则化问题是一个严格凸问题，可 以求得其全局唯一解，但其具有较大局限性。 0038 本发明的思路是利用变量分裂技术，引入辅助变量，将稀疏域含噪信号的重。
16、构 说 明 书CN 102882530 A 4/7页 6 问题转化为如下的约束l 0 -正则化问题 0039 0040 s.t.= 0041 此时即可用具有较好的性能和严格的理论保证的交替方向乘子法（ADMM）对式 (4)的约束l 0 -正则化问题进行求解，得到信号的更加精确的稀疏系数，进而得到原始信号。 0042 约束l 0 -正则化问题(4)的乘子函数为 0043 0044 其中，R M 为乘子，0,)为惩罚参数。 0045 对于变量得到： 0046 0047 0048 (6) 0049 其中，也称为乘子，I为单位矩阵。当矩阵AR MN 为随机矩阵时， (A T A+I) -1 的计算量为。
17、O(N 3 )，为了减少其计算量，根据Sherman-Morrison-Woodburg定理 得到 0050 0051 问题(6)是一个严格凸问题，其解唯一。从而得到 0052 0053 对于变量得到： 0054 0055 0056 其中hard(,)为阈值为的硬阈值函数。为了加速ADMM的收敛速度，本发明 对和进行了二次更新，其二次更新的形式分别为 0057 0058 0059 其中t 0 =1，为算法迭代更新的步长。且对乘子进行了 更新 0060 d k+1 =d k -( k+1 - k+1 ) (12) 说 明 书CN 102882530 A 5/7页 7 0061 本发明称上述的改进。
18、交替方向乘子法算法为快速ADMM，即FADMM。之所以FADMM会 加快算法的收敛速度，是因为在迭代过程中不仅仅用到当前点的信息(即算法ADMM)，而 且用到的信息，使得每步的迭代结果更准确，在终止准则不变时，使得算法的迭代次数 减少；而且FADMM的主要计算量仍在式(6)和(9)，对和进行的二次更新并无实质性 的增加计算量，即FADMM每次迭代的运算量并无实质性的增加。从而导致FADMM的运行时 间因其收敛速度快而减少。之所以取t 0 =1是为了保证对于任意的k0， 均有从而保证FADMM的超线性收敛速度。 0062 综上所述，利用FADMM求解约束l 0 -正则化含噪信号重构问题的算法步骤。
19、为 0063 1)初始化：设置初始惩罚参数0，初始辅助变量 0 ，初始乘子d 0 ，初始变量 t 0 =1，迭代次数k=0； 0064 2)更新稀疏系数变量： 0065 3)更新变量t： 0066 4)二次更新稀疏系数变量： 0067 5)更新辅助变量： 0068 6)二次更新辅助变量： 0069 7)更新乘子d：d k+1 =d k -( k+1 - k+1 )； 0070 8)判断是否满足终止条件，满足则停止；否则，令k=k+1，转步骤4)。 0071 定理1：令设为FADMM生成的序列，则对任意的 k0，有 0072 0073 其中， 0074 定理1说明FADMM具有超线性收敛性。 0。
20、075 图1显示了一种采用本发明的方法的压缩感知信号重构系统，如图所示，该系统 包括编码端和解码端，解码端对原始的含噪信号先进行小波变换，提取出原始含噪信号的 小波系数，然后对小波系数进行压缩感知采样，得到观测量；解码端采用本发明的FADMM算 说 明 书CN 102882530 A 6/7页 8 法对编码端发送的观测量进行优化求解，得到重构信号的小波系数，然后进行反小波变换， 即得到重构后的信号。 0076 为了验证本发明方法的效果，以含噪图像的压缩感知重构为例，对本发明的重构 方法与其他现有方法进行了比较。实验对象为Lena图像(大小为256256)，正交基矩阵 为小波基db3，观测矩阵A。
21、为服从高斯分布的随机矩阵，压缩比r=0.5，使用imnoise函数 对实验对象加入均值为0，方差为 2 的高斯噪声。分别采用峰值信噪比PSNR(单位为dB) 和运行时间t(单位为S)作为算法的重构性能和收敛速度的评价指标。重复进行10次实 验，各种算法的峰值信噪比PSNR和运行时间t均为10次实验的平均结果。并将本发明的 FADMM算法与SALSA、AIHT(Accelerated Iterative Hard Thresholding)算法、IHT算法、 TwIST(Two-step Iterative Shrinkage/Thresholding)算法、GPSR(Gradient Proj。
22、ection for Sparse Reconstruction)算法的性能和收敛速度进行了比较与分析。 0077 对实验对象的原图像加入均值为0，方差为 2 =0.001的高斯噪声，图2显示了实 验对象的原图像、含噪图像及各种算法重构的图像，可以看出，当压缩比r=0.5时，对原图 像加入均值为0，方差为 2 =0.001的高斯噪声，各种算法重构的图像效果均较好。 0078 表1显示了各种重构算法在不同噪声程度下，采用PSNR(单位为dB)作为衡量它 们重构性能指标的结果。 0079 表1不同噪声程度各种重构算法下的PSNR 0080 2 FADMM SALSA AIHT IHT TwIST 。
23、GPSR 0.001 32.055 32.042 29.879 29.736 31.969 32.038 0.005 28.039 27.771 25.286 25.038 27.953 27.685 0.01 25.668 24.531 22.804 22.641 25.601 25.646 0.05 19.838 18.937 16.932 16.602 19.082 18.883 0.1 17.608 15.763 14.647 14.354 17.491 17.537 0081 随着噪声方差 2 的增大，各种算法重构图像的PSNR减小。此外在相同噪声程度 下，FADMM算法重构图像的P。
24、SNR最大，IHT算法重构的图像的PSNR最小。而PSNR在一定程 度上反映算法的重建质量，同等条件下PSNR越大，算法重建质量越高。所以在相同噪声程 度下，FADMM算法重构图像的质量高于其他算法(即SALSA、AIHT算法、IHT算法、TwIST算 法、GPSR算法)重构图像的质量。 0082 表2显示了各种重构算法在不同噪声程度下，采用运行时间t(单位为S)作为衡 量它们重构性能指标的结果。 0083 表2不同噪声程度各种重构算法下的t 0084 说 明 书CN 102882530 A 7/7页 9 2 FADMM SALSA AIHT IHT TwIST GPSR 0.001 3.86。
25、0 4.695 4.561 11.767 24.031 56.375 0.005 4.249 5.454 5.079 11.924 30.352 58.239 0.01 4.203 7.218 4.960 11.999 38.281 59.996 0.05 4.405 11.852 5.454 12.010 42.399 60.925 0.1 4.842 12.406 5.639 12.313 48.420 62.817 0085 随着噪声方差 2 的增大，各种重构算法重构图像的运行时间t增大。此外相同 噪声程度下，FADMM算法的运行时间最短，GPSR算法的运行时间最长，各种算法的运行时间 快慢可以表示为t FADMM t AIHT t SALSA t IHT t TwIST t GPSR 。 0086 综合考虑峰值信噪比PSNR和运行时间t，FADMM算法的性能最好。 说 明 书CN 102882530 A 1/1页 10 图1 图2 说 明 书 附 图CN 102882530 A 10 。