On the Lyapunov equation with the state matrix in companion form

Sintesi Tecnica: Sull'equazione di Lyapunov con matrice di stato in forma companion

Enunciato del Problema
Questo articolo investiga la nonnegatività elemento per elemento dell'unica soluzione $X$ dell'equazione di Lyapunov in tempo continuo:
$$XA + A^\top X = -Q$$
dove $A \in \mathbb{R}^{n \times n}$ è una matrice companion di Hurwitz associata a un polinomio $a(s)$, e $Q = Q^\top \succeq 0$ è una matrice semidefinita positiva.

Sebbene la teoria standard di Lyapunov garantisca che $X$ sia simmetrica e semidefinita positiva (nell'ordine di Loewner), questa proprietà non implica che $X$ sia nonnegativa elemento per elemento ($X_{ij} \geq 0$ per ogni $i,j$). L'articolo affronta la Congettura 1.1, proposta dal Prof. Luigi Fortuna, la quale postula che se $A$ è una matrice companion di Hurwitz e $Q \succeq 0$, allora la soluzione $X$ è nonnegativa elemento per elemento. Questa proprietà è distinta dalla semidefinitività positiva ed è di particolare interesse nella teoria dei sistemi positivi, dove la nonnegatività componente per componente assicura la preservazione di grandezze fisiche come masse, concentrazioni o popolazioni.

Metodologia
La strategia di dimostrazione si basa sulla riduzione del problema generale a una specifica struttura algebrica che coinvolge polinomi simmetrici elementari e matrici di tipo Cauchy. La metodologia procede attraverso i seguenti passaggi:

1. Semplificazione del Problema:

   • Riduzione al Rango-1: Per la linearità dell'equazione di Lyapunov e la decomposizione di qualsiasi matrice semidefinita positiva $Q$ in una somma di matrici di rango-1 ($Q = \sum q_i q_i^\top$), gli autori assumono senza perdita di generalità che $Q = qq^\top$.
   • Diagonalizzabilità: Gli autori sostengono che la soluzione $X$ sia una funzione continua di $A$. Pertanto, è sufficiente dimostrare il risultato per matrici companion diagonalizzabili (dove gli autovalori sono distinti). I casi non diagonalizzabili possono essere approssimati da quelli diagonalizzabili con coefficienti arbitrariamente vicini.

2. Decomposizione Spettrale:

   • Assumendo che $A$ sia diagonalizzabile con autovalori reali distinti $-\lambda_i$ (ordinati $0 < \lambda_1 < \dots < \lambda_n$), gli autori utilizzano gli autovettori di $A^\top$. Tali autovettori, indicati con $v_j$, sono costruiti utilizzando i polinomi simmetrici elementari $\sigma_m^{(j)}$ degli autovalori.
   • Il vettore di forzante $q$ viene espanso in questa base propria: $q = \sum a_i v_i$.

3. Rappresentazione Integrale e Forme Quadratiche:

   • Utilizzando la formula integrale per la soluzione di Lyapunov, $X = \int_0^\infty e^{A^\top t} Q e^{At} dt$, gli autori derivano un'espressione esplicita per le entrate $X_{kl}$.
   • L'entrata $X_{kl}$ è mostrata essere una forma quadratica $a^\top M a$, dove la matrice $M$ ha elementi definiti da:
     $$M_{ij} = \frac{\sigma_{n-k}^{(i)}\sigma_{n-l}^{(j)} + \sigma_{n-k}^{(j)}\sigma_{n-l}^{(i)}}{\lambda_i + \lambda_j}$$
   • Questa matrice $M$ è fattorizzata come $M = \Delta C \Delta$, dove $\Delta$ è una matrice diagonale e $C$ è una matrice di tipo Cauchy con entrate $C_{ij} = \frac{r_i + r_j}{\lambda_i + \lambda_j}$.

4. Verifica della Positività:

   • La nonnegatività di $X_{kl}$ segue se la matrice $M$ è semidefinita positiva. Poiché $\Delta$ è nonsingolare, ciò equivale a dimostrare che $C \succeq 0$.
   • Gli autori citano risultati consolidati [8] che confermano che tali matrici di tipo Cauchy sono semidefinite positive.

Contributi Chiave e Risultati

• Teorema 1 (Risultato Principale): Il saggio dimostra la Congettura 1.1 sotto la specifica condizione che la matrice companion $A$ abbia solo autovalori reali. Sotto tale assunzione, la soluzione unica $X$ è nonnegativa elemento per elemento.
• Analisi della Total Positività: Gli autori investigano una proprietà più forte: la totale nonnegatività (tutti i minori sono nonnegativi).
  • Forniscono un controesempio che mostra come $X$ non sia generalmente totalmente nonnegativa, anche quando gli autovalori sono reali.
  • Teorema 2: Stabiliscono una condizione sufficiente per la totale positività: se il vettore di forzante $q$ ha un'espansione nonnegativa nella base propria di $A^\top$ (specificamente, $q = \sum \alpha_j v_j$ con $\alpha_j \geq 0$), allora $X$ è totalmente positiva. Ciò è dimostrato fattorizzando $X$ in un prodotto di matrici totalmente positive (che coinvolgono matrici di Vandermonde e matrici di Cauchy).

Significato e Rivendicazioni
Il saggio sostiene di risolvere una congettura di lunga data riguardante le proprietà strutturali delle soluzioni di Lyapunov per matrici companion, specificamente nell'ambito dei sistemi positivi.

• Motivazione: La nonnegatività elemento per elemento di $X$ è motivata da applicazioni nel controllo ottimo per sistemi positivi. Nello specifico, per sistemi a singolo ingresso con matrici di ingresso nonnegative, una soluzione elemento per elemento nonnegativa dell'equazione di Riccati algebrica associata (riscritta come equazione di Lyapunov) implica che il regolatore ottimo può essere implementato utilizzando hardware con soli alimentatori di tensione positivi.
• Ambito e Limiti: Gli autori sono cauti riguardo alla generalità della loro dimostrazione. Essi dichiarano esplicitamente che l'assunzione di autovalori reali è una restrizione. Il saggio lascia aperta la questione se la congettura valga per matrici companion di Hurwitz con autovalori complessi.
• Valore Teorico: Il lavoro connette la teoria di Lyapunov con la teoria delle matrici di tipo Cauchy e dei polinomi simmetrici elementari, dimostrando che la nonnegatività elemento per elemento di $X$ deriva dalla semidefinitività positiva di una specifica classe di matrici duali alle matrici di Kwong (o anti-Loewner).

In sintesi, il saggio fornisce una dimostrazione rigorosa che le soluzioni di Lyapunov per matrici companion di Hurwitz con spettri reali sono nonnegative elemento per elemento, chiarendo al contempo che la totale positività richiede vincoli aggiuntivi sul vettore di forzante.