Calcul des déterminants

Un article de Wikipédia, l'encyclopédie libre.

Il existe de nombreux procédés de calcul du déterminant d'une matrice carrée de taille n à coefficients réels, ou plus généralement à coefficients dans un corps K. La méthode la plus efficace en règle générale est l'utilisation d'une technique de pivot de Gauss.

[modifier] Cas particuliers simples

Les formules des premiers paragraphes s'appliquent dans le cas des dimensions égales à 2 ou à 3.
La règle de Sarrus est un procédé visuel, qui permet de retenir la formule de calcul des déterminants d’ordre 3. Ce n'est toutefois pas toujours la méthode la plus simple ou la plus rapide. Une approche basée sur les propriétés de linéarité du déterminant permet souvent d'effectuer moins d'opérations, ou d'obtenir une forme factorisée plus intéressante.

La règle de Sarrus consiste à écrire les trois colonnes de la matrice et à répéter dans l’ordre, les deux premières lignes en dessous de la matrice. Il suffit alors d’effectuer les produits des coefficients de chaque diagonale et d’en faire la somme si la diagonale est descendante ou la différence si la diagonale est ascendante. Plus clairement : pour calculer $\begin{vmatrix}a&b&c\\d&e&f\\g&h&i\end{vmatrix}$ , il suffit d'effectuer

a	b	c
d	e	f
g	h	i
a	b	c
d	e	f

a	b	c
d	e	f
g	h	i
a	b	c
d	e	f

affectés d'un signe positif

affectés d'un signe négatif

et le résultat est $(a\cdot e\cdot i + d\cdot h\cdot c + g\cdot b\cdot f) - (g\cdot e\cdot c +a\cdot h\cdot f +d\cdot b\cdot i)$ .

le déterminant d'une matrice triangulaire est le produit des coefficients diagonaux

$\begin{vmatrix} m_{1;1} & m_{1;2} & \cdots & m_{1;n-1} & m_{1;n} \\ 0 & m_{2;2} & \cdots & \cdots & m_{2;n} \\ \vdots & 0 & \ddots & & \vdots \\ \vdots & & \ddots & m_{n-1;n-1} & m_{n-1;n} \\ 0 & 0 & \cdots & 0 & m_{n;n} \end{vmatrix} = \prod_{i=1}^{i=n}{m_{i;i}}$

Démonstration :

On procède par récurrence. Il suffit d'appliquer la formule de Laplace à la première colonne pour se ramener d'une matrice de taille n à une matrice de taille n-1.

le déterminant d'une matrice triangulaire par blocs est le produit des déterminants des blocs diagonaux

$\det \begin{pmatrix} A & B \\ 0 & C \end{pmatrix}=\det A \det C$

Démonstration

On commence par simplifier la situation en utilisant le produit par blocs suivant

$\begin{pmatrix} A & B \\ 0 & C \end{pmatrix}=\begin{pmatrix} I & 0 \\ 0 & C \end{pmatrix}.\begin{pmatrix} A & B \\ 0 & I \end{pmatrix}$

Il suffit en effet de prouver que la première matrice a pour déterminant det C, la seconde det A. Mais pour cela on reprend la méthode de démonstration utilisée pour les matrices triangulaires. Ainsi pour la première matrice, on effectue des développements successifs par rapport aux premières lignes, qui sont les plus simples : il ne reste plus que le déterminant de C. Pour la deuxième matrice, on suit une méthode analogue avec les dernières lignes.

[modifier] Méthodes générales

1. Application directe de la formule de Leibniz

Appliquée en toute généralité, cette formule demande d'effectuer une somme de n! termes, chacun constitué d'un produit de n termes. Sans même prendre en compte le calcul de signature, le nombre de calculs nécessaires dépasse l'ordre (n+1)!.

2. Application récursive de la formule de Laplace

Le calcul d'un déterminant de taille n se ramène à une somme de n déterminants de taille n-1, chacun d'eux étant somme de n-1 déterminants de taille n-2, et ainsi de suite. Le temps de calcul croît de nouveau comme une factorielle en n.

Dans la pratique la formule de Laplace n'est intéressante que si plusieurs des coefficients a(i,j) d'une même colonne (ou d'une même ligne) sont nuls.

3. Utilisation du caractère multilinéaire

On peut suivre une méthode de type pivot de Gauss pour transformer le déterminant par des opérations élémentaires en un déterminant triangulaire dont le calcul est aisé. Le temps de calcul croît alors seulement comme n³.

La dernière méthode est donc celle qu'il convient de privilégier pour un calcul général. Mais si la matrice présente des symétries particulières, une technique adaptée peut permettre de conclure plus aisément.

[modifier] Exemple

Soit à calculer

$A = \begin{bmatrix}-2&2&-3 \\ -1& 1& 3\\ 2 &0 &-1\end{bmatrix}$

1. Par utilisation directe de la définition (ou règle de Sarrus, cela revient au même)

$\det(A)=(-2)\cdot 1 \cdot (-1) + (-3)\cdot 0 \cdot (-1) + 2\cdot 3\cdot 2 - (-3)\cdot 1 \cdot 2 - (-2)\cdot 3 \cdot 0 - 2\cdot (-1) \cdot (-1)$

$=2+0+12-(-6)-0-2 = 18.\;$

2. Par applications successives de la formule de Laplace (en commençant par la deuxième colonne, la plus avantageuse pour la disposition des zéros)

$\det(A)=(-1)^{1+2}\cdot 2 \cdot \det \begin{bmatrix}-1&3\\ 2 &-1\end{bmatrix} + (-1)^{2+2}\cdot 1 \cdot \det \begin{bmatrix}-2&-3\\ 2&-1\end{bmatrix}$

$=(-2)\cdot((-1)\cdot(-1)-2\cdot3)+1\cdot((-2)\cdot(-1)-2\cdot(-3)) = (-2)(-5)+8 = 18$