- Calcul du déterminant d'une matrice
-
Le calcul du déterminant d'une matrice est un outil nécessaire tant en algèbre linéaire pour vérifier une inversibilité ou calculer l'inverse d'une matrice qu'en analyse vectorielle avec, par exemple, le calcul d'un jacobien. S'il existe une formule générale de calcul du déterminant, sa complexité en fait une technique difficile à mettre en œuvre pour des matrices de grande taille. On lui préfère alors des méthodes de calculs plus simples comme la technique du pivot de Gauss.
Certaines matrices de forme particulière ont des déterminants déjà étudiés.
Sommaire
Présentation
Article détaillé : Déterminant (mathématiques).Le déterminant de la matrice carrée est donné par la formule de Leibniz
où désigne l'ensemble des permutations de et ε(σ) la signature de la permutation σ.
Il s'agit donc d'effectuer tous les produits possibles en prenant un élément par ligne et par colonne dans la matrice, de les multiplier tantôt par +1 tantôt par -1[1], et de faire la somme des n! termes ainsi obtenus
Soit à calculer, par exemple, le déterminant de
Il y a six produits à calculer en prenant un terme par ligne et par colonne. Le produit (-2)(1) (-1) est précédé de + car, dans toutes les paires, le terme de gauche est au dessus de celui de droite, le produit (-2)(0)(3) est précédé du signe - car il existe une seule paire, la paire {0;3}, où le terme de gauche est sous le terme de droite, le produit (-1)(2)(-1) précédé de - car il existe une seule paire, {-1;2}, où le terme de gauche est sous celui de droite, le produit (-1)(0)(-3) précédé de +à cause des paires {-1;-3} et {0;-3}, le produit (2)(2)(3) précédé de + à cause des paires {2;2}et {2;3}, et le produit (2)(1)(-3) précédé de - à cause des trois paires {2;1}{2;3} et {1;-3}
On peut aussi calculer le déterminant d'une matrice de dimension n à l'aide de n déterminants de matrices de dimension n - 1 obtenues en enlevant à la matrice de départ une ligne et une colonne. Si A est la matrice, pour tout i et j, on note Ai,j la matrice obtenue en enlevant à A sa i-ème ligne et sa j-ème colonne.
On peut alors développer le calcul du déterminant de A suivant une ligne ou une colonne.
Développement suivant la ligne i :
Le terme ( − 1)i + jdet(Ai,j) est appelé le cofacteur du terme ai,j et le terme det(Ai,j) est appelé le mineur du terme ai,j . Cette méthode porte le nom de développement suivant une ligne (ou une colonne)[2] , méthode de Laplace[3] ou méthode des cofacteurs[4] ou des mineurs[5].
Ainsi la matrice précédente se développe aisément suivant la deuxième colonne, la plus avantageuse pour la disposition des zéros.
Exemples en dimension 2 et 3
Dimension 2
Le déterminant de la matrice possède deux termes : le produit ad précédé de + car le terme le plus à gauche (a) est situé au dessus de terme le plus à droite (d) et le produit bc précédé du signe - car le terme de gauche (c) est situé sous le terme de droite (b)
Dimension 3
Il y a 6 façons de choisir trois termes un par ligne et par colonnes, il y a donc 6 produits dans un déterminant d'ordre 3, 3 sont précédé du signe + et trois sont précédés du signe -. 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 , il suffit d'effectuer
a b c d e f g h i a b c d e f et 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 .
Techniques de simplification du calcul d'un déterminant
Le calcul du déterminant d'une matrice carrée de dimension n nécessite le calcul d'autant de produits que de permutations à n éléments c'est à dire n! produits à effectuer, soit 2 pour une matrice carrée, 6 pour une matrice de dimension 3 et 24 pour une matrice de dimension 4. De plus, il s'agit de trouver la signature de chacune des permutations. Le développement suivant une ligne ou une colonne permet d'organiser plus clairement les calculs mais ne diminue en rien le nombre de produits à effectuer.
On remarque cependant que la présence d'un zéro dans une des case de la matrice permet de faire disparaitre (n-1)! calculs. L'idée est donc de trouver des techniques remplaçant le calcul du déterminant d'une matrice par celui d'une matrice contenant de nombreux zéros, dite matrice à trous. On dispose pour cela d'un certain nombre de propriétés opératoires et de quelques techniques
Propriétés opératoires élémentaires
Le déterminant est une forme n-linéaire alternée des vecteurs colonnes ou des vecteurs lignes. Cette propriété a les conséquence suivantes
- si on permute deux lignes ou deux colonnes, le déterminant change de signe
- si deux lignes ou deux colonnes sont identiques, le déterminant est nul
- si on multiplie tous les termes d'une même ligne ou d'une même colonne par un réel k, le déterminant est multiplié par k.
- on peut ajouter à une colonne (ou une ligne) le multiple d'une autre colonne (ou d'une autre ligne) sans changer la valeur du déterminant
Enfin le déterminant se comporte bien avec le produit des matrices :
- det (A \times B)=det(A).det(B)
Matrice triangulaire ou triangulaire par bloc
- le déterminant d'une matrice triangulaire est le produit des coefficients diagonaux
Déterminant d'une matrice Triangulaire - 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 de dimension n est le produit des déterminants des blocs diagonaux
- où
Déterminant d'une matrice par blocs - Démonstration par développements successifs- On commence par simplifier la situation en utilisant le produit par blocs suivant
- 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.
Déterminant d'une Matrice par blocs - Démonstration grâce aux fonctions symétriques élémentairesNotons .
Alors .
Or [[1;p]] est stable par σ et comme σ est une bijection de [[1;n]], l'ensemble [[p;n]] est aussi stable.
Donc
Méthode du pivot de Gauss
Article détaillé : Pivot de Gauss.Cette méthode consiste à remplacer la matrice par une matrice triangulaire en utilisant seulement des permutations de lignes ou colonnes et des ajouts à un ligne d'un multiple d'une autre ligne de manière à faire apparaitre un maximum de zéros.
Le principe est le suivant:
- on choisit dans la matrice un terme non nul ai,j, en général le premier terme en haut à gauche, que l'on appelle le pivot.
- si le terme choisi n'est pas a1,1, on peut, en permutant les lignes 1 et i et les colonnes 1 et j le mettre à la bonne position. On obtient alors une matrice A' telle que: : det(A) = ( − 1)i + jdet(A')
- on élimine tous les termes situés sous le pivot , a1,1 en ajoutant à la ligne k la ligne 1 multipliée par . Cette opération ne change pas la valeur du déterminant
- on recommence ensuite le même processus dans la sous matrice privée de sa première ligne et de sa première colonne
- on obtient alors à la dernière étape une matrice triangulaire dont le déterminant est égal, au signe près, au déterminant de la matrice de départ.
Ainsi, dans la matrice , on peut choisir( -2) comme premier pivot et ajouter ainsi à la seconde ligne, la première multipliée par -1/2 et ajouter à la troisième ligne la première ligne
en choisissant 2 comme second pivot et en permutant les lignes 2 et 3, ce qui conduit à multiplier par (-1) le déterminant, on obtient directement une matrice triangulaire.
Cas particuliers de déterminant
Déterminant de Vandermonde
Article détaillé : Matrice de Vandermonde.Le déterminant de Vandermonde est le déterminant d'une matrice dans laquelle chaque ligne est composée des premières puissances d'un même nombre. Si les coefficients sont dans un corps (ou un anneau intègre), ce determinant s'annule si et seulement si deux lignes sont identiques.
Déterminant circulant
Article détaillé : Matrice circulante.Un déterminant circulant droit[6] est le déterminant d'une matrice dont les lignes sont obtenues par permutations circulaires des éléments de la première ligne. Supposons donnée la famille de complexes :
Soit le polynôme dont les coefficients sont donnés par la famille :
et soit la première racine n-ième de l'unité :
Le déterminant circulant s'exprime à l'aide de et de la manière suivante :
Déterminant d'une matrice tridiagonale
Article détaillé : Matrice tridiagonale.Une matrice tridiagonale est une matrice à trous contenant des zéros sauf éventuellement sur la première diagonale ainsi que les deux sous-diagonales limitrophes supérieure et inférieure. Le déterminant d'une telle matrice se calcule par récurrence à l'aide des sous-matrices tridiagonales obtenues en ne conservant que les k premières lignes et les k premières colonnes. Si on appelle A la matrice définie par :
On peut développer le déterminant par récurrence en:
Déterminant d'une matrice de Hessenberg
Article détaillé : Matrice de Hessenberg.Une matrice de Hessenberg est une matrice quasi-triangulaire. Dans une matrice de Hessenberg supérieure, tous les termes situés sous la diagonale sont nuls sauf éventuellement ceux situés sur la première sous-diagonale. À ce titre, une matrice tridiagonale est une matrice de Hessenberg à la fois supérieure et inférieure. Le déterminant d'une matrice de Hessenberg inférieure se calcule par récurrence selon une technique voisine de celle utilisée pour le calcul du déterminant tridiagonal. En appelant les sous-matrices de Hessenberg obtenues en ne conservant que les k premières lignes et les k premières colonnes, on a[7]:
Déterminant de Sylvester
Article détaillé : Résultant.Soient Pet Q sont deux polynômes de degrés respectifs n et m tels que :
On appelle déterminant de Sylvester ou résultant des polynômes P et Q, le déterminant de la matrice de Sylvester de dimension n + m
Si on se place dans un corps dans lequels les deux polynômes sont scindés, c'est-à-dire qu'ils se décomposent en produit de polynômes du premier degré:
on a :
Déterminant de Cauchy et de Hilbert
Article détaillé : Déterminant de Cauchy.Soient et deux familles de complexes tels que, pour tout i et j, ., le déterminant de Cauchy associé à ces deux familles est le déterminant de la matrice de terme général .
Il a pour expression
En particulier, si et , le déterminant obtenu est le déterminant de Hilbert dont il existe la formule explicite suivante[8]:
avec la notation :
Calcul de déterminant et complexité
Pour des calculs par ordinateur, il est important de connaitre le coût d'un calcul, c'est à dire le nombre d'opérations nécessaires pour le réaliser. La méthode de Laplace nécessite un nombre d'opérations proportionnel à n! , on dit qu'il est de complexité O(n!)[9].
L'utilisation d'une méthode de pivot de Gauss demande la précaution de ne pas diviser par 0. Si la matrice est suffisamment régulière pour que le choix du pivot soit naturellement sur la diagonale, le nombre d'opérations est majoré[10] par un nombre proportionnel à n3. Si pour des calculs à la main, le choix se porte sur des pivots simples (proches de 1), en analyse numérique, il est souvent préférable de choisir pour pivot des nombres grands en valeur absolue pour minimiser les erreurs commises dans le calcul des quotients. Enfin, si l'on tient à donner le résultat sous forme exacte fractionnaire, il faut aussi tenir compte de la taille des nombres manipulés. Dans ce cas, d'autres méthodes se révèlent intéressantes comme la méthode de Jordan-Bareiss[11] ou la méthode de Dogson[12]
Notes et références
- Cette affectation est difficile et fait intervenir le nombre d'inversion de la permutation, c'est à dire le nombre de paires parmi les termes du produit où l’élément de gauche dans la matrice est situé plus bas que l'élément de droite. Si ce nombre est impair, le produit est multiplié par -1, sinon il est multiplié par +1
- p. 481 Stéphane Balac,Frédéric Sturm, Algèbre et analyse: cours de mathématiques de première année avec exercices,
- p.484 Arthur Adam,Francis Lousberg, Espace Math 56,
- p.44 , école polytechnique fédérale de Lausanne Lara Thomas, Algèbre Linéaire
- Calcul du déterminant d'une matrice par la méthode des mineurs M; Fouquet,
- p.325 Daniel Guinin, Bernard Joppin, Algèbre et géométrie MPSI,
- p.75 Jounaïdi Abdeljaoued, Henri Lombardi, Méthodes matricielles: introduction à la complexité algébrique,
- The Hilbert Matrix and its Determinant Robert Fossum,
- p.30 Alfio Quarteroni,Fausto Saleri,Paola Gervasio, Calcul scientifique. Cours, exercices corrigés et illustrations en matlab
- p60 Jounaïdi Abdeljaoued,Henri Lombardi, Méthodes matricielles: introduction à la complexité algébrique,
- p66 Jounaïdi Abdeljaoued,Henri Lombardi, Méthodes matricielles: introduction à la complexité algébrique,
- p71 Jounaïdi Abdeljaoued,Henri Lombardi, Méthodes matricielles: introduction à la complexité algébrique,
Voir aussi
-
Wikimedia Foundation. 2010.