Accueil🇫🇷Chercher

Équation de Sylvester

En mathématiques, dans le domaine de la théorie du contrôle, une équation de Sylvester est une équation matricielle de la forme[1] :

.

Les matrices , et sont données ; le problème est de déterminer les matrices qui satisfont cette équation. Les matrices sont supposées à coefficients complexes. Les matrices sont de tailles appropriées, par exemple elles sont toutes des matrices carrées de même taille, ou plus généralement, et sont des matrices carrées de taille et respectivement, et et sont à lignes et colonnes.

Une équation de Sylvester a une solution unique en si et seulement si et n'ont pas de valeur propre commune. Plus généralement, l'équation a été considérée comme une équation d'opérateurs bornés dans un espace de Banach (éventuellement de dimension infinie). Dans ce cas, la condition d'unicité d'une solution en est quasiment la même : il existe une solution unique en exactement lorsque les spectres de et de sont disjoints[2].

Existence et unicité des solutions

En utilisant le produit de Kronecker et l'opérateur de vectorisation (en) , on peut réécrire l'équation de Sylvester sous la forme

est de taille , est de taille , de taille et est la matrice d'identité de taille . Sous cette forme, l'équation peut être vue comme un système linéaire de taille [3].

Théorème Étant données des matrices et , l'équation de Sylvester : a une solution unique quelle que soit si et seulement si et nont pas de valeur propre commune.

La non-singularité de dans la partie (i) de la preuve ci-dessus peut également être démontrée par l' identité de Bézout pour des polynômes premiers entre eux. Soit en effet le polynôme caractéristique de . Comme et n'ont pas de valeur propre commune, et sont premiers entre eux. Il existe donc des polynômes et tel que . Par le théorème de Cayley-Hamilton, . Ainsi , ce qui implique que est non singulier.

Le théorème reste vrai pour les matrices réelles sous réserve de considérer leurs valeurs propres complexes. La preuve de la partie directe est toujours valable ; pour la réciproque, on note que et satisfont l'équation homogène , et ils ne peuvent pas être nuls simultanément.

La règle de suppression de Roth

Étant donné deux matrices carrées complexes et de taille et , et une matrice de taille , on peut se demander quand les deux matrices carrées suivantes

: et

de taille sont semblables. La réponse est que ces deux matrices sont semblables exactement lorsqu'il existe une matrice telle que , en d'autres termes, si est une solution d'une équation de Sylvester. Cette observation est appelée la règle de suppression de Roth [4].

On vérifie facilement une direction : Si alors

La règle de suppression de Roth ne se généralise pas aux opérateurs bornés de dimension infinie sur un espace de Banach[5].

Résolutions numériques

Un algorithme classique de résolution numérique de l'équation de Sylvester est l'algorithme de Bartels-Stewart, qui consiste à transformer et en forme de Schur par un algorithme QR, puis en résolvant le système triangulaire résultant par substitution. Cet algorithme, dont le coût de calcul est en opérations arithmétiques, est utilisé, entre autres, par le logiciel LAPACK et la fonction lyap dans GNU Octave[6]. Elle est proche de la fonction sylvester dans cet ensemble[7] - [8]. Dans certaines applications spécifiques de traitement d'image , la solution de l'équation de Sylvester dérivée a une forme close[9].

Articles connexes

Notes

  1. On voit aussi l'écriture équivalente .
  2. Bhatia et Rosenthal 1997.
  3. Cependant, cette réécriture de n'est pas conseillée car la solution numérique coûteuse en temps peut être mal conditionnée.
  4. F. Gerrish et A.G.B. Ward, « Sylvester's matrix equation and Roth's removal rule », The Mathematical Gazette, vol. 82, no 495, , p. 423–430 (DOI 10.2307/3619888, JSTOR 3619888).
  5. Bhatia et Rosenthal 1997, p. 3.
  6. « Function Reference: Lyap »
  7. « Functions of a Matrix (GNU Octave (version 4.4.1)) »
  8. La commande syl est obsolète depuis la version GNU Octave Version 4.0.
  9. Q. Wei, N. Dobigeon et J.-Y. Tourneret, « Fast Fusion of Multi-Band Images Based on Solving a Sylvester Equation », Transactions on Image Processing, vol. 24, no 11, , p. 4109–4121 (PMID 26208345, DOI 10.1109/TIP.2015.2458572, Bibcode 2015ITIP...24.4109W, arXiv 1502.03121, S2CID 665111).

Références

Liens externes

Cet article est issu de wikipedia. Text licence: CC BY-SA 4.0, Des conditions supplémentaires peuvent s’appliquer aux fichiers multimédias.