Cours d’Analyse numérique 1, Filière SMA, [Link] Ibn Tofail.

Université Ibn Tofail

Faculté des Sciences. Kénitra
Filière : S.M.A

Extrait du polycopié
Module : Analyse Numérique 1
7 mars 2024

(Méthodes directes de résolution des systèmes linéaires)

M
Une méthode directe de résolution d’un système linéaire est une méthode qui

RI
permet de résoudre ce système en un nombre fini d’opérations. Parmis les méthodes
on a la méthode de Gauss, la méthode de Gauss-Jordan, la factorisation LU , la fac-
torisation de Cholesky,· · · .
UK
Considérons le système linéaire AX = B
     
a1,1 a1,2 · · · a1,n x1 b1
BO
 a2,1 a2,2 · · · a2,n   x2   b2 
où , A =  .. ..  , X =  ..  et B =  .. .
     
.. . .
 . . . .  . .
an,1 an,2 · · · an,n xn bn

1- Méthode de Cramer.
f. L

On sait que lorsque la matrice A est inversible alors la solution existe et elle est
unique, cette solution est donnée par la méthode de Cramer (X = A−1 B),
Pro

a1,1 · · · a1,i−1 b1 a1,i+1 · · · a1,n

a2,1 · · · a2,i−1 b2 a2,i+1 · · · a2,n 1
ou encore : xi = .. .. ... .. × detA .
. . .
an,1 · · · an,i−1 bn an,i+1 · · · an,n
Mais cette méthode de Cramer necessite un trés grand nombre d’opérations
(+, −, ×, ÷), ce qui la rend du point de vue pratique la méthode la moins recom-
mendée voire même impossible dés que n est grand (n ≥ 20 par exemple).

A.U : 2023-2024 1 Prof. L BOUKRIM

Cours d’Analyse numérique 1, Filière SMA, [Link] Ibn Tofail.

Exercice1.
Evaluer le nombre d’opérations nécessaires pour la résolution du système linéaire
avec la méthode de Cramer (exemple n = 15, n = 25).

2- Cas de système triangulaire.

Supposons que le système cette fois-ci est triangulaire, exemple triangulaire supérieur :

a1,1 a1,2 · · · a1,n

 
.
 0 a2,2 . . a2,n 
 
A= . .. .. ..  .
 .. . . . 

M
0 0 · · · an,n

RI
Dans ce cas, on détermine la solution X en commençant par la dérnière équation
du système AX = B qui donne directement l’inconnue xn , puis en remontant, on
obtient les autres inconnues xn−1 , xn−2 , · · · , x1 . On a alors :





xn = abn,n
n UK
n
X
1
 xi = (b − ai,j xj ) (i = 1, · · · , n − 1).

ai,i i

BO

j=i+1

Le cas triangulaire inférieur donne :

x = b1

 1 a1,1

f. L



i−1
X
1
 xi = ai,i (bi − ai,j xj ) (i = 2, · · · , n).



j=1
Pro

Exercice2.
Evaluer le nombre d’opérations nécessaires pour la résolution d’un système triangu-
laire (exemple n = 15, n = 25).

Vu que le cas de système linéaire non triangulaire nécessite un trés grand nombre
d’opérations, ce qui rend la méthode de Cramer inutilisable en pratique voire im-
possible dés que n est grand (n ≥ 20 par exemple), il faudrait donc trouver d’autres
méthodes de résolution.

A.U : 2023-2024 2 Prof. L BOUKRIM

Cours d’Analyse numérique 1, Filière SMA, [Link] Ibn Tofail.

3- Méthode de Gauss.
Cette méthode permet de transformer le système régulier AX = B quelconque en
un système triangulaire équivalent A′ X = B ′ où A′ est une matrice triangulaire
supérieure.

Remarques. La solution du système linéaire reste la même si :

- On divise une ligne par un nombre non nul.
- On retranche une ligne à une autre ligne.
- On permute deux lignes.

M
Description de la méthode de Gauss.

RI

 a1,1 x1 + · · · + a1,n xn = b1 (E1 )

AX = B ⇔ .. (S)
 .

UK
 a x + · · · + a x = b (E ).
n,1 1 n,n n n n

Etape1 : On suppose a1,1 ̸= 0 (a1,1 =pivot), on transforme (S) en un autre système

équivalent dont les ai,1 = 0 , ∀i = 2, · · · , n.
La méthode est : L’équation (Ei ) est remplacée par l’équation
BO
(2) ai,1
(Ei ) = (Ei ) − (Ei ) (∀i = 2, · · · , n).
a1,1
On obtient un autre système équivalent au premier qui est de la forme :
f. L

 (2)

a1,1 x1 + a1,2 x2 + · · · + a1,n xn = b1 (E1 )
(2) (2) (2) (2)
a2,2 x2 + · · · +a2,n xn = b2

 0+ (E2 )
. (S (2) )
Pro

.. .. ..


 . . .
 (2) (2) (2) (2)
0+ an,2 x2 + · · · +an,n xn = bn (En ).

Où l’indice (2) dans le système (S (2) ) signifie l’étape 2 de la méthode de Gauss.
(2)
Si a2,2 ̸= 0, on transforme le système (S (2) ) en un autre système équivalent (S (3) ) où
(2)
l’équation (Ei ) est remplacée par par l’équation
(2)
(3) (2) ai,2 (2)
(Ei ) = (Ei ) − (2)
(Ei ) (∀i = 3, · · · , n).
a2,2

A.U : 2023-2024 3 Prof. L BOUKRIM

Cours d’Analyse numérique 1, Filière SMA, [Link] Ibn Tofail.

On obtient donc un autre système équivalent à AX = B :

a1,1 x1 + a1,2 x2 + ··· +a1,n xn = b1



 (2) (2) (2)
 0+ a2,2 x2 + ··· +a2,n xn = b2



(3) (3) (3)
0+ 0+ a3,3 x3 + · · · +a3,n xn = b3 . (S (3) )
 .. .. ..
. . .




 (3) (3) (3)
0+ 0+ an,3 x3 + · · · +an,n xn = bn .
On continue ce processus et on obtient à la fin un système triangulaire équivalent au
système de départ AX = B.

M
Exemple.
Résoudre avec la méthode de Gauss le système suivant ;

RI

 x + 2y + 3z =1
4x + 5y + 6z =2

Solution :



UK
7x + 8y + 8z = 3.

x + 2y + 3z = 1
E2 → E2 − 4E1 et E3 → E3 − 7E1 : −3y − 6z = −2
BO
−6y − 13z = −4.


x + 2y + 3z = 1
−6
E3 → E3 − ( −3 )E2 : −3y − 6z = −2
−z = 0.

 
f. L

−1/3
La solution est donc X =  2/3 .
0
Remarque.
Pro

Dans la description de la méthode de Gauss on a supposé qu’à chaque étape le

(k)
coéfficient ak,k ̸= 0, or il se peut qu’il soit non nul mais trés proche de zèro, ce qui
peut entrainer des erreurs d’arrondi importantes dans les calculs. Pour remédier à
(k) (k) (k)
ce problème, on choisit comme pivot l’élément ap,k tel que |ap,k | = max |ai,k | et on
k≤i≤n
permute les lignes p et k.

A.U : 2023-2024 4 Prof. L BOUKRIM

Cours d’Analyse numérique 1, Filière SMA, [Link] Ibn Tofail.

Exemple.
Résoudre avec la méthode de Gauss le système linéaire suivant :
 1
1016
x − y = −1
x + 2y = 1.
Autre remarque.
La méthode de Gauss sans pivotation s’applique lorsque tous les pivots sont non
nuls. Si ce n’est pas le cas, on cherche dans les équations suivantes un coéfficient non
nul. On choisit comme pivot celui le plus grand en valeur absolue et on éffectue des
permutations (de lignes et/ou colonnes) pour faire apparaı̂tre ce nouveau pivot à la

M
place du pivot nul.
On a deux types de pivot : Pivot partiel et pivot total.

RI
Méthode du pivot partiel :
(k−1) (k−1) (k−1)
tel que |ap,k | = max |ai,k | et on per-
mute les lignes k et p.
UK
On choisit comme pivot l’élément ap,k
k≤i≤n

Exemple.
BO
Résoudre avec la méthode de Gauss, , avec pivotation partielle, le système suivant ;


 2x + 4y − z + 5t = −10
x + 2y + 7t = −13


 x + y + 3z + t =4
2x + y + 2z + 4t = −5.
f. L

Méthode du pivot total :

(k−1) (k−1) (k−1)
On choisit comme pivot l’élément ap,l tel que |ap,l | = max |ai,j | et on
Pro

k≤i,j≤n
permute les lignes k et p puis les colonnes k et l.

Exemple.
Résoudre, avec pivotation totale, le système précédent.

A.U : 2023-2024 5 Prof. L BOUKRIM

Cours d’Analyse numérique 1, Filière SMA, [Link] Ibn Tofail.

Proposition. (Calcul du déterminant)

n
Y
m
detA = (−1) a(p−1)
p,p
p=1

où m est le nombre total de permutations (lignes et colonnes) effectuées.

Exemple. Calculer le déterminant de la matrice correspondante au système de

l’exemple précédent.

M
Proposition.
Le nombre d’opérations necessaires pour la résolution du système AX = B avec la
3
méthode de Gauss est de l’ordre de n3 .

RI
Autre remarque.
La méthode de Gauss devient impraticable dès que n est grand (exemple n = 100).
UK
BO
f. L
Pro

A.U : 2023-2024 6 Prof. L BOUKRIM