Méthodes directes de

I
LM
résolution des systèmes
linéaires

A
.S
(Cours N 04)
A
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 1 44
 Outline

I
LM
1. Méthode de factorisation LU

2. Application de la factorisation LU

A
2.1. Calcul du déterminant

.S
2.2. Résolution d’un système linéaire

3. Factorisation LU avec choix de pivot

A
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 2 44
 Méthode de factorisation LU
Définition :

I
LM
* La factorisation (décomposition) LU d’une matrice A d’ordre n con-
stitue un outil fondamental en analyse numérique.
* Elle consiste à transformer la matrice A en un produit de deux matrices

A
triangulaires :
A = LU
avec .S
• L une matrice triangulaire inférieure
A
• U une matrice triangulaire supérieure

A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 3 44
 Méthode de factorisation LU
    

I
a1,1 a1,2 · · · a1,n ℓ1,1 0 · · · 0 u1,1 u1,2 · · · u1,n

LM
 2,1 2,2
a a · · · a 2,n   2,1 2,2
ℓ ℓ · · · 0   0 u2,2 · · · u2,n 
 .. .. .. ..  =  .. .. .. ..   .. .. .. .. 
. . . . . . . . . . . .
an,1 an,2 · · · an,n ℓn,1 ℓn,2 · · · ℓn,n 0 0 · · · un,n

A
* Remarque :

.S
La factorisation LU impose que la diagonale de l’une des deux matrices L
ou U soit constituée d’éléments égaux à 1.
• la diagonale de L est unitaire ⇒ méthode de Doolittle
A
• la diagonale de U est unitaire ⇒ méthode de Crout
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 4 44
 Existence et unicité d’une factorisation LU

* Existence :

I
LM
• Une matrice carrée A admet une factorisation LU (sans pivot) si tous
les pivots de la méthode de Gauss sont non nuls.
• Autrement dit, les mineurs principaux de A sont non nuls.

A
det(A1 ) ̸= 0, det(A2 ) ̸= 0, . . . , det(An ) ̸= 0

* Unicité : .S
• La factorisation A = LU est unique si l’on impose que la matrice L
A
est triangulaire inférieure avec des éléments diagonaux égaux à 1.

A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 5 44
 Existence et unicité d’une factorisation LU
Exemple : Soit la matrice carrée A suivante :

I
 

LM
2 4 3
A = 2 2 −1
1 −1 2

A
- Justifier l’existence de la factorisation LU de la matrice A.
Solution :
.S
- Vérification de l’existence de LU : Une matrice carrée A admet une
factorisation LU si tous les mineurs principaux de A sont différents de zéro:
A
det(A1 ) ̸= 0, det(A2 ) ̸= 0, det(A3 ) ̸= 0
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 6 44
 Existence et unicité d’une factorisation LU

det(A1 ) = 2 ̸= 0

I
LM
2 4
det(A2 ) = = −4 ̸= 0
2 2

2 4 3

A
2 −1 2 −1 2 2
det(A3 ) = det(A) = 2 2 −1 = 2 −4 +3
−1 2 1 2 1 −1
1 −1 2
.S
= 2(4 − 1) − 4(4 + 1) + 3(−2 − 2)
= 2(3) − 4(5) + 3(−4) = 6 − 20 − 12 = −26 ̸= 0
A
⇒ factorisation LU existe.
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 7 44
 Méthode de factorisation LU
* Dans la suite, on adoptera le cas où la matrice L est unitaire – Méthode

I
de Doolittle – (car c’est le cas le plus courant en pratique).

LM
 
a1,1 a1,2 · · · a1,n
 a2,1 a2,2 · · · a2,n 
A=  ... .. .. 
. 
...
.

A
an,1 an,2 · · · an,n
   

L=  ..
1

.
0 ··· 0
 ℓ2,1 1 · · · 0 
. .
.. . . . .. 

.S U =
u1,1 u1,2 · · · u1,n
 0 u2,2 · · · u2,n 
 .
 .. .. .. 
. 
...
.
A
ℓn,1 ℓn,2 · · · 1 0 0 · · · un,n

A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 8 44
 Méthode de factorisation LU
* On suppose que la matrice A admet une factorisation A = LU . Les

I
coefficients de L et U sont déterminés par identification :

LM
X
n
aij = ℓik ukj
k=1

A
* Calcul des coefficients :
X
i−1
Pour la matrice U : .S uij = aij −
k=1
ℓik ukj (i ≤ j)

!
X
j−1
A
1
Pour la matrice L : ℓij = aij − ℓik ukj (i > j)
ujj
k=1
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 9 44
 Méthode de factorisation LU

* La méthode de factorisation LU se déroule comme suit :

I
LM
1. Initialiser L comme matrice identité.

2. Calculer les coefficients de U ligne par ligne.

A
3. Calculer les coefficients de L colonne par colonne.

.S
4. Répéter jusqu’à obtenir toutes les valeurs de L et U .
A
Important : les pivots ujj doivent être non nuls.

A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 10 44
 Méthode de factorisation LU

I
Exemple : Soit la matrice carrée A suivante :

LM
 
2 4 3
A = 2 2 −1
1 −1 2

A
- Décomposer la matrice A sous la forme A = LU (avec L unitaire).

Solution:
.S
A
Dans l’exemple précédent on a justifier l’existence de la factorisation LU
de la matrice A.
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 11 44
 Exemple

I
Factorisation: Décomposition de la matrice A sous la forme A = LU

LM
avec L unitaire :
    
2 4 3 1 0 0 u1,1 u1,2 u1,3
 2 2 −1  = ℓ2,1 1 0  0 u2,2 u2,3 
1 −1 2

A
ℓ3,1 ℓ3,2 1 0 0 u3,3
 

=
ℓ
u1,1
2,1 u 1,1 ℓ 2,1 u
.S
u1,2
1,2 + u 2,2 ℓ 2,1 u
u1,3
1,3 + u 2,3



A
ℓ3,1 u1,1 ℓ3,1 u1,2 + ℓ3,2 u2,2 ℓ3,1 u1,3 + ℓ3,2 u2,3 + u3,3

A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 12 44
 Exemple

I
Par identification LU = A, on obtient :

LM


u1,1 = 2
* de la 1ére ligne ⇒
u1,2 = 4



A
u1,3 = 3

* de la 2ème ligne ⇒

 .S
ℓ2,1 u1,1 = 2
ℓ2,1 u1,2 + u2,2 = 2 ⇒


ℓ2,1 = 1
u2,2 = −2
 
A
 
ℓ2,1 u1,3 + u2,3 = −1 u2,3 = −4

A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 13 44
 Exemple
 

I
 1

 ℓ3,1 u1,1 = 1 
 ℓ3,1 =

LM

 
 2
* de la 3ème ligne ⇒ ⇒
3
ℓ3,1 u1,2 + ℓ3,2 u2,2 = −1 ℓ3,2 =

 
 2

 

ℓ3,1 u1,3 + ℓ3,2 u2,3 + u3,3 = 2 u = 13

A
3,3
2
Matrices L et U obtenues

1 0
1 1
.S :
0
0



2 4 3


L=  U = 0 −2 −4
A
1 3 13
1 0 0
2 2 2
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 14 44
 Outline

I
LM
1. Méthode de factorisation LU

2. Application de la factorisation LU

A
2.1. Calcul du déterminant

.S
2.2. Résolution d’un système linéaire

3. Factorisation LU avec choix de pivot

A
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 15 44
 Méthode de factorisation LU

I
* Application de la méthode de factorisation LU :

LM
Cette méthode est un outil fondamental en analyse numérique, qui permet:

• Le calcul du déterminant d’une matrice.

A
• La résolution de systèmes linéaires.

.S
• La résolution de plusieurs systèmes avec la même matrice A si plusieurs
seconds membres b(1) , b(2) , . . . sont donnés.
A
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 16 44
 Outline

I
LM
1. Méthode de factorisation LU

2. Application de la factorisation LU

A
2.1. Calcul du déterminant

.S
2.2. Résolution d’un système linéaire

3. Factorisation LU avec choix de pivot

A
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 17 44
 Calcul du déterminant d’une matrice
* Le calcul du déterminant se fait en multipliant les déterminants de L et

I
LM
de U :
det(A) = det(LU )

= det(L) det(U )

A
* Comme la matrice L est triangulaire unitaire :

Alors
.S det(L) = 1
A
det(A) = det(U )

A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 18 44
 Méthode de factorisation LU
Exemple : Soit la matrice carrée A et sa décomposition LU suivante :

I
     

LM
2 4 3 1 0 0 2 4 3
   
A = 2 2 −1
 L =  1 1 0 U = 0 −2 −4
1 3 13
1 −1 2 1 0 0
2 2 2

A
- Calcul du déterminant de A.

Solution : .S
det(A) = det(L) det(U ) = det(U ) = 2×(−2)×
13
= −26 det(A) = −26
A
2
(det(L) = 1 car L est triangulaire unitaire)
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 19 44
 Outline

I
LM
1. Méthode de factorisation LU

2. Application de la factorisation LU

A
2.1. Calcul du déterminant

.S
2.2. Résolution d’un système linéaire

3. Factorisation LU avec choix de pivot

A
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 20 44
 Résolution d’un système linéaire
* La résolution d’un système linéaire par la méthode LU se déroule en

I
trois principales étapes:

LM
1. Factorisation : Décomposition de la matrice du système sous la forme

A = LU

A
Une fois la décomposition effectuée, le système linéaire :
.S
Ax = b ⇐⇒ L |{z}
Ux = b
y
A
on obtient le système Ly = b
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 21 44
 Résolution d’un système linéaire

I
2. Résolution du système triangulaire inférieur :

LM
Ly = b

par substitution avant pour déterminer y.

A
3. Résolution du système triangulaire supérieur :
.S Ux = y
A
par substitution arrière pour déterminer x.

A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 22 44
 Exemple 1
Exemple : Soit le système d’équations linéaires suivant :

I


LM

 2x1 + x2 + 2x3 = 10
6x1 + 4x2 = 26

 8x + 5x + x = 35
1 2 3

A
1. Écrire le système ci-dessus sous la forme matricielle Ax = b ?
.S
2. Justifier l’existence de la factorisation LU de la matrice A du système
Ax = b ?
A
3. Résoudre ce système par la méthode de factorisation LU ?
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 23 44
 Exemple 1
Solution :

I
LM
1. Mise du système sous la forme matricielle:
     
2 1 2 x1 10
Ax = b ⇐⇒ A = 6 4 0   
x = x2  b = 26


A
8 5 1 x3 35

2. Vérification de l’existence de LU :
.S
Une matrice carrée A admet une factorisation LU sans pivot si tous les
mineurs principaux de A sont différents de zéro :
A
det(A1 ) ̸= 0, det(A2 ) ̸= 0, det(A3 ) ̸= 0
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 24 44
 Exemple 1

det(A1 ) = 2 ̸= 0

I
LM
2 1
det(A2 ) = = 2 ̸= 0
6 4
2 1 2

A
4 0 6 0 6 4
det(A3 ) = det(A) = 6 4 0 = 2 −1 +2
5 1 8 1 8 5
8 5 1
.S
= 2(4 − 0) − 1(6 − 0) + 2(30 − 32)
= 2(4) − 6 + 2(−2) = −2 ̸= 0
A
⇒ factorisation LU existe.
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 25 44
 Exemple 1
3. Résolution du système par la méthode de factorisation LU

I
LM
a. Factorisation: Décomposition de A sous la forme A = LU (L unitaire):
    
2 1 2 1 0 0 u1,1 u1,2 u1,3
 6 4 0  = ℓ2,1 1 0  0 u2,2 u2,3 

A
8 5 1 ℓ3,1 ℓ3,2 1 0 0 u3,3
 

=

 ℓ
u1,1
2,1 u 1,1 ℓ 2,1 u
.S
u1,2
1,2 + u 2,2 ℓ 2,1 u
u1,3
1,3 + u 2,3



A
ℓ3,1 u1,1 ℓ3,1 u1,2 + ℓ3,2 u2,2 ℓ3,1 u1,3 + ℓ3,2 u2,3 + u3,3

A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 26 44
 Exemple 1

I
Par identification LU = A, on obtient :

LM


u1,1 = 2
* de la 1ére ligne ⇒
u1,2 = 1



A
u1,3 = 2

* de la 2ème ligne ⇒

 .S
ℓ2,1 u1,1 = 6
ℓ2,1 u1,2 + u2,2 = 4 ⇒


ℓ2,1 = 3
u2,2 = 1
 
A
 
ℓ2,1 u1,3 + u2,3 = 0 u2,3 = −6

A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 27 44
 Exemple 1
 

I

ℓ3,1 u1,1 = 8 
ℓ3,1 = 4

LM
* de la 3ème ligne ⇒
ℓ3,1 u1,2 + ℓ3,2 u2,2 = 5 ⇒ ℓ3,2 = 1

 

ℓ3,1 u1,3 + ℓ3,2 u2,3 + u3,3 = 1 u3,3 = −1

A
Matrices L et U obtenues :

1 0 0
L = 3 1 0
 .S 
2 1 2
U = 0 1 −6

A
4 1 1 0 0 −1

A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 28 44
 Exemple 1
b. Résolution du système Ly = b par substitution avant:

I


LM
 y1 = 10
3y1 + y2 = 26

4y1 + y2 + y3 = 35
• de la 1ère équation, on a y1 = 10

A
• de la 2ème équation, on a 3(10) + y2 = 26 ⇒ y2 = −4

La solution de Ly = b est :    
y1 10
.S
• de la 3ème équation, on a 4(10) + (−4) + y3 = 35 ⇒ y3 = 1
A
y2  = −4
y3 1
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 29 44
 Exemple 1
c. Résolution du système U x = y par substitution arrière:

I


LM
 2x1 + x2 + 2x3 = 10
x2 − 6x3 = −4

−x3 = −1
• de la 3ème équation, on a −x3 = −1 ⇒ x3 = 1

A
• de la 2ème équation, on a x2 − 6(1) = −4 ⇒ x2 = 2

La solution du système est : x  3

1
.S
• de la 1ère équation, on a 2x1 + (2) + 2(1) = 10 ⇒ 2x1 = 6 ⇒ x1 = 3
A
x2  = 2
x3 1
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 30 44
 Exemple 2
Exemple : Soit le système d’équations linéaires suivant :

I


LM
 2x1 + x2 − 2x3 = 1
4x1 + 5x2 − 3x3 = 6

−2x1 + 5x2 + 3x3 = 6

A
Résoudre ce système par la méthode de factorisation LU ?


2 1 −2


.S
Solution : Ce système est représenté par la forme matricielle:
  

x1

 
1
Ax = b ⇐⇒ A = 4 5 −3 x = x2 b = 6
A
−2 5 3 x3 6

A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 31 44
 Exemple 2

I
1. Factorisation: Décomposition de A sous la forme A = LU (L unitaire):

LM
    
2 1 −2 1 0 0 u1,1 u1,2 u1,3
 4 5 −3 = ℓ2,1 1 0  0 u2,2 u2,3 
−2 5 3 ℓ3,1 ℓ3,2 1 0 0 u3,3

A
 
u1,1 u1,2 u1,3
=

 ℓ 2,1 u 1,1 ℓ 2,1 u .S
1,2 + u 2,2 ℓ 2,1 u 1,3 + u 2,3

ℓ3,1 u1,1 ℓ3,1 u1,2 + ℓ3,2 u2,2 ℓ3,1 u1,3 + ℓ3,2 u2,3 + u3,3



A
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 32 44
 Exemple 2
Par identification LU =A, on obtient :

I

u1,1 = 2

LM
* de la 1ère ligne ⇒
u1,2 = 1


u1,3 = −2
 
 

A
ℓ2,1 u1,1 = 4 ℓ2,1 = 2
* de la 2ème ligne ⇒ ⇒
ℓ2,1 u1,2 + u2,2 = 5 u2,2 = 3
 



ℓ3,1 u1,1 = −2
.S
ℓ2,1 u1,3 + u2,3 = −3

u2,3 = 1


ℓ3,1 = −1
* de la 3 ème ⇒
A
ligne ⇒ ℓ3,2 = 2
ℓ3,1 u1,2 + ℓ3,2 u2,2 = 5

 

ℓ3,1 u1,3 + ℓ3,2 u2,3 + u3,3 = 3 u3,3 = −1
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 33 44
 Exemple 2

I
     

LM
1 0 0 u1,1 u1,2 u1,3 
u1,1 = 2 
ℓ2,1 = 2 
ℓ3,1 = −1
ℓ2,1 1 0  0 u2,2 u2,3  u1,2 = 1 u2,2 = 3 ℓ3,2 = 2

 
 

ℓ3,1 ℓ3,2 1 0 0 u3,3 u1,3 = −2 u2,3 = 1 u3,3 = −1

A
Matrices L et U obtenues :

1 0 0
.S
L =  2 1 0
 
2 1 −2
U = 0 3 1 

A
−1 2 1 0 0 −1

A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 34 44
 Exemple 2
2. Résolution du système Ly = b par substitution avant:

I


LM
 y1 = 1
2y1 + y2 = 6

−y1 + 2y2 + y3 = 6
• de la 1ère équation, on a y1 = 1

A
• de la 2ème équation, on a 2(1) + y2 = 6 ⇒ y2 = 4
• de la 3ème équation, on a −1 + 2(4) + y3 = 6 ⇒ y3 = −1
La solution de Ly = b est :    
y1 1
.S
y 2  =  4 
A
y3 −1
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 35 44
 Exemple 2
3. Résolution du système U x = y par substitution arrière:

I

 2x1 + x2 − 2x3

LM
= 1
3x2 + x3 = 4

−x3 = −1
• de la 3ème équation, on a x3 = 1

A
• de la 2ème équation, on a 3x2 + 1 = 4 ⇒ x2 = 1

La solution du système est : x  1

1
.S
• de la 1ère équation, on a 2x1 + 1 − 2 = 1 ⇒ x1 = 1

x2  = 1
A
x3 1
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 36 44
 Outline

I
LM
1. Méthode de factorisation LU

2. Application de la factorisation LU

A
2.1. Calcul du déterminant

.S
2.2. Résolution d’un système linéaire

3. Factorisation LU avec choix de pivot

A
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 37 44
 Factorisation LU avec choix de pivot
* La factorisation LU classique nécessite que les pivots soient non nuls.

I
Lorsque ce n’est pas le cas, on utilise une permutation des lignes.

LM
* On introduit alors une matrice de permutation P (obtenue en échangeant
des lignes de la matrice identité) telle que :

A
P A = LU

* Exemple : .S 
0 1 0
P = 1 0 0


0 0 1
A
* Multiplier P par A revient à échanger les lignes de A.
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 38 44
 Résolution par la méthode LU avec échange
* En utilisant la matrice de permutation P A = LU :

I
LM
P Ax = P b ⇐⇒ LU x = P b ⇐⇒ L |{z}
Ux = Pb
y

* Une fois la décomposition effectuée, la résolution s’effectue par la:

A
• Résolution du système triangulaire inférieur : Ly = P b

.S
• Résolution du système triangulaire supérieur : Ux = y

Remarque : Lorsqu’on utilise une factorisation avec pivot, il faut

A
appliquer la même permutation au second membre b.

A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 39 44
 Exemple
Exemple : Factoriser par la méthodeLU la 
matrice A:

I
2 1 1

LM
A= 4 2 3
2 5 1

Solution : Décomposition de A sous la forme A = LU avec L unitaire :

A
    
2 1 1 1 0 0 u1,1 u1,2 u1,3
4 2
2 5 1

ℓ3,1 ℓ3,2 1
u1,1
.S
3 = ℓ2,1 1 0  0 u2,2 u2,3 
0
u1,2
0 u3,3
u1,3

A
= ℓ2,1 u1,1 ℓ2,1 u1,2 + u2,2 ℓ2,1 u1,3 + u2,3 
ℓ3,1 u1,1 ℓ3,1 u1,2 + ℓ3,2 u2,2 ℓ3,1 u1,3 + ℓ3,2 u2,3 + u3,3
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 40 44
 Exemple
Par identification LU = A :

I


LM

u1,1 = 2
* 1ère ligne ⇒
u1,2 = 1


u1,3 = 1
 

A

ℓ2,1 u1,1 = 4 
ℓ2,1 = 2
* 2ème ligne ⇒

 .S
ℓ2,1 u1,2 + u2,2 = 2
ℓ2,1 u1,3 + u2,3 = 3
⇒ u2,2 = 0


u2,3 = 1
A
s Remarque : On obtient u2,2 = 0 ⇒ pivot nul ⇒ La factorisation LU
sans pivot est impossible.
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 41 44
 Exemple

* Dans ce cas, il faut utiliser une permutation P A = LU

I
LM
* On applique alors une permutation de lignes. On cherche un pivot non
(2)
nul dans la colonne 2 : a32 = 4 ̸= 0.

A
* On permute les lignes l2 et l3
   
.S
1 0 0
P = 0 0 1
2 1 1
P A = 2 5 1
0 1 0 4 2 3
A
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 42 44
 Exemple
On cherche maintenant la décomposition : P A = LU

I
    

LM
2 1 1 1 0 0 u1,1 u1,2 u1,3
2 5 1 = ℓ2,1 1 0  0 u2,2 u2,3 
4 2 3 ℓ3,1 ℓ3,2 1 0 0 u3,3
 
u1,1 u1,2 u1,3

A
= ℓ2,1 u1,1 ℓ2,1 u1,2 + u2,2 ℓ2,1 u1,3 + u2,3 

Par identification :

.S
ℓ3,1 u1,1 ℓ3,1 u1,2 + ℓ3,2 u2,2 ℓ3,1 u1,3 + ℓ3,2 u2,3 + u3,3

u1,1 = 2
* 1ère ligne ⇒
A
u1,2 = 1


u1,3 = 1
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 43 44
 Exemple
 

 ℓ u = 2 
ℓ2,1 = 1

I
2,1 1,1
* 2ème ligne ⇒
ℓ2,1 u1,2 + u2,2 = 5 ⇒ u2,2 = 4

LM

 

ℓ2,1 u1,3 + u2,3 = 1 u2,3 = 0
 

 3,1 1,1
ℓ u = 4 
ℓ3,1 = 2

A
* 3ème ligne ⇒
ℓ3,1 u1,2 + ℓ3,2 u2,2 = 2 ⇒ ℓ3,2 = 0

 


Matrices obtenues :

.S
ℓ3,1 u1,3 + ℓ3,2 u2,3 + u3,3 = 3

1 0 0
 
2 1 1

u3,3 = 1
A
L = 1 1 0 U = 0 4 0
2 0 1 0 0 1
A. SALMI (Université Mouloud Mammeri de Tizi-Ouzou) Méthodes Numériques Année universitaire : 2025–2026 44 44