Matrice définie positive - Définition

Source: Wikipédia sous licence CC-BY-SA 3.0.
La liste des auteurs de cet article est disponible ici.

Introduction

En algèbre linéaire, la notion de matrice définie positive est analogue à celle de nombre réel strictement positif : une matrice définie positive est une matrice positive inversible.

On introduit tout d'abord les notations suivantes ; si a est une matrice à éléments réels ou complexes :

  • {}^t \! a désigne la matrice transposée de a
  • a * désigne la matrice transconjuguée de a (conjuguée de la transposée)

On rappelle que :

  • \mathbb{R} désigne le corps des nombres réels
  • \mathbb{C} désigne 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 vérifie l'une des 3 propriétés équivalentes suivantes :

1. Pour toute matrice colonne non nulle \ \textbf{x} à n éléments réels, on a
{}^t \! \textbf{x} M \textbf{x} > 0.

(autrement dit, la forme quadratique définie par \,M est strictement positive pour \,\mathbf{x}\not=0 )

2. Toutes les valeurs propres de M sont strictement positives, c'est-à-dire :
\ \mathrm{sp}(M)  \subset\, ]0,\, +\infty[\, .
3. La forme bilinéaire symétrique définie par la relation
\langle \textbf{x},\textbf{y}\rangle_M = {}^t \! \textbf{x} M \textbf{y}

est un produit scalaire sur \mathbb{R}^n (identifié ici à l'espace vectoriel des matrices colonnes à n éléments réels).

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 propriété 1 signifie que M définit sur \mathbb{R}^n une forme quadratique définie positive, la propriété 2 que sur \mathbb{R}^n , vu comme espace euclidien avec le produit scalaire \langle x,y \rangle =\sum_{i=1}^nx_iy_i , M définit un opérateur auto-adjoint positif. L'équivalence entre 1 et 2 vient de cette double interprétation, à la lumière de la réduction de Gauss et du théorème spectral. Si 1 est vraie, sachant que les valeurs propres d'une matrice symétrique réelle sont réelles, on voit en appliquant 1 aux vecteurs propres que les valeurs propres sont strictement positives. Si 2 est vraie, il existe une matrice orthogonale \,Q telle que \,QAQ^{-1} soit diagonale (parce que \,A est symétrique réelle) à coefficients diagonaux strictement positifs (c'est l'hypothèse 2 sur les valeurs propres). Mais comme Q − 1 = tQ, la matrice \,A est aussi congrue à la matrice diagonale en question, donc la forme quadratique {}^t \! \textbf{x} M \textbf{x} est définie positive.

Exemple de base

Pour toute matrice réelle \,A , les matrices symétriques \,{}^tAA et \,A{}^tA sont positives ; elles sont définies positives si et seulement si \,A est inversible. Les matrices de Gram donnent un exemple de cette situation.

Plus précisément, c'est un exemple générique, puisque:

Une matrice M \in M_n( \R ) est définie positive si et seulement si il existe une matrice A \in M_n( \R ) inversible telle que M ={}^{\operatorname t}AA , c'est-à-dire si et seulement si elle est congruente à la matrice identité.

La matrice A n'est pas unique. Elle l'est si on impose qu'elle soit elle-même définie positive.

Si M ={}^{\operatorname t}AA , alors \forall x \in \R^n, {}^{\operatorname t}xMx = {}^{\operatorname t}(Ax)(Ax)  = \| Ax \|^2 \ge 0 , et si ce terme est nul, alors Ax=0~ , et si l'on suppose A inversible, alors x est nul.

Inversement, si M est définie positive, elle est diagonalisable avec une matrice de passage P orthogonale (puisque symétrique réelle), la matrice D={}^{\operatorname t}PMP ayant des valeurs propres λi strictement positives. Il suffit de définir la matrice \Delta~ comme étant la matrice diagonale dont les termes diagonaux sont les \sqrt \lambda_i , et de poser A =\Delta {}^{\operatorname t}P ~ , car alors {}^{\operatorname t}AA=M . Si l'on veut une matrice définie positive, il suffit de poser plutôt .

Exemple : matrice de Hilbert

On appelle matrice de Hilbert la matrice (symétrique d'ordre n) H = (h_{i,\, j}) , telle que h_{i,\, j} = \frac{1}{i + j  - 1} . Elle est définie positive.

En effet, soit une matrice colonne quelconque \ \textbf{x} à n éléments réels x_1, \dots, x_n .
On remarque que \forall\, i,\, \forall\, j,\, h_{i,\,j} = \int_0^1 t^{i + j - 2}\, dt . Alors, par linéarité de l'intégrale :
{}^t \! \textbf{x} H \textbf{x} = \sum_{i = 1}^n \sum_{j = 1}^n h_{i,\, j}\, x_i\, x_j = \int_0^1 \sum_{i = 1}^n \sum_{j = 1}^n t^{i + j - 2}\, x_i\, x_j\, dt =
\int_0^1 \sum_{i = 1}^n \sum_{j = 1}^n \left(t^{i - 1}\, x_i\right)\,  \left(t^{j - 1}\,  x_j\right)\,dt  = \int_0^1 \left(\sum_{i = 1}^n t^{i - 1}\, x_i\right)\, \left(\sum_{j = 1}^n t^{j - 1}\,  x_j\right)\,dt ,
d'où enfin : {}^t \! \textbf{x} H \textbf{x} = \int_0^1 \left(\sum_{i = 1}^nx_i\,  t^{i - 1}\right)^2\, dt .
Dans cette dernière intégrale, l'intégrande est continu et à valeurs positives. Par conséquent :
  • {}^t \!\textbf{x} H \textbf{x} \geq 0  ;
  • si {}^t \! \textbf{x} H \textbf{x} = 0 , alors pour tout t \in\, [0,\, 1],\, \left(\sum_{i = 1}^nx_i\,  t^{i - 1}\right)^2 = 0 .
Donc pour tout t \in\, [0,\, 1],\, \sum_{i = 1}^nx_i\,  t^{i - 1} = 0 .
Il en résulte que les \ x_i , coefficients d'un polynôme admettant une infinité de racines, sont tous nuls, c'est-à-dire \ \textbf{x} = 0 .
Ceci prouve que {}^t \! \textbf{x} H \textbf{x} > 0 pour toute matrice colonne non nulle \ \textbf{x} à n éléments réels.

Nota : ceci est un cas particulier d'une propriété des matrices de Gram. 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.

Intérêt des matrices définies positives

  • Les 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.
  • 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é.
Page générée en 0.115 seconde(s) - site hébergé chez Contabo
Ce site fait l'objet d'une déclaration à la CNIL sous le numéro de dossier 1037632
A propos - Informations légales
Version anglaise | Version allemande | Version espagnole | Version portugaise