Décomposition de Schur

décomposition algébrique de la forme M=UAU*

En algèbre linéaire, une décomposition de Schur (nommée après le mathématicien Issai Schur) d'une matrice carrée complexe M est une décomposition de la forme

M = UAU*

où U est une matrice unitaire (U*U = I) et A une matrice triangulaire supérieure.

Démonstration

modifier

On peut écrire la décomposition de Schur en termes d'applications linéaires :

Soient E un ℂ-espace vectoriel de dimension n muni d'un produit scalaire et v un endomorphisme sur E, alors ∃(f1,...,fn) une base orthonormée de E et une famille (φik)1≤i,k≤n∈ℂ telle que ∀i,v(fi)=∑k=inφikfk.


Dans le cas où v est l'application nulle, l'énoncé est directement vérifié, on peut donc se contenter de traiter le cas où v est différente de l'application nulle. On démontre par récurrence forte sur la dimension n≥1 de E le résultat énoncé. L'initialisation est triviale, pour l'hérédité on considère deux cas différents :

Si la matrice Mv associée à v dans une base quelconque est diagonalisable, alors on peut choisir fn+1 un vecteur propre normé de Mv. On pose F=Vect(fn+1)⊥ et on considère l'application linéaire v~=πF∘v|F où πF est le projecteur orthogonal sur F et v|F la restriction de v à F. Comme v~ est un endomorphisme de l'espace vectoriel F qui est de dimension n, l'hypothèse de récurrence assure l'existence d'une base orthonormée (f1,...,fn) de F dans laquelle la matrice Mv~ associée à v~ est triangulaire supérieure. Il est alors clair que (f1,...,fn,fn+1) forme une base orthonormée dans laquelle la matrice Mv associée à v est triangulaire supérieure.

Si la matrice Mv associée à v dans une base quelconque n'est pas diagonalisable on a l'inégalité 0<rang(v)<n+1. On pose alors F=ker(v)C et on considère l'application linéaire v~=πF∘v|F où πF est le projecteur orthogonal sur F et v|F la restriction de v à F. Comme 1≤dim(F)=rang(v)<n+1 on peut utiliser l'hypothèse de récurrence qui assure l'existence d'une base orthonormée (f1,...,fdim(F)) de F dans laquelle la matrice Mv~ associée à v~ est triangulaire supérieure. On complète cette base en une base orthonormée (f1,...,fdim(F),fdim(F)+1,...,fn+1) de E. Comme Vect(fdim(F)+1,...,fn+1)=ker(v), il est clair que la matrice Mv associée à v est triangulaire supérieure dans cette même base. Cela termine la récurrence.

Propriétés

modifier

Il existe une telle décomposition (non unique en général) pour toute matrice carrée complexe M[1],[2].

A étant semblable à M, elle a les mêmes valeurs propres. Et A étant triangulaire, les valeurs propres se trouvent sur sa diagonale.

Puisque A = U*MU, si M est normale (M*M = MM*) alors A aussi donc (comme elle est de plus triangulaire) elle est diagonale[3]. En particulier, si M est hermitienne (M* = M) alors A est diagonale réelle.

Cas réel

modifier

Si une matrice réelle M est trigonalisable, elle possède une décomposition de la même forme avec de plus U et A réelles, autrement dit

M = PA tP

avec P orthogonale et A réelle et triangulaire supérieure.

Si M n'est pas trigonalisable, elle a « presque » une décomposition de cette forme (avec P orthogonale et A réelle) mais où A est seulement triangulaire par blocs, avec des blocs diagonaux de polynôme caractéristique irréductible, donc d’ordre 1 ou 2.

Notes et références

modifier
(en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Schur decomposition » (voir la liste des auteurs).
  1. ↑ (en) Roger A. Horn et Charles R. Johnson, Matrix Analysis, Cambridge University Press, , 561 p. (ISBN 978-0-521-38632-6, lire en ligne), p. 79 et s.
  2. ↑ (en) Gene H. Golub et Charles F. Van Loan, Matrix Computations, Johns Hopkins University Press, , 3e éd., 694 p. (ISBN 978-0-8018-5414-9, lire en ligne), p. 313.
  3. ↑ Golub et Van Loan 1996, p. 314.