clear;
lam=600e-6;
R=3;
f=1000;
ym=R;
xm=R;
y=linspace(-ym,ym,1001);
x=linspace(-xm,xm,1001);
for m=1:1001
for n=1:1001
p=sqrt(x(m).^2+y(n).^2);
k=fix(p.^2./(lam.*f));
if p>R
I(m,n)=0;
else
if mod(k,2)==1
I(m,n)=0;
else
I(m,n)=1;
end
end
end
end
colormap(gray(255))
N=255;
Ir=I*N;
image(Ir);
评论0