Refine
Language
- English (2)
Keywords
- optimal control (2)
- BVD pencil (1)
- Cayley transformation (1)
- Hamiltonian matrix (1)
- Lagrangian subspace (1)
- disk function method (1)
- even pencil (1)
- graph subspace (1)
- matrix pencil (1)
- pencil arithmetic (1)
Project
- C29 (2) (remove)
Application Area
- C (2)
We derive a new representation of Lagrangian subspaces in the form
%
\[
{\mathrm Im}\Pi^T [I,X]^T,
\]
%
where $\Pi$ is a symplectic matrix which is the product of a permutation matrix and a real orthogonal diagonal matrix, and $X$ satisfies
%
\[
\abs{X_{ij}} \leq \begin{cases}1 & \text{if $i=j$,}\\ \sqrt{2} & \text{if $i\neq j$.} \end{cases}
\]
%
This representation allows to limit element growth in the context of doubling algorithms
for the computation of Lagrangian subspaces and the solution of Riccati equations.
It is shown that a simple doubling algorithm using this representation can reach full machine accuracy on a wide range of problems, obtaining invariant subspaces of the same quality as those computed by the state-of-the-art algorithms based on orthogonal transformations.
The same idea carries over to representations of arbitrary subspaces and can be used
for other types of structured pencils.
We propose a generalization of the Structured Doubling Algorithm (SDA) to compute invariant subspaces
of structured matrix pencils
that arise in the context of solving linear quadratic optimal control problems.
The new algorithm is
designed to attain better accuracy when the classical Riccati equation approach for the solution of the optimal control problem is not well suited because
the stable and unstable invariant subspaces are not well separated (due to eigenvalues near or on the imaginary
axis) or in the case when the Riccati solution does not exist at all. We analyze the convergence
of the method and compare the new method with the classical SDA algorithm as well as some recent structured QR-methods.