P-matrice

En mathématiques, une P-matrice ou matrice P est une matrice carrée réelle dont les mineurs principaux sont strictement positifs. Certains auteurs qualifient ces matrices de totalement strictement positives.

Les P-matrices interviennent dans l'étude des problèmes de complémentarité linéaire.

Une notion voisine est celle de [[P0-matrice|PModèle:Ind-matrice]].

Notations

On note

• $[[1,n]]:=\{1,\dots ,n\}$ l'ensemble des $n$ premiers entiers,
• $u\cdot v$ le produit de Hadamard de deux vecteurs $u$ et $v$, qui est défini par $(u\cdot v{)}_i=u_i{v}_i,$ pour tout indice $i$,
• $M_{IJ}$ la sous-matrice de $M$ formée de ses éléments avec indices de ligne dans $I$ et indices de colonne dans $J$.

Définitions

La notion de P-matrice peut se définir de différentes manières, bien sûr équivalentes.

Le nom de ces matrices a été proposé par Fiedler et Pták (1966). L'équivalence entre les définitions 1 et 2 est due à Fiedler et Pták (1962).

Propriétés immédiates

De la définition 1, on déduit que

• $M\in 𝐏$ si, et seulement si, $M^{\mathrm{\top }}\in 𝐏$,
• si $M$ est symétrique, alors $M\in 𝐏$ si, et seulement si, $M$ est définie positive,
• $𝐏\cap {ℝ}^{n\times n}$ est un ouvert de ${ℝ}^{n\times n}$.

De la définition 2, on déduit que

• si $M+M^{\mathrm{\top }}$ est définie positive, alors $M\in 𝐏.$

Autres propriétés

Complémentarité linéaire

Un problème de complémentarité linéaire consiste à trouver un vecteur $x⩾0,$ tel que $Mx+q⩾0$ et $x^{\mathrm{\top }}(Mx+q)=0.$ Dans cette définition, $M\in {ℝ}^{n\times n},$ $q\in {ℝ}^n,$ $x^{\mathrm{\top }}$ est le transposé de $x$ et les inégalités doivent se comprendre composante par composante. Ce problème est parfois noté de manière compacte comme suit

$\text{CL}(M,q)0⩽x\perp (Mx+q)⩾0.$

Existence et unicité de solution

L'importance des P-matrices dans les problèmes de complémentarité linéaire vient du résultat d'existence et d'unicité suivant.

Dès lors, si $M\notin 𝐏$, il existe un vecteur $q\in {ℝ}^n$ tel que l'une des deux situations exclusives suivantes a lieu :

• soit $\text{CL}(M,q)$ n'a pas de solution,
• soit $\text{CL}(M,q)$ a plus d'une solution.

On ne peut cependant pas affirmer que, pour une matrice $M\notin 𝐏$, même symétrique et non dégénérée, il existe un vecteur $q$ tel que la première des deux situations ci-dessus ait lieu. Ainsi

$M=\left(\begin{array}{cc}1& 2\\ 2& 1\end{array}\right)$

n'est pas une P-matrice, mais le problème $\text{CL}(M,q)$ a une solution quel que soit $q.$

Caractérisation algorithmique

La complémentarité linéaire offre une autre caractérisation des P-matrices, en termes d'une propriété d'un algorithme de résolution de ces problèmes, l'algorithme de Newton-min. Celui-ci est bien défini lorsque la matrice $M$ est non dégénérée. Pour une telle matrice et un vecteur $q$ donnés, on peut associer à un ensemble d'indices $I\subset [[1,n]]$, un nœud noté $x^{(I)}$ qui est l'unique solution $x$ du système linéaire

$x_{I^c}=0\text{et}(Mx+q{)}_I=0.$

On a noté $I^c:=[[1,n]]\setminus I$ le complémentaire de $I$ dans $[[1,n]]$. Brièvement, l'algorithme de Newton-min est celui de Newton semi-lisse pour résoudre l'équation linéaire par morceaux

$\min (x,Mx+q)=0,$

qui est équivalente au problème $\text{CL}(M,q)$. Dans la version qui importe dans le résultat ci-dessous, l'algorithme de Newton-min détermine d'abord, au point courant $x\in {ℝ}^n$, l'ensemble d'indices

$I=\{i\in [[1,n]]:x_i>(Mx+q{)}_i\}$

et calcule ensuite l'itéré suivant $x^{+}=x^{(I)}$. On a la caractérisation suivante]).

Résolution en temps polynomial ?

On ne connait pas d'algorithme permettant de résoudre le problème $\text{CL}(M,q)$ en temps polynomial lorsque $M\in 𝐏$, mais certains ont proposé des arguments en faveur de cette possibilitéModèle:,Modèle:,.

Vérifier la P-matricité

Vérifier qu'une matrice donnée dans ${\mathbb{Q}}^{n\times n}$ est une P-matrice est un problème co-NP-complet.

Une manière évidente de vérifier la P-matricité d'une matrice $M$ donnée est de calculer ses $2^n-1$ mineurs principaux (test des mineurs principaux), ce qui requiert $O(n^3{2}^n)$ opérations. Rump (2003) a montré que, quel que soit $I\subset [[1,n]]$ non vide, on peut trouver une matrice $M\in {ℝ}^{n\times n}$ telle que $\det M_{II}<0$ et $\det M_{JJ}>0$ pour tout $J\subset [[1,n]]$ non vide et différent de $I$, si bien que le test des mineurs principaux ne peut négliger aucun mineur.

Tsatsomeros et Li (2000) ont proposé un test, fondé sur le complément de Schur, qui réduit le nombre d'opérations à $72^n$. Le test requiert toujours ce nombre exponentiel d'opérations si la matrice est une P-matrice, mais peut en demander beaucoup moins si $M\notin 𝐏$, car un seul mineur négatif suffit pour montrer cette non-appartenance.

Rump (2003) a proposé le premier test qui ne demande pas nécessairement un nombre exponentiel d'opérations pour vérifier la P-matricité.

Exemples

1. Une matrice de Cauchy $C\in {ℝ}^{n\times n}$ est une matrice réelle carrée dont l'élément $(i,j)$ est donné par

Modèle:Saut$C_{ij}=\frac{1}{a_i+b_j},$Modèle:Saut

où $a_i+b_j\ne 0$. Une matrice de Cauchy est une P-matrice si $a_1+b_1>0$ et si les suites $\{a_i\}$ et $\{b_j\}$ sont strictement croissantes. En particulier, une matrice de Hilbert est une P-matrice (c'est une matrice de Cauchy avec $a_i=b_i=i-1/2$ pour tout $i\in [[1,n]]$).

1. Considérons la matrice circulante $M\in {ℝ}^{n\times n},$ $n⩾2,$ définie parModèle:Saut$M=\left(\begin{array}{cccc}1& & & \alpha \\ \alpha & \ddots & & \\ & \ddots & 1& \\ & & \alpha & 1\end{array}\right)$Modèle:Sautou de manière plus précise par $M_{ij}=1$ si $i=j$, $M_{ij}=\alpha$ si $i=(jmodn)+1$ et $M_{ij}=0$ sinon. Alors
   • si $n$ est pair, alors $M\in 𝐏$ si, et seulement si, $|\alpha |<1$,
   • si $n$ est impair, alors $M\in 𝐏$ si, et seulement si, $\alpha >-1$.
2. Considérons la matrice circulante $M\in {ℝ}^{n\times n},$ $n⩾3,$ définie parModèle:Saut$M=\left(\begin{array}{ccccc}1& & & \beta & \alpha \\ \alpha & \ddots & & & \beta \\ \beta & \ddots & 1& \\ & \ddots & \alpha & 1\\ & & \beta & \alpha & 1\end{array}\right)$Modèle:Sautou de manière plus précise parModèle:Saut$M_{ij}=\{\begin{array}{cc}1& \text{si}i=j\\ \alpha & \text{si}i=(jmodn)+1\\ \beta & \text{si}i=((j+1)modn)+1\\ 0& \text{sinon}.\end{array}$Modèle:SautAlors $M\in 𝐏$ si $|\alpha |-1<\beta <|\alpha |/4$ ou si ${\alpha }^2+8(\beta -\frac{1}{2}{)}^2<2$.

Annexes

Articles connexes

• Complémentarité linéaire
• P0-matrice

Bibliographie

• Modèle:Ouvrage.
• Modèle:Ouvrage.
• Modèle:Ouvrage.
• Modèle:Article.