azu_td_1349475_sip1_... - Arizona Campus Repository

% f3mask.m
% m.file to calculate Fresnel zone plate parameters
% for multi-level/multi-element arrays
% all #'s are in microns
% lens parameter input
clear
clc
lam=input(' what wavelength ');
mask=input(' # of masks');
fna=input(' what fiber na ');
d=input(' what lens diameter');
% constants
k=2*pi/lam;
na=fna*1.6;
f=d/(2*na);
fmin=lam/na/2 A (mask-1);
fnum=f/d;
rmax=d/2;
l=2 A mask;
% performance section
diary
clc
diffeff=(sin(pi/l)/(pi/l)) A 2;
if mask > 1
zm=(((1/lam)*(sqrt(rmax A 2+f A 2)-f)))
a=round(zm)
rm=sqrt(((lam*a) A 2)+(2*f*lam*a));
else
zm=(2/lam)*(sqrt(rmax A 2+f A 2)-f)
a=round(zm)
rm=sqrt(((a*lam/2) A 2)+(a*lam*f));
end
diary
% Gaussian beam propagation (units are mm)
fmm=f/1000;
dmm=d/1000;
lammm=lam/1000;