ATLAS Ouvrir dans l'atlas →

ExerciceMaths expertes, chapitre 9 — Chaînes de Markov · Pour aller plus loin

Exercice 22 Puissances d'une matrice de transition à deux états

DémonstrationCalcul

Soient pp et qq deux réels de ]0 ;1[]0\,;1[, et P=(1−ppq1−q)P=\begin{pmatrix} 1-p & p \\ q & 1-q \end{pmatrix}. On pose r=1−p−qr=1-p-q et : A=1p+q(qpqp)etB=1p+q(p−p−qq)A=\frac{1}{p+q}\begin{pmatrix} q & p \\ q & p \end{pmatrix} \qquad \text{et} \qquad B=\frac{1}{p+q}\begin{pmatrix} p & -p \\ -q & q \end{pmatrix}

1. Vérifier que A+B=I2A+B=I_{2}, P=A+rBP=A+rB, A2=AA^{2}=A, B2=BB^{2}=B et AB=BA=0AB=BA=0.

2. Montrer par récurrence que, pour tout n∈ℕn\in\mathbb{N}, Pn=A+rnBP^{n}=A+r^{n}B.

3. Justifier que −1<r<1-1<r<1, puis déterminer la limite de chaque coefficient de PnP^{n}. Que représentent les lignes de la matrice limite ?

4. Application : retrouver la distribution invariante de l'exemple du cours, P=(0,80,20,40,6)P=\begin{pmatrix} 0{,}8 & 0{,}2 \\ 0{,}4 & 0{,}6 \end{pmatrix}, et donner l'expression de PnP^{n}.

Source : inspiré de poly « Mathématiques : du lycée aux CPGE scientifiques » (LLG / H4), n° 7 et 11 (suites récurrentes linéaires), transposition matricielle

Corrigé

1. On calcule, en factorisant par 1p+q\frac{1}{p+q} ou 1(p+q)2\frac{1}{(p+q)^{2}} :

• A+B=1p+q(p+q00p+q)=I2A+B=\frac{1}{p+q}\begin{pmatrix} p+q & 0 \\ 0 & p+q \end{pmatrix}=I_{2} ;

• (qpqp)2=(q2+pqqp+p2q2+pqqp+p2)=(p+q)(qpqp)\begin{pmatrix} q & p \\ q & p \end{pmatrix}^{2}=\begin{pmatrix} q^{2}+pq & qp+p^{2} \\ q^{2}+pq & qp+p^{2} \end{pmatrix}=(p+q)\begin{pmatrix} q & p \\ q & p \end{pmatrix}, d'où A2=AA^{2}=A ; de même (p−p−qq)2=(p+q)(p−p−qq)\begin{pmatrix} p & -p \\ -q & q \end{pmatrix}^{2}=(p+q)\begin{pmatrix} p & -p \\ -q & q \end{pmatrix}, d'où B2=BB^{2}=B ;

• (qpqp)(p−p−qq)=(qp−pq−qp+pqqp−pq−qp+pq)=0\begin{pmatrix} q & p \\ q & p \end{pmatrix}\begin{pmatrix} p & -p \\ -q & q \end{pmatrix}=\begin{pmatrix} qp-pq & -qp+pq \\ qp-pq & -qp+pq \end{pmatrix}=0, et de même dans l'autre ordre : AB=BA=0AB=BA=0 ;

• le coefficient (1,1)(1,1) de A+rBA+rB vaut q+(1−p−q)pp+q=(p+q)−p(p+q)p+q=1−p\frac{q+(1-p-q)p}{p+q}=\frac{(p+q)-p(p+q)}{p+q}=1-p. De même pour les trois autres, on trouve P=A+rBP=A+rB.

2. Initialisation : P0=I2=A+B=A+r0BP^{0}=I_{2}=A+B=A+r^{0}B.

Hérédité : Si Pn=A+rnBP^{n}=A+r^{n}B, alors Pn+1=(A+rnB)(A+rB)=A2+rAB+rnBA+rn+1B2=A+rn+1BP^{n+1}=(A+r^{n}B)(A+rB)=A^{2}+rAB+r^{n}BA+r^{n+1}B^{2}=A+r^{n+1}B.

Conclusion : Pn=A+rnB\boxed{P^{n}=A+r^{n}B} pour tout n∈ℕn\in\mathbb{N}.

3. 0<p+q<20<p+q<2, donc −1<r=1−(p+q)<1-1<r=1-(p+q)<1, et rn→0r^{n}\to 0. Ainsi : Pn→n→+∞A=(qp+qpp+qqp+qpp+q)P^{n}\xrightarrow[n\to+\infty]{}A=\begin{pmatrix} \frac{q}{p+q} & \frac{p}{p+q} \\[1mm] \frac{q}{p+q} & \frac{p}{p+q} \end{pmatrix} Les deux lignes de la limite sont égales à la distribution invariante π=(qp+q,pp+q)\pi=\left(\frac{q}{p+q},\frac{p}{p+q}\right) trouvée dans la démonstration de la Proposition 4. Quel que soit l'état de départ, la loi de XnX_{n} tend vers π\pi.

4. Ici p=0,2p=0{,}2, q=0,4q=0{,}4 et r=0,4r=0{,}4 : Pn=(23132313)+0,4n(13−13−2323)\boxed{P^{n}=\begin{pmatrix} \frac{2}{3} & \frac{1}{3} \\[1mm] \frac{2}{3} & \frac{1}{3} \end{pmatrix}+0{,}4^{n}\begin{pmatrix} \frac{1}{3} & -\frac{1}{3} \\[1mm] -\frac{2}{3} & \frac{2}{3} \end{pmatrix}} On retrouve π=(2313)\pi=\begin{pmatrix} \frac{2}{3} & \frac{1}{3} \end{pmatrix}. Contrôle pour n=1n=1 : 23+0,4×13=2,43=0,8\frac{2}{3}+0{,}4\times\frac{1}{3}=\frac{2{,}4}{3}=0{,}8. Cette décomposition joue le rôle de la diagonalisation de l'exercice 26 du chapitre 3 : AA et BB « séparent » les deux comportements, la partie stable AA et la partie qui s'éteint, rnBr^{n}B.

Dans l'édifice

Thèmes

matrice de transition, distributions, distribution invariante, puissances de matrices, convergence, suites

Notions

matrices (espace des), mesure invariante

Ouvrir dans l'atlas interactif →