M.
Yahyaoui Analyse numérique
Analyse numérique
Solutions aux exercices du chapitre 1
Résolution des systèmes linéaire
Exercice 1.1
Résoudre les systèmes suivants par la méthode de Gauss (sans permutation des lignes) et une arithmétique de 2
chiffres significatifs avec arrondissement :
(𝑎) 4𝑥1 + 𝑥2 + 2𝑥3 = 9 (𝑏) 2𝑥1 + 4𝑥2 − 𝑥3 = −5
2𝑥1 + 4𝑥2 − 𝑥3 = −5 𝑥1 + 𝑥2 − 3𝑥3 = −9
𝑥1 + 𝑥2 − 3𝑥3 = −9 4 𝑥1 + 𝑥2 + 2𝑥3 = 9
Les solutions sont (a) 𝑥1 = 1.1, 𝑥2 = -1.1, 𝑥3 = 2.9 et (b) 𝑥1 = 1.0, 𝑥2 = -1.0, 𝑥3 = 3.0
a. Elimination :
Elimination de 𝑎21 et 𝑎31 : les pivots sont 1⁄2 = 0.5 et 1⁄4 = 0.25 et on note que :
−9 − (0.25)(9) = −11.25 → −11
4 1 2 9 4 1 2 9
(2 4 −1|−5) → (0 3.5 −2 |−9.5)
1 1 −3 −9 0 0.75 −3.5 −11
Elimination de 𝑎32 : le pivot est 0.75/3.5 = 0.2143
4 1 2 9
→ (0 3.5 −2 |−9.5)
0 0 −3.1 −9
Remontée :
𝑥3 = 2.9
−9.5 − (−2)(2.9)
𝑥2 = = −1.1
3.5
9 − (2)(2.9) − (−1.1)
𝑥1 = = 1.1
3.5
1.1
𝑥⃗ = [−1.1]
2.9
b. Elimination :
2 4 −1 −5 2 4 −1 −5 2 4 −1 −5
(1 1 −3|−9) → (0 −1 −2.5|−6.5) → (0 −1 −2.5 |−6.5)
4 1 2 9 0 −7 4 19 0 0 22 65
Remontée :
65
𝑥3 = = 3.0
22
−6.5 − (−2.5)(3.0)
𝑥2 = = −1.0
−1
−5 − (−1)(3.0) − (4)(−1)
𝑥1 = = 1.0
3.5
1.0
𝑥⃗ = [−1.0]
3.0
1
M. Yahyaoui Analyse numérique
Exercice 1.2
Résoudre les mêmes systèmes de l’exercice précédent par la méthode de Gauss-Jordan
Les solutions sont (a) 𝑥1 = 1.1, 𝑥2 = -1.0, 𝑥3 = 2.9 et (b) 𝑥1 = 1.0, 𝑥2 = -0.7, 𝑥3 = 3.0.
On note que la solution exacte est 𝑥1 = 1, 𝑥2 = -1, 𝑥3 = 3 pour les deux systèmes.
Méthode de Gauss-Jordan :
4 1 2 9 4 1 2 9 4 0 2.6 12 4 0 0 4.5
(2 4 −1|−5) → (0 3.5 −2 |−9.5) → (0 3.5 −2 |−9.5) → (0 3.5 0 |−3.6)
1 1 −3 −9 0 0.75 −3.5 −11 0 0 −3.1 −9 0 0 −3.1 −9
1.1
→ 𝑥⃗ = [−1.0]
2.9
2 4 −1 −5 2 4 −1 −5 2 0 −11 −31 2 0 0 2.0
(1 1 −3|−9) → (0 −1 −2.5|−6.5) → (0 −1 −2.5 |−6.5) → (0 −1 0 | 0.7 )
4 1 2 9 0 −7 4 19 0 0 22 65 0 0 22 65
1.0
→ 𝑥⃗ = [−0.7]
3.0
Exercice 1.3
On donne la matrice :
2 −1 −1
𝐴 = (3 3 9)
3 3 5
1. Factoriser cette matrice en un produit 𝐿𝑈 par l’approche de Doolittle.
2. Résoudre le système (𝐴)𝑥⃗ = 𝑏⃗⃗ en se basant sur la décomposition effectuée sous (1) pour
1
𝑏⃗⃗ = [−1]
1
a. Décomposition LU :
2 −1 −1
𝐴 = (3 3 9)
3 3 5
- Première ligne de 𝑨 : on multiplie la 1ère ligne de 𝐿 avec les colonnes de 𝑈 pour obtenir la 1ère ligne de 𝑈 :
⇒ 𝑢11 = 2, 𝑢12 = −1, 𝑢13 = 1
- Première colonne de A : on multiplie les reste des lignes de 𝐿 avec la 1ère colonne de 𝑈 pour déterminer le
reste de éléments de la 1ère colonne de 𝐴 :
𝑙21 𝑢11 = 𝑎21 = 3 𝑙21 = 3⁄2
𝑙31 𝑢11 = 𝑎31 = 3 ⇒ 𝑙31 = 1
1 0 0 2 −1 1
𝐿 = (3⁄2 1 0) , 𝑈 = (0 𝑢22 𝑢23 )
3⁄2 1 1 0 0 𝑢33
- Deuxième ligne de 𝑨 : on effectue les multiplications qui donne le reste de la 2ème ligne de 𝐴 (on multiplie
la 2ème ligne de 𝐿 avec les différentes colonnes de 𝑈) :
−3⁄2 + 𝑢22 = 𝑎22 = 3 ⇒ 𝑢22 = 9⁄2
3⁄2 + 𝑢23 = 𝑎23 = 9 ⇒ 𝑢23 = 15⁄2
- Multiplication qui donne 𝑎33 :
3⁄2 + 𝑢23 + 𝑢33 = 5 ⇒ 𝑢33 = −4
2
M. Yahyaoui Analyse numérique
La solution finale est :
1 0 0 2 −1 1
𝐿 = (3⁄2 1 0) , 𝑈 = (0 9⁄2 15⁄2)
3⁄2 1 1 0 0 −4
b. Résolution :
(𝐴)𝑥⃗ = 𝑏⃗⃗ ⇔ (𝐿)(𝑈)𝑥⃗ = 𝑏⃗⃗
On pose :
(𝐿)𝑧⃗ = 𝑏⃗⃗
1 0 0 𝑧1 1 𝑧1 1
(3⁄2 1 0) [𝑧2 ] = [−1] ⇒ [𝑧2 ] = [−2.5]
3⁄2 1 1 𝑧3 1 𝑧3 2
(𝑈)𝑥⃗ = 𝑧⃗
2 −1 1 𝑥1 1 𝑥1 0.89
(0 9⁄2 15⁄2) [𝑥2 ] = [−2.5] ⇒ [𝑥2 ] = [ 0.28 ]
0 0 −4 𝑥3 2 𝑥3 −0.5
