菲涅尔波带片的设计

[复制链接]

字体大小: 正常放大

李秋锐

14 主题	14 听众	76 积分

升级 74.74%

TA的每日心情

	奋斗 2015-3-29 18:44

签到天数: 17 天

[LV.4]偶尔看看III

电梯直达

1^#

发表于 2015-2-1 19:26 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

%菲涅尔波带片的设计，本程序预计运行耗时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',[0,0,0]);
imagesc(A,[0,1]);
colormap(gray);
set(gca,'xtick',[],'ytick',[]);
axis square

zan