%??本程序?qū)崿F(xiàn)圖像LENA的壓縮傳感
%??程序作者：沙威，香港大學(xué)電氣電子工程學(xué)系，wsha@eee.hku.hk
%??算法采用正交匹配法，參考文獻(xiàn)?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.
%??該程序沒(méi)有經(jīng)過(guò)任何優(yōu)化

%Wavelet_OMP

clc;clear

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

%??小波變換矩陣生成
ww=dctmtx（a）;

%??小波變換讓圖像稀疏化（注意該步驟會(huì)耗費(fèi)時(shí)間，但是會(huì)增大稀疏度）
X1=ww*sparse（X）*ww‘;
X1=full（X1）;

%??隨機(jī)矩陣生成
M=90;
R=randn（Ma）;

%??測(cè)量
Y=R*X1;

%??OMP算法
X2=zeros（ab）;??%??恢復(fù)矩陣
for?i=1:b??%??列循環(huán)???????
????rec=omp（Y（:i）Ra）;
????X2（:i）=rec;
end


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

%??變換圖像
figure（2）;
imshow（uint8（X1））;
title（‘小波變換后的圖像‘）;

%??壓縮傳感恢復(fù)的圖像
figure（3）;
X

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----

?????文件???????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