Mr. ADOUI Salah / 2022. CORRECTION DES EXERCICES Dépt Mathématiques.

L2- GROUPES : 1, 2, 4 et 5 TD : 02 Analyse numérique 2

Exercice 1
Soit le systéme linéaire dont la
matrice carrée d’ 
ordre 4 associée est :
6 2 1 −1
 2 4 1 0 
A=  1 1 4 −1


−1 0 −1 3
1 La matrice A est diagonalement dominante au sens strict car :
∀i = 1, 2, 3, 4 : | aii | ≥ ∑ j6=i | aii |, j = 1, 2, 3, 4
pour i=1 : | a11 | = 6 ≥ | a12 | + | a13 | + | a14 | = 2 + 1 + 1 = 4.
pour i=2 : | a22 | = 4 ≥ | a21 | + | a23 | + | a24 | = 2 + 1 + 0 = 3.
pour i=3 : | a33 | = 4 ≥ | a31 | + | a32 | + | a34 | = 1 + 1 + 1 = 3.
pour i=4 : | a44 | = 3 ≥ | a41 | + | a42 | + | a43 | = 1 + 0 + 1 = 2.

2 Calcul de l’inverse de A , avec la méthode de Gauss-Jordon on a :

01 TD [Link]

1
 3 Factorisation de A en produit de deux matrices L et U :

01 q3 TD [Link]
4 En déduire le déterminant de A sachant que A=L ×U :
On a det ( A) = det ( L.U ) = det ( L) × det (U ) = ∏ lii × ∏ uii ;
donc det ( A) = det (U ) = ∏ lii = 6 × 10 37 191
3 × 10 × 74 = 191.

Exercice 2
Étant donnée le systéme déquation linéaire AX=B tels que la matrice carrée :
   
3 8 9 6
A = 3 8 1 ; b = 11
  
1 2 4 20
1 Résolution approximativement le systéme par la méthode Jacobi pour un vecteur
initiale X (0) = (1; 0; 0) :

 
1
3x + 8y + 9z = 6
  x = 3 (6 − 8y − 9z)

Ax = b ⇐⇒ 3x + 8y + z = 11 ⇐⇒ y = 81 (11 − 3x − z)
z = 14 (20 − x − 2y)
 
x + 2y + 4z = 20
 

En appliquant le processus de Jacobi avec deux iterations on obtient :


(k+1) = 1 (6 − 8y(k) − 9z(k) )
x
 3
y(k+1) = 81 (11 − 3x (k) − z(k) )
= 14 (20 − x (k) − 2y(k) )
 ( k +1)

z

1ere itteration : pour k=0


(1) 1 (0) (0)
 x = 3 (6 − 8y − 9z ) = 2

y(1) = 81 (11 − 3x (0) − z(0) ) = 1 .
z = 14 (20 − x (0) − 2y(0) ) = 19
 (1)

4

2eme itteration : pour k=1


(2) 1 (1) (1) −179
 x = 3 (6 − 8y − 9z ) = 12

y(2) = 81 (11 − 3x (1) − z(1) ) = 32
1 .
z = 41 (20 − x (1) − 2y(1) ) = 4
 (2)


2
 2 Détermination les valeurs propres de A :

(3 − λ ) 8 9
λ valeur propre de A ⇐⇒ 3 (8 − λ ) 1 =0
1 2 (4 − λ )

λ1 = −0.407

⇐⇒ −(λ + 0.407)(λ − 3.224)(λ − 12.184) = 0 ⇐⇒ λ2 = 3.224

λ3 = 12.184


3 calcul du rayon spectral de la matrice A :

ρ( A) = MaxλvpA |λi | = 12.184.
4 En dduire k A k
p2 :
On a : k Ak2 = ρ( At .A) = ρ( H ) ; avec H = At .A
p

et comme A n’est pas symetrique , on va chercher les valeurs propre de H :

Exercice 3
On considére le
 systéme linéaire
 (S) :AX= b :  
2 −1 0 x 1
avec A = −1 2 −1 ;X =  y  et b = 0
0 −1 2 z 1
On cherche résoudre (S) par les méthodes indirectes de type : Xk+1 = MXk + D :
1 A est symétrique car : At = A.
A est définie positive car :
∆1 = a11 = 2 > 0,
a a 2 −1
∆2 = 11 12 = a11 a22 − a12 a21 = = 3 > 0,
a21 a22 −1 2
a11 a12 a13
∆3 = a21 a22 a23 = det( A) = 4 > 0.
a31 a32 a33
2 A est diagonalement dominante, car :
∀i = 1, 2, 3 : | aii | ≥ ∑ j6=i | aii |, j = 1, 2, 3
pour i=1 : | a11 | = 2 ≥ | a12 | + | a13 | = 1.
pour i=2 : | a22 | = 2 ≥ | a21 | + | a23 | = 2.
pour i=3 : | a33 | = 2 ≥ | a31 | + | a32 | = 1.

3 Résolution du systéme (S) avec la méthode ittérative de Jacobi en partant du vecteur de

départ initiale X (0) = (0; 0; 0) et avec 2 ittérations :
Xk+1 = JXk + D , avec :  
a13 a1n 
0 aa12 · · · · · · b1

11 a11 a11 a11
 a21 0 a23 · · · · · · a2n   b2 
 a22 a22 a22  a 
 22 
 . ..
. 0 · · · · · · a22  ; D =  ... .
a2n 

J = − . .
  
 .. .. .. ..   . 
 . 
 . . . ··· 0 .   . 
an1 an2 a13 bn
ann ann a11 · · · · · · 0
ann

3
  −1  1
0 2 0 2
Donc : J = −  −21 0 −21  ; D =  0 .
0 −21 0 1
2
1ere itteration : pour k=0
0 −21 0
    1 1
0 2 2
X1 = JX0 + D =  −21 0 −21  × 0 +  0  =  0 
0 −21 0 0 1
2
1
2
2eme itteration : pour k=1
0 −21 0
  1 1  1 
2 2 2
X2 = JX1 + D =  −21 0 −21  ×  0  +  0  =  −41 
0 −21 0 1
2
1
2
1
2
4 Récolution du systéme (S) avec la méthode ittérative de Gauss-Seidel. en partant du
vecteur de départ initiale X (0) = (0; 0; 0) et avec 2 ittérations :

1
 x = 2 (1 + y )

Ona : (S) ⇔ y = 21 ( x + z) .
z = 12 (1 + y)



On appliquant la méthode itérative de Gauss seidel :


( k +1) = 1 (1 + y ( k ) )
x
 2
(S) ⇔ y(k+1) = 12 ( x (k+1) + z(k) ) .
= 21 (1 + y(k+1) )
 ( k +1)

z
1ere itteration : pour k=0

(1) 1 (0) 1
 x = 2 (1 + y ) = 2

y(1) = 12 ( x (1) + z(0) ) = 1
4
z = 12 (1 + y(1) ) = 58
 (1)


2eme itteration : pour k=1


(2) 1 (1) 5
 x = 2 (1 + y ) = 8

y(2) = 12 ( x (2) + z(1) ) = 5
8
z = 12 (1 + y(2) ) = 13
 (2)

16

5 Quelle est celle qui converge plus rapidement , justifier ? On remarque que la solution
exacte du systéme est x (exacte)=(1,1,1) n donc si on calcul L’erreur :
itteration 1 pour k=0 :
EGS = k x (exacte) − x (0) k∞ = k(1 − 21 ; 1 − 14 ; 1 − 58 )k∞ = 21
E Jacobi = k x (exacte) − x (1) k∞ = k(1 − 12 ; 1 − 0; 1 − 21 )k∞ = 1
itteration 2 pour k=1 :
EGS = k x (exacte) − x (2) k∞ = k(1 − 58 ; 1 − 58 ; 1 − 13 16 )k∞ = 8
3

E Jacobi = k x (exacte) − x (2) k∞ = k(1 − 12 ; 1 − 14 ; 1 − 12 )k∞ = 34

On remarque que dans tout les iterations EGS < E Jacobi , donc la méthode itérative de
Gauss seidel converge plus rapide que celle de Jacobi.

4
Exercice 4
Soit le systéme d’équations linéaires :

1 1 11
x + 2 y + 4 z = 4

(S)............ −21 x + y = 32 .
1
 1
2x + 4y + z = 4

1 Résolution du systéme par la méthode itérative de Gauss seidel en partant du vecteur

de départ initiale unité. Prendre la précision : ε = 10−2 .

11 1 1
x = 4 − 2 y − 4 z

Ona : (S) ⇔ y = 32 + 12 x .
z = 4 − 12 x − 41 y



On appliquant la méthode itérative de Gauss seidel :


(k+1) = 11 − 1 y(k) − 1 z(k)
x
 4 2 4
(S) ⇔ y ( k + 1 ) 3 1 ( k ) .
= 2 + 2x
= 4 − 21 x (k) − 14 y(k)
 ( k +1)

z
1ere itteration : pour k=0

(1) 11 1 (0) 1 (0)
x = 4 − 2 y − 4 z = 2

y(1) = 32 + 21 x (1) = 52 .
z = 4 − 12 x (1) − 41 y(1) = 19
 (1)

8

L’erreur : k x (1) − x (0) k∞ = k(2 − 1; 52 − 1; 19 3 11 −2

8 − 1)k∞ = k(1; 2 ; 8 )k∞ = 1.5 > ε = 10

2eme itteration : pour k=1


(2) 11 1 (1) 1 (1) 29
 x = 4 − 2 y − 4 z = 32

y(2) = 32 + 21 x (2) = 93
64
.
z = 4 − 12 x (2) − 14 y(2) = 256
815
 (2)


29 5 35 167 −207
L’erreur : k x (1) − x (0) k∞ = k(2 − 32 ; 2 − 93 19 815 167
64 ; 8 − 256 )k∞ = k( 32 ; 64 ; 256 )k∞ = 64 >
ε = 10−2
...ect...on calcul X (3) , X (4 , X (5) ... jusqau on obtient L’erreur : k x (k) − x (k−1) k < ε = 10−2

2 Résolution du méme systéme par Jacobi sachant que la solution exacte est
X (exacte) = (1; 2; 3)

11 1 1
x = 4 − 2 y − 4 z

Ona : (S) ⇔ y = 32 + 12 x .
z = 4 − 12 x − 41 y



On appliquant la méthode itérative de Jacobi :


(k+1) = 11 − 1 y(k) − 1 z(k)
x
 4 2 4
(S) ⇔ y(k+1) = 32 + 12 x (k+1) .
= 4 − 21 x (k+1) − 14 y(k+1)
 ( k +1)

z

5
 1ere itteration : pour k=0

(1) 11 1 (0) 1 (0)
x = 4 − 2 y − 4 z = 2

y(1) = 32 + 21 x (0) = 2 .
z = 4 − 12 x (1) − 14 y(1) = 13
 (1)

4

L’erreur : k x (exacte) − x (0) k∞ = k(2 − 1; 2 − 2; 13 1 −2

4 − 3)k∞ = k(1; 0; 4 )k∞ = 1 > ε = 10
2eme itteration : pour k=1

(2) 11 1 (1) 1 (1) 15
 x = 4 − 2 y − 4 z = 16

y(2) = 32 + 21 x (1) = 52 .
z = 4 − 12 x (1) − 14 y(1) = 25
 (2)


15
L’erreur : k x (exacte) − x (2) k∞ = k( 16 − 1; 52 − 2; 25 − 3)k∞ = k( −161 ; 21 ; −21 )k∞ = 1
2 = 0.5 >
ε = 10−2
3eme itteration : pour k=2

(3) 11 1 (2) 1 (2) 17
x = 4 − 2 y − 4 z = 8

y(3) = 32 + 21 x (2) = 52 = 63
32
.
 (3) 1 (2) 1 (2) 93

z = 4 − 2 x − 4 y = 32

3 Étudier le critére de convergence de la méthode de Gauss seidel.

1 12 41
 

Comme la matrice A =  2 1 1  est diagonalement dominante car :

−1 0 1
∀i = 1, 2, 3 : | aii | ≥ ∑ j6=i | aii |, j = 1, 2, 3
pour i=1 : | a11 | = 1 ≥ | a12 | + | a13 | = 34 .
pour i=2 : | a22 | = 1 ≥ | a21 | + | a23 | = 12 .
pour i=3 : | a33 | = 1 ≥ | a31 | + | a32 | = 34 .

Donc la méthode de Gauss seidel converge vers la solution exacte du systéme (S).

6
Exercice 5
1 1
 3
8 0 4 8
On considre la matrice A et le vecteur b dfenis par : A =  0 13 0  ; b =  13 
0 −1 2 5
8
1 a Résolution du systéme AX=b par la méthode de Gauss.

2 exo [Link]
b Endéduire le déterminant de A :

1
Comme la matrice de Gauss et tria,gulaire ; donc : det (A)= a11 × a22 × a33 = 8 × 13 × 2
1
Donc det(A)= 12
c On ne peut pas utiliser la méthode de Cholesky pour retrouver la solution du
systéme AX=b car la matrice n’est pas symetrique.
2 a Le processus itératif de Jacobi correspondant au systéme AX=b :

5 [Link]

7
 5 [Link]

5 2a 3 .png
d’ou le type : Xk+1 = JXk + D , avec :  
a13
0 aa11
12
· · · · · · aa1n b1
 
a11 11 a11
 a21 0 a23 · · · · · · aa2n   b2 
 a22 a22 22 
a 
 22 
 . .. . 
. · · · · · · aa2n

J = − .
 . 0  .. .
;D = 
22 
 .. .
.. .. ..   . 
 . 
 . . ··· 0 .   . 
an1 an2 a13 bn
ann ann a11 ··· ··· 0
ann

b Soit X (0) = ( 21 ; 0; 1), calculer X (1) et X (2) . Donner l’erreur absolue commise dans
chacun des cas en norme k A k∞ :

On appliquant la méthode itérative de Jacobi :


(k+1) = 8( 3 − 1 z(k) ) = 3 − 2z(k)
x
 8 4
Ax = b ⇔ y(k+1) = 3. 13 = 1 .
= 12 ( 58 + y(k) )
 ( k +1)

z

1ere itteration : pour k=0


(1) 3 1 (0) (0)
 x = 8( 8 − 4 z ) = 3 − 2z = 1

y(1) = 3. 13 = 1 .
z = 21 ( 58 + y(0) ) = 16
5
 (1)


8
 L’erreur : k x (exacte) − x (1) k∞ = k( 11 13 5 3 1
8 − 1; 1 − 1; 16 − 16 )k∞ = k( 8 ; 0; 2 )k∞ = 0.5
2eme itteration : pour k=1

(2) 3 1 (1) (1) 19
 x = 8( 8 − 4 z ) = 3 − 2z = 8

y(2) = 3. 13 = 1 .
z = 12 ( 58 + y(1) ) = 1613
 (2)


L’erreur : k x (exacte) − x (2) k∞ = k( 11

8 −
19 13
8 ; 1 − 1; 16
13
− 16 )k∞ = k(1; 0; 0)k∞ = 1

c Etudier la convergence de cette méthode :

 
0 0 2
k J k∞ = k 0 0 0 k∞ = Max {2, 0, 21 } = 2 > 1;
0 −21 0
donc la methode de Jacobi ne converge pas.
3 On ne peut pas réécrire les systéme AX=b de maniére assurer la convergence de la
méthode de Jacobi ∀ X (0) ∈ R3 ; car si on permette les lignes de A toujours on tombe sur
aii = 0 pour un certain i ∈ {1, 2, 3}.

Exercice 6

1 Résolution par la méthode du pivot de Gauss le systéme suivant :

(
( a + 1) x + ( a − 1) y = 1
(S) ; a ∈ R∗
( a − 1) x + ( a + 1) y = 1
On a det( A) = 4a 6= 0, car : a ∈ R∗
donc le systéme admet une solution unique.
On appliquant la méthode du pivot de Gauss :
..
 
.. ! L1 → L1
( a + 1) ( a − 1) . 1 ( a + 1) ( a − 1) . 1
A= ..

1− a

→→  .. 
( a − 1) ( a + 1) . 1 L2 → L1 a+1 + L1 0 (
4a
a +1)
. 2
a +1
(
( a + 1) x + ( a − 1) y = 1 1
Donc : (S) ⇐⇒ 4a 2 ; alors :y = x = 2a
( a +1)
y = a +1

2 On pose a = 2. On Choisit la méthode itérative de GAUSS-SEIDEL pour résoudre le

systéme (S) en partant du vecteur de départ initiale ( x (0) ; y(0) ) = (0; 0) et avec deux
itérations :

( ( (
3x + y = 1 x + 3y = 1 x = 1 − 3y
(S) ⇐⇒ ⇐⇒ ⇐⇒
x + 3y = 1 y + 3x = 1 y = 1 − 3x
(
x (k+1) = 1 − 3y(k)
On appliquant la méthode de GAUSS-SEIDEL :
y(k+1) = 1 − 3x (k+1)
1ere itteration : pour k=0

9
 (
x (1) = 1 − 3y(0) = 1
.
y(1) = 1 − 3x (1) = −2

2eme itteration : pour k=1

(
x (2) = 1 − 3y(1) = 5
.
y(2) = 1 − 3x (2) = −14

10