www.pudn.com > work.rar > zaixian.m

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