Aller au contenu

L'Établi des Maths Cours

Outils pour apprendre et pratiquer

Terminale · Mathématiques Expertes · TX5

Cours et exercices corrigés de terminale : Calcul matriciel, graphes et chaînes de Markov

Matrices (lignes, colonnes, carrées), addition, multiplication matricielle, matrice inverse, puissances de matrices M^n, graphes probabilistes et état stable d'une chaîne de Markov.

Objectifs du chapitre

Le cours

1. Produit matriciel et matrice inverse

  • Le produit matriciel A×BA \times B n'est possible que si le nombre de colonnes de AA égale le nombre de lignes de BB. En général, AB≠BAAB \neq BA (non-commutativité).
  • Inverse d'une matrice 2×22 \times 2 : pour A=(abcd)A = \begin{pmatrix} a & b \\ c & d \end{pmatrix}, si det⁡(A)=ad−bc≠0\det(A) = ad - bc \neq 0, alors :

A−1=1ad−bc(d−b−ca)A^{-1} = \frac{1}{ad - bc} \begin{pmatrix} d & -b \\ -c & a \end{pmatrix}

2. Chaînes de Markov et distribution stationnaire

  • Soit PP la matrice de transition d'un graphe probabiliste. Si πn\pi_n désigne la matrice ligne des probabilités des états au rang nn, alors πn+1=πnP\pi_{n+1} = \pi_n P et πn=π0Pn\pi_n = \pi_0 P^n.
  • Un état stable (ou distribution stationnaire) est une matrice ligne π=(xy)\pi = (x \quad y) telle que πP=π\pi P = \pi avec x+y=1x + y = 1.

3. Graphes et matrice d'adjacence

  • Un graphe est formé de sommets reliés par des arêtes (orientées ou non). L'ordre est le nombre de sommets ; le degré d'un sommet est le nombre d'arêtes dont il est une extrémité (pour un graphe non orienté).
  • Une chaîne (ou chemin) de longueur nn est une suite de nn arêtes consécutives ; un graphe est connexe si deux sommets quelconques sont reliés par une chaîne ; il est complet si tous ses sommets sont adjacents deux à deux.
  • La matrice d'adjacence AA d'un graphe d'ordre pp est la matrice carrée dont le coefficient aija_{ij} est le nombre d'arêtes allant du sommet ii au sommet jj.
  • Théorème : le coefficient (i,j)(i, j) de AnA^n est le nombre de chemins de longueur nn allant de ii à jj.

4. Chaînes de Markov

  • Une chaîne de Markov sur un ensemble fini d'états est une suite de variables aléatoires (Xn)(X_n) telle que la probabilité de passer de l'état ii à l'état jj, pij=P{Xn=i}(Xn+1=j)p_{ij} = P_{\{X_n = i\}}(X_{n+1} = j), ne dépend ni de nn ni des états antérieurs.
  • La matrice de transition P=(pij)P = (p_{ij}) est stochastique : coefficients positifs, somme de chaque ligne égale à 11.
  • Si πn\pi_n est la matrice ligne de la loi de XnX_n, alors πn+1=πnP\pi_{n+1} = \pi_n P et πn=π0Pn\pi_n = \pi_0 P^n.
  • Cas de deux états : si P=(1−aab1−b)P = \begin{pmatrix} 1 - a & a \\ b & 1 - b \end{pmatrix} avec 0<a+b<20 < a + b < 2, la suite (πn)(\pi_n) converge vers l'unique distribution invariante π=(ba+baa+b)\pi = \left(\dfrac{b}{a + b} \quad \dfrac{a}{a + b}\right), quelle que soit la loi initiale.

5. Suites de matrices

  • Si Un+1=AUnU_{n+1} = AU_n (avec UnU_n matrice colonne), alors Un=AnU0U_n = A^n U_0 par récurrence.
  • Si Un+1=AUn+BU_{n+1} = AU_n + B et si CC vérifie C=AC+BC = AC + B, alors Vn=Un−CV_n = U_n - C vérifie Vn+1=AVnV_{n+1} = AV_n, donc Un=An(U0−C)+CU_n = A^n(U_0 - C) + C.
  • Puissances d'une matrice diagonale : (a00b)n=(an00bn)\begin{pmatrix} a & 0 \\ 0 & b \end{pmatrix}^n = \begin{pmatrix} a^n & 0 \\ 0 & b^n \end{pmatrix}. Si A=QDQ−1A = QDQ^{-1} avec DD diagonale, alors An=QDnQ−1A^n = QD^nQ^{-1}.

6. Démonstration : nombre de chemins de longueur nn

On note aij(n)a_{ij}^{(n)} le coefficient (i,j)(i, j) de AnA^n et on raisonne par récurrence sur n≥1n \ge 1.

  • Initialisation : pour n=1n = 1, aija_{ij} est le nombre d'arêtes, c'est-à-dire de chemins de longueur 11, de ii à jj.
  • Hérédité : supposons la propriété vraie au rang nn. Un chemin de longueur n+1n + 1 de ii à jj se décompose de façon unique en un chemin de longueur nn de ii à un sommet kk, suivi d'une arête de kk à jj. Pour kk fixé, il y en a aik(n)×akja_{ik}^{(n)} \times a_{kj} (principe multiplicatif) ; en additionnant sur tous les sommets kk (cas disjoints), le nombre de chemins est

∑k=1paik(n) akj\sum_{k=1}^{p} a_{ik}^{(n)}\, a_{kj}
C'est exactement le coefficient (i,j)(i, j) du produit An×A=An+1A^n \times A = A^{n+1}.

  • Conclusion : pour tout n≥1n \ge 1, le coefficient (i,j)(i, j) de AnA^n est le nombre de chemins de longueur nn de ii à jj.

7. Démonstration : distribution après nn transitions d'une chaîne de Markov

  • Une transition : d'après la formule des probabilités totales, pour chaque état jj,

P(Xn+1=j)=∑iP(Xn=i)×pijP(X_{n+1} = j) = \sum_{i} P(X_n = i) \times p_{ij}
ce qui s'écrit matriciellement πn+1=πnP\pi_{n+1} = \pi_n P.

  • Récurrence : π0=π0P0\pi_0 = \pi_0 P^0 et, si πn=π0Pn\pi_n = \pi_0 P^n, alors πn+1=πnP=π0Pn+1\pi_{n+1} = \pi_n P = \pi_0 P^{n+1}. Donc πn=π0Pn\pi_n = \pi_0 P^n pour tout n∈Nn \in \mathbb{N}.
  • Interprétation de PnP^n : si l'on part de l'état ii (distribution initiale avec un 11 en position ii et des 00 ailleurs), πn\pi_n est la ligne ii de PnP^n. Le coefficient (i,j)(i, j) de PnP^n est donc la probabilité de passer de l'état ii à l'état jj en nn transitions.

8. Matrices et transformations du plan

Une transformation du plan qui s'écrit (x′y′)=A(xy)\begin{pmatrix} x' \\ y' \end{pmatrix} = A\begin{pmatrix} x \\ y \end{pmatrix} est représentée par la matrice AA. Exemples :

  • rotation de centre OO et d'angle θ\theta : A=(cos⁡θ−sin⁡θsin⁡θcos⁡θ)A = \begin{pmatrix} \cos\theta & -\sin\theta \\ \sin\theta & \cos\theta \end{pmatrix} ; pour θ=π2\theta = \frac{\pi}{2}, le point (1;0)(1 ; 0) a pour image (0;1)(0 ; 1) ;
  • homothétie de centre OO et de rapport kk : A=(k00k)A = \begin{pmatrix} k & 0 \\ 0 & k \end{pmatrix} ;
  • symétrie par rapport à l'axe des abscisses : A=(100−1)A = \begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix}.

Appliquer la transformation de matrice AA, puis celle de matrice BB, revient à appliquer la transformation de matrice BABA.

Méthodes et exemples rédigés

Méthode : Résoudre un système linéaire par inversion matricielle

Exemple

Résoudre le système {3x+2y=75x+4y=11\begin{cases} 3x + 2y = 7 \\ 5x + 4y = 11 \end{cases}.

  1. Étape 1 (Forme matricielle) : AX=BAX = B avec A=(3254)A = \begin{pmatrix} 3 & 2 \\ 5 & 4 \end{pmatrix}, X=(xy)X = \begin{pmatrix} x \\ y \end{pmatrix} et B=(711)B = \begin{pmatrix} 7 \\ 11 \end{pmatrix}.

  2. Étape 2 (Inversion) : det⁡(A)=3(4)−2(5)=12−10=2≠0\det(A) = 3(4) - 2(5) = 12 - 10 = 2 \neq 0. A−1=12(4−2−53)A^{-1} = \frac{1}{2} \begin{pmatrix} 4 & -2 \\ -5 & 3 \end{pmatrix}.

  3. Étape 3 (Calcul de XX) : X=A−1B=12(4(7)−2(11)−5(7)+3(11))=12(6−2)=(3−1)X = A^{-1} B = \frac{1}{2} \begin{pmatrix} 4(7) - 2(11) \\ -5(7) + 3(11) \end{pmatrix} = \frac{1}{2} \begin{pmatrix} 6 \\ -2 \end{pmatrix} = \begin{pmatrix} 3 \\ -1 \end{pmatrix}.

  4. Conclusion : x=3x = 3 et y=−1y = -1.

Méthode : Compter des chemins dans un graphe

Exemple

Un graphe non orienté a pour sommets AA, BB, CC, DD et pour arêtes AA–BB, AA–CC, AA–DD, BB–CC, CC–DD. Combien y a-t-il de chemins de longueur 3 de BB à DD ?

  1. Étape 1 (Matrice d'adjacence) : dans l'ordre A,B,C,DA, B, C, D :
    M=(0111101011011010)M = \begin{pmatrix} 0 & 1 & 1 & 1 \\ 1 & 0 & 1 & 0 \\ 1 & 1 & 0 & 1 \\ 1 & 0 & 1 & 0 \end{pmatrix}

  2. Étape 2 (Puissances) : à la calculatrice ou à la main,
    M2=(3121121221311212)M3=(4555525255455252)M^2 = \begin{pmatrix} 3 & 1 & 2 & 1 \\ 1 & 2 & 1 & 2 \\ 2 & 1 & 3 & 1 \\ 1 & 2 & 1 & 2 \end{pmatrix} \qquad M^3 = \begin{pmatrix} 4 & 5 & 5 & 5 \\ 5 & 2 & 5 & 2 \\ 5 & 5 & 4 & 5 \\ 5 & 2 & 5 & 2 \end{pmatrix}

  3. Étape 3 (Lecture) : le coefficient ligne BB, colonne DD de M3M^3 vaut 22. Ce sont les chemins BB–AA–CC–DD et BB–CC–AA–DD.

  4. Conclusion : il y a 22 chemins de longueur 3 de BB à DD. (Remarque : la diagonale de M2M^2 donne les degrés : 3,2,3,23, 2, 3, 2.)

Méthode : Étudier une chaîne de Markov à deux états

Exemple

Chaque jour, un abonné utilise le service AA ou le service BB. S'il utilise AA, il garde AA le lendemain avec probabilité 0,90{,}9 ; s'il utilise BB, il passe à AA avec probabilité 0,30{,}3. Le jour 00, π0=(0,50,5)\pi_0 = (0{,}5 \quad 0{,}5). Calculer π2\pi_2 puis l'état stable.

  1. Étape 1 (Matrice de transition) : P=(0,90,10,30,7)P = \begin{pmatrix} 0{,}9 & 0{,}1 \\ 0{,}3 & 0{,}7 \end{pmatrix} (chaque ligne a pour somme 11).

  2. Étape 2 (Évolution) : π1=π0P=(0,45+0,150,05+0,35)=(0,60,4)\pi_1 = \pi_0 P = (0{,}45 + 0{,}15 \quad 0{,}05 + 0{,}35) = (0{,}6 \quad 0{,}4), puis π2=π1P=(0,54+0,120,06+0,28)=(0,660,34)\pi_2 = \pi_1 P = (0{,}54 + 0{,}12 \quad 0{,}06 + 0{,}28) = (0{,}66 \quad 0{,}34).

  3. Étape 3 (État stable) : on cherche π=(xy)\pi = (x \quad y) avec πP=π\pi P = \pi et x+y=1x + y = 1. La première coordonnée donne 0,9x+0,3y=x0{,}9x + 0{,}3y = x, soit 0,3y=0,1x0{,}3y = 0{,}1x, donc x=3yx = 3y. Avec x+y=1x + y = 1 : y=0,25y = 0{,}25 et x=0,75x = 0{,}75.

  4. Conclusion : π2=(0,660,34)\pi_2 = (0{,}66 \quad 0{,}34) et l'état stable est π=(0,750,25)\pi = (0{,}75 \quad 0{,}25) : à long terme, 75 %75\,\% des jours se font avec le service AA.

Méthode : Calculer les puissances d'une matrice et en déduire une suite

Exemple

Soit A=(1201)A = \begin{pmatrix} 1 & 2 \\ 0 & 1 \end{pmatrix}. Conjecturer puis démontrer l'expression de AnA^n, et en déduire UnU_n si Un+1=AUnU_{n+1} = AU_n et U0=(13)U_0 = \begin{pmatrix} 1 \\ 3 \end{pmatrix}.

  1. Étape 1 (Conjecture) : A2=(1401)A^2 = \begin{pmatrix} 1 & 4 \\ 0 & 1 \end{pmatrix} et A3=(1601)A^3 = \begin{pmatrix} 1 & 6 \\ 0 & 1 \end{pmatrix} : on conjecture An=(12n01)A^n = \begin{pmatrix} 1 & 2n \\ 0 & 1 \end{pmatrix}.

  2. Étape 2 (Récurrence) : c'est vrai pour n=0n = 0 (A0=I2A^0 = I_2). Si Ak=(12k01)A^k = \begin{pmatrix} 1 & 2k \\ 0 & 1 \end{pmatrix}, alors Ak+1=AkA=(12+2k01)=(12(k+1)01)A^{k+1} = A^k A = \begin{pmatrix} 1 & 2 + 2k \\ 0 & 1 \end{pmatrix} = \begin{pmatrix} 1 & 2(k+1) \\ 0 & 1 \end{pmatrix}.

  3. Étape 3 (Suite) : Un=AnU0=(1×1+2n×30×1+1×3)=(1+6n3)U_n = A^n U_0 = \begin{pmatrix} 1 \times 1 + 2n \times 3 \\ 0 \times 1 + 1 \times 3 \end{pmatrix} = \begin{pmatrix} 1 + 6n \\ 3 \end{pmatrix}.

  4. Conclusion : si Un=(unvn)U_n = \begin{pmatrix} u_n \\ v_n \end{pmatrix}, alors un=1+6nu_n = 1 + 6n et vn=3v_n = 3. Contrôle : la relation Un+1=AUnU_{n+1} = AU_n s'écrit un+1=un+2vn=un+6u_{n+1} = u_n + 2v_n = u_n + 6.

Pièges fréquents

Questions de révision

Essaie de répondre de tête avant d'afficher la réponse.

Quelle condition sur det⁡(A)\det(A) assure l'inversibilité d'une matrice carrée ?

Son déterminant doit être non nul : det⁡(A)≠0\det(A) \neq 0.

Quelle équation définit l'état stable π\pi d'une chaîne de Markov ?

πP=π\pi P = \pi avec la somme des composantes égale à 1.

Que représente le coefficient (i,j)(i, j) de AnA^n, où AA est la matrice d'adjacence d'un graphe ?

Le nombre de chemins de longueur nn allant du sommet ii au sommet jj.

Qu'est-ce qu'une matrice stochastique ?

Une matrice carrée à coefficients positifs dont la somme de chaque ligne vaut 11.

Donner l'inverse de (2111)\begin{pmatrix} 2 & 1 \\ 1 & 1 \end{pmatrix}.

Le déterminant vaut 11, donc l'inverse est (1−1−12)\begin{pmatrix} 1 & -1 \\ -1 & 2 \end{pmatrix}.

Pourquoi le coefficient (i,j)(i, j) de A2A^2 compte-t-il les chemins de longueur 22 de ii à jj ?

Il vaut ∑kaikakj\sum_k a_{ik}a_{kj} : pour chaque sommet intermédiaire kk, on multiplie le nombre d'arêtes de ii à kk par le nombre d'arêtes de kk à jj.

Exercices corrigés

Exercice 1 Produit matriciel

Palier 1 · découverte

Calculer le produit A×BA \times B avec A=(1234)A = \begin{pmatrix} 1 & 2 \\ 3 & 4 \end{pmatrix} et B=(2013)B = \begin{pmatrix} 2 & 0 \\ 1 & 3 \end{pmatrix}.

  1. (461012)\begin{pmatrix} 4 & 6 \\ 10 & 12 \end{pmatrix}
  2. (20312)\begin{pmatrix} 2 & 0 \\ 3 & 12 \end{pmatrix}
  3. (40312)\begin{pmatrix} 4 & 0 \\ 3 & 12 \end{pmatrix}
  4. (101246)\begin{pmatrix} 10 & 12 \\ 4 & 6 \end{pmatrix}
Un indice

Le coefficient ligne ii, colonne jj de A×BA \times B s'obtient avec la ligne ii de AA et la colonne jj de BB.

Voir la correction

Bonne réponse : A. (461012)\begin{pmatrix} 4 & 6 \\ 10 & 12 \end{pmatrix}

c11=1(2)+2(1)=4c_{11} = 1(2)+2(1) = 4, c12=1(0)+2(3)=6c_{12} = 1(0)+2(3) = 6, c21=3(2)+4(1)=10c_{21} = 3(2)+4(1) = 10, c22=3(0)+4(3)=12c_{22} = 3(0)+4(3) = 12.

Exercice 2 Matrice de transition d'un graphe probabiliste

Palier 2 · application

Un système a 2 états AA et BB. De AA, on va en AA avec proba 0,7 et en BB avec 0,3. De BB, on va en AA avec 0,4 et en BB avec 0,6. Quelle est la matrice de transition MM ?

  1. (0,70,30,40,6)\begin{pmatrix} 0{,}7 & 0{,}3 \\ 0{,}4 & 0{,}6 \end{pmatrix}
  2. (0,70,40,30,6)\begin{pmatrix} 0{,}7 & 0{,}4 \\ 0{,}3 & 0{,}6 \end{pmatrix}
  3. (0,30,70,60,4)\begin{pmatrix} 0{,}3 & 0{,}7 \\ 0{,}6 & 0{,}4 \end{pmatrix}
  4. (1001)\begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix}
Un indice

Chaque ligne correspond à l'état de départ, la somme des coefficients de chaque ligne vaut 1.

Voir la correction

Bonne réponse : A. (0,70,30,40,6)\begin{pmatrix} 0{,}7 & 0{,}3 \\ 0{,}4 & 0{,}6 \end{pmatrix}

Les lignes représentent les probabilités de transition depuis l'état courant : ligne 1 (0,7;0,3)(0{,}7 ; 0{,}3) et ligne 2 (0,4;0,6)(0{,}4 ; 0{,}6).

Exercice 3 État stable d'un graphe probabiliste

Palier 3 · approfondissement

Pour M=(0,70,30,40,6)M = \begin{pmatrix} 0{,}7 & 0{,}3 \\ 0{,}4 & 0{,}6 \end{pmatrix}, quel est l'état stable P=(xy)P = (x \quad y) vérifiant PM=PP M = P et x+y=1x + y = 1 ?

  1. P=(4737)P = \left(\frac{4}{7} \quad \frac{3}{7}\right)
  2. P=(0,50,5)P = (0{,}5 \quad 0{,}5)
  3. P=(0,70,3)P = (0{,}7 \quad 0{,}3)
  4. P=(3747)P = \left(\frac{3}{7} \quad \frac{4}{7}\right)
Un indice

0,7x+0,4y=x  ⟺  0,4y=0,3x  ⟺  3x=4y0{,}7x + 0{,}4y = x \iff 0{,}4y = 0{,}3x \iff 3x = 4y.

Voir la correction

Bonne réponse : A. P=(4737)P = \left(\frac{4}{7} \quad \frac{3}{7}\right)

0,7x+0,4y=x  ⟺  0,4y=0,3x  ⟺  y=34x0{,}7x + 0{,}4y = x \iff 0{,}4y = 0{,}3x \iff y = \frac{3}{4}x. Avec x+y=1x + y = 1, on obtient x+34x=1  ⟺  74x=1  ⟺  x=47x + \frac{3}{4}x = 1 \iff \frac{7}{4}x = 1 \iff x = \frac{4}{7}, d'où y=37y = \frac{3}{7}.

Exercice 4 Puissance n-ième d'une matrice diagonalisable

Palier 4 · défi

Si M=PDP−1M = P D P^{-1} avec D=(1000,5)D = \begin{pmatrix} 1 & 0 \\ 0 & 0{,}5 \end{pmatrix}, quelle est la matrice limite lim⁡n→+∞Dn\lim_{n \to +\infty} D^n ?

  1. (1000)\begin{pmatrix} 1 & 0 \\ 0 & 0 \end{pmatrix}
  2. (0000)\begin{pmatrix} 0 & 0 \\ 0 & 0 \end{pmatrix}
  3. (1001)\begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix}
  4. La limite n'existe pas
Un indice

Dn=(1n00(0,5)n)D^n = \begin{pmatrix} 1^n & 0 \\ 0 & (0{,}5)^n \end{pmatrix}.

Voir la correction

Bonne réponse : A. (1000)\begin{pmatrix} 1 & 0 \\ 0 & 0 \end{pmatrix}

Dn=diag(1n,0,5n)=diag(1,(0,5)n)D^n = \text{diag}(1^n, 0{,}5^n) = \text{diag}(1, (0{,}5)^n). Quand n→+∞n \to +\infty, (0,5)n→0(0{,}5)^n \to 0, donc la matrice limite est (1000)\begin{pmatrix} 1 & 0 \\ 0 & 0 \end{pmatrix}.