一种渐进式正则化自适应匹配追踪方法【专利摘要】本发明属于图像信号处理领域,具体涉及一种用于信号或者图像重建的渐进式正则化自适应匹配追踪方法。本发明包括:初始参数定义,利用余量与传感矩阵的每一列的内积计算相关系数;对J中索引值对应原子的相关系数进行正则化;更新支撑集,并进行余量更新等。本发明更合理的设计了稀疏度的迭代步长,保证了迭代能够渐进式地逼近信号的真实稀疏度,提高了迭代的收敛性及信号重建的质量。【专利说明】一种渐进式正则化自适应匹配追踪方法【
技术领域:
】[0001]本发明属于图像信号处理领域,具体涉及一种用于信号或者图像重建的渐进式正则化自适应匹配追踪方法。【
背景技术:
】[0002]近几年出现的压缩传感理论(如CANDESE.Compressivesampling[A].ProceedingsoftheInternationalCongressofMathematicians[C].Madrid,Spain,2006,3:1433-1452;D0N0H0D.L.Compressedsensing[J]·IEEETrans.OnInformationTheory.2006,52(4):1289-1306.),将信号的采样与压缩过程合二为一,在进行数据采集时,只需要获取很少的观测值,然后利用重构算法就能够精确恢复出原信号。该理论一经提出,就引起了学术界的轰动,国内外众多学者纷纷开展了相关领域的研究。[0003]信号的重建是压缩传感理论的关键内容,D0N0H0DL,ELADM,andTEMLYAK0VVN在((Stablerecoveryofsparseovercompleterepresentationsinthepresenceofnoise))IEEETransactionsonInformationTheory,2006,52(I):6-18中指出匹配追踪类方法能够有效地应用于压缩传感的信号重建并具有一定的稳定性。在此基础上引入的正交化和正则化过程,保证了迭代的最优性,减少了迭代次数。而自适应匹配追踪算法(THONGTDo,GANLu,NGUYENetal.Sparsityadaptivematchingpursuitalgorithmforpracticalcompressedsensing.AsilomarConferenceonSignals,Systems,andComputers,PacificGrove,California,2008,10:581-587.)可以在信号稀疏度未知的情况下获得较好的重建效果,且具有较快的速度。刘亚新,赵瑞珍,胡绍海等在电子与信息学报的《用于压缩感知信号重建的正则化自适应匹配追踪算法》中将上述各种思想融合起来提出了正则化自适应匹配追踪算法。正则化自适应匹配追踪算法增强了信号重建的理论依据同时不需要预先估计信号的稀疏度,但其迭代步长的设置不甚合理,跨越式的步长迭代很难保证信号重建过程的收敛性。【
发明内容】[0004]本发明的目的在于提供一种可提高信号稀疏度迭代的收敛性,并获得高质量的重建信号的渐进式正则化自适应匹配追踪方法。[0005]本发明的目的是这样实现的:[0006](I)初始参数定义,定义测量值为y,重建信号为I初始余量:Ttl=y,传感矩阵为O,初始步长size古O,初始稀疏度K。=size,索引值集合Λ=0,J=0,Λ=0,支撑集ΦΛ,迭代次数η=I,迭代阶段Stage=I,阈值I为E1、阈值2为ε2、阈值3为ε3;[0007](2)利用余量r与传感矩阵Θ的每一列的内积计算相关系数,在相关系数中寻找Ktl个最大值对应的索引值存入J中;[0008](3)对J中索引值对应原子的相关系数进行正则化,并将正则化结果存入集合Jtl中;[0009](4)更新支撑集ΦΛ,应用最小二乘法估算得到信号的估计值;,并进行余量更新;[0010](5)若IIw1II(ε:,执行步骤6,否则,令η=η+1,重新执行步骤2;[0011](6)若IIrnM2Sε2,则停止迭代,否则执行步骤7;【权利要求】1.一种渐进式正则化自适应匹配追踪方法,其特征在于:(1)初始参数定义,定义测量值为y,重建信号为i,初始余量a=y,传感矩阵为?,初始步长size关0,初始稀疏度KQ=size,索引值集合A=0,J=0,=0,支撑集①A,迭代次数n=1,迭代阶段stage=l,阈值1为£工、阈值2为e2、阈值3为e3;(2)利用余量r与传感矩阵?的每一列的内积计算相关系数,在相关系数中寻找&个最大值对应的索引值存入J中;(3)对J中索引值对应原子的相关系数进行正则化,并将正则化结果存入集合1中;(4)更新支撑集OA,应用最小二乘法估算得到信号的估计值Y,并进行余量更新;(5)若llHill(h,执行步骤6,否则,令n=n+l,重新执行步骤2;(6)若||rn||2≤e2,则停止迭代,否则执行步骤7;(7)若||rn||2≤e3,则令,size^,K0=K0+size,n=n+1,stage=stage+l,执行步骤2,否则令Kq=K0+size,n=n+1,stage=stage+l,执行步骤2。【文档编号】G06T11/00GK103489207SQ201310452181【公开日】2014年1月1日申请日期:2013年9月29日优先权日:2013年9月29日【发明者】卞红雨,吴菊,张志刚,孙慧娟申请人:哈尔滨工程大学