The nonlinear transformations incurred by the rays in an optical system can be suitably described by matrices to any desired order of approximation. In systems composed of uniform refractive-index elements each individual ray refraction or translation has an associated matrix, and a succession of transformations corresponds to the product of the respective matrices. A general method is described to find the matrix coefficients for translation and surface refraction, irrespective of the surface shape or the order of approximation. The choice of coordinates is unusual, as the orientation of the ray is characterized by the direction cosines rather than by the slopes; this is shown to greatly simplify and generalize coefficient calculation. Two examples are shown in order to demonstrate the power of the method: The first is the determination of seventh-order coefficients for spherical surfaces, and the second is the determination of third-order coefficients for a toroidal surface.
© 1999 Optical Society of America
(080.2730) Geometric optics : Matrix methods in paraxial optics
(220.1010) Optical design and fabrication : Aberrations (global)
(220.2740) Optical design and fabrication : Geometric optical design
José B. Almeida, "General method for the determination of matrix coefficients for high-order optical system modeling," J. Opt. Soc. Am. A 16, 596-601 (1999)