EMI Rabat Année 2024-2025
1ère Année Méthodes Numériques
Genie Mineral
EXERCICE 1 :
Soient α, β, γ ∈ R. On considre le systme linaire suivant d’inconnues x1 , x2 , x3
:
x1 + 2x2 − 3x3 = α
2x1 + 6x2 − 5x3 = β (1)
x1 − 2x2 + 7x3 = γ
1. crire le systme (1) sous la forme Ax = b, avec A ∈ M∋ (R), x ∈ R3 , et
b ∈ R3 , que lon explicitera.
2. Est-ce que le systme (1) admet une unique solution pour tout α, β, γ ∈ R ?
3. Montrer que A admet une unique factorisation LU.
Dans la suite, on choisit α = 1, β = −1, et γ = 2 et on va rsoudre le systme
Ax = b de plusieurs faons :
a) Rsoudre le systme (1) par l’algorithme de Gauss sans pivot.
b) Calculer la factorisation LU de A puis rsoudre le systme (1) en utilisant
cette factorisation LU.
c) Rsoudre le systme (1) par l’algorithme de Gauss avec pivot partiel.
d) Calculer la factorisation L̃U = P A (o P est la matrice produit des ma-
trices de permutations effectues dans l’algorithme de Gauss avec pivot
partiel), puis rsoudre le systme (1) en utilisant cette factorisation.
EXERCICE 2 :
Il s’agit de compléter et d’exécuter un algorithme permettant de résoudre
un système linéaire A⃗x = ⃗b grâce à la méthode de Cholesky, la matrice A
étant tridiagonale, symétrique définie positive. Soit A la matrice tridiagonale,
symétrique et définie positive suivante :
a1 c 1 0 ... 0
c 1 a2 c2 . 0
A=
. .
. . . cN −1
0 cN −1 aN
1
Notons ⃗a le vecteur dont les composantes sont (ai )N c le (N − 1)-vecteur
i=1 et ⃗
N −1
donné de composantes (ci )i=1 .
Soit L la matrice triangulaire inférieure telle que A = LLT . Alors L aura la
forme suivante :
l1 0 . . . . 0
m1 l2 0 . 0
L=
. .
. . . 0
0 mN −1 lN
Soit ⃗l le N -vecteur de coefficients (li )N ⃗ le (N − 1)-vecteur de coefficients
i=1 , m
N −1
(mi )i=1 .
1o ) Déterminer les vecteurs ⃗l et m
⃗ en fonction de ⃗a et ⃗c et compléter l’algorithme
de Cholesky suivant (le vecteur ⃗l est mémorisé dans ⃗a et m ⃗ dans ⃗c) :
Décomposition de Cholesky :
√
a1 := a1
Faire i = 1 à (N − 1)
ci
ci :=
?√
ai+1 := ai+1 −?∗?
2o ) Résoudre le système linéaire triangulaire inférieur L⃗y = ⃗b
3o ) Déduire la solution du système A⃗x = ⃗b
4o ) Compléter l’algorithme suivant :
Entrées : ⃗a,⃗b, ⃗c.
Sorties : la solution ⃗x de A⃗x = ⃗b est stockée dans le vecteur ⃗b.
Résolution (d’abord L⃗y = ⃗b, puis LT ⃗x = ⃗y ) :
b1
b1 :=
a1
Faire i = 1 à (N − 1)
bi+1 −?∗?
bi+1 :=
?
bN
bN :=
aN
2
Faire i = (N − 1) à 1 (pas de -1)
bi −?∗?
bi :=
?
EXERCICE 3 :
1o ) utiliser une décomposition convenable pour résoudre le système AX = b
avec:
0 1 2
12
A=
1 2 3
b=
11
2 3 1 2
2o ) Faire deux itérations par la méthode de Jacobi et une itération de la méthode
de Gauss-Seidel Pour approcher X, avec :
0
X =
0
0
3o ) Peut - on affirmer la convergence de ces deux méthodes dans le cas du
système considéré
4o ) Ecrire l’algorithme de Gauss-Seidel.
EXERCICE 4:
Considrons le systme non-linaire suivant:
x21 − 10x1 + x22 = −8
x1 x22 + x1 − 10x2 = −8
1. Utiliser la méthode de Newton pour résoudre le systéme considéré, faire 3
itérations avec
0 0
X =
0
2. Utiliser la méthode de Gauss-Seidel en linéarisant par A(X k )X k+1 = b, faire
3 itrations.
3. Utiliser la méthode de Newton pour calculer une approximation de la solu-
tion du système non linaire suivant:
3x1 − cos(x2 x3 ) − 0.5 = 0
x21 2
− 81(x2 + 0.1) + sin(x3 ) + 1.06 = 0
e−x1 x2 + 20x3 + 10π−3
3 = 0
On peut considérer comme approximation initiale xt0 = (0.1; 0.1; −0.1).
3
EXERCICE
11 5:
On se propose d’approcher les valeurs propres et les vecteurs propres de la
matrice suivante:
2 −1/2 0
A= −1/2 2 0
0 0 2
1. Écrire trois itérations de la méthode de la puissance itérée. Comparer avec
les valeurs exactes.
2. Utiliser la méthode de Jacobi.
EXERCICE 6 :
Considérons l’équation non linéaire
x2 − 2x − 3 = 0
1o ) Utiliser la méthode de la Newton pour calculer une approximation de la
solution (faire 2 itérations).
3o ) Considérons la méthode du point fixe, utilisant les fonctions suivantes :
√ 3 x2 − 3
g1 (x) = 2x + 3, g2 (x) = et g3 (x) = ,
x−2 2
prendre x0 = 4.
