En algèbre linéaire, une matrice circulante[1] est une matrice carrée dans laquelle on passe d'une ligne à la suivante par permutation circulaire (décalage vers la droite) des coefficients.
Une matrice circulante de taille n est donc de la forme
où les coefficients ci sont des complexes.
Une matrice circulante constitue un cas particulier de matrice de Toeplitz, de matrice de Frobenius (c'est la matrice générique de la multiplication par un élément de l'algèbre de groupe ℂ[ℤ/nℤ] et aussi un cas particulier de carré latin).
La réduction des matrices circulantes fait intervenir les formules de la transformation de Fourier discrète. En analyse numérique, les systèmes circulants peuvent être résolus très efficacement par transformée de Fourier rapide.
On parle parfois de matrice anticirculante ou circulante gauche quand on effectue un décalage à gauche des coefficients en passant d'une ligne à la suivante.
Pour alléger les notations, on désigne par C(c0, … , cn–1) la matrice circulante précédente.
on peut constater que toute matrice circulante est un polynôme en J
Réciproquement, comme Jn est la matrice identité, tout polynôme en J est une matrice circulante.
Ainsi la somme, le produit de matrices circulantes sont circulantes, et un tel produit est commutatif. L'ensemble des matrices circulantes n'est autre que l'algèbre commutative des polynômes en J.
La matrice J, vérifiant Jn = I, est diagonalisable sur ℂ avec pour valeurs propres des racines n-ièmes de l'unité.
On appelle donc , racine primitive de l'unité. On vérifie alors sans peine que pour tout k[2]:
est vecteur propre de J associé à la valeur propre ωk.
On a donc exhibé, pour k allant de 0 à n – 1, une famille de n vecteurs propres associés à des valeurs propres distinctes, soit une base propre pour J.
Par conséquent, c'est une base propre aussi pour tout polynôme en J, c'est-à-dire toute matrice circulante. Les valeurs propres de C(c0, … , cn–1) sont donc les
qui, cette fois, ne sont plus nécessairement distinctes.
On peut prendre, pour matrice de passage de la base canonique à la base propre, la matrice
Cette matrice U est unitaire (U*U = I) et les formules de passage précédentes s'écrivent, en notant Λ la matrice diagonale de coefficients les valeurs propres
Une nouvelle définition possible pour l'ensemble des matrices circulantes est l'ensemble des matrices de la forme UDU* avec D diagonale. Géométriquement, cela correspond aux endomorphismes qui admettent la base orthonormale des Xk comme base de vecteurs propres.
Le déterminant circulant est le déterminant de la matrice circulante ; comme pour toute autre matrice, il est égal au produit des valeurs propres[2]:
Toute matrice est inversible si et seulement si son déterminant est non nul, et dans le cas d'une matrice circulante sa matrice inverse est aussi une matrice circulante.
Les formules de changement de base à l'aide de la matrice U ont un intérêt particulier. La formule de passage des coefficients aux valeurs propres est la définition classique d'une transformée de Fourier discrète. On peut retrouver les coefficients à partir des valeurs propres en effectuant cette fois une transformée inverse
Soit le système circulant avec matrice circulante de taille . Ce système peut se réécrire à l'aide d'un produit de convolution discret
en notant la première colonne de la matrice et en périodisant les composantes des vecteurs , et . La transformée de Fourier discrète transforme cette équation de convolution en un produit composante par composante[4].
et ainsi
Cet algorithme de résolution est bien plus rapide que l'élimination de Gauss-Jordan, et l'est d'autant plus si l'on a recours à la transformée de Fourier rapide.