clear;
cm=0.01;um=1e-6;
N=128;
rx=0*cm*ones(N);ry=0*cm*ones(N);rz=3*ones(N);
ox=1*cm*ones(N);oy=1*cm*ones(N);oz=2*ones(N);
hx=2*cm;hy=2*cm;
dhx=hx/N;dhy=hy/N;
w=0.63*um;
k=2*pi/w;
x=dhx.*(ones(N,1)*[-N/2:N/2-1]);
y=(dhy.*(ones(N,1)*[-N/2:N/2-1]))';
z=0*ones(N);
d1=sqrt((rx-x).^2+(ry-y).^2+(rz-z).^2);
d2=sqrt((ox-x).^2+(oy-y).^2+(oz-z).^2);
ur=exp(i*k*d1)./d1;
uo=exp(i*k*d2)./d2;
p=abs(ur+uo).^2;
p=mat2gray(p);
figure;
imshow(p);
p1=p.*conj(ur);
clear;clc;
cm=0.01;um=1e-6;
N=128;
hx=1*cm;hy=1*cm;
dhx=hx/N;dhy=hy/N;
x=dhx.*(ones(N,1)*[-N/2:N/2-1]);
y=(dhy.*(ones(N,1)*[-N/2:N/2-1]))';
p1=imread('mz.jpg');
figure;imshow(p1);
z=4;
w=0.63*um;
k=2*pi/w;
[f1,dx1,dy1,x1,y1]=fresnell(p1,N,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);
rp1=f1.*conj(ur);
a=0.002;b=0.3;
rp1=pure(rp1,a,b,x1,y1);
[f2,dx2,dy2,x2,y2]=fresnell(p1,N,N,dhx,dhy,z,w);
f2=abs(f2)^2;
f2=mat2gray(f2);
figure;imshow(f2);