%??本程序實現圖像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