二维图像OMP算法代码

使用matlab代码实现的二维图像压缩感知OMP算法实现过程

function hat_x=cs_omp（yT_Matm）
% y=T_Mat*x T_Mat is n-by-m
% y - measurements
% T_Mat - combination of random matrix and sparse representation basis
% m - size of the original signal
% the sparsity is length（y）/4

n=length（y）;
s=floor（n/4）;                                     %  测量值维数
hat_x=zeros（1m）;                                 %  待重构的谱域（变换域）向量                     
Aug_t=[];                                         %  增量矩阵（初始值为空矩阵）
r_n=y;                                            %  残差值 

for times=1:s;                                  %  迭代次数（稀疏度是测量的1/4）

    product=abs（T_Mat‘*r_n）;

    [valpos]=max（product）;                       %  最大投影系数对应的位置
    Aug_t=[Aug_tT_Mat（:pos）];                   %  矩阵扩充
    T_Mat（:pos）=zeros（n1）;                      %  选中的列置零（实质上应该去掉，为了简单将其置零）
    aug_x=（Aug_t‘*Aug_t）^（-1）*Aug_t‘*y;           %  最小二乘使残差最小
    r_n=y-Aug_t*aug_x;                            %  残差
    pos_array（times）=pos;                         %  纪录最大投影系数的位置

end
hat_x（pos_array）=aug_x;                           %  重构的向量 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2015-10-21 16:25  OMP\
     文件      262308  2015-10-20 09:52  OMP\barbara.pgm
     文件        1186  2015-10-20 09:52  OMP\cs_omp.m
     文件         164  2015-10-20 09:52  OMP\dct_basis.m
     文件        1931  2015-10-21 19:56  OMP\Demo_CS_OMP.m
     文件         599  2015-10-20 09:52  OMP\GenerateProjection.m
     文件      262299  2015-10-20 09:52  OMP\goldhill.pgm
     文件      262302  2015-10-20 09:52  OMP\lenna.pgm
     文件      262218  2015-10-20 09:52  OMP\mandrill.pgm
     文件         857  2015-10-20 09:52  OMP\omp.m
     文件      262230  2015-10-20 09:52  OMP\peppers.pgm
     文件     2014489  2015-10-20 09:52  OMP\projections.1.mat
     文件        1329  2015-10-20 09:52  OMP\PSNR.m
     文件        1021  2015-10-20 09:52  OMP\test_omp.m
     文件        1473  2015-10-20 09:52  OMP\Untitled3.m