Mathematics · Livre 2 · Grades 10–12

Mathématiques du lycée

Mathématiques du lycée · Grades 10–12

30Matrices et graphes

Une matrice est un tableau rectangulaire de nombres, additionné et multiplié selon des règles conçues pour que l’algèbre matricielle représente la composition des applications linéaires. Les matrices résolvent les systèmes linéaires, pilotent les suites récurrentes couplées, et comptent les chemins dans les réseaux — les mathématiques derrière les moteurs de recherche et les algorithmes de plus court chemin.

30.1 Algèbre matricielle

Définition 30.1 (Matrice)

Une matrice m×nm \times n est un tableau de nombres réels à mm lignes et nn colonnes : A=(aij)A = (a_{ij}), où aija_{ij} est l’entrée de la ligne ii, colonne jj. Deux matrices de même taille s’additionnent entrée par entrée, et λA=(λaij)\lambda A = (\lambda a_{ij}).

Définition 30.2 (Produit matriciel)

Soit AA de taille m×nm \times n et BB de taille n×pn \times p. Le produit ABAB est la matrice m×pm \times p dont l’entrée (i,j)(i,j) est

(AB)ij=k=1naikbkj(AB)_{ij} = \sum_{k=1}^{n} a_{ik} b_{kj}

(la règle « ligne ii de AA multipliée par colonne jj de BB »).

Exemple 30.3

(1234)(0111)=(2347)\begin{pmatrix} 1 & 2\\ 3 & 4\end{pmatrix} \begin{pmatrix} 0 & 1\\ 1 & 1\end{pmatrix} = \begin{pmatrix} 2 & 3\\ 4 & 7\end{pmatrix}, tandis que (0111)(1234)=(3446)\begin{pmatrix} 0 & 1\\ 1 & 1\end{pmatrix} \begin{pmatrix} 1 & 2\\ 3 & 4\end{pmatrix} = \begin{pmatrix} 3 & 4\\ 4 & 6\end{pmatrix} : la multiplication matricielle n’est pas commutative.

Proposition 30.4 (Règles de l’algèbre matricielle)

Dès que les tailles rendent les produits significatifs :

(AB)C=A(BC),A(B+C)=AB+AC,(A+B)C=AC+BC,(AB)C = A(BC), \qquad A(B + C) = AB + AC, \qquad (A+B)C = AC + BC,

et la matrice identité InI_n (uns sur la diagonale, zéros ailleurs) satisfait ImA=AIn=AI_m A = A I_n = A pour AA de taille m×nm \times n.

Démonstration. Toutes se vérifient entrée par entrée à partir de Définition 30.2 ; l’associativité, la seule non triviale, revient à échanger deux sommes finies : ((AB)C)ij=l(kaikbkl)clj=kaik(lbklclj)=(A(BC))ij\bigl((AB)C\bigr)_{ij} = \sum_l \left(\sum_k a_{ik}b_{kl}\right) c_{lj} = \sum_k a_{ik} \left(\sum_l b_{kl} c_{lj}\right) = \bigl(A(BC)\bigr)_{ij}.

Définition 30.5 (Inverse)

Une matrice carrée AA de taille nn est inversible s’il existe une matrice BB avec AB=BA=InAB = BA = I_n ; BB est alors unique, notée A1A^{-1}.

Proposition 30.6 (Inverse d’une matrice 2×22\times2)

Soit A=(abcd)A = \begin{pmatrix} a & b\\ c & d\end{pmatrix} et detA=adbc\det A = ad - bc (le déterminant). Alors AA est inversible si et seulement si detA0\det A \neq 0, auquel cas

A1=1adbc(dbca).A^{-1} = \frac{1}{ad - bc}\begin{pmatrix} d & -b\\ -c & a\end{pmatrix}.

Démonstration. Un calcul donne A(dbca)=(dbca)A=(adbc)I2A \begin{pmatrix} d & -b\\ -c & a\end{pmatrix} = \begin{pmatrix} d & -b\\ -c & a\end{pmatrix} A = (ad - bc) I_2 ; si adbc0ad - bc \neq 0, diviser. Réciproquement, si adbc=0ad - bc = 0, les colonnes de AA sont proportionnelles, et de même les colonnes de ABAB pour toute BB ; mais les colonnes de I2I_2 ne le sont pas, donc aucune BB ne peut vérifier AB=I2AB = I_2.

Méthode 30.7 (Systèmes linéaires)

Le système {ax+by=ecx+dy=f\begin{cases} ax + by = e\\ cx + dy = f \end{cases} est l’équation matricielle AX=YAX = Y avec X=(xy)X = \begin{pmatrix} x \\ y\end{pmatrix}, Y=(ef)Y = \begin{pmatrix} e \\ f\end{pmatrix}. Si detA0\det A \neq 0, sa solution unique est X=A1YX = A^{-1}Y. Le même formalisme traite nn équations à nn inconnues.

30.2 Puissances de matrices et suites récurrentes

Définition 30.8

Pour une matrice carrée AA et kNk \in \N, Ak=A××AA^k = A \times \dots \times A (kk facteurs), avec A0=IA^0 = I.

Méthode 30.9 (Cas diagonal-plus-nilpotent et diagonalisable)

Deux façons standards de calculer AkA^k :

  • Si A=λI+NA = \lambda I + N avec N2=0N^2 = 0, la formule du binôme (valable ici car II et NN commutent) se réduit à deux termes : Ak=λkI+kλk1NA^k = \lambda^k I + k \lambda^{k-1} N.
  • Si l’on trouve une PP inversible avec A=PDP1A = PDP^{-1} et DD diagonale, alors Ak=PDkP1A^k = P D^k P^{-1}, et DkD^k se calcule entrée par entrée. (Trouver une telle PP de façon systématique est la théorie de la diagonalisation, développée à l’université ; à ce niveau PP est donnée.)

Exemple 30.10 (Suites couplées)

Soient un+1=3un+vnu_{n+1} = 3u_n + v_n et vn+1=un+3vnv_{n+1} = u_n + 3v_n. En posant Xn=(unvn)X_n = \begin{pmatrix} u_n\\ v_n \end{pmatrix} et A=(3113)A = \begin{pmatrix} 3 & 1\\ 1 & 3\end{pmatrix}, on obtient Xn+1=AXnX_{n+1} = AX_n, donc Xn=AnX0X_n = A^n X_0. Les suites auxiliaires sn=un+vns_n = u_n + v_n et dn=unvnd_n = u_n - v_n vérifient sn+1=4sns_{n+1} = 4s_n et dn+1=2dnd_{n+1} = 2d_n, donc sn=4ns0s_n = 4^n s_0, dn=2nd0d_n = 2^n d_0 et

un=4n(u0+v0)+2n(u0v0)2,vn=4n(u0+v0)2n(u0v0)2.u_n = \frac{4^n(u_0+v_0) + 2^n(u_0-v_0)}{2}, \qquad v_n = \frac{4^n(u_0+v_0) - 2^n(u_0-v_0)}{2}.

(En coulisse : (1,1)(1,1) et (1,1)(1,-1) sont des directions propres de AA.)

30.3 Graphes et chemins

Définition 30.11 (Graphe, matrice d’adjacence)

Un graphe est constitué de sommets 1,2,,n1, 2, \dots, n et d’arêtes joignant certaines paires de sommets (couples ordonnés pour un graphe orienté). Sa matrice d’adjacence est la matrice n×nn \times n MM avec mij=1m_{ij} = 1 s’il y a une arête de ii vers jj, et 00 sinon. Un chemin de longueur kk de ii vers jj est une suite de kk arêtes consécutives menant de ii à jj.

M = pmatrix 0 & 1 & 1\\ 0 & 0 & 1\\ 1 & 0 & 0 pmatrix Un graphe orienté et sa matrice d’adjacence () : m_ij = 1 exactement lorsqu’il y a une arête de i vers j.
M=(011001100)M = \begin{pmatrix} 0 & 1 & 1\\ 0 & 0 & 1\\ 1 & 0 & 0 \end{pmatrix} Un graphe orienté et sa matrice d’adjacence (Exercice 30.6) : mij=1m_{ij} = 1 exactement lorsqu’il y a une arête de ii vers jj.

Théorème 30.12 (Dénombrement des chemins)

Le nombre de chemins de longueur kk du sommet ii au sommet jj est l’entrée (i,j)(i,j) de MkM^k.

Démonstration. Récurrence sur kk. Pour k=1k = 1 c’est la définition de MM. Supposer le résultat pour kk. Un chemin de longueur k+1k+1 de ii à jj est un chemin de longueur kk de ii vers un certain sommet ll, suivi d’une arête de ll vers jj ; par les principes d’addition et de multiplication, leur nombre est

l=1n(Mk)ilmlj=(Mk+1)ij.\sum_{l=1}^{n} \bigl(M^k\bigr)_{il}\, m_{lj} = \bigl(M^{k+1}\bigr)_{ij}. \qedhere

Exemple 30.13

Pour le graphe triangle (33 sommets, toutes les paires jointes), M=(011101110)M = \begin{pmatrix} 0&1&1\\ 1&0&1\\ 1&1&0\end{pmatrix} et M2=(211121112)M^2 = \begin{pmatrix} 2&1&1\\ 1&2&1\\ 1&1&2\end{pmatrix} : depuis chaque sommet il y a 22 chemins de longueur 22 vers soi-même (via l’un ou l’autre voisin) et 11 vers chaque autre sommet.

30.4 Exercices

Exercice 30.1

Soient A=(1201)A = \begin{pmatrix} 1 & 2\\ 0 & 1 \end{pmatrix} et B=(2011)B = \begin{pmatrix} 2 & 0\\ 1 & 1 \end{pmatrix}. Calculer A+BA + B, ABAB, BABA et A2A^2.

Solution

Solution de Exercice 30.1.

A+B=(3212),AB=(4211),BA=(2413),A2=(1401).A + B = \begin{pmatrix} 3 & 2\\ 1 & 2\end{pmatrix}, \quad AB = \begin{pmatrix} 4 & 2\\ 1 & 1\end{pmatrix}, \quad BA = \begin{pmatrix} 2 & 4\\ 1 & 3\end{pmatrix}, \quad A^2 = \begin{pmatrix} 1 & 4\\ 0 & 1\end{pmatrix}.

Noter ABBAAB \neq BA.

Exercice 30.2

Déterminer si les matrices suivantes sont inversibles, et calculer les inverses lorsqu’ils existent :

A=(2513),B=(3624).A = \begin{pmatrix} 2 & 5\\ 1 & 3\end{pmatrix}, \qquad B = \begin{pmatrix} 3 & 6\\ 2 & 4\end{pmatrix}.
Solution

Solution de Exercice 30.2.

detA=65=10\det A = 6 - 5 = 1 \neq 0 : A1=(3512)A^{-1} = \begin{pmatrix} 3 & -5\\ -1 & 2 \end{pmatrix}. detB=1212=0\det B = 12 - 12 = 0 : BB n’est pas inversible.

Exercice 30.3

Résoudre par inversion matricielle le système {2x+5y=1x+3y=2.\begin{cases} 2x + 5y = 1\\ x + 3y = 2 . \end{cases}

Solution

Solution de Exercice 30.3.

Le système est AX=YAX = Y avec AA comme dans l’Exercice 30.2 et Y=(12)Y = \begin{pmatrix} 1\\ 2\end{pmatrix} :

X=A1Y=(3512)(12)=(73):x=7, y=3.X = A^{-1}Y = \begin{pmatrix} 3 & -5\\ -1 & 2\end{pmatrix} \begin{pmatrix} 1\\ 2\end{pmatrix} = \begin{pmatrix} -7\\ 3\end{pmatrix}: \qquad x = -7,\ y = 3 .

Exercice 30.4 ★★

Soit A=(2102)=2I+NA = \begin{pmatrix} 2 & 1\\ 0 & 2\end{pmatrix} = 2I + N avec N=(0100)N = \begin{pmatrix} 0 & 1\\ 0 & 0 \end{pmatrix}.

  1. Vérifier que N2=0N^2 = 0 et que II et NN commutent.
  2. En déduire AkA^k pour tout kNk \in \N et vérifier la formule pour k=2k=2 par calcul direct.
Solution

Solution de Exercice 30.4.

1. N2=(0100)(0100)=0N^2 = \begin{pmatrix} 0&1\\0&0\end{pmatrix} \begin{pmatrix} 0&1\\0&0\end{pmatrix} = 0, et II commute avec toute matrice.

2. Comme les deux termes commutent, la formule du binôme s’applique et tous les termes contenant N2N^2 s’annulent :

Ak=(2I+N)k=2kI+k2k1N=(2kk2k102k).A^k = (2I + N)^k = 2^k I + k\,2^{k-1} N = \begin{pmatrix} 2^k & k\,2^{k-1}\\ 0 & 2^k \end{pmatrix}.

Vérification pour k=2k = 2 : A2=(2102)2=(4404)A^2 = \begin{pmatrix} 2&1\\0&2\end{pmatrix}^2 = \begin{pmatrix} 4&4\\0&4\end{pmatrix}, et la formule donne 22=42^2 = 4, 2×2=42 \times 2 = 4. ✓

Exercice 30.5 ★★

Soient A=(0111)A = \begin{pmatrix} 0 & 1\\ 1 & 1\end{pmatrix} et FnF_n la suite de Fibonacci (F0=0F_0 = 0, F1=1F_1 = 1, Fn+2=Fn+1+FnF_{n+2} = F_{n+1} + F_n). Montrer par récurrence que pour n1n \geq 1,

An=(Fn1FnFnFn+1),A^n = \begin{pmatrix} F_{n-1} & F_n\\ F_n & F_{n+1}\end{pmatrix},

et en déduire l’identité Fn+1Fn1Fn2=(1)nF_{n+1}F_{n-1} - F_n^2 = (-1)^n. (Indication : les déterminants se multiplient : det(MN)=detMdetN\det(MN) = \det M \det N, ce que l’on peut vérifier pour les matrices 2×22\times2.)

Solution

Solution de Exercice 30.5.

Récurrence. Pour n=1n = 1 : A1=(0111)=(F0F1F1F2)A^1 = \begin{pmatrix} 0&1\\1&1\end{pmatrix} = \begin{pmatrix} F_0 & F_1\\ F_1 & F_2\end{pmatrix}. Supposer la formule pour nn ; alors

An+1=AnA=(Fn1FnFnFn+1)(0111)=(FnFn1+FnFn+1Fn+Fn+1)=(FnFn+1Fn+1Fn+2).A^{n+1} = A^n A = \begin{pmatrix} F_{n-1} & F_n\\ F_n & F_{n+1}\end{pmatrix} \begin{pmatrix} 0 & 1\\ 1 & 1\end{pmatrix} = \begin{pmatrix} F_n & F_{n-1} + F_n\\ F_{n+1} & F_n + F_{n+1}\end{pmatrix} = \begin{pmatrix} F_n & F_{n+1}\\ F_{n+1} & F_{n+2}\end{pmatrix}.

Identité. Pour les matrices 2×22\times2, le développement montre det(MN)=detMdetN\det(MN) = \det M \det N ; d’où det(An)=(detA)n=(1)n\det(A^n) = (\det A)^n = (-1)^n, et detAn=Fn1Fn+1Fn2\det A^n = F_{n-1}F_{n+1} - F_n^2. (C’est l’identité de Cassini.)

Exercice 30.6 ★★

Un graphe orienté sur les sommets {1,2,3}\{1, 2, 3\} a les arêtes 121\to2, 232\to3, 313\to1 et 131\to3.

  1. Écrire la matrice d’adjacence MM et calculer M2M^2 et M3M^3.
  2. Combien de chemins de longueur 33 vont de 11 à 11 ? Les lister.
Solution

Solution de Exercice 30.6.

1. En ordonnant les sommets 1,2,31, 2, 3 :

M=(011001100),M2=(101100011),M3=(111101101).M = \begin{pmatrix} 0&1&1\\ 0&0&1\\ 1&0&0\end{pmatrix}, \quad M^2 = \begin{pmatrix} 1&0&1\\ 1&0&0\\ 0&1&1\end{pmatrix}, \quad M^3 = \begin{pmatrix} 1&1&1\\ 1&0&1\\ 1&0&1 \end{pmatrix}.

2. (M3)11=1\bigl(M^3\bigr)_{11} = 1 : exactement un chemin fermé de longueur 33 au sommet 11, à savoir 12311 \to 2 \to 3 \to 1. (Le chemin 1311 \to 3 \to 1 n’a que longueur 22, et 131 \to 3 puis 313\to1 puis 131\to3 se termine en 33.)

Exercice 30.7 ★★

Une société d’autopartage déplace des véhicules entre deux villes AA et BB. Chaque semaine, 80%80\% des voitures en AA restent en AA et 20%20\% passent en BB ; 30%30\% des voitures en BB passent en AA et 70%70\% restent. Soient an,bna_n, b_n les proportions de la flotte dans chaque ville.

  1. Écrire Xn+1=MXnX_{n+1} = MX_n avec Xn=(anbn)X_n = \begin{pmatrix} a_n\\ b_n\end{pmatrix} et identifier MM.
  2. Trouver les proportions d’équilibre (résoudre MX=XMX = X avec a+b=1a + b = 1).
  3. Montrer que cn=an0.6c_n = a_n - 0.6 vérifie cn+1=0.5cnc_{n+1} = 0.5\,c_n, et conclure que la répartition de la flotte converge vers l’équilibre.
Solution

Solution de Exercice 30.7.

1. an+1=0.8an+0.3bna_{n+1} = 0.8a_n + 0.3b_n, bn+1=0.2an+0.7bnb_{n+1} = 0.2a_n + 0.7b_n : M=(0.80.30.20.7)M = \begin{pmatrix} 0.8 & 0.3\\ 0.2 & 0.7\end{pmatrix}.

2. MX=XMX = X donne 0.8a+0.3b=a0.8a + 0.3b = a, c.-à-d. 0.3b=0.2a0.3b = 0.2a, donc b=23ab = \frac23 a ; avec a+b=1a + b = 1 : a=0.6a = 0.6, b=0.4b = 0.4.

3. En utilisant bn=1anb_n = 1 - a_n : an+1=0.8an+0.3(1an)=0.5an+0.3a_{n+1} = 0.8a_n + 0.3(1 - a_n) = 0.5a_n + 0.3, donc

cn+1=an+10.6=0.5an+0.30.6=0.5(an0.6)=0.5cn.c_{n+1} = a_{n+1} - 0.6 = 0.5a_n + 0.3 - 0.6 = 0.5(a_n - 0.6) = 0.5\,c_n .

D’où cn=0.5nc00c_n = 0.5^n c_0 \to 0 : an0.6a_n \to 0.6 et bn0.4b_n \to 0.4, quelle que soit la répartition initiale.

Exercice 30.8 ★★★

Soient A=(3113)A = \begin{pmatrix} 3 & 1\\ 1 & 3 \end{pmatrix}, P=(1111)P = \begin{pmatrix} 1 & 1\\ 1 & -1 \end{pmatrix}.

  1. Calculer P1P^{-1}, puis D=P1APD = P^{-1}AP, et vérifier que DD est diagonale.
  2. En déduire une formule fermée pour AnA^n et comparer avec Exemple 30.10.
Solution

Solution de Exercice 30.8.

1. detP=2\det P = -2, donc P1=12(1111)=12(1111)P^{-1} = -\frac12\begin{pmatrix} -1 & -1\\ -1 & 1\end{pmatrix} = \frac12\begin{pmatrix} 1 & 1\\ 1 & -1\end{pmatrix}. Puis

AP=(4242),D=P1AP=12(1111)(4242)=(4002).AP = \begin{pmatrix} 4 & 2\\ 4 & -2 \end{pmatrix}, \qquad D = P^{-1}AP = \frac12\begin{pmatrix} 1&1\\1&-1\end{pmatrix} \begin{pmatrix} 4&2\\4&-2\end{pmatrix} = \begin{pmatrix} 4 & 0\\ 0 & 2\end{pmatrix}.

2. De A=PDP1A = PDP^{-1}, une récurrence immédiate donne An=PDnP1A^n = PD^nP^{-1} avec Dn=(4n002n)D^n = \begin{pmatrix} 4^n & 0\\ 0 & 2^n\end{pmatrix}, donc

An=PDnP1=(4n2n4n2n)12(1111)=12(4n+2n4n2n4n2n4n+2n).A^n = P D^n P^{-1} = \begin{pmatrix} 4^n & 2^n\\ 4^n & -2^n\end{pmatrix}\cdot \frac12\begin{pmatrix} 1&1\\1&-1\end{pmatrix} = \frac12\begin{pmatrix} 4^n + 2^n & 4^n - 2^n\\ 4^n - 2^n & 4^n + 2^n\end{pmatrix}.

Appliquer AnA^n à X0=(u0v0)X_0 = \begin{pmatrix} u_0\\v_0\end{pmatrix} reproduit exactement les formules de l’Exemple 30.10.

30.5 Problème : la matrice qui connaît Fibonacci (et la météo)

Problème 30.1

Devoir du week-end — une seule matrice 2×22 \times 2 porte tout Fibonacci, une matrice de Markov prévoit la météo à long terme, et un vecteur propre vaut un milliard

Une matrice est une machine qui avale un état et renvoie le suivant — et ses puissances contiennent donc des avenirs entiers. Ce problème s’ouvre sur l’étonnante matrice dont les puissances égrènent les nombres de Fibonacci (et démontrent leurs identités en une ligne chacune), fait ensuite tourner la météo comme une chaîne de Markov jusqu’à son état stationnaire, et se referme sur le vecteur propre sur lequel un moteur de recherche a été bâti (Théorème 30.12, Méthode 30.9).

Partie I — Aisance.

  1. Avec A=(1234)A = \begin{pmatrix} 1 & 2\\ 3 & 4\end{pmatrix} et B=(0110)B = \begin{pmatrix} 0 & 1\\ 1 & 0\end{pmatrix} : calculer ABAB et BABA. Verdict sur la commutativité ?
  2. Inverser (2153)\begin{pmatrix} 2 & 1\\ 5 & 3\end{pmatrix} (Proposition 30.6) et se servir de l’inverse pour résoudre 2x+y=42x + y = 4, 5x+3y=75x + 3y = 7.
  3. Soit N=(0100)N = \begin{pmatrix} 0 & 1\\ 0 & 0\end{pmatrix} : calculer N2N^2, et en déduire (I+N)n=I+nN(I + N)^n = I + nN pour tout nn.
  4. Le graphe triangle (trois sommets, toutes les paires reliées) : écrire sa matrice d’adjacence AA, calculer A3A^3, et interpréter les coefficients diagonaux (Théorème 30.12).
  5. Pour D=(20012)D = \begin{pmatrix} 2 & 0\\ 0 & \frac12 \end{pmatrix} : donner DnD^n et son comportement quand nn \to \infty.

Partie II — La matrice de Fibonacci. Posons F=(1110)F = \begin{pmatrix} 1 & 1\\ 1 & 0\end{pmatrix} et notons F1=F2=1,F3=2,F_1 = F_2 = 1, F_3 = 2, \dots les nombres de Fibonacci du Problème 13.1.

  1. Calculer F2F^2, F3F^3, F4F^4 et conjecturer la forme générale de FnF^n en fonction des nombres de Fibonacci.
  2. Démontrer la conjecture Fn=(Fn+1FnFnFn1)F^n = \begin{pmatrix} F_{n+1} & F_n\\ F_n & F_{n-1} \end{pmatrix} par récurrence.
  3. Prendre le déterminant des deux membres (le déterminant d’un produit est le produit des déterminants — le vérifier sur des matrices 2×22 \times 2 si vous ne l’avez jamais vu) : en déduire l’identité de Cassini Fn+1Fn1Fn2=(1)nF_{n+1}F_{n-1} - F_n^2 = (-1)^n — le moteur du carré évanoui, démontré en une ligne.
  4. À partir de Fm+n=FmFnF^{m+n} = F^m F^n, lire les coefficients en haut à droite et en déduire la formule d’addition

    Fm+n=Fm+1Fn+FmFn1.F_{m+n} = F_{m+1} F_n + F_m F_{n-1} .

    La vérifier pour m=n=3m = n = 3.

  5. Déduire de la formule d’addition (récurrence sur kk) que FnF_n divise FknF_{kn}, et le vérifier sur F3F6F_3 \mid F_6 et F3F9F_3 \mid F_9.
  6. Pour calculer F100F_{100}, nul besoin de multiplier 100100 matrices : élever au carré de façon répétée (F2,F4,F8,F^2, F^4, F^8, \dots) puis combiner. Combien de multiplications de matrices suffisent, et de quelle antique astuce de multiplication du volume précédent s’agit-il, promue aux matrices ?

Partie III — La machine à météo. Dans une certaine ville : après une journée ensoleillée, la suivante est ensoleillée avec la probabilité 0.80.8 ; après une journée pluvieuse, elle est ensoleillée avec la probabilité 0.40.4. On code la loi du jour par une colonne (psoleilppluie)\binom{p_{\text{soleil}}}{p_{\text{pluie}}} et l’évolution par

M=(0.80.40.20.6).M = \begin{pmatrix} 0.8 & 0.4\\ 0.2 & 0.6 \end{pmatrix}.
  1. Vérifier que chaque colonne de MM a pour somme 11, et dire pourquoi toute machine à météo doit posséder cette propriété.
  2. Aujourd’hui il fait soleil. Calculer la prévision pour demain, puis pour après-demain.
  3. Déterminer l’état stationnaire : la loi vv telle que Mv=vMv = v (et dont les coefficients ont pour somme 11). Quelle fraction des jours est ensoleillée à long terme ?
  4. Partir d’une journée pluvieuse, (01)\binom01, et appliquer MM quatre fois en suivant à chaque pas la distance à l’état stationnaire. Par quel facteur l’écart se réduit-il à chaque pas — et de quel type de convergence s’agit-il ?
  5. Le PageRank en miniature : trois pages, avec les liens ABA \to B, ACA \to C, BCB \to C, CAC \to A. Un internaute aléatoire suit un lien sortant uniformément au hasard. Écrire la matrice de transition, déterminer l’état stationnaire, et classer les pages.
  6. Interpréter le classement : pourquoi CC obtient-il un score aussi élevé que AA alors qu’il reçoit des liens de moins de pages — que mesure au juste l’état stationnaire ? (Le vrai PageRank ajoute un facteur d’amortissement pour les impasses et les sauts ; l’idée du vecteur propre est exactement celle-ci.)

Partie IV — Les dividendes de la diagonale.

  1. Deux quantités couplées obéissent à un+1=3un+vnu_{n+1} = 3u_n + v_n, vn+1=un+3vnv_{n+1} = u_n + 3v_n, c’est-à-dire à la matrice AA de l’Exercice 30.8. À l’aide de la diagonalisation de cet exercice (D=diag(4,2)D = \operatorname{diag}(4, 2)), donner la formule close de unu_n pour u0=1u_0 = 1, v0=0v_0 = 0, et la confronter au calcul direct pour n=1,2,3n = 1, 2, 3.
  2. En une ou deux phrases : que fait la diagonalisation à un système couplé — et en quel sens l’état stationnaire de Markov de la question 14 est-il lui aussi une histoire de vecteur propre ?
  3. Pour finir — les trois visages de la matrice ce week-end : la tenue de comptes (systèmes et inverses), le dénombrement (chemins et liens comptés par les puissances) et l’évolution (Fibonacci, la météo, le web — des avenirs lus sur les directions propres). Une phrase pour chacun, plus l’annonce : l’algèbre linéaire des volumes universitaires fait de chacun de ces visages une théorie.
Solution

Solution de Problème 30.1.

1. AB=(2143)AB = \begin{pmatrix} 2 & 1\\ 4 & 3\end{pmatrix} et BA=(3412)BA = \begin{pmatrix} 3 & 4\\ 1 & 2\end{pmatrix} : la multiplication des matrices n’est pas commutative — BB échange les colonnes à droite, les lignes à gauche.

2. Déterminant 65=16 - 5 = 1, d’où l’inverse (3152)\begin{pmatrix} 3 & -1\\ -5 & 2\end{pmatrix}. En l’appliquant à (47)\binom{4}{7} : x=127=5x = 12 - 7 = 5, y=20+14=6y = -20 + 14 = -6.

3. N2=0N^2 = 0. Alors (I+N)n=I+nN(I + N)^n = I + nN par récurrence : (I+nN)(I+N)=I+(n+1)N+nN2=I+(n+1)N(I + nN)(I + N) = I + (n+1)N + nN^2 = I + (n+1)N.

4. A=(011101110)A = \begin{pmatrix} 0&1&1\\ 1&0&1\\ 1&1&0 \end{pmatrix}, et A3A^3 a pour coefficients diagonaux 22 : depuis chaque sommet, exactement deux chemins fermés de longueur 33 (le triangle parcouru dans un sens ou dans l’autre) — le théorème de dénombrement à l’œuvre.

5. Dn=(2n002n)D^n = \begin{pmatrix} 2^n & 0\\ 0 & 2^{-n} \end{pmatrix} : une direction explose, l’autre s’éteint — sur la diagonale, les destins sont des suites géométriques indépendantes.

6. F2=(2111)F^2 = \begin{pmatrix} 2 & 1\\ 1 & 1\end{pmatrix}, F3=(3221)F^3 = \begin{pmatrix} 3 & 2\\ 2 & 1\end{pmatrix}, F4=(5332)F^4 = \begin{pmatrix} 5 & 3\\ 3 & 2\end{pmatrix} : du Fibonacci partout, d’où la conjecture énoncée.

7. Si Fn=(Fn+1FnFnFn1)F^n = \begin{pmatrix} F_{n+1} & F_n\\ F_n & F_{n-1}\end{pmatrix}, alors

Fn+1=FnF=(Fn+1+FnFn+1Fn+Fn1Fn)=(Fn+2Fn+1Fn+1Fn):F^{n+1} = F^n F = \begin{pmatrix} F_{n+1} + F_n & F_{n+1}\\ F_n + F_{n-1} & F_n \end{pmatrix} = \begin{pmatrix} F_{n+2} & F_{n+1}\\ F_{n+1} & F_n \end{pmatrix} :

hérédité ; l’initialisation en n=1n = 1 n’est autre que FF lui-même, avec la convention F0=0F_0 = 0 (qui prolonge la récurrence vers l’arrière).

8. detF=1\det F = -1, donc det(Fn)=(detF)n=(1)n\det(F^n) = (\det F)^n = (-1)^n ; et directement det(Fn)=Fn+1Fn1Fn2\det(F^n) = F_{n+1}F_{n-1} - F_n^2 : c’est Cassini, en une ligne. (La règle du produit pour les déterminants 2×22 \times 2 est un agréable développement de cinq minutes.)

9. Coefficient en haut à droite de FmFnF^m F^n : Fm+1Fn+FmFn1F_{m+1}F_n + F_m F_{n-1} ; celui de Fm+nF^{m+n} : Fm+nF_{m+n}. Pour m=n=3m = n = 3 : F4F3+F3F2=3×2+2×1=8=F6F_4 F_3 + F_3 F_2 = 3 \times 2 + 2 \times 1 = 8 = F_6.

10. Pour k=1k = 1 : c’est immédiat. Si FnFknF_n \mid F_{kn}, la formule d’addition avec m=knm = kn donne F(k+1)n=Fkn+1Fn+FknFn1F_{(k+1)n} = F_{kn+1}F_n + F_{kn}F_{n-1} : les deux termes sont des multiples de FnF_n. Donc FnFknF_n \mid F_{kn} pour tout kk : vérification, F3=2F_3 = 2 divise F6=8F_6 = 8 et F9=34F_9 = 34.

11. F100=F64F32F4F^{100} = F^{64} F^{32} F^4 : sept élévations au carré (F2,F4,,F64F^2, F^4, \dots, F^{64}) plus deux combinaisons — neuf multiplications au lieu de quatre-vingt-dix-neuf. C’est l’astuce de la table de doublements des scribes égyptiens, transposée des nombres aux matrices : écrire 100100 en binaire, puis multiplier les doublements utiles.

12. 0.8+0.2=10.8 + 0.2 = 1 et 0.4+0.6=10.4 + 0.6 = 1 : demain il fera un temps ou un autre — chaque colonne est une loi de probabilité complète, si bien que les probabilités se conservent.

13. Demain : (0.80.2)\binom{0.8}{0.2}. Après-demain : M(0.80.2)=(0.720.28)M\binom{0.8}{0.2} = \binom{0.72}{0.28}.

14. Mv=vMv = v avec v=(sr)v = \binom{s}{r} et s+r=1s + r = 1 : 0.8s+0.4r=s0.8s + 0.4r = s donne 0.4r=0.2s0.4r = 0.2s, soit s=2rs = 2r, d’où v=(2/31/3)v = \binom{2/3}{1/3}. À long terme, deux jours sur trois sont ensoleillés — quel que soit le temps qu’il fait aujourd’hui.

15. À partir de (01)\binom01 : composantes ensoleillées 0.40.4, 0.560.56, 0.6240.624, 0.64960.6496 ; écarts à 23\frac23 : 0.2670.267, 0.1070.107, 0.0430.043, 0.0170.017 — chaque pas multiplie l’écart par exactement 0.40.4 (la seconde valeur propre de la machine) : convergence géométrique vers l’état stationnaire.

16. Colonnes (issues de AA, BB, CC) : P=(00112001210)P = \begin{pmatrix} 0 & 0 & 1\\ \frac12 & 0 & 0\\ \frac12 & 1 & 0\end{pmatrix}. État stationnaire : vA=vCv_A = v_C, vB=vA2v_B = \frac{v_A}{2}, vC=vA2+vBv_C = \frac{v_A}{2} + v_B ; en sommant à 11 : v=(25,15,25)v = \left(\frac25, \frac15, \frac25\right). Classement : AA et CC à égalité en tête, BB dernière.

17. CC reçoit tout le trafic de BB et la moitié de celui de AA, et il renvoie tout vers AA : l’état stationnaire mesure où l’internaute passe du temps, et non combien de liens pointent vers la page — un lien depuis une page fréquentée pèse plus que plusieurs liens depuis des pages désertes. Cette pondération récursive est précisément l’idée fondatrice du moteur ; l’amortissement gère les pièges et les impasses.

18. An=PDnP1A^n = P D^n P^{-1} donne un=4n+2n2u_n = \frac{4^n + 2^n}{2} (et vn=4n2n2v_n = \frac{4^n - 2^n}{2}). Vérification : u1=3u_1 = 3, u2=10u_2 = 10, u3=36u_3 = 36 ; directement, (1,0)(3,1)(10,6)(36,28)(1,0) \to (3,1) \to (10,6) \to (36, 28) : cela concorde.

19. La diagonalisation change de coordonnées pour des coordonnées dans lesquelles le système couplé se disloque en suites géométriques indépendantes — chaque valeur propre court sa propre course. L’état stationnaire de Markov est le vecteur propre associé à la valeur propre 11, et la vitesse de convergence de la question 15 est la valeur propre suivante : la machine à météo était depuis le début une histoire d’éléments propres.

20. Tenue de comptes : un système est une seule équation matricielle, résolue par un seul inverse. Dénombrement : les puissances de la matrice d’adjacence comptent les chemins, les liens et les connexions. Évolution : les puissances de la machine conduisent les états vers leurs destins, et les directions propres (la direction dorée de Fibonacci, l’état stationnaire de la météo, le vecteur de classement du web) sont ces destins. L’algèbre linéaire, dans les volumes universitaires, est la science de tout cela.

Termes définis dans ce chapitre

Voir les 395 termes du glossaire