Iterative methods for the delay Lyapunov equation with T-Sylvester preconditioning

The delay Lyapunov equation is an important matrix boundary-value problem which arises as an analogue of the Lyapunov equation in the study of time-delay systems x˙(t)=A0x(t)+A1x(t−τ)+B0u(t). We propose a new algorithm for the solution of the delay Lyapunov equation. Our method is based on the fact...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:Applied numerical mathematics 2019-01, Vol.135, p.173-185
Hauptverfasser: Jarlebring, Elias, Poloni, Federico
Format: Artikel
Sprache:eng
Schlagworte:
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Beschreibung
Zusammenfassung:The delay Lyapunov equation is an important matrix boundary-value problem which arises as an analogue of the Lyapunov equation in the study of time-delay systems x˙(t)=A0x(t)+A1x(t−τ)+B0u(t). We propose a new algorithm for the solution of the delay Lyapunov equation. Our method is based on the fact that the delay Lyapunov equation can be expressed as a linear system of equations, whose unknown is the value U(τ/2)∈Rn×n, i.e., the delay Lyapunov matrix at time τ/2. This linear matrix equation with n2 unknowns is solved by adapting a preconditioned iterative method such as GMRES. The action of the n2×n2 matrix associated to this linear system can be computed by solving a coupled matrix initial-value problem. A preconditioner for the iterative method is proposed based on solving a T-Sylvester equation MX+XTN=C, for which there are methods available in the literature. We prove that the preconditioner is effective under certain assumptions. The efficiency of the approach is illustrated by applying it to a time-delay system stemming from the discretization of a partial differential equation with delay. Approximate solutions to this problem can be obtained for problems of size up to n≈1000, i.e., a linear system with n2≈106 unknowns, a dimension which is outside of the capabilities of the other existing methods for the delay Lyapunov equation.
ISSN:0168-9274
1873-5460
1873-5460
DOI:10.1016/j.apnum.2018.08.011