Comatrice

Restez toujours informé : suivez-nous sur Google (☆)

Introduction

En algèbre linéaire, la comatrice d'une matrice carrée A est une matrice introduite par une généralisation du calcul de l'inverse de A. Elle a une importance considérable pour l'étude des déterminants. Ses coefficients sont appelés cofacteurs de A, et ils permettent d'étudier les variations de la fonction déterminant.

La comatrice est aussi appelée matrice des cofacteurs.

Matrice ayant un coefficient variable

Le déterminant pour les matrices est naturellement défini comme une fonction sur les n vecteurs colonnes de la matrice. Il est cependant légitime de le considérer aussi comme une fonction qui aux n coefficients de la matrice associe un scalaire.

Quand on gèle tous les coefficients de la matrice à l'exception d'un seul, le déterminant est une fonction affine du coefficient variable. L'expression de cette fonction affine est simple à obtenir comme cas particulier de la propriété de n-linéarité ; elle fait intervenir un déterminant de taille n-1, appelé cofacteur du coefficient variable.

Ces considérations permettent d'établir une formule de récurrence ramenant le calcul d'un déterminant de taille n, à celui de n déterminants de taille n-1 : c'est la formule de Laplace.

Cofacteur

Soit A une matrice carrée de taille n. On observe l'effet d'une modification d'un des coefficients de la matrice, toutes choses égales par ailleurs. Pour cela on choisit donc deux indices i pour la ligne et j pour la colonne, et on note A(x) la matrice dont les coefficients sont les mêmes que ceux de A, sauf le terme d'indice i,j qui vaut a_i,j+x. On écrit la formule de linéarité pour la j-ème colonne

det A (x) = det A + x a_{1, 1} \dots a_{1, j - 1} 0 a_{1, j + 1} \dots a_{1, n} ⋮ ⋮ ⋮ ⋮ ⋮ a_{i - 1, 1} \dots a_{i - 1, j - 1} 0 a_{i - 1, j + 1} \dots a_{i - 1, n} a_{i, 1} \dots a_{i, j - 1} 1 a_{i, j + 1} \dots a_{i, n} a_{i + 1, 1} \dots a_{i + 1, j - 1} 0 a_{i + 1, j + 1} \dots a_{i + 1, n} ⋮ ⋮ ⋮ ⋮ ⋮ a_{n, 1} \dots a_{n, j - 1} 0 a_{n, j + 1} \dots a_{n, n} = det A + x C o f_{i, j}

Le déterminant noté Cof_i,j est appelé cofacteur d'indice i,j de la matrice A. Il admet les interprétations suivantes

augmenter de x le coefficient d'indice i,j de la matrice (toutes choses égales par ailleurs) revient à augmenter le déterminant de x fois le cofacteur correspondant
le cofacteur est la dérivée du déterminant de la matrice A(x)

Dans la pratique, on calcule les cofacteurs de la façon suivante : on appelle M(i;j) le déterminant de la sous-matrice déduite de M en ayant enlevé la ligne i et la colonne j (on parle de mineur pour un tel déterminant). Alors le cofacteur est (-1) fois M(i;j).

C o f i, j = (- 1)^{i + j} a 1, 1 \dots a_{1, j - 1} a_{1, j + 1} \dots a_{1, n} ⋮ ⋮ ⋮ ⋮ a_{i - 1, 1} \dots a_{i - 1, j - 1} a_{i - 1, j + 1} \dots a_{i - 1, n} a_{i + 1, 1} \dots a_{i + 1, j - 1} a_{i + 1, j + 1} \dots a_{i + 1, n} ⋮ ⋮ ⋮ ⋮ a_{n, 1} \dots a_{n, j - 1} a_{n, j + 1} \dots a_{n, n}

Formules de Laplace

Si n>1 et A est une matrice carrée de taille n alors on peut calculer son déterminant en fonction des coefficients d'une seule colonne et des cofacteurs correspondants. Cette formule, dite formule de Laplace, permet ainsi de ramener le calcul du déterminant à n calculs de déterminants de taille n-1.

Formule de développement par rapport à la colonne j

det A = i = 1 \sum n a_{i, j} C o f_{i, j}

On peut donner également une formule de développement par rapport à la ligne i

det A = j = 1 \sum n a_{i, j} C o f_{i, j}

Généralisation

On introduit la comatrice de A, matrice constituée des cofacteurs de A. On peut généraliser les formules de développement du déterminant par rapport aux lignes ou colonnes

A . (^{t} c o m A) = (^{t} c o m A) . A = det A \times I_{n}

La matrice transposée de la comatrice est appelée matrice complémentaire de A. Notamment si A est inversible, l'inverse de A est un multiple de la matrice complémentaire. Ce qui veut dire qu'on a obtenu une formule pour l'inverse, ne nécessitant que des calculs de déterminants

A^{- 1} = \frac{1}{det A}^{t} c o m A

Cette formule est encore valable si les matrices sont à coefficients dans un anneau A. Elle est utilisée pour démontrer que M est inversible en tant que matrice à coefficients dans A si et seulement si det(M) est inversible comme élément de A.

Elle est d'un intérêt limité pour calculer explicitement des inverses de matrices; en pratique elle est trop lourde dès que n=4 et la méthode plus élémentaire à base d'opérations élémentaires sur les lignes (inversion par pivot de Gauss) est plus efficace, aussi bien pour l'homme que pour la machine.

Propriétés de la comatrice

Nous avons

com(I_n) = I_n

pour toutes matrices M et N d'ordre n, com(MN) = com(N) com(M)

La comatrice est aussi compatible avec la transposition :

com(M) = (com(M)).

de plus,

det(com(M)) = det(M).

Si p(t) = det(M - tI_n) est le polynôme caractéristique de M et que q est le polynôme défini par q(t) = (p(0) - p(t))/t, alors

com(M) = q(M).

La comatrice apparaît dans la formule de la dérivée d'un déterminant.

Pour $A \in M_{n} (K)$ :

si A est de rang n (i.e. A inversible), Com(A) aussi. On a alors $C o m (A) = d e t (A)^{t} A^{- 1}$ et $C o m (A)^{- 1} = \frac{1}{d e t ( A )}^{t} A$ .
si A est de rang n-1, Com(A) est de rang 1.
si A est de rang au plus n-2, Com(A)=0.

Si $n \geq 3$ et $A \in M_{n} (K)$ , $C o m (C o m (A)) = d e t (A)^{n - 2} A$ (et est donc nulle si, et seulement si, A n'est pas inversible). Si n=2, on a Com(Com(A))=A pour toute matrice A (ce qu'on peut inclure dans la formule précédente avec la convention x = 1 pour tout $x \in K$ , y compris pour x=0).

Si $n \geq 3$ , les matrices $A \in M_{n} (R)$ telles que A=Com(A) sont la matrice nulle et les matrices spéciales orthogonales. Si n=2, ce sont les matrices multiples des matrices spéciales orthogonales.

Variations de la fonction déterminant

La formule de Leibniz montre que le déterminant d'une matrice A s'exprime comme somme et produit de composantes de A. Il n'est donc pas étonnant que le déterminant ait de bonnes propriétés de régularité. On suppose ici que K est le corps des réels.

Déterminant dépendant d'un paramètre

Si $t \mapsto A (t)$ est une fonction de classe $C^{k}$ à valeurs dans les matrices carrées d'ordre n, alors $t \mapsto det A (t)$ est également de classe $C^{k}$ .

La formule de dérivation s'obtient en faisant intervenir les colonnes de A

\frac{d}{d t} (det (A_{1} (t), \dots, A_{n} (t))) = i = 1 \sum n det (A_{1} (t), \dots, A_{i - 1} (t), A^{'} i (t), A i + 1 (t), \dots, A_{n} (t))

Cette formule est analogue formellement à la dérivée d'un produit de n fonctions numériques.

Le déterminant comme fonction sur l’espace des matrices

L’application qui à la matrice $A\$ associe son déterminant est continue (membre de $C^{0}$ ).

Cette propriété a des conséquences topologiques intéressantes :

ainsi le groupe $G L_{n} (R)$ est un ouvert,

le sous-groupe $S L_{n} (R)$ est un fermé.

Cette application est en fait différentiable (membre de $C^{1}$ ), et même infiniment déférentiable (membre de $C^{\infty}$ ).

En effet le calcul des cofacteurs peut être vu précisément comme un calcul de dérivée partielle

\frac{\partial det}{\partial E _{i, j}} (A) = C o f A_{i, j} .

Toutes ces dérivées partielles étant elles-mêmes des déterminants, par récurrence le déterminant est $C^{\infty}$ .

En outre on peut écrire le développement limité à l’ordre un du déterminant au voisinage de $A\$

det (A + H) = det A + t r (^{t} C o m (A) . H) + o (∣ H ∣),

c’est-à-dire que si on munit $M_{n} (R)$ de son produit scalaire canonique, l’application déterminant a pour gradient

\nabla det (A) = C o m (A) .

Notamment pour le cas où $A\$ est l’identité,

det (I + H) = 1 + t r (H) + o (∣ H ∣),

\nabla det (I) = I .

Comatrice et produit vectoriel

Si A est une matrice d'ordre trois, elle agit sur les vecteurs de l'espace à trois dimensions muni d'une base orthonormée d'orientation directe. La comatrice de A décrit alors l'interaction de A avec le produit vectoriel:

$A u \land A v = C o m A (u \land v)$ .