Chapitre 03
Chapitre 03
Chapitre 03
3-1 Introduction :
Dans la pratique, les ingénieurs trouvent souvent des difficultés dans la résolution
d’un système d’équations qui modélisent les divers éléments considérés. Ainsi, la
détermination de :
Mélange de produits chimiques
Courant et tension des réseaux électriques
Solution optimale en programmation linéaire…
Certains de ces problèmes nécessitent la résolution d’un système de quelques
milliers d’équations.
Exemple :
5𝑥 + 3𝑦 + 𝑧 = 2
{2𝑥 + 4𝑦 + 3𝑧 = 1
3𝑥 + 6𝑦 + 7𝑧 = 9
Tel que :
𝑎𝑖𝑗 sont les coefficients du système
𝑥𝑗 sont les inconnus du système
𝑏𝑖 sont les composants de la deuxième partie (second member).
Le système d’équations linéaire peut être traduit sous forme matricielle :
𝑎11 𝑎12 … 𝑎1𝑛 𝑥1 𝑏1
𝑎21 𝑎22 … 𝑎2𝑛 𝑥2 𝑏
[… … … … ] […] = [ 2] ou 𝐴𝑥 = 𝑏
…
𝑎𝑚1 𝑎𝑚2 … 𝑎𝑚𝑛 𝑥𝑛 𝑏𝑚
Qu’on peut le noter aussi sous forme de matrice augmentée : [𝐴 | 𝑏]
C.TABET HELLEL 1
Chapitre3 : Résolution des systèmes d’équations linéaires
C.TABET HELLEL 2
Chapitre3 : Résolution des systèmes d’équations linéaires
𝑙11 0 0 𝑥1 𝑏1
(𝑏2 −𝑙21 𝑥1 )
[ 𝑙21 𝑙22 0 ] [𝑥2 ] = [𝑏2 ] ⇒ 𝑙21 𝑥1 + 𝑙22 𝑥2 = 𝑏2 ⇔ 𝑥2 = 𝑙22
𝑙31 𝑙32 𝑙33 𝑥3 𝑏3
(𝑏3−𝑙31 𝑥1−𝑙32 𝑥2 )
𝑙31 𝑥1 + 𝑙32 𝑥2 + 𝑙33 𝑥3 = 𝑏3 ⇔ 𝑥3 = 𝑙33
Donc la solution d’un système avec une matrice triangulaire inférieure est :
𝑏1
𝑥1 =
𝑙11
{ 1
𝑥𝑖 = (𝑏𝑖 − ∑𝑖−1
𝑗=1 𝑙𝑖𝑗 𝑥𝑗 ) 𝑖 = 2…𝑛
𝑙𝑖𝑖
C.TABET HELLEL 3
Chapitre3 : Résolution des systèmes d’équations linéaires
𝑏3
𝑢33 𝑥3 = 𝑏3 ⇔ 𝑥3 = 𝑢33
Donc la solution d’un système avec une matrice triangulaire supérieure est :
𝑏𝑛
𝑥𝑛 =
𝑢𝑛𝑛
{ 1
𝑥𝑖 = (𝑏𝑖 − ∑𝑛𝑗=𝑖+1 𝑢𝑖𝑗 𝑥𝑗 ) 𝑖 = 𝑛 − 1 … ,1
𝑢𝑖𝑖
C.TABET HELLEL 4
Chapitre3 : Résolution des systèmes d’équations linéaires
Exemple :
𝒙𝟏 + 𝒙𝟐 = 𝟐
(S) {
−𝒙𝟏 + 𝟐𝒙𝟐 = 𝟑
1 1
⇔ (−1 2
) (𝑥𝑥1 ) = (23) ⇔ A.x=b
2
1 1
On calcule le déterminant du (S) : det (−1 2
) = 2 – (-1) = 3≠0
⇒ (S) admet une unique solution
Etape 1 :
L1
1 1
L2 (−1 2
) (𝑥𝑥1 ) = (23)
2
L1
⇔ (𝑥1 ) = (01
L2← L1 + L2𝑥
1
3
) (𝑥𝑥1 ) = (25)
2 2
Etape 2 :
𝒙 + 𝒙𝟐 = 𝟐
{ 𝟏 ⇒ 𝒙𝟐 = 𝟓⁄𝟑
𝟑𝒙𝟐 = 𝟓
𝒙𝟏 + 𝒙𝟐 = 𝟐 𝒙𝟏 + 𝟓⁄𝟑 = 𝟐 𝒙𝟏 = 𝟏⁄𝟑
⇔ { ⇔ { ⇔ {
𝒙𝟐 = 𝟓⁄𝟑 𝒙𝟐 = 𝟓⁄𝟑 𝒙𝟐 = 𝟓⁄𝟑
C.TABET HELLEL 5
Chapitre3 : Résolution des systèmes d’équations linéaires
𝟏⁄
x= (𝟓 𝟑)
⁄𝟑
Exemple :
Soit le système linéaire suivant :
𝑥1 + 2𝑥2 − 3𝑥3 = −7 L1 1 2 −3 −7
{2𝑥1 − 3𝑥2 + 5𝑥3 = 18 ⟺ L2 (24 −3
1
5
−2
18
24
)
L3
4𝑥1 + 𝑥2 − 2𝑥3 = 24
L1
1 2 −3 −7 1 2 −3 −7
L2 ← L2- 2L1 2 −3 5 18 0 −7 11 32
⟺ (4 1 −2 24
) ⟺ ( 0 −7 10 52
)
L3 ←L3-4L1
C.TABET HELLEL 6
Chapitre3 : Résolution des systèmes d’équations linéaires
L1 1 2 −3 −7 L1←L1-3L3 1 2 0 −67
L2 0 −7 11 32 L2← L2+11L3 0 −7 0 252
⟺ (0 0 −1 20
) ⟺ (0 0 −1 20
)
L3 ←L3-L2 L3
L1←2L2 +7L1 7 0 0 35
0 −7 0 252
⟺ L2 ( 0 0 −1 20
)
L3
7𝑥1 = 35 𝒙𝟏 = 5
⟺ {−7𝑥2 = 252 ⟺ {𝒙𝟐 = −36
−𝑥3 = 20 𝒙𝟑 = −20
C.TABET HELLEL 7
Chapitre3 : Résolution des systèmes d’équations linéaires
Exemple :
Soit la matrice A suivante :
𝒙𝟏
𝟐 −𝟏 𝟎 −𝟐
( 𝟒 −𝟏 𝟐 ) (𝒙𝟐 ) = ( 𝟏𝟒 )
−𝟔 𝟐 𝟎 𝒙𝟑 𝟏𝟐
A= LU
𝑳𝒚 = 𝒃
Ax=b ⇒ LUx=b ⇒ {
𝑼𝒙 = 𝒚
𝒖𝟏𝟏 = 𝟐 , 𝒖𝟏𝟐 = −𝟏 , 𝒖𝟏𝟑 = 𝟎
𝒍𝟐𝟏 . 𝒖𝟏𝟏 = 𝟒 ⇔ 𝟐𝒍𝟐𝟏 = 𝟒 ⇔ 𝒍𝟐𝟏 = 𝟐
𝒍𝟐𝟏 . 𝒖𝟏𝟐 + 𝒖𝟐𝟐 = −𝟏 ⇔ 𝟐 ∗ (−𝟏) + 𝒖𝟐𝟐 = −𝟏 ⇔ 𝒖𝟐𝟐 = 𝟏
𝒍𝟐𝟏 . 𝒖𝟏𝟑 + 𝒖𝟐𝟑 = 𝟐 ( )
⇔ 𝟐 ∗ 𝟎 + 𝒖𝟐𝟑 = 𝟐 ⇔ 𝒖𝟐𝟑 = 𝟐
𝒍𝟑𝟏 . 𝒖𝟏𝟏 = −𝟔 ⇔ 𝟐𝒍𝟑𝟏 = −𝟔 ⇔ 𝒍𝟑𝟏 = −𝟑
𝒍𝟑𝟏 . 𝒖𝟏𝟐 + 𝒍𝟑𝟐 . 𝒖𝟐𝟐 = 𝟐 ⇔ ( −𝟑 ∗ (−𝟏) + 𝒍𝟑𝟐 (𝟏) = 𝟐
) ⇔ 𝒍𝟑𝟐 = −𝟏
𝒍𝟑𝟏 . 𝒖𝟏𝟑 + 𝒍𝟑𝟐 . 𝒖𝟐𝟑 + 𝒖𝟑𝟑 = 𝟎 ⇔ ( ) ( ) ( ) ( )
−𝟑 ∗ 𝟎 + −𝟏 ∗ 𝟐 + 𝒖𝟑𝟑 = 𝟎
⇔ 𝒖𝟑𝟑 = 𝟐
C.TABET HELLEL 8
Chapitre3 : Résolution des systèmes d’équations linéaires
𝒚𝟏
1 0 0 1 0 0 −𝟐
(𝑙21 1 0) ⇔ (2 1 0) (𝒚𝟐 ) = ( 𝟏𝟒 )
𝑙31 𝑙32 1 −3 −1 1 𝒚𝟑 𝟏𝟐
𝑦1 = −2 𝑦1 = −2
⇔ { 2𝑦1 + 𝑦2 = 14 ⇔ { 𝑦2 = 18
−3𝑦1 − 𝑦2 + 𝑦3 = 12 𝑦3 = 24
2- On remplace la matrice U par ses valeurs :
2𝑥3 = 24 𝒙𝟑 = 12
⇔ { 𝑥2 + 2𝑥3 = 18 ⇔ {𝒙𝟐 = −6
2𝑥1 − 𝑥2 = −2 𝒙𝟏 = −4
Théorème : Si A est une matrice symétrique définie positive alors il existe une
matrice triangulaire inférieure L, tel que : A= L.LT.
On a : Ax=b ⇒ L.LTx=b
𝐿𝑦 = 𝑏
Donc : { 𝑇
𝐿 𝑥=𝑦
Rappel :
Une matrice A est symétrique si et seulement si : A= AT ⇔ 𝑎𝑖𝑗 = 𝑎𝑗𝑖 (i,j=1……n )
Elle est définie positive si et seulement si : det A(k) > 0 ∀ 𝑘, 1 ≤ 𝑘 ≤ 𝑛 avec
A(k) 1 ≤ k≤ n les sous matrices principales.
C.TABET HELLEL 9
Chapitre3 : Résolution des systèmes d’équations linéaires
𝑖−1
Fin
2
𝑙𝑖𝑖 = √𝑎𝑖𝑖 − ∑ 𝑙𝑖𝑘 ,
𝑘=1
Exemple :
𝑥1 + 2𝑥2 − 3𝑥3 = −1
Soit le système suivant : { 2𝑥1 + 5𝑥2 − 4𝑥3 = −3
−3𝑥1 − 4𝑥2 + 14𝑥3 = 1
1 2 −3
−1
⇒ A= (−3
2 5
−4
−4 )
14
, b=(−3)
1
C.TABET HELLEL 10
Chapitre3 : Résolution des systèmes d’équations linéaires
2
𝑙11 𝑙11𝑙21 𝑙11 𝑙31
( 𝑙 𝑙 2
𝑙21 2
+ 𝑙22 𝑙21𝑙31 + 𝑙22 𝑙32 )
21 11
2 2 2
𝑙31 𝑙11 𝑙31𝑙21 + 𝑙32 𝑙22 𝑙31 + 𝑙32 + 𝑙33
1 0 0 1 2 −3
⇔ A= ( 2 1 0 ) . (0 1 2 )
−3 2 1 0 0 1
L LT
𝐿𝑦 = 𝑏
On a : {
𝐿𝑇 𝑥 = 𝑦
𝒚𝟏
1 0 0 −1
𝑳𝒚 = 𝒃 ⇒ (2 1 0 ) (𝒚𝟐 ) = (−3)
−3 2 1 𝒚𝟑 1
2
𝑙11 = 1 ⇒ 𝑙11 = 1
𝑙11𝑙21 = -1 ⇒ 𝑙21 = -1
𝑙11 𝑙31 = 2 ⇒ 𝑙31 = 2
2 2
𝑙21 + 𝑙22 = -1 ⇒ 𝑙22 = 1
𝑙21 𝑙31 + 𝑙22 𝑙32 = -4 ⇒ 𝑙32 = 2
2 2 2
𝑙31 + 𝑙32 + 𝑙33 = 14 ⇒ 𝑙33 = 1
𝑦1 = −1 𝑦1 = −1
{ 2𝑦1 + 𝑦2 = −3 ⇒ {𝑦2 = −1
−3𝑦1 + 2𝑦2 + 𝑦3 = 1 𝑦3 = 0
𝒙𝟏
1 2 −3 −1
𝑳 𝒙=𝒚 ⇒
𝑻
(0 1 2 ) (𝒙𝟐 ) = (−1)
0 0 1 𝒙𝟑 0
𝑥1 + 2𝑥2 − 3𝑥3 = −1 𝒙𝟏 = 𝟏
{ 𝑥2 + 2𝑥3 = −1 ⇒ {𝒙𝟐 = −𝟏
𝑥3 = 0 𝒙𝟑 = 𝟎
Les méthodes directes sont utilisées surtout pour la résolution des petits
(n<20) et moyens (n<100) systèmes d’équations linéaires.
On dit que la meilleure méthode directe au sens du temps de calcul est celle
qui demande le moins d’opérations.
C.TABET HELLEL 11
Chapitre3 : Résolution des systèmes d’équations linéaires
3.7.2 La convergence :
- pour vérifier la convergence il faut que : |𝑎𝑖𝑖 | > ∑𝑛𝑗=1 |𝑎𝑖𝑗 | i=1,2….n
𝑗≠𝑖
C.TABET HELLEL 12
Chapitre3 : Résolution des systèmes d’équations linéaires
C.TABET HELLEL 13
Chapitre3 : Résolution des systèmes d’équations linéaires
𝟏
𝒙𝒌+𝟏
𝟏 = (𝒃𝟏 − 𝒂𝟏𝟐 𝒙𝒌𝟐 − ⋯ − 𝒂𝟏𝒏 𝒙𝒌𝒏 )
𝒂𝟏𝟏
𝟏
⇔ 𝒙𝒌+𝟏
𝟐 = (𝒃𝟐 − 𝒂𝟐𝟏 𝒙𝒌𝟏 − ⋯ − 𝒂𝟐𝒏 𝒙𝒌𝒏 )
𝒂𝟐𝟐
… … … … …
𝟏
𝒌+𝟏
{𝒙𝒏 = (𝒃𝒏 − 𝒂𝒏𝟏 𝒙𝒌𝟏 − ⋯ − 𝒂𝒏𝒏−𝟏 𝒙𝒌𝒏−𝟏 )
𝒂𝒏𝒏
𝟏
⇔ 𝒙𝒌+𝟏
𝒊 = (𝒃𝒊 − ∑𝒏𝒋=𝟏 𝒂𝒊𝒋𝒙𝒌𝒋 ) i=1…n
𝒂𝒊𝒊
𝒋≠𝒊
Cette méthode suppose des pivots non nuls. Dans le cas où un pivot est nul, un
simple échange des lignes suffit pour vérifier cette condition.
Test d’arrêt : |𝒙𝒌+𝟏
𝒊 − 𝒙𝒌𝒊 | ≤ 𝜺 𝜺: 𝒖𝒏𝒆 𝒑𝒓é𝒄𝒊𝒔𝒊𝒐𝒏 𝒅𝒐𝒏𝒏é𝒆
Données : A, b, 𝒙(𝟎) , 𝜺
Début
k=0
Tant que |𝒙𝒌+𝟏
𝒊 − 𝒙𝒌𝒊 | > 𝜺 faire
Pour i=1 à n faire
𝟏
𝒙𝒌+𝟏
𝒊 = 𝒂 (𝒃𝒊 − ∑𝒏𝒋=𝟏 𝒂𝒊𝒋 𝒙𝒌𝒋 )
𝒊𝒊
Fin pour
k=k+1
Fin tant que
Fin
C.TABET HELLEL 14
Chapitre3 : Résolution des systèmes d’équations linéaires
Exemple :
Résoudre le système suivant par la méthode itérative de Jacobi avec 3 décimales
exactes.
𝜺= 0.5.10-3 et le vecteur initial x(0)= (2, 3, 5)t
4𝑥1 + 0.24𝑥2 − 0.08𝑥3 = 8
{ 0.09𝑥1 + 3𝑥2 − 0.15𝑥3 = 9
0.04𝑥1 − 0.08𝑥2 + 4𝑥3 = 20
Solution :
- Vérifions d’abord la convergence :
𝑎11 = 4 > 0.24 + |−0.08| = 0.32
𝑎22 = 3 > 0.09 + |−0.15| = 0.24
𝑎33 = 4 > 0.04 + |−0.08| = 0.12
Donc le processus itératif de Jacobi est convergent.
1
𝑥1 = (8 − 0.24𝑥2 + 0.08𝑥3 )
4
1
𝑥2 = (9 − 0.09𝑥1 + 0.15𝑥3 )
3
1
𝑥3 = (20 − 0.04𝑥1 + 0.08𝑥2 )
{ 4
𝑥1 = 2 − 0.06𝑥2 + 0.02𝑥3
⇔ {𝑥2 = 3 − 0.03𝑥1 + 0.05𝑥3
𝑥3 = 5 − 0.01𝑥1 + 0.02𝑥2
En remplaçant les valeurs du vecteur x(0), on obtient :
𝑥1 = 2 − 0.06(3) + 0.02 (5) 𝒙𝟏 = 𝟏. 𝟗𝟐
{ 𝑥2 = 3 − 0.03(2) + 0.05(5) ⇔ {𝒙𝟐 = 𝟑. 𝟏𝟗 ⇔
𝑥3 = 5 − 0.01(2) + 0.02(3) 𝒙𝟑 = 𝟓. 𝟎𝟒
C.TABET HELLEL 15
Chapitre3 : Résolution des systèmes d’équations linéaires
C.TABET HELLEL 16
Chapitre3 : Résolution des systèmes d’équations linéaires
𝟏
𝒙𝒌+𝟏
𝟏 = (𝒃𝟏 − 𝒂𝟏𝟐 𝒙𝒌𝟐 − ⋯ − 𝒂𝟏𝒏 𝒙𝒌𝒏 )
𝒂𝟏𝟏
𝟏
⇔ 𝒙𝒌+𝟏
𝟐 = (𝒃𝟐 − 𝒂𝟐𝟏 𝒙𝒌+𝟏
𝟏 − ⋯ − 𝒂𝟐𝒏 𝒙𝒌𝒏 )
𝒂𝟐𝟐
… … … … …
𝟏
𝒌+𝟏
{𝒙𝒏 = (𝒃𝒏 − 𝒂𝒏𝟏 𝒙𝒌+𝟏
𝟏 − ⋯ − 𝒂𝒏𝒏−𝟏 𝒙𝒌+𝟏
𝒏−𝟏 )
𝒂𝒏𝒏
𝟏
⇔ 𝒙𝒌+𝟏
𝒊 = [𝒃𝒊 − ∑𝒊−𝟏 𝒌+𝟏
𝒋=𝟏 𝒂𝒊𝒋 𝒙𝒋 − ∑𝒏𝒋=𝒊+𝟏 𝒂𝒊𝒋𝒙𝒌𝒋 ]
𝒂𝒊𝒊
Remarque :
On constate que la différence entre cette méthode et celle de Jacobi est
l’utilisation immédiate des nouveaux estimés à l’itération (k+1).
La condition d’avoir des pivots non nuls reste nécessaire aussi pour la méthode
de Gauss-Seidel.
3.7.5.2 Algorithme de Gauss Seidel :
Données : A, b, 𝒙(𝟎) , 𝜺
Début
k=0
Tant que |𝒙𝒌+𝟏
𝒊 − 𝒙𝒌𝒊 | > 𝜺 faire
Pour i=1 à n faire
𝟏
𝒙𝒌+𝟏
𝒊 = 𝒂 [𝒃𝒊 − ∑𝒊−𝟏 𝒌+𝟏
𝒋=𝟏 𝒂𝒊𝒋 𝒙𝒋 − ∑𝒏𝒋=𝒊+𝟏 𝒂𝒊𝒋 𝒙𝒌𝒋 ]
𝒊𝒊
Fin pour
k=k+1
Fin tant que
Fin
Quelques remarques :
- La méthode de Gauss-Seidel est plus rapide en convergence que la méthode
de Jacobi. En effet, on utilise dans la même itération les valeurs déjà
estimées. Dans un certain sens, le (k+1) ième itéré est plus proche de la
limite pour la méthode de Gauss-Seidel que pour la méthode de Jacobi.
C.TABET HELLEL 17
Chapitre3 : Résolution des systèmes d’équations linéaires
Exemple :
Résoudre le système suivant par la méthode itérative de Gauss-Seidel avec 4
décimales exactes. (𝜺= 0.5.10-4)
Avec le vecteur initial x(0)= (1.2, 0, 0)t
10𝑥1 + 𝑥2 + 𝑥3 = 12
{ 2𝑥1 + 10𝑥2 + 𝑥3 = 13
2𝑥1 + 2𝑥2 + 10𝑥3 = 14
Solution :
Vérifions d’abord la convergence :
𝑎11 = 10 > 1 + 1 = 2
𝑎22 = 10 > 2 + 1 = 3
𝑎33 = 10 > 2 + 2 = 4
Ainsi, le processus itératif de Gauss-Seidel est convergent.
1
𝑥1 = (12 − 𝑥2 − 𝑥3 )
10 𝑥1 = 1.2 − 0.1𝑥2 − 0.1𝑥3 )
1
𝑥2 =
10
(13 − 2𝑥1 − 𝑥3 ) ⇔ {𝑥2 = 1.3 − 0.2𝑥1 − 0.1𝑥3 )
1 𝑥3 = 1.4 − 0.2𝑥1 − 0.2𝑥2 )
𝑥3 = (14 − 2𝑥1 − 2𝑥2 )
{ 10
C.TABET HELLEL 18
Chapitre3 : Résolution des systèmes d’équations linéaires
(2)
𝑥1 = 1.2 − 0.1 ∗ 1.06 − 0.1 ∗ 0.948) = 0.9992
{ 𝑥2(2) = 1.3 − 0.2 ∗ 0.9992 − 0.1 ∗ 0.948) = 1.00536
(2)
𝑥3 = 1.4 − 0.2 ∗ 0.9992 − 0.2 ∗ 1.00536) = 0.999098
Pour les itérations restantes on a :
𝑥 (3) = (0.9996, 1.0002, 1.0001)𝑡
𝑥 (4) = (1.0000, 1.0000, 1.0000)𝑡
𝑥 (5) = (1.0000, 1.0000, 1.0000)𝑡
Ainsi, la solution approchée du système est :
𝒙𝟏 = 𝟏. 𝟎𝟎𝟎𝟎 ± 𝟎. 𝟎𝟎𝟎𝟏
𝒙𝟐 = 𝟏. 𝟎𝟎𝟎𝟎 ± 𝟎. 𝟎𝟎𝟎𝟏
𝒙𝟑 = 𝟏. 𝟎𝟎𝟎𝟎 ± 𝟎. 𝟎𝟎𝟎𝟏
C.TABET HELLEL 19
Chapitre3 : Résolution des systèmes d’équations linéaires
C.TABET HELLEL 20