3.2.2 The InverseFree Generalized Eigenvector and Schur Methods for the Riccati Equations
The Schur, Newton, and matrix Signfunction methods (described in the following sections) for the CARE and DARE, all require the computation of . Furthermore, the Schur and Newton methods for the DARE require explicit computation of , which may lead to inaccurate solutions when these matrices are singular or nearly singular. The problems can be overcome by using the socalled inversefree generalized eigenvector and Schur methods, developed by Pappas, Laub, and Sandell (1980). The inversefree methods replace the Hamiltonian eigenproblem associated with the CARE by the generalized eigenvalue problem given by the following pencil with order
The pencil is further compressed into a pencil using the QR factorization of , without affecting the deflating subspaces (Van Dooren (1981)). The CARE is then solved by construction of either a basis of the generalized eigenvectors associated with the stable eigenvalues (the inversefree generalized eigenvector method) or an orthonormal basis of the stable deflating subspace (the inversefree generalized Schur method) of the compressed pencil.
For the DARE , the pencil considered is
This pencil is compressed into a pencil using the QR decomposition of . For details, see Datta (2003).
Inversefree methods for the Riccati equations.
When the associated matrix pencil has multiple or nearly multiple eigenvalues, the inversefree generalized Schur methods may perform better than the inversefree generalized eigenvector methods. In fact, failure of the generalized eigenvector method sometimes occurs because the pencil does not have enough eigenvectors to form a basis (illustrated in the following example).
The generalized eigenvector methods use the function GeneralizedEigensystem and the generalized Schur methods use the function GeneralizedSchurDecomposition, both in Chapter 9, to solve the appropriate generalized eigenvalue problems (see Datta (2003) for details).
Consider solving the DARE with matrices and of the following statespace system of a process control of a paper machine.
In[10]:=
Out[10]=
The weighting matrices are chosen to be identity matrices of the appropriate dimensions.
In[11]:=
This solves the DARE by the inversefree generalized Schur method.
In[12]:=
Out[12]=
This auxiliary function computes the residual norm of the solution.
In[13]:=
In[14]:=
Out[14]=
