Title: A NEW ALGORITHM FOR THE SVD OF A LONG PRODUCT OF MATRICES AND THE STABILITY OF PRODUCTS
Abstract: Lyapunov exponents can be estimated by accurately computing the singular values of long products of matrices, with perhaps 1000 or more factor matrices. These products have extremely large ratios between the largest and smallest eigenvalues. A variant of Rutishauser’s Cholesky LR algorithm for computing eigenvalues of symmetric matrices is used to obtain a new algorithm for computing the singular values and vectors of long products of matrices with small backward error in the factor matrices. The basic product SVD algorithm can also be accelerated using hyperbolic Givens’ rotations. The method is competitive with Jacobi-based methods for certain problems as numerical results indicate. Some properties of the product SVD factorization are also discussed, including uniqueness and stability. The concept of a stable product is introduced; for such products, all singular values can be computed to high relative accuracy.
Publication Year: 1997
Publication Date: 1997-01-01
Language: en
Type: article
Access and Citation
Cited By Count: 15
AI Researcher Chatbot
Get quick answers to your questions about the article from our AI researcher chatbot