|
579. LPT: Linear Predictive Single Value Decomposition of
Decaying Sinusoids for Spectral Analysis in NMR
by David Cowburn and Jewel Chaudhary, The Rockefeller University, New York, New York 10021 This program is used for the determination of a complex time-domain signal as the summation of a number of exponentially damped sinusoids. The frequencies and damping factors are determined via linear prediction. The linear prediction equations are solved with a linear least-squares procedure via diagonalization of the "normal equations" (indirect singular value decomposition). The frequencies and damping factors are determined by polynomial factorization. The amplitude and phases are obtained by a second linear least-squares procedure. There are two versions of this program, each containing a different driver. Both drivers expect an input data file of double complex data. One version requests the needed parameters which the user supplies, while the other version expects the parameters to be entered in a parameter file. FORTRAN 77 (SUN) Lines of Code: 3855 |