U NIVERSIT É DE G HARDA ÏA - D ÉPARTEMENT ST M ÉTHODES NUM ÉRIQUES - 2eme ST

Solution TD – 03 Systèmes d’équation linéaires

Exercice –1 2 Pour que le système admette une solution unique le

déterminant de A doit être différent de zéro :
1 Le système peut être mis sous forme matricielle Ax = b : 
α 6= 2
det A 6= 0 ⇒

α 0 1

x1
 
2
 α 6= 3
 α 0 2α   x2  =  1  3 La solution du système :
2α 2 −1 x3 3  
| {z }| {z } | {z }  7 2 1 
A x b

 0 α −1/5 
 
det A1  5 α 4/5  2α − 2
on déduit le déterminant de A : x1 = = =
  det A (α − 2)(α − 3) (α − 2)(α − 3)
 α 0 1   
 
 = 2α(1 − 2α)  α 7 1 
det A =  α 0 2α 
 2α 2 −1

  −4 0 −1/5 
 
det A2  1 5 4/5  α+1
x2 = = =
2 det A (α − 2)(α − 3) (α − 2)(α − 3)
 
Pour que le système admette une solution unique le  α 2 7 
 
déterminant de A doit être différent de zéro :  −4 α 0 
 
det A3  1 α 5  5α2 − 35α + 40

α 6= 0 x3 = = =
det A 6= 0 ⇒ det A (α − 2)(α − 3) (α − 2)(α − 3)
α 6= 1/2 Le deuxième système :
1 Le système sous forme matricielle Ax = b :
3 Selon la méthode de Cramer les solutions du systèmes     
det Ai α 0 1 x1 7
sont données par xi = :
det A  α 0 α   x2  =  0 
  α 2 −1 x3 5
 2 0 1 

 1 0 2α  on déduit le déterminant de A :
   
det A1  3 2 −1  2 − 8α 1 − 4α  α 0 1 
x1 = =  = =
 
 det A =  α 0 α  = 2α(1 − α)
det A  α 0 1  2α ( 1 − 2α ) α ( 1 − 2α) 

 α 0 2α 
  α 2 −1 
 
 2α 2 −1 
2 Pour que le système admette une solution unique le
  déterminant de A doit être différent de zéro :
 α 2 1  
  α 6= 0
 α 1 2α  det A 6= 0 ⇒
  α 6= 1
det A2  2α 3 −1  2α2 + 2α 1+α
x2 = =   = =
det A 
 α 0 1 
 2α ( 1 − 2α ) 1 − 2α 3 La solution du système :
 α 0 2α   
   7 2 1 
 2α 2 −1  
 0 α −1/5 
   
 α 0 2  det A1  5 α 4/5  −14α 7

0 1
 x1 = = = =
 α


 det A 2α(1 − α) 2α(1 − α) α−1
det A3  2α 2 3  2α 1  
x3 = =   = =  α 7 1 
det A  α
 0 1 
 2α ( 1 − 2α ) 1 − 2α 
 −4 0 −1/5 
 α 0 2α  
2α2 + 12α
  1 5
  det A2 4/5  6+α
 2α 2 −1  x2 = = = =
det A 2α(1 − α) 2α(1 − α) 1−α
 
 α 2 7 
Exercice –2 
 −4 α 0 



det A3  1 α 5  14α 7
Le premier système : x3 = = = =
det A 2α(1 − α) 2α(1 − α) 1−α
1 Le système sous forme matricielle Ax = b :

α 2 1

x1
 
7
 Exercice –3
 −4 α −1/5   x2  =  0  1 Le système sous forme matricielle Ax = b :
1 α 4/5 x3 5     
2 −1 3 −1 x1 −4
on déduit le déterminant de A :  −4 4 − 3 0   x2   10 
  = 
   −2 5 4 −4   x3   8 
 α 2 1  −2 3 −6 4 x4 −4

det A =  −4 α −1/5  = (α − 2)(α − 3)
 1 α 4/5  2 Triangularisation du système par la méthode de Gauss :

Page 1/5 H. Bouderba

U NIVERSIT É DE G HARDA ÏA - D ÉPARTEMENT ST M ÉTHODES NUM ÉRIQUES - 2eme ST

– Etape 1 – Etape 2
 
L1 ← L1 2 −1 3 −1 −4 L1 ← L1 2

−1 1 0

L2 ← L2 + (2) L1 
 0 2 3 −2 2 
 L2 ← L2  0 −2 1 3 
L3 ← L3 + (1) L1  0 4 7 −5 4  L3 ← L3 + (−1/2) L2 0 0 1/2 1/2
L4 ← L4 + (1) L1 0 −4 −3 3 −8
ce qui correspond au système triangulaire Ãx = b̃ :
– Etape 2
    
  1 2 1 x1 0
L1 ← L1 2 −1 3 −1 −4  0 −2 1   x2  =  3 
L2 ← L2  0
 2 3 −2 2 
 0 0 1/2 x3 1/2
L3 ← L3 + (−2) L2  0 0 1 −1 0 
L4 ← L4 + ( 2) L2 0 0 3 −1 −4
3 A partir du système Ãx = b̃ :
– Etape 3 
 x1 + 2x2 + x3 =0
  − 2 x2 + x3 =3
L1 ← L1 2 −1 3 −1 −4 1/2 x3 = 1/2

L2 ← L2  0
 2 3 −2 2 

L3 ← L3  0 0 1 −1 0 
L4 ← L4 + (−3) L3 qu’on résout par remontée triangulaire :
0 0 0 2 −4

x = −2x2 − x3 = −2(−1) − (1) = 1
ce qui correspond au système triangulaire Ãx = b̃ :  1


3 − x3 3 − (1)
x2 = = = −1
− −2
    
2 −1 3 −1 x1 −4 

 x =1 2
 0 2 3 −2   x2  
    2  3
 = 
 0 0 1 −1   x3  0 
0 0 0 2 x4 −4 4 On déduit le déterminant de A :

3 En partant du système Ãx = b̃ on obtient : det A = det à = 1 × (−2) × 1/2 = −1



 2 x1 − x2 + 3 x3 − x4 = −4 Pour le deuxième système :
2 x2 + 3 x3 − 2 x4 = 2


 x3 − x4 = 0 
x1 − 2 x2 − 2 x3 − 2 x4 = −5
−4

2 x4 =
 
−3 x1 + 7 x2 + 9 x3 + 7 x4 = 20


qu’on peut résoudre par remontée triangulaire : 

 x1 − 4 x2 − 7 x3 − 6 x4 = −16
x1 + x2 + 5 x3 + 6 x4 = 13

−4 + x2 − 3x3 + x4 −4 + 2 − 3(−2) + (−2)

 x1 = = =1


 2 2 1 La forme matricielle Ax = b :
2 − 3x3 + 2x4 2 − 3(−2) + 2(−2)



x2 = = =2
2 2 
1 −2 −2 −2

x1

−5



 x 3 = x 4 = − 2
  −3 7 9 7   x2 20 
 x = −4 = −2

  
   = 
4
2
 1 −4 −7 −6   x3   −16 
1 1 5 6 x4 13
4 On déduit le déterminant de A :
2 La triangularisation par la méthode de Gauss
det A = det à = 2 × 2 × 1 × 2 = 8
– Etape 1

Exercice –4 1

−2 −2 −2 −5

L1 ← L1
On obtient pour le système : L2 ← L2 + ( 3) L1 
 0 1 3 1 5 

 L3 ← L3 + (−1) L1  0 −2 −5 −4 −11 
 x1 + 2 x2 + x3 = 0 L4 ← L4 + (−1) L1 0 3 7 8 18
2 x1 + 2 x2 + 3 x3 = 3
− x1 − 3 x2 = 2

– Etape 2

1 l’écriture sous forme matricielle Ax = b :  

L1 ← L1 1 −2 −2 −2 −5

1 2 1

x1
 
0

L2 ← L2  0 1 3 1 5 
 
 2 2 3   x2  =  3  L3 ← L3 + ( 2) L2  0 0 1 −2 −1 
L4 ← L4 + (−3) L2 0 0 −2 5 3
−1 −3 0 x3 2

2 triangularisation par la méthode de Gauss : – Etape 3

– Etape 1
 
  L1 ← L1 1 −2 −2 −2 −5
L1 ← L1 1 2 1 0 L2 ← L2  0
 1 3 1 5 

L2 ← L2 + (−2) L1  0 −2 1 3  L3 ← L3  0 0 1 −2 −1 
L3 ← L3 + ( 1) L1 0 −1 1 2 L4 ← L4 + (2) L3 0 0 0 1 1

Page 2/5 H. Bouderba

U NIVERSIT É DE G HARDA ÏA - D ÉPARTEMENT ST M ÉTHODES NUM ÉRIQUES - 2eme ST

      
ce qui correspond au système triangulaire Ãx = b̃ : 1 2 1 x1 0 −9
 0 −2 1  x2  =  0  ⇒ v3 =  −3 
−1
    
1 −2 −2 −2 x1 −5 0 0 x3 1 1
  x2  = 
 0 1 3 1     5  | {z } | {z }
  v3 b̃3
 0 0 1 −2   x3   −1 
0 0 0 1 x4 1 donc :  
−12 20 −9
3 Le système matricielle Ãx = b̃ nous donne : A −1 =  −5 7 −3 
1 −2 1

x − 2 x2 − 2 x3 − 2 x4 = −5
 1 2 De même pour inverser la matrice :


x2 + 3 x3 + x4 = 5
− −1
 

 x3 2 x4 = 1 −1 −1
x4 = 1 A =  2 −1 −4 

−3 1 8
qu’on peut résoudre par remontée triangulaire :
 on construit le système augmenté :
x = −5 + 2x2 + 2x3 + 2x4 = −5 + 2(1) + 2(1) + 2(1) = 1
 1
  
1 −1 −1 1 0 0

x2 = 5 − 3x3 − x4 = 5 − 3(1) − 1 = 1
 2 −1 −4 0 1 0 
x = −1 + 2x4 = −1 + 2(1) = 1
 3

−3 1 8 0 0 1

x4 = 1
on le triangularise par la méthode de Gauss :
4 On déduit le déterminant de A :
– Etape 1
det A = det à = 1 × 1 × 1 × 1 = 1  
L1 ← L1 1 −1 −1 1 0 0
L2 ← L2 + (−2) L1  0 1 −2 −2 1 0 
L3 ← L3 + (3) L1 0 −2 5 3 0 1
Exercice –5
1 Pour inverser la matrice : – Etape 2
 
1 −2 3
 
L1 ← L1 1 −1 −1 1 0 0
A= 2 −3 9  L2 ← L2  0 1 −2 −2 1 0 
3 −4 16 L3 ← L3 + (2) L2 0 0 1 −1 2 1

on construit le système augmenté : la remontée triangulaire des trois systèmes nous donne :
 

1 −2 3 1 0 0
 −4 7 3
 2 −3 9 0 1 0  A −1 =  −4 5 2 
3 −4 16 0 0 1 −1 2 1

on le triangularise par la méthode de Gauss : 3 Pour la troisième matrice :

– Etape 1  
2 −1 2
  A =  1 −1 1 
L1 ← L1 1 −2 3 1 0 0 −1 1 1
L2 ← L2 + (−2) L1  0 1 3 −2 1 0 
L3 ← L3 + (−3) L1 0 2 7 −3 0 1 on construit le système augmenté :
 
2 −1 2 1 0 0
– Etape 2  1 −1 1 0 1 0 
  −1 1 1 0 0 1
1 −2 3 1 0 0
L1 ← L1  0 1 3 −2 1 0  on le triangularise par la méthode de Gauss :
L2 ← L2 0 0 1 1 −2 1
L3 ← L3 + (−2) L2 – Etape 1
b̃1 b̃2 b̃3
 
L1 ← L1 2 −1 2 1 0 0
ce qui correspond à trois systèmes triangulaires Ãvi = b̃i dont L2 ← L2 + (−1/2) L1  0 −1/2 0 −1/2 1 0 
la résolution par remontée triangulaire donne les vectreurs co- L3 ← L3 + (1/2) L1 0 1/2 2 1/2 0 1
lonnes v1 , v2 et v3 de la matrice inverse A−1 :
       – Etape 2
1 2 1 x1 1 −12  
 0 −2 1   x2  =  −2  ⇒ v1 =  −5  L1 ← L1 2 −1 2 1 0 0
0 0 −1 x3 1 1 L2 ← L2  0 −1/2 0 −1/2 1 0 
| {z } | {z } L3 ← L3 + (1) L2 0 2 0 0 1
v1 b̃1
la remontée triangulaire des trois systèmes nous donne :
      
1 2 1 x1 0 20  
 0 −2 1  x2  =  1  ⇒ v2 =  7  1 −3/2 −1/2
0 0 −1 x3 −2 −2 A −1 =  1 −2 0 
| {z } | {z } 0 1/2 1/2
v2 b̃2

Page 3/5 H. Bouderba

U NIVERSIT É DE G HARDA ÏA - D ÉPARTEMENT ST M ÉTHODES NUM ÉRIQUES - 2eme ST

Exercice –6 la multiplication de la matrice U par le vecteur x donne un

vecteur y et on arrive au système :
1 Le système :     
1 1 2 x1 y1

x1 + x2 + 2 x3 = 1  0 4 6   x2  =  y2  (2)

x1 + 5 x2 + 8 x3 = 3 0 0 1 x3 y3
2 x1 + 8 x2 + 14 x3 = 6
 | {z } | {z } | {z }
U x y

peut être mis sous forme matricielle Ax = b : on peut remplacer le produit Ux par le vecteur y dans
l’équation (1) qu’on résout par descente triangulaire :
    
1 1 2 x1 1 
1 0 0

y1
 
1
 
y1
 
1

 1 5 8   x2  =  3   1 1 0   y2  =  3  ⇒  y2  =  2  (3)
2 8 14 x3 6 2 3/2 1 y3 6 y3 1
| {z } | {z } | {z }
on procède à la décomposition de LU de A comme suit : L y b

on remplace les valeurs du vecteur y dans l’équation (2). On ar-

     rive à un système qu’on résout par remontée triangulaire pour
1 0 0 u11 u12 u13 1 1 2
 l21 1 obtenir le vecteur x :
0  0 u22 u23  =  1 5 8 
l31 l32 1 0 0 u33 2 8 14
        
1 1 2 x1 1 x1 0
| {z }| {z } | {z }  0 4 6   x2  =  2  ⇒  x2  =  −1 
L U A
0 0 1 x3 1 x3 1
On obtient les étapes suivantes : | {z } | {z } | {z }
x y
U
– Etape 1 On multiplie la ligne 1 de L par U on obtient :
3 La matrice inverse A−1 de la matrice A est une matrice
    carrée 3 × 3 elle a donc trois colonnes qu’on note v1 , v2 et v3 . Si
 u11 = 1 1 0 0 1 1 2
on note les colonnes de la matrice indentité I3 par e1 , e2 et e3 ,
u12 = 1 ⇒ A =  l21 1 0  0 u22 u23 
alors les colonnes v1 , v2 et v3 sont respectivement la solution
u13 = 2 l31 l32 1 0 0 u33

| {z }| {z } des trois systèmes : 
L U  A v 1 = e1
A v 2 = e2
– Etape 2 On multiplie la ligne 2 de L par U on obtient : 
A v 3 = e3
  qu’on peut résoudre par la décomposition LU obtenue
 l21 = 1  l21 = 1 précédement de la matrice A. Le premier système s’écrit :
l21 + u22 = 5 ⇒ u =4
 22
     
1 0 0 1 1 2 x1 1
2 l21 + u23 = 8 u23 = 6

 1 1 0   0 4 6   x2  =  0 
   2 /2 1
3 0 0 1 x3 0
1 0 0 1 1 2 | {z }| {z } | {z } | {z }
⇒A= 1 1 0  0 4 6  L U v1 e1
l31 l32 1 0 0 u33
| {z }
y
| {z }| {z }
L U la solution du système Ly = e1 par déscente triangulaire :
        
– Etape 3 On multiplie la ligne 3 de L par U on obtient : 1 0 0 y1 1 y1 1
 1 1 0   y2  =  0  ⇒  y2  =  −1 
  2 3/2 1 y3 0 y3 −1/2
 l31 = 2  l31 = 2 | {z } | {z } | {z }
l31 + 4 l32 = 8 ⇒ l = 3/2 L y e1
 32 (4)
2 l31 + 6 l32 + u33 = 14 u33 = 1

il reste à trouver v1 par remontée triangulaire :
  
1 0 0 1 1 2 
1 1 2

x1
 
1
 
3/2

⇒ A= 1 1 0  0 4 6   0 4 6  x2  =  −1  ⇒ v1 =  1/2 
2 3/2 1 0 0 1 0 0 1 x3 −1/2 −1/2
| {z }| {z } | {z }| {z } | {z }
L U v1 y
U

de même, la solution du système Ly = e2 par déscente trian-

2 Le système Ax = b :
gulaire :
             
1 1 2 x1 1 1 0 0 y1 0 y1 0
 1 5 8   x2  =  3   1 1 0   y2  =  1  ⇒  y2  =  1 
2 8 14 x3 6 2 3/2 1 y3 0 y3 −3/2
| {z } | {z } | {z }
L y e1
peut donc s’écrire :
(5)
      il reste à trouver v2 par remontée triangulaire :
1 0 0 1 1 2 x1 1
 1 1 0   0 4 6  x2  =  3  (1) 
1 1 2

x1
 
0
 
1/2

2 3/2 1 0 0 1 x3 6  0 4 6  x2  =  1  ⇒ v2 =  5/2 
| {z }| {z }| {z } | {z }
L U x b 0 0 1 x3 −3/2 −3/2
| {z } | {z }| {z } | {z }
y U v2 y

Page 4/5 H. Bouderba

U NIVERSIT É DE G HARDA ÏA - D ÉPARTEMENT ST M ÉTHODES NUM ÉRIQUES - 2eme ST

enfin, la solution du système Ly = e3 par déscente triangu- 3 Par remontée triangulaire on déduit :
laire :
         α−1 3
1 0 0 y1 0 y1 0 x1 = x2 = et x3 =
7α − 4 7α − 4
 1 1 0   y2  =  0  ⇒  y2  =  0  (6)
2 3/2 1 y3 1 y3 1
| {z
L
} | {z } | {z }
y e1
Exercice –9
reste à trouver v3 par remontée triangulaire : Inverser la matrice :
 
       1 α 2α
1 1 2 x1 0 −1/2
A= 2 2α + 1 3α 
 0 4 6  x2  =  0  ⇒ v3 =  −3/2 
−1 −α − 2 α
0 0 1 x3 1 1
| {z }| {z } | {z }
U v3 y on utilise la méthode de Gauss. On construit le système aug-
menté :
donc :  
  1 α 2α 1 0 0
3/2 1/2 −1/2  2 2α + 1 3α 0 1 0 
A −1 =  1/2 5/2 −3/2  −1 −α − 2 α 0 0 1
−1/2 −3/2 1
– Etape 1

Exercice –7  
L1 ← L1 1 α 2α 1 0 0
Voir la solution de l’exercice (5). L2 ← L2 + (−2) L1  0 1 −α −2 1 0 
L3 ← L3 + ( 1) L1 0 −2 3α 1 0 1
Exercice –8
1 Le système sous forme matricielle Ax = b : – Etape 2
    
1 3 α x1 1  
 2 5 1   x2  =  1  L1 ← L1 1 α 2α 1 0 0
L2 ← L2  0 1 −α −2 1 0 
2 2 α x3 1 L3 ← L3 + (2) L2 0 0 α −3 2 1
2 Triangularisation du système par la méthode de Gauss :
la remontée triangulaire donne les vectreurs colonnes v1 , v2 et
– Etape 1 v3 de la matrice inverse A−1 :
 
1 3 1
      
L1 ← L1 α 1 α 2α x1 1 5α + 7
L2 ← L2 + (−2) L1  0 −1 1 − 2α −1   0 1 −α   x2  =  −2  ⇒ v1 =  −5 
L3 ← L3 + (−2) L1 0 −4 −α −1 0 0 α x3 −3 −3/α
| {z } | {z }
– Etape 2 v1 b̃1
 
L1 ← L1 1 3 1
      
α 1 α 2α x1 0 −3α − 4
L2 ← L2  0 −1 1 − 2α −1   0 1 −α   x2  =  1  ⇒ v2 =  3 
L3 ← L3 + (−4) L2 0 0 7α − 4 3 0 0 α x3 2 2/α
| {z } | {z }
avec les coéfficients obtenus dans les étapes de la méthode de v2 b̃2
Gauss on déduit la décompostion LU suivante :
      
   1 α 2α x1 0 −α − 2
1 0 0 1 3 α  0 1 − α   x2  =  0  ⇒ v3 =  1 
 2 1 0   0 −1 1 − 2α 
0 0 α x3 1 1/α
2 4 1 0 0 7α − 4 | {z } | {z }
v3 b̃3
3 On déduit le déterminant de A :
donc :
det A = det U = 1 × (−1) × −4 + 7α = 4 − 7α
 
5α + 7 −3α − 4 −α − 2
pour accepter une solution unique il faut et il suffit que :
A −1 =  −5 3 1 
det A 6= 0 ⇒ α 6= 4/7 −3/α 2/α 1/α

Page 5/5 H. Bouderba