Minor and major subspace computation of large matrices

Large matrices arise in many formulations in signal processing and control. In this paper, a Rayleigh quotient iteration (RQI) method for locating the minimum eigenpair for symmetric positive definite matrix pencil has been developed. This method has a cubic convergence rate and does not require computation of matrix inversion. The core procedure is based on a modified Rayleigh quotient iteration (MRQI) which uses a line search (exact or approximate) to determine a vector of steepest descent. As a special case, the proposed algorithm is customized to solve high resolution temporal and spatial frequency tracking problems. The eigenstructure tracking algorithm has update complexity O(n/sup 2/p), where n is the data dimension and p is the dimension of the minor or major subspaces. The performance of these algorithms is tested with several examples.

[1]  F. Trench,et al.  Numerical solution of the eigenvalue problem for Hermitian Toeplitz matrices , 1989 .

[2]  Sabine Van Huffel,et al.  Total least squares problem - computational aspects and analysis , 1991, Frontiers in applied mathematics.

[3]  Friedrich L. Bauer,et al.  On certain methods for expanding the characteristic polynomial , 1959, Numerische Mathematik.

[4]  A. Ostrowski.  On the convergence of the Rayleigh quotient iteration for the computation of the characteristic roots and vectors. I , 1957 .

[5]  Fazal Noor,et al.  Recursive and iterative algorithms for computing eigenvalues of Hermitian Toeplitz matrices , 1993, IEEE Trans. Signal Process..

[6]  K. Bathe,et al.  Large Eigenvalue Problems in Dynamic Analysis , 1972 .

[7]  J. H. Wilkinson.  The algebraic eigenvalue problem , 1966 .

[8]  B. Parlett.  The Rayleigh Quotient Iteration and Some Generalizations for Nonnormal Matrices , 1974 .

[9]  G. Stewart.  Error and Perturbation Bounds for Subspaces Associated with Certain Eigenvalue Problems , 1973 .

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

[11]  David S. Watkins,et al.  Understanding the $QR$ Algorithm , 1982 .

[12]  Ricardo D. Fierro,et al.  The Total Least Squares Problem: Computational Aspects and Analysis (S. Van Huffel and J. Vandewalle) , 1993, SIAM Rev..

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

[14]  George Cybenko,et al.  Computing thr minimum eigenvalue of a symmetric positive definite Toeplitz matrix , 1984 .

[15]  P. C. Chowdhury.  The truncated Lanczos algorithm for partial solution of the symmetric eigenproblem , 1976 .