omp 用小波变换作为稀疏基

用小波变换作为稀疏基，采取OMP算法将图像重建恢复，由于算法计算量大会导致成像时间过长，程序用改进的分块处理缩短了时间，

 clc;
clear all;

for p=0.3:0.1:0.9
tic;



im=imread（‘lena256.jpg‘）;
im=im（::1）;
[m n]=size（im）;
im2=zeros（mn）;
bm=16;
bn=16;
d=bm*bn;
N=round（d*p）;
bx=ceil（m/bm）;
by=ceil（n/bn）;

rx=mod（mbm）;
ry=mod（nbn）;
imc=zeros（bx*bmby*bn）;

for i=1:bx*bm
    for j=1:by*bn
        if i<=m&&j<=n
            imc（ij）=im（ij）;
        else imc（ij）=0;
        end
    end
end
 


for x=1:bx;
    for y=1:by;
        b=imc（（1+（x-1）*bm）:（x*bm）（1+（y-1）*bn）:（y*bn））;
        xin=reshape（bd1）;
        xin=double（xin）;
       
        A=dctmtx（d）;
        Phi=randn（Nd）;   
        T=Phi*inv（A）;
        s=Phi*xin;                                   
 
%        [x2 iter]=omp（sT）;
%        [x2 iter]=fun_CMP（sT）;
%        x2=GPSR_BB（sT50‘verbose‘0）;
  

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      33951  2012-07-15 22:25  omp\lena256.jpg

     文件       1562  2012-08-02 15:34  omp\main.m

     文件        918  2012-07-16 08:54  omp\omp.m

     目录          0  2012-08-20 09:53  omp

----------- ---------  ---------- -----  ----

                36431                    4