Estimation of parameters of exponentially damped sinusoids using fast maximum likelihood estimation with application to NMR spectroscopy data

We present fast maximum likelihood (FML) estimation of parameters of multiple exponentially damped sinusoids. The FML algorithm was motivated by the desire to analyze data that have many closely spaced components, such as the NMR spectroscopy data of human blood plasma. The computational efficiency of FML lies in reducing the multidimensional search involved in ML estimation into multiple 1-D searches. This is achieved by using our knowledge of the shape of the compressed likelihood function (CLF) in the parameter space. The proposed FML algorithm is an iterative method that decomposes the original data into its constituent signal components and estimates the parameters of the individual components efficiently using our knowledge of the shape of the CLF. The other striking features of the proposed algorithm are that it provides procedures for initialization, has a fast converging iteration stage, and makes use of the information extracted in preliminary iterations to segment the data suitably to increase the effective signal-to-noise ratio (SNR). The computational complexity and the performance of the proposed algorithm are compared with other existing methods such as those based on linear prediction, KiSS/IQML, alternating projections (AP), and expectation-maximization (EM).

[1]  Ilan Ziskind,et al.  Maximum likelihood localization of multiple sources by alternating projection , 1988, IEEE Trans. Acoust. Speech Signal Process..

[2]  Louis L. Scharf,et al.  Signal processing applications of oblique projection operators , 1994, IEEE Trans. Signal Process..

[3]  D. W. Tufts,et al.  Reverberation Suppression and Modeling , 1993 .

[4]  L. Mcbride,et al.  A technique for the identification of linear systems , 1965 .

[5]  Donald W. Tufts,et al.  Determination of the dimension of a signal subspace from short data records , 1994, IEEE Trans. Signal Process..

[6]  S. Parthasarathy,et al.  Maximum-likelihood estimation of parameters of exponentially damped sinusoids , 1985, Proceedings of the IEEE.

[7]  Ehud Weinstein,et al.  Parameter estimation of superimposed signals using the EM algorithm , 1988, IEEE Trans. Acoust. Speech Signal Process..

[8]  S Umesh.  Fast maximum likelihood estimation of parameters in crowded signal environments , 1993 .

[9]  K. Steiglitz.  On the simultaneous estimation of poles and zeros in speech analysis , 1977 .

[10]  R. Kumaresan,et al.  Estimating the parameters of exponentially damped sinusoids and pole-zero modeling in noise , 1982 .

[11]  D. van Ormondt,et al.  Analysis of NMR Data Using Time Domain Fitting Procedures , 1992 .

[12]  V. U. Reddy,et al.  A Newton-based Ziskind-Wax alternating projection algorithm applied to spectral estimation , 1993, IEEE Trans. Signal Process..

[13]  Yoram Bresler,et al.  Exact maximum likelihood parameter estimation of superimposed exponential signals in noise , 1986, IEEE Trans. Acoust. Speech Signal Process..

[14]  Ramdas Kumaresan,et al.  An algorithm for pole-zero modeling and spectral analysis , 1986, IEEE Trans. Acoust. Speech Signal Process..

[15]  T. Henderson,et al.  Geometric methods for determining system poles from transient response , 1981 .

[16]  D. van Ormondt,et al.  Retrieval of frequencies, amplitudes, damping factors, and phases from time-domain signals using a linear least-squares procedure , 1985 .

[17]  J. Cadzow,et al.  Signal processing via least squares error modeling , 1990, IEEE ASSP Magazine.

[18]  Hongya Ge,et al.  Fast maximum likelihood estimation of signal parameters using the shape of the compressed likelihood function , 1993 .

[19]  Alfred O. Hero,et al.  Space-alternating generalized expectation-maximization algorithm , 1994, IEEE Trans. Signal Process..

[20]  Louis L. Scharf,et al.  On the complexity of IQML algorithms , 1992, IEEE Trans. Signal Process..

[21]  Louis L. Scharf,et al.  Matched subspace detectors , 1994, IEEE Trans. Signal Process..

[22]  Richard Travis Behrens.  Subspace signal processing structured noise , 1991 .

[23]  D. van Ormondt,et al.  Estimating the parameters of NMR signals by transforming to the frequency domain , 1990 .

[24]  Edward J. Wegman,et al.  Statistical Signal Processing , 1985 .

[25]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .