压缩传感重构算法（DCT-OMP）

压缩传感重构算法：稀疏基采用离散余弦变换的omp重构算法

%  本程序实现图像LENA的压缩传感
%  程序作者：沙威，香港大学电气电子工程学系，wsha@eee.hku.hk
%  算法采用正交匹配法，参考文献 Joel A. Tropp and Anna C. Gilbert 
%  Signal Recovery From Random Measurements Via Orthogonal Matching
%  Pursuit，IEEE TRANSACTIONS ON INFORMATION THEORY VOL. 53 NO. 12
%  DECEMBER 2007.
%  该程序没有经过任何优化

%Wavelet_OMP

clc;clear

%  读文件
X=imread（‘u.bmp‘）;
X=double（X）;
[ab]=size（X）;

%  小波变换矩阵生成
ww=dctmtx（a）;

%  小波变换让图像稀疏化（注意该步骤会耗费时间，但是会增大稀疏度）
X1=ww*sparse（X）*ww‘;
X1=full（X1）;

%  随机矩阵生成
M=90;
R=randn（Ma）;

%  测量
Y=R*X1;

%  OMP算法
X2=zeros（ab）;  %  恢复矩阵
for i=1:b  %  列循环       
    rec=omp（Y（:i）Ra）;
    X2（:i）=rec;
end


%  原始图像
figure（1）;
imshow（uint8（X））;
title（‘原始图像‘）;

%  变换图像
figure（2）;
imshow（uint8（X1））;
title（‘小波变换后的图像‘）;

%  压缩传感恢复的图像
figure（3）;
X

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1179  2013-10-25 09:26  DCT-OMP\DCT_OMP.m

     文件       1371  2013-10-24 20:08  DCT-OMP\omp.m

     文件      66614  2009-06-08 18:00  DCT-OMP\u.bmp

     目录          0  2013-10-12 09:52  DCT-OMP

----------- ---------  ---------- -----  ----

                69164                    4