On the Lyapunov equation with the state matrix in companion form

Technische Samenvatting: Over de Lyapunov-vergelijking met de Toestandmatrix in Companionvorm

Probleemstelling
Dit artikel onderzoekt de componentgewijze niet-negativiteit van de unieke oplossing $X$ van de continue-tijd Lyapunov-vergelijking:
$$XA + A^\top X = -Q$$
waarbij $A \in \mathbb{R}^{n \times n}$ een Hurwitz companionmatrix is behorend bij een polynoom $a(s)$, en $Q = Q^\top \succeq 0$ een positief semidefiniete matrix is.

Hoewel de standaard Lyapunov-theorie garandeert dat $X$ symmetrisch en positief semidefiniet is (in de Loewner-orde), impliceert deze eigenschap niet dat $X$ componentgewijs niet-negatief is ($X_{ij} \geq 0$ voor alle $i,j$). Dit artikel behandelt Conjectuur 1.1, voorgesteld door Prof. Luigi Fortuna, die stelt dat indien $A$ een Hurwitz companionmatrix is en $Q \succeq 0$, dan de oplossing $X$ componentgewijs niet-negatief is. Deze eigenschap is verschillend van positieve semidefinitie en is van bijzonder belang binnen de theorie van positieve systemen, waar componentgewijze niet-negativiteit de bewaring van fysieke grootheden zoals massa's, concentraties of populaties waarborgt.

Methodologie
De bewijsstrategie berust op het reduceren van het algemene probleem tot een specifieke algebraïsche structuur bestaande uit elementaire symmetrische polynomen en Cauchy-achtige matrices. De methodologie verloopt via de volgende stappen:

1. Vereenvoudiging van het probleem:

   • Rang-1 Reductie: Door de lineariteit van de Lyapunov-vergelijking en de decompositie van elke positief semidefiniete $Q$ in een som van rang-1 matrices ($Q = \sum q_i q_i^\top$), nemen de auteurs zonder algemeen verlies aan dat $Q = qq^\top$.
   • Diagonaliseerbaarheid: De auteurs argumenteren dat de oplossing $X$ een continue functie is van $A$. Het volstaat daarom om het resultaat te bewijzen voor diagonaliseerbare companionmatrices (waarbij de eigenwaarden verschillend zijn). Niet-diagonaliseerbare gevallen kunnen worden benaderd door diagonaliseerbare matrices met willekeurig dichtbij gelegen coëfficiënten.

2. Spectrale Decompositie:

   • Uitgaande van het geval dat $A$ diagonaliseerbaar is met verschillende reële eigenwaarden $-\lambda_i$ (geordend $0 < \lambda_1 < \dots < \lambda_n$), maken de auteurs gebruik van de eigenvectoren van $A^\top$. Deze eigenvectoren, aangeduid als $v_j$, worden geconstrueerd met behulp van elementaire symmetrische polynomen $\sigma_m^{(j)}$ van de eigenwaarden.
   • De forcingsvector $q$ wordt uitgebreid in deze eigenbasis: $q = \sum a_i v_i$.

3. Integrale Representatie en Kwadratische Vormen:

   • Met behulp van de integraalformule voor de Lyapunov-oplossing, $X = \int_0^\infty e^{A^\top t} Q e^{At} dt$, leiden de auteurs een expliciete uitdrukking af voor de elementen $X_{kl}$.
   • Het element $X_{kl}$ wordt getoond te zijn als een kwadratische vorm $a^\top M a$, waarbij de matrix $M$ elementen heeft gedefinieerd door:
     $$M_{ij} = \frac{\sigma_{n-k}^{(i)}\sigma_{n-l}^{(j)} + \sigma_{n-k}^{(j)}\sigma_{n-l}^{(i)}}{\lambda_i + \lambda_j}$$
   • Deze matrix $M$ wordt gefactoriseerd als $M = \Delta C \Delta$, waarbij $\Delta$ een diagonale matrix is en $C$ een Cauchy-achtige matrix met elementen $C_{ij} = \frac{r_i + r_j}{\lambda_i + \lambda_j}$.

4. Verificatie van Positiviteit:

   • De niet-negativiteit van $X_{kl}$ volgt indien de matrix $M$ positief semidefiniet is. Aangezien $\Delta$ niet-singulier is, is dit equivalent aan het bewijzen dat $C \succeq 0$.
   • De auteurs citeren gevestigde resultaten [8] die bevestigen dat dergelijke Cauchy-achtige matrices positief semidefiniet zijn.

Belangrijkste Bijdragen en Resultaten

• Stelling 1 (Hoofdresultaat): Het artikel bewijst Conjectuur 1.1 onder de specifieke voorwaarde dat de companionmatrix $A$ alleen reële eigenwaarden heeft. Onder deze aanname is de unieke oplossing $X$ componentgewijs niet-negatief.
• Analyse van Totale Positiviteit: De auteurs onderzoeken een sterkere eigenschap: totale niet-negativiteit (alle minor-determinanten zijn niet-negatief).
  • Zij leveren een tegenvoorbeeld dat aantoont dat $X$ over het algemeen niet totaal niet-negatief is, zelfs wanneer de eigenwaarden reëel zijn.
  • Stelling 2: Zij stellen een voldoende voorwaarde vast voor totale positiviteit: indien de forcingsvector $q$ een niet-negatieve expansie heeft in de eigenbasis van $A^\top$ (specifiek, $q = \sum \alpha_j v_j$ met $\alpha_j \geq 0$), dan is $X$ totaal positief. Dit wordt bewezen door $X$ te factoriseren in een product van totaal positieve matrices (met betrekking tot Vandermonde-matrices en Cauchy-matrices).

Betekenis en Claims
Het artikel beweert een langlopende conjectuur op te lossen met betrekking tot de structurele eigenschappen van Lyapunov-oplossingen voor companionmatrices, specifiek binnen het domein van positieve systemen.

• Motivatie: De componentgewijze niet-negativiteit van $X$ wordt gemotiveerd door toepassingen in optimale regeling voor positieve systemen. Specifiek, voor single-input systemen met niet-negatieve inputmatrices, impliceert een componentgewijs niet-negatieve oplossing van de bijbehorende algebraïsche Riccati-vergelijking (herleid tot een Lyapunov-vergelijking) dat de optimale regelaar geïmplementeerd kan worden met hardware die slechts positieve voedingsspanningen gebruikt.
• Reikwijdte en Beperkingen: De auteurs zijn bescheiden over de algemeenheid van hun bewijs. Zij stellen expliciet dat de aanname van reële eigenwaarden een restrictie is. Het artikel laat de vraag open of de conjectuur standhoudt voor Hurwitz companionmatrices met complexe eigenwaarden.
• Theoretische Waarde: Het werk verbindt de Lyapunov-theorie met de theorie van Cauchy-achtige matrices en elementaire symmetrische polynomen, en demonstreert dat de componentgewijze niet-negativiteit van $X$ voortkomt uit de positieve semidefinitie van een specifieke klasse matrices die duaal zijn aan Kwong- (of anti-Loewner-) matrices.

Samenvattend biedt het artikel een rigoureus bewijs dat de Lyapunov-oplossingen voor Hurwitz companionmatrices met een reëel spectrum componentgewijs niet-negatief zijn, terwijl het verheldert dat totale positiviteit aanvullende beperkingen op de forcingsvector vereist.