Matrice (mathématiques)

Bien que le calcul matriciel proprement dit n'apparaisse qu'au début du XIX^e siècle, les matrices, en tant que tableaux de nombres, ont une longue histoire d'applications à la résolution d'équations linéaires. Le texte chinois Les Neuf Chapitres sur l'art mathématique, écrit vers le II^e siècle av. J.-C., est le premier exemple connu de l'utilisation de tableaux pour résoudre des systèmes d'équations[1], introduisant même le concept de déterminant. En 1545, Jérôme Cardan fait connaître cette méthode en Europe en publiant son Ars Magna[2]. Le mathématicien japonais Seki Kōwa utilise indépendamment les mêmes techniques pour résoudre des systèmes d'équations en 1683[3]. Aux Pays-Bas, Johan de Witt représente des transformations géométriques à l'aide de tableaux dans son livre de 1659, Elementa curvarum linearum[2]. Entre 1700 et 1710, Leibniz montre comment utiliser les tableaux pour noter des données ou des solutions, et expérimente plus de 50 systèmes de tableaux à cet effet[2]. En 1750, Gabriel Cramer publie la règle qui porte son nom[4].

En 1850, le terme de « matrix » (qui sera traduit par matrice) est forgé (sur la racine latine mater) par James Joseph Sylvester[5], qui le voit comme un objet donnant naissance à la famille de déterminants actuellement appelés mineurs, c'est-à-dire les déterminants des sous-matrices obtenues en retirant des lignes et des colonnes. Dans un article de 1851, Sylvester précise :

« Dans des articles antérieurs, j'ai appelé matrix un tableau rectangulaire de termes à partir desquels plusieurs systèmes de déterminants peuvent être engendrés, comme issus des entrailles d'un parent commun »[6].

En 1854, Arthur Cayley publie un traité sur les transformations géométriques utilisant les matrices de façon beaucoup plus générale que tout ce qui a été fait avant lui. Il définit les opérations usuelles du calcul matriciel (addition, multiplication et division) et montre les propriétés d'associativité et de distributivité de la multiplication[2]. Jusque-là, l'utilisation des matrices s'était essentiellement limitée au calcul des déterminants ; cette approche abstraite des opérations sur les matrices est révolutionnaire. En 1858, Cayley publie son A Memoir on the Theory of Matrices[7] - [8], dans lequel il énonce et démontre le théorème de Cayley-Hamilton[2] pour les matrices 2×2.

Beaucoup de théorèmes ne sont d'ailleurs démontrés au début que pour de petites matrices : après Cauchy, Hamilton généralise le théorème aux matrices 4×4, et ce n'est qu'en 1898 que Frobenius, étudiant les formes bilinéaires, démontre le théorème en dimension quelconque. C'est aussi à la fin du XIX^e siècle que Wilhelm Jordan établit la méthode d'élimination de Gauss-Jordan (généralisant la méthode de Gauss pour les matrices échelonnées). Au début du XX^e siècle, les matrices occupent une place centrale en algèbre linéaire[9], en partie grâce au rôle qu'elles jouent dans la classification des systèmes de nombres hypercomplexes du siècle précédent.

Un mathématicien anglais du nom de Cullis est le premier, en 1913, à utiliser la notation moderne des crochets (ou des parenthèses) pour représenter les matrices, ainsi que de la notation systématique A = [a_i,j] pour représenter la matrice dont a_i,j est le terme de la i-ème ligne et de la j-ème colonne[2].

La formulation de la mécanique quantique au moyen de la mécanique matricielle, due à Heisenberg, Born et Jordan, amena à étudier des matrices comportant un nombre infini de lignes et de colonnes[10]. Par la suite, von Neumann précisa les fondements mathématiques de la mécanique quantique, en remplaçant ces matrices par des opérateurs linéaires sur des espaces de Hilbert.

L'étude théorique des déterminants vient de plusieurs sources[11]. Des problèmes de théorie des nombres amènent Gauss à relier à des matrices (ou plus précisément à leur déterminant) les coefficients d'une forme quadratique ainsi que les applications linéaires en dimension trois. Gotthold Eisenstein développe ces notions, remarquant en particulier qu'en notation moderne, le produit des matrices est non commutatif. Cauchy est le premier à démontrer des résultats généraux sur les déterminants, en utilisant comme définition du déterminant de la matrice A = [a_i,j] le résultat de la substitution, dans le polynôme ${\displaystyle a_{1}a_{2}\cdots a_{n}\prod _{i<j}(a_{j}-a_{i})\;}$, des puissances ak
j par a_jk. Il montre aussi, en 1829, que les valeurs propres d'une matrice symétrique sont réelles[12]. Jacobi étudie les « déterminants fonctionnels » (appelés par la suite jacobiens par Sylvester), utilisés pour décrire des transformations géométriques d'un point de vue infinitésimal ; les livres Vorlesungen über die Theorie der Determinanten[13] de Leopold Kronecker et Zur Determinantentheorie[14] de Karl Weierstrass, tous deux publiés en 1903, définissent pour la première fois les déterminants de manière axiomatique comme formes multilinéaires alternées.

Deux mathématiciens notables au moins ont utilisé le mot dans un sens inhabituel.

Bertrand Russell et Alfred North Whitehead, dans leur Principia Mathematica, utilisent le mot « matrice » dans le cadre de leur axiome de réductibilité (en). Cet axiome permet de réduire le type d'une fonction, les fonctions de type 0 étant identiques à leur extension (en) ; ils appellent « matrice » une fonction n'ayant que des variables libres[15]. Ainsi, par exemple, une fonction Φ(x, y) de deux variables x et y peut être réduite à une collection de fonctions d'une seule variable, par exemple y, en « considérant » la fonction pour toutes les valeurs a_i substituées à la variable x puis réduite à une « matrice » de valeurs en procédant de même pour y : ∀ b_j, ∀ a_i, Φ(a_i, b_j).

Alfred Tarski, dans son Introduction à la logique de 1946, utilise le mot « matrice » comme un synonyme de table de vérité[16].

On définit sur M_m,n(K) une loi de composition interne provenant de l'addition des scalaires :

${\displaystyle (a_{i,j})+(b_{i,j})=(c_{i,j}),\ {\textrm {avec}}\ c_{i,j}=a_{i,j}+b_{i,j}}$.

On ne peut additionner que deux matrices de même taille.

• Exemple :

${\displaystyle {\begin{array}{l}{\begin{pmatrix}0&1&2&3\\4&5&6&7\\\end{pmatrix}}+{\begin{pmatrix}0&0&1&1\\0&1&0&1\\\end{pmatrix}}={\begin{pmatrix}0&1&3&4\\4&6&6&8\\\end{pmatrix}}\end{array}}}$

Pour chaque valeur du couple (m, n), l'espace M_m,n(K) devient alors un groupe abélien, d'élément neutre la matrice nulle, celle dont tous les coefficients valent 0.

On définit aussi une opération à droite de K sur chaque espace M_m,n(K) en associant, à chaque scalaire λ dans K et à chaque matrice (a_i,j) à coefficients dans K, la matrice (a_i,j)λ = (a_i,jλ) obtenue en effectuant la multiplication à droite, dans K, de tous les coefficients de la matrice initiale par λ : c'est la multiplication par un scalaire. Lorsque l'anneau est commutatif, la multiplication peut également s'effectuer à gauche.

En reprenant toujours la matrice A du premier exemple (voir supra) :

${\displaystyle 2A={\begin{pmatrix}0&2&4&6\\8&10&12&14\\16&18&20&22\\\end{pmatrix}}}$

Les espaces M_m,n(K) ainsi obtenus ont donc une structure de K-module à droite, et plus particulièrement de K-espace vectoriel, si K est un corps commutatif.

Le K-module M_m,n(K) est libre de rang mn, c'est-à-dire qu'il possède une base de mn éléments : il suffit de considérer la base canonique (E_i,j)_{1 ≤ i ≤ m, 1 ≤ j ≤ n}. La matrice E_i,j est celle dont tous les coefficients sont nuls sauf celui d'indice (i,j), qui vaut 1.

Les coordonnées dans la base canonique d'une matrice A sont ses coefficients :

${\displaystyle A=\sum _{1\leqslant i\leqslant m \atop 1\leqslant j\leqslant n}a_{i,j}E_{i,j}}$

• Exemple :

${\displaystyle {\begin{pmatrix}0&1&2\\4&3&1\\\end{pmatrix}}=0\cdot {\begin{pmatrix}1&0&0\\0&0&0\\\end{pmatrix}}+1\cdot {\begin{pmatrix}0&1&0\\0&0&0\\\end{pmatrix}}+2\cdot {\begin{pmatrix}0&0&1\\0&0&0\\\end{pmatrix}}+4\cdot {\begin{pmatrix}0&0&0\\1&0&0\\\end{pmatrix}}+3\cdot {\begin{pmatrix}0&0&0\\0&1&0\\\end{pmatrix}}+1\cdot {\begin{pmatrix}0&0&0\\0&0&1\\\end{pmatrix}}}$

On commence par définir le produit d'une matrice ligne par une matrice colonne[18] - [19]. Soit n un nombre entier, L une matrice ligne, x_i ses coefficients, C une matrice colonne, y_i ses coefficients. On les suppose toutes deux de taille n. On définit alors le produit, considéré comme un scalaire ou une matrice de dimension (1, 1) :

${\displaystyle LC={\begin{pmatrix}x_{1}&\dots &x_{n}\end{pmatrix}}{\begin{pmatrix}y_{1}\\\vdots \\y_{n}\end{pmatrix}}=\sum _{1\leq i\leq n}x_{i}y_{i}}$

On remarque la condition de compatibilité sur les tailles des matrices (égalité du nombre de colonnes de la première avec le nombre de lignes de la deuxième). On définit maintenant plus généralement un produit entre deux matrices, la première, (x_i,j) dans M_m,n(K), la deuxième, (y_i,j) dans M_n,p(K), toujours avec une condition de compatibilité sur les tailles (et l'ordre des facteurs de la multiplication ne peut en général pas être changé). Le résultat obtenu est une matrice de M_m,p(K), dont les coefficients (z_i,j) sont obtenus par :

${\displaystyle z_{i,j}=\sum _{1\leq k\leq n}x_{i,k}y_{k,j}=x_{i,1}y_{1,j}+x_{i,2}y_{2,j}+\cdots +x_{i,n}y_{n,j}}$

À la lumière de l'exemple de la multiplication d'une matrice ligne par une matrice colonne, on peut reformuler cette définition en disant que ce coefficient est égal au produit de la ligne i de la première matrice par la colonne j de la deuxième, ce qui s'écrit de la manière suivante, si les L_i sont les lignes de la première matrice, et les C_j les colonnes de la deuxième, le produit est : ${\displaystyle {\begin{pmatrix}L_{1}\\\vdots \\L_{m}\end{pmatrix}}{\begin{pmatrix}C_{1}&\dots &C_{p}\end{pmatrix}}={\begin{pmatrix}L_{1}C_{1}&L_{1}C_{2}&\dots &L_{1}C_{p}\\L_{2}C_{1}&L_{2}C_{2}&\dots &L_{2}C_{p}\\\vdots &\vdots &\ddots &\vdots \\L_{m}C_{1}&L_{m}C_{2}&\dots &L_{m}C_{p}\\\end{pmatrix}}}$.

Le produit matriciel est associatif, distributif à droite et à gauche par rapport à l'addition matricielle. En revanche, même lorsque les dimensions permettent de donner un sens à la question et même si l'anneau des scalaires est commutatif, un produit de matrices ne commute en général pas : AB n'est pas en général égal à BA, par exemple :

${\displaystyle A={\begin{pmatrix}0&0\\1&0\end{pmatrix}}\quad B={\begin{pmatrix}1&0\\0&0\\\end{pmatrix}}\quad AB={\begin{pmatrix}0&0\\1&0\\\end{pmatrix}}\quad BA={\begin{pmatrix}0&0\\0&0\\\end{pmatrix}}.}$

Remarque : le produit de deux matrices non nulles peut être nul, comme ci-dessus ; on dit que l'anneau des matrices n'est pas intègre.

Il arrive même, selon les tailles respectives des matrices A et B, que l'un des deux produits existe et pas l'autre.

La transposition et le produit matriciel sont compatibles au sens suivant :

${\displaystyle \forall A\in M_{n,p}(K),\ \forall B\in M_{p,q}(K),\ (AB)^{\mathsf {T}}=B^{\mathsf {T}}\cdot A^{\mathsf {T}}}$

(même si l'anneau K n'est pas commutatif, en se rappelant que les matrices transposées ont leurs coefficients dans l'anneau opposé K^op).

Pour chaque nombre entier n, on note I_n la matrice carrée de taille n dont les coefficients diagonaux sont égaux à 1 et dont les autres coefficients sont nuls ; elle est appelée matrice identité de taille n.

${\displaystyle I_{1}={\begin{pmatrix}1\end{pmatrix}}\quad I_{2}={\begin{pmatrix}1&0\\0&1\end{pmatrix}}\quad I_{3}={\begin{pmatrix}1&0&0\\0&1&0\\0&0&1\end{pmatrix}}\quad I_{n}=(\delta _{i,j})_{1\leq i\leq n,1\leq j\leq n}}$

où δ_i,j désigne le symbole de Kronecker.

Sous réserve de compatibilité des tailles, les matrices I_n sont neutres à droite et à gauche pour la multiplication.

${\displaystyle \forall A\in M_{m,n}(K),\ I_{m}\,A=A\,I_{n}=A.}$

Soit A une matrice de dimension (m, n). On dit que A est inversible à droite (respectivement à gauche) s'il existe une matrice B de taille (n, m) telle que AB = I_m (respectivement BA = I_n). Elle est simplement dite inversible si elle l'est à la fois à droite et à gauche. Le sous-ensemble de M_n(K) constitué des matrices inversibles possède une structure de groupe pour le produit matriciel ; il est appelé groupe linéaire et noté GL_n(K).

Pour une matrice carrée à coefficients dans un anneau commutatif K, être inversible à droite ou à gauche ou avoir un déterminant inversible dans K (c'est-à-dire non nul si K est un corps) sont trois propriétés équivalentes.

Lorsque l'anneau K est commutatif, l'ensemble M_n(K) des matrices carrées de taille n est donc muni d'une structure de K-algèbre associative et unitaire avec l'addition matricielle, le produit par un scalaire et le produit matriciel.

On appelle matrice scalaire une matrice de la forme I_n λ où λ est un élément de l'anneau K.

${\displaystyle I_{n}a={\begin{pmatrix}a&0&\dots &0\\0&a&\dots &0\\\vdots &\vdots &\ddots &\vdots \\0&0&\dots &a\end{pmatrix}}}$

Ces matrices s'appellent matrices scalaires car elles se comportent comme des scalaires, vis-à-vis de la multiplication :

${\displaystyle \forall \lambda \in K,\ \forall A\in M_{n}(K),\ A(I_{n}\lambda )=A\lambda ~.}$

Lorsque K est commutatif, ou à défaut, lorsque λ est central dans K, c'est-à-dire lorsque λ commute avec tous les éléments de K, on a en outre :

${\displaystyle \forall \lambda \in K,\ \forall A\in M_{n}(K),(\lambda I_{n})\ A=A\lambda ~.}$

Réciproquement, toute matrice B de M_n(K) telle que ∀ A ∈ M_n(K), AB = BA est une matrice scalaire I_n λ où λ est central dans K (ceci se démontre en prenant pour A les matrices de la base canonique).

Une matrice de la forme :

${\displaystyle {\begin{pmatrix}a_{1}&0&\dots &0\\0&a_{2}&\dots &0\\\vdots &\vdots &\ddots &\vdots \\0&0&\dots &a_{n}\end{pmatrix}}}$

sera dite matrice diagonale.

Outre le déterminant, une autre fonction à noter est la trace. Toutes deux apparaissent dans un objet plus général, le polynôme caractéristique, qui à son tour permet d'obtenir certaines caractérisations des matrices diagonalisables (c'est-à-dire semblables à une matrice diagonale), ou de la trigonalisation.

Il existe plusieurs manières de faire agir le groupe linéaire GL_n(K) sur les espaces de matrices, notamment :

• action par multiplication à gauche de GL_m(K) sur M_m,n(K), qui à P et A, associe PA,
• action (à gauche) par multiplication à droite de GL_n(K) sur M_m,n(K), qui à Q ∈ GL_n(K) et A ∈ M_m,n(K), associe AQ⁻¹,
• action par conjugaison de GL_n(K) sur M_n(K), qui à P ∈ GL_n(K) et A ∈ M_n(K), associe PAP⁻¹.

On décrit maintenant les résultats classiques sur ces actions, lorsque les scalaires forment un corps commutatif. Les deux premières actions sont souvent considérées simultanément ; on s'intéresse donc à la question : deux matrices A et B de dimension (m,n) étant données, existe-t-il des matrices P ∈ GL_m(K) et Q ∈ GL_m(K) telles que A = PBQ⁻¹ ? Si tel est le cas, les deux matrices A et B sont dites équivalentes. Le résultat principal est que deux matrices sont équivalentes si et seulement si elles ont même rang, ce qui s'exprime encore en disant que le rang est un invariant complet pour les doubles classes définies par les deux actions de multiplication à gauche et à droite. Par ailleurs, une matrice étant donnée, on peut trouver d'autres matrices privilégiées (les matrices échelonnées) dans la même orbite pour une de ces actions par la méthode du pivot de Gauss.

Pour l'action par conjugaison, deux matrices carrées A et B de taille n dans la même orbite admettent une relation de la forme A = PBP⁻¹, pour une certaine matrice P inversible de taille n ; deux telles matrices sont dites semblables. 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.

Le premier point est de remarquer que le K-module Kⁿ s'identifie canoniquement à l'espace de matrices colonne M_n,1(K) : si e_i est le n-uplet de Kⁿ dont tous les coefficients sont nuls, sauf le i-ème qui vaut 1, on lui associe la i-ème matrice colonne E_i,1 de la base canonique de M_n,1(K) (celle dont tous les coefficients sont nuls sauf le i-ème qui vaut 1), et on étend l'identification par linéarité ; la matrice associée à chaque n-uplet sera appelée matrice coordonnée canonique.

D'autres identifications sont cependant possibles ; lorsqu'on peut parler de base (si l'anneau des scalaires est un corps, par exemple), on peut associer les matrices colonnes élémentaires à n'importe quelle base de l'espace Kⁿ (ou plus généralement d'un K-module libre), puis à nouveau étendre par linéarité ; les matrices associées seront appelées matrices coordonnées dans la base envisagée.

On peut juxtaposer les matrices coordonnées, dans une base fixée, de plusieurs n-uplets. On obtient ainsi la matrice coordonnée d'une famille de vecteurs. Le rang de la matrice est alors défini comme la dimension de la famille de ces vecteurs. En particulier la matrice d'une base dans une autre base est appelée matrice de passage entre ces deux bases, ou matrice de changement de base. Si X et X' sont les matrices coordonnées du même vecteur dans deux bases B et C, et que P est la matrice de passage de la base C dans la base B, on a la relation (une matrice de passage est toujours inversible) :

${\displaystyle X=PX'\quad X'=P^{-1}X}$

Soient E et F deux espaces vectoriels de dimensions respectives n et m sur un corps commutatif[20] K, B une base de E, C une base de F et φ une application linéaire de E dans F.

On appelle matrice de φ dans le couple de bases (B, C) la matrice mat_B,C(φ) de M_m,n(K) telle que pour tout vecteur x de E, si l'on note y = φ(x), X = mat_B(x) et Y = mat_C(y), alors :

${\displaystyle Y={\textrm {mat}}_{B,C}(\varphi )\times X.}$

Si ψ est une deuxième application linéaire, de F dans un troisième espace vectoriel G de base D, alors, relativement aux bases B, C, D, la matrice de la composée ψ ∘ φ est égale au produit des matrices de ψ et φ. Plus précisément :

${\displaystyle {\textrm {mat}}_{B,D}(\psi \circ \varphi )={\textrm {mat}}_{C,D}(\psi )\times {\textrm {mat}}_{B,C}(\varphi ).}$

L'application de L(E, F) dans M_m,n(K) qui à chaque φ associe sa matrice dans (B, C) est un isomorphisme d'espaces vectoriels.

Pour toute matrice M de M_m,n(K), l'application X ↦ MX du K-espace vectoriel M_n,1(K) dans le K-espace vectoriel M_m,1(K) est linéaire. C'est un point clef du lien entre algèbre linéaire et matrices. En conséquence, il arrive souvent que l'on identifie la matrice M avec cette application linéaire. On parlera alors de noyau de la matrice, d'espaces propres de la matrice, d'image de la matrice, etc.

Si B et B' sont deux bases de E, C et C' deux bases de F, P = mat_B(B') la matrice de passage de B vers B' et Q la matrice de passage de C vers C', alors les deux matrices M et M' d'une même application linéaire de E dans F, dans les couples de bases (B, C) et (B', C'), sont liées par : M' = Q⁻¹MP. On constate ainsi que deux matrices équivalentes sont deux matrices qui représentent la même application linéaire dans des bases différentes. En particulier, dans le cas d'un endomorphisme, si l'on impose B = C et B' = C', la formule précédente devient : M' = P⁻¹MP et deux matrices semblables sont deux matrices qui représentent le même endomorphisme dans des bases différentes.

Soient de nouveau E et F deux K-espaces vectoriels de dimensions finies[20], de bases respectives B et C, et φ une application linéaire de E dans F. L'application linéaire transposée ^tφ : F* → E* entre leurs duals est définie par

${\displaystyle \forall y^{*}\in F^{*},\forall x\in E,\quad \left\langle ^{\operatorname {t} }\!\varphi \left(y^{*}\right),x\right\rangle =\left\langle y^{*},\varphi \left(x\right)\right\rangle .}$

Sa matrice dans le couple de bases duales (C*, B*) est liée à celle de φ dans (B, C) par :

${\displaystyle {\textrm {mat}}_{C^{*},B^{*}}(^{\operatorname {t} }\!\varphi )={}^{\operatorname {t} }\!({\textrm {mat}}_{B,C}(\varphi )).}$

Lorsque l'anneau n'est pas commutatif, si l'on représente les vecteurs par des matrices colonne, l'algèbre linéaire n'est compatible avec le calcul matriciel que si les modules ou espaces vectoriels considérés sont à droite, comme dans les articles détaillés signalés ci-dessus, une application linéaire correspondant à la multiplication à gauche d'un vecteur colonne par la matrice qui la représente. Si l'on tient à avoir des modules ou espaces vectoriels à gauche, il faut représenter les vecteurs par des matrices ligne, une application linéaire étant cette fois représentée par la multiplication à droite d'un vecteur ligne par la matrice qui la représente[21].

En général, un système de m équations linéaires à n inconnues peut être écrit sous la forme suivante :

${\displaystyle \left\{{\begin{matrix}a_{1,1}x_{1}+a_{1,2}x_{2}+...+a_{1,n}x_{n}=b_{1}\\a_{2,1}x_{1}+a_{2,2}x_{2}+...+a_{2,n}x_{n}=b_{2}\\\vdots \\\vdots \\a_{m,1}x_{1}+a_{m,2}x_{2}+...+a_{m,n}x_{n}=b_{m}\end{matrix}}\right.}$

où x₁, ..., x_n sont les inconnues et les nombres a_i,j sont les coefficients du système.

Ce système peut s'écrire sous la forme matricielle :

${\displaystyle Ax=b}$

avec :

${\displaystyle A={\begin{pmatrix}a_{1,1}&a_{1,2}&\cdots &a_{1,n}\\a_{2,1}&a_{2,2}&\cdots &a_{2,n}\\\vdots &\vdots &\ddots &\vdots \\a_{m,1}&a_{m,2}&\cdots &a_{m,n}\end{pmatrix}};\qquad x={\begin{pmatrix}x_{1}\\x_{2}\\\vdots \\x_{n}\end{pmatrix}}\quad {\text{et}}\quad b={\begin{pmatrix}b_{1}\\b_{2}\\\vdots \\b_{m}\end{pmatrix}};}$

la théorie de la résolution des systèmes utilise les invariants liés à la matrice A (appelée matrice du système), par exemple son rang, et, dans le cas où A est inversible, son déterminant (voir l'article règle de Cramer).

Soient E un K-module libre et B = (e₁, … , e_n) une base de E.

Soit ${\displaystyle f:E\times E\to K}$ une forme bilinéaire. On définit la matrice de ${\displaystyle f}$ dans la base B par la formule suivante :

${\displaystyle {\textrm {mat}}_{B}\,f=(f(e_{i},e_{j}))_{1\leq i\leq n,\ 1\leq j\leq n}={\begin{pmatrix}f(e_{1},e_{1})&f(e_{1},e_{2})&\dots &f(e_{1},e_{n})\\f(e_{2},e_{1})&f(e_{2},e_{2})&\dots &f(e_{2},e_{n})\\\vdots &\vdots &\ddots &\vdots \\f(e_{n},e_{1})&f(e_{n},e_{2})&\dots &f(e_{n},e_{n})\\\end{pmatrix}}}$

Dans le cas particulier où K = ℝ et f est un produit scalaire, cette matrice est appelée matrice de Gram.

La matrice mat_B f est symétrique (respectivement antisymétrique) si et seulement si la forme bilinéaire ${\displaystyle f}$ 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 = mat_B f. On a alors la formule : ${\displaystyle f(x,y)={}^{\operatorname {t} }\!X~A~Y}$.

Deux formes bilinéaires sont égales si et seulement si elles ont la même matrice dans une base donnée.

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.

Soient E un K-module libre et B, C deux bases de E. Soit ${\displaystyle f:E\times E\to K}$ une forme bilinéaire.

Notons M = mat_B f la matrice de f dans la base B et M' = mat_C f la matrice de f dans la base C. Notons P = mat_{B 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) :

${\displaystyle M'={}^{\operatorname {t} }\!P~M~P}$

Deux matrices carrées A et B sont dites congruentes s'il existe une matrice inversible P telle que

${\displaystyle A={}^{\operatorname {t} }\!P~B~P~.}$

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.

• 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⁻¹DP.
• 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.

• 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.

• 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.

• 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.

Dans tout ce paragraphe, K = ℝ ou ℂ.

Une norme matricielle est une norme d'algèbre sur l'algèbre M_n(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 M_n(K), toute norme N subordonnée à une norme sur Kⁿ est une norme d'algèbre vérifiant de plus N(I_n) = 1 (la réciproque est fausse).

L'espace vectoriel M_m,n(ℝ), canoniquement isomorphe à ℝ^mn, hérite de sa structure euclidienne. Le produit scalaire se transcrit en

où ${\displaystyle \operatorname {tr} }$ désigne la trace (i.e., ${\displaystyle \textstyle \operatorname {tr} \,M=\sum _{i=1}^{n}m_{ii}}$) et les a_i,j (resp. b_i,j) désignent les éléments de A (resp. B). La norme associée à ce produit scalaire est la norme de Frobenius ou norme de Hilbert-Schmidt :

${\displaystyle \|A\|_{F}=\left(\sum _{1\leq i\leq m \atop 1\leq j\leq n}a_{i,j}^{2}\right)^{\frac {1}{2}}=\|\sigma (A)\|_{2},}$

où σ(A) est le vecteur des valeurs singulières de A et ${\displaystyle \|\cdot \|_{2}}$ est la norme euclidienne.

Si m = n > 1, il ne s'agit pas d'une norme subordonnée, puisque ${\displaystyle \|I_{n}\|_{F}={\sqrt {n}}\neq 1.}$

L'inégalité de Cauchy-Schwarz s'écrit (comme pour tout produit scalaire) :

${\displaystyle \forall A,B\in \mathbb {R} ^{m\times n},\langle A,B\rangle \leqslant \|A\|_{F}\,\|B\|_{F}.}$

Cette inégalité peut être renforcée par l'inégalité de trace de von Neumann[22] :

${\displaystyle \forall A,B\in \mathbb {R} ^{m\times n},\langle A,B\rangle \leqslant \sigma (A)^{\mathsf {T}}\sigma (B),}$

où σ(A) est le vecteur des valeurs singulières de A, rangées en ordre décroissant. Elle a la même structure que l'inégalité de Ky Fan, laquelle suppose que les matrices sont carrées et symétriques (on peut alors remplacer σ(•) par le vecteur des valeurs propres).

L'espace vectoriel M_m,n(ℂ) est muni d'une structure similaire d'espace hermitien.

• J.-M. Arnaudiès et H. Fraysse, Cours de mathématiques, Dunod, 1980
• Rached Mneimné, Réduction des endomorphismes, Calvage et Mounet, Paris, 2006 (ISBN 978-2-916352-01-5)
• P. Wira, Un rappel sur les matrices, support de cours, Université de Haute Alsace, Mulhouse, 2000

• Frédéric Brechenmacher, Les matrices : formes de représentation et pratiques opératoires (1850-1930)