clear
lam=500e-9 ;a=2e-3;
D=1;
ym=5*lam*D/a;xs=ym;
n=101;
ys=linspace(-ym,ym,n);
for i=1:n
r1=sqrt((ys(i)-a/2).^2+D^2);
r2=sqrt((ys(i)+a/2).^2+D^2);
phi=2*pi*(r2-r1)./lam;
B(i,:)=sum(4*cos(phi/2).^2);
end
N=255;
Br=(B/4.0)*N ;
subplot(1,2,1) ;
image(xs,ys,Br);
colormap(gray(N));
subplot(1,2,2) ;
plot(B,ys)
lam=500e-9 ;a=2e-3;
D=1;
ym=5*lam*D/a;xs=ym;
n=101;
ys=linspace(-ym,ym,n);
for i=1:n
r1=sqrt((ys(i)-a/2).^2+D^2);
r2=sqrt((ys(i)+a/2).^2+D^2);
phi=2*pi*(r2-r1)./lam;
B(i,:)=sum(4*cos(phi/2).^2);
end
N=255;
Br=(B/4.0)*N ;
subplot(1,2,1) ;
image(xs,ys,Br);
colormap(gray(N));
subplot(1,2,2) ;
plot(B,ys)