Algèbre linéaire 201-NYC / Maths 105 • Complément québécois de Première et cégep à Montréal

Exercices corrigés : systèmes d'équations linéaires, Gauss et Gauss-Jordan (201-NYC / Maths 105)

Voici la série d'exercices corrigés d'algèbre linéaire et géométrie vectorielle (cours 201-NYC, aussi appelé Maths 105 selon le programme visé) sur les systèmes d'équations linéaires. La partie A couvre les bases : matrice augmentée et opérations élémentaires sur les lignes, méthode de Gauss et résolution par remontée, passage à la forme échelonnée réduite par Gauss-Jordan, les trois cas possibles avec paramétrage d'un ensemble infini de solutions, et l'inversion d'une matrice par la méthode de Gauss-Jordan. La partie B monte au niveau examen : systèmes homogènes et solutions non triviales, et une discussion complète d'un système dépendant de deux paramètres.

C'est le chapitre le plus lourd de la session en volume de calcul, et paradoxalement celui où l'on perd le plus de points sur des questions qui n'en demandent aucun. Un système linéaire ne se résout pas seulement, il se DISCUTE : combien de solutions, pourquoi, et sous quelles conditions. Un étudiant qui échelonne parfaitement mais conclut « pas de solution » là où il y en a une infinité perd toute la question.

Le réflexe à installer : ne jamais mélanger les opérations élémentaires sur les LIGNES avec des manipulations de colonnes, et écrire l'opération utilisée à chaque étape. Les trois opérations autorisées sont l'échange de deux lignes, la multiplication d'une ligne par un scalaire non nul, et l'ajout à une ligne d'un multiple d'une autre ligne. Toute autre manipulation change l'ensemble des solutions, et l'erreur devient alors impossible à retrouver dans une copie de deux pages.

Série autocorrigée Tape tes réponses sous chaque question : la page te dit juste ou faux avant d'ouvrir la correction. Avec un compte, chaque bonne réponse du premier coup rapporte des points.

Ce chapitre fait partie de Algèbre linéaire et géométrie vectorielle, 201-NYC
Avant de commencer Fiche de révision : les pièges et la méthode de ce chapitre

Avant ce chapitre

Ces notions sont supposées acquises ici. Si le premier exercice résiste, le blocage vient presque toujours de l'une d'elles, pas du chapitre lui-même.

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

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

  1. 1Équations et inéquations du premier degréSecondaire 3
  2. 2Systèmes d'équations : comparaison et graphiqueSecondaire 3
  3. 3La factorisation et les expressions rationnellesSecondaire 4 SN4
  4. 4Les équations et les inéquationsSecondaire 4 SN4
  5. 5Les systèmes d'équationsSecondaire 4 SN4
  6. 6Relations métriques et trigonométrie du triangle rectangleSecondaire 4 SN4
  7. 7Les vecteurs et la géométrie vectorielleSecondaire 5 SN5
  8. 8Combinaisons linéaires, dépendance et bases
  9. 9Matrices, systèmes et déterminants

Rappel de cours

  • Matrice augmentée : pour le système AX=BAX=B, c'est la matrice [AB][A\,|\,B] obtenue en accolant la colonne des seconds membres à la matrice des coefficients. La barre verticale n'a aucun rôle calculatoire, elle sépare visuellement les coefficients du membre de droite.
  • Opérations élémentaires sur les lignes, les trois seules autorisées : échanger deux lignes (LiLjL_{i}\leftrightarrow L_{j}) ; multiplier une ligne par un scalaire NON NUL (LikLiL_{i}\leftarrow kL_{i}, k0k\neq 0) ; ajouter à une ligne un multiple d'une autre ligne (LiLi+kLjL_{i}\leftarrow L_{i}+kL_{j}). Chacune transforme le système en un système ÉQUIVALENT, c'est-à-dire ayant exactement le même ensemble de solutions.
  • Forme échelonnée (méthode de Gauss) : dans chaque ligne non nulle, le premier coefficient non nul, appelé pivot, est situé strictement à droite du pivot de la ligne précédente, et les lignes entièrement nulles sont en bas. On résout ensuite par REMONTÉE, de la dernière équation vers la première.
  • Forme échelonnée réduite (méthode de Gauss-Jordan) : en plus des conditions ci-dessus, chaque pivot vaut 11 et il est le SEUL coefficient non nul de sa colonne, donc il y a aussi des zéros AU-DESSUS de chaque pivot. La solution se lit alors directement, sans aucune remontée.
  • Les trois cas possibles pour un système linéaire, et il n'y en a jamais d'autres. Une solution unique : autant de pivots que d'inconnues. Aucune solution (système incompatible) : apparition d'une ligne du type [0 0  0c][0\ 0\ \dots\ 0\,|\,c] avec c0c\neq 0, qui traduit l'équation impossible 0=c0=c. Une infinité de solutions : le système est compatible et il y a moins de pivots que d'inconnues, les inconnues sans pivot devenant des paramètres libres.
  • Nombre de paramètres libres == nombre d'inconnues - nombre de pivots. C'est ce compte qui décide, jamais le nombre d'équations de départ.
  • Inversion par Gauss-Jordan : on forme la matrice [AI][A\,|\,I] et on applique les opérations élémentaires jusqu'à obtenir [IA1][I\,|\,A^{-1}]. Si au cours du calcul une ligne de la partie gauche devient entièrement nulle, alors AA n'est pas inversible et le procédé s'arrête.
  • Système homogène : un système AX=0AX=0, dont le second membre est nul. Il admet TOUJOURS la solution triviale X=0X=0, donc il n'est jamais incompatible. Il admet des solutions non triviales si et seulement si det(A)=0\det(A)=0.

Partie A : Les bases (/50)

Exercice 1 : Matrice augmentée, opérations élémentaires et méthode de Gauss

On considère le système {x+2yz=22x+5y+z=153x+4y7z=10\begin{cases} x+2y-z=2 \\ 2x+5y+z=15 \\ 3x+4y-7z=-10 \end{cases}.

  • a) Écrivez la matrice augmentée [AB][A\,|\,B] de ce système.
  • b) Effectuez les opérations L2L22L1L_{2}\leftarrow L_{2}-2L_{1} et L3L33L1L_{3}\leftarrow L_{3}-3L_{1}, et écrivez la matrice obtenue.
  • c) Terminez l'échelonnement pour obtenir une forme échelonnée, en indiquant l'opération utilisée.
  • d) Résolvez le système par remontée, puis vérifiez votre solution dans les trois équations d'origine.

Tape tes réponses, la page te dit juste ou faux 0/6

b)
c)
d)
Voir la correction

Réponses

  • a) Matrice augmentée 3×43\times 4
  • b) (0;1;311)(0;1;3\,|\,11) et (0;2;416)(0;-2;-4\,|\,-16)
  • c) L3L3+2L2L_3\leftarrow L_3+2L_2 : (0;0;26)(0;0;2\,|\,6)
  • d) (1; 2; 3)(1;\ 2;\ 3)

a) [AB]=[12122511534710][A\,|\,B]=\left[\begin{array}{ccc|c} 1 & 2 & -1 & 2 \\ 2 & 5 & 1 & 15 \\ 3 & 4 & -7 & -10 \end{array}\right]. Le coefficient en haut à gauche vaut 11, ce qui en fait un pivot idéal : aucune fraction ne sera introduite à la première étape.

b) L2L22L1L_{2}\leftarrow L_{2}-2L_{1} donne (22; 54; 1+2154)=(0; 1; 311)(2-2;\ 5-4;\ 1+2\,|\,15-4)=(0;\ 1;\ 3\,|\,11). Puis L3L33L1L_{3}\leftarrow L_{3}-3L_{1} donne (33; 46; 7+3106)=(0; 2; 416)(3-3;\ 4-6;\ -7+3\,|\,-10-6)=(0;\ -2;\ -4\,|\,-16).

La matrice devient [12120131102416]\left[\begin{array}{ccc|c} 1 & 2 & -1 & 2 \\ 0 & 1 & 3 & 11 \\ 0 & -2 & -4 & -16 \end{array}\right]. Attention au signe sur la troisième colonne de L2L_{2} : 12(1)=1+2=31-2(-1)=1+2=3. Soustraire un multiple d'un coefficient négatif produit une addition, et c'est là que la majorité des erreurs d'échelonnement se glissent.

c) Le pivot de la deuxième ligne vaut 11, on l'utilise pour annuler le 2-2 situé en dessous : L3L3+2L2L_{3}\leftarrow L_{3}+2L_{2} donne (0; 2+2; 4+616+22)=(0; 0; 26)(0;\ -2+2;\ -4+6\,|\,-16+22)=(0;\ 0;\ 2\,|\,6).

La forme échelonnée est [1212013110026]\left[\begin{array}{ccc|c} 1 & 2 & -1 & 2 \\ 0 & 1 & 3 & 11 \\ 0 & 0 & 2 & 6 \end{array}\right]. Les trois pivots sont 11, 11 et 22, chacun strictement à droite du précédent : la forme est bien échelonnée. Trois pivots pour trois inconnues annoncent déjà une solution unique, avant même de la calculer.

d) On remonte de la dernière ligne vers la première. Troisième ligne : 2z=62z=6, donc z=3z=3. Deuxième ligne : y+3z=11y+3z=11, soit y+9=11y+9=11, donc y=2y=2. Première ligne : x+2yz=2x+2y-z=2, soit x+43=2x+4-3=2, donc x=1x=1.

Vérification dans les trois équations d'origine, et non dans les lignes échelonnées, qui pourraient contenir une erreur propagée. Équation 1 : 1+2(2)3=1+43=21+2(2)-3=1+4-3=2, conforme. Équation 2 : 2(1)+5(2)+3=2+10+3=152(1)+5(2)+3=2+10+3=15, conforme. Équation 3 : 3(1)+4(2)7(3)=3+821=103(1)+4(2)-7(3)=3+8-21=-10, conforme. La solution est (x; y; z)=(1; 2; 3)(x;\ y;\ z)=(1;\ 2;\ 3).

Exercice 2 : De Gauss à Gauss-Jordan : la forme échelonnée réduite

On reprend la forme échelonnée obtenue à l'exercice précédent : [1212013110026]\left[\begin{array}{ccc|c} 1 & 2 & -1 & 2 \\ 0 & 1 & 3 & 11 \\ 0 & 0 & 2 & 6 \end{array}\right].

  • a) Rappelez les deux conditions supplémentaires qu'une forme échelonnée RÉDUITE doit vérifier, et dites lesquelles manquent ici.
  • b) Poursuivez l'élimination jusqu'à la forme échelonnée réduite, en indiquant chaque opération.
  • c) Lisez la solution du système sur la matrice obtenue.
  • d) Dans quel cas la méthode de Gauss suffit-elle, et dans quel cas Gauss-Jordan est-elle nettement préférable ? Justifiez par le travail réellement économisé.

Tape tes réponses, la page te dit juste ou faux 0/6

a)
b)
c)
d)
Voir la correction

Réponses

  • a) Pivots à 1 et zéros au-dessus
  • b) [I(1;2;3)][I\,|\,(1;2;3)]
  • c) x=1x=1, y=2y=2, z=3z=3
  • d) Gauss pour un système ; Gauss-Jordan pour inverser ou paramétrer

a) Une forme échelonnée réduite exige deux conditions de plus : chaque pivot doit valoir exactement 11, et chaque pivot doit être le SEUL coefficient non nul de sa colonne, donc il faut aussi des zéros au-dessus des pivots.

Ici, les deux conditions manquent. Le pivot de la troisième ligne vaut 22 et non 11. Et au-dessus des pivots subsistent des coefficients non nuls : le 22 en position (1; 2)(1;\ 2) au-dessus du pivot de la deuxième colonne, ainsi que le 1-1 et le 33 au-dessus du pivot de la troisième colonne.

b) On normalise d'abord le dernier pivot : L312L3L_{3}\leftarrow\dfrac{1}{2}L_{3} donne (0; 0; 13)(0;\ 0;\ 1\,|\,3).

On remonte ensuite colonne par colonne. L2L23L3L_{2}\leftarrow L_{2}-3L_{3} donne (0; 1; 33119)=(0; 1; 02)(0;\ 1;\ 3-3\,|\,11-9)=(0;\ 1;\ 0\,|\,2). Puis L1L1+L3L_{1}\leftarrow L_{1}+L_{3} donne (1; 2; 1+12+3)=(1; 2; 05)(1;\ 2;\ -1+1\,|\,2+3)=(1;\ 2;\ 0\,|\,5).

Il reste à nettoyer la deuxième colonne : L1L12L2L_{1}\leftarrow L_{1}-2L_{2} donne (1; 22; 054)=(1; 0; 01)(1;\ 2-2;\ 0\,|\,5-4)=(1;\ 0;\ 0\,|\,1).

La forme échelonnée réduite est [100101020013]\left[\begin{array}{ccc|c} 1 & 0 & 0 & 1 \\ 0 & 1 & 0 & 2 \\ 0 & 0 & 1 & 3 \end{array}\right]. Ordre de travail à retenir : on traite les colonnes de DROITE à GAUCHE. Nettoyer la deuxième colonne avant la troisième obligerait à y revenir, puisque les opérations sur la troisième colonne réintroduisent des termes.

c) La lecture est immédiate : la première ligne dit x=1x=1, la deuxième y=2y=2, la troisième z=3z=3. Aucune remontée n'est nécessaire, la matrice EST la solution. On retrouve bien (1; 2; 3)(1;\ 2;\ 3), ce qui confirme l'exercice précédent par un chemin différent.

d) Gauss suffit largement quand on cherche la solution d'un seul système à solution unique : la remontée est rapide, et pousser jusqu'à la forme réduite ajoute du travail pour un gain nul. C'est d'ailleurs la méthode la plus efficace en temps d'examen sur un système 3×33\times 3 ordinaire.

Gauss-Jordan devient nettement préférable dans trois situations. D'abord lorsqu'il faut inverser une matrice, puisque le procédé [AI][IA1][A\,|\,I]\to[I\,|\,A^{-1}] exige la forme réduite. Ensuite lorsque le système admet une infinité de solutions : la forme réduite donne directement l'expression de chaque inconnue principale en fonction des paramètres libres, là où la remontée oblige à une substitution pénible et propice aux erreurs de signe. Enfin lorsqu'un même système doit être résolu pour plusieurs seconds membres, que l'on peut traiter simultanément en accolant plusieurs colonnes.

Autrement dit, le surcoût de Gauss-Jordan est un investissement : il ne se rentabilise que si l'on a besoin de la structure complète de la solution, pas seulement d'un triplet de nombres.

Exercice 3 : Les trois cas : solution unique, aucune solution, une infinité

Un système linéaire n'a jamais que trois issues possibles. Cet exercice les fait apparaître toutes les trois, et demande à chaque fois de les reconnaître SUR la forme échelonnée.

  • a) Résolvez {x+2y+z=32x+4y+2z=7\begin{cases} x+2y+z=3 \\ 2x+4y+2z=7 \end{cases} et concluez.
  • b) Résolvez {x+y+z=6x+2y+3z=142x+3y+4z=20\begin{cases} x+y+z=6 \\ x+2y+3z=14 \\ 2x+3y+4z=20 \end{cases} et donnez l'ensemble des solutions sous forme paramétrique.
  • c) Dans le cas b), vérifiez votre paramétrage pour deux valeurs distinctes du paramètre.
  • d) Interprétez géométriquement les trois cas possibles pour un système de trois équations à trois inconnues.

Tape tes réponses, la page te dit juste ou faux 0/6

a)
b)
c)
d)
Voir la correction

Réponses

  • a) 0=10=1 : aucune solution
  • b) (t2; 82t; t)(t-2;\ 8-2t;\ t)
  • c) t=3t=3 : (1;2;3)(1;2;3) ; t=0t=0 : (2;8;0)(-2;8;0)
  • d) Point, rien, droite ou plan

a) Matrice augmentée [12132427]\left[\begin{array}{ccc|c} 1 & 2 & 1 & 3 \\ 2 & 4 & 2 & 7 \end{array}\right]. On applique L2L22L1L_{2}\leftarrow L_{2}-2L_{1} : (22; 44; 2276)=(0; 0; 01)(2-2;\ 4-4;\ 2-2\,|\,7-6)=(0;\ 0;\ 0\,|\,1).

La deuxième ligne s'écrit 0x+0y+0z=10x+0y+0z=1, c'est-à-dire 0=10=1, ce qui est impossible. Le système est INCOMPATIBLE, il n'admet AUCUNE solution.

Le signal à reconnaître : une ligne dont la partie gauche est entièrement nulle mais dont le second membre ne l'est pas. Notez que le membre de gauche de la deuxième équation est exactement le double de celui de la première, alors que 72×37\neq 2\times 3 : les deux équations se contredisent.

b) Matrice augmentée [11161231423420]\left[\begin{array}{ccc|c} 1 & 1 & 1 & 6 \\ 1 & 2 & 3 & 14 \\ 2 & 3 & 4 & 20 \end{array}\right]. On applique L2L2L1L_{2}\leftarrow L_{2}-L_{1}, ce qui donne (0; 1; 28)(0;\ 1;\ 2\,|\,8), puis L3L32L1L_{3}\leftarrow L_{3}-2L_{1}, ce qui donne (0; 1; 28)(0;\ 1;\ 2\,|\,8).

Les deux nouvelles lignes sont identiques, donc L3L3L2L_{3}\leftarrow L_{3}-L_{2} donne (0; 0; 00)(0;\ 0;\ 0\,|\,0). La forme échelonnée est [111601280000]\left[\begin{array}{ccc|c} 1 & 1 & 1 & 6 \\ 0 & 1 & 2 & 8 \\ 0 & 0 & 0 & 0 \end{array}\right].

Ici la ligne nulle a un second membre NUL, ce qui ne traduit aucune contradiction : c'est l'équation 0=00=0, toujours vraie, qui signale simplement une équation redondante. Le système est compatible. Il possède deux pivots (colonnes de xx et de yy) pour trois inconnues, donc 32=13-2=1 paramètre libre. L'inconnue zz, dont la colonne ne porte pas de pivot, devient le paramètre.

On pose z=tz=t. La deuxième ligne donne y+2t=8y+2t=8, donc y=82ty=8-2t. La première donne x+(82t)+t=6x+(8-2t)+t=6, soit x=68+2tt=t2x=6-8+2t-t=t-2.

L'ensemble des solutions est {(t2; 82t; t) : tR}\{(t-2;\ 8-2t;\ t)\ :\ t\in\mathbb{R}\}, une INFINITÉ de solutions. Il faut donner l'ensemble complet : répondre par une seule solution particulière, même correcte, ne vaut pas la question.

c) Pour t=3t=3 : (1; 2; 3)(1;\ 2;\ 3). Vérification : 1+2+3=61+2+3=6, conforme ; 1+4+9=141+4+9=14, conforme ; 2+6+12=202+6+12=20, conforme.

Pour t=0t=0 : (2; 8; 0)(-2;\ 8;\ 0). Vérification : 2+8+0=6-2+8+0=6, conforme ; 2+16+0=14-2+16+0=14, conforme ; 4+24+0=20-4+24+0=20, conforme. Deux valeurs distinctes du paramètre donnent deux solutions distinctes, toutes deux valides : le paramétrage est correct.

d) Chaque équation linéaire à trois inconnues représente un PLAN de l'espace, et résoudre le système revient à chercher l'intersection des trois plans.

Solution unique : les trois plans se coupent en un seul POINT, comme deux murs et un plafond dans un coin de pièce. Aucune solution : les plans n'ont aucun point commun aux trois, soit parce que deux d'entre eux sont parallèles distincts, soit parce qu'ils forment un prisme, chaque paire se coupant selon une droite mais les trois droites restant distinctes. C'est le cas a), où les deux plans sont parallèles et distincts. Une infinité de solutions : les trois plans se coupent selon une DROITE commune, ce qui est le cas b), la solution paramétrée par tt étant précisément l'équation paramétrique de cette droite. Une infinité peut aussi survenir si les trois plans sont confondus, auquel cas l'ensemble des solutions est un plan entier, à deux paramètres.

Ce pont avec la géométrie n'est pas décoratif : reconnaître que la solution de b) est une droite de vecteur directeur (1; 2; 1)(1;\ -2;\ 1), lu directement sur les coefficients de tt, relie ce chapitre à celui des droites et plans de l'espace.

solution UNIQUEelles se coupentAUCUNE solutionparallèlesune INFINITÉconfonduesun système 2 × 2 se lit comme deux droites du plan

Exercice 4 : Inversion d'une matrice par la méthode de Gauss-Jordan

On considère la matrice A=(111012124)A=\begin{pmatrix} 1 & 1 & 1 \\ 0 & 1 & 2 \\ 1 & 2 & 4 \end{pmatrix}.

On rappelle le principe : on forme [AI][A\,|\,I] et on applique les opérations élémentaires jusqu'à transformer la partie gauche en II ; la partie droite devient alors A1A^{-1}.

  • a) Écrivez la matrice [AI][A\,|\,I].
  • b) Appliquez la méthode de Gauss-Jordan pour obtenir [IA1][I\,|\,A^{-1}].
  • c) Vérifiez votre résultat en calculant AA1A\,A^{-1}.
  • d) Utilisez A1A^{-1} pour résoudre le système AX=BAX=B avec B=(6817)B=\begin{pmatrix} 6 \\ 8 \\ 17 \end{pmatrix}.

Tape tes réponses, la page te dit juste ou faux 0/7

b)
d)
Voir la correction

Réponses

  • a) [AI][A\,|\,I]
  • b) A1=(021232111)A^{-1}=\begin{pmatrix}0&-2&1\\2&3&-2\\-1&-1&1\end{pmatrix}
  • c) AA1=I3AA^{-1}=I_3
  • d) X=(1; 2; 3)X=(1;\ 2;\ 3)

a) [AI]=[111100012010124001][A\,|\,I]=\left[\begin{array}{ccc|ccc} 1 & 1 & 1 & 1 & 0 & 0 \\ 0 & 1 & 2 & 0 & 1 & 0 \\ 1 & 2 & 4 & 0 & 0 & 1 \end{array}\right]. Toutes les opérations qui suivent portent sur les lignes ENTIÈRES, partie droite comprise : oublier d'appliquer une opération à droite est l'erreur qui ruine le calcul.

b) On élimine sous le premier pivot. L3L3L1L_{3}\leftarrow L_{3}-L_{1} donne (0; 1; 31; 0; 1)(0;\ 1;\ 3\,|\,-1;\ 0;\ 1), ce qui produit [111100012010013101]\left[\begin{array}{ccc|ccc} 1 & 1 & 1 & 1 & 0 & 0 \\ 0 & 1 & 2 & 0 & 1 & 0 \\ 0 & 1 & 3 & -1 & 0 & 1 \end{array}\right].

On élimine sous le deuxième pivot : L3L3L2L_{3}\leftarrow L_{3}-L_{2} donne (0; 0; 11; 1; 1)(0;\ 0;\ 1\,|\,-1;\ -1;\ 1). La partie gauche est maintenant échelonnée avec trois pivots, ce qui garantit déjà que AA est inversible.

On remonte. L2L22L3L_{2}\leftarrow L_{2}-2L_{3} donne (0; 1; 00+2; 1+2; 02)=(0; 1; 02; 3; 2)(0;\ 1;\ 0\,|\,0+2;\ 1+2;\ 0-2)=(0;\ 1;\ 0\,|\,2;\ 3;\ -2). Puis L1L1L3L_{1}\leftarrow L_{1}-L_{3} donne (1; 1; 01+1; 0+1; 01)=(1; 1; 02; 1; 1)(1;\ 1;\ 0\,|\,1+1;\ 0+1;\ 0-1)=(1;\ 1;\ 0\,|\,2;\ 1;\ -1).

Enfin L1L1L2L_{1}\leftarrow L_{1}-L_{2} donne (1; 0; 022; 13; 1+2)=(1; 0; 00; 2; 1)(1;\ 0;\ 0\,|\,2-2;\ 1-3;\ -1+2)=(1;\ 0;\ 0\,|\,0;\ -2;\ 1).

Le tableau final est [100021010232001111]\left[\begin{array}{ccc|ccc} 1 & 0 & 0 & 0 & -2 & 1 \\ 0 & 1 & 0 & 2 & 3 & -2 \\ 0 & 0 & 1 & -1 & -1 & 1 \end{array}\right], donc A1=(021232111)A^{-1}=\begin{pmatrix} 0 & -2 & 1 \\ 2 & 3 & -2 \\ -1 & -1 & 1 \end{pmatrix}.

c) Première ligne de AA1A\,A^{-1} : (1; 1; 1)(1;\ 1;\ 1) contre la première colonne (0; 2; 1)(0;\ 2;\ -1) donne 0+21=10+2-1=1 ; contre la deuxième colonne (2; 3; 1)(-2;\ 3;\ -1) donne 2+31=0-2+3-1=0 ; contre la troisième (1; 2; 1)(1;\ -2;\ 1) donne 12+1=01-2+1=0. On obtient (1; 0; 0)(1;\ 0;\ 0).

Deuxième ligne (0; 1; 2)(0;\ 1;\ 2) : 0+22=00+2-2=0, puis 0+32=10+3-2=1, puis 02+2=00-2+2=0, soit (0; 1; 0)(0;\ 1;\ 0). Troisième ligne (1; 2; 4)(1;\ 2;\ 4) : 0+44=00+4-4=0, puis 2+64=0-2+6-4=0, puis 14+4=11-4+4=1, soit (0; 0; 1)(0;\ 0;\ 1). Le produit vaut I3I_{3}, l'inverse est validé.

d) La solution est X=A1BX=A^{-1}B, et non BA1BA^{-1} : l'ordre des facteurs n'est pas libre, et le produit BA1BA^{-1} n'est d'ailleurs même pas défini pour des dimensions 3×13\times 1 et 3×33\times 3.

X=(021232111)(6817)X=\begin{pmatrix} 0 & -2 & 1 \\ 2 & 3 & -2 \\ -1 & -1 & 1 \end{pmatrix}\begin{pmatrix} 6 \\ 8 \\ 17 \end{pmatrix}. Première composante : 0(6)2(8)+1(17)=16+17=10(6)-2(8)+1(17)=-16+17=1. Deuxième : 2(6)+3(8)2(17)=12+2434=22(6)+3(8)-2(17)=12+24-34=2. Troisième : 1(6)1(8)+1(17)=68+17=3-1(6)-1(8)+1(17)=-6-8+17=3. Donc X=(123)X=\begin{pmatrix} 1 \\ 2 \\ 3 \end{pmatrix}.

Vérification directe dans le système : 1+2+3=61+2+3=6, conforme ; 0+2+6=80+2+6=8, conforme ; 1+4+12=171+4+12=17, conforme.

Une remarque de méthode : pour résoudre UN système, calculer l'inverse est un détour coûteux, et Gauss reste plus rapide. L'inverse devient rentable quand plusieurs seconds membres BB doivent être traités avec la même matrice AA, puisqu'il suffit alors d'un produit matriciel par second membre.

Exercice 5 : Le rang et le théorème de compatibilité

Le RANG d'une matrice est le nombre de pivots de sa forme échelonnée, c'est-à-dire le nombre d'équations réellement indépendantes. Cette unique quantité décide à elle seule du nombre de solutions d'un système, et le résume dans un théorème.

  • a) Déterminez le rang de la matrice des coefficients du système {x+2yz=22x+5y+z=153x+4y7z=10\begin{cases} x+2y-z=2 \\ 2x+5y+z=15 \\ 3x+4y-7z=-10 \end{cases}, et concluez sur le nombre de solutions sans les calculer.
  • b) Énoncez le théorème de compatibilité (Rouché) : un système AX=BAX=B est compatible si et seulement si rang(A)=rang([AB])\operatorname{rang}(A)=\operatorname{rang}([A\,|\,B]). Illustrez le cas d'incompatibilité sur {x+2y+z=32x+4y+2z=7\begin{cases} x+2y+z=3 \\ 2x+4y+2z=7 \end{cases} en comparant les deux rangs.
  • c) Pour un système compatible à nn inconnues, exprimez le nombre de paramètres libres à l'aide du rang, et retrouvez ainsi le critère de solution unique.
  • d) Un système de 55 équations à 33 inconnues peut-il avoir une solution unique ? Une infinité ? Aucune ? Discutez chaque cas à l'aide du rang, en soulignant que c'est lui, et non le nombre d'équations, qui décide.

Tape tes réponses, la page te dit juste ou faux 0/5

a)
b)
c)
d)
Voir la correction

Réponses

  • a) Rang 3 : solution unique
  • b) Rangs 1 et 2 : incompatible
  • c) nrang(A)n-\operatorname{rang}(A) paramètres
  • d) Les trois cas possibles

a) On échelonne la matrice des coefficients. Avec L2L22L1L_{2}\leftarrow L_{2}-2L_{1} et L3L33L1L_{3}\leftarrow L_{3}-3L_{1}, puis L3L3+2L2L_{3}\leftarrow L_{3}+2L_{2}, on obtient la forme échelonnée [121013002]\left[\begin{array}{ccc} 1 & 2 & -1 \\ 0 & 1 & 3 \\ 0 & 0 & 2 \end{array}\right] (le calcul est celui de l'exercice 1). Elle possède trois pivots, donc rang(A)=3\operatorname{rang}(A)=3. Comme le rang égale le nombre d'inconnues, le système a une SOLUTION UNIQUE, et cela se sait avant même de la calculer. On peut confirmer par le déterminant : det(A)=20\det(A)=2\neq 0, ce qui est cohérent, un déterminant non nul signalant toujours un rang maximal.

b) Le théorème de compatibilité affirme qu'un système AX=BAX=B admet au moins une solution si et seulement si le rang de la matrice des coefficients AA est égal au rang de la matrice augmentée [AB][A\,|\,B]. Autrement dit, adjoindre la colonne des seconds membres ne doit pas faire apparaître un pivot supplémentaire. Sur le système donné, la matrice augmentée est [12132427]\left[\begin{array}{ccc|c} 1 & 2 & 1 & 3 \\ 2 & 4 & 2 & 7 \end{array}\right]. L'opération L2L22L1L_{2}\leftarrow L_{2}-2L_{1} donne [0 0 01][0\ 0\ 0\,|\,1]. La partie gauche a alors un seul pivot, donc rang(A)=1\operatorname{rang}(A)=1, tandis que la partie augmentée porte un pivot dans la colonne des seconds membres, donc rang([AB])=2\operatorname{rang}([A\,|\,B])=2. Les deux rangs diffèrent, le système est donc INCOMPATIBLE, ce qui traduit la ligne impossible 0=10=1. Le théorème formalise exactement le signal qu'on repère à l'échelonnement.

c) Pour un système compatible à nn inconnues, le nombre de paramètres libres est nrang(A)n-\operatorname{rang}(A) : chaque pivot fixe une inconnue principale, et les inconnues restantes, en nombre nrang(A)n-\operatorname{rang}(A), deviennent libres. Le critère de solution unique s'en déduit immédiatement : la solution est unique lorsqu'il n'y a AUCUN paramètre libre, c'est-à-dire lorsque rang(A)=n\operatorname{rang}(A)=n. En résumé, un système compatible a une solution unique si rang(A)=n\operatorname{rang}(A)=n, et une infinité si rang(A)<n\operatorname{rang}(A)<n, l'infinité étant de dimension nrang(A)n-\operatorname{rang}(A).

d) Les trois cas sont possibles, et c'est le rang qui les sépare, jamais le nombre d'équations. Le rang d'une matrice 5×35\times 3 ne peut pas dépasser 33, le plus petit des deux nombres de lignes et de colonnes. Solution unique : si rang(A)=rang([AB])=3\operatorname{rang}(A)=\operatorname{rang}([A\,|\,B])=3, les cinq équations se ramènent à trois équations indépendantes déterminant les trois inconnues ; les deux équations excédentaires sont des conséquences des autres, donc redondantes mais compatibles. Une infinité : si les deux rangs sont égaux mais valent moins de 33, disons 22, il reste 32=13-2=1 paramètre libre. Aucune solution : si rang(A)<rang([AB])\operatorname{rang}(A)<\operatorname{rang}([A\,|\,B]), ce qui arrive dès qu'une équation contredit les autres, le système est incompatible, situation la plus fréquente quand on impose cinq contraintes à seulement trois inconnues. La morale du chapitre est là : avoir plus d'équations que d'inconnues ne garantit rien, et avoir moins d'équations n'interdit pas l'incompatibilité ; seul le rang tranche.

COMPATIBLEun point commun aux troisINCOMPATIBLEtrois intersections distinctestrois équations à deux inconnues : une de trop

Partie B : Niveau examen (/50)

Exercice 6 : Systèmes homogènes et solutions non triviales

Un système homogène est un système AX=0AX=0, dont tous les seconds membres sont nuls. Cet exercice établit ses propriétés, puis les applique.

  • a) Démontrez qu'un système homogène n'est jamais incompatible, et précisez ce qu'on appelle sa solution triviale.
  • b) Énoncez et justifiez la condition sur det(A)\det(A) qui garantit l'existence de solutions NON triviales.
  • c) Déterminez toutes les valeurs de kk pour lesquelles le système {kx+y+z=0x+ky+z=0x+y+kz=0\begin{cases} kx+y+z=0 \\ x+ky+z=0 \\ x+y+kz=0 \end{cases} admet des solutions non triviales.
  • d) Pour chacune des valeurs trouvées en c), déterminez l'ensemble des solutions sous forme paramétrique.
  • e) Démontrez que si un système AX=BAX=B admet deux solutions distinctes, alors il en admet une infinité.

Tape tes réponses, la page te dit juste ou faux 0/6

a)
b)
c)
d)
Voir la correction

Réponses

  • a) X=0X=0 toujours solution
  • b) Non triviales si detA=0\det A=0
  • c) k=1k=1 ou k=2k=-2
  • d) k=1k=1 : plan ; k=2k=-2 : droite t(1;1;1)t(1;1;1)
  • e) X1+λ(X2X1)X_1+\lambda(X_2-X_1) : infinité

a) Quel que soit AA, le vecteur nul vérifie A0=0A\cdot 0=0. Le système homogène admet donc toujours au moins la solution X=0X=0, c'est-à-dire x=y=z=0x=y=z=0, appelée solution TRIVIALE. Un système homogène ne peut donc jamais être incompatible : les trois cas généraux se réduisent ici à deux, soit la solution triviale est la seule, soit il y en a une infinité.

Lecture sur l'échelonnement : la colonne des seconds membres étant nulle au départ, elle reste nulle à chaque opération élémentaire, puisque toute combinaison de zéros est nulle. Une ligne du type [0 0 0c][0\ 0\ 0\,|\,c] avec c0c\neq 0 ne peut donc jamais apparaître.

b) Le système homogène admet des solutions non triviales si et seulement si det(A)=0\det(A)=0.

Justification. Si det(A)0\det(A)\neq 0, alors AA est inversible, et on peut multiplier AX=0AX=0 par A1A^{-1} : X=A10=0X=A^{-1}\cdot 0=0. La solution triviale est alors la SEULE. Réciproquement, si det(A)=0\det(A)=0, la matrice n'est pas inversible, l'échelonnement produit moins de pivots que d'inconnues, donc au moins un paramètre libre apparaît ; en donnant à ce paramètre une valeur non nulle, on obtient une solution non triviale.

Lecture équivalente en termes de familles de vecteurs : det(A)=0\det(A)=0 signifie que les colonnes de AA sont linéairement dépendantes, et une relation de dépendance entre les colonnes est exactement une solution non triviale de AX=0AX=0. Les deux chapitres disent la même chose.

c) On calcule le déterminant en développant selon la première ligne : det=k(k21)1(k1)+1(1k)=k3kk+1+1k=k33k+2\det=k(k^{2}-1)-1(k-1)+1(1-k)=k^{3}-k-k+1+1-k=k^{3}-3k+2.

On cherche une racine évidente : pour k=1k=1, 13+2=01-3+2=0. On factorise donc par k1k-1 : k33k+2=(k1)(k2+k2)=(k1)(k1)(k+2)=(k1)2(k+2)k^{3}-3k+2=(k-1)(k^{2}+k-2)=(k-1)(k-1)(k+2)=(k-1)^{2}(k+2).

Le déterminant s'annule pour k=1k=1 (racine double) et k=2k=-2. Ce sont les deux seules valeurs donnant des solutions non triviales.

d) Pour k=1k=1, les trois équations deviennent identiques et se réduisent à x+y+z=0x+y+z=0. Il reste un seul pivot pour trois inconnues, donc 31=23-1=2 paramètres libres. En posant y=sy=s et z=tz=t, on obtient x=stx=-s-t, et l'ensemble des solutions est {(st; s; t) : s,tR}\{(-s-t;\ s;\ t)\ :\ s,t\in\mathbb{R}\}. Géométriquement, c'est un PLAN passant par l'origine. La racine double du déterminant se traduit ici par la perte de deux pivots au lieu d'un.

Pour k=2k=-2, le système s'écrit {2x+y+z=0x2y+z=0x+y2z=0\begin{cases} -2x+y+z=0 \\ x-2y+z=0 \\ x+y-2z=0 \end{cases}. En soustrayant la deuxième équation de la première : 3x+3y=0-3x+3y=0, donc x=yx=y. En soustrayant la troisième de la deuxième : 3y+3z=0-3y+3z=0, donc y=zy=z. Ainsi x=y=zx=y=z, et en posant z=tz=t on obtient {(t; t; t) : tR}\{(t;\ t;\ t)\ :\ t\in\mathbb{R}\}, avec un seul paramètre. Vérification sur la première équation : 2t+t+t=0-2t+t+t=0, conforme. Géométriquement, c'est une DROITE passant par l'origine, de vecteur directeur (1; 1; 1)(1;\ 1;\ 1).

Observation à retenir : dans les deux cas, l'ensemble des solutions contient l'origine. C'est une caractéristique des systèmes homogènes, et cela les distingue des systèmes généraux, dont l'ensemble de solutions est une droite ou un plan qui ne passe pas nécessairement par l'origine.

e) Soient X1X_{1} et X2X_{2} deux solutions distinctes, donc AX1=BAX_{1}=B et AX2=BAX_{2}=B avec X1X2X_{1}\neq X_{2}. Posons, pour un réel λ\lambda quelconque, Xλ=X1+λ(X2X1)X_{\lambda}=X_{1}+\lambda(X_{2}-X_{1}).

Calculons : AXλ=AX1+λ(AX2AX1)=B+λ(BB)=B+λ0=BAX_{\lambda}=AX_{1}+\lambda\big(AX_{2}-AX_{1}\big)=B+\lambda(B-B)=B+\lambda\cdot 0=B. Donc XλX_{\lambda} est solution pour TOUT réel λ\lambda.

De plus, ces solutions sont deux à deux distinctes : si Xλ=XμX_{\lambda}=X_{\mu}, alors (λμ)(X2X1)=0(\lambda-\mu)(X_{2}-X_{1})=0, et comme X2X10X_{2}-X_{1}\neq 0 par hypothèse, on doit avoir λ=μ\lambda=\mu. À chaque réel λ\lambda correspond donc une solution différente : le système en admet une infinité.

Ce résultat démontre au passage qu'il ne peut jamais y avoir exactement deux, ou exactement trois solutions. Le nombre de solutions d'un système linéaire est toujours 00, 11 ou l'infini, et cette démonstration explique pourquoi tout cas intermédiaire est impossible. Géométriquement, la formule X1+λ(X2X1)X_{1}+\lambda(X_{2}-X_{1}) est l'équation paramétrique de la droite passant par les deux solutions : dès que deux solutions existent, toute la droite qui les joint est faite de solutions.

Exercice 7 : Problème : discussion complète d'un système à deux paramètres

On considère, pour des réels aa et bb, le système {x+2y+z=12x+5y+3z=43x+7y+az=b\begin{cases} x+2y+z=1 \\ 2x+5y+3z=4 \\ 3x+7y+az=b \end{cases}.

L'objectif est de déterminer, selon les valeurs de aa et de bb, le nombre de solutions du système, puis de les calculer dans chaque cas. C'est le format d'exercice le plus fréquent en fin de chapitre, et celui où l'on perd le plus de points en oubliant un cas.

  • a) Écrivez la matrice augmentée et échelonnez-la en fonction de aa et de bb. On prendra soin de ne jamais diviser par une expression pouvant s'annuler.
  • b) Pour quelles valeurs de aa et de bb le système admet-il une solution unique ?
  • c) Pour quelles valeurs le système n'admet-il aucune solution ? Interprétez géométriquement.
  • d) Pour quelles valeurs admet-il une infinité de solutions ? Donnez alors l'ensemble des solutions sous forme paramétrique.
  • e) Résolvez complètement le système pour a=5a=5 et b=7b=7.
  • f) Calculez le déterminant de la matrice des coefficients en fonction de aa, et expliquez pourquoi il redonne exactement la condition trouvée en b).

Tape tes réponses, la page te dit juste ou faux 0/8

b)
c)
d)
e)
f)
Voir la correction

Réponses

  • a) Dernière ligne (0;0;a4b5)(0;0;a-4\,|\,b-5)
  • b) a4a\neq 4
  • c) a=4a=4, b5b\neq 5
  • d) a=4a=4, b=5b=5 : (t3; 2t; t)(t-3;\ 2-t;\ t)
  • e) (1; 0; 2)(-1;\ 0;\ 2)
  • f) detA=a4\det A=a-4

a) Matrice augmentée : [1211253437ab]\left[\begin{array}{ccc|c} 1 & 2 & 1 & 1 \\ 2 & 5 & 3 & 4 \\ 3 & 7 & a & b \end{array}\right]. Le pivot en haut à gauche vaut 11 et ne dépend d'aucun paramètre, on peut donc l'utiliser sans précaution.

L2L22L1L_{2}\leftarrow L_{2}-2L_{1} donne (0; 1; 12)(0;\ 1;\ 1\,|\,2). L3L33L1L_{3}\leftarrow L_{3}-3L_{1} donne (0; 1; a3b3)(0;\ 1;\ a-3\,|\,b-3).

Le deuxième pivot vaut 11, lui non plus ne dépend pas des paramètres. L3L3L2L_{3}\leftarrow L_{3}-L_{2} donne (0; 0; a31b32)=(0; 0; a4b5)(0;\ 0;\ a-3-1\,|\,b-3-2)=(0;\ 0;\ a-4\,|\,b-5).

La forme échelonnée est [1211011200a4b5]\left[\begin{array}{ccc|c} 1 & 2 & 1 & 1 \\ 0 & 1 & 1 & 2 \\ 0 & 0 & a-4 & b-5 \end{array}\right]. Toute la discussion se joue maintenant sur la dernière ligne, qui s'écrit (a4)z=b5(a-4)z=b-5. Point de méthode essentiel : on n'a jamais divisé par a4a-4, on a seulement soustrait. Diviser une ligne par une expression paramétrique sans discuter sa nullité est l'erreur qui fait perdre la question entière, car elle suppose implicitement a4a\neq 4 et escamote le cas le plus intéressant.

b) Si a40a-4\neq 0, c'est-à-dire a4a\neq 4, on peut diviser et obtenir z=b5a4z=\dfrac{b-5}{a-4}. Il y a alors trois pivots pour trois inconnues, donc SOLUTION UNIQUE, et ce quelle que soit la valeur de bb.

On remonte pour l'obtenir en général. La deuxième ligne donne y=2zy=2-z. La première donne x=12yz=12(2z)z=14+2zz=z3x=1-2y-z=1-2(2-z)-z=1-4+2z-z=z-3. La solution unique est donc (b5a43; 2b5a4; b5a4)\left(\dfrac{b-5}{a-4}-3;\ 2-\dfrac{b-5}{a-4};\ \dfrac{b-5}{a-4}\right).

c) Si a=4a=4, la dernière ligne devient 0z=b50z=b-5. Lorsque b5b\neq 5, cette équation s'écrit 0=b50=b-5 avec b50b-5\neq 0 : elle est impossible. Le système n'admet AUCUNE solution pour a=4a=4 et b5b\neq 5.

Interprétation géométrique : les trois plans n'ont aucun point commun aux trois. Ici les deux premiers plans se coupent bien selon une droite, mais le troisième plan est parallèle à cette droite sans la contenir. La configuration est celle d'un prisme : les plans se coupent deux à deux selon trois droites parallèles distinctes, sans point commun aux trois.

d) Si a=4a=4 ET b=5b=5, la dernière ligne devient 0=00=0, une équation toujours vraie qui n'apporte aucune information. Le système est compatible avec deux pivots pour trois inconnues, donc 32=13-2=1 paramètre libre : il admet une INFINITÉ de solutions.

On pose z=tz=t. La deuxième ligne donne y+t=2y+t=2, donc y=2ty=2-t. La première donne x+2(2t)+t=1x+2(2-t)+t=1, soit x+42t+t=1x+4-2t+t=1, donc x=t3x=t-3. L'ensemble des solutions est {(t3; 2t; t) : tR}\{(t-3;\ 2-t;\ t)\ :\ t\in\mathbb{R}\}.

Vérification pour t=0t=0, soit (3; 2; 0)(-3;\ 2;\ 0) : équation 1, 3+4+0=1-3+4+0=1, conforme ; équation 2, 6+10+0=4-6+10+0=4, conforme ; équation 3 avec a=4a=4, 9+14+0=5=b-9+14+0=5=b, conforme. Géométriquement, les trois plans se coupent selon une droite de vecteur directeur (1; 1; 1)(1;\ -1;\ 1), lu sur les coefficients de tt.

e) Pour a=5a=5 et b=7b=7, on est dans le cas a4a\neq 4, donc solution unique. La dernière ligne donne (54)z=75(5-4)z=7-5, soit z=2z=2. Puis y=2z=0y=2-z=0, et x=z3=1x=z-3=-1.

Vérification dans le système d'origine : équation 1, 1+0+2=1-1+0+2=1, conforme ; équation 2, 2+0+6=4-2+0+6=4, conforme ; équation 3, 3(1)+7(0)+5(2)=3+10=73(-1)+7(0)+5(2)=-3+10=7, conforme. La solution est (1; 0; 2)(-1;\ 0;\ 2).

f) det(A)=12125337a=1(5a21)2(2a9)+1(1415)=5a214a+181=a4\det(A)=\begin{vmatrix} 1 & 2 & 1 \\ 2 & 5 & 3 \\ 3 & 7 & a \end{vmatrix}=1(5a-21)-2(2a-9)+1(14-15)=5a-21-4a+18-1=a-4.

Le déterminant vaut exactement a4a-4, c'est-à-dire précisément le dernier pivot obtenu par échelonnement. Ce n'est pas une coïncidence : les opérations élémentaires utilisées ici, qui consistent uniquement à ajouter à une ligne un multiple d'une autre, ne modifient PAS le déterminant. Le déterminant de la matrice échelonnée est donc égal à celui de AA, et pour une matrice triangulaire il vaut le produit des coefficients diagonaux, soit 1×1×(a4)=a41\times 1\times(a-4)=a-4.

La condition det(A)0\det(A)\neq 0, c'est-à-dire a4a\neq 4, est donc rigoureusement la même que la condition de solution unique trouvée en b). Les deux approches, échelonnement et déterminant, se rejoignent nécessairement.

Ce que le déterminant ne dit PAS, en revanche, mérite d'être souligné. Il tranche la solution unique, mais lorsqu'il s'annule, il est muet sur la suite : il ne distingue pas le cas b5b\neq 5 (aucune solution) du cas b=5b=5 (une infinité), puisqu'il ne dépend même pas de bb. Seul l'échelonnement de la matrice AUGMENTÉE, qui transporte le second membre, permet de séparer ces deux cas. C'est la raison pour laquelle la règle de Cramer devient inutilisable dès que le déterminant s'annule, et pourquoi une discussion complète passe toujours par Gauss.

Exercice 8 : La structure de l'ensemble des solutions

Lorsqu'un système AX=BAX=B possède une infinité de solutions, celles-ci ne sont pas dispersées au hasard : elles s'organisent autour d'une solution particulière. Cet exercice établit la structure X=Xp+XhX=X_{p}+X_{h}, où XpX_{p} est une solution particulière et XhX_{h} parcourt les solutions du système homogène associé.

  • a) Démontrez que si XpX_{p} est une solution de AX=BAX=B et XhX_{h} une solution de AX=0AX=0, alors Xp+XhX_{p}+X_{h} est solution de AX=BAX=B. Démontrez la réciproque : toute solution de AX=BAX=B s'écrit Xp+XhX_{p}+X_{h} pour un certain XhX_{h} solution du système homogène.
  • b) On considère le système {x+y+z=6x+2y+3z=142x+3y+4z=20\begin{cases} x+y+z=6 \\ x+2y+3z=14 \\ 2x+3y+4z=20 \end{cases}. Vérifiez que Xp=(2; 8; 0)X_{p}=(-2;\ 8;\ 0) en est une solution, puis résolvez le système homogène associé.
  • c) Écrivez la solution générale sous la forme Xp+tXhX_{p}+t\,X_{h}, et confrontez-la au paramétrage (t2; 82t; t)(t-2;\ 8-2t;\ t) obtenu par échelonnement direct.
  • d) Interprétez géométriquement : pourquoi l'ensemble des solutions du système homogène passe-t-il toujours par l'origine, alors que celui du système complet, en général, non ? Que représente XpX_{p} dans cette image ?

Tape tes réponses, la page te dit juste ou faux 0/3

b)
c)
d)
Voir la correction

Réponses

  • a) A(Xp+Xh)=BA(X_p+X_h)=B et réciproque
  • b) Xh=t(1; 2; 1)X_h=t(1;\ -2;\ 1)
  • c) (2;8;0)+t(1;2;1)(-2;8;0)+t(1;-2;1) : même ensemble
  • d) Droite translatée par XpX_p

a) Sens direct. Supposons AXp=BAX_{p}=B et AXh=0AX_{h}=0. Alors, par linéarité du produit matriciel, A(Xp+Xh)=AXp+AXh=B+0=BA(X_{p}+X_{h})=AX_{p}+AX_{h}=B+0=B. Donc Xp+XhX_{p}+X_{h} est bien solution de AX=BAX=B. Réciproque. Soit XX une solution quelconque de AX=BAX=B. Posons Xh=XXpX_{h}=X-X_{p}. Alors AXh=A(XXp)=AXAXp=BB=0AX_{h}=A(X-X_{p})=AX-AX_{p}=B-B=0, donc XhX_{h} est solution du système homogène. En écrivant X=Xp+XhX=X_{p}+X_{h}, on a exprimé XX sous la forme voulue. Les deux sens ensemble montrent que l'ensemble des solutions de AX=BAX=B est EXACTEMENT l'ensemble des Xp+XhX_{p}+X_{h} lorsque XhX_{h} parcourt les solutions homogènes.

b) Vérification de Xp=(2; 8; 0)X_{p}=(-2;\ 8;\ 0) : équation 1, 2+8+0=6-2+8+0=6, conforme ; équation 2, 2+16+0=14-2+16+0=14, conforme ; équation 3, 4+24+0=20-4+24+0=20, conforme. C'est bien une solution. Système homogène associé : {x+y+z=0x+2y+3z=02x+3y+4z=0\begin{cases} x+y+z=0 \\ x+2y+3z=0 \\ 2x+3y+4z=0 \end{cases}. Avec L2L2L1L_{2}\leftarrow L_{2}-L_{1} on obtient y+2z=0y+2z=0, et L3L32L1L_{3}\leftarrow L_{3}-2L_{1} donne y+2z=0y+2z=0, ligne identique. Il reste deux pivots pour trois inconnues, donc un paramètre libre. En posant z=tz=t, on a y=2ty=-2t et x=yz=2tt=tx=-y-z=2t-t=t. L'ensemble des solutions homogènes est {t(1; 2; 1) : tR}\{t(1;\ -2;\ 1)\ :\ t\in\mathbb{R}\}.

c) La solution générale est donc X=Xp+tXh=(2; 8; 0)+t(1; 2; 1)=(2+t; 82t; t)X=X_{p}+t\,X_{h}=(-2;\ 8;\ 0)+t(1;\ -2;\ 1)=(-2+t;\ 8-2t;\ t). Le paramétrage obtenu par échelonnement direct était (t2; 82t; t)(t-2;\ 8-2t;\ t) : c'est rigoureusement le même ensemble, puisque 2+t=t2-2+t=t-2. Les deux approches se rejoignent, mais la forme Xp+tXhX_{p}+t\,X_{h} est plus parlante : elle sépare visiblement la partie constante, une solution particulière, de la partie mobile, qui décrit toutes les directions le long desquelles on peut se déplacer sans quitter l'ensemble des solutions.

d) Le système homogène AX=0AX=0 admet toujours la solution triviale X=0X=0 : son ensemble de solutions contient donc l'origine, et c'est ici une droite vectorielle, de vecteur directeur (1; 2; 1)(1;\ -2;\ 1), passant par OO. L'ensemble des solutions du système complet, lui, est cette même droite TRANSLATÉE par le vecteur XpX_{p} : elle a la même direction (1; 2; 1)(1;\ -2;\ 1), mais elle est décalée et ne passe plus par l'origine, sauf dans le cas particulier où B=0B=0. C'est le sens géométrique de la décomposition : XhX_{h} donne la DIRECTION de l'ensemble des solutions, un objet vectoriel passant par l'origine, et XpX_{p} donne la TRANSLATION qui le place au bon endroit. On dit que l'ensemble des solutions d'un système compatible est un sous-espace AFFINE, c'est-à-dire un sous-espace vectoriel déplacé. Choisir une autre solution particulière ne change pas cet ensemble : cela revient seulement à repérer le même objet à partir d'un autre de ses points.

-4-3-2-1123456-4-3-2-1123456solution particulièresolutions du système HOMOGÈNEl'ensemble des solutions est une DROITE,translatée de la droite homogènesolution générale = particulière + homogène

Exercice 9 : Interpolation : la parabole qui passe par trois points

Faire passer une courbe par des points imposés se traduit par un système linéaire dont les inconnues sont les coefficients de la courbe. C'est le principe de l'interpolation, omniprésent en calcul numérique. On cherche ici la parabole y=ax2+bx+cy=ax^{2}+bx+c.

  • a) Écrivez le système linéaire d'inconnues aa, bb, cc traduisant le passage de la parabole par les points (1; 2)(1;\ 2), (2; 3)(2;\ 3) et (3; 6)(3;\ 6).
  • b) Résolvez ce système et donnez l'équation de la parabole.
  • c) Le déterminant de la matrice des coefficients est un déterminant de Vandermonde. Calculez-le et montrez qu'il est non nul dès que les trois abscisses sont distinctes. Que garantit ce fait quant à l'existence et l'unicité de la parabole ?
  • d) On impose maintenant le passage par (0; 1)(0;\ 1), (1; 3)(1;\ 3) et (2; 5)(2;\ 5), trois points ALIGNÉS. Résolvez et interprétez la valeur trouvée pour aa. Le système reste-t-il de solution unique ?

Tape tes réponses, la page te dit juste ou faux 0/8

b)
c)
d)
Voir la correction

Réponses

  • a) a+b+c=2a+b+c=2, 4a+2b+c=34a+2b+c=3, 9a+3b+c=69a+3b+c=6
  • b) y=x22x+3y=x^2-2x+3
  • c) Déterminant 2-2 : parabole unique
  • d) y=2x+1y=2x+1 : a=0a=0, solution unique

a) Écrire que la parabole y=ax2+bx+cy=ax^{2}+bx+c passe par un point (x0; y0)(x_{0};\ y_{0}), c'est imposer ax02+bx0+c=y0ax_{0}^{2}+bx_{0}+c=y_{0}. Les trois points donnent donc, en substituant chaque abscisse : {a(1)2+b(1)+c=2a(2)2+b(2)+c=3a(3)2+b(3)+c=6\begin{cases} a(1)^{2}+b(1)+c=2 \\ a(2)^{2}+b(2)+c=3 \\ a(3)^{2}+b(3)+c=6 \end{cases}, c'est-à-dire {a+b+c=24a+2b+c=39a+3b+c=6\begin{cases} a+b+c=2 \\ 4a+2b+c=3 \\ 9a+3b+c=6 \end{cases}. Les inconnues sont les coefficients aa, bb, cc, et non les points, qui sont des données.

b) On élimine cc par différences. La deuxième équation moins la première : 3a+b=13a+b=1. La troisième moins la deuxième : 5a+b=35a+b=3. En soustrayant ces deux relations : 2a=22a=2, donc a=1a=1. Alors b=13a=13=2b=1-3a=1-3=-2, et c=2ab=21+2=3c=2-a-b=2-1+2=3. La parabole est y=x22x+3y=x^{2}-2x+3. Vérification sur les trois points : en x=1x=1, 12+3=21-2+3=2 ; en x=2x=2, 44+3=34-4+3=3 ; en x=3x=3, 96+3=69-6+3=6. Les trois passages sont conformes.

c) La matrice des coefficients est (111421931)\begin{pmatrix} 1 & 1 & 1 \\ 4 & 2 & 1 \\ 9 & 3 & 1 \end{pmatrix}, dont les lignes sont bâties sur les puissances des abscisses x1=1x_{1}=1, x2=2x_{2}=2, x3=3x_{3}=3 : c'est une matrice de Vandermonde. Son déterminant se factorise, au signe près, en produit des différences d'abscisses : les colonnes étant rangées de x2x^{2} à 11, dans l'ordre inverse de la matrice de Vandermonde usuelle, det=(x2x1)(x3x1)(x3x2)=(21)(31)(32)=2\det=-(x_{2}-x_{1})(x_{3}-x_{1})(x_{3}-x_{2})=-(2-1)(3-1)(3-2)=-2. Contrôle direct : 1(23)1(49)+1(1218)=21(2-3)-1(4-9)+1(12-18)=-2. Ce produit est non nul dès que les trois abscisses sont DISTINCTES, car chaque facteur est alors non nul. Un déterminant non nul garantit que le système a une solution unique : il existe donc une parabole, et une seule, passant par trois points d'abscisses distinctes. C'est le résultat fondamental de l'interpolation polynomiale : n+1n+1 points d'abscisses distinctes déterminent un unique polynôme de degré au plus nn.

d) Le système devient {c=1a+b+c=34a+2b+c=5\begin{cases} c=1 \\ a+b+c=3 \\ 4a+2b+c=5 \end{cases}. De la première, c=1c=1. La deuxième donne a+b=2a+b=2, la troisième 4a+2b=44a+2b=4, soit 2a+b=22a+b=2. En soustrayant, a=0a=0, puis b=2b=2. La solution est a=0a=0, b=2b=2, c=1c=1, ce qui donne y=2x+1y=2x+1 : une DROITE. La valeur a=0a=0 indique que la parabole dégénère en droite, ce qui était prévisible puisque les trois points sont alignés sur la droite de pente 22. Le système reste de solution UNIQUE : les abscisses 00, 11, 22 sont distinctes, donc le déterminant de Vandermonde (10)(20)(21)=2-(1-0)(2-0)(2-1)=-2 reste non nul, et l'unicité est garantie. Trois points d'abscisses distinctes déterminent toujours un unique polynôme de degré au plus 22 ; ce polynôme est une vraie parabole si et seulement si les points ne sont pas alignés, auquel cas le coefficient dominant aa s'annule.

-112345246810(1 ; 2)(2 ; 3)(3 ; 6)y = x² - 2x + 3trois points, trois équations, trois inconnues :l'interpolation EST un système linéaire

Exercice 10 : Problème : équilibrer une équation chimique par un système homogène

Équilibrer une réaction chimique, c'est trouver les coefficients qui conservent chaque type d'atome. Le problème se traduit par un système linéaire HOMOGÈNE, dont l'ensemble des solutions révèle pourquoi une réaction s'équilibre toujours, et à une constante multiplicative près. On équilibre la combustion du propane : x1C3H8+x2O2x3CO2+x4H2Ox_{1}\,\mathrm{C_{3}H_{8}}+x_{2}\,\mathrm{O_{2}}\rightarrow x_{3}\,\mathrm{CO_{2}}+x_{4}\,\mathrm{H_{2}O}.

  • a) En écrivant la conservation du carbone, de l'hydrogène et de l'oxygène, établissez le système homogène de trois équations vérifié par x1x_{1}, x2x_{2}, x3x_{3}, x4x_{4}.
  • b) Résolvez ce système. Combien de paramètres libres possède-t-il, et pourquoi était-ce prévisible ?
  • c) Déterminez les plus petits coefficients entiers positifs et écrivez l'équation équilibrée. Vérifiez la conservation de chaque atome.
  • d) Expliquez pourquoi le système admet nécessairement des solutions non triviales, et quel sens physique porte le paramètre libre.

Tape tes réponses, la page te dit juste ou faux 0/6

b)
c)
d)
Voir la correction

Réponses

  • a) 3x1=x33x_1=x_3, 8x1=2x48x_1=2x_4, 2x2=2x3+x42x_2=2x_3+x_4
  • b) t(1; 5; 3; 4)t(1;\ 5;\ 3;\ 4) : un paramètre
  • c) C3H8+5O23CO2+4H2OC_3H_8+5\,O_2\to 3\,CO_2+4\,H_2O
  • d) Plus d'inconnues que d'équations ; échelle

a) On écrit la conservation atome par atome. Carbone : le propane C3H8\mathrm{C_{3}H_{8}} apporte 3x13x_{1} atomes de carbone, le CO2\mathrm{CO_{2}} en emporte x3x_{3}, d'où 3x1=x33x_{1}=x_{3}. Hydrogène : C3H8\mathrm{C_{3}H_{8}} apporte 8x18x_{1} atomes, l'eau H2O\mathrm{H_{2}O} en emporte 2x42x_{4}, d'où 8x1=2x48x_{1}=2x_{4}. Oxygène : le O2\mathrm{O_{2}} apporte 2x22x_{2} atomes, répartis entre 2x32x_{3} dans le CO2\mathrm{CO_{2}} et x4x_{4} dans l'eau, d'où 2x2=2x3+x42x_{2}=2x_{3}+x_{4}. Le système homogène, tous termes ramenés à gauche, est {3x1x3=08x12x4=02x22x3x4=0\begin{cases} 3x_{1}-x_{3}=0 \\ 8x_{1}-2x_{4}=0 \\ 2x_{2}-2x_{3}-x_{4}=0 \end{cases}.

b) On exprime tout en fonction de x1x_{1}. La première équation donne x3=3x1x_{3}=3x_{1}. La deuxième donne x4=4x1x_{4}=4x_{1}. La troisième donne 2x2=2x3+x4=6x1+4x1=10x12x_{2}=2x_{3}+x_{4}=6x_{1}+4x_{1}=10x_{1}, donc x2=5x1x_{2}=5x_{1}. En posant x1=tx_{1}=t, on obtient la famille de solutions (x1; x2; x3; x4)=(t; 5t; 3t; 4t)=t(1; 5; 3; 4)(x_{1};\ x_{2};\ x_{3};\ x_{4})=(t;\ 5t;\ 3t;\ 4t)=t(1;\ 5;\ 3;\ 4). Il y a UN seul paramètre libre. C'était prévisible : le système compte quatre inconnues pour seulement trois équations, donc au moins 43=14-3=1 inconnue reste libre. Un système homogène ayant plus d'inconnues que d'équations a toujours des solutions non triviales.

c) La famille t(1; 5; 3; 4)t(1;\ 5;\ 3;\ 4) donne des coefficients entiers et positifs pour tt entier positif ; les plus petits correspondent à t=1t=1, soit (x1; x2; x3; x4)=(1; 5; 3; 4)(x_{1};\ x_{2};\ x_{3};\ x_{4})=(1;\ 5;\ 3;\ 4). L'équation équilibrée est C3H8+5O23CO2+4H2O\mathrm{C_{3}H_{8}}+5\,\mathrm{O_{2}}\rightarrow 3\,\mathrm{CO_{2}}+4\,\mathrm{H_{2}O}. Vérification atome par atome : carbone, 33 à gauche (3×13\times 1) et 33 à droite (3×13\times 1) ; hydrogène, 88 à gauche (8×18\times 1) et 88 à droite (4×24\times 2) ; oxygène, 1010 à gauche (5×25\times 2) et 1010 à droite (3×2+4×13\times 2+4\times 1). Chaque type d'atome est conservé : l'équation est équilibrée.

d) Le système est homogène, donc il admet toujours la solution triviale x1=x2=x3=x4=0x_{1}=x_{2}=x_{3}=x_{4}=0, qui n'a évidemment aucun sens chimique. Mais il admet AUSSI des solutions non triviales, et ce nécessairement, pour une raison structurelle : avec quatre inconnues et trois équations seulement, la forme échelonnée ne peut avoir au plus que trois pivots, laissant au moins une inconnue libre ; le rang est strictement inférieur au nombre d'inconnues, ce qui garantit une infinité de solutions. Autrement dit, une réaction met en jeu plus d'espèces que d'éléments distincts, et cette inégalité assure l'existence d'un équilibrage non trivial. Le sens physique du paramètre libre est celui de l'ÉCHELLE de la réaction : multiplier tous les coefficients par un même facteur décrit la même réaction menée en plus grande ou plus petite quantité. La chimie fixe les proportions, pas les quantités absolues, et c'est exactement ce qu'exprime la droite de solutions t(1; 5; 3; 4)t(1;\ 5;\ 3;\ 4). On choisit par convention les plus petits entiers, mais tous les points de cette droite décrivent la même transformation.

Partie C : les classiques (/50)

Exercice 11 : Pivot nul et échange de lignes

On considère le système {y+2z=8x+y+z=62x+yz=1\begin{cases} y+2z=8 \\ x+y+z=6 \\ 2x+y-z=1 \end{cases}.

  • a) Écrivez la matrice augmentée. Pourquoi ne peut-on pas commencer l'élimination avec le coefficient en haut à gauche ?
  • b) Échangez les deux premières lignes, puis échelonnez la matrice.
  • c) Résolvez le système et vérifiez la solution.
  • d) Calculez le déterminant de la matrice des coefficients avant et après l'échange. Que change l'échange, et que ne change-t-il pas ?
  • e) Si toute la première colonne sous la diagonale avait été nulle, avec un zéro en haut, qu'aurait-on dû conclure sur l'inconnue xx ?

Tape tes réponses, la page te dit juste ou faux 0/9

a)
b)
c)
d)
e)
Voir la correction

Réponses

  • a) Pivot nul
  • b) L1L2L_1\leftrightarrow L_2 ; dernière ligne (0;0;13)(0;0;-1\,|\,-3)
  • c) (1; 2; 3)(1;\ 2;\ 3)
  • d) 11 puis 1-1 : signe changé
  • e) xx paramètre libre

a) [012811162111]\left[\begin{array}{ccc|c} 0 & 1 & 2 & 8 \\ 1 & 1 & 1 & 6 \\ 2 & 1 & -1 & 1 \end{array}\right]. Le coefficient en haut à gauche est nul : aucun multiple de la première ligne ne peut annuler le 11 ou le 22 placés en dessous. Un pivot doit être NON NUL.

b) L1L2L_{1}\leftrightarrow L_{2} donne [111601282111]\left[\begin{array}{ccc|c} 1 & 1 & 1 & 6 \\ 0 & 1 & 2 & 8 \\ 2 & 1 & -1 & 1 \end{array}\right]. Puis L3L32L1L_{3}\leftarrow L_{3}-2L_{1} donne (0; 1; 311)(0;\ -1;\ -3\,|\,-11), et L3L3+L2L_{3}\leftarrow L_{3}+L_{2} donne (0; 0; 13)(0;\ 0;\ -1\,|\,-3).

c) z=3-z=-3, donc z=3z=3 ; y+6=8y+6=8, donc y=2y=2 ; x+2+3=6x+2+3=6, donc x=1x=1. Vérification : 2+6=82+6=8, 1+2+3=61+2+3=6 et 2+23=12+2-3=1.

d) Avant : 0(11)1(12)+2(12)=32=10(-1-1)-1(-1-2)+2(1-2)=3-2=1. Après l'échange, le déterminant vaut 1-1 : échanger deux lignes change son SIGNE. L'ensemble des solutions, lui, ne change pas, puisqu'on a seulement réécrit les mêmes équations dans un autre ordre.

e) Une colonne sans aucun pivot possible signifie que l'inconnue xx n'apparaît dans aucune équation : si le système est compatible, xx est un paramètre libre, et il n'y a pas de solution unique.

Le fil : un pivot nul se contourne par un échange de lignes, qui ne change pas les solutions mais change le signe du déterminant.

Exercice 12 : Un système de quatre équations par Gauss

On considère le système {x+y+z+w=10xy+zw=22x+yz+w=5x+2y+3zw=10\begin{cases} x+y+z+w=10 \\ x-y+z-w=-2 \\ 2x+y-z+w=5 \\ x+2y+3z-w=10 \end{cases}.

  • a) Éliminez xx des trois dernières équations en utilisant la première. Écrivez les trois lignes obtenues.
  • b) Poursuivez l'échelonnement jusqu'à une forme triangulaire.
  • c) Résolvez le système par remontée.
  • d) Vérifiez la solution dans les quatre équations d'origine.
  • e) Quel est le rang de la matrice des coefficients ? Combien de solutions le système possède-t-il ?

Tape tes réponses, la page te dit juste ou faux 0/9

a)
b)
c)
e)
Voir la correction

Réponses

  • a) (0;2;0;212)(0;-2;0;-2\,|\,-12), (0;1;3;115)(0;-1;-3;-1\,|\,-15), (0;1;2;20)(0;1;2;-2\,|\,0)
  • b) Pivots 11, 11, 3-3, 9-9
  • c) (1; 2; 3; 4)(1;\ 2;\ 3;\ 4)
  • d) Quatre équations vérifiées
  • e) Rang 4 : solution unique

a) L2L2L1L_{2}\leftarrow L_{2}-L_{1} : (0; 2; 0; 212)(0;\ -2;\ 0;\ -2\,|\,-12). L3L32L1L_{3}\leftarrow L_{3}-2L_{1} : (0; 1; 3; 115)(0;\ -1;\ -3;\ -1\,|\,-15). L4L4L1L_{4}\leftarrow L_{4}-L_{1} : (0; 1; 2; 20)(0;\ 1;\ 2;\ -2\,|\,0).

b) On simplifie L2L_{2} par 2-2 : (0; 1; 0; 16)(0;\ 1;\ 0;\ 1\,|\,6). Puis L3L3+L2L_{3}\leftarrow L_{3}+L_{2} : (0; 0; 3; 09)(0;\ 0;\ -3;\ 0\,|\,-9) et L4L4L2L_{4}\leftarrow L_{4}-L_{2} : (0; 0; 2; 36)(0;\ 0;\ 2;\ -3\,|\,-6). Enfin L43L4+2L3L_{4}\leftarrow 3L_{4}+2L_{3} : (0; 0; 0; 936)(0;\ 0;\ 0;\ -9\,|\,-36).

c) 9w=36-9w=-36, donc w=4w=4. 3z=9-3z=-9, donc z=3z=3. y+w=6y+w=6, donc y=2y=2. x+y+z+w=10x+y+z+w=10, donc x=1x=1. La solution est (1; 2; 3; 4)(1;\ 2;\ 3;\ 4).

d) 1+2+3+4=101+2+3+4=10 ; 12+34=21-2+3-4=-2 ; 2+23+4=52+2-3+4=5 ; 1+4+94=101+4+9-4=10. Les quatre équations sont vérifiées.

e) La forme échelonnée a quatre pivots, 11, 11, 3-3 et 9-9 : le rang vaut 44, égal au nombre d'inconnues, et la solution est unique. La méthode est la même qu'en dimension 3, avec une étape de plus ; multiplier une ligne par 33 pour éviter les fractions est permis, à condition de ne jamais multiplier par zéro.

Le fil : quelle que soit la taille, Gauss élimine une inconnue par colonne, et le nombre de pivots décide du nombre de solutions.

Exercice 13 : La décomposition LU

On considère A=(211433879)A=\begin{pmatrix} 2 & 1 & 1 \\ 4 & 3 & 3 \\ 8 & 7 & 9 \end{pmatrix}. Échelonner AA sans échange de lignes revient à l'écrire A=LUA=LU, avec UU triangulaire supérieure, la forme échelonnée, et LL triangulaire inférieure, de diagonale formée de 11, qui garde les multiplicateurs utilisés.

  • a) Échelonnez AA en notant chaque multiplicateur, et donnez UU.
  • b) Écrivez LL en plaçant chaque multiplicateur à la position du coefficient qu'il a annulé, puis vérifiez que LU=ALU=A sur la dernière ligne.
  • c) Déduisez-en detA\det A.
  • d) Résolvez AX=BAX=B avec B=(41024)B=\begin{pmatrix} 4 \\ 10 \\ 24 \end{pmatrix} en deux temps : LY=BLY=B par descente, puis UX=YUX=Y par remontée.
  • e) Pourquoi cette décomposition est-elle avantageuse quand on doit résoudre de nombreux systèmes avec la même matrice AA ?

Tape tes réponses, la page te dit juste ou faux 0/10

a)
b)
c)
d)
e)
Voir la correction

Réponses

  • a) U=(211011002)U=\begin{pmatrix}2&1&1\\0&1&1\\0&0&2\end{pmatrix}
  • b) L=(100210431)L=\begin{pmatrix}1&0&0\\2&1&0\\4&3&1\end{pmatrix}
  • c) detA=4\det A=4
  • d) Y=(4;2;2)Y=(4;2;2), X=(1; 1; 1)X=(1;\ 1;\ 1)
  • e) Échelonnement fait une fois

a) L2L22L1L_{2}\leftarrow L_{2}-2L_{1} donne (0; 1; 1)(0;\ 1;\ 1), multiplicateur 22. L3L34L1L_{3}\leftarrow L_{3}-4L_{1} donne (0; 3; 5)(0;\ 3;\ 5), multiplicateur 44. L3L33L2L_{3}\leftarrow L_{3}-3L_{2} donne (0; 0; 2)(0;\ 0;\ 2), multiplicateur 33. Donc U=(211011002)U=\begin{pmatrix} 2 & 1 & 1 \\ 0 & 1 & 1 \\ 0 & 0 & 2 \end{pmatrix}.

b) L=(100210431)L=\begin{pmatrix} 1 & 0 & 0 \\ 2 & 1 & 0 \\ 4 & 3 & 1 \end{pmatrix}. Dernière ligne de LULU : (42; 41+31; 41+31+12)=(8; 7; 9)(4\cdot 2;\ 4\cdot 1+3\cdot 1;\ 4\cdot 1+3\cdot 1+1\cdot 2)=(8;\ 7;\ 9), c'est bien la dernière ligne de AA.

c) detA=detL×detU=1×(2×1×2)=4\det A=\det L\times\det U=1\times(2\times 1\times 2)=4 : les déterminants de matrices triangulaires se lisent sur leur diagonale.

d) LY=BLY=B : y1=4y_{1}=4, 2y1+y2=102y_{1}+y_{2}=10 donc y2=2y_{2}=2, 4y1+3y2+y3=244y_{1}+3y_{2}+y_{3}=24 donc y3=2y_{3}=2. UX=YUX=Y : 2z=22z=2 donc z=1z=1, y+z=2y+z=2 donc y=1y=1, 2x+y+z=42x+y+z=4 donc x=1x=1. La solution est (1; 1; 1)(1;\ 1;\ 1), et en effet A(111)=(41024)A\begin{pmatrix}1\\1\\1\end{pmatrix}=\begin{pmatrix}4\\10\\24\end{pmatrix}.

e) L'échelonnement, la partie coûteuse, n'est fait qu'une fois. Chaque nouveau second membre ne demande plus qu'une descente et une remontée, deux calculs triangulaires très rapides : c'est ainsi que les logiciels résolvent les grands systèmes, plutôt qu'en calculant A1A^{-1}.

Le fil : la matrice LL garde la mémoire de l'échelonnement ; résoudre AX=BAX=B devient deux systèmes triangulaires.

Exercice 14 : Problème : la température d'équilibre d'une plaque

Une plaque carrée est maintenue à 100 C100\ ^{\circ}\mathrm{C} sur son bord supérieur, à 0 C0\ ^{\circ}\mathrm{C} sur son bord inférieur et à 50 C50\ ^{\circ}\mathrm{C} sur ses bords gauche et droit. On étudie quatre points intérieurs disposés en carré : T1T_{1} en haut à gauche, T2T_{2} en haut à droite, T3T_{3} en bas à gauche et T4T_{4} en bas à droite.

À l'équilibre, la température de chaque point intérieur est la moyenne des températures de ses quatre voisins : les deux points du bord les plus proches et les deux points intérieurs adjacents.

  • a) Écrivez l'équation vérifiée par T1T_{1}, puis les trois autres.
  • b) Mettez le système sous la forme AT=BAT=B et échelonnez-le.
  • c) Déterminez les quatre températures.
  • d) Retrouvez le résultat en utilisant la symétrie gauche-droite de la plaque.
  • e) Calculez la moyenne des quatre températures. Pourquoi était-elle prévisible ?

Tape tes réponses, la page te dit juste ou faux 0/5

a)
c)
d)
e)
Voir la correction

Réponses

  • a) 4T1=150+T2+T34T_1=150+T_2+T_3, etc.
  • b) Quatre pivots : solution unique
  • c) T1=T2=62,5T_1=T_2=62{,}5 °C, T3=T4=37,5T_3=T_4=37{,}5 °C
  • d) 3ab=1503a-b=150, a+3b=50-a+3b=50
  • e) 50 °C

a) T1T_{1} a pour voisins le bord supérieur, 100100, le bord gauche, 5050, puis T2T_{2} et T3T_{3} : 4T1=150+T2+T34T_{1}=150+T_{2}+T_{3}. De même 4T2=150+T1+T44T_{2}=150+T_{1}+T_{4}, 4T3=50+T1+T44T_{3}=50+T_{1}+T_{4} et 4T4=50+T2+T34T_{4}=50+T_{2}+T_{3}.

b) [41101501401150104150011450]\left[\begin{array}{cccc|c} 4 & -1 & -1 & 0 & 150 \\ -1 & 4 & 0 & -1 & 150 \\ -1 & 0 & 4 & -1 & 50 \\ 0 & -1 & -1 & 4 & 50 \end{array}\right]. Sans fraction : L24L2+L1L_{2}\leftarrow 4L_{2}+L_{1} donne (0; 15; 1; 4750)(0;\ 15;\ -1;\ -4\,|\,750) et L34L3+L1L_{3}\leftarrow 4L_{3}+L_{1} donne (0; 1; 15; 4350)(0;\ -1;\ 15;\ -4\,|\,350) ; puis L315L3+L2L_{3}\leftarrow 15L_{3}+L_{2} donne (0; 0; 224; 646000)(0;\ 0;\ 224;\ -64\,|\,6\,000) et L415L4+L2L_{4}\leftarrow 15L_{4}+L_{2} donne (0; 0; 16; 561500)(0;\ 0;\ -16;\ 56\,|\,1\,500) ; enfin L414L4+L3L_{4}\leftarrow 14L_{4}+L_{3} donne (0; 0; 0; 72027000)(0;\ 0;\ 0;\ 720\,|\,27\,000). Quatre pivots non nuls : la solution est unique.

c) Remontée : 720T4=27000720T_{4}=27\,000, donc T4=37,5T_{4}=37{,}5 ; 224T3=6000+64×37,5=8400224T_{3}=6\,000+64\times 37{,}5=8\,400, donc T3=37,5T_{3}=37{,}5 ; 15T2=750+37,5+4×37,5=937,515T_{2}=750+37{,}5+4\times 37{,}5=937{,}5, donc T2=62,5T_{2}=62{,}5 ; 4T1=150+62,5+37,54T_{1}=150+62{,}5+37{,}5, donc T1=62,5T_{1}=62{,}5. Contrôle : 4×62,5=250=150+62,5+37,54\times 62{,}5=250=150+62{,}5+37{,}5 et 4×37,5=150=50+62,5+37,54\times 37{,}5=150=50+62{,}5+37{,}5.

d) La plaque et ses conditions de bord sont symétriques par rapport à l'axe vertical, et la solution est unique : elle est donc symétrique, T1=T2=aT_{1}=T_{2}=a et T3=T4=bT_{3}=T_{4}=b. Le système se réduit à 3ab=1503a-b=150 et a+3b=50-a+3b=50. La première donne b=3a150b=3a-150, puis a+9a450=50-a+9a-450=50, soit a=62,5a=62{,}5 et b=37,5b=37{,}5.

e) 62,5+62,5+37,5+37,54=50 C\frac{62{,}5+62{,}5+37{,}5+37{,}5}{4}=50\ ^{\circ}\mathrm{C}, la moyenne des quatre bords, 100+0+50+504\frac{100+0+50+50}{4}. Ce n'est pas un hasard : en additionnant les quatre équations, 4Ti=400+2Ti4\sum T_{i}=400+2\sum T_{i}, donc Ti=200\sum T_{i}=200. La méthode, répétée sur des grilles de milliers de points, est celle qu'emploient les logiciels de calcul thermique.

Le fil : une loi locale, la moyenne des voisins, produit un système linéaire ; la symétrie et l'unicité permettent de réduire le nombre d'inconnues.

Exercice 15 : Problème : un système mal conditionné

On considère le système {x+y=2x+1,001y=2,001\begin{cases} x+y=2 \\ x+1{,}001\,y=2{,}001 \end{cases}, qui modélise deux mesures presque redondantes.

  • a) Résolvez le système.
  • b) Calculez le déterminant de la matrice des coefficients.
  • c) Une erreur de mesure change le second membre 2,0012{,}001 en 2,0022{,}002. Résolvez le nouveau système.
  • d) Le second membre a varié de combien de pour cent, au centième ? Comparez à la variation de la solution.
  • e) Interprétez géométriquement, et dites ce qui arrive si le coefficient 1,0011{,}001 est arrondi à 1,001{,}00.

Tape tes réponses, la page te dit juste ou faux 0/7

a)
b)
c)
d)
e)
Voir la correction

Réponses

  • a) (1; 1)(1;\ 1)
  • b) 0,0010{,}001
  • c) (0; 2)(0;\ 2)
  • d) 0,050{,}05 % contre une solution bouleversée
  • e) Droites presque parallèles ; arrondi : aucune solution

a) En soustrayant la première équation de la seconde : 0,001y=0,0010{,}001\,y=0{,}001, donc y=1y=1 et x=1x=1.

b) det=1×1,0011×1=0,001\det=1\times 1{,}001-1\times 1=0{,}001 : non nul, la solution est unique, mais le déterminant est minuscule.

c) 0,001y=0,0020{,}001\,y=0{,}002, donc y=2y=2 et x=0x=0. La solution passe de (1; 1)(1;\ 1) à (0; 2)(0;\ 2).

d) 0,0012,0010,0005\frac{0{,}001}{2{,}001}\approx 0{,}0005, soit environ 0,05 %0{,}05\ \%. La solution, elle, a changé du tout au tout : xx passe de 11 à 00, yy double. Une erreur minuscule sur les données produit une erreur énorme sur le résultat : le système est MAL CONDITIONNÉ.

e) Les deux équations représentent deux droites presque parallèles : leurs pentes, 1-1 et 11,001-\frac{1}{1{,}001}, diffèrent à peine. Déplacer très légèrement l'une d'elles fait glisser leur point d'intersection très loin le long des droites. Avec 1,001{,}00, les droites deviennent parallèles : x+y=2x+y=2 et x+y=2,001x+y=2{,}001 n'ont plus aucune solution, et le système change de nature pour un arrondi.

Le fil : un déterminant proche de zéro, relativement aux coefficients, annonce une solution très sensible aux erreurs de mesure.

Chapitre précédent Matrices, systèmes et déterminants Chapitre suivant Déterminants, cofacteurs et règle de Cramer

Ce chapitre resservira dans

Les chapitres qui le réclament en amont, plus tard dans l'année ou dans les années suivantes.

Voir aussi

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

Contactez-moi pour une première séance. On travaille le 201-NYC et le complément québécois de Première au niveau réel des évaluations, de l'échelonnement jusqu'à la discussion complète d'un système à paramètres.

Site par Studio Squalli