Abstract
This paper addresses some numerical issues that arise in computing a basis for the stable invariant subspace of a Hamiltonian matrix. Such a basis is required in solving the algebraic Riccati equation using the well-known method due to Laub. Two algorithms based on certain properties of Hamiltonian matrices are proposed as viable alternatives to the conventional approach.
Original language | American English |
---|---|
Journal | SIAM Journal on Matrix Analysis and Applications |
Volume | 15 |
DOIs | |
State | Published - Jan 1 1994 |
Keywords
- HAMILTONIAN MATRICES; EIGENVALUES; INVARIANT SUBSPACES; ALGEBRAIC RICCATI EQUATION
Disciplines
- Electrical and Computer Engineering
- Engineering