Matrice stochastique
# Matrice stochastique

En mathématiques, une matrice stochastique (aussi appelée matrice de Markov) est une matrice carrée (finie ou infinie) dont chaque élément est un réel positif et dont la somme des éléments de chaque ligne vaut 1. Cela correspond, en théorie des probabilités, à la matrice de transition d'une chaîne de Markov.

## Définitions

Une matrice $M\in {\mathcal {M}}_{n}(\mathbb {R} )$ est dite stochastique si toutes ses entrées sont positives (ou nulles) et si, pour tout $i=1,...,n$, on a $\sum _{j=1}^{n}m_{ij}=1$, c'est-à-dire que la somme des coordonnées de chaque ligne vaut 1.

Une matrice stochastique est dite régulière s'il existe un entier $k>0$ tel que la matrice $M^{k}$ ne contient que des réels strictement positifs.

Une matrice est dite bistochastique (ou doublement stochastique) si la somme des éléments de chaque ligne et de chaque colonne vaut 1, autrement si $M$ et sa transposée $M^{t}$ sont stochastiques.

## Propriétés

Une autre caractérisation des matrices stochastiques est donnée par :

* $M$ est une matrice stochastique si et seulement si $M\geq 0$ (ses coefficients sont positifs ou nuls) et $Me=e$, où $e$ désigne le vecteur de $\mathbb {R} ^{n}$ dont toutes les coordonnées valent 1.
* $M$ est bistochastique si $M\geq 0$, $Me=e$ et $e^{*}M=e^{*}$, où $e^{*}$ est le vecteur transposé de $e$.

D'après la propriété précédente, puisque 1 est une valeur propre de $M$ avec comme vecteur propre à droite le vecteur colonne dont toutes les coordonnées valent 1 :

* Si $M$ est une matrice stochastique, on appelle vecteur stable pour $M$ un vecteur ligne non nul $v$ tel que $vM=v$ , autrement dit : un vecteur propre à gauche pour la valeur propre 1 (et $M$ possède toujours au moins un vecteur stable).

Une caractérisation du rayon spectral d'une matrice stochastique est donnée par :

* Si $M$ est une matrice stochastique, alors $||Mx||_{\infty }\leq ||x||_{\infty }$ pour tout $x\in \mathbb {C} ^{n}$, de sorte que le rayon spectral $\rho (M)\leq 1$. Or, comme $Me=e$, on a en fait $\rho (M)=1$. Ainsi, le rayon spectral d'une matrice stochastique vaut précisément 1.

D'autres résultats sont donnés par :

* Le produit de deux matrices stochastiques est stochastique.
* Toute matrice stochastique indexée par E×E opère sur l'espace des applications bornées de E dans $\mathbb {R}$ .
* Si $M$ est une matrice stochastique et si $p$ est une probabilité alors $pM$ est une probabilité.

## Exemple

La matrice suivante est stochastique mais pas bistochastique :

    $M={\begin{pmatrix}0{,}5&0{,}3&0{,}2\\0{,}2&0{,}8&0\\0{,}3&0{,}3&0{,}4\end{pmatrix}}.$

Le vecteur ${\begin{pmatrix}3&6&1\end{pmatrix}}$ est stable pour M.

La matrice stochastique M est régulière car

    $M^{2}={\begin{pmatrix}0{,}37&0{,}45&0{,}18\\0{,}26&0{,}70&0{,}04\\0{,}33&0{,}45&0{,}22\end{pmatrix}}.$

## Théorème des matrices stochastiques

Le théorème des matrices stochastiques stipule que, si A est une matrice stochastique régulière, alors

* la chaîne de Markov de matrice de transition A est irréductible ;
* A possède un unique vecteur stable t dont la somme des coordonnées vaut 1 (autrement dit : le sous-espace vectoriel des vecteurs stables est de dimension 1 — les vecteurs stables sont tous colinéaires) ;
* les coordonnées de t sont toutes strictement positives.

De plus, si $x_0$ est une loi initiale quelconque (i.e. est un vecteur à coordonnées positives ou nulles et de somme 1), et si $x_{k+1}$ = $x_kA$ pour k = 0, 1, 2, … alors la chaîne de Markov ${x_k}$ converge vers t quand $k\to \infty$. C’est-à-dire :

$\lim _{k\to \infty }\mathbf {x} _{0}A^{k}={\textbf {t}}.$

## Quelques autres résultats

Le rôle des matrices stochastiques est important, notamment dans l'étude des chaînes de Markov. Une caractéristique importante des matrices doublement stochastiques (ou bistochastiques) est fourni par les matrices de permutation $P(\sigma )$, $\sigma \in {\mathcal {S}}^{n}$, dont les coefficients valent $p_{ij}=\delta _{\sigma (i)}^{j}$, avec $\delta$ le symbole de Kronecker.

Le théorème de Birkhoff montre ce rôle central qu'ont les matrices de permutations dans la caractérisation des matrices bistochastiques :

Théorème de Birkhoff — Une matrice $M\in {\mathcal {M}}_{n}(\mathbb {R} )$ est doublement stochastique si et seulement si elle est barycentre de matrices de permutations.

Une conséquence de théorème est donnée par le résultat suivant :

Corollaire — Soit $||.||$ une norme sur $\mathbb {R} ^{n}$, invariante par permutation des coordonnées. Alors $||M||=1$ pour toute matrice doublement stochastique.

Deux autres résultats sur les matrices bistochastiques utilisent la relation décrite par le symbole $\prec$, défini par : Soient $a=(a_{1},...,a_{n})$ et $b=(b_{1},...,b_{n})$ deux suites de $n$ nombres réels. On dit que b majore a, et on note $a\prec b$ si :

* $a_{1}+...+a_{k}\leq b_{1}+...+b_{k}$ pour tout $k=1,...,n-1$ ;
* $a_{1}+...+a_{n}=b_{1}+...+b_{n}$.

Il s'agit d'une relation d'ordre partielle.

Les deux théorèmes sont :

Théorème — Une matrice $M$ est doublement stochastique si et seulement si $x\prec Mx$ pour tout $x\in \mathbb {R} ^{n}$.

Théorème — Soient $x,y\in \mathbb {R} ^{n}$. ALors $x\prec y$ si et seulement s'il existe une matrice $M$, doublement stochastique, telle que $y=Mx$.

