论坛 › MATLAB论坛 › 菲涅尔波带片的设计

李秋锐 发表于 2015-2-1 19:26

菲涅尔波带片的设计

%菲涅尔波带片的设计，本程序预计运行耗时25s
r=6.3;
points=10000;
lamb=0.5*10E-3;%光波长
R=1000;%物距
b=250;%像距
X=linspace(-r,r,points);
Y=linspace(r,-r,points);
A=ones(points,points);
for m=1:points
for n=1:points
if((sqrt(X(n)^2+Y(m)^2+b^2)-b)/(lamb/2)<30)
if(mod(floor((sqrt(X(n)^2+Y(m)^2+b^2)-b)/(lamb/2)),2)==0) A(m,n)=0;
end
end
end
end
figure('color',);
imagesc(A,);
colormap(gray);
set(gca,'xtick',[],'ytick',[]);
axis square

查看完整版本: 菲涅尔波带片的设计