L’ALGORITHME DE GAUSS-JORDAN
1. La forme échelonnée (réduite) d’une matrice augmentée
Considérer le système linéaire suivant.
a11 x1 + · · · + a1n xn =b1
..
.
am1 x1 + · · · + amn xn =bm
Sa matrice augmentée est la matrice de taille m × (n + 1)
a11 ··· a1n b1
.. .. .. ..
[A.~b] := . . .
am1 ··· amn bm .
.
La matrices augmentée [A..~b] est en forme échelonnée si les trois conditions suiv-
antes sont satisfaites.
.
(1) Si ~`i [A..~b] 6= ~0, alors son premier coefficient non-nul depuis la gauche est
.
un 1, que l’on appelle le pivot de ~`i [A..~b] .
. .
(2) Si ~`i [A..~b] = ~0, alors ~`j [A..~b] = ~0 pour tout j > i. Autrement dit, toutes
les lignes de zéros se trouvent ensemble en bas de la matrice.
. . . .
(3) Si ~`i [A..~b] 6= ~0 et ~`i+1 [A..~b] 6= ~0, avec pivots [A..~b])ij = 1 et [A..~b])i+1k =
1, alors j < k. Autrement dit le pivot de la ième -ligne se trouve à gauche
de celui de la (i + 1)ème -ligne.
Elle est en forme échelonnée réduite si, de plus, la condition suivante est vérifiée.
. . .
(4) Si ~cj [A..~b] contient un pivot [A..~b])ij = 1, alors [A..~b])kj = 0 pour tout
k 6= i. Autrement dit, un pivot est le seul coefficient non-nul dans sa colonne.
L’intérêt des matrices augmentées en forme échelonnée réside dans la facilité avec
laquelle on trouve les solutions du système linéaire y associé.
Exemples.
(1) Une matrice en forme échelonnée et son système linéaire associé:
1 2 3 4 5
0 0 1 2 3
0 0 0 0 0
1
2 L’ALGORITHME DE GAUSS-JORDAN
et
x1 + 2x2 + 3x3 + 4x4 = 5
x3 + 2x4 = 3
0=0
On voit facilement que les solutions de ce système sont les vecteurs
(−4 − 2s + 2t, s, 3 − 2t, t)
pour tous s, t ∈ F.
(2) Une matrice en forme échelonnée réduite et son système linéaire associé:
1 0 2 0 3
0 1 −1 0 −2
0 0 0 1 −7
et
x1
+ 2x3 =3
x2 − x3 = −2
x4 = −7
On voit encore plus facilement que les solutions de ce système sont les
vecteurs
(3 − 2s, −2 + s, s, −7)
pour tout s ∈ F.
2. L’algorithme de Gauss-Jordan
Nous résumons ici un algorithme qui permet de transformer toute matrice en
une matrice en forme échelonnée (réduite). Cet algorithme utilise les opérations
suivantes.
Définition. Soit A une matrice de taille m × n, à coefficients dans F. Il y a trois
types d’opérations élémentaires sur les lignes de la matrice A.
(1) Opération Ii,j : Permuter ~`i (A) et ~`j (A), pour obtenir une nouvelle matrice
A0 de même taille et telle que
~` (A) : k = j
i
~`k (A0 ) = ~`j (A) : k = i
~
`k (A) : sinon.
(2) Opération IIi,λ : Multiplier ~`i (A) par λ ∈ F, pour obtenir une nouvelle
matrice A00 de même taille et telle que
λ · ~`i (A) : k = i
(
~`k (A ) =
00
~`k (A) : sinon.
(3) Opération IIIi,j,λ : Remplacer ~`i (A) par ~`i (A) + λ · ~`j (A) pour obtenir une
nouvelle matrice A000 de même taille et telle que
~`i (A) + λ · ~`j (A) : k = i
(
~`k (A000 ) =
~`k (A) : sinon.
L’ALGORITHME DE GAUSS-JORDAN 3
L’algorithme. Soit A une matrice m × n.
(1) Soit j0 = min{j | ~cj (A) 6= ~0}. Autrement dit, ~cj0 (A) est la première colonne
non-nulle de A depuis la gauche. S’il n’y a pas de colonne non-nulle, alors
la matrice A est la matrice zéro, qui est en forme échelonné réduite, et
l’algorithme s’arrête.
(2) Si (A)1j0 = 0, appliquer à A l’opération I1,i0 , où (A)i0 j0 6= 0, pour obtenir
une nouvelle matrice A0 telle que (A0 )1j0 6= 0. (Il existe forcément un tel i0 ,
puisque ~cj0 (A) 6= ~0.)
(3) Appliquer à A0 l’opération IIi,λ , où λ = (A0 )−1
1j0 , pour obtenir une nouvelle
00 00
matrice A telle que (A )1j0 = 1.
(4) Pour tout i > 1, appliquer à A00 l’opération IIIi,1,λi , où λi = −(A00 )ij0 , pour
obtenir une nouvelle matrice A000 telle que (A000 )ij0 = 0 quelque soit i > 1.
(5) Soit B la matrice de taille (m − 1) × (n − j0 ) obtenue de A000 en enlevant les
colonnes ~c1 (A000 ), ..., ~cj0 (A000 ) et la ligne ~`1 (A000 ). Appliquer les étapes (1)-(4)
de l’algorithme à B et ensuite l’étape (5) à la matrice B 000 qui en résulte.
Après au plus m circuits à travers l’algorithme, il s’arrête, et la matrice  obtenue
à la fin est en forme échelonnée. Si l’on veut une forme échelonnée réduite, il faut
suivre les étapes ci-dessous.
(6) Soit
j1 = max{j | ~cj (Â) contient un pivot et au moins un autre coefficient non-nul}.
Si toute colonne contenant un pivot n’a aucun autre coefficient non-nul,
l’algorithme s’arrête. Soit (Â)i1 j1 = 1 le pivot.
Pour tout i < i1 , appliquer à Â l’opération IIIi,i1 ,µi , où µi = −(Â)ij1 ,
pour obtenir une nouvelle matrice Â0 telle que (Â0 )ij1 = 0 pour tout i 6= i1 .
(7) Appliquer l’étape (6) à la matrice Â0 .
Cette deuxième partie de l’algorithme s’arrête après au plus n étapes. La matrice
Ǎ obtenue à la fin est en forme échelonnée réduite.
Une illustration détaillée de l’application de cet algorithme sera fournie au cours.
