matlab模拟方孔衍射,matlab 矩孔衍射程序

下面的是矩孔衍射的程序，我每次在网吧做的时候程序运行总是得不到图形，请哪位方便的帮我调试运行一下！！！谢谢各位啦！！！

请把运行结果得到图片发到我的邮箱里，谢谢了！！！

我的邮箱：lxiuyue@126.com

程序一：

clear

lam=500e-9

a=1e-3

b=1e-3

f=1

m=500;

ym=8000*lam*f;

ys=linspace(-ym,ym,m);

xs=ys;

n=255;

col(m,m)=0

for i=1:m

sinth1=xs(i)/sqrt(xs(i)^2+f^2);

sinth2=ys./sqrt(ys.^2+f^2);

angleA=pi*a*sinth1/lam;

angleB=pi*b*sinth2./lam;

col(:,i)=(sin(angleA).^2.*sin(angleB).^2.*5000./(angleA.^2.*angleB.^2));

end

subplot(1,2,1)

image(xs,ys,col)

colormap(gray(n))

subplot(1,2,2)

plot(xs,col(m/2,:))

程序二：

clear

lam=500e-9<