TD1 - Analyse Numérique
TD1 - Analyse Numérique
TD1 - Analyse Numérique
Cycle Ingénieur
Analyse Numérique
Fiche de TD N°1
Exercice N°1
Ax = b (1)
où
6 1 2 −1
3 −8 −2 2
A =
2 0 5 2
−2 3 2 −9
1. Déterminer la décomposition LU de la matrice A.
2. En déduire det(A), ainsi que le nombre de solutions du système (1)
3. On donne b = (30, −14, 20, −4)T . Résoudre (1) suivant la décomposition LU de A.
Exercice N°2
Attention : cette matrice nécessite une permutation des lignes pour la décomposition LU
1. Écrire ce système sous forme matricielle.
On notera A la matrice qui en découle.
2. La matrice A est-elle à diagonale dominante stricte ? Justifier votre réponse.
3. Déterminer la décomposition LU de A.
4. En déduire le déterminant de la matrice A, ainsi que le nombre de solutions du
système (2).
5. Résoudre le système (1) en utilisant la décomposition LU de B.
Exercice N°3
N.B : On donnera les résultats avec cinq chiffres significatifs. On considère le système
12x1 − 2x2 + 3x3 = 29
−2x1 + 7x2 + x3 = 19 (3)
2x1 − 3x2 + 8x3 = 52
1
1. Écrire le système (3) sous la forme Bx = b. où B est la matrice qui en découle.
2. Étudier la convergence des méthodes de Jacobi et de Gauss-Seidel.
3. Écrire B sous la forme B = D − E − F , où
— D : est une matrice diagonale ;
— E : est une matrice triangulaire inférieure ;
— F : est une matrice triangulaire supérieure.
4. Déterminer les matrices de Jacobi BJ de Gauss-Seidel BGS .
5. On donne x(0) = (0, 0, 0)T , en déduire x(3) suivant la méthode de Jacobi et celle de
Gauss-Seidel.
6. Peut-on estimer la solution de (3) ? Si oui, vérifier le système (3) avec la solution
estimée x∗ .
7. On note rJ = ∥x∗ −x(3) ∥2 et rGS = ∥x∗ −x(3) ∥2 les erreurs commises respectivement
suivant les méthodes de Jacobi et de Gauss-Seidel.
a) Calculer rJ et rGS .
b) En déduire la meilleur des deux méthodes .