University of Aboubakr Belkaïd – Tlemcen

Faculté des Sciences de L’ingénieur

Département de Génie Civil

Résolution des systèmes

d’équations linéaires

Nom de l'enseignant 1
Système d’équation linéaire

E1 a11 x1 + a12 x2 ............a1n xn = b1

E2 a21 x1 + a22 x2 ............a2 n xn = b2

.
.
En an1 x1 + an 2 x2 ............ann xn = bn

Nom de l'enseignant 2
 Opérations sur un système
d’équation

„ Trois opérations permises:

(λEi) équivalente à (Ei)

(λEij + Ej) équivalent à Ei
Ej ↔ Ei

Nom de l'enseignant 3
 Exemple:
1° Etape
E1 : x1 + x 2 + ....... + 3x 4 = 4 (E 2 − 2E1 ) → E 2
(1)
E2 : 2x1 + x 2 − x 3 + x 4 = 1 (2) (E 3 − 3E1 ) → E 3
E3 : 3x1 − x 2 − x 3 + 2x 4 = 3 (E 4 + E1 ) → E 4
E4 : − x1 + 2x 2 + 3x 3 − x 4 = 4

E1 : x1 + x 2 + 3x 4 = +4
E2 : − x 2 − x 3 − 5x 4 = −7
E3 : − 4x 2 − x 3 − 7x 4 = −15
E4 : − 3x 2 + 3x 3 + 2x 4 = +8

Nom de l'enseignant 4
 Exemple:
2° Étape

(E 3 − 4E 2 ) → E 3
(E 4 + 3E 2 ) → E 4

E1 : x1 + x 2 + + 3x 4 = 4
E2 : − x 2 − x 3 − 5x 4 = −7
(3)
E3 : + 3x 3 − 13x 4 = 13
E4 : − 13x 4 = 8

Nom de l'enseignant 5
 Résultat Final
− 13
x4 = = −1
13
x 3 13 − 13x 4 ) = (13 − 13) = 0
1( 1
3 3
x 2 = −(− 7 + 5x 4 + x 3 ) = −(− 7 + 5 + 0) = +2
x1 = 4 − 3x 4 − x 2 = 4 − 3 − 2 = −1

(x1=-1, x2=2, x3=0, x4=1) est la solution unique du système (3) et on peut
alors vérifier que c’est aussi la solution du système d’équation (1) proposé.

Nom de l'enseignant 6
 Remarque importante

Inconnues x1, x2, x3 et x4 non utilisées

Elles ont gardés leurs places initiales

Plus commode d’utiliser le calcul Matricielle

Nom de l'enseignant 7
Méthode de résolution directe
d’élimination de Gauss

Triangularisé la matrice [A]

Résoudre le système équivalent

Nom de l'enseignant 8
 Qu’est ce qu’une matrice
triangulaire?
„ C’est une matrice dont les termes aij (i # j) du triangle
inférieur ou supérieur sont tous nuls.

⎡a11 a 12 a 1n ⎤ ⎡a 11 0 0 ⎤
⎢ ⎥ ⎢ ⎥
⎢0 a 22 a 2n ⎥ a
⎢ 21 a 22 0 ⎥
⎢............................. ⎥ ⎢............................. ⎥
⎢ ⎥ ⎢ ⎥
⎢............................. ⎥ ⎢............................. ⎥
⎢0 0.................a nn ⎥ ⎢a n1 a n2 ..............a nn ⎥
⎣ ⎦ ⎣ ⎦

Aij = 0 Aij = 0
si i<j matrice triangulaire supérieure si i<j matrice triangulaire supérieure

Nom de l'enseignant 9
 Triangularisation

Pour résoudre [A]{x}={b} on commence par triangularisé la

matrice [A]
Si a11<>0 pivot sur la ligne n° 1

((EJ-aJ1/a11Ei) Æ Ej) est conduite pour tous les j=2,3…..n

afin d’annuler les aj1

Nom de l'enseignant 10
Triangularisation
La matrice résultante sera alors de la forme.

⎡a11 a12 ..........

...a1n a1,n+1⎤
⎢ ⎥
⎢0 a22 ..........
...a2n a2,n+1 ⎥
⎢. ⎥
⎢ ⎥
⎢. ⎥
⎢ ⎥
⎣0......0 .ann an,n+1 ⎦

Nom de l'enseignant 11
 Résolution du système équivalent

a n, n +1
xn =
a nn
x n −1 = (a n −1, n +1 − a n −1, n x n )/a n −1, n −1

Nom de l'enseignant 12
 Résolution du système équivalent

a i , x +1 − a ix x x − a i , x −1 x x −1 − ...... − a i ,i +1 x i +1
xi =
a ii

x
ai , x +1 − aix x x − ∑a ij xj
j =i +1
xi =
aii

Tous les xi peuvent être définie pour : i= n-1, …………..1 connaissant xn.

Nom de l'enseignant 13
 Méthode d’élimination de Gauss

„ [A] triangulaire supérieure est obtenu par une série de

transformation [A(1)], [A(2)]…..A(k).
„ Avec
„ [A(0)] est la matrice initiale du système.
„

„ [A(1)] : matrice transformé obtenu de [A(0)] en annulant a1j

„

„ [A(k)] : matrice transformé obtenu de [A(k-1)] en annulant akj

Nom de l'enseignant 14
 Méthode d’élimination de Gauss
⎧
⎪ (k−1)
⎪aij pour i = 1, 2...k−1 et j = 1,2...n+1
(k) ⎪⎪
=
aij ⎨⎪0 pour i = k, k +1,....n et j = 1,2...k −1
(k−1)
⎪ (k−1) ai,k−1 (k−1)
⎪aij − (k−1) ak−1,j pour i = k, k +1......n et j = k,k +1....n +1
⎪⎩ ak−1,k−1

Nom de l'enseignant 15
Matrice à la transformation k-1

(0) () () () () ()
⎡a11 a120 ............a1,0k−1 a1k0 ..........a1n0 a1,0n+1 ⎤
⎢ (2) (2) (2) (2) (2)
⎥
⎢0 a22 ............a2,k−1 a2k ..........a2n a2,n+1 ⎥
⎢ ⎥
⎢..........
...................................................................... ..........
.... ⎥
[A~( )]
k-1
= ⎢0
⎢
( −)
akk−1,1k−1
( −) ( −)
akk−1,1k .....akk−1,1n akk−1,1k+1⎥
(−)
⎥
⎢........................................0 ( ) ( ) ( )
akkk ........aknk ak,kx+1 ⎥
⎢ ⎥
⎢.................................................................................................... ⎥
⎢ (k) (k ) (k) ⎥
⎣0......... ...............0.... ...........................a nk .........a nn a n,n +1 ⎦

Nom de l'enseignant 16
 Permutation d’équations
Si l’on reconsidère maintenant la transformation :

⎛ a (k )
⎞
⎜ Ei - , k E k ÷ → E i
i
⎜ (k ) ÷
⎝ a kk ⎠

Si akk=0 alors on doit chercher une valeur apk<>0 tel que k<p<n
Si p est trouvé alors on permute l’équation p avec l’équation k
Sinon il n’y pas de solution unique

Nom de l'enseignant 17
 Programmation de la méthode
d’élimination de Gauss
Etape 1: Lire [A] et le vecteur {B} pour i = 1 à n et j = 1 à n
Etape 2: Triangularisation de [A]
a., Si apk =0 avec p=k alors on
b., Chercher p = p+1 tel que apk <>0
Si p trouver et si p#k alors permutation
équation p avec équation k
sinon ‘pas de solution unique’; Halt

Nom de l'enseignant 18
 Programmation de la méthode
d’élimination de Gauss

Etape 2: Triangularisation de A
c., Calculer aij=aij-aik/akk*akj
pour k=1 to n-1
i=k+1 to n
j=k+1 to n+1

Nom de l'enseignant 19
 Programmation de la méthode
d’élimination de Gauss

Etape 3: Résolution du système

a., Si ann= 0 alors ‘Pas de solution unique’; Halt
Sinon xn = an,n+1/ann
b., xi=ai,n+1-Σaij*xj pour i = n-1 to 1

Nom de l'enseignant 20