Exercice 1

1. Méthode itérative
(a) Résolvons le système Ax = B avec

2 −1 0 0
 
−1 2 −1 0 
A= 
 0 −1 2 −1
 

0 0 1 2
 T
et b = 3 −5 7 −6

- La méthode de LU :

Soit : L = (lij )1≤i,j≤n tels que lij = 0 si i < j et

U = (uij )1≤i,j≤n tels que uij = 0 si i > j et uij = 1 si i = j

l11 0 0 0 1 u12 u13 u14

   
l
21 l22 0 0 
0 1 u23 u24 
 
L=  et U = 
l31 l32 l33 0  0 0 1 u34 
 

l41 l42 l43 l44 0 0 0 1

X j−1
X
aij = lik ukj , 1 ≤ i, j ≤ n; lij = aij − lik ukj , i < j; i = 1, ..., n
k=1

i−1
1 X
uij = [aij − lik ukj ], i < j, i = 1, ..., n
lij k=1

Application numérique :

l11 = a11 = 2; l21 = a21 = −1; l22 = 32 ; l32 = −1; l33 = 34 ; l43 = 1; l44 = 11
4

l41 = l31 = l42 = 0

u12 = − 12 ; u23 = − 32 ; u34 = − 34

u13 = u14 = u24 = 0

2 0 0 0 1 − 21 0 0
   
−1 3 2
2
0 0   0 1 − 3
0 
Nous obtenons : L =   et U = 
  
 0 −1 43 0 0 0 1 − 34 

11
0 0 1 4
0 0 0 1

1
 
Ly =b
Ax = b ⇐⇒ (LU )x = b ⇐⇒
U x =y
2 0 0 0 y1 3
    
−1 3
2
0 0  y2  −5
   
Ly = b =⇒    = 
 
 0 −1 43 0  y3   7 

11
0 0 1 4
y4 −6




2y1
=3
3

−y
1 + 2 y2 = −5

=⇒




−y2 + 34 y3 = 7
y3 + 11 y = −6


4 4


= 32


y1
7

2 = −3
y

=⇒




y3 = 27
y4 = − 38


11

1 − 12 0 0 x1 3
    
2
0 1 − 32 0   x2 
   −7 
 3
U x = y =⇒    = 

0 0 1 − 34  x3   72 

38
0 0 0 1 x4 − 11

− 21 x2 = 32


x1
2 7

2 − 3 x3 = − 3
x

=⇒




x3 − 43 x4 = 72
x4 = − 38


11


7
= 11


x1
19

2 = − 11
x

=⇒




x3 = 10
11
x4 = − 38


11
7
Alors x1 = 11
; x2 = − 19
11
; x3 = 10
11
: x4 = − 38
11

2
- La méthode de Choleski :

Soit A une matrice donnée il existe une matrice triangulaire inférieure

telle que A = LLT si seulement si A est une matrice symétrique définie
positive.
Dans notre cas A est non symétrique car AT ̸= A , alors nous allons multi-
plier le système par la transposé de A.
Ax = B ⇐⇒ AT Ax = AT b.

5 −4 1 0
 
−4 6 −4 1
La matrice AT A =   est symétrique.

 1 −4 6 0


0 1 0 5
5 −4 1 0
5 −4 1
5 −4 −4 6 −4 1
De plus 5 > 0 ; = 14 > 0 ; −4 6 −4 = 30 > 0 ; =
−4 6 1 −4 6 0
1 −4 6
0 1 0 5
121 > 0
La matrice AT A est bien symétrique et définie positive. Il suffit de poser
AT Ax = AT b pour résoudre le système par la méthode de Choleski.

l11 0 0 0
 
v
i−1 Pi−1
aij − k=1
u
 21 l22 0
l 0  u X
2 lik ljk
Soit L =   ; lii = aii −
t lik ; lji =
l31 l32 l33 0  lii
k=1
l41 l42 l43 l44
√ √ √ √ √ √
l11 =√ 5; l21 = √− 4 5 5 ; l22 = √ 570 ; l31 = 5
5
; l32 = − 8 3570 ; l33 = 105
7
; l41 =
0; l42 = 1470 ; l43 = 8 105
105
; l44 = 113030
 √ 
5 0 0 0
 4√5 √
70
− 0 0

 √5 5√

Alors L =  √ 
5
 5
−√8 3570 √7
105
0
√


70 11 308 105
0 14 30 105
T
En résolvant d’abord Ly = A b ; on trouve y1 = 4, 9193; y2 = −6, 6933; y3 =
−1, 3663; y4 = −6, 9378

Par la suite en résolvant LT x = AT b ; on trouve comme solution x1 =

7
11
= − 19
; x2 11
; x3 = 10
11
38
: x4 = − 11

3
(b) Reprenons le (a) par la méthode de Jacobi et de Gauss-Seidel
- La méthode de jacobi



 x1 − 12 x2 = 32
2 7

2 − 3 x3 = − 3
x

D’après la résolution précédente on a : 



x3 − 34 x4 = 72
x4 = − 38


11
Le système récursif pour la méthode de jacobi s’écrit :


2 −1 0 0 x1
  
3
 

 2x1 − x2 = 3

−1 2 −1 0  x  −5 −x + 2x − x = −5

  2 1 2 3
  =   =⇒
  
 0 −1 2 −1 x3   7 



 −x2 + 2x3 − x4 = 7
0 0 1 2 x4 −6

x3 + 2x4 = −6



3 + x2

x1 =



2


−5 + x1 + x3



2 =
x


=⇒  2
7 + x2 + x4

x3 =

 2
−6 − x3



x4 =


2

3 + xk2

xk+1 =


1

2



−5 + xk1 + xk3


xk+1

=


2
=⇒  2
7 + xk2 + xk4
xk+1 =


 3


2
−6 − xk3



xk+1


4 =
2
- La méthode de Gauss-Seidel



 x1 − 12 x2 = 32
2 7

2 − 3 x3 = − 3
x

D’après la résolution précédente on a : 

 x − 34 x4 = 72
 3
x4 = − 38


11
Le système récursif pour la méthode de Gauss-Seidel s’écrit :


2 −1 0 0 x1
  
3
 

 2x1 − x2 = 3

−1 2 −1 0  x  −5 −x + 2x − x = −5

  2 1 2 3
  =   =⇒
  
 0 −1 2 −1 x3   7 



 −x2 + 2x3 − x4 = 7
0 0 1 2 x4 −6

x3 + 2x4 = −6



4
 3 + x2

x1 =



2


−5 + x1 + x3



2 =
x


=⇒ 2
7 + x2 + x4
x3 =




 2
−6 − x3



x4 =


2

3 + xk2

xk+1 =


 1


 2
−5 + xk+1 + xk3


1
xk+1

=


2
=⇒  2
k+1 7 + xk+1
2 + xk4
x =


 3


 2k+1
−6 − x


xk+1 =

 3
4
2
(c) Faisons un tableau contenant les 6 premières itérations des deux méthodes avec les
erreurs
- La méthode de jacobi

⃗ (0) = [0 0 0 0]T , et Erreur = ∥XLU − Xjacobi ∥ On trouve :

A partir de X

1 2 3 4 5 6
x1 1,5 0,25 1,5 0,5 0,90 0,58
x2 -2,5 0,00 -2 -1,19 -1,84 -1,48
x3 3,5 0,75 1,12 0,81 1,12 0,87
x4 -3 -4,75 -3,37 -3,56 -3,40 -3,56
Erreurs 2,874 2,199 0,934 0,575 0,368 0,274

Table 1 – Itérations et erreurs pour la méthode de jacobi

- La méthode de Gauss-seidel

⃗ (0) = [0 0 0 0]T , et Erreur = ∥XLU − XGauss ∥ On trouve :

A partir de X

1 2 3 4 5 6
x1 1,5 0,62 1,06 0,74 0,68 0,65
x2 -1,75 -0,87 -1,51 -1,62 -1,69 -1,71
x3 2,62 0,90 1,01 0,93 0,92 0,91
x4 -4,31 -3,45 -3,50 -3,47 -3,46 -3,46
Erreurs 2,104 0,852 0,490 0,152 0,068 0,026

Table 2 – Itérations et erreurs pour La méthode de Gauss-seidel

5
 Nous pouvons remarquer que les erreurs diminuent au fil des itérations entre
les vrais valeurs et les itérations peut importe la méthode utilisé (Jacobi et
Gauss-seidel).