Computation of Stable Invariant Subspaces of Hamiltonian Matrices

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.

[1]  G. W. Stewart,et al.  Algorithm 506: HQR3 and EXCHNG: Fortran Subroutines for Calculating and Ordering the Eigenvalues of a Real Upper Hessenberg Matrix [F2] , 1976, TOMS.

[2]  C. Loan,et al.  A Schur decomposition for Hamiltonian matrices , 1981 .

[3]  A. Laub.  A schur method for solving algebraic Riccati equations , 1978, 1978 IEEE Conference on Decision and Control including the 17th Symposium on Adaptive Processes.

[4]  A. Bunse-Gerstner,et al.  A symplectic QR like algorithm for the solution of the real algebraic Riccati equation , 1986 .

[5]  Ralph Byers,et al.  Hamiltonian and symplectic algorithms for the algebraic riccati equation , 1983 .

[6]  Edward J. Davison,et al.  Further discussion on the calculation of transmission zeros , 1978, Autom..

[7]  A. MacFarlane.  An Eigenvector Solution of the Optimal Linear Regulator Problem , 1963 .

[8]  V. Mehrmann,et al.  On Hamiltonian and symplectic Hessenberg forms , 1991 .

[9]  James Hardy Wilkinson,et al.  Rounding errors in algebraic processes , 1964, IFIP Congress.

[10]  A. Bunse-Gerstner,et al.  Numerical Methods for Algebraic Riccati Equations , 1989 .

[11]  Rajnikant V. Patel,et al.  Multivariable System Theory and Design , 1981 .

[12]  Jack J. Dongarra,et al.  Numerical Considerations in Computing Invariant Subspaces , 1992, SIAM J. Matrix Anal. Appl..

[13]  C. Loan.  A Symplectic Method for Approximating All the Eigenvalues of a Hamiltonian Matrix , 1982 .

[14]  A. Laub.  Invariant Subspace Methods for the Numerical Solution of Riccati Equations , 1991 .

[15]  R. Byers.  A Hamiltonian $QR$ Algorithm , 1986 .