function?CoSa

clc;clear

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

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

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

%??隨機矩陣生成
M=190;
R=randn（Ma）;%傳感矩陣

%??測量
Y=R*X1;

%??OMP算法
X2=zeros（ab）;??%??恢復(fù)矩陣

tic
for?i=1:b??%??列循環(huán)???????
????rec=CS_CoSaMP（Y（:i）Ra）;
????X2（:i）=rec;
end
toc

%??原始圖像
figure（1）;
imshow（uint8（X））;
title（‘原始圖像‘）;
%??變換圖像
figure（2）;
imshow（uint8（X1））;
title（‘小波變換后的圖像‘）;
%??壓縮傳感恢復(fù)的圖像
figure（3）;
X3=ww‘*sparse（X2）*ww;??%??小波反變換
X3=full（X3）;
imshow（uint8（X3））;
title（‘恢復(fù)的圖像‘）;
%??誤差（PSNR）
errorx=sum（sum（abs（X3-X）.^2））;????????%??MSE誤差
%psnr=10*log10（255*255/（errorx/a/b））???%??PSNR

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????830??2018-03-01?17:03??cs\CoSa.m

?????文件???????1220??2018-03-12?15:11??cs\cs.m

?????文件????????917??2018-03-12?15:51??cs\cs64.m

?????文件???????2170??2018-03-01?15:34??cs\CS_CoSaMP.m

?????文件???????1115??2018-02-28?20:11??cs\DWT.m

?????文件???????3652??2018-03-08?09:40??cs\lena128.jpg

?????文件??????66614??2009-06-08?18:00??cs\lena256.bmp

?????文件???????1388??2018-03-08?14:37??cs\lena64.jpg

?????文件???????1346??2018-03-02?10:19??cs\omp.m

?????文件????????842??2018-03-05?15:55??cs\PartHadamardMtx.m

?????文件????????444??2018-03-07?16:26??cs\SparseRandom.m

?????文件????????444??2018-03-07?16:30??cs\SparseRandomMtx.m

?????目錄??????????0??2018-03-12?17:03??cs

-----------?---------??----------?-----??----

????????????????80982????????????????????13