Résolution numérique de systèmes d’équations linéaires

Méthode de Gauss-Seidel

AN - 3A&B - A.U. 2022/2023

 2

La Méthode de Gauss-Seidel
Comme pour la méthode de Jacobi, la méthode de Gauss-Seidel pour la résolution
d’un système d’équations linéaires (S) : AX = b, consiste en premier lieu à
décomposer A sous la forme:
A = D − E − F,
où D est une matrice diagonale, E est une matrice triangulaire inférieure et F est
une matrice triangulaire supérieure.

@UP-Maths Résolution numérique des systèmes d’équations linéaires MNI

 2

La Méthode de Gauss-Seidel
Comme pour la méthode de Jacobi, la méthode de Gauss-Seidel pour la résolution
d’un système d’équations linéaires (S) : AX = b, consiste en premier lieu à
décomposer A sous la forme:
A = D − E − F,
où D est une matrice diagonale, E est une matrice triangulaire inférieure et F est
une matrice triangulaire supérieure.
     
a1,1 0 ··· 0 0 ··· ··· 0 0 −a1,2 ··· −a1,n
 .. ..
  ..  . .. ..  ..
. . . . . .
0   −a2,1
0  0
 . .
  
A= −  − .
 
.. .. .. 
 ..
 .. .. . .. 

 . . . 0 
  . . . 0
 . . −an−1,n 

0 · · · 0 an,n −an,1 · · · −an,n−1 0 0 ··· ··· 0
| {z } | {z } | {z }
D E F

@UP-Maths Résolution numérique des systèmes d’équations linéaires MNI

 3

Considérons, par exemple, le cas où n = 3. On a

 
a1,1 a1,2 a1,3
 
A = 
 a2,1 a2,2 a2,3


a3,1 a3,2 a3,3
     
a1,1 0 0 0 0 0 0 −a1,2 −a1,3
     
= 
 0 a2,2 0  − −a
  2,1 0 0  −
 
0 0 −a 
2,3 

0 0 a3,3 −a3,1 −a3,2 0 0 0 0

| {z } | {z } | {z }
D E F

Le système (S) : AX = b est équivalent alors à

(D − E)X − F X = b
⇐⇒ (D − E)X = F X + b

@UP-Maths Résolution numérique des systèmes d’équations linéaires MNI

 4

4
 
(k)
x1
 (k) 
Soit (X (k) )k≥0 la suite de vecteurs dans R3 définie par X (k) = x ,
 2  vérifiant
(k)
x3

(D − E) X (k+1) = |{z}
F X (k) + b
| {z }
M N

Si A est à diagonale strictement dominante, alors les coefficients diagonaux de A

sont non nuls. Par conséquent, M est inversible. Dans ce cas
X (k+1) = M −1 N X (k) + M −1 b

@UP-Maths Résolution numérique des systèmes d’équations linéaires MNI

 4

4
 
(k)
x1
 (k) 
Soit (X (k) )k≥0 la suite de vecteurs dans R3 définie par X (k) = x ,
 2  vérifiant
(k)
x3

(D − E) X (k+1) = |{z}
F X (k) + b
| {z }
M N

Si A est à diagonale strictement dominante, alors les coefficients diagonaux de A

sont non nuls. Par conséquent, M est inversible. Dans ce cas
X (k+1) = M −1 N X (k) + M −1 b

Remarque
Si la suite (Xk )k≥0 est convergente, alors

lim X (k) = X,
k→+∞

avec X l’unique solution du système (S).

@UP-Maths Résolution numérique des systèmes d’équations linéaires MNI
 5

Les composantes du vecteur X (k+1) s’écrivent en fonction des composantes du

vecteur X (k) comme suit:
   (k+1)     (k)   
a1,1 0 0 x1 0 −a1,2 −a1,3 x b1
   (k+1)     1(k)   
 = 0
 x2  + b2 
−a2,3 
a2,1 a2,2 0   x2 0
    
  
(k+1) (k)
a3,1 a3,2 a3,3 x3 0 0 0 x3 b3
Ainsi, on en déduit que  (k) (k)
(k+1) (k) (k) (k+1) b1 −a1,2 x2 −a1,3 x3
x1 =

a1,1 x1 + a1,2 x2 + a1,3 x3 = b1

a1,1
 

 (k+1) (k)
(k+1) (k+1) (k) b2 −a2,1 x1 −a2,3 x3
a2,1 x1 + a2,2 x2 + a2,3 x3 = b2 ⇐⇒ x(k+1)
2 = a2,2
(k+1) (k+1) (k+1)
  (k+1) (k+1)
a3,1 x1 + a3,2 x2 + a3,3 x3 = b3 x(k+1) = b3 −a3,1 x1 −a3,2 x2
 

3 a3,3

@UP-Maths Résolution numérique des systèmes d’équations linéaires MNI

 6

(k+1)
Dans Rn , les composantes xi (i ∈ {0, 1, · · · , n}) du vecteur
(k+1) (k)
X s’écrivent en fonction des composantes xi du vecteur X (k)
comme suit:
i−1 n
!
(k+1) 1 X (k+1)
X (k)
xi = bi − aij xj − aij xj , ∀i ∈ {1, · · · , n}
aii j=1 j=i+1

@UP-Maths Résolution numérique des systèmes d’équations linéaires MNI

 7

Convergence de la méthode de Gauss-Seidel

Question: Existe-il une condition sur la matrice A assurant la convergence de la
suite (X (k) )k≥0 est convergente?

@UP-Maths Résolution numérique des systèmes d’équations linéaires MNI

 7

Convergence de la méthode de Gauss-Seidel

Question: Existe-il une condition sur la matrice A assurant la convergence de la
suite (X (k) )k≥0 est convergente?

Théorème
Soit A une matrice à diagonale strictement dominante alors la méthode de
Gauss-Seidel appliquée au système (S) : AX = b est convergente vers la solution de
(S) pour tout X (0) ∈ Mn,1 (R).

@UP-Maths Résolution numérique des systèmes d’équations linéaires MNI

 7

Convergence de la méthode de Gauss-Seidel

Question: Existe-il une condition sur la matrice A assurant la convergence de la
suite (X (k) )k≥0 est convergente?

Théorème
Soit A une matrice à diagonale strictement dominante alors la méthode de
Gauss-Seidel appliquée au système (S) : AX = b est convergente vers la solution de
(S) pour tout X (0) ∈ Mn,1 (R).

Remarque
Comme pour la méthode de Jacobi, on peut considérer le critère d’arrêt suivant pour
la méthode de Gauss-Seidel:
||AX (k) − b|| ≤ ε, avec la tolérence ε assez petite.
@UP-Maths Résolution numérique des systèmes d’équations linéaires MNI
 8

Étude d’un exemple

On considère un système d’équations linéaires (S), telle que;

(S) ⇔ AX = b
avec :      
5 2 −1 x1 6
     
1 6 −3 , X = x2  et b = 4
A=     

2 1 4 x3 7

1 Montrer qu’il existe une unique solution de (S) dans R3 .

2 Etudier la convergence de la méthode de Gauss-Seidel pour la résolution de (S).
3 Donner le schéma itératif de la méthode de Gauss-Seidel associé à (S).
4 Calculer les trois premiers itérés par la méthode de Gauss-Seidel.

@UP-Maths Résolution numérique des systèmes d’équations linéaires MNI

 9

• Existence d’une unique solution de (S) dans R3 :

On a det(A) = 126 6= 0 ⇒ ∃! X ∈ R3 / AX = b.

@UP-Maths Résolution numérique des systèmes d’équations linéaires MNI

 10

10

• Etude de la convergence de la méthode de Gauss-Seidel pour la résolution de (S):

On a: 
|a1,1 | > |a1,2 | + |a1,3 | (car |5| > |2| + | − 1|)

|a2,2 | > |a2,1 | + |a2,3 | (car |6| > |1| + | − 3|)

|a3,3 | > |a3,1 | + |a3,2 | (car |4| > |2| + |1|)


⇒ A est une matrice à diagonale strictement dominante

⇒ La méthode de Gauss-Seidel est convergente.

@UP-Maths Résolution numérique des systèmes d’équations linéaires MNI

 11

11

• Schéma itératif associé à (S) avec la méthode de Gauss-Seidel:

   (k+1)     (k)   
5 0 0 x1 0 −2 1 x 6
   (k+1)     1(k)   
1 6 0 x  = 0 0 3 x  + 4
  2    2   
(k+1) (k)
2 1 4 x3 0 0 0 x3 7

@UP-Maths Résolution numérique des systèmes d’équations linéaires MNI

 11

11

• Schéma itératif associé à (S) avec la méthode de Gauss-Seidel:

   (k+1)     (k)   
5 0 0 x1 0 −2 1 x 6
   (k+1)     1(k)   
1 6 0 x  = 0 0 3 x  + 4
  2    2   
(k+1) (k)
2 1 4 x3 0 0 0 x3 7
 (k) (k) (k) (k)
(k+1) b1 −a1,2 x2 −a1,3 x3 6−2x2 +x3
 x = =
 1 a1,1 5


(k+1) (k) (k+1) (k)
(k+1) b2 −a2,1 x1 −a2,3 x3 4−x1 +3x3
⇒ x2 = a2,2 = 6
 (k+1) (k+1) (k+1) (k+1)
x(k+1) = b3 −a3,1 x1 −a3,2 x2 7−2x1 −x2

=

3 a3,3 4

@UP-Maths Résolution numérique des systèmes d’équations linéaires MNI

 12

12

• Application de la méthode de Gauss-Seidel avec 3 itérations:

 
0
 
Considérons par exemple un vecteur initial X (0) = 
0.


@UP-Maths Résolution numérique des systèmes d’équations linéaires MNI

 12

12

• Application de la méthode de Gauss-Seidel avec 3 itérations:

 
0
 
Considérons par exemple un vecteur initial X (0) = 
0.


0
 
6/5
 
I Itération 1 : X (1) =  7/15 
 
31/30

@UP-Maths Résolution numérique des systèmes d’équations linéaires MNI

 12

12

• Application de la méthode de Gauss-Seidel avec 3 itérations:

 
0
 
Considérons par exemple un vecteur initial X (0) = 
0.


0
 
6/5
 
I Itération 1 : X (1) =  7/15 
 
31/30
 
1, 22
 
I Itération 2 : X (2) =  0.98 
 
0, 895

@UP-Maths Résolution numérique des systèmes d’équations linéaires MNI

 12

12

• Application de la méthode de Gauss-Seidel avec 3 itérations:

 
0
 
Considérons par exemple un vecteur initial X (0) = 
0.


0
 
6/5
 
I Itération 1 : X (1) =  7/15 
 
31/30
 
1, 22
 
I Itération 2 : X (2) =  0.98 
 
0, 895
 
1, 014
 
I Itération 3 : X (3) =  1.014 
 
0.989

@UP-Maths Résolution numérique des systèmes d’équations linéaires MNI