Algorithms for separable nonlinear least squares with application to modelling time-resolved spectra

Katharine Mullen, Ivo van Stokkum, Mikas Vengris

Research output: Contribution to JournalArticleAcademicpeer-review

69 Downloads (Pure)

Abstract

The multiexponential analysis problem of fitting kinetic models to time-resolved spectra is often solved using gradient-based algorithms that treat the spectral parameters as conditionally linear. We make a comparison of the two most-applied such algorithms, alternating least squares and variable projection. A numerical study examines computational efficiency and linear approximation standard error estimates. A new derivation of the Fisher information matrix under the full Golub-Pereyra gradient allows a numerical comparison of parameter precision under variable projection variants. Under the criteria of efficiency, quality of standard error estimates and parameter precision, we conclude that the Kaufman variable projection technique performs well, while techniques based on alternating least squares have significant disadvantages for application in the problem domain.
Original languageEnglish
Pages (from-to)201-213
JournalJournal of Global Optimization
Volume38
Issue number2
DOIs
Publication statusPublished - 2007

Fingerprint

Dive into the research topics of 'Algorithms for separable nonlinear least squares with application to modelling time-resolved spectra'. Together they form a unique fingerprint.

Cite this