Simple, effective computation of principal eigenvectors and their eigenvalues and application to high-resolution estimation of frequencies

We present the results of an investigation of the Prony-Lanczos (P-L) method [14], [38] and the power method [39] for simple computation of approximations to a few eigenvectors and eigenvalues of a Hermitian matrix. We are motivated by realization of high-resolution signal processing in an integrated circuit. The computational speeds of the above methods are analyzed. They are completely dependent on the speed of a matrix-vector product operation. If only a few eigenvalues or eigenvectors are needed, the suggested methods can substitute for methods of the LINPACK or EISPACK subroutine libraries. The accuracies of the suggested methods are evaluated using matrices formed from simulated data consisting of two sinusoids plus Gaussian noise. Comparisons are made to the corresponding eigenvalues and eigenvectors obtained using LINPACK as a reference. Also, the accuracies of frequency estimates obtained from the eigenvectors are compared.

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

[2]  Kenji Nakayama.  Permuted difference coefficient digital filters , 1981, ICASSP.

[3]  C. Lanczos.  Applied Analysis , 1961 .

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

[5]  R. Kumaresan,et al.  Data-Adaptive Detection of a Weak Signal , 1983, IEEE Transactions on Aerospace and Electronic Systems.

[6]  R. Kumaresan,et al.  Data adaptive signal estimation by singular value decomposition of a data matrix , 1982, Proceedings of the IEEE.

[7]  Georges Bienvenu,et al.  Adaptivity to background noise spatial coherence for high resolution passive methods , 1980, ICASSP.

[8]  R. Frisch.  Correlation and scatter in statistical variables , 1951 .

[9]  L. Scharf,et al.  A Prony method for noisy data: Choosing the signal components and selecting the order in exponential signal models , 1984, Proceedings of the IEEE.

[10]  S. Lang,et al.  Frequency estimation with maximum entropy spectral estimators , 1980 .

[11]  R. Kumaresan,et al.  Improved spectral resolution III: Efficient realization , 1980, Proceedings of the IEEE.

[12]  F. B. Hildebrand,et al.  Introduction To Numerical Analysis , 1957 .

[13]  S.M. Kay,et al.  Spectrum analysis—A modern perspective , 1981, Proceedings of the IEEE.

[14]  Karl Pearson F.R.S..  LIII. On lines and planes of closest fit to systems of points in space , 1901 .

[15]  G. Boudreaux-Bartels,et al.  Discrete Fourier transform using summation by parts , 1987, ICASSP '87. IEEE International Conference on Acoustics, Speech, and Signal Processing.

[16]  I. Kirsteins,et al.  On the probability density of signal-to-noise ratio in an improved adaptive detector , 1984, ICASSP '85. IEEE International Conference on Acoustics, Speech, and Signal Processing.

[17]  Donald W. Tufts,et al.  Simple, effective computation of principal eigenvectors and their eigenvalues and application to high-resolution estimation of frequencies , 1985, ICASSP '85. IEEE International Conference on Acoustics, Speech, and Signal Processing.

[18]  Ramdas Kumaresan,et al.  Improved spectral resolution II , 1980, ICASSP.

[19]  S.S. Reddi,et al.  Multiple Source Location-A Digital Approach , 1979, IEEE Transactions on Aerospace and Electronic Systems.

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

[21]  Gene H. Golub,et al.  Matrix computations , 1983 .

[22]  Georges Bienvenu,et al.  ADAPTIVE HIGH RESOLUTION SPATIAL DISCRIMINATION OF PASSIVE SOURCES , 1981 .

[23]  Gloria Faye Boudreaux-Bartels,et al.  Time-frequency signal processing algorithms: analysis and synthesis using wigner distributions , 1984 .

[24]  Albert H Nuttall,et al.  Spectral Analysis of a Univariate Process with Bad Data Points, via Maximum Entropy and Linear Predictive Techniques , 1976 .

[25]  R. Kumaresan,et al.  Estimation of frequencies of multiple sinusoids: Making linear prediction perform like maximum likelihood , 1982, Proceedings of the IEEE.

[26]  Chandler Davis.  The rotation of eigenvectors by a perturbation , 1963 .

[27]  T. Ulrych,et al.  Time series modeling and maximum entropy , 1976 .

[28]  J. E. Hudson.  Adaptive Array Principles , 1981 .

[29]  T. Claasen,et al.  THE WIGNER DISTRIBUTION - A TOOL FOR TIME-FREQUENCY SIGNAL ANALYSIS , 1980 .

[30]  W. Kahan,et al.  The Rotation of Eigenvectors by a Perturbation. III , 1970 .

[31]  R. Bracewell.  The fast Hartley transform , 1984, Proceedings of the IEEE.

[32]  A. Cantoni,et al.  Resolving the directions of sources in a correlated field incident on an array , 1980 .

[33]  Calyampudi R. Rao.  The use and interpretation of principal component analysis in applied research , 1964 .

[34]  C. Eckart,et al.  The approximation of one matrix by another of lower rank , 1936 .

[35]  R. R. Hocking,et al.  Selection of the Best Subset in Regression Analysis , 1967 .

[36]  K. Nakayama.  Permuted difference coefficient realization of FIR digital filters , 1982 .

[37]  Ramdas Kumaresan,et al.  Accurate parameter estimation of noisy speech-like signals , 1982, ICASSP.