- Formule de Laplace
-
Comatrice
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, ou encore, matrice adjointe (par exemple dans le logiciel Maple), à ne pas confondre avec la matrice adjointe au sens de transposée de la matrice conjuguée.
Sommaire
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 n2 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 ai,j+x. On écrit la formule de linéarité pour la j-ème colonne
Le déterminant noté Cofi,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)i+j fois M(i;j).
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
- On peut donner également une formule de développement par rapport à la ligne i
Démonstration de la formule de LaplaceQuitte à transposer la matrice, il suffit de prouver la formule du développement par rapport à une colonne.
On considère la matrice M0 obtenue en remplaçant la colonne j de la matrice A par une colonne de 0. Le déterminant de M0 est nul.
On passe de M0 à A en modifiant successivement les différents coefficients de la colonne j. On fait d'abord passer le premier coefficient de 0 à a1;j, puis le second de 0 à a2;j et ainsi de suite. La première opération revient à ajouter a1;j multiplié par son cofacteur, la seconde a1;j multiplié par son cofacteur.
Au bout du compte, pour passer du déterminant de M0 (nul) à celui de A, on a ajouté successivement tous les termes intervenant dans la formule de Laplace.
Remarque : pour la correction de la preuve, les cofacteurs à chaque étape sont bien les mêmes que ceux de la matrice A, puisqu'on ne modifie que la colonne 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
DémonstrationOn prend de nouveau pour matrice M0 la matrice obtenue en remplaçant la colonne j de la matrice A par une colonne de 0. Le déterminant de M0 est nul.
Mais cette fois au lieu de faire apparaître en colonne j les coefficients ai;j, on fait apparaître successivement les coefficients ai;k d'une autre colonne. La matrice finale obtenue admet alors deux colonnes identiques, donc est elle aussi de déterminant nul. Ceci signifie
En effet par rapport à la démonstration précédente, si les coefficients ont changé, les cofacteurs non.
Si on ajoute à cette étude le cas j=k qui a été examiné dans le paragraphe précédent, on arrive à
qui est exactement, composante par composante, la formule souhaitée.
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
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(In) = In
et
- pour toutes matrices d'ordre n M et N, com(MN) = com(M) com(N)
La comatrice est aussi compatible avec la transposition :
- com(tM) = t (com(M)).
de plus,
- det(com(M)) = det(M)n-1.
Si p(t) = det(M - tIn) 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
- tcom(M) = q(M).
La comatrice apparaît dans la formule de la dérivée d'un déterminant.
Pour :
- si A est de rang n (i.e. A inversible), Com(A) aussi. On a alors et .
- 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 et , (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 x0 = 1 pour tout , y compris pour x=0).
Si , les matrices 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 est une fonction de classe à valeurs dans les matrices carrées d'ordre n, alors est également de classe .
La formule de dérivation s'obtient en faisant intervenir les colonnes de A
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 associe son déterminant est continue (membre de ).
- Cette propriété a des conséquences topologiques intéressantes :
- ainsi le groupe est un ouvert,
- le sous-groupe est un fermé.
- Cette application est en fait différentiable (membre de ), et même infiniment déférentiable (membre de ).
- En effet le calcul des cofacteurs peut être vu précisément comme un calcul de dérivée partielle
- Toutes ces dérivées partielles étant elles-mêmes des déterminants, par récurrence le déterminant est .
- En effet le calcul des cofacteurs peut être vu précisément comme un calcul de dérivée partielle
- En outre on peut écrire le développement limité à l’ordre un du déterminant au voisinage de
-
- c’est-à-dire que si on munit de son produit scalaire canonique, l’application déterminant a pour gradient
- Notamment pour le cas où est l’identité,
-
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:
- .
Voir aussi
- Portail des mathématiques
Catégories : Déterminant | Matrice
Wikimedia Foundation. 2010.