Perturbation of monic matrix polynomials

1. Problématique et Contexte

L'article étudie la stabilité spectrale des polynômes matriciels soumis à des perturbations structurées de leurs coefficients. Soit un polynôme matriciel de degré $d$ défini par :
$$P_u(\lambda) = A_d(u)\lambda^d + A_{d-1}(u)\lambda^{d-1} + \dots + A_0(u)$$
où $\lambda \in \mathbb{C}$ et les coefficients matriciels $A_k(u) \in \mathbb{C}^{n \times n}$ dépendent continûment et de manière semi-algébrique d'un vecteur de paramètres $u \in \mathbb{C}^p$.

Les auteurs se concentrent sur le cas monique, c'est-à-dire où le coefficient dominant est l'identité ($A_d(u) \equiv I_n$), bien que les résultats soient extensibles au cas où $A_d(u)$ est localement non singulier.

L'objectif principal est d'analyser la continuité (spécifiquement la continuité Hölderienne) des objets spectraux suivants en fonction du paramètre $u$ :

1. Le spectre $\sigma(u) = \{ \lambda \in \mathbb{C} \mid \det P_u(\lambda) = 0 \}$.
2. Le $\epsilon$-pseudospectre $\Lambda_\epsilon(u)$, défini par les perturbations admissibles du vecteur de paramètres $u$ (perturbations structurées), par opposition aux perturbations arbitraires de norme bornée.
3. La valeur numérique (numerical range) $W(u) = \{ \lambda \in \mathbb{C} \mid x^* P_u(\lambda) x = 0, \|x\|=1 \}$.
4. La valeur numérique conjointe (joint numerical range) $JW(u)$ associée à la $(d+1)$-tuplet des matrices coefficients.
5. Les paires de Jordan (Jordan pairs) et le rayon spectral.

2. Méthodologie

La démarche méthodologique repose sur l'application rigoureuse de la géométrie semi-algébrique et de l'analyse réelle complexe.

• Cadre Semi-Algébrique : Les auteurs exploitent le fait que les coefficients $A_k(u)$ sont des fonctions semi-algébriques. Cela implique que les ensembles définis par des équations et inégalités polynomiales (comme le spectre, le pseudospectre, etc.) sont eux-mêmes des ensembles semi-algébriques.
• Théorème de Tarski-Seidenberg : Ce théorème fondamental garantit que l'image d'un ensemble semi-algébrique par une projection (comme la projection du spectre sur le plan complexe) reste semi-algébrique. Cela permet d'affirmer que les objets spectraux étudiés possèdent une structure géométrique bien définie (nombre fini de composantes connexes, frontières algébriques).
• Inégalité de Łojasiewicz : Pour établir la continuité, les auteurs utilisent l'inégalité de Łojasiewicz. Celle-ci stipule que pour toute application semi-algébrique continue définie sur un compact, la fonction est Hölderienne. Cela permet de quantifier la stabilité des ensembles spectraux face aux variations du paramètre $u$.
• Décomposition Cellulaire : L'article utilise le théorème de décomposition cellulaire pour partitionner l'espace des paramètres $\mathbb{C}^p$ en une union finie de sous-variétés analytiques semi-algébriques. Sur chaque variété, la structure spectrale devient régulière (analytique).

3. Résultats Principaux

Les résultats sont présentés en plusieurs étapes, allant de la continuité globale à la régularité analytique locale.

A. Continuité Hölderienne des Applications à Valeur Ensemble

Sous l'hypothèse que les coefficients sont continus et semi-algébriques, les auteurs démontrent que les applications à valeur ensemble suivantes sont Hölder-continues par rapport au paramètre $u$ (localement) :

• Spectre $\sigma(u)$ : L'application $u \mapsto \sigma(u)$ est Hölder-continue.
• $\epsilon$-Pseudospectre $\Lambda_\epsilon(u)$ : L'application est Hölder-continue, même lorsque $\epsilon$ varie. De plus, la frontière $\partial \Lambda_\epsilon(u)$ est contenue dans une courbe algébrique et ne possède qu'un nombre fini de singularités (points de rebroussement ou d'auto-intersection).
• Valeur Numérique $W(u)$ et sa frontière $\partial W(u)$ : Ces applications sont également Hölder-continues. La frontière $\partial W(u)$ est une courbe algébrique avec un nombre fini de points non différentiables.
• Valeur Numérique Conjointe $JW(u)$ : Cette application est Hölder-continue.

B. Continuité du Rayon Spectral

Le rayon spectral $r(u) = \max \{ |\lambda| : \lambda \in \sigma(u) \}$ est une fonction continue et semi-algébrique, donc Hölder-continue sur les compacts.

C. Régularité Analytique Générique (Résultat Fort)

Un résultat central de l'article est la décomposition générique de l'espace des paramètres. Il existe un ouvert dense et semi-algébrique $U \subset \mathbb{C}^p$ tel que :

1. Sur $U$, les applications $\sigma, \Lambda_\epsilon, W, \partial W$ sont Hölder-continues avec des exposants explicites (dépendant du degré $d$ et de la taille $n$).
2. L'application $JW$ est Lipschitzienne sur $U$.
3. Le rayon spectral $r(u)$ est une fonction analytique sur $U$.

D. Stabilité des Paires de Jordan

Pour les polynômes moniques, l'espace des paramètres $\mathbb{C}^p$ peut être décomposé en une union finie de sous-variétés analytiques semi-algébriques $M_1, \dots, M_N$. Sur chaque variété $M_t$ :

• Le nombre de blocs de Jordan et leurs tailles sont constants (indépendants de $u$).
• Les valeurs propres distinctes $\lambda_k(u)$ sont des fonctions analytiques de $u$.
• Les blocs de la paire de Jordan $X_u$ et $J_u$ dépendent analytiquement de $u$.

Cela signifie que les singularités (changements de structure de Jordan) ne se produisent que sur des ensembles de mesure nulle (les frontières entre les variétés).

4. Contributions Clés

1. Généralisation aux Perturbations Structurées : Contrairement à la littérature classique qui étudie souvent les perturbations de norme arbitraire, ce travail se concentre sur les perturbations induites par la variation d'un vecteur de paramètres physique ou de conception ($u$). Cela rend les résultats plus pertinents pour les applications en ingénierie et en contrôle.
2. Unification par la Géométrie Semi-Algébrique : L'article fournit un cadre unifié pour analyser la stabilité du spectre, du pseudospectre et des valeurs numériques en utilisant les propriétés de stabilité des ensembles semi-algébriques sous projection et intersection.
3. Quantification de la Stabilité : L'obtention d'exposants de Hölder explicites (bien que dépendant de la dimension) et la preuve de la régularité analytique générique offrent des garanties quantitatives sur la stabilité des systèmes.
4. Analyse de la Frontière : La caractérisation géométrique des frontières du pseudospectre et de la valeur numérique (courbes algébriques, nombre fini de singularités) est une contribution significative pour la compréhension visuelle et computationnelle de ces ensembles.

5. Signification et Impact

Ce travail est significatif pour plusieurs domaines :

• Théorie Spectrale : Il renforce la compréhension de la stabilité des polynômes matriciels au-delà du cas linéaire (matrices simples).
• Analyse Numérique et Contrôle : La notion de pseudospectre structuré est cruciale pour l'analyse de la sensibilité des systèmes dynamiques et des problèmes de valeurs propres généralisés. Les résultats garantissent que de petites variations dans les paramètres de conception (géométrie, matériaux) ne provoquent pas de changements catastrophiques dans le spectre, sauf sur des ensembles de paramètres très spécifiques.
• Géométrie Algébrique Appliquée : L'article démontre la puissance des outils de la géométrie semi-algébrique (Tarski-Seidenberg, décomposition cellulaire) pour résoudre des problèmes d'analyse fonctionnelle et de stabilité des systèmes.

En résumé, l'article établit que, pour les polynômes matriciels moniques à coefficients semi-algébriques, les objets spectraux sont non seulement stables (Hölder-continus), mais qu'ils présentent une régularité analytique exceptionnelle sur un ensemble de paramètres générique, permettant une analyse prédictive fiable des systèmes perturbés.