Résolution des systèmes linéaires : méthodes directes- Part 1
I) Les matrices réelles
On note par A = (aij )1 i n; 1 j m une matrice à éléments dans R ayant
n lignes et m colonnes.
On note par Mn;m (R) l’ensemble des matrices réelles du type (n; m) ou
de taille n m formées de n lignes et de m colonnes.
La matrice nulle de Mn;m (R) est designée par 0n m .
Si n = m, alors on note simplement Mn (R) : c’est l’ensemble des matrices
carrées d’ordre n. La matrice unité (ou identité) de Mn (R) est designée par
In .
Dé…nitions
Soit A = (aij )1 i; j n 2 Mn (R)
1) La diagonale principale de A est faii = i = 1; :::; ng:
P
n
2) La trace tr de A est tr(A) = aii :
i=1
3) A
0 est diagonale1 si aij = 0 pour tous i 6= j:
1 0 0 0
B 0 2 0 0 C
B C
@ 0 0 3 0 A
0 0 0 4
4) A est
0 triangulaire 1supérieure si aij = 0 pour tous j < i:
1 7 8 1
B 0 2 6 1 C
B C
@ 0 0 3 5 A
0 0 0 4
5) A est
0 triangulaire inférieure
1 si aij = 0 pour tous j > i:
3 0 0 0
B 4 1 0 0 C
B C
@ 5 4 1 0 A
2 3 9 3
6) A est symétrique si aij = aji pour tous i; j:
7) La transposée de A est notée t A = (a0ij )1 i;j n où a0ij = aji pour tous
i; j.
8) La matrice A est inversible s’il existe A0 2 Mn (R) telle que AA0 =
A0 A = In . Si A0 existe, on écrit alors A0 = A 1 .
Exemple de transposée
1
0 1 0 1
1 5 3 1 1 0
A = @ 1 6 2 A ! tA = @ 5 6 4 A
0 4 7 3 2 7
(les lignes de A sont respectivement les colonnes de t A).
On a toujours t (t A) = A:
Propriétés
Soient A; B 2 Mn (R)
1) La matrice A est symétrique si et seulement si t A = A
2) t (A + B) = t A + t B et t (AB) = t B t A
4) Si A 1 existe alors t (A 1 ) = (t A) 1
5) Si A 1 et B 1 existent alors (AB) 1 = B 1 A 1
6) Si A et B sont triangulaires supérieures (resp. inférieures) alors AB
est triangulaire
0 supérieure (resp.
1 0inférieure). 1 0 1
1 7 8 1 5 7 4 1 5 14 38 2
B 0 2 6 1 C B 0 1 6 1 C B 4 C
exemple : B C B C=B 0 2 6 C
@ 0 0 3 5 A @ 0 0 1 2 A @ 0 0 3 24 A
0 0 0 4 0 0 0 6 0 0 0 24
1
7) Supposons que A existe. Si A est trinagulaire supérieure (resp. tri-
nagulaire inférieure) alors A 1 est trinagulaire supérieure (resp. inférieure).
II) Les déterminants
La notion du déterminant concerne seulement les matrices carrées.
Soit A = (aij )1 i; j n 2 Mn (R), le déterminant de A est le nombre réel
a b a b
noté det(A) ou jAj; il est bien connu que det = = ad cb.
c d c d
On désigne par Ai; j la matrice d’ordre n 1, obtenue à partir de A en
supprimant la i eme ligne et la j eme colonne de A.
Le déterminant de A est donné par :
P
n
1) Pour tout j 2 f1; :::; ng, on a det(A) = ( 1)i+j aij det(Ai; j ) : c’est
i=1
le développement de det(A) selon la j eme colonne de A.
P
n
2) Pour tout i 2 f1; :::; ng, on a det(A) = ( 1)i+j aij det(Ai; j ) : c’est
j=1
le développement de det(A) selon la i eme ligne de A.
Noter que la valeur de det(A) est indépendante du choix de i ou de j.
Exemples
1)
2
0 1
1 3 4
A=@ 2 2 1 A
2 3 2
On développe par rapport à la 1ere ligne
1 3 4
2 1 2 1
det A = 2 2 1 = 1( 1)1+1 + 3( 1)1+2 +
3 2 2 2
2 3 2
2 2
4( 1)1+3
2 3
det A = 7 6 40 = 53
2) 0 1
1 3 4 2
B 1 0 2 4 C
A=B @ 2 4 3
C
3 A
3 0 5 2
On développe par rapport à la 2ème colonne
1 3 4 2
1 2 4 1 4 2
1 0 2 4 1+2 2+2
det A = = 3( 1) 2 3 3 +0( 1) 2 3 3
2 4 3 3
3 5 2 3 5 2
3 0 5 2
1 4 2 1 4 2
3+2 4+2
+4( 1) 1 2 4 + 0( 1) 1 2 4
3 5 2 2 3 3
1 2 4 1 4 2
1+2 3+2
det A = 3( 1) 2 3 3 + 4( 1) 1 2 4
3 5 2 3 5 2
il faut encore calculer les 2 déterminants 3x3 suivant une colonne ou une
ligne de chacun d’eux.
Remarque : choisir la ligne ou la colonne qui contient le plus de 0.
Propriétés
Soient A = (aij )1 i; j n 2 Mn (R) et 2 R
1) La valeur de det(A) reste inchangée si on ajoute à une ligne (resp.
colonne) de A une combinaison linéaire formée de lignes (resp. colonnes)
restantes de A.
2) La valeur de det(A) est égale au déterminant de la matrice obtenue
à partir de A en multipliant l’une des lignes ou colonnes de A par .
3) Si on permute 2 lignes ou 2 colonnes de A alors le déterminant de la
matrice obtenue est ègal à ( det(A)).
3
4) Si A est une matrice diagonale, triangulaire supérieure ou inférieure
Q
n
alors det(A) = aii :
i=1
det In = 1:
Exemple
2 3 1 5 2 5 3 5 1 5 2 3 1 2 3 1
5 5 2 4 = 5 2 4 = 5 5 2 4 = 5 5 5 ( 2) 5 4
3 2 3 3 2 3 5 3 2 3 3 2 3
6 2 6 3 6 1 6 2 6 3 6 1 2 3 1
3 5 3 ( 2) 3 4 = 3 5 ( 2) 4 = 3 6 5 ( 2) 4
3 2 3 3 2 3 3 2 3
Propriétés Soient A; B 2 Mn (R)
1) A est inversible (i.e A 1 existe) si et seulement si det(A) 6= 0:
2) det(t A) = det(A)
3) det(AB) = det(A) det(B)
1
4) Si A 1 existe alors det(A 1 ) =
det(A)
t
1
5) Si A 1 existe alors A 1 = (( 1)i+j det(Ai; j ))1 i; j n où Ai; j
det(A)
est la matrice d’ordre n 1, obtenue à partir de A en supprimant la i eme
ligne et la j eme colonne de A.
1
a b 1 d b
6) Si ad cb 6= 0 alors =
c d ad cb c a
1 4 5 5
3 5 1 4 5 2 2
2 2
= = 2 3 = 3
2 4 2 2 3 2 2
1 2
0 1 1 0 1 1
a11 :: 0 a11 :: 0
7) Si a11 a22 :::ann 6= 0 alors @ :: :: :: A = @ :: :: :: A
0 :: ann 0 :: ann1
0 1 1 0 1
1 0 0 0 1 0 0 0
B 0 2 0 0 C B 0 1 0 0 C
B C B 2 C
@ 0 0 3 0 A =@ 0 0 1 0 A
3
0 0 0 4 0 0 0 41
(l’inverse d’une matrice diagonale inversible est une matrice diagonale)
III) Décomposition par blocs
Soit A une matrice du type (n; m). On peut exprimer la matrice A comme
un tableau de matrices comme suit :
4
0 1
A11 ::: ::: A1M
B ::: ::: C
B C
B
A=B AIJ C avec 1 I N et 1 J M
C
@ ::: ::: A
AN 1 ::: ::: AN M
P
N P
M
où AIJ est une matrice du type (nI ; mJ ) et on a n = nI et m = mJ .
I=1 J=1
Exemple
0 1
2 1 3 6 8 7 1
B 6 8 5 2 5 1 0 C
B C
A=B B 1 4 5 0 9 8 7 CC du type (5,7)
@ 1 2 8 9 0 2 3 A
1 2 2 8 4 5 5
1) 1ère décomposition : on décompose A en 4 blocs, comme suit
lignes : 5=2+3
colonnes : 7=3+4
alors
A11 A12
A=
A21 A22
0 1
2 1 3 6 8 7 1
B 6 8 5 2 5 1 0 C
2 1 3 6 8 7 1 B C
avec A11 = ; A12 = A=B
B 1 4 5 0 9 8 7 C
C
6 8 5 2 5 1 0 @ A
1 2 8 9 0 2 3
0 1 0 1 1 2 2 8 4 5 5
1 4 5 0 9 8 7
A21 = @ 1 2 8 A ; A22 = @ 9 0 2 3 A
1 2 2 8 4 5 5
2) 2ème décomposition : on décompose A en 6 blocs, comme suit
lignes : 5=3+2
colonnes : 7=2+3+2
A11 A12 A13
alors A =
0A21 A1 22 A23
0 1 0 1
2 1 3 6 8 7 1
avec A11 = @ 6 8 A ; A12 = @ 5 2 5 A ; A13 = @ 1 0 A
1 4 5 0 9 8 7
1 2 8 9 0 2 3
A21 = ; A22 = ; A23 =
1 2 2 8 4 5 5
5
Remarque : on peut décomposer une matrice en blocs de matrices de
plusiseurs façons.
L’interêt de la décomposition par blocs de matrices, réside dans le fait
que certaines opératrions valables sur les éléments, restent valables sur les
blocs : par exemple la somme et le produit de matrices.
– Le produit par blocs : soient A = (aij ) du type (n; p), B = (bij ) du
type (p; m) et C = AB.
P
N PP
On décompose A et B en blocs comme suit : n = nI , p = pk et
I=1 K=1
P
M
m= mJ
J=1
P
p
Les éléments de C sont donnés par : cij = aik bkj avec 1 i n et
k=1
1 j m.
P
P
Les blocs de C sont donnés par : CIJ = AIK BKJ avec 1 I N et
K=1
1 J M.
Théorème 1
Soit A une matrice carrée diagonale (resp. triangulaire supérieure),(resp.
trinagulaire inférieure) par blocs telle que les blocs diagonaux Aii ; i = 1; :::; r
Qr
de A sont carrés alors det(A) = det(Aii ).
i=1
Théorème 2
Soit A une matrice carrée0inversible et diagonale1par blocs carrés, notée
A111 0 ::: 0
B 0 A221 ::: 0 C
Aii ; i = 1; :::; r, alors A 1 = B
@ :::
C:
::: ::: 0 A
0 0 ::: Arr1
Exemple
11 ) = ja11 j = a11
det (a0 1
1 3 0 0 0 0 0
B 2 4 0 0 0 0 0 C
B C
B 0 0 1 2 0 0 0 C
B C
A=B B 0 0 1 4 0 0 0 C
C matrice 7 7 diagonale par blocs
B 0 0 0 0 4 2 0 C
B C
@ 0 0 0 0 1 3 0 A
0 0 0 0 0 0 5
carrés
6
1 3 1 2 4 2
det A = ( 5) = ( 2)(2)(10)( 5) = 200
2 4 1 4 1 3
0 1
1 3 0 0 0 0 0
B 2 4 0 0 0 0 0 C
B C
B 0 0 1 2 0 0 0 C
B C
si A = B
B 0 0 1 4 0 0 0 C C
B 0 0 0 0 4 2 0 C
B C
@ 0 0 0 0 1 3 0 A
0 0 0 0 0 0 5
0 1
4= 2 3= 2 0 0 0 0 0
B 2= 2 1= 2 0 0 0 0 0 C
B C
B 0 0 4=2 2=2 0 0 0 C
B C
alors A 1
=B
B 0 0 1=2 1=2 0 0 0 C
C
B 0 0 0 0 3=10 2=10 0 C
B C
@ 0 0 0 0 1=10 4=10 0 A
1
0 0 0 0 0 0 5