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