Show simple item record

dc.contributor.advisor Embree, Mark
dc.contributor.advisor Cox, Steven J.
dc.creatorHokanson, Jeffrey
dc.date.accessioned 2014-09-11T16:24:58Z
dc.date.available 2014-09-11T16:24:58Z
dc.date.created 2013-12
dc.date.issued 2013-12-06
dc.date.submitted December 2013
dc.identifier.citation Hokanson, Jeffrey. "Numerically Stable and Statistically Efficient Algorithms for Large Scale Exponential Fitting." (2013) Diss., Rice University. https://hdl.handle.net/1911/77161.
dc.identifier.urihttps://hdl.handle.net/1911/77161
dc.description.abstract The exponential fitting problem appears in diverse applications such as magnetic resonance spectroscopy, mechanical resonance, chemical reactions, system identification, and radioactive decay. In each application, the exponential fitting problem decomposes measurements into a sum of exponentials with complex coefficients plus noise. Although exponential fitting algorithms have existed since the invention of Prony's Method in 1795, the modern challenge is to build algorithms that stably recover statistically optimal estimates of these complex coefficients while using millions of measurements in the presence of noise. Existing variants of Prony's Method prove either too expensive, most scaling cubically in the number of measurements, or too unstable. Nonlinear least squares methods scale linearly in the number of measurements, but require well-chosen initial estimates lest these methods converge slowly or find a spurious local minimum. We provide an analysis connecting the many variants of Prony's Method that have been developed in different fields over the past 200 years. This provides a unified framework that extends our understanding of the numerical and statistical properties of these algorithms. We also provide two new algorithms for exponential fitting that overcome several practical obstacles. The first algorithm is a modification of Prony's Method that can recover a few exponential coefficients from measurements containing thousands of exponentials, scaling linearly in the number of measurements. The second algorithm compresses measurements onto a subspace that minimizes the covariance of the resulting estimates and then recovers the exponential coefficients using an existing nonlinear least squares algorithm restricted to this subspace. Numerical experiments suggest that small compression spaces can be effective; typically we need fewer than 20 compressed measurements per exponential to recover the parameters with 90% efficiency. We demonstrate the efficacy of this approach by applying these algorithms to examples from magnetic resonance spectroscopy and mechanical vibration. Finally, we use these new algorithms to help answer outstanding questions about damping in mechanical systems. We place a steel string inside vacuum chamber and record the free response at multiple pressures. Analyzing these measurements with our new algorithms, we recover eigenvalue estimates as a function of pressure that illuminate the mechanism behind damping.
dc.format.mimetype application/pdf
dc.language.iso eng
dc.subjectExponential fitting
Exponential sums
Exponentially damped sinusoids
Nonlinear least squares
Parameter estimation
Prony's method
Damping
Compression
Compressed nonlinear least squares
Large scale nonlinear least squares
Big data
dc.title Numerically Stable and Statistically Efficient Algorithms for Large Scale Exponential Fitting
dc.contributor.committeeMember Antoulas, Athanasios C.
dc.contributor.committeeMember Heinkenschloss, Matthias
dc.date.updated 2014-09-11T16:24:58Z
dc.type.genre Thesis
dc.type.material Text
thesis.degree.department Computational and Applied Mathematics
thesis.degree.discipline Engineering
thesis.degree.grantor Rice University
thesis.degree.level Doctoral
thesis.degree.name Doctor of Philosophy


Files in this item

Thumbnail

This item appears in the following Collection(s)

Show simple item record