Rapid accurate frequency estimation of multiple resolved exponentials in noise

The estimation of the frequencies of the sum of multiple resolved exponentials in noise is an important problem due to its application in diverse areas from engineering to chemistry. Yet to date, no low cost Fourier-based algorithm has been successful at obtaining unbiased estimates that achieve the CramrRao lower bound (CRLB) over a wide range of signal-to-noise ratios. In this work, we achieve precisely this goal, proposing a fast yet accurate estimator that combines an iterative frequency-domain interpolation step with a leakage subtraction scheme. By analysing the asymptotic performance and the convergence behaviour of the estimator, we show that the estimate of each frequency converges to the asymptotic fixed point. Thus, the estimator is asymptotically unbiased and the variance is extremely close to the CRLB. We verify the theoretical analysis by extensive simulations, and demonstrate that the proposed algorithm is capable of obtaining more accurate estimates than state-of-the-art high resolution methods while requiring significantly less computational effort. HighlightsAn efficient frequency estimator for multiple resolved complex exponentials in noise is proposed.The proposed estimator is Fourier-based with no singular value decomposition or matrix inversion involved.The variance of the estimates obtained by the proposed method is extremely close to the CramrRao bound.Simulation results show that the proposed estimator can outperform state-of-the-art estimation approaches.

[1]  Luzhou Xu,et al.  IAA Spectral Estimation , 2011 .

[2]  Bernard Mulgrew,et al.  Iterative frequency estimation by interpolation on Fourier coefficients , 2005, IEEE Transactions on Signal Processing.

[3]  Jian Li,et al.  Source Localization and Sensing: A Nonparametric Iterative Adaptive Approach Based on Weighted Least Squares , 2010, IEEE Transactions on Aerospace and Electronic Systems.

[4]  Zhijiang Xie,et al.  Frequency estimation of the weighted real tones or resolved multiple tones by iterative interpolation DFT algorithm , 2015, Digit. Signal Process..

[5]  Jiufei Luo,et al.  Interpolated DFT algorithms with zero padding for classic windows , 2016 .

[6]  Elias Aboutanios,et al.  An algorithm for the parameter estimation of multiple superimposed exponentials in noise , 2015, 2015 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).

[7]  Sudhakar M. Pandit,et al.  Cramer-Rao lower bounds for a damped sinusoidal process , 1995, IEEE Trans. Signal Process..

[8]  K. Duda,et al.  Frequency and Damping Estimation Methods - An Overview , 2011 .

[9]  Thomas Kailath,et al.  ESPRIT-estimation of signal parameters via rotational invariance techniques , 1989, IEEE Trans. Acoust. Speech Signal Process..

[10]  Tapan K. Sarkar,et al.  Matrix pencil method for estimating parameters of exponentially damped/undamped sinusoids in noise , 1990, IEEE Trans. Acoust. Speech Signal Process..

[11]  S. Vanhuffel,et al.  Algorithm for time-domain NMR data fitting based on total least squares , 1994 .

[12]  Weize Sun,et al.  Efficient parameter estimation of multiple damped sinusoids by combining subspace and weighted least squares techniques , 2012, 2012 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).

[13]  J. Capon High-resolution frequency-wavenumber spectrum analysis , 1969 .

[14]  Krzysztof Duda,et al.  Interpolated DFT for $\sin^{\alpha}(x)$ Windows , 2014, IEEE Transactions on Instrumentation and Measurement.

[15]  George-Othon Glentis,et al.  Computationally Efficient Capon- and APES-Based Coherence Spectrum Estimation , 2012, IEEE Transactions on Signal Processing.

[16]  Elias Aboutanios,et al.  Estimating the parameters of sinusoids and decaying sinusoids in noise , 2011, IEEE Instrumentation & Measurement Magazine.

[17]  Tomasz P. Zielinski,et al.  DFT-based Estimation of Damped Oscillation Parameters in Low-Frequency Mechanical Spectroscopy , 2011, IEEE Transactions on Instrumentation and Measurement.

[18]  Gang Wei,et al.  A Noniterative Frequency Estimator With Rational Combination of Three Spectrum Lines , 2011, IEEE Transactions on Signal Processing.

[19]  Jian Li,et al.  Efficient mixed-spectrum estimation with applications to target feature extraction , 1995, Conference Record of The Twenty-Ninth Asilomar Conference on Signals, Systems and Computers.

[20]  Jian Li,et al.  IAA Spectral Estimation: Fast Implementation Using the Gohberg-Semencul Factorization , 2011, IEEE Trans. Signal Process..

[21]  J. Rissanen A UNIVERSAL PRIOR FOR INTEGERS AND ESTIMATION BY MINIMUM DESCRIPTION LENGTH , 1983 .

[22]  Jian Li,et al.  An adaptive filtering approach to spectral estimation and SAR imaging , 1996, IEEE Trans. Signal Process..

[23]  George-Othon Glentis,et al.  A Fast Algorithm for APES and Capon Spectral Estimation , 2008, IEEE Transactions on Signal Processing.

[24]  Barry G. Quinn,et al.  Estimating frequency by interpolation using Fourier coefficients , 1994, IEEE Trans. Signal Process..

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

[26]  Frankie K. W. Chan,et al.  Efficient Approach for Sinusoidal Frequency Estimation of Gapped Data , 2010, IEEE Signal Processing Letters.

[27]  Daniel Belega,et al.  Sine-wave parameter estimation by interpolated DFT method based on new cosine windows with high interference rejection capability , 2014, Digit. Signal Process..

[28]  R. O. Schmidt,et al.  Multiple emitter location and signal Parameter estimation , 1986 .

[29]  Barry G. Quinn,et al.  The Estimation and Tracking of Frequency , 2001 .

[30]  Petre Stoica,et al.  Introduction to spectral analysis , 1997 .

[31]  Jesus Selva An efficient Newton-type method for the computation of ML estimators in a uniform linear array , 2005, IEEE Transactions on Signal Processing.

[32]  H. Akaike A new look at the statistical model identification , 1974 .

[33]  J. Selva Efficient maximum likelihood estimation of a 2-D complex sinusoidal based on barycentric interpolation , 2011, 2011 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).

[34]  E. J. Hannan,et al.  The maximum of the periodogram , 1983 .

[35]  Elias Aboutanios,et al.  Estimating Parameters of Multiple Damped Complex Sinusoids with Model Order Estimation , 2016, 2016 IEEE International Workshop on Signal Processing Systems (SiPS).

[36]  Peter T. Gough,et al.  A fast spectral estimation algorithm based on the FFT , 1994, IEEE Trans. Signal Process..

[37]  Qingfeng Meng,et al.  An Interpolation Algorithm for Discrete Fourier Transforms of Weighted Damped Sinusoidal Signals , 2014, IEEE Transactions on Instrumentation and Measurement.

[38]  Srinivasan Umesh,et al.  Estimation of parameters of exponentially damped sinusoids using fast maximum likelihood estimation with application to NMR spectroscopy data , 1996, IEEE Trans. Signal Process..

[39]  Josef A. Nossek,et al.  Unitary ESPRIT: how to obtain increased estimation accuracy with a reduced computational burden , 1995, IEEE Trans. Signal Process..