%…………………………………………………………
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);