zaixian.m

来自「基于幻方加密算发的全数字全息图制作与加密与解密」· M 代码 · 共 22 行

mim=setstr(mim);
 if strcmp(mima,mim)==0
    handle=errordlg('解密错误，不能正确再现，请重新解密：');
 else
     wei=size(ur);fff=fff(1:wei(1),1:wei(2));
rp1=fff.*conj(ur);
%……………………………………高通滤波……………………………………………………
[a,b]=size(rp1);
n=4;d0=10;n1=(a/2);n2=(b/2);
for u=1:a
     for t=1:b
         d=sqrt((u-n1)^2+(t-n2)^2);
         h=1/(2+0.614*(d0/d)^(2*n));
         rp1(u,t)=h*rp1(u,t);
     end
end 
%………………………………………………………………………………
[f2,dx2,dy2,x2,y2]=fresnell(rp1,M,N,dx1,dy1,z,w);
f2=abs(f2).^2;
f2=mat2gray(f2);
figure;imshow(f2);
 end