matlab 光学衍射模拟

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）;