Mathematics · Book 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.