Matrice définie positive

En algèbre linéaire, une matrice définie positive est une matrice positive inversible.

Définitions

Notations

Soit $A$ une matrice à éléments réels ou complexes, par la suite on notera :

$A^{\mathsf {T}}$ la matrice transposée de $A$ ;
$A^{*}$ la matrice transconjuguée de $A$ (conjuguée de la transposée) ;
$\mathbb {R}$ le corps des nombres réels ;
$\mathbb {C}$ le corps des nombres complexes.

Matrice symétrique réelle définie positive

Soit $M$ une matrice symétrique réelle d'ordre $n$ . Elle est dite définie positive si elle est positive et inversible, autrement dit si elle vérifie l'une des quatre propriétés équivalentes suivantes :

Pour toute matrice colonne non nulle ${\textbf {x}}$ à $n$ éléments réels, on a : ${\textbf {x}}^{\mathsf {T}}M{\textbf {x}}>0$ . Autrement dit, la forme quadratique définie par $M$ est strictement positive pour $\mathbf {x} \neq 0$ .
Toutes les valeurs propres de $M$ (qui sont nécessairement réelles) sont strictement positives.
La forme bilinéaire symétrique $\mathbb {R} ^{n}\times \mathbb {R} ^{n}\to \mathbb {R} ,\quad ({\textbf {x}},{\textbf {y}})\mapsto {\textbf {x}}^{\mathsf {T}}M{\textbf {y}}$ est un produit scalaire sur $\mathbb {R} ^{n}$ .
Il existe une matrice $N\in {\mathcal {M}}_{n}(\mathbb {R} )$ inversible telle que $M=N^{\mathsf {T}}\,N$ (autrement dit : $M$ est congruente à la matrice identité).

Une matrice symétrique réelle est dite définie négative si son opposée (symétrique elle aussi) est définie positive.

La caractérisation 4 ci-dessus peut se justifier ainsi :

pour toute matrice carrée réelle $N$ , telle que son noyau se réduit au singleton 0, la matrice symétrique $N^{\mathsf {T}}\,N$ est positive ;
réciproquement, toute matrice réelle symétrique positive est de cette forme (la matrice $N$ n'est pas unique ; elle l'est si l'on impose qu'elle soit elle-même positive) ;
or si $M=N^{\mathsf {T}}\,N$ (avec $N$ carrée) alors $M$ est inversible si et seulement si $N$ l'est.

Elle permet de montrer que la matrice de Gram d'une famille de $n$ vecteurs d'un espace préhilbertien (réel ou complexe) est définie positive si et seulement si la famille est libre. Par exemple, toute matrice de Hilbert est définie positive.

Matrice hermitienne définie positive

On étend les propriétés et définitions précédentes aux matrices complexes.

Soit $M$ une matrice carrée complexe d'ordre $n$ . Elle est dite définie positive si elle vérifie l'une des quatre propriétés équivalentes suivantes :

Pour toute matrice colonne non nulle ${\textbf {z}}$ à $n$ éléments complexes, le nombre complexe ${\textbf {z}}^{*}M{\textbf {z}}$ est un réel strictement positif.
$M$ est hermitienne et toutes ses valeurs propres sont strictement positives.
La forme sesquilinéaire $\mathbb {C} ^{n}\times \mathbb {C} ^{n}\to \mathbb {C} ,\quad ({\textbf {x}},{\textbf {y}})\mapsto {\textbf {x}}^{*}M{\textbf {y}}$ est un produit scalaire sur $\mathbb {C} ^{n}$ (au sens : forme hermitienne définie positive).
Il existe une matrice $N\in {\mathcal {M}}_{n}(\mathbb {C} )$ inversible telle que $M=N^{*}\,N$ .

Une matrice $M$ est dite définie négative si son opposée est définie positive.

Intérêt des matrices définies positives

Beaucoup de problèmes de résolution de systèmes linéaires les plus faciles à traiter numériquement sont ceux dont les matrices sont symétriques définies positives^[1] : on dispose d'algorithmes numériquement stables et rapides pour l'inversion^[2] et la diagonalisation des matrices définies positives.

Toute matrice symétrique réelle positive est limite d'une suite de matrices symétriques réelles définies positives, ce qui est à la base de nombreux raisonnements par densité^[3].

Propriétés et critères

Propriétés

La matrice inverse d'une matrice définie positive est définie positive.
Si $M$ est définie positive et si $r$ est un réel strictement positif, alors $rM$ est définie positive.
Si $M$ et $N$ sont positives et si l'une des deux est inversible, alors $M+N$ est définie positive.
Une matrice positive est définie positive si et seulement si sa racine carrée positive est inversible. Cette propriété est utilisée pour la décomposition polaire (voir infra).
Inégalité de Hadamard : le déterminant d'une matrice définie positive est inférieur ou égal au produit de ses éléments diagonaux.

Critère de Sylvester

Pour qu'une matrice $A=\left(a_{ij}\right)_{1\leqslant i,j\leqslant n}$ , réelle symétrique ou complexe hermitienne, soit définie positive, il faut et suffit que les $n$ matrices $A_{p}=\left(a_{ij}\right)_{1\leqslant i,j\leqslant p}$ pour $p$ de 1 à $n$ , aient leur déterminant strictement positif, autrement dit que les $n$ mineurs principaux dominants soient strictement positifs.

Remarques

Pour $n=2$ , le critère de Sylvester est essentiellement le critère de positivité du trinôme du second degré.
En fait, sur un corps (commutatif) quelconque, cette condition de non-nullité des mineurs principaux est une condition nécessaire et suffisante pour qu'il existe une matrice $Q$ triangulaire supérieure telle que $Q^{\mathsf {T}}AQ$ soit diagonale et de rang maximum (il suffit d'adapter la démonstration qui suit).

Démonstration

Notons $q$ la forme quadratique associée à $A$ , définie par $q(\mathbf {x} )=\sum _{1\leqslant i,j\leqslant n}a_{ij}x_{i}x_{j}$ .

La condition est nécessaire. On remarque d'abord que si $q$ est définie positive, alors $\det(A)>0$ . En effet, par rapport à une base orthogonale pour cette forme quadratique (il en existe, d'après la réduction de Gauss), la matrice de $q$ s'écrit $\mathrm {diag} (c_{1},\cdots ,c_{n})$ les $c_{i}$ étant tous strictement positifs. Alors $c_{1}\cdots c_{n}=(\det A)(\det Q)^{2}$ ( $Q$ étant la matrice de passage), donc $\det A>0$ . Le résultat s'ensuit, en appliquant le même raisonnement à la restriction de $q$ aux sous-espaces $\mathbb {R} ^{k}\times \{0\}^{n-k}$ , pour $1\leqslant k\leqslant n-1$ .

Montrons maintenant que la condition est suffisante. On procède par récurrence sur la dimension.

Pour $n=0$ c'est évident puisqu'en dimension 0 l'ensemble des vecteurs non nuls est vide. Supposons la propriété vraie pour $n-1$ et notons $E=\mathbb {R} ^{n-1}\times \{0\}$ . Par hypothèse de récurrence, $q_{\vert E}$ est définie positive. De plus, $q$ est non dégénérée (parce que le déterminant de $A$ est non nul) donc :

$\,\mathbb {R} ^{n}=E\oplus E^{\perp }\quad {\textrm {avec}}\quad \dim E^{\perp }=1$

Soient $e$ un vecteur non nul de $E^{\perp }$ et $a=q(e)$ . Alors $\det A$ et $a\det A_{n-1}$ ont même signe d'après le même argument que dans la première partie (qui met implicitement en jeu le discriminant), or par hypothèse $\det A$ et $\det A_{n-1}$ sont strictement positifs. Donc $a>0$ , si bien que la restriction de $q$ à $E^{\perp }$ est, elle aussi, définie positive, ce qui montre que $q$ est définie positive.

Dans le cas complexe, plus général, la preuve est analogue, en considérant la forme hermitienne définie par la matrice.

Une autre méthode est d'utiliser le théorème d'entrelacement de Cauchy^[4].

Décomposition polaire

Article détaillé : Décomposition polaire.

Toute matrice inversible réelle se décompose de façon unique en produit d'une matrice orthogonale et d'une matrice symétrique définie positive.
Toute matrice inversible complexe se décompose de façon unique en produit d'une matrice unitaire et d'une matrice hermitienne définie positive.

La notion de matrice définie positive est donc analogue à celle de nombre réel strictement positif dans la forme polaire d'un nombre complexe.

Notes et références

↑ Philippe Ciarlet, Introduction à l'analyse numérique matricielle et à l'optimisation, Paris, Dunod, 1998, p. 26.
↑ Laboratoire bordelais de recherche en informatique (LaBRI), « Matrices symétriques définies positives et leur inversion » [PDF], sur labri.fr (consulté le 12 mars 2022).
↑ Jean Voedts, Cours de mathématiques : MP-MP*, Paris, Ellipses, 2002, 1119 p. (ISBN 2-7298-0666-0, OCLC 470114804, BNF 38853625, SUDOC 061031119), p. 634.
↑ (en) Roger A. Horn et Charles R. Johnson, Matrix Analysis, Cambridge University Press, 2013, 2^e éd. (1^re éd. 1985), 643 p. (ISBN 978-0-521-83940-2, lire en ligne), p. 439.

Articles connexes

Matrice à diagonale dominante
Factorisation de Cholesky - Pour toute matrice symétrique définie positive $A$ , il existe une matrice triangulaire inférieure $L$ telle que $A=LL^{\mathsf {T}}$ .
Complément de Schur