Numerical Solution of Algebraic Matrix Riccati Equations

摘要

Numerical issues related to the computational solution of the algebraic matrix Riccati equation are studied. The approach uses the generalized eigenproblem formulation for the solution of general forms of applications. These general forms result from control and filtering problems for systems in generalized (or implicit or descriptor) state space form. A Newton-type iterative refinement procedure for the generalized Riccati solution is derived. The issue of numerical condition of the Riccati problem is addressed. Balancing to improve numerical condition is discussed. An overview of a software package coded in FORTRAN is given. Results of numerical experiments are reported.