Performance of ESPRIT for Estimating Mixtures of Complex Exponentials Modulated by Polynomials

High-resolution (HR) methods are known to provide accurate frequency estimates for discrete spectra. The polynomial amplitude complex exponentials (PACE) model, also called quasi-polynomial model in the literature, was presented as the most general model tractable by HR methods. A subspace-based estimation scheme was recently proposed, derived from the classical ESPRIT algorithm. In this paper, we focus on the performance of this estimator. We first present some asymptotic expansions of the estimated parameters, obtained at the first order under the assumption of a high signal-to-noise ratio (SNR). Then the performance of the generalized ESPRIT algorithm for estimating the parameters of this model is analyzed in terms of bias and variance, and compared to the Cramer-Rao bounds (CRB). This performance is studied in an asymptotic context, and it is proved that the efficiency of undamped single poles estimators is close to the optimality. Moreover, our results show that the best performance is obtained for a proper dimensioning of the data. To illustrate the practical capabilities of the generalized ESPRIT algorithm, we finally propose an application to ARMA filter synthesis, in the context of system conversion from continuous time to discrete time.

[1]  Roland Badeau,et al.  Méthodes à haute résolution pour l'estimation et le suivi de sinusoïdes modulées. Application aux signaux de musique , 2005 .

[2]  Sarangarajan Parthasarathy,et al.  The statistical performance of state-variable balancing and Prony's method in parameter estimation , 1987, ICASSP '87. IEEE International Conference on Acoustics, Speech, and Signal Processing.

[3]  Ronald L. Graham,et al.  Concrete mathematics - a foundation for computer science , 1991 .

[4]  Thomas Kailath,et al.  Linear Systems , 1980 .

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

[6]  Steven Roman The Umbral Calculus , 1984 .

[7]  Petre Stoica,et al.  Second-order properties of MUSIC and ESPRIT estimates of sinusoidal frequencies in high SNR scenarios , 1993 .

[8]  Zhonggang Zeng,et al.  On Ill-conditioned Eigenvalues, Multiple Roots of Polynomials, and Their Accurate Computations , 1998 .

[9]  Tapan K. Sarkar,et al.  Perturbation analysis of TK method for harmonic retrieval problems , 1988, IEEE Trans. Acoust. Speech Signal Process..

[10]  Roland Badeau,et al.  Yet another subspace tracker , 2005, Proceedings. (ICASSP '05). IEEE International Conference on Acoustics, Speech, and Signal Processing, 2005..

[11]  H. Zeiger,et al.  Approximate linear realizations of given dimension via Ho's algorithm , 1974 .

[12]  A. Swindlehurst,et al.  Subspace-based signal analysis using singular value decomposition , 1993, Proc. IEEE.

[13]  A.M. Sabatini,et al.  Correlation receivers using Laguerre filter banks for modelling narrowband ultrasonic echoes and estimating their time-of-flights , 1997, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[14]  Benjamin Friedlander,et al.  On the accuracy of the Kumaresan-Tufts method for estimating complex damped exponentials , 1987, IEEE Trans. Acoust. Speech Signal Process..

[15]  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..

[16]  Cramér-Rao bound for the estimates of frequencies and damping factors of quasipolynomials in noise , 1992 .

[17]  Petre Stoica,et al.  Study of the statistical performance of the Pisarenko harmonic decomposition method , 1988 .

[18]  Charles R. Johnson,et al.  Matrix analysis , 1985, Statistical Inference for Engineers and Data Scientists.

[19]  Roland Badeau,et al.  High-resolution spectral analysis of mixtures of complex exponentials modulated by polynomials , 2006, IEEE Transactions on Signal Processing.

[20]  Thomas Kailath,et al.  ESPRIT-A subspace rotation approach to estimation of parameters of cisoids in noise , 1986, IEEE Trans. Acoust. Speech Signal Process..

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

[22]  Roland Badeau,et al.  A new perturbation analysis for signal enumeration in rotational invariance techniques , 2006, IEEE Transactions on Signal Processing.

[23]  Roland Badeau,et al.  Sliding window adaptive SVD algorithms , 2004, IEEE Transactions on Signal Processing.

[24]  Roland Badeau,et al.  CramÉr–Rao Bounds for Multiple Poles and Coefficients of Quasi-Polynomials in Colored Noise , 2008, IEEE Transactions on Signal Processing.

[25]  J. L. Guiñón,et al.  The spectrum of atomic hydrogen , 1989 .

[26]  Torsten Söderström,et al.  Statistical analysis of MUSIC and subspace rotation estimates of sinusoidal frequencies , 1991, IEEE Trans. Signal Process..

[27]  Tapan K. Sarkar,et al.  On SVD for estimating generalized eigenvalues of singular matrix pencil in noise , 1991, IEEE Trans. Signal Process..

[28]  Bin Yang,et al.  Projection approximation subspace tracking , 1995, IEEE Trans. Signal Process..

[29]  Philippe Loubaton,et al.  A subspace algorithm for certain blind identification problems , 1997, IEEE Trans. Inf. Theory.

[30]  Thomas Kailath,et al.  Detection of signals by information theoretic criteria , 1985, IEEE Trans. Acoust. Speech Signal Process..

[31]  V. Pisarenko The Retrieval of Harmonics from a Covariance Function , 1973 .

[32]  Roland Badeau,et al.  Fast approximated power iteration subspace tracking , 2005, IEEE Transactions on Signal Processing.

[33]  Laura Marcu,et al.  A Nonparametric Method for Analysis of Fluorescence Emission in Combined Time and Wavelength Dimensions , 2005, Annals of Biomedical Engineering.

[34]  S. Kay Fundamentals of statistical signal processing: estimation theory , 1993 .

[35]  A. Medvedev,et al.  L/sup 2/ time delay estimation by means of Laguerre functions , 1999, Proceedings of the 1999 American Control Conference (Cat. No. 99CH36251).

[36]  D. Filipović EXPONENTIAL-POLYNOMIAL FAMILIES AND THE TERM STRUCTURE OF INTEREST RATES , 2000 .