login

A Generalized Eigenvalue Approach for Solving Riccati Equations

SIAM Journal on Scientific and Statistical ComputingPublished 1 June 1981
Paul Van Dooren
Citations380

TL;DR

A numerically stable algorithm is derived to compute orthonormal bases for any deflating subspace of a regular pencil $\lambda B - A$ based on an update of the QZ-algorithm, in order to obtain any desired ordering of eigenvalues in the quasitriangular forms constructed by this algorithm.

Abstract

A numerically stable algorithm is derived to compute orthonormal bases for any deflating subspace of a regular pencil $\lambda B - A$. The method is based on an update of the $QZ$-algorithm, in order to obtain any desired ordering of eigenvalues in the quasitriangular forms constructed by this algorithm. As applications we discuss a new approach to solve Riccati equations arising in linear system theory. The computation of deflating subspaces with specified spectrum is shown to be of crucial importance here.

Keywords

Computer ScienceMathematicsPhysics and Astronomy