Abstract
This paper outlines a method for the determination of the unknown wave-front aberration function of an optical system from noisy measurements of the corresponding point-spread function. The problem is cast as a nonlinear least-squares estimation problem for the values of the wave-front aberration function at N points over the slit aperture, from measurements of the point-spread function at M points with . Newton’s method is used to replace the nonlinear minimization problem with a sequence of linear problems. Each such problem requires the inversion of the Hessian matrix of the error metric that is shown to be both singular (with rank ) and ill-conditioned. To overcome singularity, the pseudoinverse is used; to overcome ill-conditioning the pseudoinverse is calculated using singular value decomposition and the singular values then filtered. Attention is drawn to difficulties such as nonuniqueness, sensitivity of algorithms to initial guess, etc.; the ancillary mathematical details being set out in appendices. Some illustrative numerical results are presented and analyzed.
© 1980 Optical Society of America
Full Article | PDF ArticleMore Like This
Richard Barakat and Barbara H. Sandler
J. Opt. Soc. Am. A 9(10) 1715-1723 (1992)
W.H. Southwell
J. Opt. Soc. Am. 70(8) 998-1006 (1980)
Richard Barakat
J. Opt. Soc. Am. 70(6) 739-742 (1980)