压缩感知算法

压缩感知去算法代码并由参数评估算法，实物图和效果图

function G = BGI（gwNrc）
I1I2=zeros（rc）;
I1=zeros（rc）;
I2=0;
i = sqrt（-1）;
tic
for mn=1:N
g0=（normrnd（01rc））.*exp（1i*（randn（rc））*2*pi）;
lamde=633e-9;           			%波长
boshuk=2*pi/lamde;   
d=350; 
A=zeros（rc）;
for m=1:r
    for n=1:c
 A（mn）=exp（（-1i*boshuk*（m^2+n^2））/（2*d*1e-3））;
   end
end
z=zeros（rc）;
z1=zeros（rc）;
%Ir的矩阵
z（1:r1:c）=g0.*A;%见专利[0047]中对应的公式光源强度和光源散斑的乘积
z1（1:r1:c）=conj（z）.*z; %Ir矩阵w的虚部和实部的平方差
t1=zeros（rc）;
t1=（double（gw））.*z1（1:r1:c）;%t1即密文
R0=0;
  for m=1:r
    for n=1:c
        R0=R0+t1（mn）;%密文对应矩阵中所有点的和
    end
  end
  I1I2=I1I2+z1.*R0;
  I1=I1+z1;
  I2=I2+R0;
end
G=I1I2/N-（（I1/N）.*（I2/N））;%关联计算公式
toc

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

     文件       1118  2016-01-06 13:47  Demo_CS_BP\Demo_CS_GBP.m

     文件       1065  2015-08-26 11:21  Demo_CS_BP\DWT.m

     文件        740  2016-01-05 16:59  Demo_CS_BP\BGI.m

     文件       4289  2015-06-04 14:17  Demo_CS_BP\ssim.m

     文件       6145  2016-01-06 13:40  Demo_CS_BP\cs_gbp.m

     文件       1092  2016-01-06 13:22  Demo_CS_BP\main.m

     文件        720  2015-06-04 09:13  Demo_CS_BP\Image\PZL.png

     文件        555  2014-08-28 08:47  Demo_CS_BP\Image\G1.png

     文件      66614  2009-06-08 18:00  Demo_CS_BP\Image\lena256.bmp

     文件     263224  2015-04-09 20:12  Demo_CS_BP\Image\lena512.bmp

     目录          0  2016-01-06 11:06  Demo_CS_BP\result

     目录          0  2016-01-06 11:06  Demo_CS_BP\Image

     目录          0  2015-01-12 17:15  Demo_CS_BP

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

               345562                    13