%??讀文件
X=imread（‘lena.bmp‘）;
X=double（X）;
[ab]=size（X）;
size_kuai=16*4;
X2=zeros（size_kuai）;??%??恢復矩陣
X3=zeros（ab）;??%??恢復矩陣
%??小波變換矩陣生成
ww=DWT（size_kuai）;
%??隨機矩陣生成
M=80;
R=randn（Msize_kuai）;
tic
for?i_x=1:ceil（a/size_kuai）
????for?i_y=1:ceil（b/size_kuai）
????????XX=X（（i_x-1）*size_kuai+1:i_x*size_kuai（i_y-1）*size_kuai+1:i_y*size_kuai）;
????????%??小波變換讓圖像稀疏化（注意該步驟會耗費時間，但是會增大稀疏度）
????????X1=ww*sparse（XX）*ww‘;
????????X1=full（X1）;
????????%??測量
????????Y=R*X1;
????????%??OMP算法
????????for?i=1:size_kuai??%??列循環(huán)???????
????????????rec=omp_fenkuai（Y（:i）Rsize_kuai）;
????????????X2（:i）=rec;
????????end
????????X3（（i_x-1）*size_kuai+1:i_x*size_kuai（i_y-1）*size_kuai+1:i_y*size_kuai）=ww‘*sparse（X2）*ww;??%??小波反變換
????end
end
X3=full（X3）;
use_time=toc

%??原始圖像
figure（1）;
imshow（uint8（X））;
title（‘原始圖像‘）;
%??壓縮傳感恢復的圖像
figure（2）;
imshow（uint8（X3））;
title（‘分塊恢復的圖像‘）;

%??誤差（PSNR）
errorx=sum（sum（abs（X3-X）.^2））;????????%??MSE誤差
psnr=10*log10（255*255/（errorx/a/b））??%??PSNR

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

?????文件???????1115??2013-10-24?21:02??gezhongsuanfa\block_omp\blk_omp.m

?????文件????????821??2013-10-18?10:53??gezhongsuanfa\block_omp\DWT.m

?????文件??????66614??2004-06-30?00:47??gezhongsuanfa\block_omp\lena.bmp

?????文件???????2217??2013-10-18?10:54??gezhongsuanfa\block_omp\omp_fenkuai.m

?????文件????????519??2013-10-14?16:00??gezhongsuanfa\BP\BP.m

?????文件??????66614??2004-06-30?00:47??gezhongsuanfa\BP\lena.bmp

?????文件??????54596??2006-07-30?14:13??gezhongsuanfa\BP\pdco.m

?????文件??????10186??2006-07-30?14:13??gezhongsuanfa\BP\pdcoSet.m

?????文件???????4483??2007-08-18?20:54??gezhongsuanfa\BP\SolveBP.m

?????文件??????66614??2009-06-08?18:00??gezhongsuanfa\OMP\lena.bmp

?????文件????????516??2013-10-14?16:10??gezhongsuanfa\OMP\omp.m

?????文件???????5451??2013-10-14?16:11??gezhongsuanfa\OMP\SolveOMP.m

?????文件????????525??2013-10-14?16:07??gezhongsuanfa\StOMP\fdrthresh.m

?????文件??????66614??2004-06-30?00:47??gezhongsuanfa\StOMP\lena.bmp

?????文件????????475??2013-10-14?16:07??gezhongsuanfa\StOMP\reverse.m

?????文件???????4889??2013-10-14?16:06??gezhongsuanfa\StOMP\SolveStOMP.m

?????文件????????532??2013-10-14?16:06??gezhongsuanfa\StOMP\StOMP.m

?????目錄??????????0??2013-10-18?11:03??gezhongsuanfa\block_omp

?????目錄??????????0??2013-10-18?11:00??gezhongsuanfa\BP

?????目錄??????????0??2013-10-24?20:41??gezhongsuanfa\OMP

?????目錄??????????0??2013-10-18?10:59??gezhongsuanfa\StOMP

?????目錄??????????0??2013-10-18?11:03??gezhongsuanfa

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

???????????????352781????????????????????22