压缩感知重构算法
语音压缩感知及其重构算法

误差 , 不 同的噪 声 强度 下 具有 较高 的鲁棒 性 . 在
关 键词 :次梯 度投 影 ; 正交 小波 ; 双 压缩感 知 ; 疏重 构 稀 中图分 类号 : N9 2 T 1 文 献标 志码 : A 文 章编 号 : 0 1 0 0 ( 0 1 0 -0 10 10 — 5 5 2 1 ) 1 0 -5 0
第4 卷第 1 1 期
21 0 1年 1月
东 南 大 学 学 报 (自然科 学版 )
J R L OF S T E T U VE S T ( aua S in eE io ) OU NA OU H AS NI R I Y N tr l ce c dt n i
V o . N o. 141 1 J n. 2 1 a 01
( S h o fIfr t n S in ea d E gn eig o te s U ies y Naj g2 0 9 , hn ) c o l o mai ce c n n ier ,S uh a t nv ri o n o n t ni 10 6 C ia n
立 日
感 知 ( s 语音 信 号重构 的 方案. P 算 法 首先根 据 压缩感 知 重构模 型 建 立包含 稀疏 重 构信 号 c ) AS M
并具有 随机 属 性 的 凸集 , 后运用 次梯 度投 影 的 思想将 该 凸集 的投 影 转化 为对 多个 闭合半 平 面 然 的投影 , 最后 将 更新 后 的稀 疏 重构 信 号投影 到 限定 集合 上 . 时 , 算 法设 计 了 自适应 调 节 膨胀 同 该 系数 的机制 以获 得 快速 收敛性 . 理论 分 析和 仿 真结果 表 明 , 算法 具有 快速 收敛 性和 较低 的重构 该
v r poet no t o v x st t rjc o no m lpec sd hl p c s ial , h p a d e rjci noc n e es o poet n it ut l l e a sae.Fn l teu dt t o i i o f y e
基于图像结构模型的压缩感知图像重构方法

基于图像结构模型的压缩感知图像重构方法随着数字图像在各个领域的广泛应用,如视频监控、医学影像以及移动通信等,对图像传输和存储的要求也越来越高。
然而,由于图像数据量庞大,传输和存储的成本也随之增加。
为了解决这一问题,压缩感知(Compressed Sensing,CS)技术被提出并逐渐得到应用。
压缩感知技术的基本思想是在图像采集中对图像进行压缩,即不直接采集完整的图像数据,而是对其进行稀疏采样,然后通过稀疏的采样数据来重构完整的图像。
而图像结构模型就是其中一种常用的重构方法之一。
图像结构模型是一种基于图像自身的特性进行建模和重构的方法。
它利用图像的边缘、纹理和结构等特征来提取图像信息,从而实现更加准确和高质量的图像重构。
下面将介绍基于图像结构模型的压缩感知图像重构方法的具体步骤和原理。
一、图像结构模型的建立在压缩感知图像重构过程中,首先需要建立图像结构模型。
这个步骤涉及到对图像的稀疏表示,常用的方法有小波变换、稀疏表示字典以及图像分割等。
小波变换是一种常用的图像分析和压缩方法,通过将图像进行小波变换来提取图像的频域信息,进而实现图像的稀疏表示。
稀疏表示字典则是通过提前建立一个字典,将图像的局部结构进行编码,从而实现图像的稀疏表示。
图像分割是将图像划分为若干个小块,每个小块可以看做是具有相似结构的局部区域,从而实现图像的稀疏表示。
二、图像重构算法建立好图像结构模型后,下一步就是利用稀疏采样数据对图像进行重构。
常用的图像重构算法有基于最小二乘法的估计(Least Squares,LS)、基于迭代阈值法的估计(Iterative Shrinkage-Thresholding Algorithm,ISTA)以及基于广义估计最小二乘法的估计(Generalized Estimation of Signal and Noise,GESPAR)等。
LS方法是一种常见的图像重构算法,它通过将图像重构问题转换成一个最小二乘问题,通过最小化重构图像与原始图像之间的欧式距离来进行重构。
压缩感知重构算法综述

收稿日期: 2013 - 01 - 14 ;
修订日期: 2013 - 02 - 16
基金项目: 陕西省自然科学基金( 2012JM8021 ) 作者简介: 李珅( 1980 - ) , 女, 助理研究员, 博士, 主要从事压缩感知和图像超分辨率分析方面的研究。Email: waterblue_333@ opt. ac. cn
( Compressed Sensing,CS ) 为解 决该问题 提供了 契机 。该理论指出: 对于稀 疏或可压 缩 的 信 号, 能够以远低 于奈奎 斯特频率 对其 进行采 样 , 并通过设计重 构算 法 来 精确 的 恢 复 该 信 号。 介绍 了 压 缩 感 知 理 论的 基 本 框架并讨论了 该理论 关于信 号压 缩 的采样过程 ; 综 述了压缩感 知理论的 重 构算 法 , 其 中 着 重介绍 了 最 优 化 算法和贪婪 算法并比较了各 种算 法 之间的 优 劣 , 最 后 探 讨 了 压 缩 感 知 理 论 重 构算 法 未 来 的 研究 重 点 。 通 为压 缩感知 重构算 法的 改进和 应 用 提 供 了理 论 依 据。 过对压缩感知理论重 构算法 较为 系 统的介绍 和比较 , 关键词: 压缩感 知 ; 信号 采样 ; 稀疏 ; 重构算法 中图分类号: TP301. 6 文献标志码: A 文章编号: 1007 - 2276 ( 2013 ) S1 - 0225 - 08
,
稀疏( sparsity) 和不相关( incoher - ence ) , 稀疏关于感 兴趣的信号, 它所表达的意思为: 连续时间信号的信息 率可能比根据其带宽所建议的小得多, 离散时间信号 所依赖自由度的数量比它的长度少得多 。 可以说, 许 多自然界的信号在某种程度上都是稀疏的或可压 缩 的, 当以合适的基 Ψ 来表示时, 信号可以有很多简练 的表达式 。不相关表达了一种含义即以 Ψ 稀疏表示 的信号一定可以在其所需要的域中展开 。 基于这两个 原则, 压缩感知理论指出, 长度为 n 的信号 X 在某组正
基于压缩感知的信号重构算法研究共3篇

基于压缩感知的信号重构算法研究共3篇基于压缩感知的信号重构算法研究1基于压缩感知的信号重构算法研究随着信息技术的发展以及现代通信系统的广泛应用,人们对于信号重构算法的研究也越来越深入。
其中,基于压缩感知的信号重构算法受到了广泛关注。
本文将从以下四个方面来探讨该算法的研究。
一、压缩感知的基本原理压缩感知的核心思想是将一个高维信号(如图像、音频等)映射到一个较低维的空间中,然后再通过一个线性投影方式将数据压缩。
利用测量矩阵可以将压缩后的数据重构到原来的高维空间中,并且能够利用未知信号的稀疏性完成恢复过程。
这种低维的表示方式可以使数据占用的空间大大减小,因此压缩感知成为了高效的信号采样方式。
二、常见的压缩感知算法常见的压缩感知算法包括OMP算法、CoSaMP算法、MPCP算法等。
其中OMP算法是一种迭代算法,用于寻找稀疏表示向量。
CoSaMP算法考虑到了噪声的影响,能够更准确地进行稀疏重构。
MPCP算法则是多向量压缩感知的拓展,用于处理多个信号的联合稀疏性问题。
三、压缩感知在图像压缩方面的应用基于压缩感知的信号重构算法在图像压缩方面的应用也是较为广泛的。
传统的JPEG和PNG等图像压缩算法虽然能够将图像进行压缩,但是重构后的图像质量较差,并且对于稀疏性较强的图像处理能力有限。
基于压缩感知的算法能够更好地处理稀疏性强的图像,同时也能够提高图像的显示效果。
四、压缩感知在音频处理方面的应用除了在图像处理方面的应用,基于压缩感知的信号重构算法在音频处理方面也具有广泛的应用前景。
例如在音频采样、去噪、提取声音等方面都有着极为广泛的应用。
此外,利用压缩感知的技术,人们还可以用较小的存储空间存储大量音乐等高质量音频数据。
综上所述,基于压缩感知的信号重构算法是一种高效且优越的信号处理方法,具有较广泛的应用前景。
在未来的研究中,我们可以结合更多的数据处理技术来提高算法的效率和精度基于压缩感知的信号重构算法在信号处理中具有广泛应用前景,能够更好地处理稀疏性较强的信号,并提高信号质量。
几种压缩感知算法

.1 压缩感知部分压缩感知算法主要可分为三类:贪婪迭代算法、凸凸优化(或最优化逼近方法)和基于贝叶斯框架提出的重构算法。
由于第三类方法注重信号的时间相关性,不适合图像处理问题,故目前的研究成果主要集中在前两类中。
目前已实现6中算法,分别为正交匹配追踪法(OMP)、迭代硬阈值法(IHT)、分段正交匹配追踪法(StOMP)、分段弱正交匹配追踪法(SwOMP)、广义正交匹配追踪(GOMP)、基追踪法(BP)。
1.1 正交匹配追踪法(OMP)在正交匹配追踪OMP中,残差是总与已经选择过的原子正交的。
这意味着一个原子不会被选择两次,结果会在有限的几步收敛。
OMP的算法如下(1)用x表示你的信号,初始化残差e0=x;(2)选择与e0内积绝对值最大的原子,表示为φ1;(3)将选择的原子作为列组成矩阵Φt,定义Φt列空间的正交投影算子为通过从e0减去其在Φt所张成空间上的正交投影得到残差e1;(4)对残差迭代执行(2)、(3)步;其中I为单位阵。
需要注意的是在迭代过程中Φt为所有被选择过的原子组成的矩阵,因此每次都是不同的,所以由它生成的正交投影算子矩阵P每次都是不同的。
(5)直到达到某个指定的停止准则后停止算法。
OMP减去的Pem是em在所有被选择过的原子组成的矩阵Φt所张成空间上的正交投影,而MP减去的Pem是em在本次被选择的原子φm所张成空间上的正交投影。
经OMP算法重构后的结果如下所示:算法的使用时间如下:1.2 迭代硬阈值法(IHT)目标函数为这里中的M应该指的是M-sparse,S应该指的是Surrogate。
这里要求:之后我们利用式对目标函数进行变形。
接着便是获得极值点:利用该式进行迭代可以得到极值点,我们需要的是最小值。
此时目标函数的最小值就得到了。
此时便得到我们需要的公式:我们要保证向量y的稀疏度不大于M,即,为了达到这一目标,要保留最大的M项(因为是平方,所以要取绝对值absolute value),剩余的置零(注意这里有个负号,所以要保留最大的M项)。
压缩感知重构算法——SP算法

压缩感知重构算法——SP算法SP(subspace pursuit)算法是压缩感知中⼀种⾮常重要的贪婪算法,它有较快的计算速度和较好的重构概率,在实际中应⽤较多。
本⽂给出了SP算法的matlab代码,以及相应的测试函数。
参考⽂献:Dai W, Milenkovic O. Subspace pursuit for compressive sensing signal reconstruction[J]. Information Theory, IEEE Transactions on, 2009, 55(5): 2230-2249.⽂献下载地址:Matlab代码:SP_paper.mfunction x=SP_paper(Phi,y,K)%SP算法%获取Phi矩阵的⾏数和列数[M,N]=size(Phi);%初始化步骤%将Phi的每列与y做相关,得到⼀个N*1的矩阵(列向量)correlation=Phi'*y;%对correlation取绝对值后排序,按从⼤到⼩的顺序[var,pos] = sort(abs(correlation),'descend');%声明⼀个空集T,⽤于记录Phi的列数标值T=[];T=union(T,pos(1:K));y_r=resid_paper(y,Phi(:,T));%迭代%使⽤如下形式的do---while结构% while(1)% if(condition)% break;% end% endcount=1;while(1)%根据残差计算待增加的列数,得到T_addcorrelation=Phi'*y_r;[var,pos] = sort(abs(correlation),'descend');T_add=union([],pos(1:K));%合并已有的T和T_addT=union(T,T_add);%x_p=((Phi(:,T)'*Phi(:,T))\eye(length(T)))*Phi(:,T)'*y;%proj_paper(y,Phi(:,T));%更新下标记录T[var,pos] = sort(abs(x_p),'descend');%取前K个最⼤值T=union([],T(pos(1:K)));%计算新的残差y_r_n=resid_paper(y,Phi(:,T));%判断是否退出循环,且置为最⼤迭代100次if(norm(y_r_n)>=norm(y_r) || count>100)break;end%若不退出循环,进⾏新⼀轮的迭代y_r=y_r_n;count=count+1;end%退出循环后,做最后的数据输出x=zeros(N,1);x(T)=((Phi(:,T)'*Phi(:,T))\eye(length(T)))*Phi(:,T)'*y;endfunction y_r=resid_paper(y,Phi)%计算y在Phi上的投影残差%获取矩阵Phi的⾏数和列数,M没有⽤[M,N]=size(Phi);%判断矩阵(Phi'*Phi)是否可逆if(rank(Phi'*Phi)~=N)error('矩阵不可逆');endy_p=Phi*((Phi'*Phi)\eye(N))*Phi'*y;y_r=y-y_p;end%%%%%%%%%%%%%%%%%%%%%%%%dataGen.mfunction [y,Phi,x]=dataGen(M,N,K)% 产⽣贪婪算法所需要的数据%⽣成-1/+1的原始信号xx = zeros(N,1);q = randperm(N); %y=randperm(n),是把1到n这些数随机打乱得到的⼀个数字序列。
压缩感知的重构算法

压缩感知的重构算法算法的重构是压缩感知中重要的一步,是压缩感知的关键之处。
因为重构算法关系着信号能否精确重建,国内外的研究学者致力于压缩感知的信号重建,并且取得了很大的进展,提出了很多的重构算法,每种算法都各有自己的优缺点,使用者可以根据自己的情况,选择适合自己的重构算法,大大增加了使用的灵活性,也为我们以后的研究提供了很大的方便。
压缩感知的重构算法主要分为三大类:1.组合算法2.贪婪算法3.凸松弛算法每种算法之中又包含几种算法,下面就把三类重构算法列举出来。
组合算法:先是对信号进行结构采样,然后再通过对采样的数据进行分组测试,最后完成信号的重构。
(1) 傅里叶采样(Fourier Representaion)(2) 链式追踪算法(Chaining Pursuit)(3) HHS追踪算法(Heavy Hitters On Steroids)贪婪算法:通过贪婪迭代的方式逐步逼近信号。
(1) 匹配追踪算法(Matching Pursuit MP)(2) 正交匹配追踪算法(Orthogonal Matching Pursuit OMP)(3) 分段正交匹配追踪算法(Stagewise Orthogonal Matching Pursuit StOMP)(4) 正则化正交匹配追踪算法(Regularized Orthogonal Matching Pursuit ROMP)(5) 稀疏自适应匹配追踪算法(Sparisty Adaptive Matching Pursuit SAMP)凸松弛算法:(1) 基追踪算法(Basis Pursuit BP)(2) 最小全变差算法(Total Variation TV)(3) 内点法(Interior-point Method)(4) 梯度投影算法(Gradient Projection)(5) 凸集交替投影算法(Projections Onto Convex Sets POCS)算法较多,但是并不是每一种算法都能够得到很好的应用,三类算法各有优缺点,组合算法需要观测的样本数目比较多但运算的效率最高,凸松弛算法计算量大但是需要观测的数量少重构的时候精度高,贪婪迭代算法对计算量和精度的要求居中,也是三种重构算法中应用最大的一种。
正交半迭代硬阈值(OSIHT)压缩感知重构算法

正交半迭代硬阈值(OSIHT)压缩感知重构算法杨爱萍;刘华平;何宇清;栗改【摘要】信号重构是压缩感知过程中的重要环节,迭代硬阈值(IHT)算法因具有较好的重构性能被广泛应用,但其收敛速度比较慢.近期提出的半迭代硬阈值算法(SIHT)虽然可实现快速收敛,但对测量矩阵的尺度缩放非常敏感,依赖性强,大大限制了其应用范围.受OMP对MP算法改进启发,对SIHT算法进行改进,提出了正交半迭代硬阈值(OSIHT)重构算法.该算法不仅取消了对测量矩阵的依赖性,还有效改善了图像重构质量,减少运行时间.【期刊名称】《计算机工程与应用》【年(卷),期】2016(052)005【总页数】5页(P79-83)【关键词】压缩感知;重构算法;迭代阈值;半迭代法【作者】杨爱萍;刘华平;何宇清;栗改【作者单位】天津大学电子信息工程学院,天津300072;天津大学电子信息工程学院,天津300072;天津大学电子信息工程学院,天津300072;天津大学电子信息工程学院,天津300072【正文语种】中文【中图分类】TP391近年来,Candès和Donoho提出的压缩感知(CS)理论[1-4],突破了传统采样定理的束缚,在获取信号的同时对数据进行适当的压缩,避免了传统采样压缩的资源浪费,缓解了信号采样、传输和存储过程中所面临的越来越大的压力,为信号采集技术带来了革命性的进展。
压缩感知理论主要包括:信号的稀疏表示、测量矩阵设计和信号的重构算法。
信号重构是由测量值及投影矩阵重构原始信号,是压缩感知理论的必备手段,是压缩感知处理过程最重要的环节。
信号压缩感知重构算法研究已经取得了不少的成果,其中典型的算法有凸松弛方法[5-6]、贪婪方法[7-9]、迭代阈值算法[10-12]等。
其中,迭代阈值算法因迭代简单、可单分量处理且具有较好的重构性能而被广泛应用。
但是传统的迭代阈值类算法具有较慢的收敛速度,文献[12]基于半迭代加速技术提出的半迭代硬阈值算法(Semi-Iterative Hard Thresholding algorithm,SIHT),可实现快速收敛,但该算法对测量矩阵的尺度缩放非常敏感,有很强的依赖性,在很大程度上限制了其应用范围。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
4.2.1 高斯白噪声影响
高斯白噪声均值为0,方差为0.01
图A 受高斯白噪声影响OMP算法重构
图B 受高斯白噪声影响信号CoSAMP算法重构
4.2.2 椒盐噪声影响
椒盐噪声噪声密度为0.05
图A 受椒盐噪声影响信号OMP算法重构
图B 受椒盐噪声影响信号IRLS算法重构
4.2.3 受噪声影响结果比较
4
2
0
SAMP
-2 -4
-6
0
50
100
150
200
250
300
2.7.2 各算法仿真实验分析
重构算法时间都较短,误差也都较小
2.7.2 各算法性能分析
随着稀疏度K的增加,到达一定的临界值以后,信号重构的概率逐步降低
3.1 改进算法
基于变步长的正则化回溯自适应追踪算法 2018年,王欣等人针对压缩感知重构过程中稀疏度 未知以及步长大小固定的情况,基于SAMP 算法 存在的而不足,提出一种基于变步长的正则化回
每次迭代时选择多列
2.4 压缩采样匹配追踪CoSaMP
每次迭代选择多个原子, 选择的原子在下次迭代中可能会被抛弃
2.5 压缩采样匹配追踪SP
后向追踪的方式,在每次迭代中将找出的错 误的原子剔出,再把新找的原子放入支撑集 中,最后进行逼近
2.6 压缩采样匹配追踪SAMP
稀疏度自适应匹配追踪算法(SAMP)解决了需己知稀疏度才能准确重 构信号的要求。该算法沿用SP算法的回溯思想,加入步长和分阶段选 择的思想,通过迭代次数的增加,支撑集增加,从而达到每次选择更 多的原子,更快的重构信号。但是,如何选择初始步长及固定步长对 该算法有一定影响。
3.3 迭代硬阈值算法IHT
文献:T. Blumensath and M. Davies, "Iterative
Hard Thresholding for Compressed Sensing,"
2008.
基本思想:给定一个初始的X0,然后通过以下的 阈值公式不断地迭代。
3.4 迭代加权最小二乘算法(IRLS)
6
4
2
5
0
0
-2
-4
-5
-6
-8
0
50
100
150
200
250
300
-10
0
50
100
150
200
250
300
GOMP
StOMP
2.7.1 仿真实验
StOMP算 法 重 构 信 号 8 Recovery Original 6
4
2
0
SP
-2
-4
-6
0
50
100
150
200
250
300
SP算 法 重 构 信 号 8 Recovery Original 6
压缩感知重构算法
日期:2018.1Biblioteka .27目录1 2 3
4 5
压缩感知理论
经典贪婪迭代算法理论和实验 改进的重构算法
二维信号仿真及抗噪声分析
结论和展望
1 压缩感知理论
1.1 压缩感知重构算法
0范式重构算法:贪婪迭代算法(MP,OMP,ROMP等) 迭代算法(IHT,NIHT,AIHT) 1范式重构算法:凸优化算法(ADMM,Interior point,Gradient Projection,GPSR,Homotopy
Method) 结构化稀疏:Model based CS,SBL,BSBL
2.1.1 经典重构算法---OMP
1
2.1.2 程序角度了解OMP
2.2 分段正交匹配追踪StOMP
每次迭代可以选择多个原子。 此算法的输入参数中没有信号稀疏度K, 通过门限来确定候选原子集
2.3 正则化正交匹配追踪ROMP
2.7.1 仿真实验
CoSaMp算 法 重 构 信 号 8 Recovery Original 6 4
2
0
-2
-4
-6
-8
-10
0
50
100
150
200
250
300
ROMP
ROMP算 法 重 构 信 号 10 Recovery Original 8
15
CoSaMP
Recovery Original 10
控 制 迭 代 参 数
3.1算法流程图
3.2 RBAPVS算法性能分析
RBAPVS算法在S = 4 时的迭代次数是最少的,明显优于其他算法; 无论S取何值,重构误差都远低于SAMP 算法和SAMPRB算法,所以该 算法在重构误差的优化上具有很大的优势;
3.2 RBAPVS算法性能分析
在运行时间上,RBAPVS 算法 无论S取何值,都优于原始的 SAMP 算法,但与SAMP-RB 算 法相比较,运行时间略长,但 是差距不大.对RBAPVS 算法中 的S 分别取2,3,4,进行比较, 综合来说,当S 取4 时,所达到 的性能最优.
峰值信噪比的下降顺序仍然都为: IRLS,GBP,OMP,SP,CoSaMP,IHT.
5 结论和展望
结论:介绍了现在已经存在的几种经典贪婪算法,包括 OMP、OMP、 ROMP、GOMP、CoSaMP、StOMP、SP和SAMP 算法。
介绍了一种改进算法并且分析它的性能。 对算法进行仿真和参数分析,抗噪声实验分析。 展望:贪婪算法由于其结构简单、运行时间短而且重构效果比较好等 优点得到了广泛的应用。但是,它不能直接求解原始的优化问题,因 而在重构质量上还是次于基于1范数的算法。后续需要继续学习其它 算法达到更好的重构效果。
4.1.1 二维信号仿真---OMP
4.1.2 二维信号仿真---CoSaMP
4.1.3 二维信号仿真---SP
4.1.4 二维信号仿真---GBP
4.1.5 二维信号仿真---IHT
4.1.6 二维信号仿真---IRLS
4.1.7算法性能分析
重构算法效果参数PSNR(DB)比较表
由分析表格可知,在本实验背景下,相同压缩比条 件下,峰值信噪比的下降顺序为: IRLS,GBP,OMP,SP,CoSaMP,IHT。
谢谢!
溯自适应追踪算法( RBAPVS) .
该算法充分融合了已经存在的ROMP 算法的正则 化思想和SP 算法的回溯思想,并加入了变步长的 设计,使其重构精度和速度与同类算法相比能达 到更好效果.
3.1 RBAPVS算法
稀 疏 度 估 计
正 则 化 二 次 筛 选
3.1 RBAPVS算法
回 溯 思 想 再 筛 选