MATLAB模拟菲涅尔半波带程序.doc
需积分: 45 61 浏览量
2020-04-18
14:58:01
上传
评论 1
收藏 12KB DOC 举报
clear;
tic
lam=600e-6; %波长 600nm
R=3; %菲涅尔波带片半径
为 3mm
f=1000; %焦距 1m
ym=R;xm=R;theta=0:pi/100:2*pi;
x=linspace(-xm,xm,1001); %设置程序中 x 分为
1001 等分
for n=1:1001
k=fix(x(n).^2./(lam.*f)); %求各点半波数点
if mod(k,2)==1 %判断半波带数奇偶，
为奇数则涂黑
x1(:,n)=x(n)*cos(theta);y1(:,n)=x(n)*sin(theta);
end
end
figure(1)
plot(x1,y1,'-k')
toc
资源评论
