Estimating Frequency, Amplitude and Phase of Two Sinusoids with Very Close Frequencies

This paper presents an algorithm to estimate the parameters of two closely spaced sinusoids, providing a frequency resolution that is more than 800 times greater than that obtained by using the Discrete Fourier Transform (DFT). The strategy uses a highly optimized grid search approach to accurately estimate frequency, amplitude and phase of both sinusoids, keeping at the same time the computational effort at reasonable levels. The proposed method has three main characteristics: 1) a high frequency resolution; 2) frequency, amplitude and phase are all estimated at once using one single package; 3) it does not rely on any statistical assumption or constraint. Potential applications to this strategy include the difficult task of resolving coincident partials of instruments in musical signals.

• Jayme G. A. Barbedo
• Amauri Lopes

• Closely spaced sinusoids
• high-resolution parameter estimation
• optimized grid search.

[1] L. Benaroya, F. Bimbot, R. Gribonval, "Audio source separation with a
single sensor," IEEE Tran. on Audio, Speech, and Language Proc., vol.
14, no. 1, pp. 191-199, Jan. 2006.
[2] S. Kay and S. L. Marple, "Spectrum analysis- A modern perspective,"
Proc. IEEE, vol. 9, pp. 1380-1419, Nov. 1981.
[3] Proceedings of IEEE, Special Issue on Spectral Estimation, vol. 70,
Sept. 1982.
[4] S. M. Kay, Modern Spectral Estimation: Theory and Application,
Prentice-Hall Signal Processing Series, Englewood Cliffs, NJ, USA,
1988.
[5] S. L. Marple, Digital Spectral Analysis with Applications, Prentice-Hall
Signal Processing Series, Englewood Cliffs, NJ, USA, 1987.
[6] P. Stoica, R. Moses, Introduction to Spectral Analysis, Prentice-Hall,
Upper Saddle River, NJ, USA, 1997.
[7] K. Kitsios, A. Spanias, B. Welfert, P. Loizou, "An Adaptive Modified
Covariance Algorithm For Spectral Analysis," in Proc. of 8th IEEE
Signal Processing Workshop on Statistical Signal and Array Processing,
pp. 56-59, Corfu, Greece, Jun 1996.
[8] [8] Y. Bresler and A. Macovski, "Exact maximum likelihood parameter
estimation of superimposed exponential signals in noise," IEEE Trans.
Acoust., Speech and Signal Processing, vol. ASSP-34, n. 5, pp. 1081-
1089, Oct. 1986.
[9] P. Stoica, R. L. Moses, B. Friedlander and T. Söderström, "Maximum
Likelihood estimation of the parameters of multiple sinusoids from noisy
measurements," IEEE Trans. Acoust., Speech and Signal Processing,
vol. 37, n. 3, pp. 378-391, March 1989.
[10] P. Stoica and C. Sharman "Novel eigenanalysis method for direction
estimation," IEE Proceedings, vol. 137, No. 1, pp. 19-26, Feb. 1990.
[11] D. C. Rife, R. R. Boorstyn, "Multiple tone parameter estimation from
discrete-time observation," Bell System Technical Journal, vol. 55, no.
9, pp. 1389 - 1410, Nov. 1976.
[12] P. Stoica, P. H├ñndel, T. Söderström, "Approximate maximum likelihood
frequency estimation," Automatica, vol. 30, no. 1, pp. 131 - 145, January
1994.
[13] M. D. Macleod, "Joint detection and high resolution ML estimation of
multiple sinusoids in noise," in Proc. of 2001 IEEE International
Conference on Acoustics, Speech, and Signal Processing (ICASSP), Vol.
5, pp. 3065 - 3068, Salt Lake City, May 2001.
[14] J.-K. Hwang, Y.-C. Chen, "Superresolution Frequency Estimation by
Alternating Notch Periodogram," IEEE Transaction on Signal
Processing, vol. 41, no. 2, pp. 727 - 741, February 1993.
[15] P. Stoica, K.C. Sharman, "Maximum Likelhood Methods for Directionof-
Arrival Estimation", IEEE Trans. Acoust., Speech and Signal
Processing, vol. 38, N. 7 pp. 1132-1143, July 1990.
[16] J. P. Burg, "Maximum Entropy Spectral Analysis," in Proceedings of the
37th Annual International Meeting of the Society Exploration
Geophysicists, Oklahoma City, USA, 1967.
[17] P. Depalle, T. Helie, "Extraction of spectral peak parameters using a
short-time Fourier transform modeling and no sidelobe windows," in
Proc.of the IEEE Workshop on Applications of Signal Processing to
Audio and Acoustics, pp. 19 - 22, New Paltz, USA, October 1997.
[18] D. Tufts and R. Kumaresan, "Estimation of Frequencies of Multiple
Sinusoids: Making Linear Prediction Perform Like Maximum
Likelihood", Proceedings of the IEEE, vol. 70, No. 9, pp. 975-989, Sept.
1982.
[19] R. O. Schmidt, "A signal subspace approach to multiple emitter location
and spectral estimation," Ph.D. dissertation, Stanford University,
Stanford, CA, 1981.
[20] H. Wang and M. Kaveh, "On the Performance of Signal-Subspace
Processing-Part I: Narrow-Band Systems", IEEE Trans. Acoust., Speech
and Signal Processing, vol. 34, N. 5 pp. 1201- 1209, Oct. 1986.
[21] J. A. Cadzow, "Multiple Source Location-The Signal Subspace
Approach", IEEE Trans. Acoust., Speech and Signal Processing, vol. 38,
N. 7 pp. 1110-1125, July 1990.
[22] R. Roy, T. Kailath, "ESPRIT-Estimation of Signal Parameters Via
Rotational Invariance Techniques," IEEE Trans. on Acoustics, Speech
and Signal Processing, vol. 37, no. 7, July 1989.
[23] M. Viberg and B. Ottersten, "Sensor array processing based on subspace
fitting", IEEE Trans. on Signal Processing, vol. 39, No. 5, pp. 1110-
1121, May. 1991.
[24] M. Viberg, B. Ottersten and T. Kailath, "Detection and estimation in
sensor arrays using weighted subspace fitting", IEEE Trans. on Signal
Processing, vol. 39, No. 11, pp. 2436-2449, Nov. 1991.
[25] D. M. Wilkes, J. A. Cadzow, "The effects of phase on high-resolution
frequency estimators," IEEE Transaction on Signal Processing, vol. 41,
no. 3, pp. 1319 - 1330, March 1993.
[26] G. Feyh, "Inverse Eigenvalue Problem for Sinusoidal Frequency
Estimation," in Proc. of 1997 IEEE International Conference on
Acoustics, Speech, and Signal Processing (ICASSP), vol. 5, pp. 3985 -
3988, Munich, Germany, April 1997.
[27] F.-X. Ge, H. Zhang, J. Yang, Y.-N. Peng, "Super-resolution frequency
estimation of the sinusoidal signals with unknown lowpass envelopes,"
in Proc. of 2003 IEEE Radar Conference, pp. 273 - 273, Huntsville,
USA, May 2003.
[28] M. G. Christensen, A. Jakobsson, S. H. Jensen, "Joint High-Resolution
Fundamental Frequency and Order Estimation," IEEE Transactions on
Audio, Speech, and Language Processing, Vol. 15, No. 5, pp. 1635 -
1644, July 2007.
[29] J. Capon, "High-resolution frequency-wavenumber spectrum analysis,"
Proc. IEEE, Vol. 57, pp. 1408 - 1418, August 1969.
[30] S. Nishimura, "Adaptive detection and enhancement of closely spaced
sinusoids using multirate techniques," in Proc. of 1993 IEEE
International Conference on Acoustics, Speech, and Signal Processing
(ICASSP), vol. 3, pp. 400 - 403, Minneapolis, USA, April 1993.
[31] C. I. Byrnes, T. T. Georgiou, A. Lindquist, "A New Approach to
Spectral Estimation: A Tunable High-Resolution Spectral Estimator,"
IEEE Transactions on Signal Processing, Vol. 48, no. 11, pp. 3189 -
3205, November 2000.
[32] C. B. Lambrecht, M. Karrakchou, "Wavelet packets-based highresolution
spectral estimation," Signal Processing, vol. 47, no. 2, pp. 135
- 144, November 1995.
[33] H.-T. Li, P. M. Djuric, "A novel approach to detection of closely spaced
sinusoids," Signal Processing, Vol. 51, No. 2, pp. 93 - 104, June 1996.
[34] P. O-Shea, "A High-Resolution Spectral Analysis Algorithm for Power-
System Disturbance Monitoring," IEEE Transactions on Power Systems,
Vol. 17, no. 3, August 2002.
[35] S. Hainsworth, M. Macleod, "On Sinusoidal Parameter Estimation," in
Proc. of the 6th International Conference on Digital Audio Effects,
London, UK, September 2003.
[36] I. Sarkar, A. T. Fam, "The interlaced chirp Z transform," Signal
Processing, vol. 86, no. 9, pp. 2221 - 2232, September 2006.
[37] O. Besson, P. Stoica, "Nonlinear least-squares approach to frequency
estimation and detection for sinusoidal signals with arbitrary envelope,"
Digital Signal Processing, vol. 9, no. 1, pp. 45-56, January 1999.
[38] T. Tolonen "Methods for Separation of Harmonic Sound. Sources using
Sinusoidal Modeling," presented at the 106th AES Convention, paper
4958, Munich, Germany, 1999.
[39] R. H. Roy, "ESPRIT-Estimation of signal parameters via rotational
invariance techniques," Ph.D. dissertation, Stanford University,
Stanford, CA, 1987.
[40] P. Stoica, T. Soderstrom, "Statistical analysis of MUSIC and ESPRIT
estimates of sinusoidal frequencies," in Proc. International Conference
on Acoustics, Speech, and Signal Processing, vol. 5, pp. 3273 - 3276,
Toronto, Canada, April 1991.
[41] O. Besson, P. Stoica, "Analysis of MUSIC and ESPRIT frequency
estimates for sinusoidal signals with lowpass envelopes," IEEE
Transactions on Signal Processing, vo1.44, no. 9, pp. 2359-2364,
September 1996.
[42] P. Stoica, H. Li, J. Li, "Amplitude estimation of sinusoidal signals:
survey, new results, and an application," IEEE Transactions on Signal
Processing, Vol. 48, No. 2, pp. 338 - 352, February 2000.
[43] S. M. Kay, Fundamentals of Statistical Signal Processing - Estimation
Theory, Prentice-Hall Signal Processing Series, Englewood Cliffs, NJ,
USA, 1993.