Algèbre linéaire 201-NYC / Maths 105 • Sciences humaines et cégep à Montréal

Fiche de révision : modèles matriciels, Markov, Leontief et Leslie

Cette fiche traite les trois modèles matriciels que le cours d'algèbre linéaire du cégep utilise pour montrer à quoi servent les matrices : les chaînes de Markov, le modèle entrées-sorties de Leontief et les matrices de Leslie. Elle ne redonne pas le calcul matriciel, elle donne les huit endroits où un produit juste répond à la mauvaise question.

Elle s'adresse aux étudiants de cégep qui suivent le cours 201-NYC, aussi appelé Maths 105, en sciences de la nature comme en sciences humaines. La série d'exercices corrigés du même chapitre reprend chaque modèle, du graphe de transition jusqu'à la lecture d'un tableau entrées-sorties.

Le fil du chapitre

Une matrice de modèle n'est pas un tableau de nombres : c'est une machine qui transporte un vecteur d'état d'une étape à la suivante. Presque toutes les fautes du chapitre viennent d'ignorer DANS QUEL SENS elle le transporte.

Ce chapitre fait partie de Algèbre linéaire et géométrie vectorielle, 201-NYC

Avant ce chapitre

Cette fiche suppose ces notions acquises. Si une méthode ci-dessous reste opaque, c'est presque toujours l'une d'elles qui manque, pas la fiche.

Remonter plus loin : la chaîne complète (4 chapitres) ↓

Le chemin de remédiation, du plus ancien au plus proche. Un élève qui reprend ce chapitre de zéro le reprend dans cet ordre.

  1. 1Matrices, systèmes et déterminants
  2. 2Systèmes linéaires : Gauss et Gauss-Jordan
  3. 3Déterminants, cofacteurs et règle de Cramer
  4. 4Diagonalisation et valeurs propres

L'essentiel

La convention, à déclarer avant tout calcul

  • Convention COLONNE : pijp_{ij} se lit « de l'état jj VERS l'état ii », et chaque COLONNE somme à 11.
  • En convention ligne, tout est transposé : pijp_{ij} va de ii vers jj et ce sont les LIGNES qui somment à 11.
  • Une étape : xn+1=Pxnx_{n+1}=Px_{n}. Après nn étapes : xn=Pnx0x_{n}=P^{n}x_{0}.
  • L'entrée (i,j)(i,j) de P2P^{2} additionne AUTOMATIQUEMENT tous les chemins de jj à ii en deux étapes : c'est ce qui rend la matrice utile.
convention COLONNE : de j vers iAB0,300,10reste : 0,70reste : 0,90
La flèche du haut porte pBA=0,30p_{BA}=0{,}30 : elle part de la colonne AA et arrive à la ligne BB, et le complément 0,700{,}70 reste sur place.

Contrôle immédiat de la convention : additionner une colonne. Si elle donne 11, on est en convention colonne et le produit s'écrit PxPx.

L'état stationnaire, et quand il existe

  • On cherche xx tel que Px=xPx=x, donc on résout (PI)x=0\left(P-I\right)x=0 et JAMAIS Px=0Px=0.
  • Le système est de rang déficient : il a une infinité de solutions, toutes proportionnelles. Il faut ensuite NORMALISER par la somme voulue.
  • Une chaîne RÉGULIÈRE, c'est-à-dire dont une puissance a toutes ses entrées strictement positives, converge vers un unique état stationnaire, quelle que soit la répartition de départ.
  • Les puissances d'une matrice stochastique restent stochastiques : les colonnes de PnP^{n} somment encore à 11.
les trois programmes convergent vers un état stablen = 0n = 1n = 2limite
Les trois barres se rapprochent d'étape en étape et finissent à la répartition limite : c'est le point fixe de la matrice, atteint quel que soit le départ.

Toute matrice stochastique admet 11 comme valeur propre, et l'état stationnaire est justement le vecteur propre associé.

Leontief et Leslie, deux matrices qui se lisent autrement

  • LEONTIEF : x=Ax+dx=Ax+d, donc (IA)x=d\left(I-A\right)x=d et x=(IA)1dx=\left(I-A\right)^{-1}d. La COLONNE jj de AA dit ce que consomme UNE UNITÉ produite par le secteur jj.
  • Chaque entrée de (IA)1\left(I-A\right)^{-1} est la production totale d'un secteur exigée par une unité de demande finale d'un autre; la somme d'une COLONNE est le MULTIPLICATEUR de ce secteur.
  • (IA)1=I+A+A2+\left(I-A\right)^{-1}=I+A+A^{2}+\cdots : effet direct, puis effets indirects de tous les tours.
  • LESLIE : première LIGNE de fécondités, SOUS-DIAGONALE de survies, tout le reste nul. Les premières périodes oscillent, c'est l'onde démographique.

Inverse d'une matrice 2×22\times 2 : échanger la diagonale, changer le signe des deux autres, diviser par le déterminant.

Les pièges qui coûtent des points

Les erreurs ci-dessous sont celles que je corrige le plus souvent en séance. Chacune coûte des points sur une copie, même quand le raisonnement est juste.

1. Confondre la convention ligne et la convention colonne

toute la question, chaque résultat étant transposé

Ce qu'il ne faut pas écrire

« La ligne AA somme à 11, donc j'écris xn+1=Pxnx_{n+1}=Px_{n}. »

Ce qu'il faut écrire

« Les LIGNES somment à 11 : je suis en convention ligne, donc soit j'écris xn+1=xnPx_{n+1}=x_{n}P avec un vecteur ligne, soit je TRANSPOSE la matrice et je garde PxPx. »

Pourquoi : Les deux conventions se lisent exactement de la même façon sur le papier. La seule chose qui les distingue est le sens de sommation, et c'est pour cela qu'on la vérifie AVANT tout calcul.

2. Écrire le produit dans le mauvais ordre

toute la question, et souvent une erreur de dimension visible

Ce qu'il ne faut pas écrire

« x1=x0Px_{1}=x_{0}P, avec x0x_{0} en colonne. »

Ce qu'il faut écrire

« Un vecteur COLONNE se multiplie à DROITE de la matrice : x1=Px0x_{1}=Px_{0}. L'autre ordre n'est même pas défini pour ces dimensions. »

Pourquoi : Le produit matriciel n'est pas commutatif, et les dimensions servent de garde-fou : une matrice 3×33\times 3 suivie d'un vecteur 3×13\times 1 donne bien un vecteur 3×13\times 1, l'inverse ne se calcule pas.

3. Élever chaque entrée à la puissance n

2 points, et un résultat qui n'est même plus stochastique

Ce qu'il ne faut pas écrire

« P2P^{2} : je mets chaque probabilité au carré. »

Ce qu'il faut écrire

« P2=P×PP^{2}=P\times P, produit matriciel : l'entrée (i,j)(i,j) y additionne tous les chemins de jj à ii en DEUX étapes. »

Pourquoi : Élever les entrées au carré détruit la somme des colonnes : 0,72+0,22+0,12=0,540{,}7^{2}+0{,}2^{2}+0{,}1^{2}=0{,}54, ce qui n'est plus une distribution de probabilités.

4. Chercher l'état stationnaire en résolvant le mauvais système

toute la question

Ce qu'il ne faut pas écrire

« L'état stationnaire vérifie Px=0Px=0. »

Ce qu'il faut écrire

« Il vérifie Px=xPx=x, donc (PI)x=0\left(P-I\right)x=0. Résoudre Px=0Px=0 donnerait le noyau de PP, qui est presque toujours réduit au vecteur nul. »

Pourquoi : Un état stationnaire est un point FIXE, pas un vecteur envoyé sur zéro. Retrancher l'identité est précisément ce qui transforme « PP ne change rien » en un système homogène.

5. Oublier de normaliser l'état stationnaire

1,5 point, et une réponse sans unité

Ce qu'il ne faut pas écrire

« Le système donne x=(10 ;11 ;8)x=(10\ ;11\ ;8) : voilà la répartition finale. »

Ce qu'il faut écrire

« Le système est de rang déficient, donc toutes les solutions sont proportionnelles. Avec 10001\,000 étudiants, je normalise : 100029(10 ;11 ;8)(345 ;379 ;276)\dfrac{1000}{29}(10\ ;11\ ;8)\approx(345\ ;379\ ;276). »

Pourquoi : La direction du vecteur propre est unique, sa longueur non. C'est la contrainte de somme, donnée par l'énoncé, qui la fixe.

6. Multiplier par IAI-A au lieu d'inverser

toute la question, et une production largement sous-estimée

Ce qu'il ne faut pas écrire

« (IA)x=d\left(I-A\right)x=d, donc x=(IA)dx=\left(I-A\right)d. »

Ce qu'il faut écrire

« On INVERSE : x=(IA)1dx=\left(I-A\right)^{-1}d. Avec A=(0,20,30,40,1)A=\begin{pmatrix}0{,}2&0{,}3\\0{,}4&0{,}1\end{pmatrix} et d=(100 ;200)d=(100\ ;200), on obtient x=(250 ;333,33)x=(250\ ;333{,}33). »

agricultureindustrieagriculture0,200,30industrie0,400,10la COLONNE dit ce que consomme une unité
La colonne « industrie » se lit ainsi : produire une unité industrielle consomme 0,300{,}30 d'agriculture et 0,100{,}10 d'industrie.

Pourquoi : (IA)1=I+A+A2+\left(I-A\right)^{-1}=I+A+A^{2}+\cdots contient tous les tours de consommation intermédiaire. Multiplier par (IA)\left(I-A\right) ne garde que l'effet direct, à l'envers.

7. Lire la matrice de Leontief en lignes

1,5 point, et une interprétation économique inversée

Ce qu'il ne faut pas écrire

« La première ligne donne ce que l'agriculture consomme : 0,200{,}20 d'elle-même et 0,300{,}30 d'industrie. »

Ce qu'il faut écrire

« La première ligne dit ce que l'agriculture FOURNIT aux autres. Ce qu'un secteur CONSOMME se lit dans sa COLONNE. »

Pourquoi : La convention est la même que pour Markov : l'indice de colonne désigne l'origine, celui de ligne la destination. La matrice se lit toujours « de la colonne vers la ligne ».

8. Prendre les premières périodes de Leslie pour la structure stable

2 points, et la conclusion démographique

Ce qu'il ne faut pas écrire

« Après une période, la classe des jeunes explose : c'est la structure d'équilibre de la population. »

Ce qu'il faut écrire

« Ce sont les oscillations de l'ONDE DÉMOGRAPHIQUE. La structure stable est le vecteur propre de la valeur propre DOMINANTE, atteint seulement après plusieurs périodes. »

Pourquoi : La seconde valeur propre gouverne la vitesse de convergence : l'écart à la limite est multiplié par λ2\left|\lambda_{2}\right| à chaque étape, ce qui prend d'autant plus de temps que ce module est proche de 11.

Quelle méthode choisir

Quel modèle l'énoncé décrit-il ?

On regarde ce que la matrice contient et ce que le vecteur représente.

  • Si les colonnes somment à 11 et les entrées sont des probabilités de passage CHAÎNE DE MARKOV : xn+1=Pxnx_{n+1}=Px_{n}, état stationnaire par (PI)x=0\left(P-I\right)x=0

    Exemple : la mobilité entre trois programmes

  • Si les entrées disent ce qu'un secteur consomme pour produire une unité, et une demande finale est donnée LEONTIEF : x=(IA)1dx=\left(I-A\right)^{-1}d

    Exemple : un tableau entrées-sorties à deux secteurs

  • Si la première ligne porte des fécondités et la sous-diagonale des taux de survie LESLIE : croissance donnée par la valeur propre dominante

    Exemple : une population par classes d'âge

  • Si la matrice n'a aucune de ces structures c'est une application linéaire ordinaire : on revient aux outils généraux du cours

    Exemple : une matrice de rotation

Le premier geste est toujours le même : additionner une colonne. Une somme de 11 annonce Markov, une somme inférieure à 11 annonce Leontief.

Que demande la question, une étape ou la limite ?

On repère l'HORIZON de la question dans l'énoncé.

  • Si « après une année », « à la prochaine étape » un seul produit Px0Px_{0}

    Exemple : de (600 ;300 ;100)(600\ ;300\ ;100) à (490 ;330 ;180)(490\ ;330\ ;180)

  • Si « après nn étapes », avec nn petit itérer le produit, ou calculer PnP^{n} puis multiplier

    Exemple : deux ou trois itérations à la main

  • Si « à long terme », « à l'équilibre », « la répartition finale » état STATIONNAIRE : résoudre (PI)x=0\left(P-I\right)x=0 puis normaliser

    Exemple : (345 ;379 ;276)(345\ ;379\ ;276) sur mille étudiants

  • Si « quel taux de croissance à long terme » VALEUR PROPRE dominante, et son vecteur propre donne la structure stable

    Exemple : une matrice de Leslie

Si l'énoncé demande la limite, itérer le produit ne rapporte presque rien : le correcteur attend le système, pas dix multiplications.

La rédaction attendue

Le correcteur coche des étapes. Les voici dans l'ordre, avec la phrase de conclusion qu'il attend mot pour mot.

Trouver et rédiger un état stationnaire

Quand l'utiliser : Dès que l'énoncé demande la répartition à long terme d'une chaîne de Markov.

  1. 1 Vérifier et ANNONCER la convention : les colonnes somment à 11, donc xn+1=Pxnx_{n+1}=Px_{n}.
  2. 2 Écrire le système (PI)x=0\left(P-I\right)x=0 en toutes lettres, ligne par ligne.
  3. 3 Résoudre, en signalant que le système est de rang déficient et admet donc une infinité de solutions proportionnelles.
  4. 4 NORMALISER selon la contrainte de l'énoncé, effectif total ou somme égale à 11.
  5. 5 Vérifier en calculant PxPx : on doit retrouver xx.

Phrase de conclusion

« Le système (PI)x=0\left(P-I\right)x=0 donne les solutions proportionnelles à (10 ;11 ;8)(10\ ;11\ ;8); normalisées sur 10001\,000 étudiants, elles valent (345 ;379 ;276)(345\ ;379\ ;276), et l'on vérifie que Px=xPx=x. »

Le piège : S'arrêter au vecteur proportionnel : sans normalisation, la réponse ne s'exprime dans aucune unité de l'énoncé.

Barème : 1 point pour la convention, 2 points pour le système, 1 point pour la normalisation, 1 point pour la vérification.

Résoudre un problème de Leontief

Quand l'utiliser : Quand l'énoncé donne une matrice technologique et une demande finale.

  1. 1 Écrire AA et rappeler que sa COLONNE jj décrit la consommation d'une unité du secteur jj.
  2. 2 Former IAI-A et calculer son déterminant, en vérifiant qu'il n'est pas nul.
  3. 3 Inverser, par la formule du 2×22\times 2 ou par la méthode de Gauss-Jordan.
  4. 4 Calculer x=(IA)1dx=\left(I-A\right)^{-1}d et vérifier que Ax+d=xAx+d=x.
  5. 5 Interpréter les MULTIPLICATEURS, c'est-à-dire les sommes des colonnes de l'inverse.

Phrase de conclusion

« La production requise est x=(250 ;333,33)x=(250\ ;333{,}33), et l'on vérifie Ax+d=xAx+d=x. La somme de la première colonne de (IA)1\left(I-A\right)^{-1}, soit 2,172{,}17, est le multiplicateur de l'agriculture : chaque dollar de demande finale agricole mobilise 2,172{,}17 dollars de production totale. »

Le piège : Oublier la vérification Ax+d=xAx+d=x : c'est elle qui attrape une inversion fautive, et elle prend trente secondes.

Barème : 1 point pour IAI-A, 2 points pour l'inversion, 1 point pour le produit, 1 point pour la vérification, 1 point pour l'interprétation.

Vérifier avant de rendre

Cinq minutes de vérification récupèrent plus de points qu'un exercice de plus commencé à la hâte.

L'exercice type décortiqué

La mobilité entre trois programmes, de la première année à la limite

Chaque année, les étudiants changent de programme selon la matrice de transition, en convention colonne : P=(0,70,20,10,20,60,30,10,20,6)P=\begin{pmatrix}0{,}7&0{,}2&0{,}1\\0{,}2&0{,}6&0{,}3\\0{,}1&0{,}2&0{,}6\end{pmatrix}.

Au départ, les effectifs sont x0=(600 ;300 ;100)x_{0}=(600\ ;300\ ;100), soit 10001\,000 étudiants.

On demande la répartition après un an, après deux ans, et la répartition à long terme.

Étape 1

Contrôle de convention : 0,7+0,2+0,1=10{,}7+0{,}2+0{,}1=1, 0,2+0,6+0,2=10{,}2+0{,}6+0{,}2=1, 0,1+0,3+0,6=10{,}1+0{,}3+0{,}6=1. Les COLONNES somment à 11, donc xn+1=Pxnx_{n+1}=Px_{n}.

Pourquoi

Cette vérification tient en une ligne et décide de tout le reste : en convention ligne, chaque produit du problème serait transposé.

Étape 2

x1=Px0=(0,7(600)+0,2(300)+0,1(100) ;)=(490 ;330 ;180)x_{1}=Px_{0}=(0{,}7(600)+0{,}2(300)+0{,}1(100)\ ;\ldots)=(490\ ;330\ ;180), de somme 10001\,000.

Pourquoi

Chaque composante mélange les trois provenances : 490490 étudiants sont en programme 1 parce que 420420 y sont restés, 6060 sont venus du 2 et 1010 du 3.

Étape 3

x2=Px1=(427 ;350 ;223)x_{2}=Px_{1}=(427\ ;350\ ;223), de somme 10001\,000.

Pourquoi

Le premier programme continue de perdre et le troisième de gagner, mais de moins en moins vite : la chaîne s'approche de son point fixe.

Étape 4

État stationnaire : (PI)x=0\left(P-I\right)x=0 donne 0,3x1+0,2x2+0,1x3=0-0{,}3x_{1}+0{,}2x_{2}+0{,}1x_{3}=0 et 0,1x1+0,2x20,4x3=00{,}1x_{1}+0{,}2x_{2}-0{,}4x_{3}=0, d'où x(10 ;11 ;8)x\propto(10\ ;11\ ;8).

Pourquoi

Le système est de rang 22 pour trois inconnues : il ne peut donner qu'une DIRECTION, jamais des effectifs. C'est normal et il faut l'écrire.

Étape 5

Normalisation sur 10001\,000 étudiants : 100029(10 ;11 ;8)(344,8 ;379,3 ;275,9)\dfrac{1000}{29}(10\ ;11\ ;8)\approx(344{,}8\ ;379{,}3\ ;275{,}9).

Pourquoi

La contrainte de somme vient de l'énoncé, pas de la matrice. Sans elle, la réponse resterait un rapport entre trois nombres.

Étape 6

Vérification : 0,7(344,8)+0,2(379,3)+0,1(275,9)344,80{,}7(344{,}8)+0{,}2(379{,}3)+0{,}1(275{,}9)\approx 344{,}8.

Pourquoi

Retrouver la première composante confirme d'un coup la résolution du système et la normalisation, sans refaire aucun calcul.

Conclusion rédigée

« Après un an la répartition est (490 ;330 ;180)(490\ ;330\ ;180), après deux ans (427 ;350 ;223)(427\ ;350\ ;223), et la chaîne étant régulière elle converge vers (345 ;379 ;276)(345\ ;379\ ;276), quelle qu'ait été la répartition de départ. »

L'erreur classique sur cet exercice : Itérer le produit dix fois pour approcher la limite : le correcteur attend le système (PI)x=0\left(P-I\right)x=0, et l'itération ne rapporte presque aucun point.

À savoir par cœur

  • Convention COLONNE : pijp_{ij} va de jj VERS ii, et chaque colonne somme à 11. On l'annonce avant tout calcul.
  • Une étape, xn+1=Pxnx_{n+1}=Px_{n}; nn étapes, xn=Pnx0x_{n}=P^{n}x_{0}, et PnP^{n} est un PRODUIT matriciel, pas des puissances d'entrées.
  • État stationnaire : résoudre (PI)x=0\left(P-I\right)x=0, jamais Px=0Px=0, puis NORMALISER.
  • LEONTIEF : x=(IA)1dx=\left(I-A\right)^{-1}d, et la somme d'une colonne de l'inverse est le MULTIPLICATEUR du secteur.
  • La colonne d'une matrice technologique dit ce qu'un secteur CONSOMME; sa ligne dit ce qu'il FOURNIT.
  • LESLIE : première ligne de fécondités, sous-diagonale de survies, et les premières périodes OSCILLENT.
  • La valeur propre DOMINANTE donne le taux de croissance à long terme, son vecteur propre la structure stable.

Questions fréquentes

Comment savoir si une matrice de transition est en convention colonne ?

On additionne une colonne : si le total vaut un, la convention est celle des colonnes et le vecteur d'état se multiplie à droite de la matrice. Si ce sont les lignes qui somment à un, la matrice est transposée par rapport à cette convention, et il faut soit la transposer, soit travailler avec un vecteur ligne multiplié à gauche.

Comment trouver l'état stationnaire d'une chaîne de Markov ?

On résout le système obtenu en retranchant l'identité à la matrice de transition et en égalant à zéro. Ce système est de rang déficient, donc il donne une famille de solutions proportionnelles : il faut ensuite normaliser selon la contrainte de l'énoncé, effectif total ou somme égale à un. On vérifie enfin que la matrice laisse ce vecteur inchangé.

Que représente l'inverse de la matrice de Leontief ?

Chacune de ses entrées donne la production totale qu'un secteur doit fournir pour satisfaire une unité de demande finale d'un autre secteur, effets indirects compris. La somme d'une colonne est le multiplicateur de ce secteur : elle dit combien de production totale mobilise chaque dollar de demande finale qui lui est adressé.

Pourquoi une population de Leslie oscille-t-elle au début ?

Parce que la structure de départ n'est pas celle d'équilibre. Une classe d'âge trop nombreuse produit une génération nombreuse une période plus tard, qui en produit une autre encore, et ce mouvement met plusieurs périodes à s'amortir. C'est l'onde démographique, et sa vitesse d'extinction dépend de la deuxième valeur propre.

Peut-on élever une matrice à une puissance entrée par entrée ?

Non. La puissance d'une matrice est un produit matriciel répété, et c'est justement ce qui lui donne son sens : l'entrée d'une puissance deuxième additionne automatiquement tous les chemins en deux étapes. Élever les entrées séparément détruit d'ailleurs la somme des colonnes, qui ne vaut plus un.

Passer à la pratique

Exercices corrigés : Modèles matriciels : Markov, Leontief, Leslie

Une méthode se prouve sur une copie, pas sur une fiche. La série du même chapitre reprend chacun de ces pièges dans un exercice, avec le corrigé rédigé étape par étape.

  • 15 exercices corrigés
  • 150 points
  • 225 minutes
Faire les exercices
Fiche précédente Diagonalisation et valeurs propres Fiche suivante Nombres complexes et De Moivre

Voir aussi

Vous cherchez un tuteur en algèbre linéaire à Montréal ?

Contactez-moi pour une première séance. Les modèles matriciels sont la partie du cours où l'algèbre linéaire cesse d'être abstraite, et c'est souvent là que le déclic se produit.

Site par Studio Squalli