Code covered by the BSD License

Highlights from
Optical_bench

opt_absorption(glass,waveleng... OPT_ABSORPTION - Calculate the absorption coefficient
opt_all_refrindx(wavelength) OPT_REFRINDX refractive index calculations.
opt_aperture(opt_type,opt_arg... OPT_APERTURE Circular aperture, iris
opt_build(file) OPT_BUILD - Builds optical system specified in FILE
opt_cbar(cax)
opt_elem(type) OPT_ELEM - Default opt_elem structure
opt_f2r(f,n) OPT_F2R Thin lens focal width to lens curvature radius.
opt_f_of_r1nr2zc(r1,n,r2,zc) OPT_R2_OF_FNR1_THICK - second radii of curvature given F, N R1 and ZC, thick lens
opt_fcn(opt_type,opt_args) OPT_FCN optical surface with general shape
opt_findstr(H,s) OPT_FINDSTR - find string S in character array H,
opt_fresnel(opt_type,opt_spec) OPT_LENS - Spherical lens single glass.
opt_fresnellens(opt_type,opt_... OPT_FRESNELLENS - Fresnel lens conical groves.
opt_grid(opt_type,opt_args) OPT_GRID - Simple difractive grating
opt_intersection(optelem,rayi... OPT_INTERSECTION - Determine the impact point of an optical
opt_lens(opt_type,opt_spec) OPT_LENS - Spherical lens single glass.
opt_plotoptics(opt_el,plotmod... OPT_PLOTOPTICS - Plot the optical system.
opt_prism(opt_type,opt_spec) OPT_LENS - Spherical lens single glass.
opt_project_direction(e,I,wav... OPT_PROJECT_POINT - Project light from direction E with intensity I throug OPT_EL
opt_project_point(r,I,wavelen... OPT_PROJECT_POINT - Project light from point R with intensity I through optics OPT_EL
opt_r2_of_fnr1_thick(f,n,r1,z... r2 = opt_r2_of_fnr1_thick(f,n,r1,zc)
opt_r2_of_fnr1_thin(f,n,r1) OPT_R2_OF_FNR1_THIN - second radii of curvature given F, N and R1, thin lens
opt_ray() OPT_RAY - RAY creator for OPT_TOOLS
opt_refraction(r_int,ray,opte... OPT_REFRACTION - calculation of optical refraction - Snells law.
opt_refrindx(glass,wavelength) OPT_REFRINDX refractive index calculations.
opt_rot(in_v,rot_v,theta) OPT_ROT - rotate vector IN_V THEDA radians around ROT_V
opt_screen(opt_type,opt_args) OPT_SCREEN - Screen - imaging detector.
opt_slit(opt_type,opt_args) OPT_SLIT - optical slit.
opt_sphereintersection(rl,e_i... OPT_SPHEREINTERSECTION - intersection between a ray and a spherical lens surface
opt_trace(optelements,rays,op... OPT_TRACE - ray tracing through optical systems.
opt_typical_ops() OPT_TYPICAL_OPS - typical standard options for the ray tracer
opt_vini(r1,r2,ng,zc) OPT_VINI - Rear nodal point distance from rear vertex of thick lens
optf_coshyp(r,s_or_n_or_p,arg... OPTF_HYPERBOLC defines a cosh lens surface
point_on_line(r_0,e_l,l); POINT_ON_LINE - calculates the vector to a point
Contents.m
aberation_spherical.m
dGauss_cell.m
grating_fig.m
opt_exempel.m
opt_exempel2.m
opt_exempel3.m
opt_exempel5.m
opt_exempel6.m
opt_exmpl_astigmatic.m
opt_exmpl_chromatic_ab.m
opt_exmpl_coma.m
opt_exmpl_field_curv.m
opt_exmpl_sphere_coma.m
opt_pinhole1.m
rover_cell.m
View all files

from Optical_bench by Bjorn Gustavsson
Optical_bench - a simple ray-tracing tool for optical systems.

Description of opt_sphereintersection

opt_sphereintersection

PURPOSE

OPT_SPHEREINTERSECTION - intersection between a ray and a spherical lens surface

SYNOPSIS

function [rinter,en] = opt_sphereintersection(rl,e_in,rsf,curvature,lensradius,en)

DESCRIPTION

 OPT_SPHEREINTERSECTION - intersection between a ray and a spherical lens surface 
   
 CALLING: 
  [rinter,en] = opt_sphereintersection(rl,e_in,rsf,curvature,lensradius,en)
 INPUT:

CROSS-REFERENCE INFORMATION

This function calls:

point_on_line POINT_ON_LINE - calculates the vector to a point

This function is called by:

opt_plotoptics OPT_PLOTOPTICS - Plot the optical system.

SOURCE CODE

0001 function [rinter,en] = opt_sphereintersection(rl,e_in,rsf,curvature,lensradius,en)
0002 % OPT_SPHEREINTERSECTION - intersection between a ray and a spherical lens surface
0003 %
0004 % CALLING:
0005 %  [rinter,en] = opt_sphereintersection(rl,e_in,rsf,curvature,lensradius,en)
0006 % INPUT:
0007 %
0008 ex = e_in(1);
0009 ey = e_in(2);
0010 ez = e_in(3);
0011 xl = rl(1);
0012 yl = rl(2);
0013 zl = rl(3);
0014 x0 = rsf(1);
0015 y0 = rsf(2);
0016 z0 = rsf(3);
0017 R = curvature;
0018 
0019 % I'd like to say that I did this by myself...
0020 
0021 l(2) = [ 1/2/(ex^2+ey^2+ez^2)*(2*ez*z0-2*zl*ez-2*yl*ey+2*ex*x0+2*ey*y0-2*xl*ex+2*(2*ex^2*zl*z0+2*ez*z0*ex*x0-ex^2*y0^2-ex^2*z0^2-ex^2*yl^2+ex^2*R^2-ex^2*zl^2-ey^2*xl^2-ey^2*z0^2-ey^2*x0^2+ey^2*R^2-ey^2*zl^2-ez^2*xl^2-ez^2*y0^2-ez^2*x0^2-ez^2*yl^2+ez^2*R^2+2*ex^2*yl*y0+2*ey^2*zl*z0+2*ey^2*xl*x0+2*ez^2*xl*x0+2*ez^2*yl*y0-2*ez*z0*yl*ey+2*ez*z0*ey*y0-2*ez*z0*xl*ex+2*zl*ez*yl*ey-2*zl*ez*ex*x0-2*zl*ez*ey*y0+2*zl*ez*xl*ex-2*yl*ey*ex*x0+2*yl*ey*xl*ex+2*ex*x0*ey*y0-2*ey*y0*xl*ex)^(1/2))];
0022 
0023 % ...But who am I trying to fool
0024 
0025 l(1) = [ 1/2/(ex^2+ey^2+ez^2)*(2*ez*z0-2*zl*ez-2*yl*ey+2*ex*x0+2*ey*y0-2*xl*ex-2*(2*ex^2*zl*z0+2*ez*z0*ex*x0-ex^2*y0^2-ex^2*z0^2-ex^2*yl^2+ex^2*R^2-ex^2*zl^2-ey^2*xl^2-ey^2*z0^2-ey^2*x0^2+ey^2*R^2-ey^2*zl^2-ez^2*xl^2-ez^2*y0^2-ez^2*x0^2-ez^2*yl^2+ez^2*R^2+2*ex^2*yl*y0+2*ey^2*zl*z0+2*ey^2*xl*x0+2*ez^2*xl*x0+2*ez^2*yl*y0-2*ez*z0*yl*ey+2*ez*z0*ey*y0-2*ez*z0*xl*ex+2*zl*ez*yl*ey-2*zl*ez*ex*x0-2*zl*ez*ey*y0+2*zl*ez*xl*ex-2*yl*ey*ex*x0+2*yl*ey*xl*ex+2*ex*x0*ey*y0-2*ey*y0*xl*ex)^(1/2))];
0026 
0027 % Where would I be without the symbolic toolbox?
0028 
0029 [qwe,ii] = min(abs(l));
0030 l = l(ii);
0031 
0032 rinter = point_on_line(rl,e_in,l);
0033 
0034 r1 = rinter - rsf;
0035 r2 = cross(r1,en);
0036 
0037 en = [];
0038 if norm(r2)<lensradius
0039   en = r1/norm(r1);
0040 else
0041   rinter = [];
0042 end
0043