Maths Experte Graphes & Matrices Markovtransitionétatsprobabilitésétat stationnaire

Chaînes de Markov

Matrices de transition, états, évolution d'un système et état stationnaire.

Définition

Une chaîne de Markov est une suite de variables aléatoires X0,X1,X_0, X_1, \ldots prenant des valeurs dans un ensemble fini d’états {1,2,,n}\{1, 2, \ldots, n\}, avec la propriété de Markov :

P(Xn+1=jXn=i)=pij(ne deˊpend pas du passeˊ)P(X_{n+1} = j \mid X_n = i) = p_{ij} \quad \text{(ne dépend pas du passé)}


Matrice de transition

La matrice de transition T=(pij)T = (p_{ij}) contient les probabilités de passage d’un état à l’autre.

Propriétés :

  • pij0p_{ij} \geq 0
  • jpij=1\displaystyle\sum_j p_{ij} = 1 pour chaque ligne ii (matrice stochastique)

Distribution de probabilité

L’état à l’instant nn est représenté par un vecteur ligne de probabilités :

πn=(πn(1),,πn(k))avec jπn(j)=1\vec{\pi}_n = (\pi_n(1), \ldots, \pi_n(k)) \quad \text{avec } \sum_j \pi_n(j) = 1

Évolution :

πn+1=πnT\vec{\pi}_{n+1} = \vec{\pi}_n \cdot T

πn=π0Tn\vec{\pi}_n = \vec{\pi}_0 \cdot T^n


Exemple

Deux états : Soleil (S) et Pluie (P).

T=(0,70,30,40,6)T = \begin{pmatrix} 0{,}7 & 0{,}3 \\ 0{,}4 & 0{,}6 \end{pmatrix}

Si π0=(1,0)\vec{\pi}_0 = (1, 0) (soleil aujourd’hui) :

π1=(0,7,  0,3)\vec{\pi}_1 = (0{,}7,\; 0{,}3) π2=(0,7×0,7+0,3×0,4,  )=(0,61,  0,39)\vec{\pi}_2 = (0{,}7 \times 0{,}7 + 0{,}3 \times 0{,}4,\; \ldots) = (0{,}61,\; 0{,}39)


État stationnaire

L’état stationnaire π\vec{\pi}^* vérifie :

π=πTetjπ(j)=1\vec{\pi}^* = \vec{\pi}^* \cdot T \quad \text{et} \quad \sum_j \pi^*(j) = 1

Méthode : Résoudre le système linéaire.

Exemple : π(S)=0,7π(S)+0,4π(P)\vec{\pi}^*(S) = 0{,}7\,\pi^*(S) + 0{,}4\,\pi^*(P) avec π(S)+π(P)=1\pi^*(S) + \pi^*(P) = 1

0,3π(S)=0,4π(P)    π(S)=470{,}3\,\pi^*(S) = 0{,}4\,\pi^*(P) \implies \pi^*(S) = \dfrac{4}{7}, π(P)=37\pi^*(P) = \dfrac{3}{7}


Convergence

Pour une chaîne de Markov irréductible (tous les états communicants) et apériodique, la distribution converge vers l’unique état stationnaire quel que soit l’état initial :

πnn+π\vec{\pi}_n \xrightarrow[n\to+\infty]{} \vec{\pi}^*