Année universitaire: 2024 − 2025
École Nationale Supérieure des Mines De
Rabat Analyse Numérique: Tronc commun
TD1: Résolution des Systèmes Linéaires et Non-linéaires
Exercice 1 :
Soit le système linéaire Ax = b suivant :
1 2 3 4 x1 10
2 3 4 1 x 10
2
=
3 4 1 2 x3 10
4 1 2 3 x4 10
1. Résoudre le systèmes linéaire par la méthode du pivot de Gauss.
2. Déduire le déterminant de A.
3. Donner les matrices de passage Mk vérifiant A(k+1) = Mk A(k) où A(k) est la matrice à
l’étape k avec A(1) = A.
4. Déduire les matrices L et U avec L une matrice triangulaire inférieur à diagonale unitaire
et A = LU .
5. Utiliser maintenant l’algorithme de décomposition LU pour retrouver les matrices L et U .
Exercice 2 :
Considérons la matrice symétrique défnie positive A et le vecteur b suivants :
9 −3 3 3
A = −3 10 −4 b = −4
3 −4 18 18
1. Calculer la factorisation la plus adéquate de A. Justifier votre choix.
2. Utiliser cette factorisation pour calculer la solution du système AX = b.
3. Les méthodes de Gauss-Seidel et Jacobi appliquées au système Ax = b seront-t-elles con-
vergentes si le point de départ x0 est bien choisi ?
On donne X (0) = [1, 1, 0]T .
4. Effectuer deux itérations de la méthode de Jacobi pour estimer la solution X.
ENSMR— RABAT 1ere année Tronc Commun I. MEDARHRI & [Link] 1
5. Effectuer deux itérations de la méthode de Gauss-Siedel pour estimer la solution X.
6. A chaque itération et pour chaque méthode, calculer l’écart entre la solution exacte X et
son approximée. Commenter le résultat..
Exercice 3 :
On considère le système linéaire A⃗x = ⃗b suivant:
" #" # " #
0, 780 0, 563 x1 0, 217
=
0, 913 0, 659 x2 0, 254
dont la solution exacte est ⃗x = [1, −1]T . En résolvant le système à l’aide d’une méthode
numérique, on a trouvé une solution approximative ⃗x∗ dont le résidu est donné par
" #
⃗b − A⃗x∗ = 0, 000001
0, 0
1. Sachant que A−1 ∞ = 1, 693 × 106 , est-ce que la solution approximative ⃗x∗ est accept-
able? Justifier votre réponse.
Exercice 4 :
Considérons le calcul numérique de la racine positive de la fonction :
f (x) = x3 − 3x − 2
1. Donner la formule itérative de Newton.
2. Calculer les solutions approchées xk , k = 1; 2; 3; 4, pour x0 = 2, 1.
Exercice 5 :
Soit la fonction f définie par f (x) = x3 − x2 + 8x − 8. Pour calculer ses racines on utilise 4
méthodes de point fixe différentes décrites par les fonctions d’itération suivantes:
8 − x3 1 3 1 1 4 2x3 − x2 + 8
φ1 (x) = −x3 + x2 − 7x + 8 φ2 (x) = φ3 (x) = − x + x2 + x + φ4 (x) =
8−x 10 10 5 5 3x2 − 2x + 8
Le tableau ci-dessous présente les suites des itérées obtenues par ces quatre méthodes.
Méth Méthode B Méthode C Méthode D
x0 0.5000000000000000 0.500000000000000 0.5000000000000000 0.5000000000000000
x1 0.9125000000000001 1.032258064516129 4.625000000000000 1.050000000000000
x2 0.9897857421875000 1.000235245684712 −101.9160156250000 0.9845143884892086
x3 0.9989578145726552 1.000000012299503 1.069697123778202 × 106 1.004312677086027
x4 0.9998955643403695 1.000000000000000 −1.224001861234915 × 1018 0.9987590594698483
x5 0.9999895542527895 1.000000000000000 1.833775789385161 × 1054 1.000353832012369
x6 0.9999989554034564 1.000000000000000 −6.166499545700052 × 10162 0.9998988463640411
1. Montrer que ℓ = 1 est l’unique racine réelle de f .
ENSMR— RABAT 1ere année Tronc Commun I. MEDARHRI & [Link] 2
2. Associer chaque méthode à sa fonction d’itération (justifier chaque réponse). Commenter
le résultat.
Exercice 6 :
Pour chaque valeur du paramètre s, on a une fonction :
gs (x) = s3 + x(1 − s) + (x − s2 )3
1. Montrer que pour chaque valeur de s, x = s2 est un point fixe de gs .
2. La méthode des points fixes converge vers le point fixe x = s2 pour quelles valeurs de
s? si la convergence est linéaire, donner le taux de convergence. Sinon donner l’ordre de
convergence.
Exercice 7 :
1. Obtenir tous les points fixes de la fonction :
g (x) = λx(1 − x)
où λ est un paramètre (λ ̸= 0).
2. Déterminer pour chaque point fixe trouvé en 1) les valeurs de λ pour lesquelles ces points
fixes sont attractifs.
3. Déterminer pour chaque point fixe trouvé en 1) la valeur de λ pour laquelle la convergence
de la méthode des points fixes sera quadratique.
Exercice 8 :
On considère la matrice:
0 3 1
A = 4 1 1
2 2 4
1. Montrer que la décomposition LU de la matrice obtenue en permutant les lignes 1 et 2 de
la matrice A s’écrit P A = LU , où P est une matrice de permutation. Donner P , L et U .
2. Résoudre Ax = b, où b = (1, 5, 6)T à l’aide de la décomposition P A = LU .
ENSMR— RABAT 1ere année Tronc Commun I. MEDARHRI & [Link] 3