Décomposition d'une matrice en éléments propres

En algèbre linéaire, la décomposition d'une matrice en éléments propres est la factorisation de la matrice en une forme canonique où les coefficients matriciels sont obtenus à partir des valeurs propres et des vecteurs propres.

Aspects théoriques de la détermination des valeurs propres et vecteurs propres d'une matriceModifier

Un vecteur non nul v à N lignes est un vecteur propre d'une matrice carrée A à N lignes et N colonnes si et seulement si il existe un scalaire λ tel que :

 

λ est appelé valeur propre associée à v. Cette dernière équation est appelée « équation aux valeurs propres ».

Ces valeurs propres sont les solutions de l'équation :

 

On appelle p(λ) le polynôme caractéristique de A, et cette équation, l'équation caractéristique, est une équation polynomiale de degré N dont λ est l'inconnue. Cette équation admet Nλ solutions distinctes, avec 1 ≤ NλN. L'ensemble des solutions, i. e. des valeurs propres, est appelé le spectre de A.

On peut factoriser p :

 

avec

 

Pour chaque valeur propre λi, on a une équation particulière :

 

qui admet mi vecteurs solutions linéairement indépendants, formant une base de l'espace de toutes les solutions (le sous-espace propre associé à la valeur propre λi). Il est important de remarquer que cette multiplicité géométrique mi peut être égale ou pas à la multiplicité algébrique ni, mais qu'on a toujours : 1 ≤ mini. Le cas le plus simple est évidemment mi = ni = 1.

Le nombre de vecteurs propres indépendants de la matrice, noté ici Nv, est égal à la somme :   Les vecteurs propres peuvent alors être indexés par leurs valeurs propres respectives, avec un double indice : on appellera alors vi,j le j-ième vecteur propre associé à la i-ième valeur propre. Les vecteurs propres peuvent aussi être notés plus simplement, avec un seul indice : vk, avec k = 1, 2, ... , Nv.

Décomposition d'une matrice en éléments propresModifier

Soit A une matrice carrée (N lignes et N colonnes) admettant N vecteurs propres linéairement indépendants,   Alors, A peut s'écrire sous la forme :

 

Où la matrice de passage Q est une matrice carrée (à N lignes et N colonnes) dont la i-ième colonne est le vecteur propre   de A et Λ est la matrice diagonale dont les coefficients diagonaux sont les valeurs propres, i.e.,  .

Les vecteurs propres   sont souvent choisis unitaires, mais pas toujours.

Inversion d'une matrice via sa décomposition en éléments propresModifier

Si une matrice carrée A est diagonalisable et que toutes ses valeurs propres sont non nulles, alors A est inversible, et son inverse vaut :

 

Or, Λ étant diagonale, les coefficients de son inverse se calculent trivialement :

 

Conséquences sur le calcul des puissancesModifier

La décomposition en éléments simples permet de calculer facilement les fonctions polynomiales de matrices. Si f(x) s'écrit

 

alors, on sait que :

 

et Λ étant diagonale, un polynôme en Λ est très facile à calculer :

 

Les coefficients non diagonaux de f(Λ) sont nuls ; f(Λ) est donc également diagonale. Le calcul de f(A) revient donc à calculer l'image par f de chaque valeur propre.

Exemples
  •  
  •  

Une technique similaire s'applique plus généralement au calcul fonctionnel holomorphe, en utilisant la formule

 

ci-dessus. On trouve encore

 

Cas particuliers de décomposition en éléments simplesModifier

Matrices symétriques réellesModifier

Toute matrice symétrique réelle à N lignes et N colonnes admet N vecteurs propres linéairement indépendants. De plus, ces vecteurs peuvent être choisis de façon à être orthogonaux deux à deux et unitaires. Donc, toute matrice symétrique réelle A peut s'écrire sous la forme :

 

Q est une matrice orthogonale, et Λ est une matrice diagonale réelle.

Matrices normalesModifier

De la même façon, une matrice normale complexe admet une base orthonormale de vecteurs propres, et peut donc s'écrire sous la forme :

 

U est une matrice unitaire. De plus, si A est hermitienne, la matrice diagonale Λ a tous ses coefficients réels, et si A est unitaire, les coefficients diagonaux de Λ ont tous pour module 1.

Articles connexesModifier

(en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Eigendecomposition of a matrix » (voir la liste des auteurs).