On the Lyapunov equation with the state matrix in companion form

Technische Zusammenfassung: Über die Lyapunov-Gleichung mit der Zustandsmatrix in Begleiterform

Problemstellung
Diese Arbeit untersucht die eintragsweise Nichtnegativität der eindeutigen Lösung $X$ der kontinuierlichen Zeit-Lyapunov-Gleichung:
$$XA + A^\top X = -Q$$
wobei $A \in \mathbb{R}^{n \times n}$ eine Hurwitz-Begleitmatrix (Companion Matrix) eines Polynoms $a(s)$ ist und $Q = Q^\top \succeq 0$ eine positiv semidefinite Matrix ist.

Obwohl die Standard-Lyapunov-Theorie garantiert, dass $X$ symmetrisch und positiv semidefinit (im Sinne der Loewner-Ordnung) ist, impliziert diese Eigenschaft nicht, dass $X$ eintragswiese nichtnegativ ist ($X_{ij} \geq 0$ für alle $i,j$). Die Arbeit adressiert Konjektur 1.1, vorgeschlagen von Prof. Luigi Fortuna, die postuliert, dass falls $A$ eine Hurwitz-Begleitmatrix und $Q \succeq 0$ ist, dann die Lösung $X$ eintragswiese nichtnegativ ist. Diese Eigenschaft unterscheidet sich von der positiven Semidefinitheit und ist von besonderem Interesse in der Theorie positiver Systeme, in denen die komponentenweise Nichtnegativität die Erhaltung physikalischer Größen wie Massen, Konzentrationen oder Populationen sicherstellt.

Methodik
Die Beweisstrategie beruht auf der Reduktion des allgemeinen Problems auf eine spezifische algebraische Struktur unter Verwendung elementarer symmetrischer Polynome und Cauchy-ähnlicher Matrizen. Die Methodik erfolgt durch die folgenden Schritte:

1. Vereinfachung des Problems:

   • Rang-1-Reduktion: Durch die Linearität der Lyapunov-Gleichung und die Zerlegung einer beliebigen positiv semidefiniten Matrix $Q$ in eine Summe von Rang-1-Matrizen ($Q = \sum q_i q_i^\top$), nehmen die Autoren ohne Verlust der Allgemeingültigkeit an, dass $Q = qq^\top$ gilt.
   • Diagonalisierbarkeit: Die Autoren argumentieren, dass die Lösung $X$ eine stetige Funktion von $A$ ist. Daher genügt es zu beweisen, dass das Ergebnis für diagonalisierbare Begleitmatrizen (bei denen die Eigenwerte verschieden sind) gilt. Nicht-diagonalisierbare Fälle können durch diagonalisierbare Matrizen mit beliebig nahen Koeffizienten approximiert werden.

2. Spektrale Zerlegung:

   • Unter der Annahme, dass $A$ diagonalisierbar mit distinkten reellen Eigenwerten $-\lambda_i$ (geordnet $0 < \lambda_1 < \dots < \lambda_n$) ist, nutzen die Autoren die Eigenvektoren von $A^\top$. Diese Eigenvektoren, bezeichnet als $v_j$, werden unter Verwendung elementarer symmetrischer Polynome $\sigma_m^{(j)}$ der Eigenwerte konstruiert.
   • Der erregende Vektor $q$ wird in dieser Eigenbasis expandiert: $q = \sum a_i v_i$.

3. Integraldarstellung und quadratische Formen:

   • Unter Verwendung der Integralformel für die Lyapunov-Lösung, $X = \int_0^\infty e^{A^\top t} Q e^{At} dt$, leiten die Autoren einen expliziten Ausdruck für die Einträge $X_{kl}$ ab.
   • Der Eintrag $X_{kl}$ wird als quadratische Form $a^\top M a$ dargestellt, wobei die Matrix $M$ Elemente besitzt, die definiert sind durch:
     $$M_{ij} = \frac{\sigma_{n-k}^{(i)}\sigma_{n-l}^{(j)} + \sigma_{n-k}^{(j)}\sigma_{n-l}^{(i)}}{\lambda_i + \lambda_j}$$
   • Diese Matrix $M$ wird als $M = \Delta C \Delta$ faktorisiert, wobei $\Delta$ eine Diagonalmatrix und $C$ eine Cauchy-ähnliche Matrix mit Einträgen $C_{ij} = \frac{r_i + r_j}{\lambda_i + \lambda_j}$ ist.

4. Verifizierung der Positivität:

   • Die Nichtnegativität von $X_{kl}$ folgt, wenn die Matrix $M$ positiv semidefinit ist. Da $\Delta$ nicht singulär ist, ist dies äquivalent zum Beweis, dass $C \succeq 0$.
   • Die Autoren zitieren etablierte Ergebnisse [8], die bestätigen, dass solche Cauchy-ähnlichen Matrizen positiv semidefinit sind.

Wesentliche Beiträge und Ergebnisse

• Theorem 1 (Hauptresultat): Die Arbeit beweist Konjektur 1.1 unter der spezifischen Bedingung, dass die Begleitmatrix $A$ nur reelle Eigenwerte besitzt. Unter dieser Annahme ist die eindeutige Lösung $X$ eintragswiese nichtnegativ.
• Analyse der totalen Positivität: Die Autoren untersuchen eine stärkere Eigenschaft: totale Nichtnegativität (alle Minorante sind nichtnegativ).
  • Sie liefern ein Gegenbeispiel, das zeigt, dass $X$ im Allgemeinen nicht total nichtnegativ ist, selbst wenn die Eigenwerte reell sind.
  • Theorem 2: Sie etablieren eine hinreichende Bedingung für totale Positivität: Wenn der erregende Vektor $q$ eine nichtnegative Expansion in der Eigenbasis von $A^\top$ besitzt (speziell $q = \sum \alpha_j v_j$ mit $\alpha_j \geq 0$), dann ist $X$ total positiv. Dies wird durch die Faktorisierung von $X$ in ein Produkt total positiver Matrizen (unter Beteiligung von Vandermonde-Matrizen und Cauchy-Matrizen) bewiesen.

Bedeutung und Ansprüche
Die Arbeit behauptet, eine langjährige Vermutung bezüglich der strukturellen Eigenschaften von Lyapunov-Lösungen für Begleitmatrizen zu klären, speziell im Bereich der positiven Systeme.

• Motivation: Die eintragswiese Nichtnegativität von $X$ wird durch Anwendungen in der optimalen Steuerung positiver Systeme motiviert. Speziell für Single-Input-Systeme mit nichtnegativen Eingangsmatrizen impliziert eine eintragswiese nichtnegative Lösung der zugehörigen algebraischen Riccati-Gleichung (umgeschrieben als Lyapunov-Gleichung), dass der optimale Regler mit Hardware implementiert werden kann, die nur positive Spannungsversorgungen nutzt.
• Umfang und Einschränkungen: Die Autoren äußern sich bescheiden hinsichtlich der Allgemeingültigkeit ihres Beweises. Sie geben explizit an, dass die Annahme reeller Eigenwerte eine Einschränkung darstellt. Die Arbeit lässt die Frage offen, ob die Konjektur auch für Hurwitz-Begleitmatrizen mit komplexen Eigenwerten gilt.
• Theoretischer Wert: Die Arbeit verbindet die Lyapunov-Theorie mit der Theorie der Cauchy-ähnlichen Matrizen und der elementaren symmetrischen Polynome und zeigt auf, dass die eintragswiese Nichtnegativität von $X$ aus der positiven Semidefinitheit einer spezifischen Klasse von Matrizen resultiert, die dual zu Kwong- (oder Anti-Loewner-) Matrizen sind.

Zusammenfassend liefert die Arbeit einen rigorosen Beweis dafür, dass Lyapunov-Lösungen für Hurwitz-Begleitmatrizen mit reellem Spektrum eintragswiese nichtnegativ sind, während sie gleichzeitig klärt, dass totale Positivität zusätzliche Bedingungen an den erregenden Vektor stellt.