Numerical computation of deflating subspaces of Skew-Hamiltonian/Hamiltonian pencils

We discuss the numerical solution of structured generalized eigenvalue problems that arise from linear-quadratic optimal control problems, $H_{\infty}$ optimization, multibody systems, and many other areas of applied mathematics, physics, and chemistry. The classical approach for these problems requ...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:SIAM journal on matrix analysis and applications 2002-01, Vol.24 (1), p.165-190
Hauptverfasser: BENNER, Peter, BYERS, Ralph, MEHRMANN, Volker, HONGGUO XU
Format: Artikel
Sprache:eng
Schlagworte:
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Beschreibung
Zusammenfassung:We discuss the numerical solution of structured generalized eigenvalue problems that arise from linear-quadratic optimal control problems, $H_{\infty}$ optimization, multibody systems, and many other areas of applied mathematics, physics, and chemistry. The classical approach for these problems requires computing invariant and deflating subspaces of matrices and matrix pencils with Hamiltonian and/or skew-Hamiltonian structure. We extend the recently developed methods for Hamiltonian matrices to the general case of skew-Hamiltonian/Hamiltonian pencils. The algorithms circumvent problems with skew-Hamiltonian/Hamiltonian matrix pencils that lack structured Schur forms by embedding them into matrix pencils that always admit a structured Schur form. The rounding error analysis of the resulting algorithms is favorable. For the embedded matrix pencils, the algorithms use structure-preserving unitary matrix computations and are strongly backwards stable, i.e., they compute the exact structured Schur form of a nearby matrix pencil with the same structure.
ISSN:0895-4798
1095-7162
DOI:10.1137/S0895479800367439