Aller au contenu
Wikipédia l'encyclopédie libre

Matrice de Hankel

Un article de Wikipédia, l'encyclopédie libre.

En algèbre linéaire, une matrice de Hankel, du nom du mathématicien Hermann Hankel, est une matrice carrée dont les valeurs sont constantes le long des diagonales ascendantes, c'est-à-dire dont les indices vérifient la relation a i , j = a i 1 , j + 1 {\displaystyle a_{i,j}=a_{i-1,j+1}} {\displaystyle a_{i,j}=a_{i-1,j+1}}

Par exemple une matrice de Hankel de taille 5 s'écrit sous la forme

( a b c d e b c d e f c d e f g d e f g h e f g h i ) {\displaystyle {\begin{pmatrix}a&b&c&d&e\\b&c&d&e&f\\c&d&e&f&g\\d&e&f&g&h\\e&f&g&h&i\end{pmatrix}}} {\displaystyle {\begin{pmatrix}a&b&c&d&e\\b&c&d&e&f\\c&d&e&f&g\\d&e&f&g&h\\e&f&g&h&i\end{pmatrix}}}.

Les matrices de Toeplitz ont, elles, des valeurs constantes sur les diagonales descendantes.

Sur un espace de Hilbert muni d'une base hilbertienne, on peut définir plus généralement un opérateur de Hankel. Ce dernier admet pour représentation une matrice de Hankel infinie, c'est-à-dire que le coefficient a i , j = ( e i | a ( e j ) ) {\displaystyle a_{i,j}=(e_{i}|a(e_{j}))} {\displaystyle a_{i,j}=(e_{i}|a(e_{j}))}, dépend seulement de i + j {\displaystyle i+j} {\displaystyle i+j}.

Déterminant et transformation de Hankel

[modifier | modifier le code ]

À toute suite ( b n ) n 0 {\displaystyle (b_{n})_{n\geq 0}} {\displaystyle (b_{n})_{n\geq 0}} on peut associer la suite des déterminants ( h n ) n 0 {\displaystyle (h_{n})_{n\geq 0}} {\displaystyle (h_{n})_{n\geq 0}} des matrices de Hankel successives[1]  :

h n = | b 0 b 1 b 2 b n b 1 b 2 b 3 b n + 1 b 2 b 3 b 4 b n + 2 b n b n + 1 b n + 2 b 2 n | {\displaystyle h_{n}={\begin{vmatrix}b_{0}&b_{1}&b_{2}&\dots &b_{n}\\b_{1}&b_{2}&b_{3}&\dots &b_{n+1}\\b_{2}&b_{3}&b_{4}&\dots &b_{n+2}\\\vdots &\vdots &\vdots &&\vdots \\b_{n}&b_{n+1}&b_{n+2}&\dots &b_{2n}\\\end{vmatrix}}} {\displaystyle h_{n}={\begin{vmatrix}b_{0}&b_{1}&b_{2}&\dots &b_{n}\\b_{1}&b_{2}&b_{3}&\dots &b_{n+1}\\b_{2}&b_{3}&b_{4}&\dots &b_{n+2}\\\vdots &\vdots &\vdots &&\vdots \\b_{n}&b_{n+1}&b_{n+2}&\dots &b_{2n}\\\end{vmatrix}}} pour tout n 0 {\displaystyle n\geq 0} {\displaystyle n\geq 0}.

Propriétés

[modifier | modifier le code ]
  • Une suite est nulle si et seulement si sa transformée de Hankel est nulle.
  • Une suite satisfait à une relation de récurrence linéaire à coefficients constants si et seulement si la transformée de Hankel est nulle à partir d'un certain rang.

Éléments propres dans un cas particulier

[modifier | modifier le code ]

Un cas particulier de matrices de Hankel est celui de matrices anticirculantes, lorsque b n + j = b j 1 {\displaystyle b_{n+j}=b_{j-1}} {\displaystyle b_{n+j}=b_{j-1}}.

Dans ce cas, on peut aisément diagonaliser les matrices de Hankel : On considère la base de diagonalisation des matrices circulantes de taille n {\displaystyle n} {\displaystyle n} ; on note P {\displaystyle P} {\displaystyle P} la matrice de passage associée d'éléments génériques w ( i 1 ) ( j 1 ) {\displaystyle w^{(i-1)(j-1)}} {\displaystyle w^{(i-1)(j-1)}} w = exp ( 2 i π n ) {\displaystyle w=\exp \left({\frac {2i\pi }{n}}\right)} {\displaystyle w=\exp \left({\frac {2i\pi }{n}}\right)} est une racine n {\displaystyle n} {\displaystyle n}-ième, primitive, de l'unité. On note C j {\displaystyle C_{j}} {\displaystyle C_{j}} ses colonnes.

Si M {\displaystyle M} {\displaystyle M} est une matrice de Hankel, on note M = P 1 M P {\displaystyle M'=P^{-1}MP} {\displaystyle M'=P^{-1}MP}. On remarque que M C 1 = L C 1 {\displaystyle M'C_{1}=LC_{1}} {\displaystyle M'C_{1}=LC_{1}}, où L = b 0 + b 1 + b 2 + + b n 1 {\displaystyle L=b_{0}+b_{1}+b_{2}+\dots +b_{n-1}} {\displaystyle L=b_{0}+b_{1}+b_{2}+\dots +b_{n-1}} est la somme des coefficients de chaque ligne.

Si n {\displaystyle n} {\displaystyle n} est pair, -1 est racine n-ième de l'unité, le vecteur colonne C n + 2 2 {\displaystyle C_{\frac {n+2}{2}}} {\displaystyle C_{\frac {n+2}{2}}} est ( 1 , 1 , 1 , , 1 ) {\displaystyle (1,-1,1,\ldots ,-1)} {\displaystyle (1,-1,1,\ldots ,-1)} et

M C n + 2 2 = L C n + 2 2 {\displaystyle M'C_{\frac {n+2}{2}}=L'C_{\frac {n+2}{2}}} {\displaystyle M'C_{\frac {n+2}{2}}=L'C_{\frac {n+2}{2}}}, où L = b 0 b 1 + b 2 + b n 1 {\displaystyle L'=b_{0}-b_{1}+b_{2}+\dots -b_{n-1}} {\displaystyle L'=b_{0}-b_{1}+b_{2}+\dots -b_{n-1}} est la somme alternée des coefficients de chaque ligne.

Si j { 1 , n + 2 2 } {\displaystyle j\notin \left\{1,{\frac {n+2}{2}}\right\}} {\displaystyle j\notin \left\{1,{\frac {n+2}{2}}\right\}}, on voit que M C j = L j C n + 2 j {\displaystyle M'C_{j}=L_{j}C_{n+2-j}} {\displaystyle M'C_{j}=L_{j}C_{n+2-j}} L j = b 0 + b 1 w j 1 + b 2 w 2 ( j 1 ) + b n 1 w ( n 1 ) ( j 1 ) {\displaystyle L_{j}=b_{0}+b_{1}w^{j-1}+b_{2}w^{2(j-1)}+\dots -b_{n-1}w^{(n-1)(j-1)}} {\displaystyle L_{j}=b_{0}+b_{1}w^{j-1}+b_{2}w^{2(j-1)}+\dots -b_{n-1}w^{(n-1)(j-1)}} ; la valeur en w j 1 {\displaystyle w^{j-1}} {\displaystyle w^{j-1}} du polynôme associé : P ( X ) = k = 0 n b k X k {\displaystyle P(X)=\sum _{k=0}^{n}b_{k}X^{k}} {\displaystyle P(X)=\sum _{k=0}^{n}b_{k}X^{k}}.

On remarque que le plan V e c t ( C j , C n + 2 j ) {\displaystyle \mathrm {Vect} (C_{j},C_{n+2-j})} {\displaystyle \mathrm {Vect} (C_{j},C_{n+2-j})} est stable sous l'action de l'endomorphisme canoniquement associé à M {\displaystyle M} {\displaystyle M} ; la restriction à cet espace, de l'endomorphisme (complexe) associé, a pour matrice dans la base des deux colonnes ( C j , C n + 2 j ) {\displaystyle (C_{j},C_{n+2-j})} {\displaystyle (C_{j},C_{n+2-j})} : A = ( 0 L j L n + 2 j 0 ) {\displaystyle A={\begin{pmatrix}0&L_{j}\\L_{n+2-j}&0\end{pmatrix}}} {\displaystyle A={\begin{pmatrix}0&L_{j}\\L_{n+2-j}&0\end{pmatrix}}} ; on en déduit M {\displaystyle M'} {\displaystyle M'}, puis une diagonalisation de M {\displaystyle M} {\displaystyle M}.

Notes et références

[modifier | modifier le code ]
  1. (en) John W. Layman, « The Hankel transform and some of its properties », Journal of Integer Sequences , vol. 4,‎ (lire en ligne).

Voir aussi

[modifier | modifier le code ]
v · m
Forme
Transformée
Relation
Propriété
Famille
Associée
Résultats
Décompositions
Articles liés

AltStyle によって変換されたページ (->オリジナル) /