%本程序實現圖像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.
%??該程序沒有經過任何優化

%function?Wavelet_OMP

clc;clear

%??讀文件
%?X=imread（‘C:\Users\zkdn\Desktop\MATLABdata\I30\I30.jpg‘）;
%?X=double（X）;
%?[ab]=size（X）;
load?C:\Users\zkdn\Desktop\MATLABdata\I1\I1.mat
X=I1;
background=imopen（Xstrel（‘disk‘15））;%獲取背景信息
X=imsubtract（Xbackground）;%利用函數去除背景

X=double（X）;
[ab]=size（X）;
%??小波變換矩陣生成
ww=dwt2（a‘haar‘）;
%??小波變換讓圖像稀疏化（注意該步驟會耗費時間，但是會增大稀疏度）
X1=ww*sparse（X）*ww‘;
%?X1=X;
X1=full（X1）;

%??隨機矩陣生成
M=256;
R=randn（Ma）;
%?R=randn（a）;
%?R=mapminmax（R0255）;
%?R=ro