UIT. FSK.

Département de Mathématiques Filière SMI : S4

Exercice 93 page 246

On suppose que ab et bn sont les proportions de colza des espèces A et B du champs l’année n et
qu’elles vérifient 
an+1 = 0.7an − 0.2bn + 0.07
bn+1 = −0.4an + 0.7bn + 0.1
 
an
On note Un la matrice colonne
bn
Partie A: détermination de l’état stable

1. déterminer les matrices A et C telles que pour tout n ∈, Un+1 = AUn + C.

2. trouver l’état Stable S qui vérifie S = AS + C, après avoir justifié son existence

Partie B: Etude de la convergence de la suite (Un+1 )

1. déterminer les matrices A et C telles que pour tout

2. trouver l’état Stable S qui vérifie S = AS + C, après avoir justifié son existence

M. H. EL ALJ 1 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

Module Analyse Numérique

Série N° 1 : Systèmes linéaires

Exercice 1 décomposition LU par la méthode du pivot de Gauss

soit la matrice A définie par:  
1 1 1
A=  2 3 4 
2 5 7
1. Décrire la suite des opérations sur les lignes pour mettre en place la méthode du pivot sur A

2. En déduire la matrice U de la décomposition LU

3. Résoudre le système Ax = b avec b = (1.5, 3, 4)t .

Exercice 2 Donner la décomposition en LU des matrices suivantes

 
  2 1 3 1
1 2 3  −4 −1 −4 −3 
A =  3 2 1 , B =   0 −1 −3

2 
2 1 1
2 3 2 −1

Puis résoudre le système Ax = b avec b = (4, 8, 4)t et Bx = c avec c = (10, −13, −12, 9)t .

Exercice 3 soit la matrice A définie par:

 
2 −1 4 0
 4 −1 5 1 
A=
 −2

2 −2 3 
0 3 −9 4

1. Écrire la décomposition de la matrice LU de A

2. Resoudre les systèmes linéaires Aui = ei 1 ≤ i ≤ 4 oú ei est le vecteur de la base

canonique de R4 .

Exercice 4 pour n ∈ N on définit la matrice tridiagonale d’ordre n

 
2 −1
 −1 2 −1 
 
An = 
 .. .. .. 
. . . 
 
 −1 2 −1 
−1 2

1. Calculer la factorisation LU de A3
2. Calculer la factorisation LU de A4
3. Calculer la factorisation LU de An .

M. H. EL ALJ 2 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

Exercice 5 décomposition De Cholesky

1. Donner la décomposition De Cholesky RRT des matrices A et B suivantes
 
  4 0 12 −6
1 −2 0  0 1 2 1 
A =  −2 8 −6  ; B =   12 2 49 −4 

0 −6 25
−6 1 −4 51

2. Donner la décomposition en LU des matrices A et B.

3. Retrouver les matrices R à partir des décompositions LU précedentes

Exercice 6 on pose  
−2 1 −1 1
 2 0 4 1 
A=
 −4 −1 −12

9 
−2 1 1 −4
1. donner la décomposition LU de la matrice A
2. En déduire la solution du système Ax = b oú b = (1.5, 4, −14, −6.5)t .
3. soit B = U t ALt , sans calcul supplémentaire donner une décomposition LU de la matrice B

Exercice 7 1. réaliser la décomposition LU de la matrice

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

3 1 8 13

2. En déduire la solution du système Ax = b oú b = (0, 2, −1, 5)t .

3. sans calculer A2 résoudre le système A2 x = b

Exercice 8 pour n ∈ N on définit la matrice tridiagonale d’ordre n

 
2 −1
 −1 2 −1 
 
An = 
 . . . . . . 
. . . 
 
 −1 2 −1 
−1 2

1. Déterminer la factorisation Cholesky des matrices A3 et A4

2. Déterminer la factorisation Cholesky de la matrice An .

Exercice 9 Méthodes itératives

1. Reprendre par la méthode de Jacobi et de Gauss seidel l’exemple du cours avec comme vecteur
initial X 0 = (1, 1, 1)t
2. Résoudre le système A4 x = b de l’Exercice ?? o’u b = (3, −5, 7, −6)t par la méthode LU

3. Reprendre le 2. (A4 x = b) par la méthode de Jacobi et de Gauss seidel

M. H. EL ALJ 3 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

Module Analyse Numérique

Corrigé de la série N° 1 :
Solution de l’exercice 1 Méthode de Gauss
1. Opérations sur les lignes
 
1 1 1 (l1 )
pour la matrice A =  2 3 4  de lignes (l2 ) on effectue les transformations suivantes
2 5 7 (l3 )
 
pour la 1re étape, on effectue (l1 ) ← (l1 ) 1 1 1
(l2 ) ← (l2 ) − 2(l1 ) , ce qui donne A1 =  0 1 2 
(l3 ) ← (l3 ) − 2(l1 ) 0 3 5
 
pour la 2me étape, on effectue (l1 ) ← (l1 ) 1 1 1
(l2 ) ← (l2 ) , ce qui donne A2 =  0 1 2 =U
(l3 ) ← (l3 ) − 3(l2 ) 0 0 −1

2. Décomposition LU
    
1 0 0 1 1 1 1 1 1
on a L1 A =  −2 1 0   2 3 4  =  0 1 2  = A1
−2 0 1 2 5 7 0 3 5
    
1 0 0 1 1 1 1 1 1
(1)
et L2 A =  0 1 0   0 1 2  =  0 1 2  = A2 = U
0 −3 1 0 3 5 0 0 −1
et on en déduit que L, qui ne nécessite aucun calcul, est donnée par
    
1 0 0 1 0 0 1 0 0
−1 −1
L = L1 L2 =  2 1 0   0 1 0 = 2 1 0 
2 0 1 0 3 1 2 3 1

et on vérifie que
    
1 0 0 1 1 1 1 1 1
LU =  2 1 0   0 1 2 = 2 3 4 =A
2 3 1 0 0 −1 2 5 7

3. Résolution de Ax = b

• résolution de Ly = b
     
1 0 0 y1 1.5  y1 = 1.5
Ly =  2 1 0   y2  =  3  = b ⇒ 2y1 + y2 = 3
2 3 1 y3 4 2y1 + 3y2 + y3 = 4



 y1 = 1.5
ainsi, par la méthode
y2 = 3 − 2y1 = 0
de la descente, on a
y3 = 4 − 4y1 − 3y2 = 1


M. H. EL ALJ 4 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

• résolution de U x = y
     
1 1 1 x1 1.5  x1 + x2 + x3 = 1.5
Ux =  0 1 2   x2  =  0  =y ⇒ x2 + 2x3 = 0
0 0 −1 x3 1 −x3 = 1



 x3 = −1
ainsi, par la méthode
x2 = −2x3 = 2
de la remontée, on a
x1 = 1.5 − x2 − x3 = 0.5

    
1 1 1 0.5 1.5
on vérifie aisément que Ax =  2 3 4   2  =  3  = b
2 5 7 −1 4

Solution de l’exercice 2 Décomposition LU

• pour la matrice A
– Décomposition LU
    
1 0 0 1 2 3 1 2 3
on a L1 A =  −3 1 0   3 2 1 = 0 −4 −8  = A1
−2 0 1 2 1 1 0 −3 −5
    
1 0 0 1 2 3 1 2 3
1
et L2 A =  0 1 0   0 −4 −8  =  0 −4 −8  = A2 = U
3
0 −4 1 0 −3 −5 0 0 1
ainsi L est donnée par
    
1 0 0 1 0 0 1 0 0
L = L−1 −1
1 L2 =  3 1 0  0 1 0  =  3 1 0 
2 0 1 0 34 1 2 34 1

et on vérifie que
    
1 0 0 1 2 3 1 2 3
LU =  3 1 0   0 −4 −8  =  3 2 1  = A
2 34 1 0 0 1 2 1 1

– Résolution de Ax = b
∗ résolution de Ly = b
     
1 0 0 y1 4  y1 = 4
Ly =  3 1 0   y2  =  8  =b ⇒ 3y1 + y2 = 8
3 3
2 4 1 y3 4 2y1 + 4 y2 + y3 = 4



 y1 = 4
ainsi, par la méthode
y2 = 8 − 3y1 = −4
de la descente, on a
y3 = 4 − 2y1 − 43 y2 = −1


M. H. EL ALJ 5 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

∗ résolution de U x = y
     
1 2 3 x1 4  x1 + 2x2 + 3x3 = 4
Ux =  0 −4 −8   x2  =  −4  =y ⇒ −4x2 − 8x3 = −4
0 0 1 x3 1 x3 = −1



 x3 = −1
ainsi, par la méthode
x2 = − 14 (−4 + 8x3 ) = 3
de la remontée, on a
x1 = 4 − 2x2 − 3x3 = 1

    
1 2 3 1 4
on vérifie aisément que Ax =  3 2 1   3  =  8  = b
2 1 1 −1 4

• de même pour la matrice B

– Décomposition LU
    
1 0 0 0 2 1 3 1 2 1 3 1
 2 1 0 0   −4 −1 −4 −3
    0 1 2 −1 
 = B1
L1 B =  =
 0 0 1 0   0 −1 −3 2   0 −1 −3 2 
−1 0 0 1 2 3 2 −1 0 2 −1 −2
    
1 0 0 0 2 1 3 1 2 1 3 1
 0 1 0 0  0
  1 2 −1   0 1 2 −1 
L2 B 1 =  =  = B2
 0 1 1 0  0 −1 −3 2   0 0 −1 1 
0 −2 0 1 0 2 −1 −2 0 0 −5 0
    
1 0 0 0 2 1 3 1 2 1 3 1
 0 1 0 0  0
  1 2 −1   0
  1 2 −1 
L3 B 2 =  =  = B3 = U
 0 0 1 0  0 0 −1 1   0 0 1 0 
0 0 −5 1 0 0 −5 0 0 0 0 −5

et on en déduit que L = L−1 −1 −1

1 L2 L3 est donnée par
     
1 0 0 0 1 0 0 0 1 0 0 0 1 0 0 0
 −2 1 0 0   0 1 0 0   0 1 0 0   −2 1 0 0 
L=  0 0 1 0  0
  = 
1 1 0  0 0 1 0   0 1 1 0 
1 0 0 1 0 −2 0 1 0 0 5 1 1 −2 5 1

et on vérifie que
    
1 0 0 0 2 1 3 1 2 1 3 1
 −2 1 0 0  0 1
  2 −1   −4 −1 −4 −3
  
LU =  = =B
 0 1 1 0  0 0 1 0   0 −1 −3 2 
1 −2 5 1 0 0 0 −5 2 3 2 −1

M. H. EL ALJ 6 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

– Résolution de Bx = c avec c = (10, −13, −12, 9)t

∗ résolution de Ly = c
     
1 0 0 0 y1 10 
 y1 = 10
 −2 1 0 0   y2   −13 
 
−2y1 + y2 = −13
Ly = 
 0 −1 1 0  
 = =b ⇒
y3   −12  
 −y2 + y3 = −12
1 2 5 1 y4 9 y1 + 2y2 + 5y3 + y4 = 9




 y1 = 10
ainsi, par la méthode y2 = −13 + 2y1 = 7

de la descente, on a 
 y3 = −12 + y2 = −5
y4 = 9 − y1 − 2y2 − 5y3 = 9 − 10 − 14 + 25 = 10


∗ résolution de U x = y de même
     
2 1 3 1 x1 10 
 2x1 + x2 + 3x3 + x4 = 10
 0 1 2 −1   x2   7  x2 + 2x3 − x4 = 7

Ux =   = =y ⇒
 0 0 −1 1   x3   −5   −x3 + x4 = −5

0 0 0 −5 x4 10 −5x4 = 10




 x4 = −2
ainsi, par la méthode x3 = 5 + x4 = 3

de la remontée, on a 
 x2 = 7 − 2x3 + x4 = 7 − 6 + 2 = −1
x1 = 21 (10 − x2 − 3x3 − x4 ) = 12 (10 + 1 − 9 + 2) = 2

    
2 1 3 1 2 10
 −4 −1 −4 −3   −1   −13 
on vérifie aisément que Ax =  0 −1 −3
 = =b
2   3   −12 
2 3 2 −1 −2 9

Solution de l’exercice 3 Décomposition LU

1. Décomposition LU
    
1 0 0 0 2 −1 4 0 2 −1 4 0
 −2 1 0 0   4 −1 5 1  1 −3 1 
= 0

L1 A =    = A1
 1 0 1 0   −2 2 −2 3   0 1 2 3 
0 0 0 1 0 3 −9 4 0 3 −9 4
    
1 0 0 0 2 −1 4 0 2 −1 4 0
 0 1 0 0  0 1 −3 1 
= 0
 1 −3 1 
L2 A1 =  0 −1
  = A2 = A3 = U
1 0  0 1 2 3   0 0 5 2 
0 −3 0 1 0 3 −9 4 0 0 0 1

et on en déduit que L = L−1 −1

1 L2 est donnée par
    
1 0 0 0 1 0 0 0 1 0 0 0
 2 1 0 0  0 1 0 0   2 1
  0 0 
L=  −1 0 1 0
 = 
 0 1 1 0   −1 1 1 0 
0 0 0 1 0 3 0 1 0 3 0 1

M. H. EL ALJ 7 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

et on vérifie que
    
1 0 0 0 2 −1 4 0 2 −1 4 0
 2 1 0 0  0
  1 −3 1   4 −1 5 1 
LU =   −1 1
=
  −2
=A
1 0  0 0 5 2 2 −2 3 
0 3 0 1 0 0 0 1 0 3 −9 4

2. Résolution des systèmes Aui = ei pour i = 1, . . . , 4

nous ne traiterons ici en détail que la résolution du cas oú i = 2 c’est à dire Au2 = e2

• résolution de Ly = e2
    
1 0 0 0 y1 0 
 y1 = 0
 2 1 0 0   y2
    1  
2y1 + y2 = 1
Ly =   −1 1 1
= =b ⇒
0   y3   0  
 −y1 + y2 + y3 = 0
0 3 0 1 y4 0 3y2 + y4 = 0




 y1 = 0
ainsi, par la méthode y2 = 1 − 2y1 = 1

de la descente, on a  y3 =
 y1 − y2 = 0 − 1 = −1
y4 = −3y2 = −3


• résolution de U x = y de même
   
 
2 −1 4 0 x1 0 
 2x1 − x2 + 4x3 = 0
 0 1 −3 1   x2   1  x2 − 3x3 + x4 = 1
  
Ux =    −1  = y ⇒ 
= 
 0 0 5 2   x3  5x3 + 2x4 = −1
0 0 0 1 x4 −3 x4 = −3



 x4 = −3
= 15 (−1 − 2x4 ) = 1

 x3

ainsi, par la méthode
de la remontée, on a x2 = 1 + 3x3 − x4 = 1 + 3 + 3 = 7
3


 x1 = 12 (x2 − 4x3 ) = 21 (7 − 4) =

2
 3   
2 −1 0

4 0
 4 −1 5 1   2  
0 

 7  
on vérifie aisément que Ax =  = =b
 

 −2 2 −2 3  1     1 
0 3 −9 4 −3 0

on obtient ainsi le tableau récapitulatif suivant:

i=1 on a y = (1, −2, 3, 6)t et u1 = (−2.6, −13.4, −1.8, 6)t

i=2 on a y = (0, 1, −1, 3)t et u2 = (1.5, 7, 1, −3)t
i=3 on a y = (0, 0, 1, 0)t et u3 = (−0.1, 0.6.4, 0.2, 0)t
i=4 on a y = (0, 0, 0, 1)t et u4 = (−0.3, −2.2, −0.4, 1)t

M. H. EL ALJ 8 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

Solution de l’exercice 4 Décomposition LU

1. Décomposition LU de A3
    
1 0 0 2 −1 0 2 −1 0
on a L1 A3 =  1
2
1 0   −1 2 −1  =  0 3
2
−1  = A13
0 0 1 0 −1 2 0 −1 2
    
1 0 0 2 −1 0 2 −1 0
1
et L2 A3 =  0 1 0   0 3
2
−1  =  0 3
2
−1  = A23 = U
0 23 1 0 −1 2 0 0 4
3
et on en déduit que L, est donnée par
    
1 0 0 1 0 0 1 0 0
−1 −1
L = L1 L2 =  − 12 1 0  0 1 0  =  − 12 1 0 
2
0 0 1 0 −3 1 0 − 23 1
et on vérifie que
    
1 0 0 2 −1 0 2 −1 0
LU =  − 21 1 0  0 3
2
−1  =  −1 2 −1  = A3
2 4
0 −3 1 0 0 3
0 −1 2

2. Décomposition LU de A4
    
1 0 0 0 2 −1 0 0 2 −1 0 0
 1 1 0 0   −1 2 −1 0  3
−1
= 0 0 

L1 A4 =  2  2  = A14
 0 0 1 0   0 −1 2 −1   0 −1 2 −1 
0 0 0 1 0 0 −1 2 0 0 −1 2
    
1 0 0 0 2 −1 0 0 2 −1 0 0
3 3
 0 1 0 0  0 −1 0 = 0
 −1 0 
L2 A14 =  2

  0 −1
2 2
4
 = A24
 0
3
1 0 2 −1   0 0 3
−1 
0 0 0 1 0 0 −1 2 0 0 −1 2
    
1 0 0 0 2 −1 0 0 2 −1 0 0
3 3
 0 1 0 0  0
  −1 0   0
  −1 0 
L3 B 2 =  2
4 = 2
4
 = A34 = U
 0 0 1 0  0 0 3
−1   0 0 3
−1 
3 5
0 0 4 1 0 0 −1 2 0 0 0 4

et on en déduit que L = L−1 −1 −1

1 L2 L3 est donnée par
     
1 0 0 0 1 0 0 0 1 0 0 0 1 0 0 0
 −1 1 0 0   0 1 0 0   0 1 0 0   −1 1 0 0 
L= 2   = 2 
 0 0 1 0   0 −2 1 0  0 0 1 0   2
0 −3 1 0 
3
0 0 0 1 0 0 0 1 0 0 − 34 1 0 0 − 43 1
et on vérifie que
    
1 0 0 0 2 −1 0 0 2 −1 0 0
 −1 1 0 0  0
  3
−1 0   −1 2 −1 0 
LU =  2 2 =  = A4
 0 −2 1 0   0 0 4
−1   0 −1 2 −1 
3 3
0 0 − 43 1 0 0 0 5
4
0 0 −1 2

M. H. EL ALJ 9 Mars 2022

 UIT. FSK. Département de Mathématiques Filière SMI : S4

3. Décomposition LU de An par extrapolation par rapport aux résultats obtenus pour n = 3 et

n = 4 les matrices Ln et Un sont définies comme suit:

1 0 ... ... ... 0 2 −1 0 ... ... ... 0

   
 −1 .. 3 ..
 2 1 . 

 0
 2
−1 .


 0 ..
 . ... ..  
 0 ... ... ..

 . 
 

.
Ln =  ... ..  . ..
− i−1  et Un =  .. i+1
  
 . i
1 . i
−1  .
 . ... ... ..   .
 . ... ...
 ..
 . .  . .
 
 
 .  .
 ..  ..
 
− n−2
n−1
1 0  n
n−1
−1 
n−1 n+1
0 ... ... ... 0 − n 1 0 ... ... ... ... 0 n

En effet si l’on note An = (aij )i,j=1,...,n oú An = Ln Un , on vérifie aisément

i−2


 ai−1i = − × 0 + 1 × (−1) = −1

 i−1
i−1 i+1

aii = − × (−1) + 1 × =2
 i i
i−1


 aii+1 = − × 0 + 1 × (−1) = −1 ∀i = 2, . . . , n − 1

i

Solution de l’exercice 5 Décomposition de Cholesky

1. Résolution directe On rappelle que les formules directes de la méthode de Cholesky sont
données par :
√
r11 = a11 (1)
al1
rl1 = pour l = 2, ..., n (2)
a
v11 !
u
u Xi−1
rii = t aii − 2
rik pour i = 2, ..., n (3)
k=1
i−1
!
1 X
rli = ali − rlk rik pour i = 2, ..., n − 1 et l > i (4)
rii k=1

on calcule ainsi tous les termes de la matrice R colonne par colonne, la première colonne en
commençant par r11 , puis la seconde et ainsi de suite en déterminant tout d’abord rii pour
la ième colonne.

• Matrice A Le produit matriciel A = RRt s’écrit :

    
1 −2 0 r11 0 0 r11 r21 r31
A =  −2 8 −6  = RRt =  r21 r22 0   0 r22 r32 
0 −6 25 r31 r32 r33 0 0 r33

M. H. EL ALJ 10 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

Les Formules établies précédemment entrainent :

Première Colonne
√ √
- r11 = a11 = 1 = 1 d’après (??).

a21 −2
- r21 = = = −2 d’après (??).
r11 1
a31 0
- r31 = = =0 d’après (??).
r11 1

Deuxième Colonne
p
2
p √
- r22 = (a22 − r21 ) = 8 − (−2)2 = 4 = 2 d’après (??).

a32 − r31 r21 1

- r32 = = (−6 − 0 × (−2)) = −3 d’après (??).
r22 2

Troisième Colonne
p
2 2
p √
- r33 = (a33 − r31 − r32 ) = 25 − (O)2 − (−3)2 = 16 = 4 d’après (??).

et on vérifie que
    
1 0 0 1 −2 0 1 −2 0
RRt =  −2 2 0  0 2 −3  =  −2 8 −6  = A
0 −3 4 0 0 4 0 −6 25

• Matrice B Le produit matriciel B = RRt s’écrit :

    
4 0 12 −6 r11 0 0 0 r11 r21 r31 r41
 0 1 2 1   r21 r22 0 0   0 r22 r32 r42 
 t
  
B= 12 2 49 −4  = RR =  r31 r32 r33

0  0 0 r33 r43 
−6 1 −4 51 r41 r42 r43 r44 0 0 0 r44

Les Formules établies précédemment entrainent:

Première Colonne
√ √
- r11 = a11 = 4 = 2 d’après (??).
a21 0
- r21 = = =0 d’après (??).
r11 2
a31 12
- r31 = = =6 d’après (??).
r11 2
a41 −6
- r41 = = = −3 d’après (??).
r11 2
Deuxième Colonne
p p
2
- r22 = (a22 − r21 ) = 1 − (0)2 = 1 d’après (??).
a32 − r31 r21 (2 − 6 × 0)
- r32 = = =2 d’après (??).
r22 1
a42 − r41 r21 (1 − (−3) × 0)
- r42 = = =1 d’après (??).
r22 1

M. H. EL ALJ 11 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

Troisième Colonne
p p
2 2
- r33 = (a33 − r31 − r32 ) = 49 − (6)2 − (2)2 = 3 d’après (??).

a43 − r41 r31 − r42 r32 (−4 − (−3) × 6 − 1 × 2)

- r43 = = =4 d’après (??).
r33 3

Quatrième Colonne
p p
2 2 2
- r44 = (a44 − r41 − r42 − r43 ) = 51 − (−3)2 − (1)2 − (4)2 = 5 d’après (??).

et on vérifie que
    
2 0 0 0 2 0 6 −3 4 0 12 −6
 0 1 0 0  0 1
  2 1   0
  1 2 1 
RRt =  = =B
 6 2 3 0  0 0 3 4   12 2 49 −4 
−3 1 4 5 0 0 0 5 −6 1 −4 51

2. Décomposition LU

• Matrice A
    
1 0 0 1 −2 0 1 −2 0
on a L1 A =  2 1 0   −2 8 −6  =  0 4 −6  = A1
0 0 1 0 −6 25 0 −6 25
    
1 0 0 1 −2 0 1 −2 0
1
et L2 A =  0 1 0   0 4 −6  =  0 4 −6  = A2 = U
0 32 1 0 −6 25 0 0 16
et on vérifie que
    
1 0 0 1 −2 0 1 −2 0
LU =  −2 1 0  0 4 −6  =  −2 8 −6  = A
3
0 −2 1 0 0 16 0 −6 25

• Matrice B
    
1 0 0 0 4 0 12 −6 4 0 12 −6
 0 1 0 0  0 1 2 1 
  0 1 2 1 
 1
L1 B = 
 −3 0 1 0   12 2 49 −4  =  0 2 13 14  = B
 
3
2
0 0 1 −6 1 −4 51 0 1 14 42
    
1 0 0 0 4 0 12 −6 4 0 12 −6
 0 1 0 0   0 1 2 = 0 1 2
1   1 
L2 B 1 = 
 0 −2 1
  = B2
0  0 2 13 14   0 0 9 12 
0 −1 0 1 0 1 14 42 0 0 12 41
    
1 0 0 0 4 0 12 −6 4 0 12 −6
 0 1 0 0  0 1 2 1 
= 0 1 2
 1 
L3 B 2 =  0 0
  = B3 = U
1 0  0 0 9 12   0 0 9 12 
0 0 − 43 1 0 0 12 41 0 0 0 25

M. H. EL ALJ 12 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

et on en déduit que L = L−1 −1 −1

1 L2 L3 est donnée par
     
1 0 0 0 1 0 0 0 1 0 0 0 1 0 0 0
 0 1 0 0  0 1 0
  0  0 1 0
  0   0
  1 0 0 
L=  3 0 = 
1 0  0 2 1 0  0 0 1 0   3 2 1 0 
− 23 0 0 1 0 1 0 1 0 0 43 1 − 32 1 34 1

et on vérifie que
    
1 0 0 0 4 0 12 −6 4 0 12 −6
 0 1 0 0  0 1 2 1 
  0 1 2 1 

LU =   3 2 1 0   0 0 9 12  =  12 2 49 −4  = B
 

− 32 1 43 1 0 0 0 25 −6 1 −4 51

3. Résolution à partir de la décomposition LU

• Matrice A
   
1 0 0 1 −2 0
√
Oon note D = diag ( uii )i=1,...,3 = 0 2 0  avec U =  0 4 −6 
0 0 4 0 0 16
    
1 0 0 1 0 0 1 0 0
Alors LD =  −2 1 0   0 2 0  =  −2 2 0 =R
3
0 −2 1 0 0 4 0 −3 4
    
1 0 0 1 −2 0 1 −2 0
−1
et D U =  1
0 2 0   0 4 −6  =  0 2 −3  = Rt
0 0 41 0 0 16 0 0 4
• Matrice B
   
2 0 0 0 4 0 12 −6
√  0 1 0 0   0 1 2 1 
Oon note D = diag ( uii )i=1,...,3 =
 0
 avec U =  
0 3 0   0 0 9 12 
0 0 0 5 0 0 0 25
    
1 0 0 0 2 0 0 0 2 0 0 0
 0 1 0 0  0 1 0 0   0
  1 0 0 
Alors LD =  
 3 2 1 0  0 0 = =R
3 0   6 2 3 0 
− 23 1 43 1 0 0 0 5 −3 1 4 5

1
    
0 0 0
2
4 0 12 −6 2 0 6 −3
0 1 0 0 0 1 2 1   0 1 2 1 
et D−1 U = 
    
 = Rt
 0 0 1 0   0 0 9 12  =  0 0

3
3 4 
0 0 0 51 0 0 0 25 0 0 0 5

M. H. EL ALJ 13 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

Solution de l’exercice 6 Décomposition LU

1. Décomposition LU de la matrice A
    
1 0 0 0 −2 1 −1 1 −2 1 −1 1
 1 1 0 0  2 0 4 1   0
  1 3 2 
L1 A =   =  = A1
 −2 0 1 0   −4 −1 −12 9   0 −3 −10 7 
−1 0 0 1 −2 1 1 −4 0 0 2 −5
    
1 0 0 0 −2 1 −1 1 −2 1 −1 1
 0 1 0 0   0 1 3 2 = 0
  1 3 2 
L2 A 1 =  
 0 3 1 0   0 −3 −10
 = A2
7   0 0 −1 13 
0 0 0 1 0 0 2 −5 0 0 2 −5
    
1 0 0 0 −2 1 −1 1 −2 1 −1 1
 0 1 0 0  0 1
  3 2   0 1
  3 2 
L3 A2 =  0 0 1 0   0 0 −1 13  =  0 0
 = A3 = U
−1 13 
0 0 2 1 0 0 2 −5 0 0 0 21

et on en déduit que L = L−1 −1 −1

1 L2 L3 est donnée par
     
1 0 0 0 1 0 0 0 1 0 0 0 1 0 0 0
 −1 1 0 0  0
  1 0 0  0
  1 0 0   −1
  1 0 0 
L= = 
 2 0 1 0   0 −3 1 0  0 0 1 0   2 −3 1 0 
1 0 0 1 0 0 0 1 0 0 −2 1 1 0 −2 1

et on vérifie que
    
1 0 0 0 −2 1 −1 1 −2 1 −1 1
 −1 1 0 0  0 1
  3 2   2 0 4 1 
LU =   2 −3
= =A
1 0   0 0 −1 13   −4 −1 −12 9 
1 0 −2 1 0 0 0 21 −2 1 1 −4

2. Résolution du système Ax = b avec b = (1.5, 4, −16, −6.5)

• résolution de Ly = b
  
y
  
1 0 0 0 1 0 
 y1 = 1.5
 −1  y  
1 0 0  0  −y 1 + y2 = 4

2  =b ⇒
Ly =  2 −3
 =
1 0  y3   1 2y1 − 3y2 + y3 = −14
   

1 0 −2 1 0 y1 − 2y3 + y4 = −6.5

y4


 y1 = 1.5
ainsi, par la méthode y2 = 4 + y1 = 5.5

de la descente, on a 
 y3 = −14 − 2y1 + 3y2 = −14 − 3 + 16.5 = −0.5
y4 = −6.5 − y1 + 2y3 = −6.5 − 1.5 − 1 = −9


• résolution de U x = y
     
−2 1 −1 1 x1 1.5 
 −2x1 + x2 − x3 + x4 = 1.5
 0 1 3 2   x2   5.5  
x2 + 3x3 + 2x4 = 5.5
Ux =   = =y ⇒
 0 0 −1 13   x3   −0.5   −x3 + 13x4 = −0.5

0 0 0 21 x4 −9 21x4 = −9


M. H. EL ALJ 14 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

3

9

 x4 = − 21 =−

 7

 1 39 71
 x3 = (0.5 + 13x4 ) = − =−

ainsi, par la méthode 2 7
14
de la remontée, on a 213 6 151

 x2 = 5.5 − 3x3 − 2x4 = 5.5 + 14 + 7
=
7


173


 x1 = − 1 (1.5 − x2 + x3 − x4 ) = − 1 (1.5 − 151 −
 71
) =
2 2 7 14
14

− 173
    
−2 1 −1 1 14
1.5
151
 2 0 4 1  − 7   4 
on vérifie aisément que Ax =  = =b
    
71
−4 −1 −12 − 14 −14

 9    
−2 1 1 −4 − 37 −6.5

on obtient ainsi le tableau récapitulatif suivant:

3. Décomposition LU de la matrice B = U t ALt

′ ′
D’après 1. A = LU donc B = U t LU Lt = (U t LD−1 ) (DU Lt ) = L U
   1 
−2 0 0 0 −2 0 0 0
 0 1 0 0  −1  0 1 0 0 
 
oú D = diag(uii )i=1,...,4 = 
 0 0 −1 0  et D =  0

0 −1 0 
1
0 0 0 21 0 0 0 21

′
Alors L = U t LD−1 est donnée par
   1   
2 0 0 0 1 0 0 0 −2 0 0 0 1 0 0 0
′  1 1 0 0   −2
  1 0 0  0 1
  0 0   0 1 0 0 
L =
 −1
= 
3 −1 0   2 −3 1 0   0 0 −1 0   3 6 1 0 
1
1 2 13 21 1 0 2 1 0 0 0 21
−23 −37 29 1
′
de même U = DU Lt est donné par:
     
−2 0 0 0 −2 1 −1 1 1 −2 2 1 4 −6 16 −2
′  0 1 0 0  0 1
  3 2  0 1 −3 0   0
  1 0 −4 
U = 0 0 −1 0   0 0 −1 13
 = 
 0 0 1 2   0 0 1 −15 
0 0 0 21 0 0 0 21 0 0 0 1 0 0 0 441

En effet, on a B = U t LU Lt soit :
     
1 0 0 0 −2 1 −1 1 1 −2 2 1 4 −6 16 −4
 −2 1 0 0  2
  0 4 1  0 1 −3 0  
  0 1 0 −3 
B=
 2 −3 1
 = 
0   −4 −1 −12 9  0 0 1 2   12 −12 49 −45 
1 0 2 1 −2 1 1 −4 0 0 0 1 −92 101 −339 209

et on vérifie que
    
1 0 0 0 4 −6 16 −2 4 −6 16 −4
′ ′  0 1 0 0  0 1 0 −4   0 1 0 −3 
LU =  = =B
 3 6 1 0  0 0 1 −15   12 −12 49 −45 
−23 −37 29 −1 0 0 0 441 −92 101 −339 209

M. H. EL ALJ 15 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

Solution de l’exercice 7 Décomposition LU

1. Décomposition LU
    
1 0 0 0 −1 1 −3 0 −1 1 −3 0
 1 1 0 0   1 1 3 8   0 2 0 8 
L1 A =   =  = A1
 −2 0 1 0   −2 2 −5 −1   0 0 1 1 
3 0 0 1 3 1 8 13 0 4 −1 13
    
1 0 0 0 −1 1 −3 0 −1 1 −3 0
 0 1 0 0  0 2
  0 8   0
  2 0 8 
L2 A1 = 
 0 =  = A2
0 1 0  0 0 1 1   0 0 1 1 
0 −2 0 1 0 4 −1 13 0 0 −1 −3
    
1 0 0 0 −1 1 −3 0 −1 1 −3 0
 0 1 0 0  0 2
  0 8   0 2
  0 8 
L3 A2 = 
 0 0 1 =  = A3 = U
0  0 0 1 1   0 0 1 1 
0 0 1 1 0 0 −1 −3 0 0 0 −2

et on en déduit que L = L−1 −1 −1

1 L2 L3 est donnée par
     
1 0 0 0 1 0 0 0 1 0 0 0 1 0 0 0
 −1 1 0 0   0 1 0 0   0 1 0 0   −1 1 0 0 
L=  2 0 1
  = 
0  0 0 1 0  0 0 1 0   2 0 1 0 
−3 0 0 1 0 2 0 1 0 0 −1 1 −3 2 −1 1

et on vérifie que
    
1 0 0 0 −1 1 −3 0 −1 1 −3 0
 −1 1 0 0   0 2 0 8   1 1 3 8 
LU =   =
  −2
=A
 2 0 1 0  0 0 1 1 2 −5 −1 
−3 2 −1 1 0 0 0 −2 3 1 8 13

2. Résolution de Ax = b avec b = (0, 2, −1, 5)t

• résolution de Ly = b
   
y 0
 
1 0 0 0 1 
 y1 = 0
 −1 1 0 0 
 y2   2
   
−y1 + y2 = 2
Ly =   = =b ⇒

 2 0 1 0  y3   −1 2y1 + y3 = −1
 
 

−3 2 −1 1 −3y1 + 2y2 − y3 + y4 = 5

y4 5


 y1 = 0
ainsi, par la méthode y2 = 2 + y1 = 2

de la descente, on a 
 y3 = −1 − 2y1 = −1 − 2 × 0 = −1
y4 = 5 + 3y1 − 2y2 + y3 = −5 + 3 × 0 − 2 × 2 − 1 = 0


• résolution de U x = y de même
     
−1 1 3 0 x1 0 
 −x1 + x2 + 3x3 = 0
 0 2 0 8   x2   2
     
2x2 + 8x4 = 2
Ux =  = =y ⇒
 0 0 1 1   x3   −1  
 x3 + x4 = −1
0 0 0 −2 x4 0 −2x4 = 0


M. H. EL ALJ 16 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4



 x4 = 0
ainsi, par la méthode x3 = (−1 − x4 ) = −1

de la remontée, on a 
 x2 = 12 (2 − 8x4 ) = 21 (2 − 8 × 0) = 1
x1 = x2 − 3x3 = 1 − 3 × −1 = 4


    
−1 1 −3 0 4 0
 1 1 3 8  1
    2 
on vérifie aisément que Ax = 
 −2
=
  −1  = b

2 −5 −1   −1
3 1 8 13 0 5

3. Résolution de A2 x = b avec b = (0, 2, −1, 5)t

′
La résolution du système A2 x = b est équivalente à la résolution du système Ax = b oú
′
b = (4, 1, −1, 0)t est la solution du 3.
′
• résolution de Ly =b
     
1 0 0 0 y1 4 
 y1 = 4
 −1 1 0 0   y2   1  
−y 1 + y2 = 1
Ly =   = =b ⇒
 2 0 1 0 
 y3   −1  
 2y1 + y3 = −1
−3 2 −1 1 y4 0 −3y1 + 2y2 − y3 + y4 = 0




 y1 = 4
ainsi, par la méthode y2 = 1 + y1 = 1 + 4 = 5

de la descente, on a 
 y3 = −1 − 2y1 = −1 − 2 × 4 = −9
y4 = 0 + 3y1 − 2y2 + y3 = −5 + 3 × 4 − 2 × 5 − 9 = −7


• résolution de U x = y de même
     
−1 1 3 0 x1 4 
 −x1 + x2 + 3x3 = 4
 0 2 0 8 
  x2  =  5  2x2 + 8x4 = 5
    
Ux =  =y ⇒
 0 0 1 1   x3   −9  
 x3 + x4 = −9
0 0 0 −2 x4 −7 −2x4 = −7


−7


 x4 = −2
= 3.5
ainsi, par la méthode x3 = (−9 − x4 ) = −9 − 3.5 = −12.5

1
de la remontée, on a x = (5 − 8x4 ) = 21 (5 − 8 × 72 0) = −11.5
 2 2


x1 = −4 + x2 − 3x3 = −4 − 11.5 − 3 × −12.5 = 22

La solution est donc donnée par x = (22, −11.5, −12.5, 3.5)t , en effet:
    
−1 1 −3 0 −1 1 −3 0 8 −6 21 5
 1 1 3 8  1 1
  3 8   18 16
  49 115 
A2 = AA =   −2 2 −5 −1   −2 2 −5 −1  =  17 −9

45 24 
3 1 8 13 3 1 8 13 21 33 58 185
    
8 −6 21 5 22 0
 18 16 49 115   −11.5
    2 
on vérifie aisément que A2 x = 
 17 −9 45 24   −12.5
=
  −1  = b


21 33 58 185 3.5 5

M. H. EL ALJ 17 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

Solution de l’exercice 8 Factorisation de Cholesky A = RRt

1. Factorisation de cholesky de A3

• Méthode directe Première Colonne

√ √
- r11 = a11 = 2 d’après (??).

a21 −1
- r21 = =√ d’après (??).
r11 2
a31 0
- r31 = = √ =0 d’après (??).
r11 2

Deuxième Colonne
p q q
−1 2
- r22 = (a22 − r21 ) = 2 − ( 2 ) = 32
2 √ d’après (??).

a32 − r31 r21 q 2  −1

 q
2
- r32 = = 3 −1 − 0 × ( √ 2
) = − 3
d’après (??).
r22

Troisième Colonne
s r 2
p
2 2 2 2
- r33 = (a33 − r31 − r32 )= 2− (0)2 − (− ) =√ d’après (??).
3 3

et on vérifie que
 √  √
2 q0 0 2 − √12 0
 
2 −1 0

1 3
q q
− √2 0 
 
R3 R3t =  0 3
− 23  =  −1 2 −1 
 = A3
 
 q2  2
0 − 23 √2 0 0 √2 0 −1 2
3 3

• à partir de la décomposition LU
D’après L’exercice 4, on a:
    
1 0 0 2 −1 0 2 −1 0
L3 U3 =  − 12 1 0  0 3
2
−1  =  −1 2 −1  = A3
2 4
0 −3 1 0 0 3
0 −1 2
 √ √1 0 0
  
2 q0 0 2 q
Oon note D3 =  0 3
0  d’oú D3−1 =  0 2
0 
   
2 3 √
0 0 √2 0 0 23
3

  √2   √2 0 0


1 0 0 q 0 0 q
  − √1 3
0  = R3
1 3

Alors L3 D3 = −2 1 0  0 0 =

2
   2
2
0 − 23 1
q 
2 2 2
0 0 √
3 0 − 3 √3
 √1
0 0  2 −1
   √2 − √1 0

2 q 0 q2 q
et D3−1 U3 =  0 2
0 0 3
−1 =  0 3
− 23  = R3t
    
3  2 2
√ 4
0 0 23 0 0 3 0 0 √2
3

M. H. EL ALJ 18 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

2. Factorisation de cholesky de A4

• Méthode directe Première Colonne

√ √
- r11 = a11 = 2 d’après (??).
a21 −1
- r21 = =√ d’après (??).
r11 2
a31 0
- r31 = = √ =0 d’après (??).
r11 2
a41 0
- r41 = = √ =0 d’après (??).
r11 2
Deuxième Colonne
p q q
2 −1 2 3
- r22 = (a22 − r21 ) = 2 − (√ 2
) = 2
d’après (??).
a32 − r31 r21 q 
−1
 q
- r32 = = 3 −1 − 0 × 2 = − 23
2 √ d’après (??).
r22
a42 − r41 r21 q 2  −1

- r42 = = 3 0−0× √ 2
=0 d’après (??).
r22
Troisième Colonne
r q
p
- r33 = (a33 − r31 − r32 ) = 2 − (0)2 − (− 23 )2 =
2 2 √2
3
d’après (??).

a43 − r41 r31 − r42 r32 √  q  √

3 2 3
- r43 = = 2
−1 − 0 × 0 − 0 × − 3
=− 2
d’après (??).
r33

Quatrième Colonne
p q √ √
2
- r44 = (a44 − r41 2
− r42 2
− r43 ) = 2 − (0)2 − (0)2 − (− 23 )2 = 2
5
d’après (??).

et on vérifie que
 √  √
2 − √12

2 q0 0 0 0 0  
 √1 q q 2 −1 0 0
3
 − 2 0 0  3
− 23
 
q2 0 0   −1 2 −1 0 
R4 R4t =   2 √  =   = A4

0 − 23 √2 0 

0 0 √2 3
−2 
  0 −1 2 −1 
3 3

√ √ √ 0 0 −1 2
3 5 5
0 0 − 2 2
0 0 0 2

• à partir de la décomposition LU
D’après l’exercice 4, on a
    
1 0 0 0 2 −1 0 0 2 −1 0 0
 −1 3
1 0 0   =  −1 2 −1
 0 2 0   0 
L4 U4 =  2   2  = A4
 0 −2 1 0   0 0 4
−1   0 −1 2 −1 
3 3
3 5
0 0 −4 1 0 0 0 4
0 0 −1 2
 √   1 
2 q0 0 0 √ 0 0 0
 3
  2 q2 
 0 0 0  −1
 0 0 0 
Oon note D4 =  2 et D = 3
 4
 √ 
 0 2
0 √
3
0 
 0

0 2 3
0 

√
 
0 0 0 25 0 0 0 √25

M. H. EL ALJ 19 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

 √   √ 
  2 q0 0 0 2 q0 0 0
1 0 0 0    √1 3
3   − 2 0 0 

 −1 1 0 0  0 0 0 q2

Alors L4 D4 =  2  2 =  = R4
 0 −2 1 0  √2
0 0 0 0 − 23 √2 0 
  
3 3 √ √3
0 − 34
  
0 1 5 √
0 0 0 2 0 0 − 23 5
2
 √
2 − √12
  
√1 0 0 0   0 0
2 q  2 −1 0 0 q q
3
− 23
 2 3
 
0 0 0   0 0
−1  0 2 0 
  
et D4 U4 =  3 √ 2
4
 =  2 √  = R4t

0 0 23 0  0
  0 3 −1 
 0

0 √2 − 23

5 3
 
0 0 0 √
0 0 0 √25 4 0 0 0 5
2

3. Factorisation de cholesky de An

• Méthode directe  q
i+1
 rii = i = 1, . . . , n

qi

Pour n fixé, soit P(i) ri+1i = − i+1 i
i = 1, . . . , n − 1


rli = 0 i = 1, . . . , n et l = i + 1, . . . , n

Montrons que P(i) est vraie pour tout i ≤ n
i La propriété est vraie pour i = 1 car
√ √ q
∗ r11 = a11 = 2 = 21
q
∗ r21 = r22 = − 2 = − 12
a21 √1

∗ rl1 = ra22
l1
= − √02 = 0
ii On suppose que la propriété est vraie jusqu’à i − 1 alors
v
u i−1 q q
∗ rii =
u X
2
p
t aii − rik 2
= aii − rii−1 = 2 − i−1i
= i+1
i
k=1
i−1
!
X q q
∗ ri+1i = 1
rii
ai+1i − ri+1k rik = i
i+1
(−1 i
− 0 × ...) = − i+1
k=1
i−1
!
X q
∗ rli = 1
rii
ali − rlk rik = i
i+1
(0 − 0 × ...) = 0
k=1

la propriété est donc vraie pour tout i = 1, . . . , n, et la matrice Rn est donnée par:
 √ 
2 q0 ... ... ... ... 0
.. 
 − √12 3

2
. 
 
 ... ... .. 
 0 . 
. .. 
 q q 
 .
. − i−1 i+1
.
i
q i q
 
Rn =  .. .
 
. − i+1i i+2 .
. 

i+1

 
 .. .. .. .. 
 . . . . 
..
 q 
n−2 n
p
. − n−1 0 
 

 n−1 
q q
0 ... ... ... ... 0 − n−1 n
n+1
n

M. H. EL ALJ 20 Mars 2022

 UIT. FSK. Département de Mathématiques Filière SMI : S4

• à partir de la décomposition LU
D’après la décomposition An = Ln Un de l’Exercice 4. on pose
 √   1

2 q0 . . . ... ... ... 0 sqrt2
0 ... ... ... ... 0
3 ..  q .. 
0 .  2
 
 2  0 3
. 
... ... .. 
   

0 . 
 ... ... .. 
  0 . 
 .. ..
q
i+1 ..   
.. ..
q .. 

Dn = 
 . . i
.  et Dn−1 =  . . i
. 
i+1

 .. .. .. ..  
.. .. .. .. 


 . . . . 
 
 . . . . 
 .. p n  
.. q 
 . n−1
0  
. n−1
0 

 q   n
n+1 p n
0 ... ... ... ... 0 n 0 ... ... ... ... 0 n+1

alors la matrice Rn = Ln Dn est donnée par

 √ 
2 q0 . . . ... ... ... 0
1 0 ... ... ... ... 0
 
3 ..
1 ..   0 .

 −2 1 .   2 
.. .. ..
 
. .. ..  . .
0 ..  0 .
  
 . .   .


.. . .
q .. 
Ln Dn = 

. − i−1 1 ..   ..
  .. i+1
i
. 

i  .
 .. .. .. . . . .. 
. . . ..    .. .. .. .
  
 
..  .
.
 
n
− n−2
p
 . n−1
1 0   . n−1 q
0 

n−1
0 ... ... ... 0 − n 1 0 ... ... ... ... 0 n+1
n

de même la matrice Rnt = Dn−1 Un est donnée par:

 1 
sqrt2 q
0 ... ... ... ... 0  2 −1 0 ... ... ... 0

 2 ..  ..
 0 3
.   0
 3
−1 .
 2 
 . .
.. .. .
..   0
  ... ... ..
 0  .
.. ..   ..
  
... ..
q
Dn−1 Un =  . i
.  . i+1
 
i+1  . i
−1  .

.. . . .  . ... ...
 ..

 . . . . . .
.  .

 .
 .
..   ..
 q 
n

. n−1
0 n−1
−1 
 n 
n+1
0 ... ... ... ... 0
p n 0 ... ... ... ... 0 n
n+1

M. H. EL ALJ 21 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

Rappels de Cours
 
a11 a12 ... ... a1n
. .. ..
a21 . . . .
 
 
 .. . . ... .. 
A= . . aii . =D−E−F
.. .. ..
 
. . .
 
 an−1n 
an1 . . . . . . ann−1 ann

     
a11 0 ... ... 0 0 ... ... ... 0 0 −a12 . . . . . . −a1n
.. .. .. .. .. .. . . .. ..
0 . . .   −a21 . . . . . .
     
   
.. .. .. .   . ... .. .. .. ... ... ..
. .. − ..
   
= . . aii . . − . . 
.. .. ..   . .. .. .. .. ..
     
. . . 0   .. . . . . . −an−1n
   
   
0 . . . . . . 0 ann −an1 . . . . . . −ann−1 0 0 ... ... ... 0

où D : diagonale de A
E : triangulaire inférieure avec des 0 sur la diagonale
F : triangulaire supérieure avec des 0 sur la diagonale

. Méthode de Jacobi
La matrice d’itération de la méthode de Jacobi associée à Ax = b est donnée par
MJ = D−1 (E + F ) = I − D−1 A

Et L’itération de la méthode de Jacobi s’écrit:

xk+1 = MJ xk + bJ avec bJ = D−1 b

Cette itération est déterminée pratiquement par la formule

Dxk+1 = (E + F )xk + b

dont l’écriture matricielle est:

  0 −a12 ... ... −a1n

a11 xk+1 xk1
   
1 b1
..   −a21 . . . ..  .   .. 
. .   ..   . 
 
.. ... ... ..
    k  
k+1 k+1 
Dx =  aii xi = . .   xi  +  bi 
 
..   .. .. .. ..   ..  .. 
     
 . . . . .  .   . 
ann xk+1n −an1 . . . . . . −ann−1 −an−1n xkn bn

et pour une donnée x0 choisie on calcule xk+1 par

n
!
1 X
xk+1
i = bi − aij xkj , i = 1, . . . , n (5)
aii j=1,j̸=i

L’algorithme de Jacobi nécessite le stockage des deux vecteurs xki et xk+1

i

M. H. EL ALJ 22 Mars 2022

 UIT. FSK. Département de Mathématiques Filière SMI : S4

. Méthode de Gauss-Seidel
L’itération de la méthode de Gauss-Seidel s’écrit:
MG = (D − E)−1 F

x k+1
= MG x + bGk
avec
bG = (D − E)−1 b

Cette itération est déterminée pratiquement par la formule

Dxk+1 = Exk+1 + F xk + b

dont l’écriture matricielle est:

 
a11 0 ... ... 0 xk+1
1
 
a11 xk+1 1

.. ..  .   .
0 . .   ..   ..
 

.. ..
   k+1   
Dxk+1 = . .  x  =  aii xk+1

aii  i i 
.. .. ..   ..  ..
  
. . . . .
   
 
k+1
0 ... ... 0 ann xn ann xk+1 n

et Exk+1 + F xk + b
    
0 0 ... ... 0 xk+1
1
0 −a12 . . . ... −a1n xk1
 
b1

.. .. .   .. ..  .   .. 
 −a21 . .  ..   0 . .   ..   . 

 . ... ... .. .. . . . ..
 .. .
  k+1    k  
.
 x + . . . .   xi + bi 

 i
 . ..
..   .. 
.. .. .. ..   ..  .. 
  
 ..

.. . . . . . .   . 

−an−1n
   
−an1 . . . . . . −ann−1 0 x k+1
n 0 . . . . . . −ann−1 0 xkn bn

Ainsi, pour une valeur initiale x0 donnée on détermine xk+1 par les formules
i−1 n
!
1 X X
xk+1
i = bi − aij xk+1
j − aij xkj , i = 1, . . . , n (6)
aii j=1 j=i+1


1

 xk+1
1 = ( b1 −a12 xk2 ... ... −a1i xki ... ... −ain xkn )
a11



 1
xk+1 ( b2 −a21 xk+1

= −a23 xk3 ... −a2i xki ... ... −ain xkn )

 2 1
a22






..


.






1

 xk+1
i = ( bi −ai1 xk+1
1 ... ... −aii−1 xk+1
i−1 −aii+1 xki+1 ... −ain xkn )
aii






..


.








1


 xk+1 = ( bn −an1 xk+1 ... ... ... −ani xk+1 ... −ann−1 xk+1 )


n 1 i n
ann

L’algorithme de Gauss-Seidel ne nécessite qu’un vecteur de stockage xki et il est en général plus rapide.

M. H. EL ALJ 23 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

. Application Reprise de l’exemple du cours, soit le système Ax = b suivant:

    
3 1 −1 x1 2
Ax =  1 5 2   x2  = b =  17  (7)
2 −1 −6 x3 −18

Pour trouver la solution exacte, rappelons que La décomposition LU de A est donnée par
    
1 0 0 3 1 −1 3 1 −1
LU =  13 1 0   0 14 3
7 
3
= 1 5 2 =A
2 5 9
3
− 14 1 0 0 −2 2 −1 −6

  
1 0 0 y1
 
2


 y1 = 2
 1  
 49
1 0  y2  =  17 y2 = 17 − 31 y1 =
  
Ly =  =b ⇒

3  3
2 5
   
 
−18

 27
− 1 y3  y3 = −18 − 2 × 2 +
 5
× 49
= −
3 14 4 14 3
2

de même U x = y
 

31 −1

x1
 2   
x1 1
14 7 49
 x2   2  = Xe solutionn exacte
    
Ux =  0 3
 x      
 2  =  =y ⇒

 3
9    3     du système
0 0 − x3
 29 
x3 3
2 −
2

alors la décomposition A = D − E − F de la matrice A est donnée par

     
3 0 0 0 0 0 0 −1 1
D= 0 5 0  E =  −1 0 0  et F =  0 0 −2 
0 0 −6 −2 1 0 0 0 0

et l’itération de Jacobi xk+1 = D−1 [(E + F )xk + b] du système (??) s’écrit:

 k+1   1    k   
x1 3
0 0 0 −1 1 x 1 2
xk+1 =  xk+1
2
= 0 1 0   −1 0 −2   xk2  +  17 
5
1
xk+1
3 0 0 −6 −2 1 0 xk3 −18

ce qui conduit pour x0 = (x01 , x03 , x03 )t donné, aux formules itératives suivantes:

1

 xk+1
1 = ( 2 −xk2 +xk3 )
3


1


k+1
x2 = ( 17 −xk1 −2xk3 ) (8)
 5
 1
 xk+1

= − ( −18 −2xk1 +xk2 )

 3
6

Notons que les formules (??) ci-dessus ne sont autre que l’application directes de celles
données par (??). Et leurs applications donnent les résultats numériques suivant:

M. H. EL ALJ 24 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

Et pour X 0 = (0, 0, 0)t la première itération donne


1 2

 x11 = ( 2 −0 +0 )= = 0.666667
3 3




1 17

x12 = ( 17 −0 −2 × 0 ) = = 3.400000

 5 5

 1 18
 x13 = − ( −18 −2 × 0 +0 )= = 3.000000


6 6

et la seconde itération donne


1 17 8

 x21 = ( 2 − +3 )== 0.533333
3 5 15




1 2 31

x22 = ( 17 − −2 × 3 ) = = 2.066667

 5 3 15

 2 1 2 17 239
 x3 = − ( −18 −2 × + )= = 2.655556


6 3 15 90

de façon plus détaillé, nous avons

On retrouve les deux premières itérations sur les colonnes X 1 et X 2 .

X 0 est la valeur initiale de départ pour déterminer les élément de la suite (X n )n∈N .
et Xe représente la solution exacte
L’erreur est calculée par la formule suivante

Erreurk = Max |Xik − Xie | i = 1, . . . , 3

Dans le cas oú la solution X e n’est pas connue, l’incertitude oú la précision peut être déter-
minée par
Ek = Max |Xik − Xik−1 | i = 1, . . . , 3

De même les formules itératives de la méthode de Gauss-Seidel données par (??) s’écrivent
également comme suit:

1

 xk+1
1 = ( 2 −xk2 +xk3 )
3


1


xk+1
2 = ( 17 −xk+1
1 −2xk3 ) (9)
 5
 1
 xk+1 ( −18 −2xk+1 +xk+1

= − )

 3 1 2
6

M. H. EL ALJ 25 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

Et pour X 0 = (0, 0, 0)t la première itération donne


1 2

 x11 = ( 2 −0 +0 )= = 0.666667
3 3




1 2 49

x12 = ( 17 − −2 × 0 ) = = 3.266667

 5 3 15
 x3 = − 1 2 49 241

 1
( −18 −2 × + )= = 2.677778

6 3 15 90


et la seconde itération donne


1 49 241 127

 x21 = ( 2 − + )= = 0.470373
3 5 90 270




1 127 241 3017

x22 = ( 17 − −2 × )= = 2.234815

 5 270 90 1350

 1 127 3017 22553
 x23 = − ( −18 −2 × + )= = 2.784321


6 270 1350 8100

de façon plus précise on a les résultats numériques suivants:

Solution de l’exercice 9 Méthodes itératives

1. Reprise de l’exemple du Cours avec X 0 = (1, 1, 1)t

   
Soit : 3 1 −1 2
A= 1 5 2  et b =  17 
2 −1 −6 −18
L’itération de Jacobi s’écrit:

1

 xk+1
1 = ( 2 −xk2 +xk3 )
3




1

xk+1
2 = ( 17 −xk1 −2xk3 )

 5
 1
 xk+1

= − ( −18 −2xk1 +xk2 )

 3
6

d’oú: 1 1 2

 x 1 = ( 2 −1 +1 ) = = 0.6666667
3 3




1 14

x12 = ( 17 −1 −2 × 1 ) = = 3.4

 5 5

 1 19
 x13 = − ( −18 −2 × 1 +1 )= = 3.166667


6 6

M. H. EL ALJ 26 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

et la seconde itération donne


1 71

 x21 = ( 2 − 14
5
+ 19
6
)= = 0.788889
3 90




1 30

x22 = ( 17 − 23 −2 × 19
6
)= = 2.0

 5 15

 1 2 248
 x23 = − ( −18 −2 × + 14 )= = 2.755556


3 5
6 90

De façon pius précise, nous avons les résultats suivants

L’itération de Jacobi s’écrit:


 xk+1 1

1 = ( 2 −xk2 +xk3 )
3




1

xk+1
2 = ( 17 −xk+1
1 −2xk3 )

 5
 1
 xk+1 ( −18 −2xk+1 +xk+1

= − )

 3 1 2
6

d’oú: pour X 0 = (1, 1, 1)t , la première itération donne


 x11 1 2
 = ( 2 −1 +1 )= = 0.666667
3 3




1 2 43

x12 = ( 17 − −2 × 1 ) = = 2.866667

 5 3 15

 1 2 43 247
 x13 = − ( −18 −2 × + )= = 2.744444


6 3 15 90

et la seconde itération donne


1 43 169

 x21 = ( 2 − + 247
90
)= = 0.625926
3 15 270




1 169 247 2939

x22 = ( 17 − −2 × )= = 2.177037

 5 270 90 1350

 1 169 2939 23051
 x23 = − ( −18 −2 × + )= = 2.845802


6 2703 1350 8100

M. H. EL ALJ 27 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

De façon pius précise, nous avons les résultats suivants

2. Décomposition LU de A4 x = b
    
On a: 2 −1 0 0 x1 3
 −1 2 −1 0   x2   −5 
Ax = 
 0 −1   x3  = b =  7 
    (10)
2 −1
0 0 −1 2 x4 −6

Pour trouver la solution exacte, rappelons que La décomposition LU de A est donnée par
    
1 0 0 0 2 −1 0 0 2 −1 0 0
 −1 3
−1
1 0 0  0  =  −1
0  2 −1 0 
 
LU =  2 2  = A4
 0 − 2 4
3
1 0  0 0 3
−1   0 −1 2 −1 
0 0 − 43 1 0 0 0 5
4
0 0 −1 2
 
y
 
1 0 0 0 1

3

 y1 = 3
 1 1 7

  

 − 1 0 0   y2   −5 
  
 y2 = −5 + y1 = −


 2  2 2
 0 −2
Ly =   = =b ⇒ 2 14
  
y
1 0  3   7  y3 = 7 + 3 y2 =
 
 
 3      3
 3 
−6

 3 5
y
 
0 0 − 1 y4 = −6 + y3 = −

4 
4 4 2
et      4 5
2 −1 0 0 x1 3  x4 = ×− = −2
5 2

3 7


 x3 = 3 ( 14 + x4 )
     
 0 −1 0   x2   −  

= 2
 2    2 
4 3
Ux =  4   = 14 =y ⇒ 2 7
 0 0 −1   x3   
 3    3  
 x 2 = (− + x3 ) = −1
3 2

 5  x4
  5  
 x = 1 (3 + x )

0 0 0 −


4 2 1 2 = 1
2
La solution exacte du système A4 x = b est donc X e = (1, −1, 2, −2)t

3. Méthodes itératives de A4 x = b
Par ailleurs la décomposition A = D − E − F de la matrice A est donnée par
     
2 0 0 0 0 0 0 0 0 1 0 0
 0 2 0 0   1 0 0 0   0 0 1 0 
D=  0 0 2 0  E= 0 1
   et F =  
0 0   0 0 0 1 
0 0 0 2 0 0 1 0 0 0 0 0

M. H. EL ALJ 28 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

• Méthode de Jacobi
l’itération de Jacobi xk+1 = D−1 [b + (E + F )xk ] s’écrit:
 k+1   1      k 
x1 2
0 0 0 3 0 1 0 0 x1
k+1
 xk+1
2
  0 1 0 0   −5   1 0 1 0   xk2 
x =  k+1  = 
   2   +  
x3 0 0 12 0   7   0 1 0 1   xk3 
xk+1
4 0 0 0 12 −6 0 0 1 0 xk4

on a donc les formules itératives de Jacobi suivantes:

1


 xk+1
1 = ( 3 +xk2 )



 2
 xk+1 = 1


( −5 +xk1 +xk3 )

2
2 (11)
k+1 1
+xk2 +xk4


 x 3 = ( 7 )
2




 k+1 1
 x = ( −6 +xk3 )

3
2
Ces formules (??) sont l’application directes des formules (??).
D’oú pour X 0 = (0, 0, 0)t , la première itération donne
1 3


 x11 = ( 3 +0 ) = = 1.5



 2 2
 x1 = 1 ( −5 +0 5


+0 ) = − = −2.5

2
2 2
1 7
x13 =


 ( 7 +0 +0 ) = = 3.5
2 2



 x1 = 1 ( −6 6


+0 ) = − = −3

3
2 2
La Deuxième itération donne
1 5 1


 x21 = ( 3 − ) = = 0.25



 2 2 4
 x2 = 1 ( −5 + 3 7 0


+ ) = = 0

2
2 2 2 2
 2 1 21

 x 3 = ( 7 +0 +0 ) = = 2.625
2 8



 x2 = 1 ( −6 27


+0 ) = − = −1.6875

3
2 6
Leurs applications donnent les résultats numériques suivants:

M. H. EL ALJ 29 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

Et pour X 0 = (1, 1, 1)t on a

• Méthode de Gauss Seidel

l’itération de Jacobi Dxk+1 = [b + Exk+1 + F xk ] s’écrit:

2xk+1
       k+1    k 
1 3 0 0 0 0 x1 0 1 0 0 x1
 k+1 
2x2   −5   1 0 0 0   xk+1
     2
  0 0 1 0   xk2 
 = + k+1  + 
   
 2xk+1
3
  7   0 1 0 0   x 3 0 0 0 1   xk3 
2xk+1
4 −6 0 0 1 0 xk+1
4 0 0 0 0 xk4

ce qui conduit pour x0 = (x01 , x03 , x03 )t donné, aux formules itératives suivantes:
1

 xk+1

1 = ( 3 −xk2 )



 2


 xk+1 = 1
( −5 −xk+1 −xk3 )

2 1
2 (12)
k+1 1
−xk+1 −xk4 )


 x 3 = − ( 7 2
2



 xk+1 = − 1 ( −6


−xk+1 )

3 3
2
Ces formules (??) sont l’application directes des formules (??).

D’oú pour X 0 = (0, 0, 0)t , la première itération donne

1 3


 x11 = ( 3 +0 ) = = 1.5



 2 2
 x1 = 1 ( −5 + 3 7


+0 ) = − = −1.75

2
2 2 4
 1 1 7 21
 x3 =

 ( 7 − +0 ) = = 2.625

 2 4 8
 x1 = 1 ( −6 21 27


+ ) = − = −1.6875

3
2 8 16

M. H. EL ALJ 30 Mars 2022

UIT. FSK. Département de Mathématiques Filière SMI : S4

La Deuxième itération donne

1 7 5

 x21 =
 ( 3 − ) = = 0.625



 2 4 8
 x2 = 1 ( −5 + 5 21 7


+ ) = − = −0.875

2
2 8 8 8
 2 1 7 27 71

 x 3 = ( 7 − − ) = = 2.21875
2 8 16 32



 x2 = 1 ( −6 7 121


+ ) = − = −1.8906

3
2 32 64
Leurs applications donnent les résultats numériques suivants:

De même pour X 0 = (1, 1, 1)t on a les résultats suivants :

M. H. EL ALJ 31 Mars 2022