Algèbre linéaire

Vecteurs, matrices et systèmes d’équations : élimination, LU, inverses, modèles entrées–sorties, les quatre sous-espaces fondamentaux, projections et moindres carrés, déterminants, valeurs propres, chaînes de Markov et SVD, avec six TP. Gratuit, à votre rythme, et lié au GPS, à Landsat, aux satellites de niveau de la mer, aux modèles météo, au recyclage et à la cartographie minérale.

NASA/Kim Shiflett (KSC-20210927-PH-KLS04_0013, Sept. 27, 2021)NASA/JPL-Caltech (PIA22431)CNES (PIA25623)NASA Blue Marble

Cours › MATH 230

Télécharger le guide d’étude en PDF (anglais)

À propos de ce cours

Ceci est le plan d’un cours gratuit en autonomie rédigé par Saving Planets. Ce n’est pas un cours universitaire : il ne donne aucun crédit académique et aucun établissement n’y participe.

Son contenu suit un cours d’algèbre linéaire d’un semestre de deuxième année. Il s’inspire des thèmes et du niveau de MIT 18.06 Linear Algebra enseigné par Gilbert Strang (gratuit sur MIT OpenCourseWare) et des cours d’algèbre linéaire de deuxième année de Georgia Tech et de Purdue, en s’appuyant sur le manuel gratuit de Georgia Tech Interactive Linear Algebra de Dan Margalit et Joseph Rabinoff et sur le manuel gratuit Linear Algebra de Jim Hefferon, que nous citons quand ils sont utiles. Aucun de ces auteurs ou établissements ne cautionne cette page.

Ce que nous ajoutons : chaque unité est reliée à un vrai travail sur la planète avec une source publique citée (GPS, Landsat 9, satellites de niveau de la mer de la NASA, supercalculateurs météo de la NOAA, données entrées–sorties et recyclage de l’EPA, études de bâtiments de l’USGS, spectromètre minéral EMIT de la NASA), chaque nombre de la page est recalculé par un script de vérification avant sa mise en ligne, et chaque TP se fait dans un tableur gratuit ou sur papier.

Suivez-le dans l’ordre. Chaque unité comporte des objectifs, des notions clés, des équations, des exemples pas à pas, des exercices corrigés en détail, des erreurs fréquentes, des cartes mémoire et un quiz. Votre progression est enregistrée uniquement dans ce navigateur.

Comment ce cours aide à sauver la planète

L’algèbre linéaire est le langage des données sur la planète. Un pixel satellite est un vecteur, une position GPS est la solution d’un système linéaire, une prévision météo, ce sont des millions d’équations résolues encore et encore, et la pollution d’une économie peut se suivre avec une seule matrice.

Les mêmes quelques idées, élimination, projection, valeurs propres et SVD, transforment des mesures brutes en cartes de forêts, en tendances du niveau de la mer et en cartes de poussières minérales. Chaque unité se termine par une courte section sourcée. Les valeurs qui ne proviennent pas d’une source sont indiquées comme valeurs d’exemple.

Haut de page

Avant de commencer

Prérequis

  • Algèbre du lycée : résoudre deux équations à deux inconnues et manipuler des fonctions.
  • Analyse I aide pour les limites et les dérivées dans quelques exemples, mais aucune analyse n’est nécessaire pour la plus grande partie du cours.
  • Être à l’aise avec un tableur ; Introduction à la programmation aide pour les versions Python facultatives des TP.
  • Analyse III et Équations différentielles vont bien avec ce cours mais ne sont pas requis.

Comment étudier

  • Faites d’abord chaque calcul à la main sur de petits exemples (2 × 2 et 3 × 3), puis vérifiez-le avec un tableur ou une calculatrice.
  • Pensez à Ax de deux façons : comme les produits scalaires des lignes avec x, et comme une combinaison des colonnes de A.
  • Vérifiez chaque réponse : remplacez x dans Ax = b, multipliez A par son inverse, ou testez Av = λv.
  • Surveillez les dimensions : une matrice m × n fois une matrice n × p donne m × p, et des dimensions qui ne concordent pas signalent une erreur.
  • Faites des dessins en 2D et en 3D : droites, plans, projections et étirements rendent les règles faciles à retenir.
  • Faites chaque exercice sur papier avant d’ouvrir la réponse, et refaites les quiz quelques jours plus tard.

Matériel de TP (pour les TP en laboratoire)

  • Un ordinateur avec un tableur gratuit (ses fonctions MMULT, MINVERSE, MDETERM et LINEST font le calcul matriciel).
  • Du papier millimétré, une règle, un rapporteur et un crayon.
  • Facultatif : Python gratuit avec NumPy pour les versions plus grandes des TP.
  • Une ficelle et des punaises ou du ruban adhésif pour le TP de trilatération, et une carte imprimée ou une feuille de papier.
  • Une photo en niveaux de gris ou une petite grille de nombres pour le TP de compression d’image.

Haut de page

Programme suggéré sur 15 semaines

Environ 6 à 9 heures par semaine : lecture, exercices et, les semaines de TP, le TP.

SemaineUnitéThèmesLecture
1Vecteurs : combinaisons linéaires, normes et produits scalairesVecteurs, combinaisons linéaires, espace engendré, normeMIT 18.06 L1; ILA 2.1–2.2
2Vecteurs : combinaisons linéaires, normes et produits scalairesProduits scalaires, angles, orthogonalité ; TP A ; début de l’éliminationILA 6.1
3Systèmes d’équations linéaires et élimination de GaussÉlimination, formes échelonnées, ensembles de solutions ; TP BMIT 18.06 L2; ILA 1.1–1.3
4Algèbre matricielle : multiplication, inverses et factorisation LUMultiplication matricielle, inverses, A = LU ; TP CMIT 18.06 L3–4; ILA 3.4–3.5
5Modèles linéaires : économie entrées–sorties (Leontief), réseaux et équilibrageModèles entrées–sorties, réseaux, équilibrageHefferon Topics; EPA USEEIO
6Espaces vectoriels, sous-espaces, indépendance, base, rang et les quatre sous-espaces fondamentauxSous-espaces, espace des colonnes, noyau, Ax = 0MIT 18.06 L5–7; ILA 2.6
7Espaces vectoriels, sous-espaces, indépendance, base, rang et les quatre sous-espaces fondamentauxIndépendance, base, dimension, rang, quatre sous-espaces ; entraînement de mi-parcoursMIT 18.06 L8–10; ILA 2.5, 2.7, 2.9
8Orthogonalité, projections, moindres carrés et Gram–SchmidtOrthogonalité et projectionsMIT 18.06 L14–15; ILA 6.3
9Orthogonalité, projections, moindres carrés et Gram–SchmidtMoindres carrés, Gram–Schmidt, A = QR ; TP DMIT 18.06 L16–17; ILA 6.4–6.5
10Déterminants : propriétés, cofacteurs, règle de Cramer et volumeDéterminants, cofacteurs, règle de Cramer, volumeMIT 18.06 L18–20; ILA 4.1–4.3
11Valeurs propres et vecteurs propresValeurs propres, vecteurs propres, méthode de la puissance itérée, modesMIT 18.06 L21; ILA 5.1–5.2, 5.5
12Diagonalisation, puissances de matrices, chaînes de Markov et systèmes dynamiques linéairesDiagonalisation, chaînes de Markov, systèmes dynamiques ; TP EMIT 18.06 L22–25; ILA 5.4, 5.6
13La décomposition en valeurs singulières, l’analyse en composantes principales et les matrices définies positivesMatrices définies positives et SVDMIT 18.06 L27–29
14La décomposition en valeurs singulières, l’analyse en composantes principales et les matrices définies positivesApproximation de rang faible, ACP, conditionnement ; TP FMIT 18.06 L29 ; site du livre de Strang
15RévisionRévision et série d’exercices finaleToutes les unités

Haut de page

Votre progression

Cochez les unités et les TP à mesure que vous les terminez, et faites les quiz. La progression et les scores sont enregistrés uniquement dans ce navigateur et ne nous sont jamais envoyés.

Activez JavaScript pour suivre votre progression.

Haut de page

Unité 1 · Semaines 1–2

Vecteurs : combinaisons linéaires, normes et produits scalaires

Vecteurs en deux, trois et n dimensions, addition et multiplication par un scalaire, combinaisons linéaires et espace engendré, norme, vecteurs unitaires, produit scalaire, angle entre vecteurs, orthogonalité et inégalités de Cauchy–Schwarz et triangulaire.

Objectifs d’apprentissage

À la fin de cette unité, vous saurez :

  • Additionner, multiplier et combiner des vecteurs à la main et en dessin.
  • Calculer des normes, des vecteurs unitaires et des produits scalaires.
  • Trouver l’angle entre deux vecteurs et tester l’orthogonalité.
  • Décrire l’espace engendré par un, deux ou trois vecteurs.

Notions clés

  • Un vecteur de ℝⁿ est une liste ordonnée de n nombres. En 2D et 3D c’est une flèche ; en dimension plus grande c’est une liste, comme les 11 valeurs de bandes d’un pixel satellite.
  • Une combinaison linéaire c₁v₁ + c₂v₂ + … multiplie des vecteurs par des scalaires et les additionne. L’espace engendré est l’ensemble de toutes les combinaisons linéaires : une droite, un plan ou tout l’espace.
  • Norme ‖v‖ = √(v·v). Diviser par la norme donne un vecteur unitaire u = v/‖v‖ de norme 1.
  • Produit scalaire v·w = v₁w₁ + v₂w₂ + … = ‖v‖‖w‖cos θ. Si v·w = 0, les vecteurs sont orthogonaux (perpendiculaires).
  • Cauchy–Schwarz : |v·w| ≤ ‖v‖‖w‖. Inégalité triangulaire : ‖v + w‖ ≤ ‖v‖ + ‖w‖.

Équations clés

∥v∥=v12+v22+⋯+vn2\|v\|=\sqrt{v_1^2+v_2^2+\cdots+v_n^2}
Norme
v⋅w=v1w1+⋯+vnwnv\cdot w=v_1w_1+\cdots+v_nw_n
Produit scalaire
cos⁡θ=v⋅w∥v∥ ∥w∥\cos\theta=\frac{v\cdot w}{\|v\|\,\|w\|}
Angle entre vecteurs
∣v⋅w∣≤∥v∥ ∥w∥|v\cdot w|\le\|v\|\,\|w\|
Inégalité de Cauchy–Schwarz

Exemples résolus pas à pas

Exemple 1.1 Trouvez la norme de v = (3, 4) et le vecteur unitaire de même direction.

Afficher les étapes
  1. ‖v‖ = √(3² + 4²) = √25 = 5.
  2. u = v/5.

Réponse : ‖v‖ = 5, u = (0.6, 0.8).

Exemple 1.2 Trouvez l’angle entre v = (3, 4) et w = (1, 2).

Afficher les étapes
  1. v·w = 3 + 8 = 11.
  2. ‖v‖ = 5, ‖w‖ = √5 = 2.236.
  3. cos θ = 11/(5 × 2.236) = 0.9839.

Réponse : θ = 10.30°.

Exemple 1.3 Écrivez 2(1, 0, 1) + 3(0, 1, 1) sous forme d’un seul vecteur. (2, 3, 4) appartient-il à l’espace engendré par ces deux vecteurs ?

Afficher les étapes
  1. 2(1, 0, 1) + 3(0, 1, 1) = (2, 3, 5).
  2. Toute combinaison c₁(1, 0, 1) + c₂(0, 1, 1) = (c₁, c₂, c₁ + c₂), donc la troisième composante doit égaler la somme des deux premières.

Réponse : (2, 3, 5) ; (2, 3, 4) n’est pas dans l’espace engendré car 2 + 3 ≠ 4.

Exercices

Faites chacun d’abord sur papier, puis vérifiez.

  1. Trouvez ‖(1, 2, 2)‖.

    Afficher la réponse

    √(1 + 4 + 4) = 3.

  2. Trouvez (1, 2, 3)·(4, −5, 6).

    Afficher la réponse

    4 − 10 + 18 = 12.

  3. (2, −1) et (1, 2) sont-ils orthogonaux ?

    Afficher la réponse

    Yes: 2 − 2 = 0.

  4. Trouvez le vecteur unitaire de même direction que (6, 8).

    Afficher la réponse

    ‖(6, 8)‖ = 10, donc (0.6, 0.8).

Idées reçues fréquentes

  • Idée reçue : Le produit scalaire est un vecteur.

    En réalité : C’est un seul nombre (un scalaire).

  • Idée reçue : ‖v + w‖ = ‖v‖ + ‖w‖ toujours.

    En réalité : Seulement quand v et w pointent dans le même sens ; sinon c’est plus petit.

  • Idée reçue : Les vecteurs n’ont de sens qu’en 2D ou en 3D.

    En réalité : Les vecteurs de données peuvent avoir des centaines de composantes ; les mêmes règles s’appliquent.

Cartes mémoire

  • Combinaison linéaire

    c₁v₁ + c₂v₂ + …

  • Espace engendré

    Toutes les combinaisons linéaires d’un ensemble de vecteurs.

  • Norme

    ‖v‖ = √(v·v).

  • Produit scalaire

    v·w = Σ vᵢwᵢ = ‖v‖‖w‖cos θ.

  • Orthogonal

    v·w = 0.

  • Vecteur unitaire

    v/‖v‖.

Quiz d’autocontrôle

1. ‖(3, 4)‖ vaut…

Afficher la réponse

5. √(9 + 16) = 5.

2. v·w = 0 signifie que les vecteurs sont…

Afficher la réponse

orthogonaux. cos θ = 0, donc θ = 90°.

3. L’espace engendré par deux vecteurs non parallèles de ℝ³ est…

Afficher la réponse

un plan. Deux directions indépendantes remplissent un plan.

4. Doubler v change ‖v‖ d’un facteur…

Afficher la réponse

2. ‖cv‖ = |c|‖v‖.

Comment cela aide à sauver la planète

Chaque pixel de Landsat 9 est un vecteur. Selon l’USGS, l’Operational Land Imager 2 de Landsat 9 enregistre 9 bandes spectrales à 30 m (15 m pour la bande panchromatique) et son Thermal Infrared Sensor 2 enregistre 2 bandes thermiques à 100 m, si bien que chaque point au sol devient une liste de 11 nombres. Les données sont sur 14 bits, 2¹⁴ = 16384 niveaux par bande, et Landsat 9 acquiert environ 750 scènes par jour, toutes gratuites.

La santé des forêts et des cultures se lit dans la bande rouge (bande 4, 0.64–0.67 µm) et la bande proche infrarouge (bande 5, 0.85–0.88 µm). L’indice de végétation NDVI = (NIR − red)/(NIR + red) compare deux composantes du vecteur pixel. Avec des réflectances d’exemple de 0.40 (NIR) et 0.08 (rouge), NDVI = 0.32/0.48 = 0.667, typique d’une végétation en bonne santé ; un sol nu donne des valeurs proches de zéro.

Sur Saving Planets : Observation de la Terre · Déforestation

Sources : USGS: Landsat 9 · Interactive Linear Algebra, 2.1 Vectors

Ressources d’apprentissage gratuites

Haut de page

Unité 2 · Semaines 2–3

Systèmes d’équations linéaires et élimination de Gauss

Systèmes linéaires vus par lignes et par colonnes, matrices augmentées, élimination jusqu’à la forme échelonnée, remontée, forme échelonnée réduite, pivots et variables libres, et quand un système a une, aucune ou une infinité de solutions.

Objectifs d’apprentissage

À la fin de cette unité, vous saurez :

  • Écrire un système sous la forme Ax = b et comme matrice augmentée.
  • Résoudre des systèmes par élimination et remontée.
  • Trouver la forme échelonnée réduite et identifier pivots et variables libres.
  • Décider si un système a une, aucune ou une infinité de solutions.

Notions clés

  • Vue par lignes : chaque équation est une droite (2D) ou un plan (3D) et la solution est leur intersection. Vue par colonnes : b est une combinaison des colonnes de A, avec les poids x.
  • Trois opérations sur les lignes ne changent pas les solutions : échanger des lignes, multiplier une ligne par un nombre non nul, et ajouter à une ligne un multiple d’une autre.
  • L’élimination produit une forme triangulaire supérieure (échelonnée) ; la remontée résout ensuite de bas en haut. Les pivots sont les premiers coefficients non nuls de chaque ligne.
  • Une colonne sans pivot donne une variable libre et une infinité de solutions. Une ligne 0 = c avec c ≠ 0 signifie aucune solution (système incompatible).
  • Des équations de distance au carré comme (x − a)² + (y − b)² = d² ne sont pas linéaires, mais en soustraire une d’une autre élimine x² et y² et laisse des équations linéaires. C’est ainsi que les problèmes de positionnement deviennent des systèmes linéaires.

Équations clés

Ax=bA\mathbf{x}=\mathbf{b}
Forme matricielle
Ri←Ri−ℓijRj,  ℓij=aij/ajjR_i\leftarrow R_i-\ell_{ij}R_j,\ \ \ell_{ij}=a_{ij}/a_{jj}
Étape d’élimination et multiplicateur
solutions: one, none, or infinitely many\text{solutions: one, none, or infinitely many}
Les trois seuls cas
d=c td=c\,t
Distance à partir du temps de trajet du signal (GPS)

Exemples résolus pas à pas

Exemple 2.1 Résolvez x + y + z = 6, 2x + y − z = 1, x − y + 2z = 5.

Afficher les étapes
  1. R₂ − 2R₁: −y − 3z = −11. R₃ − R₁: −2y + z = −1.
  2. R₃ − 2R₂ : 7z = 21, donc z = 3.
  3. Remontée : −y − 9 = −11 donne y = 2 ; puis x = 6 − 2 − 3 = 1.

Réponse : (x, y, z) = (1, 2, 3).

Exemple 2.2 Des balises sont en (0, 0), (10, 0) et (0, 10) km. Un récepteur mesure des distances au carré de 25, 65 et 45 km². Trouvez sa position.

Afficher les étapes
  1. Équation 1 moins équation 2 : x² + y² − [(x − 10)² + y²] = 25 − 65, donc 20x − 100 = −40.
  2. Équation 1 moins équation 3 : 20y − 100 = −20.
  3. Résolvez les deux équations linéaires.

Réponse : (x, y) = (3, 4) km.

Exemple 2.3 Combien de solutions ont x + 2y = 3, 2x + 4y = 6 ? Et x + 2y = 3, 2x + 4y = 7 ?

Afficher les étapes
  1. R₂ − 2R₁ donne 0 = 0 dans le premier système : un pivot, une variable libre.
  2. Dans le second, R₂ − 2R₁ donne 0 = 1.

Réponse : Une infinité ; aucune.

Exercices

Faites chacun d’abord sur papier, puis vérifiez.

  1. Résolvez x + y = 5, x − y = 1.

    Afficher la réponse

    En additionnant, 2x = 6 : x = 3, y = 2.

  2. Un signal de type GPS met 0.00505783 s à arriver. Quelle distance a-t-il parcourue à 299792458 m/s ?

    Afficher la réponse

    d = ct = 1516299 m, environ 1516 km.

  3. Trouvez la forme échelonnée réduite de [[1, 2, 3], [2, 5, 8]] (augmentée).

    Afficher la réponse

    R₂ − 2R₁ = [0, 1, 2]; R₁ − 2R₂ = [1, 0, −1]: x = −1, y = 2.

  4. Un système 3 × 5 a 3 pivots. Combien de variables libres ?

    Afficher la réponse

    5 − 3 = 2.

Idées reçues fréquentes

  • Idée reçue : Multiplier une ligne par 0 est une opération sur les lignes permise.

    En réalité : Cela détruit de l’information ; le multiplicateur doit être non nul.

  • Idée reçue : Un système avec plus d’inconnues que d’équations n’a pas de solution.

    En réalité : S’il est compatible, il en a une infinité ; il peut aussi être incompatible.

  • Idée reçue : Une ligne nulle signifie aucune solution.

    En réalité : Seule une ligne 0 = c avec c ≠ 0 le signifie ; 0 = 0 indique juste une variable libre.

Cartes mémoire

  • Matrice augmentée

    [A | b].

  • Pivot

    Premier coefficient non nul d’une ligne en forme échelonnée.

  • Variable libre

    Variable dont la colonne n’a pas de pivot.

  • RREF

    Pivots égaux à 1, zéros au-dessus et en dessous.

  • Incompatible

    Une ligne 0 = c, c ≠ 0.

  • Remontée

    Résoudre le système triangulaire de bas en haut.

Quiz d’autocontrôle

1. Une ligne 0 0 0 | 5 en forme échelonnée signifie…

Afficher la réponse

aucune solution. 0 = 5 est impossible.

2. Laquelle n’est PAS une opération sur les lignes ?

Afficher la réponse

élever une ligne au carré. Élever au carré n’est pas linéaire.

3. Une colonne sans pivot donne…

Afficher la réponse

une variable libre. Sa variable peut être choisie librement.

4. Vu par colonnes, Ax = b demande si b est…

Afficher la réponse

une combinaison des colonnes de A. Ax = x₁a₁ + x₂a₂ + …

Comment cela aide à sauver la planète

Le GPS est une machine d’algèbre linéaire. GPS.gov indique qu’au moins 24 satellites opérationnels sont disponibles 95 % du temps et que 31 satellites opérationnels sont en service depuis plus d’une décennie, à environ 20200 km d’altitude, chacun faisant le tour de la Terre deux fois par jour. Ils sont répartis sur six plans orbitaux également espacés, si bien que les utilisateurs voient au moins quatre satellites depuis pratiquement n’importe quel point de la planète.

L’activité de trilatération de GPS.gov explique l’idée : distance = vitesse × temps, la vitesse étant celle de la lumière, 299792458 m/s, et les distances à au moins 3 satellites fixent une position 3D. L’horloge d’un vrai récepteur est aussi inconnue, donc dans le modèle classique des manuels les inconnues sont x, y, z et le décalage d’horloge, quatre inconnues qui exigent quatre satellites ; soustraire les équations de distance, comme dans l’exemple résolu, en fait un système linéaire. Une erreur d’horloge d’une seule microseconde fausserait une distance de 299.8 m. Les positions précises guident les tracteurs qui épandent moins d’engrais, les navires sur des routes efficaces et les stations de mesure qui suivent le niveau de la mer et la glace.

Sur Saving Planets : Observation de la Terre · Espace

Sources : GPS.gov: Space segment · GPS.gov: Trilateration (NOAA classroom activity)

Ressources d’apprentissage gratuites

Haut de page

Unité 3 · Semaine 4

Algèbre matricielle : multiplication, inverses et factorisation LU

Produits matrice–vecteur et matrice–matrice, règles (et non-règles) de l’algèbre matricielle, transposées, matrices d’élimination, inverse et Gauss–Jordan, formule de l’inverse 2 × 2, A = LU et coût de l’élimination.

Objectifs d’apprentissage

À la fin de cette unité, vous saurez :

  • Multiplier des matrices en vérifiant que les dimensions concordent.
  • Trouver des inverses par la formule 2 × 2 et par Gauss–Jordan.
  • Factoriser A = LU et s’en servir pour résoudre des systèmes.
  • Estimer le coût de l’élimination quand n augmente.

Notions clés

  • (AB)ᵢⱼ = produit scalaire de la ligne i de A par la colonne j de B. Une matrice m × n fois une matrice n × p donne m × p. En général AB ≠ BA.
  • L’inverse A⁻¹ vérifie AA⁻¹ = A⁻¹A = I. Elle existe exactement quand A a n pivots (det A ≠ 0). (AB)⁻¹ = B⁻¹A⁻¹ et (AB)ᵀ = BᵀAᵀ.
  • Gauss–Jordan : réduire [A | I] en [I | A⁻¹].
  • L’élimination sans échange de lignes enregistre ses multiplicateurs dans une matrice triangulaire inférieure L avec des uns sur la diagonale, ce qui donne A = LU. Résolvez Lc = b, puis Ux = c ; avec L et U conservées, chaque nouveau second membre coûte peu.
  • L’élimination sur une matrice n × n coûte environ n³/3 multiplications–additions ; résoudre avec L et U conservées coûte environ n² par second membre.

Équations clés

(AB)ij=∑kaikbkj(AB)_{ij}=\sum_k a_{ik}b_{kj}
Produit matriciel
[abcd]−1=1ad−bc[d−b−ca]\begin{bmatrix}a&b\\c&d\end{bmatrix}^{-1}=\frac{1}{ad-bc}\begin{bmatrix}d&-b\\-c&a\end{bmatrix}
Inverse 2 × 2
A=LUA=LU
Factorisation LU
cost≈n3/3\text{cost}\approx n^3/3
Nombre d’opérations de l’élimination

Exemples résolus pas à pas

Exemple 3.1 Trouvez AB et BA pour A = [[1, 2], [3, 4]] et B = [[0, 1], [1, 0]].

Afficher les étapes
  1. AB : ligne 1 · colonnes = (2, 1) ; ligne 2 = (4, 3).
  2. BA échange les lignes de A : [[3, 4], [1, 2]].

Réponse : AB = [[2, 1], [4, 3]], BA = [[3, 4], [1, 2]], donc AB ≠ BA.

Exemple 3.2 Inversez A = [[4, 7], [2, 6]].

Afficher les étapes
  1. ad − bc = 24 − 14 = 10.
  2. A⁻¹ = (1/10)[[6, −7], [−2, 4]].

Réponse : A⁻¹ = [[0.6, −0.7], [−0.2, 0.4]].

Exemple 3.3 Factorisez A = [[2, 1], [6, 8]] sous la forme LU.

Afficher les étapes
  1. Multiplicateur ℓ₂₁ = 6/2 = 3.
  2. R₂ − 3R₁ = [0, 5].

Réponse : L = [[1, 0], [3, 1]], U = [[2, 1], [0, 5]].

Exercices

Faites chacun d’abord sur papier, puis vérifiez.

  1. Quelle est la taille d’une matrice 3 × 2 fois une matrice 2 × 4 ?

    Afficher la réponse

    3 × 4.

  2. Inversez [[2, 1], [5, 3]].

    Afficher la réponse

    det = 1, inverse [[3, −1], [−5, 2]].

  3. Combien d’opérations environ l’élimination demande-t-elle pour n = 1000 ?

    Afficher la réponse

    n³/3 ≈ 3.33 × 10⁸.

  4. Si A et B sont inversibles, que vaut (AB)⁻¹ ?

    Afficher la réponse

    B⁻¹A⁻¹.

Idées reçues fréquentes

  • Idée reçue : AB = BA.

    En réalité : La multiplication matricielle est en général non commutative.

  • Idée reçue : (AB)⁻¹ = A⁻¹B⁻¹.

    En réalité : L’ordre s’inverse : B⁻¹A⁻¹.

  • Idée reçue : Les ordinateurs résolvent Ax = b en calculant A⁻¹.

    En réalité : Ils factorisent A = LU (ou similaire) ; former A⁻¹ est plus lent et moins précis.

Cartes mémoire

  • Identité

    AI = IA = A.

  • Inverse

    AA⁻¹ = I.

  • Inverse 2 × 2

    Échanger a et d ; changer le signe de b et c ; diviser par ad − bc.

  • Règle de la transposée

    (AB)ᵀ = BᵀAᵀ.

  • LU

    A = LU issue de l’élimination.

  • Coût

    Environ n³/3 pour l’élimination.

Quiz d’autocontrôle

1. Une matrice 2 × 3 fois une matrice 3 × 5 donne…

Afficher la réponse

2 × 5. Les dimensions extérieures.

2. Une matrice est inversible exactement quand elle a…

Afficher la réponse

n pivots. Un jeu complet de pivots.

3. Doubler n multiplie le coût de l’élimination par environ…

Afficher la réponse

8. n³ croît de 2³ = 8 fois.

4. L dans A = LU contient…

Afficher la réponse

les multiplicateurs. Les multiplicateurs sous la diagonale, des uns sur la diagonale.

Comment cela aide à sauver la planète

Les prévisions météo et climatiques sont, au fond, d’énormes calculs d’algèbre linéaire répétés plusieurs fois par jour. La NOAA a annoncé en août 2023 que ses deux supercalculateurs opérationnels jumeaux, Dogwood à Manassas (Virginie) et Cactus à Phoenix (Arizona), avaient été améliorés de 20 % pour atteindre 14.5 pétaflops chacun, si bien qu’ensemble le système de prévision peut effectuer 29 quadrillions de calculs par seconde, et qu’avec ses supercalculateurs de recherche sa capacité totale est de 49 pétaflops.

La règle n³/3 montre pourquoi la vitesse compte. À titre d’exemple, un système dense à n = 100000 inconnues demande environ 3.33 × 10¹⁴ opérations ; à 14.5 × 10¹⁵ opérations par seconde en continu cela fait environ 0.023 s, mais doubler n multiplie le travail par 8. Les vrais modèles utilisent des matrices creuses et des solveurs astucieux, mais c’est cette même croissance cubique que la puissance supplémentaire doit financer selon la NOAA : des modèles à plus haute résolution qui améliorent prévisions et alertes.

Sur Saving Planets : Introduction à la programmation · Solutions

Sources : NOAA: NOAA completes upgrade to weather and climate supercomputer system (Aug. 10, 2023) · MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010)

Ressources d’apprentissage gratuites

Haut de page

Unité 4 · Semaine 5

Modèles linéaires : économie entrées–sorties (Leontief), réseaux et équilibrage

Construire des modèles linéaires à partir de systèmes réels : le modèle entrées–sorties de Leontief x = Ax + d et sa solution x = (I − A)⁻¹d, les modèles étendus à l’environnement qui convertissent les dépenses en émissions, les flux dans les réseaux et l’équilibrage des équations chimiques.

Objectifs d’apprentissage

À la fin de cette unité, vous saurez :

  • Poser et résoudre un modèle entrées–sorties de Leontief.
  • Interpréter les coefficients de (I − A)⁻¹.
  • Étendre un modèle entrées–sorties avec des intensités d’émission pour obtenir une empreinte.
  • Écrire des problèmes d’équilibre de flux et d’équilibre chimique comme des systèmes linéaires.

Notions clés

  • Dans un modèle entrées–sorties, la colonne j de la matrice de consommation A indique ce que le secteur j achète à chaque secteur pour produire un dollar de production.
  • La production totale x doit couvrir ce qu’utilisent les industries (Ax) plus la demande finale d : x = Ax + d, donc (I − A)x = d et x = (I − A)⁻¹d.
  • Le coefficient (i, j) de l’inverse de Leontief (I − A)⁻¹ est la production totale du secteur i nécessaire, directement et via toutes les chaînes d’approvisionnement, pour un dollar de demande finale adressée au secteur j.
  • Les modèles étendus à l’environnement multiplient la production par une ligne d’intensités d’émission e (par exemple kg CO₂ par dollar) : empreinte = e·x.
  • Les lois de conservation donnent aussi des équations linéaires : le flux entrant égale le flux sortant en chaque nœud d’un réseau, et les atomes de chaque élément s’équilibrent dans une réaction chimique.

Équations clés

x=Ax+d\mathbf{x}=A\mathbf{x}+\mathbf{d}
Équilibre entrées–sorties
x=(I−A)−1d\mathbf{x}=(I-A)^{-1}\mathbf{d}
Solution de Leontief
footprint=e⋅x\text{footprint}=\mathbf{e}\cdot\mathbf{x}
Modèle étendu à l’environnement
∑in=∑out\sum \text{in}=\sum \text{out}
Équilibre des flux en un nœud

Exemples résolus pas à pas

Exemple 4.1 Une économie à deux secteurs a A = [[0.2, 0.3], [0.1, 0.4]] et une demande finale d = (100, 50) (valeurs d’exemple, en millions de $). Trouvez la production totale.

Afficher les étapes
  1. I − A = [[0.8, −0.3], [−0.1, 0.6]], det = 0.48 − 0.03 = 0.45.
  2. (I − A)⁻¹ = (1/0.45)[[0.6, 0.3], [0.1, 0.8]].
  3. x = (1/0.45)(60 + 15, 10 + 40).

Réponse : x = (166.7, 111.1) millions de $.

Exemple 4.2 Avec des intensités d’émission d’exemple e = (0.5, 2.0) kg CO₂ par dollar, trouvez l’empreinte de cette demande finale.

Afficher les étapes
  1. empreinte = 0.5 × 166.7 + 2.0 × 111.1 (en kg par $ × millions de $).

Réponse : 305.6 millions de kg = 305.6 kt CO₂ (exemple).

Exemple 4.3 Équilibrez x CH₄ + y O₂ → z CO₂ + w H₂O.

Afficher les étapes
  1. C: x = z. H: 4x = 2w. O: 2y = 2z + w.
  2. Set x = 1: z = 1, w = 2, y = 2.

Réponse : CH₄ + 2 O₂ → CO₂ + 2 H₂O.

Exercices

Faites chacun d’abord sur papier, puis vérifiez.

  1. Dans le modèle résolu, quelle production du secteur 1 faut-il au total pour 1 $ de demande finale adressée au secteur 2 ?

    Afficher la réponse

    Coefficient (1, 2) de (I − A)⁻¹ = 0.3/0.45 = 0.667.

  2. Une économie à un seul secteur consomme 0.25 de sa propre production. Quelle production totale satisfait une demande de 300 ?

    Afficher la réponse

    x = 300/(1 − 0.25) = 400.

  3. Un nœud reçoit des flux entrants de 40 et 25, et a un flux sortant f₁ plus un flux sortant de 30. Trouvez f₁.

    Afficher la réponse

    f₁ = 65 − 30 = 35.

  4. Équilibrez H₂ + O₂ → H₂O.

    Afficher la réponse

    2 H₂ + O₂ → 2 H₂O.

Idées reçues fréquentes

  • Idée reçue : La production totale est égale à la demande finale.

    En réalité : Les industries consomment aussi la production les unes des autres, donc x est plus grand que d.

  • Idée reçue : L’empreinte provient uniquement du secteur auquel on achète.

    En réalité : L’inverse de Leontief additionne tous les fournisseurs en amont.

  • Idée reçue : N’importe quelle A donne une économie sensée.

    En réalité : Les sommes des colonnes doivent être inférieures à 1 pour que (I − A)⁻¹ ait des coefficients positifs ou nuls.

Cartes mémoire

  • Matrice de consommation

    Intrants par dollar de production.

  • Modèle de Leontief

    x = Ax + d.

  • Inverse de Leontief

    (I − A)⁻¹, besoins totaux.

  • Modèle étendu

    empreinte = e·x.

  • Équilibre des flux

    Entrée = sortie en chaque nœud.

  • Système homogène

    Ax = 0.

Quiz d’autocontrôle

1. La solution de Leontief est…

Afficher la réponse

x = (I − A)⁻¹d. À partir de x = Ax + d.

2. La colonne j de A décrit…

Afficher la réponse

ce qu’achète le secteur j par dollar de production. Les intrants par unité de production.

3. Dans un modèle étendu, une empreinte est…

Afficher la réponse

e·x. L’intensité fois la production, sommée.

4. Équilibrer une équation chimique est un…

Afficher la réponse

système linéaire homogène. Le décompte des atomes donne Ax = 0.

Comment cela aide à sauver la planète

L’EPA des États-Unis utilise exactement ce modèle pour suivre la pollution à travers l’économie. Selon l’EPA, ses modèles USEEIO combinent des données sur les transactions économiques entre 389 secteurs industriels avec des données environnementales sur ces secteurs, couvrant l’utilisation des sols, de l’eau, de l’énergie et des minéraux, la pollution de l’air, les nutriments et les substances toxiques, afin d’évaluer les impacts des biens et services. L’EPA les décrit comme utiles pour le calcul de l’empreinte carbone, et ils soutiennent son travail de gestion durable des matériaux.

Une matrice de consommation à 389 secteurs a 389² = 151321 coefficients, bien trop pour un calcul à la main mais facile pour les méthodes d’élimination et LU du module 3. C’est l’inverse de Leontief qui permet à une ville, une entreprise ou un ménage de voir les émissions cachées dans ses chaînes d’approvisionnement, et pas seulement à la sortie de l’usine.

Sur Saving Planets : Solutions · Chimie générale I

Sources : U.S. EPA: US Environmentally-Extended Input-Output (USEEIO) models · Interactive Linear Algebra, 1.1 Systems of linear equations

Ressources d’apprentissage gratuites

Haut de page

Unité 5 · Semaines 6–7

Espaces vectoriels, sous-espaces, indépendance, base, rang et les quatre sous-espaces fondamentaux

Espaces vectoriels et sous-espaces, espace des colonnes et noyau, résolution complète de Ax = 0 et Ax = b, indépendance linéaire, base et dimension, rang, théorème du rang et les quatre sous-espaces fondamentaux de Strang.

Objectifs d’apprentissage

À la fin de cette unité, vous saurez :

  • Tester si un ensemble est un sous-espace.
  • Trouver des bases de l’espace des colonnes et du noyau d’une matrice.
  • Tester l’indépendance linéaire de vecteurs.
  • Utiliser le théorème du rang et les dimensions des quatre sous-espaces fondamentaux.

Notions clés

  • Un sous-espace contient le vecteur nul et est stable par addition et multiplication par un scalaire. Les droites et plans passant par l’origine sont des sous-espaces ; ceux qui ne passent pas par l’origine n’en sont pas.
  • Espace des colonnes C(A) : tous les Ax, l’ensemble des b pour lesquels Ax = b a une solution. Noyau N(A) : tous les x tels que Ax = 0.
  • Des vecteurs sont linéairement indépendants si la seule combinaison donnant zéro est celle dont tous les coefficients sont nuls. Une base est une famille génératrice indépendante ; toutes les bases d’un espace ont le même nombre de vecteurs, sa dimension.
  • Rang r = nombre de pivots = dim C(A). Théorème du rang : r + dim N(A) = n, le nombre de colonnes.
  • Quatre sous-espaces fondamentaux d’une matrice m × n de rang r : C(A) (dim r) et N(Aᵀ) (dim m − r) dans ℝᵐ ; C(Aᵀ) (dim r) et N(A) (dim n − r) dans ℝⁿ. N(A) est orthogonal à l’espace des lignes.

Équations clés

C(A)={Ax},  N(A)={x:Ax=0}C(A)=\{A\mathbf{x}\},\ \ N(A)=\{\mathbf{x}:A\mathbf{x}=0\}
Espace des colonnes et noyau
r+dim⁡N(A)=nr+\dim N(A)=n
Théorème du rang
x=xp+xn\mathbf{x}=\mathbf{x}_p+\mathbf{x}_n
Solution complète de Ax = b
dim⁡N(AT)=m−r\dim N(A^T)=m-r
Noyau à gauche

Exemples résolus pas à pas

Exemple 5.1 Trouvez le rang et une base du noyau de A = [[1, 2, 3], [2, 4, 6], [1, 0, 1]].

Afficher les étapes
  1. R₂ − 2R₁ = 0 ; R₃ − R₁ = [0, −2, −2]. Deux pivots, donc r = 2.
  2. Ax = 0 : x + z = 0 et −2y − 2z = 0, avec z libre.
  3. Prenez z = 1 : x = −1, y = −1.

Réponse : r = 2, N(A) engendré par (−1, −1, 1) ; dim N(A) = 3 − 2 = 1.

Exemple 5.2 (1, 2, 3), (4, 5, 6), (7, 8, 9) sont-ils indépendants ?

Afficher les étapes
  1. Échelonnez la matrice ayant ces colonnes : R₂ − 2R₁ et R₃ − 3R₁ donnent les lignes [0, −3, −6] et [0, −6, −12], qui sont proportionnelles.
  2. Seulement deux pivots.

Réponse : No: (7, 8, 9) = 2(4, 5, 6) − (1, 2, 3).

Exemple 5.3 Une matrice 5 × 8 est de rang 3. Donnez les dimensions de ses quatre sous-espaces fondamentaux.

Afficher les étapes
  1. C(A) et C(Aᵀ) : r = 3.
  2. N(A): n − r = 8 − 3; N(Aᵀ): m − r = 5 − 3.

Réponse : 3, 3, 5 et 2.

Exercices

Faites chacun d’abord sur papier, puis vérifiez.

  1. La droite y = 2x + 1 est-elle un sous-espace de ℝ² ?

    Afficher la réponse

    Non : elle ne contient pas (0, 0).

  2. Une matrice 4 × 4 a un noyau de dimension 2. Quel est son rang ?

    Afficher la réponse

    4 − 2 = 2.

  3. Quelle est la dimension de l’espace des matrices 2 × 2 ?

    Afficher la réponse

    4.

  4. 4 vecteurs de ℝ³ peuvent-ils être indépendants ?

    Afficher la réponse

    Non : au plus 3.

Idées reçues fréquentes

  • Idée reçue : Tout plan de ℝ³ est un sous-espace.

    En réalité : Seulement les plans passant par l’origine.

  • Idée reçue : Le rang compte les coefficients non nuls.

    En réalité : Il compte les pivots (colonnes indépendantes).

  • Idée reçue : Le noyau vit dans le même espace que b.

    En réalité : N(A) est dans ℝⁿ (côté x) ; b est dans ℝᵐ.

Cartes mémoire

  • Sous-espace

    Contient 0 ; stable par + et par multiplication scalaire.

  • Espace des colonnes

    All Ax.

  • Noyau

    Tous les x tels que Ax = 0.

  • Base

    Famille génératrice indépendante.

  • Rang

    Nombre de pivots.

  • Théorème du rang

    r + dim N(A) = n.

Quiz d’autocontrôle

1. Le théorème du rang dit que r + dim N(A) est égal au…

Afficher la réponse

n. Nombre de colonnes.

2. Ax = b a une solution exactement quand b est dans…

Afficher la réponse

C(A). b doit être une combinaison des colonnes.

3. Le noyau est orthogonal à…

Afficher la réponse

l’espace des lignes. Chaque ligne scalaire x vaut 0.

4. Une base de ℝ³ contient combien de vecteurs ?

Afficher la réponse

3. Dimension 3.

Comment cela aide à sauver la planète

La dimension compte les inconnues, et de meilleures prévisions en demandent davantage. La NOAA a indiqué que la mise à niveau de ses supercalculateurs en 2023 permettrait au Global Forecast System des États-Unis d’améliorer sa résolution horizontale de 13 km à 9 km, ce qui aide à voir et à simuler des phénomènes de plus petite échelle.

Une grille de pas h compte environ (1/h)² colonnes par unité de surface, donc passer de 13 km à 9 km multiplie le nombre de colonnes de la grille par (13/9)² = 2.09. Chacune de ces colonnes porte la température, le vent, l’humidité et la pression à de nombreux niveaux, si bien que le vecteur d’état d’une prévision vit dans un espace de dimension bien plus grande, et les matrices qui le mettent à jour grossissent d’autant.

Sur Saving Planets : Analyse III · Équations différentielles

Sources : NOAA: NOAA completes upgrade to weather and climate supercomputer system (Aug. 10, 2023) · Interactive Linear Algebra, 2.7 Basis and dimension

Ressources d’apprentissage gratuites

Haut de page

Unité 6 · Semaines 8–9

Orthogonalité, projections, moindres carrés et Gram–Schmidt

Vecteurs et sous-espaces orthogonaux, projection sur une droite et sur un sous-espace, équations normales AᵀAx̂ = Aᵀb, ajustement de droites et de courbes à des données par moindres carrés, bases orthonormées, Gram–Schmidt et A = QR.

Objectifs d’apprentissage

À la fin de cette unité, vous saurez :

  • Projeter un vecteur sur une droite et sur un espace des colonnes.
  • Ajuster une droite à des données avec les équations normales.
  • Construire une base orthonormée par Gram–Schmidt.
  • Expliquer pourquoi les moindres carrés sont une projection.

Notions clés

  • La projection de b sur la droite passant par a est p = (a·b/a·a) a ; l’erreur e = b − p est orthogonale à a.
  • Quand Ax = b n’a pas de solution (plus d’équations que d’inconnues, données bruitées), on choisit x̂ pour rendre ‖b − Ax̂‖ aussi petit que possible. L’erreur doit être orthogonale à C(A), ce qui donne les équations normales AᵀAx̂ = Aᵀb.
  • Pour une droite y = C + Dt passant par des points (tᵢ, yᵢ), les colonnes de A sont une colonne de uns et les valeurs tᵢ.
  • Des vecteurs orthonormés qᵢ vérifient qᵢ·qⱼ = 0 (i ≠ j) et ‖qᵢ‖ = 1. Si Q a des colonnes orthonormées, QᵀQ = I et les projections deviennent p = QQᵀb.
  • Gram–Schmidt retranche à chaque nouveau vecteur ses projections sur les précédents, puis normalise. On obtient A = QR avec R triangulaire supérieure.

Équations clés

p=a⋅ba⋅a a\mathbf{p}=\frac{\mathbf{a}\cdot\mathbf{b}}{\mathbf{a}\cdot\mathbf{a}}\,\mathbf{a}
Projection sur une droite
ATAx^=ATbA^TA\hat{\mathbf{x}}=A^T\mathbf{b}
Équations normales
P=A(ATA)−1ATP=A(A^TA)^{-1}A^T
Matrice de projection
A=QRA=QR
Factorisation de Gram–Schmidt

Exemples résolus pas à pas

Exemple 6.1 Projetez b = (2, 3) sur a = (1, 1).

Afficher les étapes
  1. a·b = 5, a·a = 2.
  2. p = 2.5(1, 1) ; e = b − p = (−0.5, 0.5), et e·a = 0.

Réponse : p = (2.5, 2.5).

Exemple 6.2 Ajustez y = C + Dt à (0, 1), (1, 2), (2, 4) par moindres carrés.

Afficher les étapes
  1. A a pour lignes (1, 0), (1, 1), (1, 2) ; AᵀA = [[3, 3], [3, 5]], Aᵀb = (7, 10).
  2. det = 6 ; C = (5 × 7 − 3 × 10)/6, D = (3 × 10 − 3 × 7)/6.

Réponse : C = 0.833, D = 1.5: y = 0.833 + 1.5t.

Exemple 6.3 Appliquez Gram–Schmidt à a₁ = (1, 1, 0), a₂ = (1, 0, 1).

Afficher les étapes
  1. q₁ = (1, 1, 0)/√2.
  2. a₂ − (a₂·q₁)q₁ = (1, 0, 1) − ½(1, 1, 0) = (0.5, −0.5, 1), de longueur √1.5 = 1.225.

Réponse : q₁ = (0.707, 0.707, 0), q₂ = (0.408, −0.408, 0.816).

Exercices

Faites chacun d’abord sur papier, puis vérifiez.

  1. Projetez (3, 4) sur (1, 0).

    Afficher la réponse

    (3, 0).

  2. Écrivez les équations normales pour ajuster y = Dt (sans ordonnée à l’origine) à (1, 2), (2, 3).

    Afficher la réponse

    Σt² D = Σty: 5D = 8, D = 1.6.

  3. Q = [[0, 1], [1, 0]] est-elle orthogonale ?

    Afficher la réponse

    Yes: QᵀQ = I.

  4. Si b est déjà dans C(A), quelle est l’erreur des moindres carrés ?

    Afficher la réponse

    Zéro : p = b.

Idées reçues fréquentes

  • Idée reçue : Les moindres carrés résolvent Ax = b exactement.

    En réalité : Ils trouvent le Ax̂ le plus proche quand il n’existe pas de solution exacte.

  • Idée reçue : N’importe quelle base rend les projections faciles.

    En réalité : Les bases orthonormées oui ; sinon il faut (AᵀA)⁻¹.

  • Idée reçue : Un ajustement par une droite signifie que les données sont linéaires.

    En réalité : Vérifiez les résidus ; une courbe peut mieux convenir (ajoutez une colonne t²).

Cartes mémoire

  • Projection sur une droite

    p = (a·b/a·a)a.

  • Équations normales

    AᵀAx̂ = Aᵀb.

  • Matrice de projection

    P = A(AᵀA)⁻¹Aᵀ.

  • Orthonormé

    qᵢ·qⱼ = 0, ‖qᵢ‖ = 1.

  • Gram–Schmidt

    Retrancher les projections, normaliser.

  • QR

    A = QR.

Quiz d’autocontrôle

1. Les équations normales sont…

Afficher la réponse

AᵀAx̂ = Aᵀb. Erreur orthogonale à C(A).

2. Le vecteur d’erreur des moindres carrés est orthogonal…

Afficher la réponse

l’espace des colonnes. C’est ce qui le rend le plus court.

3. Pour Q à colonnes orthonormées, QᵀQ =…

Afficher la réponse

I. Produits scalaires 0 ou 1.

4. Gram–Schmidt produit…

Afficher la réponse

A = QR. Q orthonormée, R triangulaire.

Comment cela aide à sauver la planète

La montée du niveau de la mer se mesure en ajustant des courbes à des données satellitaires bruitées. La NASA indique que les altimètres satellitaires, TOPEX/Poseidon, Jason-1, -2 et -3 et désormais Sentinel-6 Michael Freilich, suivent le niveau moyen mondial de la mer depuis 1993, qu’il a monté de 114.4 mm depuis, et que le rythme de la hausse est passé de 0.08 pouce par an (0.20 cm par an) en 1993 à 0.17 pouce par an (0.44 cm par an) aujourd’hui, soit environ 2.2 fois plus vite.

Un ajustement linéaire par moindres carrés donne un rythme moyen ; ajouter une colonne t² à A ajuste une courbe dont la pente variable mesure l’accélération décrite par la NASA. Le TP D ajuste des données d’exemple ayant l’allure de la série de la NASA avec exactement les équations normales de ce module.

Sur Saving Planets : Océans · Observation de la Terre

Sources : NASA: Global mean sea level (Vital Signs) · Interactive Linear Algebra, 6.5 The method of least squares

Ressources d’apprentissage gratuites

Haut de page

Unité 7 · Semaine 10

Déterminants : propriétés, cofacteurs, règle de Cramer et volume

Le déterminant à partir de ses propriétés, son calcul par élimination et par développement en cofacteurs, det(AB) = det A det B, le critère d’inversibilité, la règle de Cramer, et le déterminant comme facteur d’échelle d’aire ou de volume d’une application linéaire.

Objectifs d’apprentissage

À la fin de cette unité, vous saurez :

  • Calculer des déterminants 2 × 2 et 3 × 3 par cofacteurs et par élimination.
  • Utiliser les propriétés des déterminants sous les opérations sur les lignes.
  • Résoudre de petits systèmes par la règle de Cramer.
  • Interpréter |det A| comme un facteur d’échelle d’aire ou de volume.

Notions clés

  • det A ≠ 0 exactement quand A est inversible. Le déterminant d’une matrice triangulaire est le produit de ses coefficients diagonaux, donc det A = ± (produit des pivots).
  • Opérations sur les lignes : un échange change le signe, multiplier une ligne par c multiplie det par c, et ajouter à une ligne un multiple d’une autre laisse det inchangé.
  • det(AB) = det A · det B, det(Aᵀ) = det A et det(A⁻¹) = 1/det A.
  • Développement en cofacteurs selon une ligne : det A = Σ aᵢⱼCᵢⱼ avec Cᵢⱼ = (−1)ⁱ⁺ʲ fois le mineur. Règle de Cramer : xᵢ = det(Bᵢ)/det A, où Bᵢ a b en colonne i.
  • Une application linéaire x ↦ Ax multiplie les aires (2D) ou les volumes (3D) par |det A| ; un déterminant négatif inverse aussi l’orientation.

Équations clés

det⁡[abcd]=ad−bc\det\begin{bmatrix}a&b\\c&d\end{bmatrix}=ad-bc
Déterminant 2 × 2
det⁡(AB)=det⁡A det⁡B\det(AB)=\det A\,\det B
Règle du produit
xi=det⁡Bidet⁡Ax_i=\frac{\det B_i}{\det A}
Règle de Cramer
area(AS)=∣det⁡A∣ area(S)\text{area}(AS)=|\det A|\,\text{area}(S)
Facteur d’échelle des aires

Exemples résolus pas à pas

Exemple 7.1 Trouvez det [[1, 2, 0], [0, 1, 3], [2, 0, 1]] par cofacteurs selon la ligne 1.

Afficher les étapes
  1. 1 × (1 × 1 − 3 × 0) = 1.
  2. − 2 × (0 × 1 − 3 × 2) = −2 × (−6) = 12.
  3. + 0 × (…) = 0.

Réponse : det = 13.

Exemple 7.2 Résolvez 2x + y = 5, x + 3y = 10 par la règle de Cramer.

Afficher les étapes
  1. det A = 6 − 1 = 5.
  2. x = det[[5, 1], [10, 3]]/5 = (15 − 10)/5.
  3. y = det[[2, 5], [1, 10]]/5 = (20 − 5)/5.

Réponse : x = 1, y = 3.

Exemple 7.3 Trouvez l’aire du parallélogramme de côtés (3, 0) et (1, 2).

Afficher les étapes
  1. Aire = |det [[3, 1], [0, 2]]| = |6 − 0|.

Réponse : 6 unités d’aire.

Exercices

Faites chacun d’abord sur papier, puis vérifiez.

  1. det [[3, 1], [2, 4]] ?

    Afficher la réponse

    12 − 2 = 10.

  2. det [[2, 0, 1], [1, 3, 2], [1, 1, 1]] ?

    Afficher la réponse

    2(3 − 2) − 0 + 1(1 − 3) = 0, elle est donc singulière.

  3. Si det A = 4 pour une matrice A 3 × 3, que vaut det(2A) ?

    Afficher la réponse

    2³ × 4 = 32.

  4. Si det A = 5, que vaut det(A⁻¹) ?

    Afficher la réponse

    1/5 = 0.2.

Idées reçues fréquentes

  • Idée reçue : det(A + B) = det A + det B.

    En réalité : Faux en général ; seul det(AB) se factorise.

  • Idée reçue : det(2A) = 2 det A.

    En réalité : Pour n × n, det(2A) = 2ⁿ det A.

  • Idée reçue : Le développement en cofacteurs est la méthode rapide.

    En réalité : Il demande environ n! étapes ; l’élimination environ n³/3.

Cartes mémoire

  • det 2 × 2

    ad − bc.

  • Inversibilité

    det A ≠ 0.

  • Règle du produit

    det(AB) = det A det B.

  • Échange de lignes

    Change le signe.

  • Règle de Cramer

    xᵢ = det Bᵢ/det A.

  • Facteur d’aire

    |det A|.

Quiz d’autocontrôle

1. Échanger deux lignes multiplie det par…

Afficher la réponse

−1. Le signe change.

2. det A = 0 signifie que A est…

Afficher la réponse

singulière. Non inversible.

3. Le det d’une matrice triangulaire est…

Afficher la réponse

le produit des coefficients diagonaux. Développez le long de la diagonale.

4. Une application de det = 3 multiplie les aires par…

Afficher la réponse

3. |det A|.

Comment cela aide à sauver la planète

Convertir des pixels en surface au sol est un déterminant. L’USGS indique pour Landsat 9 des bandes multispectrales à 30 m et une bande panchromatique à 15 m. L’application qui envoie les indices de pixel (i, j) sur les coordonnées au sol est, localement, une matrice ; pour une grille alignée de 30 m c’est diag(30, 30) m, de déterminant 900, donc chaque pixel couvre 900 m². La grille panchromatique de 15 m a pour déterminant 225 m², si bien que quatre pixels panchromatiques tiennent dans un pixel multispectral.

Compter les pixels classés comme forêt et multiplier par ce facteur d’échelle, c’est ainsi que les relevés satellitaires mesurent la perte de forêt en hectares (1 ha = 10000 m², environ 11.1 pixels Landsat). Quand une image est tournée ou cisaillée pour coller à une carte, le déterminant de la nouvelle matrice garde l’aire exacte.

Sur Saving Planets : Déforestation · Observation de la Terre

Sources : USGS: Landsat 9 · Interactive Linear Algebra, 4.3 Determinants and volumes

Ressources d’apprentissage gratuites

Haut de page

Unité 8 · Semaine 11

Valeurs propres et vecteurs propres

Les vecteurs propres comme directions qu’une matrice ne fait qu’étirer, l’équation caractéristique det(A − λI) = 0, les vérifications par la trace et le déterminant, les valeurs propres de matrices particulières, valeurs propres complexes et rotations, la méthode de la puissance itérée, et les modes propres des structures en vibration.

Objectifs d’apprentissage

À la fin de cette unité, vous saurez :

  • Trouver les valeurs propres via det(A − λI) = 0 et les vecteurs propres dans le noyau de A − λI.
  • Vérifier les valeurs propres avec la trace et le déterminant.
  • Reconnaître des valeurs propres réelles, multiples et complexes.
  • Utiliser la puissance itérée pour trouver la valeur propre dominante.

Notions clés

  • Av = λv avec v ≠ 0 : le vecteur propre v garde sa direction et est multiplié par la valeur propre λ.
  • Les valeurs propres sont les racines de det(A − λI) = 0, un polynôme de degré n. Pour chaque λ, les vecteurs propres remplissent le noyau de A − λI.
  • Vérifications : la somme des valeurs propres est la trace (somme de la diagonale) et leur produit est det A.
  • Les matrices symétriques ont des valeurs propres réelles et des vecteurs propres orthogonaux. Les matrices de rotation ont des valeurs propres complexes. Les matrices triangulaires affichent leurs valeurs propres sur la diagonale.
  • Puissance itérée : multipliez sans cesse un vecteur de départ par A et renormalisez ; il s’aligne sur le vecteur propre du plus grand |λ|. Pour les vibrations, Kv = ω²Mv donne les pulsations propres ω et les déformées modales v.

Équations clés

Av=λvA\mathbf{v}=\lambda\mathbf{v}
Équation aux valeurs propres
det⁡(A−λI)=0\det(A-\lambda I)=0
Équation caractéristique
λ1+⋯+λn=tr⁡A,  λ1⋯λn=det⁡A\lambda_1+\cdots+\lambda_n=\operatorname{tr}A,\ \ \lambda_1\cdots\lambda_n=\det A
Vérifications par la trace et le déterminant
Kv=ω2MvK\mathbf{v}=\omega^2M\mathbf{v}
Modes propres d’une structure

Exemples résolus pas à pas

Exemple 8.1 Trouvez les valeurs propres et vecteurs propres de A = [[2, 1], [1, 2]].

Afficher les étapes
  1. det(A − λI) = (2 − λ)² − 1 = λ² − 4λ + 3 = 0.
  2. λ = 3 : A − 3I = [[−1, 1], [1, −1]] donne v = (1, 1). λ = 1 : v = (1, −1).
  3. Vérification : 3 + 1 = 4 = trace ; 3 × 1 = 3 = det.

Réponse : λ = 3 avec (1, 1) ; λ = 1 avec (1, −1).

Exemple 8.2 Trouvez les valeurs propres de [[4, 1], [2, 3]].

Afficher les étapes
  1. λ² − 7λ + 10 = 0 (trace 7, det 10).
  2. (λ − 5)(λ − 2) = 0.

Réponse : λ = 5 et λ = 2.

Exemple 8.3 Un bâtiment d’exemple à deux niveaux a des masses de plancher égales m et une rigidité d’étage k, avec K/m = (k/m)[[2, −1], [−1, 1]] et k/m = 100 s⁻². Trouvez la fréquence fondamentale.

Afficher les étapes
  1. Valeurs propres de [[2, −1], [−1, 1]] : λ = (3 ± √5)/2 = 0.382 et 2.618.
  2. ω₁² = 100 × 0.382 = 38.2 s⁻², ω₁ = 6.18 rad/s.
  3. f₁ = ω₁/(2π), T₁ = 1/f₁.

Réponse : f₁ = 0.984 Hz, période T₁ = 1.02 s (valeurs d’exemple).

Exercices

Faites chacun d’abord sur papier, puis vérifiez.

  1. Valeurs propres de [[5, 0], [0, −2]] ?

    Afficher la réponse

    5 et −2 (diagonale).

  2. Une matrice 3 × 3 a pour valeurs propres 1, 2, 3. Trouvez sa trace et son déterminant.

    Afficher la réponse

    Trace 6, det 6.

  3. Valeurs propres de la rotation de 90° [[0, −1], [1, 0]] ?

    Afficher la réponse

    λ² + 1 = 0: λ = ±i.

  4. Si Av = 4v, que vaut A²v ?

    Afficher la réponse

    16v.

Idées reçues fréquentes

  • Idée reçue : Le vecteur nul est un vecteur propre.

    En réalité : Les vecteurs propres doivent être non nuls.

  • Idée reçue : L’échelonnement conserve les valeurs propres.

    En réalité : Il les modifie ; utilisez det(A − λI).

  • Idée reçue : Toute matrice réelle a des valeurs propres réelles.

    En réalité : Les rotations en donnent des complexes ; les matrices symétriques sont le cas sûr.

Cartes mémoire

  • Vecteur propre

    Av = λv, v ≠ 0.

  • Équation caractéristique

    det(A − λI) = 0.

  • Vérification par la trace

    Σλ = trace.

  • Vérification par le déterminant

    Πλ = det A.

  • Puissance itérée

    Des multiplications répétées trouvent le vecteur propre dominant.

  • Modes propres

    Kv = ω²Mv.

Quiz d’autocontrôle

1. Les valeurs propres d’une matrice triangulaire sont…

Afficher la réponse

ses coefficients diagonaux. det(A − λI) est un produit de termes diagonaux.

2. La somme des valeurs propres est égale…

Afficher la réponse

la trace. Vérification par la trace.

3. La puissance itérée trouve…

Afficher la réponse

la valeur propre dominante. Le plus grand |λ| l’emporte.

4. Les matrices symétriques ont…

Afficher la réponse

des valeurs propres réelles et des vecteurs propres orthogonaux. Théorème spectral.

Comment cela aide à sauver la planète

Les ingénieurs parasismiques lisent les bâtiments à travers leurs valeurs propres. Une étude de l’USGS de Mehmet Çelebi (2022) a utilisé les périodes fondamentales et les pourcentages d’amortissement critique de 41 immeubles de grande hauteur dans plusieurs pays. Elle a trouvé une nette tendance linéaire entre la hauteur du bâtiment et la période fondamentale pour les bâtiments en acier et en béton armé, avec toutefois une grande dispersion due aux différences de conception structurelle, et un amortissement critique nettement inférieur à 3 %, en accord avec d’autres études récentes sur les tours.

La période fondamentale est fixée par la plus petite valeur propre du problème rigidité–masse, comme dans l’exemple du bâtiment à deux niveaux. Un faible amortissement signifie qu’une tour peut vibrer longtemps quand des ondes sismiques de longue période correspondent à cette période ; c’est pourquoi les bâtiments instrumentés et leurs modes mesurés guident des villes plus sûres et plus durables.

Sur Saving Planets : Équations différentielles · Dynamique

Sources : USGS: Preliminary study of correlation of natural periods and damping percentages of tall buildings in several countries (Çelebi, 2022) · Interactive Linear Algebra, 5.1 Eigenvalues and eigenvectors

Ressources d’apprentissage gratuites

Haut de page

Unité 9 · Semaine 12

Diagonalisation, puissances de matrices, chaînes de Markov et systèmes dynamiques linéaires

Diagonaliser A = PDP⁻¹, calculer rapidement les puissances Aᵏ, matrices symétriques et théorème spectral A = QΛQᵀ, matrices de Markov (stochastiques) et états stationnaires, et le comportement à long terme de xₖ₊₁ = Axₖ.

Objectifs d’apprentissage

À la fin de cette unité, vous saurez :

  • Diagonaliser une matrice ayant n vecteurs propres indépendants.
  • Calculer Aᵏ et prévoir le comportement à long terme de xₖ₊₁ = Axₖ.
  • Trouver l’état stationnaire d’une chaîne de Markov.
  • Utiliser le théorème spectral pour les matrices symétriques.

Notions clés

  • Si A a n vecteurs propres indépendants (les colonnes de P), alors A = PDP⁻¹ avec les valeurs propres sur la diagonale de D, et Aᵏ = PDᵏP⁻¹.
  • Écrivez x₀ = c₁v₁ + … + cₙvₙ. Alors xₖ = Aᵏx₀ = c₁λ₁ᵏv₁ + … + cₙλₙᵏvₙ ; le plus grand |λ| domine quand k augmente.
  • Une matrice de Markov a des coefficients positifs ou nuls et des colonnes de somme 1. Elle a toujours λ = 1, et les autres valeurs propres vérifient |λ| ≤ 1. L’état stationnaire est le vecteur propre associé à λ = 1, normalisé pour que ses coefficients aient pour somme 1.
  • Le deuxième plus grand |λ| fixe la vitesse à laquelle une chaîne de Markov oublie son départ : l’écart décroît comme |λ₂|ᵏ.
  • Théorème spectral : une matrice symétrique réelle peut s’écrire A = QΛQᵀ avec des vecteurs propres orthonormés dans Q.

Équations clés

A=PDP−1,  Ak=PDkP−1A=PDP^{-1},\ \ A^k=PD^kP^{-1}
Diagonalisation et puissances
xk=c1λ1kv1+⋯+cnλnkvn\mathbf{x}_k=c_1\lambda_1^k\mathbf{v}_1+\cdots+c_n\lambda_n^k\mathbf{v}_n
Solution de x_{k+1} = Ax_k
Pq=q,  ∑qi=1P\mathbf{q}=\mathbf{q},\ \ \textstyle\sum q_i=1
État stationnaire de Markov
A=QΛQTA=Q\Lambda Q^T
Théorème spectral

Exemples résolus pas à pas

Exemple 9.1 Utilisez la diagonalisation pour trouver A¹⁰ pour A = [[2, 1], [1, 2]].

Afficher les étapes
  1. D’après le module 8 : λ = 3, 1 avec v = (1, 1), (1, −1). P = [[1, 1], [1, −1]], P⁻¹ = ½[[1, 1], [1, −1]].
  2. A¹⁰ = P diag(3¹⁰, 1) P⁻¹ ; chaque coefficient vaut (3¹⁰ ± 1)/2.

Réponse : A¹⁰ = [[29525, 29524], [29524, 29525]].

Exemple 9.2 Une matrice de Markov est P = [[0.8, 0.3], [0.2, 0.7]] (valeurs d’exemple). Trouvez l’état stationnaire.

Afficher les étapes
  1. (P − I)q = 0 : −0.2q₁ + 0.3q₂ = 0, donc q₁ = 1.5q₂.
  2. q₁ + q₂ = 1 donne q₂ = 0.4.

Réponse : q = (0.6, 0.4).

Exemple 9.3 En partant de x₀ = (1, 0), trouvez x₂ pour la même P, ainsi que xₖ en général.

Afficher les étapes
  1. x₁ = (0.8, 0.2); x₂ = (0.64 + 0.06, 0.16 + 0.14).
  2. Les valeurs propres sont 1 et 0.5 (trace 1.5, det 0.5), donc xₖ = (0.6, 0.4) + 0.5ᵏ(0.4, −0.4).

Réponse : x₂ = (0.70, 0.30); xₖ → (0.6, 0.4).

Exercices

Faites chacun d’abord sur papier, puis vérifiez.

  1. [[0.5, 0.2], [0.5, 0.8]] est-elle une matrice de Markov ?

    Afficher la réponse

    Oui : coefficients positifs ou nuls, colonnes de somme 1.

  2. Trouvez son état stationnaire.

    Afficher la réponse

    0.5q₁ = 0.2q₂, donc q = (2/7, 5/7) = (0.286, 0.714).

  3. Si A = PDP⁻¹ avec D = diag(0.9, 0.5), que devient Aᵏx pour k grand ?

    Afficher la réponse

    Il tend vers zéro, dominé par 0.9ᵏ.

  4. [[1, 1], [0, 1]] est-elle diagonalisable ?

    Afficher la réponse

    Non : λ = 1 double mais un seul vecteur propre indépendant.

Idées reçues fréquentes

  • Idée reçue : Toute matrice est diagonalisable.

    En réalité : Il faut n vecteurs propres indépendants ; [[1, 1], [0, 1]] n’en a qu’un.

  • Idée reçue : L’état stationnaire d’une chaîne de Markov dépend du départ.

    En réalité : Pour une chaîne régulière, non ; tout départ converge vers le même q.

  • Idée reçue : Aᵏ revient à élever chaque coefficient à la puissance k.

    En réalité : Seules les matrices diagonales fonctionnent ainsi ; c’est pourquoi on diagonalise.

Cartes mémoire

  • Diagonalisation

    A = PDP⁻¹.

  • Puissance de matrice

    Aᵏ = PDᵏP⁻¹.

  • Matrice de Markov

    Coefficients positifs ou nuls, colonnes de somme 1.

  • État stationnaire

    Pq = q, coefficients de somme 1.

  • Théorème spectral

    A = QΛQᵀ.

  • Matrice non diagonalisable

    Trop peu de vecteurs propres pour diagonaliser.

Quiz d’autocontrôle

1. Toute matrice de Markov a la valeur propre…

Afficher la réponse

1. Colonnes de somme 1.

2. Aᵏ pour A = PDP⁻¹ vaut…

Afficher la réponse

PDᵏP⁻¹. Les P⁻¹P intérieurs se simplifient.

3. La vitesse de convergence d’une chaîne de Markov est fixée par…

Afficher la réponse

|λ₂|. L’écart décroît comme |λ₂|ᵏ.

4. Les matrices symétriques se diagonalisent sous la forme…

Afficher la réponse

QΛQᵀ. Théorème spectral.

Comment cela aide à sauver la planète

Le devenir des matériaux est une chaîne de Markov. L’EPA indique que les États-Unis ont produit 292.4 millions de tonnes de déchets ménagers et assimilés en 2018, soit 4.9 livres par personne et par jour. Environ 69 millions de tonnes ont été recyclées et 25 millions de tonnes compostées, soit un taux combiné de 32.1 % ; environ 35 millions de tonnes (11.8 %) ont été incinérées avec valorisation énergétique ; et plus de 146 millions de tonnes (50 %) ont été mises en décharge.

Si les parts de chaque année restaient fixes, un matériau quitterait l’usage vers ces destinations avec ces probabilités, et seule la part recyclée reviendrait dans l’économie. Le TP E construit une matrice de transition d’exemple avec une part de valorisation de 32.1 % et montre, avec les puissances de matrices et le vecteur propre de λ = 1, comment augmenter le recyclage change la durée pendant laquelle le matériau reste en usage et la quantité qui finit en décharge.

Sur Saving Planets : Solutions

Sources : U.S. EPA: National overview, facts and figures on materials, wastes and recycling · Interactive Linear Algebra, 5.6 Stochastic matrices

Ressources d’apprentissage gratuites

Haut de page

Unité 10 · Semaines 13–14

La décomposition en valeurs singulières, l’analyse en composantes principales et les matrices définies positives

Matrices définies positives et AᵀA, la SVD A = UΣVᵀ, valeurs singulières et géométrie d’une matrice comme rotation–étirement–rotation, meilleure approximation de rang faible et compression d’images, conditionnement, et analyse en composantes principales de données.

Objectifs d’apprentissage

À la fin de cette unité, vous saurez :

  • Tester si une matrice symétrique est définie positive.
  • Calculer les valeurs singulières à partir des valeurs propres de AᵀA.
  • Construire une approximation de rang k et mesurer ce qu’elle conserve.
  • Réaliser une petite analyse en composantes principales.

Notions clés

  • Une matrice symétrique est définie positive quand xᵀAx > 0 pour tout x ≠ 0 ; de façon équivalente, toutes ses valeurs propres (ou tous ses pivots) sont positives. AᵀA est toujours symétrique et semi-définie positive.
  • Toute matrice m × n a une SVD A = UΣVᵀ : V et U ont des colonnes orthonormées, et Σ contient les valeurs singulières σ₁ ≥ σ₂ ≥ … ≥ 0, racines carrées des valeurs propres de AᵀA.
  • A = σ₁u₁v₁ᵀ + σ₂u₂v₂ᵀ + …. Garder les k premiers termes donne la meilleure approximation de rang k (Eckart–Young) ; l’erreur en norme 2 est σₖ₊₁.
  • Le conditionnement σ₁/σₙ mesure à quel point des erreurs sur b peuvent être amplifiées en résolvant Ax = b.
  • ACP : centrez les données (retranchez les moyennes des colonnes) ; les vecteurs singuliers à droite (vecteurs propres de la matrice de covariance) sont alors les directions principales, et σᵢ² mesure la variance expliquée par chacune.

Équations clés

A=UΣVTA=U\Sigma V^T
Décomposition en valeurs singulières
σi=λi(ATA)\sigma_i=\sqrt{\lambda_i(A^TA)}
Valeurs singulières
Ak=∑i=1kσiuiviTA_k=\sum_{i=1}^{k}\sigma_i\mathbf{u}_i\mathbf{v}_i^T
Meilleure approximation de rang k
κ=σ1/σn\kappa=\sigma_1/\sigma_n
Conditionnement

Exemples résolus pas à pas

Exemple 10.1 Trouvez les valeurs singulières de A = [[3, 0], [4, 5]].

Afficher les étapes
  1. AᵀA = [[25, 20], [20, 25]], valeurs propres 45 et 5.
  2. σ₁ = √45, σ₂ = √5. Vérification : σ₁σ₂ = 15 = |det A|.

Réponse : σ₁ = 6.708, σ₂ = 2.236.

Exemple 10.2 Quelle part de l’« énergie » Σσ² de AᵀA l’approximation de rang 1 de cette A conserve-t-elle ?

Afficher les étapes
  1. σ₁²/(σ₁² + σ₂²) = 45/50.

Réponse : 90%.

Exemple 10.3 Une image en niveaux de gris de 1000 × 1000 est stockée sous forme de SVD de rang 50. Combien de nombres sont stockés, par rapport à l’image complète ?

Afficher les étapes
  1. Chaque terme demande un u (1000), un v (1000) et un σ : 2001 nombres.
  2. 50 × 2001 = 100050 contre 1000000.

Réponse : 100050 nombres, environ 10.0 % de l’original.

Exercices

Faites chacun d’abord sur papier, puis vérifiez.

  1. [[2, 1], [1, 2]] est-elle définie positive ?

    Afficher la réponse

    Oui : les valeurs propres 3 et 1 sont positives.

  2. Trouvez le conditionnement de diag(10, 0.1).

    Afficher la réponse

    10/0.1 = 100.

  3. Des données 2D centrées ont pour covariance [[2, 1], [1, 2]]. Quelle part de la variance la première composante principale explique-t-elle ?

    Afficher la réponse

    Valeurs propres 3 et 1 : 3/4 = 75 %, selon (1, 1).

  4. Quelles sont les valeurs singulières d’une matrice orthogonale ?

    Afficher la réponse

    Toutes égales à 1.

Idées reçues fréquentes

  • Idée reçue : Les valeurs singulières sont les valeurs propres de A.

    En réalité : Ce sont les racines carrées des valeurs propres de AᵀA, et elles ne sont jamais négatives.

  • Idée reçue : Seules les matrices carrées ont une SVD.

    En réalité : Toutes les matrices en ont une.

  • Idée reçue : L’ACP fonctionne sur les données brutes.

    En réalité : Centrez d’abord les données, et mettez à l’échelle les variables mesurées dans des unités différentes.

Cartes mémoire

  • SVD

    A = UΣVᵀ.

  • Valeur singulière

    σᵢ = √λᵢ(AᵀA).

  • Définie positive

    xᵀAx > 0, tous les λ > 0.

  • Eckart–Young

    La SVD tronquée est la meilleure approximation de rang k.

  • Conditionnement

    σ₁/σₙ.

  • PCA

    SVD des données centrées.

Quiz d’autocontrôle

1. Les valeurs singulières sont…

Afficher la réponse

positives ou nulles. Racines carrées des valeurs propres de AᵀA, ≥ 0.

2. La meilleure approximation de rang k conserve…

Afficher la réponse

les k plus grands termes singuliers. Eckart–Young.

3. Un grand conditionnement signifie…

Afficher la réponse

que les erreurs peuvent être amplifiées. Presque singulière.

4. En ACP, la première direction principale est…

Afficher la réponse

le premier vecteur singulier à droite des données centrées. La direction de plus grande variance.

Comment cela aide à sauver la planète

Le spectromètre imageur EMIT du JPL de la NASA, à bord de la Station spatiale internationale, cartographie les minéraux des régions sources de poussière de la Terre. Le JPL indique qu’il recueille plus de 100000 spectres par seconde, soit 3.6 × 10⁸ spectres par heure, et que sa mission a été prolongée au moins jusqu’en 2026. La poussière compte pour le climat : le JPL note que les minéraux sombres absorbent la lumière du soleil et réchauffent l’air, tandis que les minéraux clairs la réfléchissent et refroidissent.

Chaque spectre est un long vecteur de réflectances à de nombreuses longueurs d’onde, et les longueurs d’onde voisines sont fortement corrélées. L’ACP et la SVD compressent ces données en quelques composantes, séparent le signal du bruit et isolent les signatures de minéraux individuels : les mêmes mathématiques que le TP de compression d’images.

Sur Saving Planets : Observation de la Terre · Espace

Sources : NASA JPL: EMIT, Earth Surface Mineral Dust Source Investigation · MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010)

Ressources d’apprentissage gratuites

Haut de page

Le laboratoire (Computation lab)

Le labo est la partie pratique du cours. Ici, il est gratuit et sert uniquement à s’entraîner, sans crédits.

Chaque TP se fait dans un tableur gratuit ou sur papier ; les versions Python sont facultatives. Faites des pauses régulières loin de l’écran, et manipulez épingles et ciseaux avec précaution.

Référence de sécurité : OSHA: Computer workstations eTool

TP A : vecteurs de pixels et indice de végétation

Semaine 2

Objectif : Traiter les pixels satellitaires comme des vecteurs : calculer des longueurs, des angles entre pixels et le NDVI, et voir que végétation, sol et eau pointent dans des directions différentes.

cos⁡θ=p⋅q∥p∥∥q∥,NDVI=NIR−redNIR+red\cos\theta=\frac{\mathbf{p}\cdot\mathbf{q}}{\|\mathbf{p}\|\|\mathbf{q}\|},\quad \text{NDVI}=\frac{\text{NIR}-\text{red}}{\text{NIR}+\text{red}}
Angle entre pixels et NDVI

Mode opératoire

  1. Dans un tableur, saisissez trois pixels d’exemple comme lignes de réflectances (vert, rouge, NIR) : végétation (0.10, 0.08, 0.40), sol (0.20, 0.25, 0.30), eau (0.06, 0.04, 0.02).
  2. Calculez la longueur de chaque pixel avec SQRT(SUMSQ(…)) et le produit scalaire de chaque paire avec SUMPRODUCT.
  3. Trouvez l’angle entre chaque paire avec ACOS et DEGREES.
  4. Calculez le NDVI de chaque pixel. Facultatif : téléchargez une scène Landsat gratuite et recommencez avec de vrais pixels.

Tableau de données

PixelVertRedNIRNDVI
Végétation0.100.080.400.667
Sol0.200.250.300.091
Eau0.060.040.02−0.333

Réflectances d’exemple choisies pour ressembler à une végétation, un sol nu et une eau typiques.

Exemple de calcul

  1. NDVI végétation = (0.40 − 0.08)/(0.40 + 0.08) = 0.667.
  2. NDVI sol = 0.05/0.55 = 0.091 ; NDVI eau = −0.02/0.06 = −0.333.
  3. Angle végétation–sol : produit scalaire = 0.160, longueurs 0.420 et 0.4387, cos θ = 0.8683, θ = 29.7°.

Sécurité

  • Faites une courte pause d’écran toutes les demi-heures.

Alternative à la maison ou virtuelle

  • Tout le TP se fait dans un tableur gratuit.

PhET : Vector Addition · USGS: Landsat 9

Lien avec la planète : L’USGS indique les bandes rouge (bande 4) et proche infrarouge (bande 5) de Landsat 9, et ses données sont gratuites, donc chacun peut cartographier la végétation ainsi. USGS: Landsat 9

TP B : trilatération, des ficelles à un système linéaire

Semaine 3

Objectif : Trouver une position à partir de distances, d’abord avec des ficelles sur une carte comme dans l’activité de GPS.gov, puis en transformant les équations de cercles en un système linéaire résolu par élimination.

2(a2−a1)x+2(b2−b1)y=d12−d22+a22−a12+b22−b122(a_2-a_1)x+2(b_2-b_1)y=d_1^2-d_2^2+a_2^2-a_1^2+b_2^2-b_1^2
Différence de deux équations de cercles

Mode opératoire

  1. Imprimez la carte de trilatération de GPS.gov (sans mise à l’échelle) ou marquez trois « satellites » sur du papier millimétré en (0, 0), (10, 0) et (0, 10) cm.
  2. Convertissez les temps de signal en distances avec d = ct (c = 299792458 m/s) ; pour la version papier, utilisez les distances d’exemple ci-dessous.
  3. Tracez des cercles avec une ficelle et un crayon et marquez leur point d’intersection.
  4. Soustrayez les équations des cercles deux à deux pour obtenir deux équations linéaires, résolvez-les par élimination dans un tableur, et comparez avec votre dessin.

Tableau de données

BalisePosition (cm)Distance (cm)d²
1(0, 0)5.00025
2(10, 0)8.06265
3(0, 10)6.70845

Distances d’exemple pour un récepteur en (3, 4) cm. L’activité de GPS.gov utilise quatre satellites sur une carte au 1:21283839.

Exemple de calcul

  1. Eq. 1 − Eq. 2: 20x − 100 = −40, x = 3. Eq. 1 − Eq. 3: 20y − 100 = −20, y = 4.
  2. Le temps A de GPS.gov = 0.00505783 s donne d = 1516299 m ; au 1:21283839 cela fait 7.12 cm sur la carte.

Sécurité

  • Manipulez punaises et ciseaux avec précaution.

Alternative à la maison ou virtuelle

  • La version ficelle et papier est le TP à faire chez soi.

GPS.gov: Trilateration activity · Interactive Linear Algebra, 1.2 Row reduction

Lien avec la planète : GPS.gov a conçu cette activité avec la NOAA ; les vrais récepteurs résolvent le même type de système, avec une inconnue d’horloge, plusieurs fois par seconde. GPS.gov: Trilateration (NOAA classroom activity) · GPS.gov: Space segment

TP C : les matrices comme transformations du plan

Semaine 4

Objectif : Voir les matrices agir sur des formes : étirer, cisailler, tourner et réfléchir un carré sur papier millimétré, et vérifier que les aires sont multipliées par |det A|.

area(AS)=∣det⁡A∣ area(S)\text{area}(AS)=|\det A|\,\text{area}(S)
Facteur d’échelle des aires

Mode opératoire

  1. Tracez le carré unité de sommets (0, 0), (1, 0), (1, 1), (0, 1).
  2. Multipliez chaque sommet par A = [[2, 1], [0, 1]] (un cisaillement et un étirement), par R = [[0, −1], [1, 0]] (une rotation de 90°) et par F = [[1, 0], [0, −1]] (une réflexion).
  3. Tracez chaque image et trouvez son aire en comptant les carreaux ou avec la formule du lacet.
  4. Comparez chaque aire avec |det| et notez quelles applications inversent l’orientation (det < 0). Vérifiez ensuite det(AR) = det A det R.

Tableau de données

MatricedetAire de l’imageOrientation
A = [[2, 1], [0, 1]]22conservée
R (rotation de 90°)11conservée
F (réflexion)−11inversée
AR22conservée

Valeurs exactes pour le carré unité.

Exemple de calcul

  1. det A = 2 × 1 − 1 × 0 = 2 ; l’image est un parallélogramme de base 2 et de hauteur 1.
  2. det(AR) = det A × det R = 2 × 1 = 2.

Sécurité

  • Aucun au-delà d’un travail de bureau normal.

Alternative à la maison ou virtuelle

  • Du papier millimétré et un crayon suffisent.

Interactive Linear Algebra, 3.1 Matrix transformations · Interactive Linear Algebra, 4.3 Determinants and volumes

Lien avec la planète : Les projections cartographiques et le recalage d’images utilisent ces mêmes matrices ; le déterminant garantit des mesures d’aire honnêtes pour les forêts et les inondations. USGS: Landsat 9

TP D : ajustement par moindres carrés de données du niveau de la mer

Semaine 9

Objectif : Ajuster une droite, puis une parabole, à des données du niveau de la mer avec les équations normales, et comparer avec l’ajustement intégré du tableur.

ATAx^=ATb,A=[ 1  t ]A^TA\hat{\mathbf{x}}=A^T\mathbf{b},\quad A=[\,\mathbf{1}\ \ \mathbf{t}\,]
Équations normales pour y = C + Dt

Mode opératoire

  1. Saisissez les données d’exemple ci-dessous (ou téléchargez la série de la NASA sur sa page Vital Signs) : t en années depuis 1993 et h en mm.
  2. Construisez A avec une colonne de uns et une colonne de t. Calculez AᵀA et Aᵀb avec MMULT et TRANSPOSE, puis x̂ = MINVERSE(AᵀA)·Aᵀb.
  3. Vérifiez avec LINEST ou une courbe de tendance de graphique. Tracez les résidus b − Ax̂.
  4. Ajoutez une colonne t² et refaites l’ajustement. Un coefficient de t² positif signifie que la hausse s’accélère.

Tableau de données

t (années depuis 1993)08162432
h (mm)0224878114

Valeurs d’exemple ayant l’allure de la série de la NASA (qui indique 114.4 mm de hausse depuis 1993) ; utilisez les vraies données pour une vraie réponse.

Exemple de calcul

  1. t̄ = 16, h̄ = 52.4; Σ(t − t̄)(h − h̄) = 2272, Σ(t − t̄)² = 640.
  2. D = 2272/640 = 3.55 mm/yr; C = 52.4 − 3.55 × 16 = −4.4 mm.
  3. Les résidus sont positifs aux deux extrémités et négatifs au milieu : une courbe convient mieux, comme le suggère le rythme croissant indiqué par la NASA.

Sécurité

  • Faites une courte pause d’écran toutes les demi-heures.

Alternative à la maison ou virtuelle

  • Tout le TP se fait dans un tableur gratuit.

NASA: Global mean sea level · Interactive Linear Algebra, 6.5 Least squares

Lien avec la planète : La NASA indique que le rythme de hausse est passé de 0.20 cm par an en 1993 à 0.44 cm par an aujourd’hui, une accélération détectable avec une colonne t². NASA: Global mean sea level (Vital Signs)

TP E : un modèle de Markov des matériaux et du recyclage

Semaine 12

Objectif : Construire une matrice de transition pour la matière en usage et perdue en décharge ou par incinération, utiliser les puissances de matrices pour la suivre dans le temps, et voir comment un taux de recyclage plus élevé garde la matière en usage plus longtemps.

xk+1=Pxk,P=[1−f(1−r)0f(1−r)1]\mathbf{x}_{k+1}=P\mathbf{x}_k,\quad P=\begin{bmatrix}1-f(1-r)&0\\f(1-r)&1\end{bmatrix}
Chaîne à deux états : en usage, perdue

Mode opératoire

  1. Choisissez un taux de mise au rebut d’exemple f = 0.2 par an et un taux de valorisation r = 0.321 (la part recyclage plus compostage de l’EPA en 2018).
  2. Construisez P dans un tableur et partez de x₀ = (1, 0) : toute la matière en usage.
  3. Utilisez MMULT pour avancer de 10 ans, et comparez avec la formule (1 − f(1 − r))ᵏ pour la part encore en usage.
  4. Recommencez avec r = 0.5 et comparez la durée moyenne d’usage en années, 1/(f(1 − r)).

Tableau de données

Taux de valorisation rProbabilité de rester en usageEn usage après 10 ansDurée moyenne d’usage (années)
0.3210.86420.2327.36
0.500.90000.34910.0

Modèle d’exemple : le taux de mise au rebut f = 0.2 est une valeur d’exemple ; r = 0.321 est la part EPA 2018.

Exemple de calcul

  1. Probabilité de rester = 1 − 0.2 × 0.679 = 0.8642 ; 0.8642¹⁰ = 0.232.
  2. Durée moyenne d’usage = 1/0.1358 = 7.36 ans ; avec r = 0.5 elle vaut 1/0.1 = 10.0 ans.
  3. λ = 1 correspond à l’état « perdue », qui absorbe tout à long terme ; λ = 0.8642 fixe la vitesse.

Sécurité

  • Faites une courte pause d’écran toutes les demi-heures.

Alternative à la maison ou virtuelle

  • Tout le TP se fait dans un tableur gratuit.

Interactive Linear Algebra, 5.6 Stochastic matrices · EPA: National overview of materials and waste

Lien avec la planète : L’EPA indique un taux de recyclage et de compostage de 32.1 % en 2018 et 50 % mis en décharge ; des modèles comme celui-ci montrent ce que vaut une hausse de ce taux. U.S. EPA: National overview, facts and figures on materials, wastes and recycling

TP F : puissance itérée et compression d’images par SVD

Semaine 14

Objectif : Trouver une valeur propre dominante par multiplications répétées, puis compresser une petite image en niveaux de gris avec une SVD tronquée et voir avec combien peu de nombres on garde l’image.

xk+1=Axk∥Axk∥,Ak=∑i=1kσiuiviT\mathbf{x}_{k+1}=\frac{A\mathbf{x}_k}{\|A\mathbf{x}_k\|},\quad A_k=\sum_{i=1}^k\sigma_i\mathbf{u}_i\mathbf{v}_i^T
Puissance itérée et SVD tronquée

Mode opératoire

  1. Dans un tableur, partez de x₀ = (1, 0) et multipliez dix fois par A = [[2, 1], [1, 2]], en renormalisant à chaque fois ; observez le rapport des longueurs tendre vers 3.
  2. Créez une grille 8 × 8 de niveaux de gris (0–255) pour une image simple, ou chargez une photo dans Python (gratuit) avec NumPy.
  3. Calculez la SVD (numpy.linalg.svd) et reconstruisez l’image avec k = 1, 2 et 4 termes.
  4. Notez le stockage k(8 + 8 + 1) et la part de Σσ² conservée, et jugez quand l’image paraît correcte.

Tableau de données

kNombres stockés (image 8 × 8)Part de 64
11726.6%
23453.1%
468106.3%

Nombres exacts de valeurs stockées. Pour une image 8 × 8, k = 4 stocke déjà plus que les 64 nombres bruts ; les gains apparaissent pour les grandes images (module 10 : le rang 50 d’une image 1000 × 1000 conserve 10.0 %).

Exemple de calcul

  1. Puissance itérée : (1, 0) → (2, 1) → (5, 4) → (14, 13) ; ‖(14, 13)‖/‖(5, 4)‖ = 19.10/6.403 = 2.98, proche de λ = 3.
  2. Stockage pour k = 1 : 8 + 8 + 1 = 17 nombres, 17/64 = 26.6 %.

Sécurité

  • Faites une courte pause d’écran toutes les demi-heures.

Alternative à la maison ou virtuelle

  • Les versions tableur et Python gratuit sont toutes deux des TP à faire chez soi.

MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010) · Interactive Linear Algebra, 5.1 Eigenvectors

Lien avec la planète : EMIT du JPL de la NASA renvoie plus de 100000 spectres par seconde ; compresser et séparer ces données utilise exactement ces outils. NASA JPL: EMIT, Earth Surface Mineral Dust Source Investigation

Haut de page

Séries d’entraînement aux examens

Chronométrez-vous : environ 10 minutes par problème. Comparez ensuite avec les solutions complètes.

Série d’exercices de mi-parcours (unités 1–5)

  1. Trouvez l’angle entre (1, 0, 1) et (0, 1, 1).

    Afficher la solution complète
    1. Produit scalaire = 1 ; longueurs √2 chacune ; cos θ = 1/2.

    Réponse : 60°.

  2. Résolvez x + 2y = 4, 3x − y = 5.

    Afficher la solution complète
    1. R₂ − 3R₁: −7y = −7.

    Réponse : x = 2, y = 1.

  3. Inversez [[3, 2], [1, 1]].

    Afficher la solution complète
    1. det = 1.

    Réponse : [[1, −2], [−1, 3]].

  4. Une économie à un seul secteur consomme 0.4 de sa propre production. Quelle production satisfait une demande de 120 ?

    Afficher la solution complète
    1. x = 120/0.6.

    Réponse : 200.

  5. Une matrice 6 × 9 est de rang 4. Trouvez dim N(A) et dim N(Aᵀ).

    Afficher la solution complète
    1. n − r et m − r.

    Réponse : 5 et 2.

Série d’entraînement finale (toutes les unités)

  1. Ajustez y = C + Dt à (0, 0), (1, 2), (2, 2) par moindres carrés.

    Afficher la solution complète
    1. AᵀA = [[3, 3], [3, 5]], Aᵀb = (4, 6) ; det 6.

    Réponse : C = 0.333, D = 1.

  2. Trouvez det [[2, 1, 0], [1, 2, 1], [0, 1, 2]].

    Afficher la solution complète
    1. 2(4 − 1) − 1(2 − 0) + 0.

    Réponse : 4.

  3. Trouvez les valeurs propres de [[6, 2], [2, 3]].

    Afficher la solution complète
    1. Trace 9, det 14 : λ² − 9λ + 14 = 0.

    Réponse : 7 et 2.

  4. Trouvez l’état stationnaire de la matrice de Markov [[0.9, 0.2], [0.1, 0.8]].

    Afficher la solution complète
    1. 0.1q₁ = 0.2q₂.

    Réponse : (0.667, 0.333).

  5. Trouvez les valeurs singulières de [[2, 0], [0, −3]].

    Afficher la solution complète
    1. AᵀA = diag(4, 9).

    Réponse : 3 et 2.

  6. Combien de nombres une SVD de rang 20 d’une image 500 × 800 stocke-t-elle, et quelle part de l’image complète cela représente-t-il ?

    Afficher la solution complète
    1. 20(500 + 800 + 1).

    Réponse : 26020, environ 6.5 %.

Haut de page

Fiche de formules

∥v∥=v⋅v, cos⁡θ=v⋅w∥v∥∥w∥\|v\|=\sqrt{v\cdot v},\ \cos\theta=\frac{v\cdot w}{\|v\|\|w\|}
Ax=b,  A=LUA\mathbf{x}=\mathbf{b},\ \ A=LU
[abcd]−1=1ad−bc[d−b−ca]\begin{bmatrix}a&b\\c&d\end{bmatrix}^{-1}=\frac{1}{ad-bc}\begin{bmatrix}d&-b\\-c&a\end{bmatrix}
(AB)−1=B−1A−1, (AB)T=BTAT(AB)^{-1}=B^{-1}A^{-1},\ (AB)^T=B^TA^T
x=(I−A)−1d\mathbf{x}=(I-A)^{-1}\mathbf{d}
r+dim⁡N(A)=nr+\dim N(A)=n
p=a⋅ba⋅aa,  ATAx^=ATb\mathbf{p}=\frac{a\cdot b}{a\cdot a}a,\ \ A^TA\hat{x}=A^Tb
A=QR, QTQ=IA=QR,\ Q^TQ=I
det⁡(AB)=det⁡Adet⁡B, xi=det⁡Bi/det⁡A\det(AB)=\det A\det B,\ x_i=\det B_i/\det A
det⁡(A−λI)=0, ∑λ=tr⁡A, ∏λ=det⁡A\det(A-\lambda I)=0,\ \sum\lambda=\operatorname{tr}A,\ \prod\lambda=\det A
A=PDP−1, Ak=PDkP−1, A=QΛQTA=PDP^{-1},\ A^k=PD^kP^{-1},\ A=Q\Lambda Q^T
A=UΣVT, σi=λi(ATA), κ=σ1/σnA=U\Sigma V^T,\ \sigma_i=\sqrt{\lambda_i(A^TA)},\ \kappa=\sigma_1/\sigma_n

Constantes

ConstanteValeur
Vitesse de la lumière (utilisée pour les distances GPS)299,792,458 m/s
1 mile1.609344 km
1 pied0.3048 m
1 pouce25.4 mm
Flops pour résoudre un système n × n par éliminationenviron n³/3 multiplications–additions (Strang)
Matrice identité IDes uns sur la diagonale, des zéros ailleurs : AI = IA = A

Source : NIST Guide to the SI, Appendix B.9: conversion factors

Haut de page

Textes et cours gratuits

Utilisez ces manuels, cours et outils gratuits en parallèle du cours.

Haut de page

Glossaire

Base
Une famille linéairement indépendante qui engendre un espace.
Espace des colonnes
Toutes les combinaisons Ax des colonnes de A.
Conditionnement
σ₁/σₙ, de combien les erreurs peuvent croître en résolvant Ax = b.
Déterminant
Facteur d’échelle signé des volumes de A ; nul exactement quand A est singulière.
Diagonalisation
A = PDP⁻¹ avec les vecteurs propres dans P.
Dimension
Nombre de vecteurs de n’importe quelle base.
Valeur propre
λ tel que Av = λv pour un v non nul.
Factorisation LU
A = LU issue de l’élimination.
Moindres carrés
Le x̂ qui minimise ‖b − Ax‖.
Inverse de Leontief
(I − A)⁻¹, production totale nécessaire par unité de demande finale.
Indépendance linéaire
Seule la combinaison nulle donne le vecteur nul.
Matrice de Markov
Coefficients positifs ou nuls, colonnes de somme 1.
Noyau
Tous les x tels que Ax = 0.
Orthonormé
Vecteurs unitaires deux à deux orthogonaux.
Pivot
Premier coefficient non nul d’une ligne d’une forme échelonnée.
Définie positive
Symétrique avec xᵀAx > 0 pour x ≠ 0.
Projection
Point le plus proche de b dans un sous-espace.
Rang
Nombre de pivots ; dimension de l’espace des colonnes.
Valeur singulière
σᵢ = √λᵢ(AᵀA).
Espace engendré
Toutes les combinaisons linéaires d’un ensemble de vecteurs.
Sous-espace
Un ensemble contenant 0, stable par addition et multiplication scalaire.

Haut de page

Grille d’autoévaluation

Aucune note ni aucun crédit ne sont attribués. Utilisez cette grille pour juger honnêtement où vous en êtes dans chaque unité.

NiveauÀ quoi cela ressemble
4 · MaîtriseChoisit le bon outil (élimination, projection, valeurs propres ou SVD), calcule les petits cas à la main, vérifie les réponses (Ax = b, trace et déterminant), et explique le sens par des images et en contexte.
3 · CompétentRésout les exercices d’entraînement avec quelques erreurs de calcul occasionnelles.
2 · En progressionSuit les exemples corrigés mais confond lignes et colonnes ou les quatre sous-espaces.
1 · DébutantManipule des symboles sans vérifier dimensions ni réponses ; doit revoir le module.

Pondération suggérée si vous voulez une note

PartiePoids
Exercices25%
Travaux pratiques20%
Série d’entraînement de mi-parcours25%
Série d’entraînement finale30%

Visez le niveau 3 ou mieux dans chaque unité avant de passer au cours suivant.

Haut de page

Aider / Faire un don