On the Lyapunov equation with the state matrix in companion form

Résumé Technique : Sur l'équation de Lyapunov avec la matrice d'état sous forme compagnon

Énoncé du problème
Cet article étudie la nonnégativité composante par composante de l'unique solution $X$ de l'équation de Lyapunov en temps continu :
$$XA + A^\top X = -Q$$
où $A \in \mathbb{R}^{n \times n}$ est une matrice compagnon de Hurwitz associée à un polynôme $a(s)$, et $Q = Q^\top \succeq 0$ est une matrice semi-définie positive.

Bien que la théorie de Lyapunov standard garantisse que $X$ est symétrique et semi-définie positive (dans l'ordre de Loewner), cette propriété n'implique pas que $X$ est nonnégative composante par composante ($X_{ij} \geq 0$ pour tout $i,j$). L'article traite de la Conjecture 1.1, proposée par le Professeur Luigi Fortuna, qui postule que si $A$ est une matrice compagnon de Hurwitz et $Q \succeq 0$, alors la solution $X$ est nonnégative composante par composante. Cette propriété est distincte de la semi-définie positivité et présente un intérêt particulier dans la théorie des systèmes positifs, où la nonnégativité composante garantit la préservation des quantités physiques telles que les masses, les concentrations ou les populations.

Méthodologie
La stratégie de preuve repose sur la réduction du problème général à une structure algébrique spécifique impliquant des polynômes symétriques élémentaires et des matrices de type Cauchy. La méthodologie procède à travers les étapes suivantes :

1. Simplification du problème :

   • Réduction au rang 1 : Par la linéarité de l'équation de Lyapunov et la décomposition de toute matrice semi-définie positive $Q$ en une somme de matrices de rang 1 ($Q = \sum q_i q_i^\top$), les auteurs supposent sans perte de généralité que $Q = qq^\top$.
   • Diagonalisabilité : Les auteurs soutiennent que la solution $X$ est une fonction continue de $A$. Il suffit donc de prouver le résultat pour les matrices compagnon diagonalisables (où les valeurs propres sont distinctes). Les cas non diagonalisables peuvent être approchés par des cas diagonalisables avec des coefficients arbitrairement proches.

2. Décomposition spectrale :

   • En supposant que $A$ est diagonalisable avec des valeurs propres réelles distinctes $-\lambda_i$ (ordonnées $0 < \lambda_1 < \dots < \lambda_n$), les auteurs utilisent les vecteurs propres de $A^\top$. Ces vecteurs propres, notés $v_j$, sont construits à l'aide de polynômes symétriques élémentaires $\sigma_m^{(j)}$ des valeurs propres.
   • Le vecteur de forçage $q$ est développé dans cette base propre : $q = \sum a_i v_i$.

3. Représentation intégrale et formes quadratiques :

   • En utilisant la formule intégrale pour la solution de Lyapunov, $X = \int_0^\infty e^{A^\top t} Q e^{At} dt$, les auteurs dérivent une expression explicite pour les entrées $X_{kl}$.
   • L'entrée $X_{kl}$ est montrée comme étant une forme quadratique $a^\top M a$, où la matrice $M$ possède des éléments définis par :
     $$M_{ij} = \frac{\sigma_{n-k}^{(i)}\sigma_{n-l}^{(j)} + \sigma_{n-k}^{(j)}\sigma_{n-l}^{(i)}}{\lambda_i + \lambda_j}$$
   • Cette matrice $M$ est factorisée en $M = \Delta C \Delta$, où $\Delta$ est une matrice diagonale et $C$ est une matrice de type Cauchy dont les entrées sont $C_{ij} = \frac{r_i + r_j}{\lambda_i + \lambda_j}$.

4. Vérification de la positivité :

   • La nonnégativité de $X_{kl}$ découle du fait que la matrice $M$ est semi-définie positive. Puisque $\Delta$ est non singulière, cela équivaut à prouver que $C \succeq 0$.
   • Les auteurs citent des résultats établis [8] confirmant que de telles matrices de type Cauchy sont semi-définies positives.

Contributions clés et résultats

• Théorème 1 (Résultat principal) : L'article prouve la Conjecture 1.1 sous la condition spécifique que la matrice compagnon $A$ possède uniquement des valeurs propres réelles. Sous cette hypothèse, la solution unique $X$ est nonnégative composante par composante.
• Analyse de la positivité totale : Les auteurs étudient une propriété plus forte : la nonnégativité totale (tous les mineurs sont nonnégatifs).
  • Ils fournissent un contre-exemple montrant que $X$ n'est pas généralement totalement nonnégative, même lorsque les valeurs propres sont réelles.
  • Théorème 2 : Ils établissent une condition suffisante pour la positivité totale : si le vecteur de forçage $q$ possède une expansion nonnégative dans la base propre de $A^\top$ (spécifiquement, $q = \sum \alpha_j v_j$ avec $\alpha_j \geq 0$), alors $X$ est totalement positive. Ceci est prouvé en factorisant $X$ en un produit de matrices totalement positives (impliquant des matrices de Vandermonde et de Cauchy).

Signification et prétentions
L'article prétend résoudre une conjecture de longue date concernant les propriétés structurelles des solutions de Lyapunov pour les matrices compagnon, spécifiquement dans le domaine des systèmes positifs.

• Motivation : La nonnégativité composante de $X$ est motivée par des applications en contrôle optimal pour les systèmes positifs. Plus précisément, pour les systèmes à commande unique avec des matrices d'entrée nonnégatives, une solution nonnégative composante de l'équation de Riccati algébrique associée (réécrite comme une équation de Lyapunov) implique que le régulateur optimal peut être implémenté via du matériel ne nécessitant que des alimentations de tension positives.
• Portée et limites : Les auteurs font preuve de modestie quant à la généralité de leur preuve. Ils déclarent explicitement que l'hypothèse de valeurs propres réelles est une restriction. L'article laisse ouverte la question de savoir si la conjecture est vérifiée pour les matrices compagnon de Hurwitz possédant des valeurs propres complexes.
• Valeur théorique : Le travail relie la théorie de Lyapunov à la théorie des matrices de type Cauchy et des polynômes symétriques élémentaires, démontrant que la nonnégativité composante de $X$ provient de la semi-définie positivité d'une classe spécifique de matrices duales aux matrices de Kwong (ou anti-Loewner).

En résumé, l'article fournit une preuve rigoureuse que les solutions de Lyapunov pour les matrices compagnon de Hurwitz à spectre réel sont nonnégatives composante par composante, tout en précisant que la positivité totale nécessite des contraintes supplémentaires sur le vecteur de forçage.