clear?all;clc;
cm?=?1e-2;
mm?=?1e-3;
um?=?1e-6;
nm?=?1e-9;%單位準備部分%

Lamda=638*nm;
M=512;?N=512;
dx=10.0*um;?dy=10.0*um;?


x?=?linspace（-255256N）;???%?衍射屏和觀察屏
screen?=?zeros（MN）;??????%?狹縫
[mn]=meshgrid（x）;
l=30?;h=50;
a=10;
D=find（（abs（m）<=l）&（abs（n）<=h）&（abs（m）>=a））;
screen（D）?=?1;??%?長為2h，高為l的矩形
figureimshow（screen[]）;

?for?Dist=5*cm:1*cm:7*cm
x=linspace（-M*dx/2（M-1）*dx/2N）;
y=linspace（-N*dy/2（N-1）*dy/2M）;