www.pudn.com > work.rar > tiaoshi.m
%¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡ picturesiz=size(p1); N=picturesiz(1); M=picturesiz(2); cm=0.001;um=1e-6; hx=0.06*cm;hy=0.06*cm; dhx=hx/M;dhy=hy/N; x=dhx.*(ones(N,1)*[-M/2:M/2-1]); y=(dhy.*(ones(M,1)*[-N/2:N/2-1]))'; z=0.00002; w=0.63*um; k=2*pi/w; %R1=rand(1,1);p1=p1.*exp(j*2*pi*R1); [f,dx1,dy1,x1,y1]=fresnell(p1,M,N,dhx,dhy,z,w); sit=pi/6; ur=exp(-j*k*y*sin(sit)); f1=f+ur; f1=abs(f1).^2; f1=mat2gray(f1); figure;imshow(f1);