Title :
On the matrix sign function method for the computation of invariant subspaces
Author :
Byers, R. ; He, Chunyang ; Mehrmann, Volker
Author_Institution :
Dept. of Math., Kansas Univ., Lawrence, KS, USA
Abstract :
There is some concern about the numerical stability of algorithms that use the matrix sign function to solve Riccati and Lyapunov equations and to find bases of invariant subspaces. In this paper, we demonstrate that evaluating the matrix sign function is a more ill-conditioned computational problem than the problem of finding bases of the two invariant subspaces. Nevertheless, we also give perturbation and error analyses, which show that the accuracy of the Newton iteration with the scaling for the computation of the invariant subspaces in most circumstances is competitive with conventional methods
Keywords :
Lyapunov matrix equations; Newton method; Riccati equations; convergence of numerical methods; Lyapunov equations; Newton iteration; Riccati equations; error analysis; ill-conditioned computational problem; invariant subspace computation; matrix sign function method; numerical stability; perturbation analysis; Computer architecture; Contracts; Eigenvalues and eigenfunctions; Error analysis; Helium; Mathematics; Matrix decomposition; Numerical stability; Riccati equations; Snow;
Conference_Titel :
Computer-Aided Control System Design, 1996., Proceedings of the 1996 IEEE International Symposium on
Conference_Location :
Dearborn, MI
Print_ISBN :
0-7803-3032-3
DOI :
10.1109/CACSD.1996.555200
