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.
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.
Programme suggéré sur 15 semaines
Environ 6 à 9 heures par semaine : lecture, exercices et, les semaines de TP, le TP.
| Semaine | Unité | Thèmes | Lecture |
|---|---|---|---|
| 1 | Vecteurs : combinaisons linéaires, normes et produits scalaires | Vecteurs, combinaisons linéaires, espace engendré, norme | MIT 18.06 L1; ILA 2.1–2.2 |
| 2 | Vecteurs : combinaisons linéaires, normes et produits scalaires | Produits scalaires, angles, orthogonalité ; TP A ; début de l’élimination | ILA 6.1 |
| 3 | Systèmes d’équations linéaires et élimination de Gauss | Élimination, formes échelonnées, ensembles de solutions ; TP B | MIT 18.06 L2; ILA 1.1–1.3 |
| 4 | Algèbre matricielle : multiplication, inverses et factorisation LU | Multiplication matricielle, inverses, A = LU ; TP C | MIT 18.06 L3–4; ILA 3.4–3.5 |
| 5 | Modèles linéaires : économie entrées–sorties (Leontief), réseaux et équilibrage | Modèles entrées–sorties, réseaux, équilibrage | Hefferon Topics; EPA USEEIO |
| 6 | Espaces vectoriels, sous-espaces, indépendance, base, rang et les quatre sous-espaces fondamentaux | Sous-espaces, espace des colonnes, noyau, Ax = 0 | MIT 18.06 L5–7; ILA 2.6 |
| 7 | Espaces vectoriels, sous-espaces, indépendance, base, rang et les quatre sous-espaces fondamentaux | Indépendance, base, dimension, rang, quatre sous-espaces ; entraînement de mi-parcours | MIT 18.06 L8–10; ILA 2.5, 2.7, 2.9 |
| 8 | Orthogonalité, projections, moindres carrés et Gram–Schmidt | Orthogonalité et projections | MIT 18.06 L14–15; ILA 6.3 |
| 9 | Orthogonalité, projections, moindres carrés et Gram–Schmidt | Moindres carrés, Gram–Schmidt, A = QR ; TP D | MIT 18.06 L16–17; ILA 6.4–6.5 |
| 10 | Déterminants : propriétés, cofacteurs, règle de Cramer et volume | Déterminants, cofacteurs, règle de Cramer, volume | MIT 18.06 L18–20; ILA 4.1–4.3 |
| 11 | Valeurs propres et vecteurs propres | Valeurs propres, vecteurs propres, méthode de la puissance itérée, modes | MIT 18.06 L21; ILA 5.1–5.2, 5.5 |
| 12 | Diagonalisation, puissances de matrices, chaînes de Markov et systèmes dynamiques linéaires | Diagonalisation, chaînes de Markov, systèmes dynamiques ; TP E | MIT 18.06 L22–25; ILA 5.4, 5.6 |
| 13 | La décomposition en valeurs singulières, l’analyse en composantes principales et les matrices définies positives | Matrices définies positives et SVD | MIT 18.06 L27–29 |
| 14 | La décomposition en valeurs singulières, l’analyse en composantes principales et les matrices définies positives | Approximation de rang faible, ACP, conditionnement ; TP F | MIT 18.06 L29 ; site du livre de Strang |
| 15 | Révision | Révision et série d’exercices finale | Toutes les unités |
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.
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
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
- ‖v‖ = √(3² + 4²) = √25 = 5.
- 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
- v·w = 3 + 8 = 11.
- ‖v‖ = 5, ‖w‖ = √5 = 2.236.
- 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
- 2(1, 0, 1) + 3(0, 1, 1) = (2, 3, 5).
- 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.
Trouvez ‖(1, 2, 2)‖.
Afficher la réponse
√(1 + 4 + 4) = 3.
Trouvez (1, 2, 3)·(4, −5, 6).
Afficher la réponse
4 − 10 + 18 = 12.
(2, −1) et (1, 2) sont-ils orthogonaux ?
Afficher la réponse
Yes: 2 − 2 = 0.
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
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
- Interactive Linear Algebra, 2.1 Vectors: Texte gratuit avec des figures interactives.
- Interactive Linear Algebra, 6.1 Dot products and orthogonality: Normes, angles et orthogonalité.
- MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010): Cours 1 : la géométrie des équations linéaires.
- PhET : Vector Addition: Interactif.
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
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
- R₂ − 2R₁: −y − 3z = −11. R₃ − R₁: −2y + z = −1.
- R₃ − 2R₂ : 7z = 21, donc z = 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
- Équation 1 moins équation 2 : x² + y² − [(x − 10)² + y²] = 25 − 65, donc 20x − 100 = −40.
- Équation 1 moins équation 3 : 20y − 100 = −20.
- 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
- R₂ − 2R₁ donne 0 = 0 dans le premier système : un pivot, une variable libre.
- Dans le second, R₂ − 2R₁ donne 0 = 1.
Réponse : Une infinité ; aucune.
Exercices
Faites chacun d’abord sur papier, puis vérifiez.
Résolvez x + y = 5, x − y = 1.
Afficher la réponse
En additionnant, 2x = 6 : x = 3, y = 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.
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.
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
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
- Interactive Linear Algebra, 1.1 Systems of linear equations: Section de texte gratuit.
- Interactive Linear Algebra, 1.2 Row reduction: Formes échelonnées pas à pas.
- MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010): Cours 2 : l’élimination avec des matrices.
- Jim Hefferon, Linear Algebra (free text): Chapitre un : systèmes linéaires.
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
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
- AB : ligne 1 · colonnes = (2, 1) ; ligne 2 = (4, 3).
- 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
- ad − bc = 24 − 14 = 10.
- 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
- Multiplicateur ℓ₂₁ = 6/2 = 3.
- 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.
Quelle est la taille d’une matrice 3 × 2 fois une matrice 2 × 4 ?
Afficher la réponse
3 × 4.
Inversez [[2, 1], [5, 3]].
Afficher la réponse
det = 1, inverse [[3, −1], [−5, 2]].
Combien d’opérations environ l’élimination demande-t-elle pour n = 1000 ?
Afficher la réponse
n³/3 ≈ 3.33 × 10⁸.
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
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
- Interactive Linear Algebra, 3.4 Matrix multiplication: Section de texte gratuit.
- Interactive Linear Algebra, 3.5 Matrix inverses: Inverses et Gauss–Jordan.
- MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010): Cours 3–4 : multiplication, inverses et A = LU.
- MIT OpenCourseWare 18.06SC Linear Algebra (self-paced, Fall 2011): Vidéos de travaux dirigés sur LU.
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
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
- I − A = [[0.8, −0.3], [−0.1, 0.6]], det = 0.48 − 0.03 = 0.45.
- (I − A)⁻¹ = (1/0.45)[[0.6, 0.3], [0.1, 0.8]].
- 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
- 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
- C: x = z. H: 4x = 2w. O: 2y = 2z + w.
- 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.
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.
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.
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.
É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
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
- U.S. EPA: US Environmentally-Extended Input-Output (USEEIO) models: Les modèles entrées–sorties étendus à l’environnement de l’EPA.
- Jim Hefferon, Linear Algebra (free text): Thème : analyse entrées–sorties de Leontief.
- MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010): Cours et séries d’exercices d’applications.
- Interactive Linear Algebra, 1.2 Row reduction: Résoudre les systèmes.
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
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
- R₂ − 2R₁ = 0 ; R₃ − R₁ = [0, −2, −2]. Deux pivots, donc r = 2.
- Ax = 0 : x + z = 0 et −2y − 2z = 0, avec z libre.
- 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
- Échelonnez la matrice ayant ces colonnes : R₂ − 2R₁ et R₃ − 3R₁ donnent les lignes [0, −3, −6] et [0, −6, −12], qui sont proportionnelles.
- 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
- C(A) et C(Aᵀ) : r = 3.
- 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.
La droite y = 2x + 1 est-elle un sous-espace de ℝ² ?
Afficher la réponse
Non : elle ne contient pas (0, 0).
Une matrice 4 × 4 a un noyau de dimension 2. Quel est son rang ?
Afficher la réponse
4 − 2 = 2.
Quelle est la dimension de l’espace des matrices 2 × 2 ?
Afficher la réponse
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
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
- Interactive Linear Algebra, 2.6 Subspaces: Section de texte gratuit.
- Interactive Linear Algebra, 2.5 Linear independence: Tests d’indépendance.
- Interactive Linear Algebra, 2.9 The rank theorem: Rang et dimension du noyau.
- MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010): Cours 5–10 : espaces, rang et les quatre sous-espaces.
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
Exemples résolus pas à pas
Exemple 6.1 Projetez b = (2, 3) sur a = (1, 1).
Afficher les étapes
- a·b = 5, a·a = 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
- A a pour lignes (1, 0), (1, 1), (1, 2) ; AᵀA = [[3, 3], [3, 5]], Aᵀb = (7, 10).
- 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
- q₁ = (1, 1, 0)/√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.
Projetez (3, 4) sur (1, 0).
Afficher la réponse
(3, 0).
É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.
Q = [[0, 1], [1, 0]] est-elle orthogonale ?
Afficher la réponse
Yes: QᵀQ = I.
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
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
- Interactive Linear Algebra, 6.3 Orthogonal projection: Section de texte gratuit.
- Interactive Linear Algebra, 6.5 The method of least squares: Les moindres carrés en images.
- Interactive Linear Algebra, 6.4 Orthogonal sets: Gram–Schmidt.
- MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010): Cours 14–17 : orthogonalité, projections, moindres carrés, Gram–Schmidt.
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
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 − 3 × 0) = 1.
- − 2 × (0 × 1 − 3 × 2) = −2 × (−6) = 12.
- + 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
- det A = 6 − 1 = 5.
- x = det[[5, 1], [10, 3]]/5 = (15 − 10)/5.
- 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
- 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.
det [[3, 1], [2, 4]] ?
Afficher la réponse
12 − 2 = 10.
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.
Si det A = 4 pour une matrice A 3 × 3, que vaut det(2A) ?
Afficher la réponse
2³ × 4 = 32.
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
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
- Interactive Linear Algebra, 4.1 Determinants: definition: Section de texte gratuit.
- Interactive Linear Algebra, 4.3 Determinants and volumes: Sens géométrique.
- MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010): Cours 18–20 : déterminants, cofacteurs, règle de Cramer.
- Jim Hefferon, Linear Algebra (free text): Chapitre quatre : déterminants.
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
Exemples résolus pas à pas
Exemple 8.1 Trouvez les valeurs propres et vecteurs propres de A = [[2, 1], [1, 2]].
Afficher les étapes
- det(A − λI) = (2 − λ)² − 1 = λ² − 4λ + 3 = 0.
- λ = 3 : A − 3I = [[−1, 1], [1, −1]] donne v = (1, 1). λ = 1 : v = (1, −1).
- 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
- λ² − 7λ + 10 = 0 (trace 7, det 10).
- (λ − 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
- Valeurs propres de [[2, −1], [−1, 1]] : λ = (3 ± √5)/2 = 0.382 et 2.618.
- ω₁² = 100 × 0.382 = 38.2 s⁻², ω₁ = 6.18 rad/s.
- 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.
Valeurs propres de [[5, 0], [0, −2]] ?
Afficher la réponse
5 et −2 (diagonale).
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.
Valeurs propres de la rotation de 90° [[0, −1], [1, 0]] ?
Afficher la réponse
λ² + 1 = 0: λ = ±i.
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
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
- Interactive Linear Algebra, 5.1 Eigenvalues and eigenvectors: Section de texte gratuit.
- Interactive Linear Algebra, 5.2 The characteristic polynomial: Trouver des valeurs propres.
- Interactive Linear Algebra, 5.5 Complex eigenvalues: Rotation et homothétie.
- MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010): Cours 21 : valeurs propres et vecteurs propres.
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
Exemples résolus pas à pas
Exemple 9.1 Utilisez la diagonalisation pour trouver A¹⁰ pour A = [[2, 1], [1, 2]].
Afficher les étapes
- D’après le module 8 : λ = 3, 1 avec v = (1, 1), (1, −1). P = [[1, 1], [1, −1]], P⁻¹ = ½[[1, 1], [1, −1]].
- 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
- (P − I)q = 0 : −0.2q₁ + 0.3q₂ = 0, donc q₁ = 1.5q₂.
- 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
- x₁ = (0.8, 0.2); x₂ = (0.64 + 0.06, 0.16 + 0.14).
- 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.
[[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.
Trouvez son état stationnaire.
Afficher la réponse
0.5q₁ = 0.2q₂, donc q = (2/7, 5/7) = (0.286, 0.714).
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ᵏ.
[[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
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
- Interactive Linear Algebra, 5.4 Diagonalization: Section de texte gratuit.
- Interactive Linear Algebra, 5.6 Stochastic matrices: Chaînes de Markov et états stationnaires.
- MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010): Cours 22–25 : diagonalisation, matrices de Markov, matrices symétriques.
- MIT OpenCourseWare 18.06SC Linear Algebra (self-paced, Fall 2011): Travaux dirigés sur les matrices de Markov.
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
Exemples résolus pas à pas
Exemple 10.1 Trouvez les valeurs singulières de A = [[3, 0], [4, 5]].
Afficher les étapes
- AᵀA = [[25, 20], [20, 25]], valeurs propres 45 et 5.
- σ₁ = √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
- σ₁²/(σ₁² + σ₂²) = 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
- Chaque terme demande un u (1000), un v (1000) et un σ : 2001 nombres.
- 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.
[[2, 1], [1, 2]] est-elle définie positive ?
Afficher la réponse
Oui : les valeurs propres 3 et 1 sont positives.
Trouvez le conditionnement de diag(10, 0.1).
Afficher la réponse
10/0.1 = 100.
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).
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
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
- MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010): Cours 27–29 : matrices définies positives et SVD.
- MIT OpenCourseWare 18.06SC Linear Algebra (self-paced, Fall 2011): Séance sur la décomposition en valeurs singulières.
- Gilbert Strang, Introduction to Linear Algebra (book site): Notes de Strang sur la SVD et les données.
- Jim Hefferon, Linear Algebra (free text): Thèmes sur les matrices orthogonales et les projections.
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.
Mode opératoire
- 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).
- Calculez la longueur de chaque pixel avec SQRT(SUMSQ(…)) et le produit scalaire de chaque paire avec SUMPRODUCT.
- Trouvez l’angle entre chaque paire avec ACOS et DEGREES.
- 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
| Pixel | Vert | Red | NIR | NDVI |
|---|---|---|---|---|
| Végétation | 0.10 | 0.08 | 0.40 | 0.667 |
| Sol | 0.20 | 0.25 | 0.30 | 0.091 |
| Eau | 0.06 | 0.04 | 0.02 | −0.333 |
Réflectances d’exemple choisies pour ressembler à une végétation, un sol nu et une eau typiques.
Exemple de calcul
- NDVI végétation = (0.40 − 0.08)/(0.40 + 0.08) = 0.667.
- NDVI sol = 0.05/0.55 = 0.091 ; NDVI eau = −0.02/0.06 = −0.333.
- 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.
Mode opératoire
- 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.
- 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.
- Tracez des cercles avec une ficelle et un crayon et marquez leur point d’intersection.
- 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
| Balise | Position (cm) | Distance (cm) | d² |
|---|---|---|---|
| 1 | (0, 0) | 5.000 | 25 |
| 2 | (10, 0) | 8.062 | 65 |
| 3 | (0, 10) | 6.708 | 45 |
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
- Eq. 1 − Eq. 2: 20x − 100 = −40, x = 3. Eq. 1 − Eq. 3: 20y − 100 = −20, y = 4.
- 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|.
Mode opératoire
- Tracez le carré unité de sommets (0, 0), (1, 0), (1, 1), (0, 1).
- 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).
- Tracez chaque image et trouvez son aire en comptant les carreaux ou avec la formule du lacet.
- 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
| Matrice | det | Aire de l’image | Orientation |
|---|---|---|---|
| A = [[2, 1], [0, 1]] | 2 | 2 | conservée |
| R (rotation de 90°) | 1 | 1 | conservée |
| F (réflexion) | −1 | 1 | inversée |
| AR | 2 | 2 | conservée |
Valeurs exactes pour le carré unité.
Exemple de calcul
- det A = 2 × 1 − 1 × 0 = 2 ; l’image est un parallélogramme de base 2 et de hauteur 1.
- 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.
Mode opératoire
- 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.
- 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.
- Vérifiez avec LINEST ou une courbe de tendance de graphique. Tracez les résidus b − Ax̂.
- 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) | 0 | 8 | 16 | 24 | 32 |
|---|---|---|---|---|---|
| h (mm) | 0 | 22 | 48 | 78 | 114 |
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
- t̄ = 16, h̄ = 52.4; Σ(t − t̄)(h − h̄) = 2272, Σ(t − t̄)² = 640.
- D = 2272/640 = 3.55 mm/yr; C = 52.4 − 3.55 × 16 = −4.4 mm.
- 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.
Mode opératoire
- 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).
- Construisez P dans un tableur et partez de x₀ = (1, 0) : toute la matière en usage.
- Utilisez MMULT pour avancer de 10 ans, et comparez avec la formule (1 − f(1 − r))ᵏ pour la part encore en usage.
- 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 r | Probabilité de rester en usage | En usage après 10 ans | Durée moyenne d’usage (années) |
|---|---|---|---|
| 0.321 | 0.8642 | 0.232 | 7.36 |
| 0.50 | 0.9000 | 0.349 | 10.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
- Probabilité de rester = 1 − 0.2 × 0.679 = 0.8642 ; 0.8642¹⁰ = 0.232.
- Durée moyenne d’usage = 1/0.1358 = 7.36 ans ; avec r = 0.5 elle vaut 1/0.1 = 10.0 ans.
- λ = 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.
Mode opératoire
- 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.
- 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.
- Calculez la SVD (numpy.linalg.svd) et reconstruisez l’image avec k = 1, 2 et 4 termes.
- 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
| k | Nombres stockés (image 8 × 8) | Part de 64 |
|---|---|---|
| 1 | 17 | 26.6% |
| 2 | 34 | 53.1% |
| 4 | 68 | 106.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
- 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.
- 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
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)
Trouvez l’angle entre (1, 0, 1) et (0, 1, 1).
Afficher la solution complète
- Produit scalaire = 1 ; longueurs √2 chacune ; cos θ = 1/2.
Réponse : 60°.
Résolvez x + 2y = 4, 3x − y = 5.
Afficher la solution complète
- R₂ − 3R₁: −7y = −7.
Réponse : x = 2, y = 1.
Inversez [[3, 2], [1, 1]].
Afficher la solution complète
- det = 1.
Réponse : [[1, −2], [−1, 3]].
Une économie à un seul secteur consomme 0.4 de sa propre production. Quelle production satisfait une demande de 120 ?
Afficher la solution complète
- x = 120/0.6.
Réponse : 200.
Une matrice 6 × 9 est de rang 4. Trouvez dim N(A) et dim N(Aᵀ).
Afficher la solution complète
- n − r et m − r.
Réponse : 5 et 2.
Série d’entraînement finale (toutes les unités)
Ajustez y = C + Dt à (0, 0), (1, 2), (2, 2) par moindres carrés.
Afficher la solution complète
- AᵀA = [[3, 3], [3, 5]], Aᵀb = (4, 6) ; det 6.
Réponse : C = 0.333, D = 1.
Trouvez det [[2, 1, 0], [1, 2, 1], [0, 1, 2]].
Afficher la solution complète
- 2(4 − 1) − 1(2 − 0) + 0.
Réponse : 4.
Trouvez les valeurs propres de [[6, 2], [2, 3]].
Afficher la solution complète
- Trace 9, det 14 : λ² − 9λ + 14 = 0.
Réponse : 7 et 2.
Trouvez l’état stationnaire de la matrice de Markov [[0.9, 0.2], [0.1, 0.8]].
Afficher la solution complète
- 0.1q₁ = 0.2q₂.
Réponse : (0.667, 0.333).
Trouvez les valeurs singulières de [[2, 0], [0, −3]].
Afficher la solution complète
- AᵀA = diag(4, 9).
Réponse : 3 et 2.
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
- 20(500 + 800 + 1).
Réponse : 26020, environ 6.5 %.
Fiche de formules
Constantes
| Constante | Valeur |
|---|---|
| Vitesse de la lumière (utilisée pour les distances GPS) | 299,792,458 m/s |
| 1 mile | 1.609344 km |
| 1 pied | 0.3048 m |
| 1 pouce | 25.4 mm |
| Flops pour résoudre un système n × n par élimination | environ n³/3 multiplications–additions (Strang) |
| Matrice identité I | Des uns sur la diagonale, des zéros ailleurs : AI = IA = A |
Source : NIST Guide to the SI, Appendix B.9: conversion factors
Textes et cours gratuits
Utilisez ces manuels, cours et outils gratuits en parallèle du cours.
- MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010): Les cours vidéo complets, séries d’exercices et examens de Strang.
- MIT OpenCourseWare 18.06SC Linear Algebra (self-paced, Fall 2011): Version en autonomie avec vidéos de travaux dirigés et problèmes corrigés.
- Interactive Linear Algebra (Margalit and Rabinoff, Georgia Tech): Manuel interactif gratuit utilisé à Georgia Tech.
- Jim Hefferon, Linear Algebra (free text): Manuel complet gratuit avec réponses.
- Gilbert Strang, Introduction to Linear Algebra (book site): Le site du livre de Strang avec notes et figures.
- PhET : Vector Addition: Vecteurs et composantes interactifs.
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.
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îtrise | Choisit 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étent | Résout les exercices d’entraînement avec quelques erreurs de calcul occasionnelles. |
| 2 · En progression | Suit les exemples corrigés mais confond lignes et colonnes ou les quatre sous-espaces. |
| 1 · Débutant | Manipule des symboles sans vérifier dimensions ni réponses ; doit revoir le module. |
Pondération suggérée si vous voulez une note
| Partie | Poids |
|---|---|
| Exercices | 25% |
| Travaux pratiques | 20% |
| Série d’entraînement de mi-parcours | 25% |
| Série d’entraînement finale | 30% |
Visez le niveau 3 ou mieux dans chaque unité avant de passer au cours suivant.



