yangshishuangfengganshe.m

一个杨氏双孔成像干涉模拟

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)