?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