En algèbre linéaire, une matrice définie positive est une matrice positive inversible.
Définitions
Notations
Soit une matrice à éléments réels ou complexes, par la suite on notera :
- la matrice transposée de ;
- la matrice transconjuguée de (conjuguée de la transposée) ;
- le corps des nombres réels ;
- le corps des nombres complexes.
Matrice symétrique réelle définie positive
Soit une matrice symétrique réelle d'ordre . 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 à éléments réels, on a : . Autrement dit, la forme quadratique définie par est strictement positive pour .
- Toutes les valeurs propres de (qui sont nécessairement réelles) sont strictement positives.
- La forme bilinéaire symétrique est un produit scalaire sur .
- Il existe une matrice inversible telle que (autrement dit : 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 , telle que son noyau se réduit au singleton 0, la matrice symétrique est positive ;
- réciproquement, toute matrice réelle symétrique positive est de cette forme (la matrice n'est pas unique ; elle l'est si l'on impose qu'elle soit elle-même positive) ;
- or si (avec carrée) alors est inversible si et seulement si l'est.
Elle permet de montrer que la matrice de Gram d'une famille de 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 une matrice carrée complexe d'ordre . 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 à éléments complexes, le nombre complexe est un réel strictement positif.
- est hermitienne et toutes ses valeurs propres sont strictement positives.
- La forme sesquilinéaireest un produit scalaire sur (au sens : forme hermitienne définie positive).
- Il existe une matrice inversible telle que .
Une matrice 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 est définie positive et si est un réel strictement positif, alors est définie positive.
- Si et sont positives et si l'une des deux est inversible, alors 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 , réelle symétrique ou complexe hermitienne, soit définie positive, il faut et suffit que les matrices pour de 1 à , aient leur déterminant strictement positif, autrement dit que les mineurs principaux dominants soient strictement positifs.
- Remarques
- Pour , 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 triangulaire supérieure telle que soit diagonale et de rang maximum (il suffit d'adapter la démonstration qui suit).
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
- 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, , p. 26.
- ↑ Laboratoire bordelais de recherche en informatique (LaBRI), « Matrices symétriques définies positives et leur inversion » [PDF], sur labri.fr (consulté le ).
- ↑ Jean Voedts, Cours de mathématiques : MP-MP*, Paris, Ellipses, , 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, , 2e éd. (1re é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 , il existe une matrice triangulaire inférieure telle que .
- Complément de Schur