Calculer la probabilité à long terme
Énoncé
Soit une chaîne de Markov avec P = pmatrix 0{,}4 & 0{,}6 0{,}2 & 0{,}8 pmatrix et ^{(0)} = (1, 0).
Calculer ^{(10)} en utilisant la distribution stationnaire.
Indice : Trouve d'abord la distribution stationnaire . Pour un grand n, ^{(n)} (distribution stationnaire).
Correction
- Étape 1 : Calculer P^{10} à la main est hors de question. L'idée est de contourner : quand n grandit, ^{(n)} tend vers la distribution stationnaire, quelle que soit la distribution de départ. On calcule donc , puis on regarde si n = 10 est assez grand.
- Étape 2 : On résout = P. La première composante donne :
_1 = 0{,}4\,_1 + 0{,}2\,_2 0{,}6\,_1 = 0{,}2\,_2 _2 = 3\,_1
- Étape 3 : Avec la condition de somme :
_1 + 3_1 = 1 _1 = 1{4} _2 = 3{4}
- Étape 4 : Mais « n est grand » n'est pas une justification — il faut savoir À QUELLE VITESSE la convergence a lieu. Pour une chaîne à deux états, l'écart à l'équilibre est multiplié à chaque étape par un facteur qui se lit directement sur la matrice :
= 1 - p_{1,2} - p_{2,1} = 1 - 0{,}6 - 0{,}2 = 0{,}2
- Étape 5 : L'écart initial vaut 1 - 1{4} = 3{4}, puisqu'on part de ^{(0)} = (1\,;\,0). Après n étapes il ne reste donc que :
^{(n)}_1 = 1{4} + 3{4} (0{,}2)^n
- Étape 6 : Application à n = 10. Comme (0{,}2)^{10} 1{,}02 10^{-7} :
^{(10)}_1 = 0{,}25 + 0{,}75 1{,}02 10^{-7} 0{,}250\,000\,08
- Étape 7 : L'approximation est donc excellente — l'écart vaut 8 cent-millionièmes :
^{(10)} (1{4}\ ;\ 3{4}) = (0{,}25\ ;\ 0{,}75)
- Étape 8 : La convergence est très rapide, et le tableau le montre : n = 1 donne 0{,}400, n = 2 donne 0{,}280, n = 3 donne 0{,}256, n = 5 donne 0{,}250\,24. Dès la cinquième étape, on est à trois décimales de l'équilibre — parce que 0{,}2 est petit.
« n est grand » ne se dit pas sans mesurer :
Si valait 0{,}99, il faudrait des centaines d'étapes pour approcher l'équilibre, et ^{(10)} serait encore très loin de .
C'est , et lui seul, qui dit si n = 10 suffit. Répondre « ^{(10)} car 10 est grand » n'est pas une justification.
Le point le plus important du chapitre :
La limite ne dépend PAS de ^{(0)}. En partant de (0\,;\,1) ou de (0{,}5\,;\,0{,}5), on aboutirait au même (1{4}\,;\,3{4}).
La chaîne oublie son point de départ, et mesure la vitesse de cet oubli.