Polynôme caractéristique

En mathématiques, et plus particulièrement en algèbre linéaire, à toute matrice carrée à coefficients dans un anneau commutatif ou à tout endomorphisme d'un espace vectoriel de dimension finie est associé un polynôme appelé polynôme caractéristique. Il renferme d'importantes informations sur la matrice ou sur l'endomorphisme, comme ses valeurs propres, son déterminant et sa trace. Le théorème de Cayley-Hamilton assure que toute matrice carrée annule son polynôme caractéristique.

Motivation

Étant donné une matrice carrée Modèle:Mvar d'ordre n, on cherche un polynôme dont les racines sont précisément les valeurs propres de Modèle:Mvar.

Si Modèle:Mvar est une matrice diagonale ou plus généralement une matrice triangulaire, alors les valeurs propres de Modèle:Mvar sont les coefficients diagonaux λ₁, …, λ_n de Modèle:Mvar et nous pouvons définir le polynôme caractéristique comme étant

<math>(X-\lambda_1)(X-\lambda_2)\ldots(X-\lambda_n)\qquad (1)</math>

Ce polynôme est le déterminant Modèle:Math où Modèle:Mvar est la matrice identité.

Pour une matrice quelconque Modèle:Mvar, si λ est une valeur propre de Modèle:Mvar, alors il existe un vecteur colonne propre Modèle:Mvar non nul tel que MV = λV, soit (λI_n – M)V = 0 (où Modèle:Mvar est la matrice unité.) Puisque V est non nul, cela implique que la matrice λI_n – M est singulière, et donc a son déterminant nul. Cela montre que les valeurs propres de Modèle:Mvar sont des zéros de la fonction Modèle:Math ou des racines du polynôme <math>\det(XI_n-M)</math>.

Définition formelle

Soit Modèle:Mvar une matrice carrée d'ordre n à coefficients dans un anneau commutatif. Le polynôme caractéristique de Modèle:Mvar, noté p_M(X), est<ref>Modèle:Ouvrage.</ref> le polynôme défini par

<math>p_M(X):=\det(XI_n-M) = \sum_{\sigma\in\Sigma_n} \varepsilon(\sigma) a_{1\sigma(1)}\dots a_{n\sigma(n)} \qquad (2)</math>

où Modèle:Math est le déterminant des matrices, Modèle:Mvar désigne la matrice identité d'ordre n, a_{i j} = le polynôme δ_{i j} X − m_{i j}, coefficient d'indice (i, j) de la matrice XI_n – M ; la somme de droite (prise sur l'ensemble des permutations des indices) donne une expression explicite du polynôme caractéristique.

Remarque. Au lieu de l'expression (2), certains auteurs définissent le polynôme caractéristique comme étant <math>\det(M-XI_n)</math>. Avec cette définition, on a l'équation <math>p_M(0)=\det(M)</math>. Ceci n'est pas le cas pour la définition (2) lorsque l'ordre Modèle:Mvar est impair et <math>\det(M)\ne 0</math>, puisque l'on a : <math>\det(M-XI_n) = (\!-\,1)^n ~\det(XI_n-M)</math>. La définition (2) présente l'« avantage » de rendre le polynôme caractéristique unitaire.

Coefficients

Le développement du polynôme caractéristique Modèle:Math d'une matrice carrée Modèle:Mvar d'ordre n est donné par

<math>\det(XI_n-M)=X^n-f_1(M)X^{n-1}+f_2(M)X^{n-2}-\dots+(-1)^nf_n(M)</math>

où Modèle:Math est une fonction polynomiale<ref>Modèle:Ouvrage.</ref> en les coefficients de la matrice Modèle:Mvar.

Un développement explicite du déterminant de Modèle:Math donne<ref>Modèle:Harvsp.</ref> :

<math>f_k(M)=\sum_{I\subset\{1,\dots,n\},\;\operatorname{card}(I)=k}\det \left(M_{I,I}\right)</math>

c'est-à-dire la somme des mineurs principaux d'ordre Modèle:Mvar. En particulier, le coefficient constant Modèle:Math est égal à (–1)ⁿ fois le déterminant de Modèle:Mvar, et le coefficient de X^n–1 est égal à l'opposé de la trace de M.

La propriété la plus importante des polynômes caractéristiques est que les valeurs propres de Modèle:Mvar sont exactement les racines du polynôme Modèle:Math. En notant <math>(\lambda_1,\dots,\lambda_n)</math> les racines de Modèle:Mvar prises avec multiplicité,

<math>f_k(M)=s_k(\lambda_1,\dots,\lambda_n)</math>

où s_k désigne le k-ième polynôme symétrique élémentaire. (Ici, les racines sont prises dans une extension finie L de K lorsque K n'est pas algébriquement clos ; ainsi, Modèle:Mvar est trigonalisable sur L. C'est ce qui permet, pour démontrer la formule ci-dessus, de se ramener au cas décrit dans le paragraphe Motivation).

Lorsque le corps de base K est de caractéristique nulle (par exemple, <math>K = \R</math> ou <math>K = \Complex</math>), grâce aux identités de Newton, les coefficients Modèle:Math s'expriment comme des fonctions polynomiales des sommes de Newton des valeurs propres :

<math>\sum_{i=1}^n \lambda_i^j=\operatorname{Tr}(M^j).</math>

Toute fonction <math>M\mapsto f(M)</math> polynomiale en les coefficients de la matrice Modèle:Mvar et invariante par similitude est une fonction polynomiale en les coefficients du polynôme caractéristique. Cette propriété est par exemple utilisée dans la définition et la classification des classes caractéristiques en géométrie différentielle, ce qui dépasse de loin le niveau de cet article.

Exemples

• Pour une matrice Modèle:Mvar d'ordre 2, le polynôme caractéristique s'exprime simplement commeModèle:Retraitmais peut aussi se calculer directement à partir de la définition.
  Déterminons par exemple le polynôme caractéristique Modèle:Math de la matriceModèle:RetraitC'est le déterminant de la matriceModèle:RetraitOn a doncModèle:RetraitOn en déduit que 1 est une valeur propre double de la matrice.
• Pour une matrice Modèle:Mvar d'ordre 3, le polynôme caractéristique s'exprime commeModèle:RetraitoùModèle:Retraitavec Modèle:Math l'élément en position Modèle:Math dans la matrice Modèle:Mvar.

Propriétés

• Le polynôme p_M(X) d'une matrice carrée Modèle:Mvar d'ordre n est unitaire (son coefficient dominant est égal à 1) et son degré est égal à n.
• La matrice Modèle:Mvar et sa transposée ont le même polynôme caractéristique.
• Le théorème de Cayley-Hamilton affirme qu'en remplaçant X par Modèle:Mvar dans p_M(X), on obtient la matrice nulle : p_M(M) = 0, autrement dit le polynôme caractéristique est un polynôme annulateur de Modèle:Mvar. Ceci est équivalent à l'affirmation que le polynôme minimal divise le polynôme caractéristique de Modèle:Mvar.
• Deux matrices semblables ont même polynôme caractéristique. Autrement dit, pour toute matrice inversible P, <math>\det(XI_n-P^{-1}MP)=\det(XI_n-M)</math>. La réciproque est fausse : par exemple, les deux matrices <math>\begin{pmatrix}0&1\\0&0\end{pmatrix}</math> et <math>\begin{pmatrix}0&0\\0&0\end{pmatrix}</math> ont même polynôme caractéristique mais ne sont pas semblables.
• Une matrice Modèle:Mvar est semblable à une matrice triangulaire si et seulement si son polynôme caractéristique peut être complètement décomposé en produit de facteurs de degré 1 à coefficients dans K. En fait, Modèle:Mvar est même semblable à une matrice de Jordan dans ce cas.
• Par trigonalisation comme dans le point précédent (ce qui est toujours possible quitte à remplacer K par une clôture algébrique), le polynôme caractéristique d'un polynôme d'endomorphisme P(u) est le polynôme unitaire dont les racines sont les images par P de celles du polynôme caractéristique de u (répétées en cas de multiplicité). Sur C, cette propriété des polynômes en u se démontre de même, plus généralement, pour les fonctions entières de u.
• Propriété de commutation<ref>Voir la page d'exercices corrigés du chapitre « Valeurs et vecteurs propres - Polynôme caractéristique » sur Wikiversité.</ref> : pour toutes matrices <math>A\in\mathrm M_{m,n}(K)</math> et <math>B\in\mathrm M_{n,m}(K)</math>, où <math>K</math> est un anneau unitaire :<math>X^np_{AB}(X)=X^mp_{BA}(X)</math>.

Matrice compagnon

Soit <math>p(X)=X^n-\sum_{k=0}^{n-1} a_{k}X^{k}</math> un polynôme à coefficients dans K. La matrice d'ordre n

<math>M=\begin{pmatrix}

0 & 1 & 0 & \ldots & \ldots & 0 \\ \vdots & \ddots & \ddots & \ddots & & \vdots \\ \vdots & & \ddots & \ddots & \ddots & \vdots \\ \vdots & & & \ddots & \ddots & 0 \\ 0 & \ldots & \ldots & \ldots & 0 & 1 \\ a_0 & a_1 & a_2 & \ldots &a_{n-2} & a_{n-1} \end{pmatrix} </math> qui admet p(X) comme polynôme caractéristique (et polynôme minimal), est appelée matrice compagnon du polynôme (ou selon certains ouvrages, sa transposée). Une des méthodes utilisées en calcul numérique pour calculer des valeurs approchées des racines d'un polynôme est d'en construire la matrice compagnon puis de calculer des valeurs approchées des valeurs propres de cette matrice à l'aide d'une méthode itérative.

Matrice triangulaire

Dans le cas d'une matrice triangulaire (supérieure) d'ordre <math>n</math>, matrice de la forme :

<math>T=\begin{pmatrix}

t_{1,1} & t_{1,2} & \ldots & \ldots & t_{1,n} \\ 0 & t_{2,2} & \ldots & \ldots & t_{2,n} \\ \vdots & \ddots & \ddots & & \vdots \\ \vdots & & \ddots & \ddots & \vdots \\ 0 & \ldots & \ldots & 0 & t_{n,n} \end{pmatrix} </math> le déterminant <math>p_T(X)=\det(XI_n-T)</math> qui exprime le polynôme caractéristique se factorise :

<math>

p_T(X)=(X-t_{1,1}) (X-t_{2,2}) \ldots (X-t_{n,n}) </math>

Le même raisonnement s'applique bien sûr au cas d'une matrice triangulaire inférieure. D'une façon générale, les valeurs propres d'une matrice triangulaire coïncident donc effectivement avec ses éléments diagonaux, comme annoncé au début.

Généralisations et applications

Polynôme caractéristique d'un endomorphisme

Si f est un endomorphisme d'un espace vectoriel E de dimension finie, on appelle polynôme caractéristique de f le polynôme caractéristique de la matrice représentant f dans une base de E. Les polynômes caractéristiques de deux matrices semblables étant égaux, cette définition ne dépend pas de la base choisie ; on a <math> p_f(X)= \det(X{\rm Id}_E-f)</math>.

Équations caractéristiques

La recherche des valeurs propres d'une matrice (ou d'un endomorphisme) revient à déterminer les zéros de son polynôme caractéristique, et donc à résoudre l'équation <math>p_A(X)=0</math>. Lorsque cette matrice apparaît comme outil de résolution d'un problème, tel que la recherche des solutions d'une équation différentielle, ou d'une formule explicite pour une suite définie par récurrence, par exemple, on dit que l'équation précédente est l’équation caractéristique de ce problème. Ainsi, pour résoudre l'équation différentielle <math>y+by'+cy=0</math>, on construit le système différentiel <math>y'=z</math>, <math>z'=-cy-bz</math>, de matrice <math>M=\begin{pmatrix} 0 & 1\\ -c&-b \end{pmatrix}</math> ; le polynôme caractéristique de Modèle:Mvar est <math>X^2+bX+c</math>, et on retrouve bien l'équation caractéristique au sens d'Euler.

Polynôme caractéristique d'un graphe

On appelle polynôme caractéristique du graphe G le polynôme caractéristique de la matrice d'adjacence de G. L'étude de ce polynôme et de ses racines est l'objet de la théorie spectrale des graphes.

Voir aussi

Modèle:Autres projets Algorithme de Faddeev-Leverrier : algorithme permettant de calculer le polynôme caractéristique d'une matrice.