TD1 An
TD1 An
TD1 An
Année universitaire :
Département : IEEE Semestre : 2
Module d’enseignement : Analyse Numérique Date :
Durée : 1h
Enseignant : Mme. Khaoula ASSADI EL ANDOLSI Lieu : Visio-conférence
Type de séance : Travaux dirigés TD N°1
Nombre de pages : 12 p
Enoncés
Exercice 1
1 1 1 0
Soit A et B Calculer le rayon spectrale de A et B
1 1 1 1
Exercice 2
1 0
A. On considére la matrice A 6
0 10
5
1 5 x1 6
B. Soit le systéme linéaire x de solution exacte x 1
1.0001 5 2 6.0005
5
1. Trouver la nouvelle solution du système ( x x) lorsqu’on remplace le membre droite du système par
5
6
(b b) et x 1
6
5
b
x
2. Calculer les valeurs de et que peut-on conclure ?
b
x
Page 1 sur 15
3. Calculer le conditionnement de ce système pour la norme ∞et vérifier la relation des majorations des
perturbations
Exercice 3
2 x1 3x2 4 x3 3
2 x1 2 x2 2 x3 4
2 x 5 x 10 x 1
1 2 3
x1 2 x2 3x3 1
2 x1 2 x2 2 x3 2
x 4x 7x 5
1 2 3
Exercice 4
y z
x 2 3 1
x y z
0
2 3 4
x y z
3 4 4 0
Exercice 5
1 1 4 x1 3
2 3 4 x2 14
0 2 8 x 0
3
1 1 1 b1
0, A 0 1 1 et b b2
1 0 b
3
4- Indiquer à quelle condition nécessaire sur b, le système A x b admet une solution lorsque 0
Exercice 7
Vérifier que la matrice suivante est définie positive et calculer sa décomposition de Cholesky
2 1 0
A 1 2 1
0 1 2
Page 3 sur 15
Correction
Exercice 1
1 1 1 0
Soit A et B Calculer le rayon spectrale de A et B
1 1 1 1
1 1 0
PA ( )
1 1 0
1 1
2 2 ( 2)
1 1
0
PA ( ) 0 1
2 2
( A) 2
1 0 0
PB ( )
1 1 0
1 0
2 2 1
1 1
b
racine double 1 2 1
2a
( A) 1
Exercice 2
1 0
A. On considére la matrice A 6
0 10
Page 4 sur 15
Détermination de cond2(A),
Indication : A 2 ( At A)
On pose B At A
1 0
B 12
0 10
PB ( ) ( 1)(1012 )
2 1012
1 0
1 0, 1 1
A 2 1
Détermination de cond1(A),
A 1 1
A1 ?
1
1 0
A1 6
0 10
A1 106
1
Détermination de cond00(A),
A 1
A1 ?
1 0
A1 6
0 10
A1 106
5
1 5 x1 6
B. Soit le systéme linéaire x de solution exacte x 1
1.0001 5 2 6.0005
5
Page 5 sur 15
1. Trouver la nouvelle solution du système ( x x) lorsqu’on remplace le membre droite du système par
5
6
(b b) et x 1
6
5
A( x x) b b
det( A) ?
1 5
5 5.0005 5 104 0
1.0001 5
Aest inversible
( x x) A1 (b b)
0
1 1 6 0
6
0.2 0.2 6 1.2
5
b
x
2. Calculer les valeurs de et que peut-on conclure ?
b
x
6 6 0
b
6 6.0005 0.0005
b 0
b
6.0005
b 0
0
b
6.0005
0 5 5
x
1.2 0.2 1
x 5
x
5
x 5
1
x
5
On voit que la variation relative de x est beaucoup plus grande que celle de b autrement dit une petite
variation sur le second membre peut produire une grande variation sur la solution.
3. Calculer le conditionnement de ce système pour la norme ∞et vérifier la relation des majorations des
perturbations
Page 6 sur 15
1 5
A
1.0001 5
A 6.0001
1 1
Soit A1
0.2 0.2
A1 2
x b
cond ( A)
12.0002 0 0
x
b
Exercice 3
2 x1 3x2 4 x3 3
2 x1 2 x2 2 x3 4 (1)
2 x 5 x 10 x 1
1 2 3
det( A j )
xj avec j 1,.., n
det( A)
2 3 4
2 2 2 2 2 2
det( A) 2 2 2 2 3 4
5 10 2 10 2 5
2 5 10
Page 7 sur 15
3 3 4
2 2 4 2 4 2
1 det( A ) 4 2 2 3
1
3 4
5 10 1 10 1 5
1 5 10
2 3 3
2 4 2 4 2 2
3 det( A ) 2 2 4 2
3
3 3
5 1 2 1 2 5
2 5 1
1 8
x1 2
det( A) 4
2
2 4
x2 1 X 1
det( A) 4 1
3 4
x3 1
det( A) 4
2 3 4 3 x1
A 2 2 2 ; b 4 , X x2
2 5 10 1 x
3
2 3 4 3 2 3 4 3
2 2 2 4 0 1 2 1
2
5 10 1 0 2 6 4
A(1) ; b(1) A(2) ; b (2)
Page 8 sur 15
2éme étape : a22(2) 1 0
2 3 4 3 2 3 4 3
0 1 2 1 0 1 2 1
0 6 4 0 0 2 2
2
A(2) ; b(2) A(3) ; b (3)
2 3 4 x1 3
A(3) x b(3) 0 1 2 x2 1
0 0 2 x 2
3
x1 1
Ce qui donne x2 1
x 2
3
x1 2 x2 3x3 1
2 x1 2 x2 2 x3 2
x 4x 7x 5
1 2 3
1 2 3
2 2 2 2 2 2
det( A) 2 2 2 2 3
4 7 1 7 1 4
1 4 7
On obtient det(A )= 0 (matrice singuliére) ce qui ne permet pas d’utiliser les formules de Cramer meme
aussi par la méthode de Gauss (A est non inversible)
Exercice 4
Page 9 sur 15
y z
x 2 3 1
x y z
0
2 3 4
x y z
3 4 4 0
1 1
1 2 3
Soit A
1 1 1
2 3 4
1 1 1
3 4 4
1 1
1
2 3
1 1 1 1
det( A) 0.004 0
2 3 4 216
1 1 1
3 4 4
A est inversible
1 1
1 2 3
1 x
A
1 1 1
; b 0 , X y
2 3 4
0 z
1 1 1
3 4 4
Page 10 sur 15
1 1 1 1
1 2 3
1 1
2 3
1
1 1 1
0 0
1 1
1
2 3 4 12 12 2
1 1 1 1 5 1
0 0
3 4 4 12 36 3
A ; b A ; b
(1) (1) (2) (2)
1
2éme étape : a22 (2) 0
12
3 L3 L2
L(3) (2) (3)
1 5 1 1 1 1
0 0
12 36 3 12 12 2
1 1
0 0
18 6
1 1 1 1
1 2 3
1 1
2 3
1
0 1 1
0
1 1 1
1
12 12 2 12 12 2
0 1 5 1
0 0
1 1
12 36 3 18 6
A(2) ; b (2) A(3) ; b (3)
1 1
1 2 3 1
x
1 1 1
A X b 0
(3) (3)
y
12 12 2
z
0 0 1 1
18 6
1 6 3
(3) 1
( A ) 0 12 18
0 0 18
(3) 1
X (A ) b
4.5
9
3
Page 11 sur 15
Exercice 5
1 1 4 x1 3
2 3 4 x2 14
0 2 8 x 0
3
1 1 4 3 1 1 4 3
2 3 4 14 0 1 12 8
0 2 8 0 0 2 8 0
A(1) ; b(1) A(2) ; b (2)
1 1 4 3 1 0 8 5
0 1 12 8 0 1 12 8
0 2 8 0 0 0 32 16
A(2) ; b(2) A(3) ; b(3)
1 0 8 5 1 0 0 1
0 1 12 8 0 1 0 2
0 0 32 16
0 0 1 1
2
A ; b A ; b
(3) (3) (4) (4)
Page 12 sur 15
1 0 0 x1 1
A(3) X b (3) 0 1 0 x2 2
0 0 1 x 1
3
2
X b (3)
1
2
1
2
Exercice 6
1 1 1 b1
0, A 0 1 1 et b b2
1 0 b
3
1 0 0 1 1 1
Soit A LU et par identification : L 0 1 0 etU 0 1 1
1 1 1 0 0
y t (b1 b2 b3 b2 b1 )
Ly b
A x b t (b3 b2 b1 ) (b3 b2 b1 ) 2 2
Ux y x (b1 b2 b2 ) 1 2
Page 13 sur 15
1 1 1
1 0 0
Avec : L1 0 1 0
1 1
etU 0 1
1 1 1
1
0 0
4- Indiquer à quelle condition nécessaire sur b, le système A x b admet une solution lorsque 0
Exercice 7
Vérifier que la matrice suivante est définie positive et calculer sa décomposition de Cholesky
2 1 0
A 1 2 1
0 1 2
1 a11 2 0
a11 a12 2 1
2 3 0
a21 a22 1 2
2 1 0
2 1 1 1
3 1 2 1 2 4 0
1 2 0 2
0 1 2
Page 14 sur 15
Par identification
En résume
Page 15 sur 15