TD1 An

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 15

Ecole Centrale Supérieure Polytechnique Privée de Tunis

3 Avenue Mohamed V 2015 Le Kram Tunisie

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

Classe : ING_1 Filière :

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 

1. Calculer cond2(A), cond1(A) et cond00(A)

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

1- Résoudre le système ci-dessous :

2 x1  3x2  4 x3  3

2 x1  2 x2  2 x3  4
2 x  5 x  10 x  1
 1 2 3

En utilisant les formules de Cramer puis la méthode de GAUSS

2- Est-il possible de résoudre par les formules de Cramer le système suivant ?

 x1  2 x2  3x3  1

2 x1  2 x2  2 x3  2
x  4x  7x  5
 1 2 3

Exercice 4

On considère le système linéaire Ax  b donné par :

 y z
x  2  3  1

x y z
   0
2 3 4
x y z
3  4  4  0

1. Vérifier que la matrice A est inversible

2. Résoudre le système en utilisant la méthode de GAUSS

Exercice 5

Soit le systéme décrit ci-dessous :

 1 1 4  x1   3 
    
 2 3 4  x2   14 
 0 2 8  x   0 
  3   

Résoudre le système en utilisant la méthode de GAUSS Jordan


Page 2 sur 15
Exercice 6

On considère le système linéaire A x  b avec:

1 1 1  b1 
  0, A   0 1 1  et b   b2 
 
1 0   b 
   3

1- Dans quel cas A est une matrice inversible ?

2- Montrer que dans ce cas A admet une décomposition LU

3- Résoudre le système A x  b pour bt  1 2 1

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

Indication : Le rayon spectrale de A est :  ( A)  max   ,   S p ( A)

Détermination des valeurs propres de A :

P( )  det( A   I 2 ) : polynome caractéristique

 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 

1. Calculer cond2(A), cond1(A) et cond00(A)

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)(1012   )
  2    1012
 1 0
1  0, 1  1
A 2 1

Détermination de cond1(A),

A 1 1
A1  ?
1

1 0 
A1   6
 0 10 
A1  106
1

cond1 ( A)  A1  A 1  106


1

Détermination de cond00(A),

A  1
A1 ?

1 0 
A1   6
 0 10 
A1  106

cond  ( A)  A1  A   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  104  0
1.0001 5
 Aest inversible
( x  x)  A1 (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 A1   
 0.2 0.2 
A1 2

cond  ( A)  A1  A   12.0002


On a bien la relation des majorations des perturbations :

x b

 cond ( A) 
 12.0002  0  0
x 
b 

Exercice 3

1- Soit le système d’équations suivant :

2 x1  3x2  4 x3  3

2 x1  2 x2  2 x3  4 (1)
2 x  5 x  10 x  1
 1 2 3

a- Résolution par la méthode de Cramer

Indication :La régle de Cramer

det( A j )
xj  avec j  1,.., n
det( A)

En appliquant la régle de Cramer sur le systéme décrit ci-dessus :

  
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

 2(10)  3(16)  4(6)   4

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

 3(10)  3(42)  4(22)  8


  
2 3 4
4 2 2 2 2 4
 2  det( A )  2 4 2  2
2
3 4
1 10 2 10 2 1
2 1 10

 2(42)  3(16)  4(10)   4

  
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

 2(22)  3(10)  3(6)   4

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

b- Résolution par la méthode de GAUSS

Le systéme (1) s’écrit encore Ax  b avec

 2 3 4  3  x1 
     
A   2 2 2  ; b   4  , X   x2 
 2 5 10   1  x 
     3

1 ere étape : a11(1)  2  0

  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) 

Résolution de A(3) x  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  

2- Est-il possible de résoudre par les formules de Cramer le système suivant ?

 x1  2 x2  3x3  1

2 x1  2 x2  2 x3  2
x  4x  7x  5
 1 2 3

En appliquant la régle de Cramer sur le systéme décrit ci-dessus :

  
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

 (6)  2(12)  3(6)   0

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

On considère le système linéaire Ax  b donné par :

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. Vérifier que la matrice A est inversible

 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

2. Résoudre le système en utilisant la méthode de GAUSS

 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

1 ere étape : a11(1)  1  0

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

Soit le systéme décrit ci-dessous :

 1 1 4  x1   3 
    
 2 3 4  x2   14 
 0 2 8  x   0 
  3   

Résoudre le système en utilisant la méthode de GAUSS Jordan

1 ere étape : a11(1)  1  0

 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) 

2éme étape : a22(2)  1  0

 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) 

3éme étape : a33(3)  32  0

 
 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)

Résolution du systéme A(3) X  b(3)

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

On considère le système linéaire A x  b avec:

1 1 1  b1 
  0, A   0 1 1  et b   b2 
 
1 0   b 
   3

1- Dans quel cas A est une matrice inversible ?

det( A )   donc A est inversible si   0

2- Montrer que dans ce cas A admet une décomposition LU

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  
   

3- Résoudre le système A x  b pour bt  1 2 1

 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 : L1  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

Lorsque   0 la seule condition pour que x soit finie est b3  b2  b1  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. Vérifier si A est définie positive

D’après le cours , on va vérifie si les mineurs de la matrice sont strictement positifs :

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

La matrice A est définie positif,donc A admet une décomposition de Cholesky

Soit A  LLT avec L est une matrice triangulaire inférieure :

 l11 0 0  l11 l21 l31 


  
A   l21 l22 0  0 l22 l32 
l l33  
 31 l32  0 0 l33 
 l112 l11l21 l11l31 
 
  l11l21 l21  l22 2
2
l21l31  l22l32 
 l11l31 l21l31  l22l32 l312  l32 2  l332 

Page 14 sur 15
Par identification

l112  2  l11  2  1.4142


1
l11l21  1  l21   0.7071
2
l11l31  0  l31  0
1 2 3 3
l212  l22 2  2  l22 2  2  ( )   l22   1.2247
2 2 2
1
l21l31  l22l32  1  l32   0.8165
1.2247
l312  l22 2  l332  2  l332  2  l312  l22 2  2  1.5  0.5  l33  0.7071

En résume

 1.4142 0 0  1.4142 0.7071 0 


   
L   0.7071 1.2247 0  et L   0
T
1.2247 0 
 0 0.8165 0.7071  0 0.8165 0.7071
 

Page 15 sur 15

Vous aimerez peut-être aussi