TD Méthodes Numériques
Série 2
Système d’équations linéaires : Méthodes directes II.
Exercice 1 :
On pose :
−2 1 −1 1 𝑥1 1.5
2 0 4 −3 𝑥2
[ ][ ] = [ 4 ]
−4 −1 −12 9 𝑥3 −14
−2 1 1 −4 4 𝑥 −6.5
1. Trouver la décomposition LU de la matrice.
2. Résoudre le système avec la décomposition LU.
Solution :
−2 1 −1 1
𝐴=[ 2 0 4 −3 ]
−4 −1 −12 9
−2 1 1 −4
1 0 0 0 1 0 0 0
𝐸1 = [ 1 1 0 0] 𝐿1 = [−1 1 0 0]
−2 0 1 0 2 0 1 0
−1 0 0 1 1 0 0 1
−2 1 −1 1
𝐸1 𝐴 = [ 0 1 3 −2 ]
0 −3 −10 7
0 0 2 −5
1 0 0 0 1 0 0 0
𝐸2 = [0 1 0 0] 𝐿2 = [ 0 1 0 0]
0 3 1 0 0 −3 1 0
0 0 0 1 0 0 0 1
−2 1 −1 1
𝐸2 𝐸1 𝐴 = [ 0 1 3 −2]
0 0 −1 1
0 0 2 −5
FGMGM/USTHB Page 1
TD Méthodes Numériques
1 0 0 0 1 0 0 0
𝐸3 = [0 1 0 0] 𝐿3 = [0 1 0 0]
0 0 1 0 0 0 1 0
0 0 2 1 0 0 −2 1
−2 1 −1 1 1 0 0 0
𝑈 = 𝐸3 𝐸2 𝐸1 𝐴 = [ 0 1 3 −2] 𝐿 = 𝐿 𝐿 𝐿 = [−1 1 0
1 2 3
0]
0 0 −1 1 2 −3 1 0
0 0 0 −3 1 0 −2 1
1.5 1.5
𝐿𝑦 = 𝑏 = [ 4 ] ⇒ 𝑦 = [ 5.5 ]
−14 −0.5
−6.5 −9
1.5 −0.5
𝑈𝑥 = 𝑦 = [ 5.5 ] ⇒ 𝑥 = [ 1 ]
−0.5 3.5
−9 3
Exercice 2 :
Donner la factorisation de Cholesky des matrices suivantes :
4 0 12 −6
1 −2 0
𝐴1 = [−2 8 −6] et 𝐴2 = [ 0 1 2 1 ]
12 2 49 −4
0 −6 25 −6 1 −4 51
Solution :
𝑏11 0 0 𝑏11 𝑏21 𝑏31
𝑇 𝑇
𝐴1 = 𝐵𝐵 𝐵 = [𝑏21 𝑏22 0 ] 𝐵 =[ 0 𝑏22 𝑏32 ]
𝑏31 𝑏32 𝑏33 0 0 𝑏33
2
𝑏11 = 1 ⇒ 𝑏11 = 1
𝑏11 𝑏21 = −2 ⇒ 𝑏21 = −2
𝑏11 𝑏31 = 0 ⇒ 𝑏31 = 0
2 2
𝑏21 + 𝑏22 = 8 ⇒ 𝑏22 = 2
𝑏21 𝑏31 + 𝑏22 𝑏32 = −6 ⇒ 𝑏32 = −3
2 2 2
𝑏31 + 𝑏32 + 𝑏33 = 25 ⇒ 𝑏33 = 4
1 0 0
𝐵 = [−2 2 0]
0 −3 4
FGMGM/USTHB Page 2
TD Méthodes Numériques
On fait la même démarche pour la matrice 𝐴2 c’est-à-dire que :
𝑏11 0 0 0 𝑏11 𝑏21 𝑏31 𝑏41
𝑏 𝑏22 0 0 0 𝑏22 𝑏32 𝑏42
𝐴2 = 𝐵𝐵𝑇 𝐵 = [ 21 ] 𝐵𝑇 = [ ]
𝑏31 𝑏32 𝑏33 0 0 0 𝑏33 𝑏43
𝑏41 𝑏42 𝑏43 𝑏44 0 0 0 𝑏44
et on trouve :
2 0 0 0
𝐵=[ 0 1 0 0]
6 2 3 0
−3 1 4 5
Exercice 3 :
Résoudre les systèmes suivants avec la méthode TDMA.
4 0 0 0 8
1 2 0 2
𝐴 = [1 4 2 0 ] 𝐴𝑥 = [2] et 𝐵 = [2 1 −1] 𝐵𝑥 = [−1]
0 2 9 4 8
0 0 −4 5 3 0 −1 −1 1
Solution :
Pour le premier système nous avons :
𝑐 𝑦 8
𝛾2 = 𝑏1 = 0 𝑧2 = 𝑏1 = 4 = 2
1 1
𝑐 𝑎 𝑧 −𝑦
𝛾𝑖+1 = − 𝑎 𝛾 𝑖−𝑏 𝑧𝑖+1 = 𝑎𝑖 𝛾𝑖 −𝑏𝑖 𝑥𝑖 + 𝛾𝑖+1 𝑥𝑖+1 = 𝑧𝑖+1
𝑖 𝑖 𝑖 𝑖 𝑖 𝑖
𝑐2 2 𝑎2 𝑧2 −𝑦2 2−2
𝛾3 = − =− = 0.5 𝑧3 = = =0
𝑎2 𝛾2 −𝑏2 0−4 𝑎2 𝛾2 −𝑏2 0−4
𝑐3 4 𝑎 𝑧 −𝑦 0−8
𝛾4 = − 𝑎 = − 1−9 = 0.5 𝑧4 = 𝑎3 𝛾3 −𝑏3 = 1−9 = 1
3 𝛾3 −𝑏3 3 3 3
𝑎 𝑧 −𝑦 −4−3
𝑥4 = 𝑧5 = 𝑎4 𝛾4 −𝑏4 = −2−5 = 1
4 4 4
𝑥3 = 1 − 0.5 = 0.5
𝑥2 = 0 − 0.25 = −0.25
𝑥1 = 2 − 0 = 2
FGMGM/USTHB Page 3