Les matrices sont l’outil de calcul privilégié de l’Algèbre Linéaire . Elles permettent de représenter les applications linéaires par des tableaux de nombres et de ramener les problèmes abstraits à des calculs concrets.
Matrices et opérations, matrice d’une application linéaire, changement de base, inversibilité, rang, systèmes linéaires (pivot de Gauss), déterminants, cofacteurs, formule de Cramer.
1. Matrices : définitions et opérations§
1.1 Définition§
Une matrice A ∈ M n , p ( K ) A \in \mathcal{M}_{n,p}(K) A ∈ M n , p ( K ) est un tableau de n n n lignes et p p p colonnes à coefficients dans K K K :
A = ( a i , j ) 1 ≤ i ≤ n 1 ≤ j ≤ p = ( a 1 , 1 a 1 , 2 ⋯ a 1 , p a 2 , 1 a 2 , 2 ⋯ a 2 , p ⋮ ⋮ ⋱ ⋮ a n , 1 a n , 2 ⋯ a n , p ) A = (a_{i,j})_{\substack{1 \leq i \leq n \ 1 \leq j \leq p}} = \begin{pmatrix} a_{1,1} & a_{1,2} & \cdots & a_{1,p} \ a_{2,1} & a_{2,2} & \cdots & a_{2,p} \ \vdots & \vdots & \ddots & \vdots \ a_{n,1} & a_{n,2} & \cdots & a_{n,p} \end{pmatrix} A = ( a i , j ) 1 ≤ i ≤ n 1 ≤ j ≤ p = ( a 1 , 1 a 1 , 2 ⋯ a 1 , p a 2 , 1 a 2 , 2 ⋯ a 2 , p ⋮ ⋮ ⋱ ⋮ a n , 1 a n , 2 ⋯ a n , p )
M n , p ( K ) \mathcal{M}_{n,p}(K) M n , p ( K ) est un K K K -espace vectoriel de dimension n p np n p .
Notation : M n ( K ) = M n , n ( K ) \mathcal{M}_n(K) = \mathcal{M}_{n,n}(K) M n ( K ) = M n , n ( K ) (matrices carrées).
1.2 Opérations§
Somme : ( A + B ) i , j = a i , j + b i , j (A + B)_{i,j} = a_{i,j} + b_{i,j} ( A + B ) i , j = a i , j + b i , j
Produit par un scalaire : ( λ A ) i , j = λ ⋅ a i , j (\lambda A)_{i,j} = \lambda \cdot a_{i,j} ( λ A ) i , j = λ ⋅ a i , j
Produit matriciel : Pour A ∈ M n , p ( K ) A \in \mathcal{M}_{n,p}(K) A ∈ M n , p ( K ) et B ∈ M p , q ( K ) B \in \mathcal{M}_{p,q}(K) B ∈ M p , q ( K ) :
( A B ) i , j = ∑ k = 1 p a i , k ⋅ b k , j (AB)_{i,j} = \sum_{k=1}^{p} a_{i,k} \cdot b_{k,j} ( A B ) i , j = k = 1 ∑ p a i , k ⋅ b k , j
Le produit A B AB A B est défini si et seulement si le nombre de colonnes de A A A égale le nombre de lignes de B B B .
Attention — Le produit matriciel n’est PAS commutatif En général, A B ≠ B A AB \neq BA A B = B A . Même si A A A et B B B sont carrées, A B AB A B et B A BA B A peuvent être différentes, et l’une peut exister sans l’autre (si les tailles diffèrent).
Transposée : ( t ! A ) i , j = a j , i ({}^t!A)_{i,j} = a_{j,i} ( t ! A ) i , j = a j , i . Propriétés :
t ! ( A + B ) = t ! A + t ! B {}^t!(A + B) = {}^t!A + {}^t!B t ! ( A + B ) = t ! A + t ! B
t ! ( λ A ) = λ , t ! A {}^t!(\lambda A) = \lambda ,{}^t!A t ! ( λ A ) = λ , t ! A
t ! ( A B ) = t ! B ⋅ t ! A {}^t!(AB) = {}^t!B \cdot {}^t!A t ! ( A B ) = t ! B ⋅ t ! A (inversion de l’ordre !)
t ! ( t ! A ) = A {}^t!({}^t!A) = A t ! ( t ! A ) = A
1.3 Matrices particulières§
Type Définition Matrice identité I n I_n I n ( δ i , j ) (\delta_{i,j}) ( δ i , j ) (symbole de Kronecker)Matrice diagonale a i , j = 0 a_{i,j} = 0 a i , j = 0 si i ≠ j i \neq j i = j Matrice triangulaire supérieure a i , j = 0 a_{i,j} = 0 a i , j = 0 si i > j i > j i > j Matrice triangulaire inférieure a i , j = 0 a_{i,j} = 0 a i , j = 0 si i < j i < j i < j Matrice symétrique t ! A = A {}^t!A = A t ! A = A Matrice antisymétrique t ! A = − A {}^t!A = -A t ! A = − A Matrice nilpotente ∃ , p ∈ N ∗ , ; A p = 0 \exists, p \in \mathbb{N}^*,; A^p = 0 ∃ , p ∈ N ∗ ,; A p = 0
Décomposition symétrique / antisymétrique Toute matrice carrée A A A s’écrit de manière unique :
A = A + t ! A 2 ⏟ sym e ˊ trique + A − t ! A 2 ⏟ antisym e ˊ trique A = \underbrace{\frac{A + {}^t!A}{2}}_{\text{symétrique}} + \underbrace{\frac{A - {}^t!A}{2}}_{\text{antisymétrique}} A = sym e ˊ trique 2 A + t ! A + antisym e ˊ trique 2 A − t ! A
2. Matrice d’une application linéaire§
2.1 Construction§
Soient E E E de base B = ( e 1 , … , e p ) \mathcal{B} = (e_1, \ldots, e_p) B = ( e 1 , … , e p ) et F F F de base C = ( f 1 , … , f n ) \mathcal{C} = (f_1, \ldots, f_n) C = ( f 1 , … , f n ) . La matrice de f ∈ L ( E , F ) f \in \mathcal{L}(E, F) f ∈ L ( E , F ) dans les bases B \mathcal{B} B et C \mathcal{C} C est la matrice M = Mat B , C ( f ) ∈ M n , p ( K ) M = \text{Mat}_{\mathcal{B}, \mathcal{C}}(f) \in \mathcal{M}_{n,p}(K) M = Mat B , C ( f ) ∈ M n , p ( K ) dont la j j j -ème colonne contient les coordonnées de f ( e j ) f(e_j) f ( e j ) dans la base C \mathcal{C} C :
f ( e j ) = ∑ i = 1 n m i , j f i f(e_j) = \sum_{i=1}^{n} m_{i,j} f_i f ( e j ) = i = 1 ∑ n m i , j f i
2.2 Lien avec les coordonnées§
Si X = Mat B ( u ) X = \text{Mat}_\mathcal{B}(u) X = Mat B ( u ) est le vecteur colonne des coordonnées de u u u dans B \mathcal{B} B et Y = Mat C ( f ( u ) ) Y = \text{Mat}_\mathcal{C}(f(u)) Y = Mat C ( f ( u )) celles de f ( u ) f(u) f ( u ) dans C \mathcal{C} C , alors :
Y = M X Y = MX Y = M X
2.3 Composition et produit§
Mat ( g ∘ f ) = Mat ( g ) × Mat ( f ) \text{Mat}(g \circ f) = \text{Mat}(g) \times \text{Mat}(f) Mat ( g ∘ f ) = Mat ( g ) × Mat ( f )
L’isomorphisme L ( E , F ) → ∼ M n , p ( K ) \mathcal{L}(E, F) \xrightarrow{\sim} \mathcal{M}_{n,p}(K) L ( E , F ) ∼ M n , p ( K ) est un isomorphisme d’espaces vectoriels (et d’algèbres quand E = F E = F E = F ).
3. Changement de base§
3.1 Matrice de passage§
Définition — Matrice de passage Soient B = ( e 1 , … , e n ) \mathcal{B} = (e_1, \ldots, e_n) B = ( e 1 , … , e n ) et B ′ = ( e 1 ′ , … , e n ′ ) \mathcal{B}' = (e_1', \ldots, e_n') B ′ = ( e 1 ′ , … , e n ′ ) deux bases de E E E . La matrice de passage de B \mathcal{B} B à B ′ \mathcal{B}' B ′ est :
P = P B → B ′ = Mat B ( Id E ) lu de B ′ vers B P = P_{\mathcal{B} \to \mathcal{B}'} = \text{Mat}_\mathcal{B}(\text{Id}_E) \text{ lu de } \mathcal{B}' \text{ vers } \mathcal{B} P = P B → B ′ = Mat B ( Id E ) lu de B ′ vers B
La j j j -ème colonne de P P P contient les coordonnées de e j ′ e_j' e j ′ dans la base B \mathcal{B} B .
P P P est toujours inversible, et P − 1 = P B ′ → B P^{-1} = P_{\mathcal{B}' \to \mathcal{B}} P − 1 = P B ′ → B .
Coordonnées d’un vecteur : Si X X X et X ′ X' X ′ sont les coordonnées de u u u dans B \mathcal{B} B et B ′ \mathcal{B}' B ′ :
X = P X ′ X = PX' X = P X ′
Matrice d’un endomorphisme : Si M = Mat B ( f ) M = \text{Mat}_\mathcal{B}(f) M = Mat B ( f ) et M ′ = Mat B ′ ( f ) M' = \text{Mat}_{\mathcal{B}'}(f) M ′ = Mat B ′ ( f ) :
M ′ = P − 1 M P \boxed{M' = P^{-1} M P} M ′ = P − 1 M P
Deux matrices liées par cette relation sont dites semblables .
Matrice d’une application linéaire : Si f : E → F f : E \to F f : E → F avec changement B → B ′ \mathcal{B} \to \mathcal{B}' B → B ′ dans E E E (passage P P P ) et C → C ′ \mathcal{C} \to \mathcal{C}' C → C ′ dans F F F (passage Q Q Q ) :
M ′ = Q − 1 M P M' = Q^{-1} M P M ′ = Q − 1 M P
4. Matrices inversibles§
4.1 Définition§
A ∈ M n ( K ) A \in \mathcal{M}_n(K) A ∈ M n ( K ) est inversible s’il existe B ∈ M n ( K ) B \in \mathcal{M}_n(K) B ∈ M n ( K ) telle que A B = B A = I n AB = BA = I_n A B = B A = I n . On note B = A − 1 B = A^{-1} B = A − 1 .
Le groupe des matrices inversibles est noté G L n ( K ) GL_n(K) G L n ( K ) .
4.2 Critères d’inversibilité§
A A A est inversible ⟺ \iff ⟺ l’une des conditions équivalentes suivantes est vérifiée :
det ( A ) ≠ 0 \det(A) \neq 0 det ( A ) = 0
Les colonnes de A A A forment une famille libre (ou une base) de K n K^n K n
Le système A X = 0 AX = 0 A X = 0 admet uniquement la solution triviale X = 0 X = 0 X = 0
L’endomorphisme associé est bijectif
rg ( A ) = n \text{rg}(A) = n rg ( A ) = n
4.3 Calcul de l’inverse par la méthode du pivot§
On forme le tableau augmenté ( A ∣ I n ) (A \mid I_n) ( A ∣ I n ) et on effectue des opérations élémentaires sur les lignes pour obtenir ( I n ∣ A − 1 ) (I_n \mid A^{-1}) ( I n ∣ A − 1 ) .
A = ( 1 2 3 7 ) A = \begin{pmatrix} 1 & 2 \ 3 & 7 \end{pmatrix} A = ( 1 2 3 7 )
( A ∣ I 2 ) = ( 1 2 1 0 3 7 0 1 ) → L 2 ← L 2 − 3 L 1 ( 1 2 1 0 0 1 − 3 1 ) → L 1 ← L 1 − 2 L 2 ( 1 0 7 − 2 0 1 − 3 1 ) (A \mid I_2) = \left(\begin{array}{cc|cc} 1 & 2 & 1 & 0 \ 3 & 7 & 0 & 1 \end{array}\right) \xrightarrow{L_2 \leftarrow L_2 - 3L_1} \left(\begin{array}{cc|cc} 1 & 2 & 1 & 0 \ 0 & 1 & -3 & 1 \end{array}\right) \xrightarrow{L_1 \leftarrow L_1 - 2L_2} \left(\begin{array}{cc|cc} 1 & 0 & 7 & -2 \ 0 & 1 & -3 & 1 \end{array}\right) ( A ∣ I 2 ) = ( 1 2 1 0 3 7 0 1 ) L 2 ← L 2 − 3 L 1 ( 1 2 1 0 0 1 − 3 1 ) L 1 ← L 1 − 2 L 2 ( 1 0 7 − 2 0 1 − 3 1 )
Donc A − 1 = ( 7 − 2 − 3 1 ) A^{-1} = \begin{pmatrix} 7 & -2 \ -3 & 1 \end{pmatrix} A − 1 = ( 7 − 2 − 3 1 ) . Vérification : det ( A ) = 7 − 6 = 1 \det(A) = 7 - 6 = 1 det ( A ) = 7 − 6 = 1 , et A − 1 = 1 1 ( 7 − 2 − 3 1 ) A^{-1} = \frac{1}{1}\begin{pmatrix} 7 & -2 \ -3 & 1 \end{pmatrix} A − 1 = 1 1 ( 7 − 2 − 3 1 ) .
5. Rang d’une matrice§
Le rang d’une matrice A ∈ M n , p ( K ) A \in \mathcal{M}_{n,p}(K) A ∈ M n , p ( K ) est la dimension du sous-espace vectoriel engendré par ses colonnes (ou de manière équivalente, par ses lignes) :
rg ( A ) = dim ( Im ( f A ) ) \text{rg}(A) = \dim(\text{Im}(f_A)) rg ( A ) = dim ( Im ( f A ))
où f A f_A f A est l’application linéaire canoniquement associée.
Propriétés :
rg ( A ) = rg ( t ! A ) \text{rg}(A) = \text{rg}({}^t!A) rg ( A ) = rg ( t ! A )
rg ( A B ) ≤ min ( rg ( A ) , rg ( B ) ) \text{rg}(AB) \leq \min(\text{rg}(A), \text{rg}(B)) rg ( A B ) ≤ min ( rg ( A ) , rg ( B ))
Si P P P ou Q Q Q est inversible : rg ( P A Q ) = rg ( A ) \text{rg}(PAQ) = \text{rg}(A) rg ( P A Q ) = rg ( A )
Le rang est le nombre de pivots non nuls après échelonnement
6. Systèmes linéaires§
6.1 Écriture matricielle§
Le système de n n n équations à p p p inconnues :
{ a 1 , 1 x 1 + ⋯ + a 1 , p x p = b 1 ⋮ a n , 1 x 1 + ⋯ + a n , p x p = b n \begin{cases} a_{1,1}x_1 + \cdots + a_{1,p}x_p = b_1 \ \vdots \ a_{n,1}x_1 + \cdots + a_{n,p}x_p = b_n \end{cases} { a 1 , 1 x 1 + ⋯ + a 1 , p x p = b 1 ⋮ a n , 1 x 1 + ⋯ + a n , p x p = b n
s’écrit A X = B AX = B A X = B avec A ∈ M n , p ( K ) A \in \mathcal{M}_{n,p}(K) A ∈ M n , p ( K ) , X ∈ M p , 1 ( K ) X \in \mathcal{M}_{p,1}(K) X ∈ M p , 1 ( K ) , B ∈ M n , 1 ( K ) B \in \mathcal{M}_{n,1}(K) B ∈ M n , 1 ( K ) .
6.2 Structure de l’ensemble des solutions§
Théorème — Structure affine L’ensemble des solutions de A X = B AX = B A X = B est :
Vide si le système est incompatible
Un sous-espace affine x 0 + ker ( A ) x_0 + \ker(A) x 0 + ker ( A ) sinon, où x 0 x_0 x 0 est une solution particulière
En particulier, les solutions du système homogène A X = 0 AX = 0 A X = 0 forment un sev de dimension p − rg ( A ) p - \text{rg}(A) p − rg ( A ) .
6.3 Méthode du pivot de Gauss§
flowchart TD
A["Écrire le système sous<br>forme matricielle augmentée<br>(A | B)"] --> B["Choisir un pivot non nul<br>dans la colonne courante"]
B --> C{"Pivot trouvé ?"}
C -- Non --> D["Passer à la<br>colonne suivante"]
C -- Oui --> E["Échanger les lignes<br>si nécessaire"]
E --> F["Éliminer les coefficients<br>sous le pivot par<br>L_i ← L_i - (a_ij/pivot) × L_pivot"]
F --> G{"Dernière colonne<br>traitée ?"}
G -- Non --> B
G -- Oui --> H["Système échelonné obtenu"]
H --> I{"Vérifier la<br>compatibilité"}
I -- "Ligne 0 = c ≠ 0" --> J["Système<br>incompatible<br>∅"]
I -- Compatible --> K["Remontée :<br>résolution par<br>substitution arrière"]
K --> L["Paramétrer les<br>variables libres"]
L --> M["Écrire l'ensemble<br>des solutions"]
style A fill:#3498db,stroke:#333,color:#fff
style J fill:#e74c3c,stroke:#333,color:#fff
style M fill:#2ecc71,stroke:#333,color:#fff
{ x + 2 y + z = 1 2 x + 5 y + 3 z = 3 x + 3 y + 2 z = 2 \begin{cases} x + 2y + z = 1 \ 2x + 5y + 3z = 3 \ x + 3y + 2z = 2 \end{cases} { x + 2 y + z = 1 2 x + 5 y + 3 z = 3 x + 3 y + 2 z = 2
Tableau augmenté :
( 1 2 1 1 2 5 3 3 1 3 2 2 ) → L 3 ← L 3 − L 1 L 2 ← L 2 − 2 L 1 ( 1 2 1 1 0 1 1 1 0 1 1 1 ) → L 3 ← L 3 − L 2 ( 1 2 1 1 0 1 1 1 0 0 0 0 ) \left(\begin{array}{ccc|c} 1 & 2 & 1 & 1 \ 2 & 5 & 3 & 3 \ 1 & 3 & 2 & 2 \end{array}\right) \xrightarrow[\substack{L_3 \leftarrow L_3 - L_1}]{L_2 \leftarrow L_2 - 2L_1} \left(\begin{array}{ccc|c} 1 & 2 & 1 & 1 \ 0 & 1 & 1 & 1 \ 0 & 1 & 1 & 1 \end{array}\right) \xrightarrow{L_3 \leftarrow L_3 - L_2} \left(\begin{array}{ccc|c} 1 & 2 & 1 & 1 \ 0 & 1 & 1 & 1 \ 0 & 0 & 0 & 0 \end{array}\right) ( 1 2 1 1 2 5 3 3 1 3 2 2 ) L 2 ← L 2 − 2 L 1 L 3 ← L 3 − L 1 ( 1 2 1 1 0 1 1 1 0 1 1 1 ) L 3 ← L 3 − L 2 ( 1 2 1 1 0 1 1 1 0 0 0 0 )
Le rang est 2 2 2 , il y a 3 − 2 = 1 3 - 2 = 1 3 − 2 = 1 variable libre. On pose z = t ∈ R z = t \in \mathbb{R} z = t ∈ R :
y = 1 − t y = 1 - t y = 1 − t
x = 1 − 2 ( 1 − t ) − t = − 1 + t x = 1 - 2(1-t) - t = -1 + t x = 1 − 2 ( 1 − t ) − t = − 1 + t
Solutions : ( − 1 + t , ; 1 − t , ; t ) ∣ t ∈ R = ( − 1 , 1 , 0 ) + Vect ( 1 , − 1 , 1 ) {(-1 + t,; 1 - t,; t) \mid t \in \mathbb{R}} = (-1, 1, 0) + \text{Vect}(1, -1, 1) ( − 1 + t ,; 1 − t ,; t ) ∣ t ∈ R = ( − 1 , 1 , 0 ) + Vect ( 1 , − 1 , 1 ) .
7. Déterminants§
7.1 Définition par les permutations§
Le déterminant d’une matrice carrée A = ( a i , j ) ∈ M n ( K ) A = (a_{i,j}) \in \mathcal{M}_n(K) A = ( a i , j ) ∈ M n ( K ) est :
det ( A ) = ∑ σ ∈ S n ε ( σ ) ∏ i = 1 n a i , σ ( i ) \det(A) = \sum_{\sigma \in \mathfrak{S}_n} \varepsilon(\sigma) \prod_{i=1}^{n} a_{i,\sigma(i)} det ( A ) = σ ∈ S n ∑ ε ( σ ) i = 1 ∏ n a i , σ ( i )
où S n \mathfrak{S}_n S n est le groupe des permutations de 1 , … , n {1, \ldots, n} 1 , … , n et ε ( σ ) \varepsilon(\sigma) ε ( σ ) la signature de σ \sigma σ .
Pour n = 2 n = 2 n = 2 :
det ( a b c d ) = a d − b c \det\begin{pmatrix} a & b \ c & d \end{pmatrix} = ad - bc det ( a b c d ) = a d − b c
Pour n = 3 n = 3 n = 3 (règle de Sarrus) :
det ( a b c d e f g h i ) = a e i + b f g + c d h − c e g − b d i − a f h \det\begin{pmatrix} a & b & c \ d & e & f \ g & h & i \end{pmatrix} = aei + bfg + cdh - ceg - bdi - afh det ( a b c d e f g h i ) = a e i + b f g + c d h − ce g − b d i − a f h
7.2 Propriétés fondamentales§
Le déterminant est caractérisé comme l’unique application det : M n ( K ) → K \det : \mathcal{M}_n(K) \to K det : M n ( K ) → K vérifiant :
n n n -linéarité : linéaire par rapport à chaque colonne (ou ligne)
Alternée : s’annule dès que deux colonnes sont égales
Normalisation : det ( I n ) = 1 \det(I_n) = 1 det ( I n ) = 1
det ( t ! A ) = det ( A ) \det({}^t!A) = \det(A) det ( t ! A ) = det ( A )
det ( A B ) = det ( A ) ⋅ det ( B ) \det(AB) = \det(A) \cdot \det(B) det ( A B ) = det ( A ) ⋅ det ( B ) (morphisme multiplicatif)
det ( λ A ) = λ n det ( A ) \det(\lambda A) = \lambda^n \det(A) det ( λ A ) = λ n det ( A )
det ( A − 1 ) = 1 det ( A ) \det(A^{-1}) = \frac{1}{\det(A)} det ( A − 1 ) = d e t ( A ) 1 (si A A A inversible)
Échange de deux lignes (ou colonnes) : change le signe
Ajout d’un multiple d’une ligne à une autre : ne change pas le déterminant
Multiplication d’une ligne par λ \lambda λ : multiplie le déterminant par λ \lambda λ
Théorème — Critère d’inversibilité
A inversible ⟺ det ( A ) ≠ 0 A \text{ inversible} \iff \det(A) \neq 0 A inversible ⟺ det ( A ) = 0
7.3 Déterminant de matrices triangulaires§
det ( a 1 , 1 ∗ ⋯ ∗ 0 a 2 , 2 ⋯ ∗ ⋮ ⋱ ⋱ ⋮ 0 ⋯ 0 a n , n ) = ∏ i = 1 n a i , i \det\begin{pmatrix} a_{1,1} & * & \cdots & * \ 0 & a_{2,2} & \cdots & * \ \vdots & \ddots & \ddots & \vdots \ 0 & \cdots & 0 & a_{n,n} \end{pmatrix} = \prod_{i=1}^n a_{i,i} det ( a 1 , 1 ∗ ⋯ ∗ 0 a 2 , 2 ⋯ ∗ ⋮ ⋱ ⋱ ⋮ 0 ⋯ 0 a n , n ) = i = 1 ∏ n a i , i
Le déterminant d’une matrice triangulaire (supérieure ou inférieure) est le produit de ses coefficients diagonaux .
8. Développement selon une ligne ou colonne§
8.1 Cofacteurs§
Le cofacteur C i , j C_{i,j} C i , j est :
C i , j = ( − 1 ) i + j M i , j C_{i,j} = (-1)^{i+j} M_{i,j} C i , j = ( − 1 ) i + j M i , j
où M i , j M_{i,j} M i , j est le mineur d’ordre ( i , j ) (i,j) ( i , j ) , c’est-à-dire le déterminant de la matrice obtenue en supprimant la i i i -ème ligne et la j j j -ème colonne.
Théorème — Développement par rapport à la i i i -ème ligne
det ( A ) = ∑ j = 1 n a i , j C i , j = ∑ j = 1 n ( − 1 ) i + j a i , j M i , j \det(A) = \sum_{j=1}^{n} a_{i,j} C_{i,j} = \sum_{j=1}^{n} (-1)^{i+j} a_{i,j} M_{i,j} det ( A ) = j = 1 ∑ n a i , j C i , j = j = 1 ∑ n ( − 1 ) i + j a i , j M i , j
De même, développement par rapport à la j j j -ème colonne :
det ( A ) = ∑ i = 1 n a i , j C i , j \det(A) = \sum_{i=1}^{n} a_{i,j} C_{i,j} det ( A ) = i = 1 ∑ n a i , j C i , j
Développer par rapport à la ligne ou colonne contenant le plus de zéros pour minimiser les calculs.
det ( 2 0 1 3 1 0 1 0 4 ) \det\begin{pmatrix} 2 & 0 & 1 \ 3 & 1 & 0 \ 1 & 0 & 4 \end{pmatrix} det ( 2 0 1 3 1 0 1 0 4 )
Développement selon la 2e colonne (deux zéros) :
= 0 ⋅ C 1 , 2 + 1 ⋅ C 2 , 2 + 0 ⋅ C 3 , 2 = ( − 1 ) 2 + 2 det ( 2 1 1 4 ) = 8 − 1 = 7 = 0 \cdot C_{1,2} + 1 \cdot C_{2,2} + 0 \cdot C_{3,2} = (-1)^{2+2} \det\begin{pmatrix} 2 & 1 \ 1 & 4 \end{pmatrix} = 8 - 1 = 7 = 0 ⋅ C 1 , 2 + 1 ⋅ C 2 , 2 + 0 ⋅ C 3 , 2 = ( − 1 ) 2 + 2 det ( 2 1 1 4 ) = 8 − 1 = 7
9.1 Comatrice§
La comatrice (ou matrice des cofacteurs) de A A A est Com ( A ) = ( C i , j ) 1 ≤ i , j ≤ n \text{Com}(A) = (C_{i,j})_{1 \leq i,j \leq n} Com ( A ) = ( C i , j ) 1 ≤ i , j ≤ n .
A ⋅ t ! Com ( A ) = t ! Com ( A ) ⋅ A = det ( A ) ⋅ I n A \cdot {}^t!\text{Com}(A) = {}^t!\text{Com}(A) \cdot A = \det(A) \cdot I_n A ⋅ t ! Com ( A ) = t ! Com ( A ) ⋅ A = det ( A ) ⋅ I n
En particulier, si det ( A ) ≠ 0 \det(A) \neq 0 det ( A ) = 0 :
A − 1 = 1 det ( A ) , t ! Com ( A ) A^{-1} = \frac{1}{\det(A)} ,{}^t!\text{Com}(A) A − 1 = det ( A ) 1 , t ! Com ( A )
Théorème — Formules de Cramer Pour le système A X = B AX = B A X = B avec A ∈ G L n ( K ) A \in GL_n(K) A ∈ G L n ( K ) :
x j = det ( A j ) det ( A ) x_j = \frac{\det(A_j)}{\det(A)} x j = det ( A ) det ( A j )
où A j A_j A j est la matrice obtenue en remplaçant la j j j -ème colonne de A A A par B B B .
Les formules de Cramer sont théoriquement élégantes mais coûteuses en calcul pour n n n grand. En pratique, on préfère le pivot de Gauss.
Exemple — Cramer en dimension 2
{ 3 x + 2 y = 7 x + 4 y = 9 \begin{cases} 3x + 2y = 7 \ x + 4y = 9 \end{cases} { 3 x + 2 y = 7 x + 4 y = 9
det ( A ) = det ( 3 2 1 4 ) = 10 \det(A) = \det\begin{pmatrix} 3 & 2 \ 1 & 4 \end{pmatrix} = 10 det ( A ) = det ( 3 2 1 4 ) = 10
x = det ( 7 2 9 4 ) 10 = 28 − 18 10 = 1 , y = det ( 3 7 1 9 ) 10 = 27 − 7 10 = 2 x = \frac{\det\begin{pmatrix} 7 & 2 \ 9 & 4 \end{pmatrix}}{10} = \frac{28-18}{10} = 1, \quad y = \frac{\det\begin{pmatrix} 3 & 7 \ 1 & 9 \end{pmatrix}}{10} = \frac{27-7}{10} = 2 x = 10 det ( 7 2 9 4 ) = 10 28 − 18 = 1 , y = 10 det ( 3 7 1 9 ) = 10 27 − 7 = 2
10. Déterminants classiques§
10.1 Déterminant de Vandermonde§
Théorème — Déterminant de Vandermonde
V n ( x 1 , … , x n ) = det ( 1 1 ⋯ 1 x 1 x 2 ⋯ x n x 1 2 x 2 2 ⋯ x n 2 ⋮ ⋮ ⋱ ⋮ x 1 n − 1 x 2 n − 1 ⋯ x n n − 1 ) = ∏ 1 ≤ i < j ≤ n ( x j − x i ) V_n(x_1, \ldots, x_n) = \det\begin{pmatrix} 1 & 1 & \cdots & 1 \ x_1 & x_2 & \cdots & x_n \ x_1^2 & x_2^2 & \cdots & x_n^2 \ \vdots & \vdots & \ddots & \vdots \ x_1^{n-1} & x_2^{n-1} & \cdots & x_n^{n-1} \end{pmatrix} = \prod_{1 \leq i < j \leq n} (x_j - x_i) V n ( x 1 , … , x n ) = det ( 1 1 ⋯ 1 x 1 x 2 ⋯ x n x 1 2 x 2 2 ⋯ x n 2 ⋮ ⋮ ⋱ ⋮ x 1 n − 1 x 2 n − 1 ⋯ x n n − 1 ) = 1 ≤ i < j ≤ n ∏ ( x j − x i )
On considère V n V_n V n comme un polynôme en x n x_n x n . En effectuant C n ← C n − C i C_n \leftarrow C_n - C_i C n ← C n − C i pour chaque i < n i < n i < n , on montre que x n − x i x_n - x_i x n − x i divise V n V_n V n pour tout i < n i < n i < n . Par récurrence et argument de degré, on obtient la formule.
Conséquence : V n ≠ 0 ⟺ V_n \neq 0 \iff V n = 0 ⟺ les x i x_i x i sont deux à deux distincts.
10.2 Déterminant par blocs§
Si A A A est triangulaire par blocs :
det ( A 1 ∗ 0 A 2 ) = det ( A 1 ) ⋅ det ( A 2 ) \det\begin{pmatrix} A_1 & * \ 0 & A_2 \end{pmatrix} = \det(A_1) \cdot \det(A_2) det ( A 1 ∗ 0 A 2 ) = det ( A 1 ) ⋅ det ( A 2 )
11. Exercices types corrigés§
Exercice 1 : Calcul de déterminant par opérations§
Calculer D = det ( 1 1 1 1 1 2 3 4 1 3 6 10 1 4 10 20 ) D = \det\begin{pmatrix} 1 & 1 & 1 & 1 \ 1 & 2 & 3 & 4 \ 1 & 3 & 6 & 10 \ 1 & 4 & 10 & 20 \end{pmatrix} D = det ( 1 1 1 1 1 2 3 4 1 3 6 10 1 4 10 20 ) .
Solution : Opérations C j ← C j − C j − 1 C_j \leftarrow C_j - C_{j-1} C j ← C j − C j − 1 pour j = 4 , 3 , 2 j = 4, 3, 2 j = 4 , 3 , 2 :
D = det ( 1 0 0 0 1 1 1 1 1 2 3 4 1 3 6 10 ) D = \det\begin{pmatrix} 1 & 0 & 0 & 0 \ 1 & 1 & 1 & 1 \ 1 & 2 & 3 & 4 \ 1 & 3 & 6 & 10 \end{pmatrix} D = det ( 1 0 0 0 1 1 1 1 1 2 3 4 1 3 6 10 )
On recommence : C j ← C j − C j − 1 C_j \leftarrow C_j - C_{j-1} C j ← C j − C j − 1 :
D = det ( 1 0 0 0 1 1 0 0 1 2 1 1 1 3 3 4 ) D = \det\begin{pmatrix} 1 & 0 & 0 & 0 \ 1 & 1 & 0 & 0 \ 1 & 2 & 1 & 1 \ 1 & 3 & 3 & 4 \end{pmatrix} D = det ( 1 0 0 0 1 1 0 0 1 2 1 1 1 3 3 4 )
Encore une fois :
D = det ( 1 0 0 0 1 1 0 0 1 2 1 0 1 3 3 1 ) = 1 × 1 × 1 × 1 = 1 D = \det\begin{pmatrix} 1 & 0 & 0 & 0 \ 1 & 1 & 0 & 0 \ 1 & 2 & 1 & 0 \ 1 & 3 & 3 & 1 \end{pmatrix} = 1 \times 1 \times 1 \times 1 = 1 D = det ( 1 0 0 0 1 1 0 0 1 2 1 0 1 3 3 1 ) = 1 × 1 × 1 × 1 = 1
(Matrice triangulaire inférieure.)
Exercice 2 : Matrice d’une application linéaire et changement de base§
Soit f : R 2 → R 2 f : \mathbb{R}^2 \to \mathbb{R}^2 f : R 2 → R 2 définie par f ( x , y ) = ( 2 x + y , x + 2 y ) f(x,y) = (2x+y, x+2y) f ( x , y ) = ( 2 x + y , x + 2 y ) . Écrire la matrice de f f f dans la base canonique, puis dans la base B ′ = ( ( 1 , 1 ) , ( 1 , − 1 ) ) \mathcal{B}' = ((1,1), (1,-1)) B ′ = (( 1 , 1 ) , ( 1 , − 1 )) .
Solution :
Base canonique B = ( e 1 , e 2 ) \mathcal{B} = (e_1, e_2) B = ( e 1 , e 2 ) :
f ( e 1 ) = f ( 1 , 0 ) = ( 2 , 1 ) f(e_1) = f(1,0) = (2,1) f ( e 1 ) = f ( 1 , 0 ) = ( 2 , 1 ) et f ( e 2 ) = f ( 0 , 1 ) = ( 1 , 2 ) f(e_2) = f(0,1) = (1,2) f ( e 2 ) = f ( 0 , 1 ) = ( 1 , 2 ) .
M = ( 2 1 1 2 ) M = \begin{pmatrix} 2 & 1 \ 1 & 2 \end{pmatrix} M = ( 2 1 1 2 )
Changement de base : La matrice de passage P P P de B \mathcal{B} B à B ′ \mathcal{B}' B ′ a pour colonnes les coordonnées de ( 1 , 1 ) (1,1) ( 1 , 1 ) et ( 1 , − 1 ) (1,-1) ( 1 , − 1 ) dans B \mathcal{B} B :
P = ( 1 1 1 − 1 ) , P − 1 = 1 − 2 ( − 1 − 1 − 1 1 ) = ( 1 / 2 1 / 2 1 / 2 − 1 / 2 ) P = \begin{pmatrix} 1 & 1 \ 1 & -1 \end{pmatrix}, \quad P^{-1} = \frac{1}{-2}\begin{pmatrix} -1 & -1 \ -1 & 1 \end{pmatrix} = \begin{pmatrix} 1/2 & 1/2 \ 1/2 & -1/2 \end{pmatrix} P = ( 1 1 1 − 1 ) , P − 1 = − 2 1 ( − 1 − 1 − 1 1 ) = ( 1/2 1/2 1/2 − 1/2 )
M ′ = P − 1 M P = ( 1 / 2 1 / 2 1 / 2 − 1 / 2 ) ( 2 1 1 2 ) ( 1 1 1 − 1 ) M' = P^{-1}MP = \begin{pmatrix} 1/2 & 1/2 \ 1/2 & -1/2 \end{pmatrix}\begin{pmatrix} 2 & 1 \ 1 & 2 \end{pmatrix}\begin{pmatrix} 1 & 1 \ 1 & -1 \end{pmatrix} M ′ = P − 1 M P = ( 1/2 1/2 1/2 − 1/2 ) ( 2 1 1 2 ) ( 1 1 1 − 1 )
P − 1 M = ( 3 / 2 3 / 2 1 / 2 − 1 / 2 ) , M ′ = ( 3 / 2 3 / 2 1 / 2 − 1 / 2 ) ( 1 1 1 − 1 ) = ( 3 0 0 1 ) P^{-1}M = \begin{pmatrix} 3/2 & 3/2 \ 1/2 & -1/2 \end{pmatrix}, \quad M' = \begin{pmatrix} 3/2 & 3/2 \ 1/2 & -1/2 \end{pmatrix}\begin{pmatrix} 1 & 1 \ 1 & -1 \end{pmatrix} = \begin{pmatrix} 3 & 0 \ 0 & 1 \end{pmatrix} P − 1 M = ( 3/2 3/2 1/2 − 1/2 ) , M ′ = ( 3/2 3/2 1/2 − 1/2 ) ( 1 1 1 − 1 ) = ( 3 0 0 1 )
Dans la base B ′ \mathcal{B}' B ′ , la matrice est diagonale ! Les vecteurs ( 1 , 1 ) (1,1) ( 1 , 1 ) et ( 1 , − 1 ) (1,-1) ( 1 , − 1 ) sont des vecteurs propres de f f f , de valeurs propres 3 3 3 et 1 1 1 .
Exercice 3 : Système linéaire paramétrique§
Discuter et résoudre selon m ∈ R m \in \mathbb{R} m ∈ R :
{ x + y + z = 1 x + m y + z = m x + y + m z = m 2 \begin{cases} x + y + z = 1 \ x + my + z = m \ x + y + mz = m^2 \end{cases} { x + y + z = 1 x + m y + z = m x + y + m z = m 2
Solution :
det ( A ) = det ( 1 1 1 1 m 1 1 1 m ) \det(A) = \det\begin{pmatrix} 1 & 1 & 1 \ 1 & m & 1 \ 1 & 1 & m \end{pmatrix} det ( A ) = det ( 1 1 1 1 m 1 1 1 m )
C 2 ← C 2 − C 1 C_2 \leftarrow C_2 - C_1 C 2 ← C 2 − C 1 , C 3 ← C 3 − C 1 C_3 \leftarrow C_3 - C_1 C 3 ← C 3 − C 1 :
= det ( 1 0 0 1 m − 1 0 1 0 m − 1 ) = ( m − 1 ) 2 = \det\begin{pmatrix} 1 & 0 & 0 \ 1 & m-1 & 0 \ 1 & 0 & m-1 \end{pmatrix} = (m-1)^2 = det ( 1 0 0 1 m − 1 0 1 0 m − 1 ) = ( m − 1 ) 2
Cas m ≠ 1 m \neq 1 m = 1 : det ( A ) = ( m − 1 ) 2 ≠ 0 \det(A) = (m-1)^2 \neq 0 det ( A ) = ( m − 1 ) 2 = 0 . Système de Cramer, solution unique. Par Cramer ou substitution :
L 2 − L 1 L_2 - L_1 L 2 − L 1 : ( m − 1 ) y = m − 1 (m-1)y = m - 1 ( m − 1 ) y = m − 1 , d’où y = 1 y = 1 y = 1
L 3 − L 1 L_3 - L_1 L 3 − L 1 : ( m − 1 ) z = m 2 − 1 = ( m − 1 ) ( m + 1 ) (m-1)z = m^2 - 1 = (m-1)(m+1) ( m − 1 ) z = m 2 − 1 = ( m − 1 ) ( m + 1 ) , d’où z = m + 1 z = m + 1 z = m + 1
L 1 L_1 L 1 : x = 1 − y − z = 1 − 1 − ( m + 1 ) = − ( m + 1 ) = − m − 1 x = 1 - y - z = 1 - 1 - (m+1) = -(m+1) = -m - 1 x = 1 − y − z = 1 − 1 − ( m + 1 ) = − ( m + 1 ) = − m − 1
Solution : ( x , y , z ) = ( − m − 1 , 1 , m + 1 ) (x,y,z) = (-m-1, 1, m+1) ( x , y , z ) = ( − m − 1 , 1 , m + 1 ) .
Cas m = 1 m = 1 m = 1 : Le système devient x + y + z = 1 x + y + z = 1 x + y + z = 1 (trois fois la même équation). Solutions : ( 1 − s − t , s , t ) ∣ s , t ∈ R {(1-s-t, s, t) \mid s, t \in \mathbb{R}} ( 1 − s − t , s , t ) ∣ s , t ∈ R , sous-espace affine de dimension 2 2 2 .
Liens§
Commentaires