Matrice carrée
matrice avec autant de lignes que de colonnes
From Wikipedia, the free encyclopedia
En mathématiques, une matrice carrée est une matrice avec autant de lignes que de colonnes.
Lorsque les coefficients sont pris dans un anneau unitaire R, l’ensemble des matrices carrées de taille n (c’est-à-dire avec n lignes et n colonnes) constitue également un anneau non commutatif, noté , dont la matrice identité est l’élément unité. L’ensemble des matrices diagonales constitue un sous-anneau isomorphe à Rn.
En plus du rang et de la transposition existant sur toutes les matrices, chaque matrice carrée a une trace égale à la somme de ses éléments diagonaux.
Si R est commutatif, l’anneau admet même une structure d’algèbre associative, munie de la forme multilinéaire du déterminant qui permet notamment de définir le polynôme caractéristique d’une matrice et de caractériser les matrices inversibles.
Si les coefficients sont pris dans un corps commutatif k, les matrices de représentent les endomorphismes de kn ou de n’importe quel espace vectoriel sur k de dimension n muni d’une base. Leurs valeurs propres et vecteurs propres permettent d'en obtenir une forme réduite (voire une diagonalisation) par conjugaison avec une matrice inversible, ce qui facilite l’expression des puissances voire de leur exponentielle.
De multiples classes de matrices peuvent être définies parmi les matrices carrées, notamment lorsque les coefficients sont entiers, réels ou complexes : matrices triangulaires, symétriques ou antisymétriques, orthogonales, hermitiennes, unitaires… Différentes normes d’algèbre peuvent être définies pour enrichir la structure de ou .
Certains résultats peuvent s’exprimer dans le contexte plus large de coefficients pris dans un pseudo-anneau (anneau sans unité) ou un demi-anneau (sans symétrique pour la première loi, par exemple une algèbre de Boole).
Vocabulaire de forme
Dans une matrice carrée, la taille décrit le nombre de lignes qui est égal au nombre de colonnes.
Les coefficients diagonaux sont ceux dont le numéro de ligne est égal au numéro de colonne. Lorsque la matrice est présentée comme un tableau, ils sont lisibles sur la diagonale principale qui rejoint le coin supérieur gauche au coin inférieur droit.
Les coefficients sur-diagonaux sont ceux dont le numéro de ligne est strictement inférieur au numéro de colonne et sont lisibles au-dessus de la diagonale principale. Les coefficients sous-diagonaux sont à l’inverse ceux dont le numéro de ligne est strictement supérieur au numéro de colonne, lisibles en-dessous de la diagonale principale.
Une matrice diagonale est une matrice carrée dont tous les coefficients sont nuls en dehors de la diagonale principale. On la note à l’aide de l’opérateur Diag et la liste de ses coefficients diagonaux :
Une matrice scalaire est une matrice diagonale dont les tous les coefficients diagonaux sont égaux. En particulier, la matrice identité (de taille n) est la matrice scalaire dont tous les coefficients diagonaux valent 1, notée In ou I.
Une matrice triangulaire supérieure est une matrice carrée dont tous les coefficients sous-diagonaux sont nuls. Une triangulaire inférieure a tous ses coefficients sur-diagonaux nuls.
Structure algébrique
Addition et multiplication
Si R est un corps comme ℝ ou ℂ, ou plus généralement un anneau unitaire quelconque, et si n un entier strictement positif, l’ensemble des matrices carrées de taille n à coefficients dans R est muni d’une opération d’addition des coefficients terme à terme, et d'une opération de multiplication interne induite par le produit matriciel. La matrice identité I est neutre pour cette multiplication : pour toute matrice carrée A de taille n,
- .
Cette multiplication est associative et distributive par rapport à l’addition, mais même si R est commutatif, dès lors que n > 1 la multiplication des matrices carrées n’est pas commutative : il existe des matrices A et B telles que AB ≠ BA.
Plus formellement, l’addition et la multiplication induisent une structure d’anneau unitaire qui contient l’ensemble des matrices scalaires comme sous-anneau isomorphe à R, et même l’ensemble des matrices diagonales comme sous-anneau isomorphe à Rn.
Inverse
Les matrices carrées de même taille à coefficients dans un même anneau R peuvent toutes être multipliées entre elles mais cela n’implique pas qu’on puisse diviser par une matrice. En effet, il est possible d’avoir une égalité AB = AC avec des matrices carrées de même taille et A non nulle sans que B soit égal à C.
Si deux matrices carrées A et B de même taille vérifient les égalités AB = BA = I, on dit qu'elles sont inverses l'une de l'autre. Comme une matrice ne peut pas avoir plusieurs inverses distinctes[1], on peut noter B = A−1. À défaut de pouvoir diviser par A, on peut alors multiplier par son inverse, ce qui permet de résoudre certaines équations linéaires
La détermination de l'inversibilité d'une matrice et le calcul de son inverse le cas échéant sont obtenus par des techniques spécifiques dépendant de la structure de l'anneau de coefficients, de la taille de la matrice et éventuellement le degré de précision souhaité.
L'ensemble des matrices inversibles dans forme un groupe appelé groupe linéaire et noté .
Multiplication scalaire
La multiplication scalaire (à gauche ou à droite) d'une matrice M par un élément λ de R revient à multiplier M (du même côté) par la matrice scalaire dont les coefficients diagonaux valent tous λ (et qui peut se noter λ.I) :
- à gauche, ;
- à droite, .
Si l’anneau de coefficients est commutatif, ou plus généralement si l’élément λ est central dans R, les deux multiplications scalaires à gauche et à droite coïncident.
Réciproquement, si une matrice commute avec toutes les autres alors elle est de la forme λ.I avec λ central dans R.
L’addition et la multiplication scalaire induisent une structure de R-module (voire de bimodule) libre isomorphe à .
Si l'anneau R est commutatif, ces opérations font de une R-algèbre associative unitaire.
Puissances et polynômes
Une matrice carrée peut être multipliée par elle-même et la répétition de ces multiplications définit la suite de ses puissances par récurrence :
Pour une même matrice A, toutes ses puissances commutent entre elles.
Si l'anneau de coefficients R est commutatif, par combinaison linéaire de ces puissances on peut définir n'importe quel polynôme de matrice P(A) avec et . Dans ce cas, tous les polynômes d'une même matrice commutent entre eux.
Plus généralement, si deux matrices commutent, tout polynôme de l'une commute avec tout polynôme de l'autre, et dans ce cas les identités remarquables du second degré usuelles sont satisfaites, ainsi que leurs généralisation pour des exposants supérieurs, qui sont la formule du binôme de Newton et la formule de Bernoulli. En particulier, puisque toute matrice commute avec une matrice scalaire :
Un polynôme annulateur d'une matrice A est un polynôme P tel que P(A) = 0.
Représentation d'endomorphisme
Toute matrice carrée représente un endomorphisme de R-module à droite (voire d'espace vectoriel si R est un corps) par
Réciproquement, tout endomorphisme de R-module libre à droite est ainsi représenté par une matrice dont les colonnes correspondent aux images des éléments de la base canonique de la forme (0, ... , 0, 1, 0, ... , 1).
La multiplication matricielle traduit alors la composition des endomorphismes.
Une matrice carrée est idempotente (c'est-à-dire égale à son carré) si et seulement si elle représente un projecteur.
Elle est involutive (c'est-à-dire égale à sa propre inverse) si elle représente une symétrie vectorielle.
Une matrice carrée est dite nilpotente si l'une de ses puissances est la matrice nulle.
Similitude et invariants
Matrices semblables
Deux matrices A et B de sont dites semblables s'il existe une matrice inversible P telle que
- .
Cette relation de similitude définit une relation d'équivalence sur les matrices. Elle exprime le fait que les matrices A et B représentent le même endomorphisme de Rn dans des bases différentes. La matrice P correspond alors à une matrice de passage de la deuxième base vers la première.
Deux matrices semblables sont forcément équivalentes et ont donc le même rang (lorsque ce dernier est bien défini, par exemple si l'anneau de coefficients est commutatif ou noethérien) mais la réciproque est fausse en général.
La description d'un système complet d'invariants (caractérisant des matrices semblables) est plus délicate. On appelle ces invariants les invariants de similitude. D'un point de vue algorithmique, la réduction d'une matrice quelconque à une matrice sous une forme privilégiée se fait par un algorithme inspiré de celui du pivot de Gauss, voir théorème des facteurs invariants.
Opérateurs caractéristiques
Dans toute cette section, l'anneau R des coefficients sera supposé commutatif.
Trace
La trace d'une matrice carrée est la somme de ses coefficients diagonaux : Si ,
- .
Cet opérateur linéaire a la particularité d'être invariant par interversion des facteurs d'un produit de deux matrices (qu'elles soient carrées ou rectangulaires) :
- .
Déterminant
Le déterminant d'une matrice carrée peut être défini par une somme indexée par le groupe symétrique et utilisant la fonction signature ε des permutations sur l'ensemble :
La fonction déterminant est multiplicative, c'est-à-dire qu'elle vérifie pour toutes les matrices A et B :
- .
En outre, le déterminant de l'identité vaut 1, donc une matrice inversible dans aura toujours un déterminant inversible dans R avec
- .
Réciproquement, toute matrice avec un déterminant inversible dans R est inversible dans et son inverse peut s'écrire avec la transposée de la comatrice :
- .
Polynôme caractéristique
Le déterminant permet de définir le polynôme caractéristique de degré n et unitaire dans R[X] par
Ce polynôme est un invariant de similitude : si alors
- .
Par conséquent, tous ses coefficients sont également invariants, et en particulier son coefficient constant qui est (au signe près) le déterminant de la matrice et son coefficient sous-dominant (en degré n−1) qui est l'opposé de la trace.
Le polynôme caractéristique constitue aussi un polynôme annulateur de la matrice en vertu du théorème de Cayley-Hamilton.
Valeurs propres, vecteurs propres et espaces propres
Soit A est une matrice carrée de taille n à coefficients dans un anneau commutatif R, avec X un vecteur colonne non nulle de taille n et λ dans R. On dit que X est un vecteur propre de A pour la valeur propre λ si
- .
En particulier, cela signifie que le vecteur X appartient au noyau de la matrice , qu'on appelle espace propre associé à la valeur propre λ, et qu'on note :
- .
Si R est un corps commutatif, le noyau est non trivial si et seulement si le déterminant est égal à zéro, ce qui revient à dire que λ est une racine du polynôme caractéristique. Il y a donc au maximum n valeurs propres, et chaque vecteur propre n'est associé qu'à une seule valeur propre. Toute famille de vecteurs propres associés à des valeurs propres deux à deux distinctes est libre et on montre même que les espaces propres sont en somme directe.
Sans cette hypothèse, le noyau est non trivial si et seulement si le déterminant est égal à zéro ou est un diviseur de zéro[2]. Il peut exister plus de n valeurs propres, et un même vecteur propre peut être associé à plusieurs valeurs propres. Par exemple, avec et , la matrice a 4 valeurs propres :
- 2 pour tous les vecteurs non nuls ;
- 5 pour les vecteurs 2 et 4 ;
- 0 et 4 pour le vecteur 3 uniquement.
dont le spectre ne contient aucun entier,
contient une seule valeur propre réelle
et deux valeurs propres supplémentaires dans ℂ.
Dans tous les cas, l'ensemble des valeurs propres d'une matrice carrée A forme son spectre, noté Sp(A) ou σ(A), qui est invariant par similitude mais dépend de l'anneau de coefficients choisi : pour une même matrice à coefficients entiers, le spectre peut contenir des valeurs propres réelles non entières et même des valeurs propres supplémentaires dans l'ensemble des nombres complexes.
Réduction
- Réduire une matrice, c'est trouver une matrice qui lui est semblable la plus simple possible.
- Une matrice diagonalisable est une matrice semblable à une matrice diagonale : A est diagonalisable s'il existe une matrice inversible P et une matrice diagonale D telles que A = P−1DP.
- Sur un corps algébriquement clos, on dispose de la réduction de Jordan qui est optimale et il existe des décompositions intermédiaires comme la décomposition de Dunford qui utilise les sous-espaces caractéristiques ou celle de Frobenius qui utilise les sous-espaces cycliques.
- Les polynômes d'endomorphismes jouent un rôle crucial dans les techniques de réduction.
Matrice normale
Une matrice normale est une matrice à coefficients réels qui commute avec sa transposée, ou plus généralement une matrice carrée à coefficients complexes qui commute avec son adjointe, qui est la conjuguée de sa transposée :
- ;
Symétrique ou antisymétrique
Une matrice carrée à coefficients réels est dite symétrique si elle est égale à sa transposée. Elle est dite antisymétrique si elle est opposée à sa transposée.
L'ensemble des matrices symétriques Sn(ℝ) et l'ensemble des matrices antisymétriques An(ℝ) forment deux sous-espaces supplémentaires orthogonaux de , ce qui implique notamment que toute matrice carrée réelle se décompose d'une unique manière en somme d'une matrice symétrique et d'une matrice antisymétrique
- .
Orthogonales
Une matrice carrée à coefficients réels est dite orthogonale si elle est l'inverse de sa transposée.
L'ensemble des matrices orthogonales de forme un sous-groupe du groupe linéaire, appelé groupe orthogonal et noté .
L'ensemble des matrices orthogonales de déterminant 1 forme un sous-groupe du précédent, appelé groupe spécial orthogonal et noté .
Hermitiennes et unitaires
Une matrice carrée à coefficients complexes est dite hermitienne si elle est égale à sa matrice adjointe.
Elle est dite unitaire si elle est inverse de sa matrice adjointe.
L'ensemble des matrices unitaires de taille n forme un groupe appelé groupe unitaire et noté Un.
Interprétations bilinéaires
Dans ce paragraphe, l'anneau K des scalaires sera supposé commutatif. Dans la plupart des applications, ce sera un corps commutatif.
Le cas non commutatif existe aussi mais il faut prendre quelques précautions et les notations deviennent trop lourdes pour cet article.
Matrice d'une forme bilinéaire
Soient E un K-module libre et B = (e1, … , en) une base de E.
Soit une forme bilinéaire. On définit la matrice de dans la base B par la formule suivante :
Dans le cas particulier où K = ℝ et f est un produit scalaire, cette matrice est appelée matrice de Gram.
La matrice matB f est symétrique (respectivement antisymétrique) si et seulement si la forme bilinéaire est symétrique (respectivement antisymétrique).
Soient x et y deux vecteurs de E. Notons X et Y leurs coordonnées dans la base B et A = matB f. On a alors la formule : .
Deux formes bilinéaires sont égales si et seulement si elles ont la même matrice dans une base donnée.
Matrice d'une forme quadratique
Lorsque K est un corps de caractéristique différente de 2, on appelle matrice d'une forme quadratique la matrice de la forme bilinéaire symétrique dont est issue la forme quadratique.
Formule de changement de base
Soient E un K-module libre et B, C deux bases de E. Soit une forme bilinéaire.
Notons M = matB f la matrice de f dans la base B et M' = matC f la matrice de f dans la base C. Notons P = matB C la matrice de passage. On a alors la formule de changement de base pour une forme bilinéaire (à ne pas confondre avec celle pour une application linéaire) :
Matrices congruentes
Deux matrices carrées A et B sont dites congruentes s'il existe une matrice inversible P telle que
Deux matrices congruentes sont deux matrices qui représentent la même forme bilinéaire dans deux bases différentes.
Lorsque K est un corps de caractéristique différente de 2, toute matrice symétrique est congruente à une matrice diagonale. L'algorithme utilisé s'appelle réduction de Gauss à ne pas confondre avec le pivot de Gauss.
Décomposition d'une matrice
On utilise abusivement le terme décomposition d'une matrice, qu'il s'agisse d'une véritable décomposition (en somme) comme dans la décomposition de Dunford ou d'une factorisation comme dans la plupart des autres décompositions.
Décomposition LU
- C'est une factorisation en produit de deux matrices triangulaires.
- En lien avec le pivot de Gauss, c'est une méthode qui permet d'inverser une matrice.
Décomposition QR
- C'est un résultat sur les matrices à coefficients réels ou à coefficients complexes.
- C'est une factorisation en produit d'une matrice orthogonale et d'une matrice triangulaire.
- C'est une traduction matricielle du procédé de Gram-Schmidt.
Décomposition polaire
- C'est un résultat sur les matrices à coefficients réels ou à coefficients complexes.
- C'est une factorisation en produit d'une matrice orthogonale et d'une matrice symétrique définie positive dans le cas réel, en produit d'une matrice unitaire et d'une matrice hermitienne définie positive dans le cas complexe.
- On peut décomposer à droite ou à gauche.
- On a unicité de la factorisation pour les matrices inversibles.
Normes d'algèbre
Dans tout ce paragraphe, K = ℝ ou ℂ.
Une norme matricielle est une norme d'algèbre sur l'algèbre Mn(K), c'est-à-dire une norme d'espace vectoriel qui est de plus sous-multiplicative.
Le rayon spectral d'une matrice carrée A à coefficients complexes est le plus grand module de ses valeurs propres. Il est égal à la borne inférieure des normes matricielles de A.
Sur Mn(K), toute norme N subordonnée à une norme sur Kn est une norme d'algèbre vérifiant de plus N(In) = 1 (la réciproque est fausse).
Extension de fonction d'une matrice
On peut étendre aux matrices carrées réelles voire complexes certaines fonctions d'une variable réelle qui peuvent être définies comme des séries entières, avec ou non, selon les fonctions, des limitations concernant la convergence de la série.
Soit A ∈ Mn(ℂ), soit N une norme d'algèbre et une série entière de rayon de convergence R.
Alors si N(A) < R, la série est absolument convergente. (On le montre en utilisant que N(An) ≤ N(A)n.)
Exponentielle d'une matrice
Pour toute matrice carrée réelle ou complexe A, on peut définir son exponentielle par :
- , une série qui est toujours convergente.
Le calcul effectif de cette exponentielle se fait par réduction de la matrice.
L'exponentielle joue un rôle central dans l'étude des systèmes linéaires d'équations différentielles.