0% ont trouvé ce document utile (0 vote)
504 vues15 pages

TD1 An

Ce document contient les énoncés de plusieurs exercices d'analyse numérique portant sur la résolution de systèmes linéaires, le calcul du rayon spectral de matrices et le conditionnement de matrices. Les exercices sont corrigés.

Transféré par

bouzidmohamedali119
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd
0% ont trouvé ce document utile (0 vote)
504 vues15 pages

TD1 An

Ce document contient les énoncés de plusieurs exercices d'analyse numérique portant sur la résolution de systèmes linéaires, le calcul du rayon spectral de matrices et le conditionnement de matrices. Les exercices sont corrigés.

Transféré par

bouzidmohamedali119
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd
Vous êtes sur la page 1/ 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