azu_td_1349475_sip1_... - Arizona Campus Repository

theta=4*lammm/pi/dmm;
z=0:100;
dg=dmm*sqrt(1+((theta.*z)/dmm). A 2);
div=(dg)./dmm;
mxdiv=max(div);
73
% output performance
diary
disp('
Lens Parameters')
disp("')
disp(sprintf(' wavelength %5.4f um'.lam))
disp(sprintf(' diameter %7.3f um',rm*2))
disp(sprintf(" NA %5.4f ',na))
disp(sprintf(' focal length %6.4f urn ',f))
disp(sprintf(' F# %4.1f \fnum))
disp(sprintf(' #masks %2.0f '.mask))
disp(sprintf(' diffrac eff. %4.3f \ diffeff))
zmax=( 1 /lam)*(sqrt(rmax A 2+f A 2)-f);
for n=1 :mask
int=a*n;
fmin=lam/na/2 A (n);
disp(sprintf('mask %2.0f has %3.0f zones w/ min feat. %6.4f um',n,int,fmin))
end;
dispC')
disp(sprintf('beam divergence ratio at 100mm is %6.3f',mxdiv))
disp('')
% mask generation section
if mask > 1
for n=1 :mask
disp(sprintf(the radii for mask %2.0f are \n));
a=round(zmax)*n:-1:1;
r=sqrt(((lam*a/(2 A (n-1))). A 2)+(2Tlam*a/(2. A (n-1))))
end;
else
disp(sprintf('the radi for mask %2.0f are ',n));
zmax=(2/lam)*(sqrt(rmax A 2+f A 2)-f);
a=round(zmax)*n:-1:1;
r=sqrt(((a*lam/2). A 2)+(a*lam*f))
end;
diary