Movatterモバイル変換


[0]ホーム

URL:


Jump to content
WikipediaThe Free Encyclopedia
Search

Hankel matrix

From Wikipedia, the free encyclopedia
Square matrix in which each ascending skew-diagonal from left to right is constant

Inlinear algebra, aHankel matrix (orcatalecticant matrix), named afterHermann Hankel, is a rectangular matrix in which each ascending skew-diagonal from left to right is constant. For example,

[abcdebcdefcdefgdefghefghi].{\displaystyle \qquad {\begin{bmatrix}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{bmatrix}}.}

More generally, aHankel matrix is anyn×n{\displaystyle n\times n}matrixA{\displaystyle A} of the form

A=[a0a1a2an1a1a2a2a2n4a2n4a2n3an1a2n4a2n3a2n2].{\displaystyle A={\begin{bmatrix}a_{0}&a_{1}&a_{2}&\ldots &a_{n-1}\\a_{1}&a_{2}&&&\vdots \\a_{2}&&&&a_{2n-4}\\\vdots &&&a_{2n-4}&a_{2n-3}\\a_{n-1}&\ldots &a_{2n-4}&a_{2n-3}&a_{2n-2}\end{bmatrix}}.}

In terms of the components, if thei,j{\displaystyle i,j} element ofA{\displaystyle A} is denoted withAij{\displaystyle A_{ij}}, and assumingij{\displaystyle i\leq j}, then we haveAi,j=Ai+k,jk{\displaystyle A_{i,j}=A_{i+k,j-k}} for allk=0,...,ji.{\displaystyle k=0,...,j-i.}

Properties

[edit]

Hankel operator

[edit]

Given aformal Laurent seriesf(z)=n=Nanzn,{\displaystyle f(z)=\sum _{n=-\infty }^{N}a_{n}z^{n},}the correspondingHankel operator is defined as[2]Hf:C[z]z1C[[z1]].{\displaystyle H_{f}:\mathbf {C} [z]\to \mathbf {z} ^{-1}\mathbf {C} [[z^{-1}]].}This takes apolynomialgC[z]{\displaystyle g\in \mathbf {C} [z]} and sends it to the productfg{\displaystyle fg}, but discards all powers ofz{\displaystyle z} with a non-negative exponent, so as to give an element inz1C[[z1]]{\displaystyle z^{-1}\mathbf {C} [[z^{-1}]]}, theformal power series with strictly negative exponents. The mapHf{\displaystyle H_{f}} is in a natural wayC[z]{\displaystyle \mathbf {C} [z]}-linear, and its matrix with respect to the elements1,z,z2,C[z]{\displaystyle 1,z,z^{2},\dots \in \mathbf {C} [z]} andz1,z2,z1C[[z1]]{\displaystyle z^{-1},z^{-2},\dots \in z^{-1}\mathbf {C} [[z^{-1}]]} is the Hankel matrix[a1a2a2a3a3a4].{\displaystyle {\begin{bmatrix}a_{1}&a_{2}&\ldots \\a_{2}&a_{3}&\ldots \\a_{3}&a_{4}&\ldots \\\vdots &\vdots &\ddots \end{bmatrix}}.}Any Hankel matrix arises in this way. Atheorem due toKronecker says that therank of this matrix is finite precisely iff{\displaystyle f} is arational function, that is, a fraction of two polynomialsf(z)=p(z)q(z).{\displaystyle f(z)={\frac {p(z)}{q(z)}}.}

Approximations

[edit]

We are often interested in approximations of the Hankel operators, possibly by low-order operators. In order to approximate the output of the operator, we can use the spectral norm (operator 2-norm) to measure the error of our approximation. This suggestssingular value decomposition as a possible technique to approximate the action of the operator.

Note that the matrixA{\displaystyle A} does not have to be finite. If it is infinite, traditional methods of computing individual singular vectors will not work directly. We also require that the approximation is a Hankel matrix, which can be shown withAAK theory.

Hankel matrix transform

[edit]
Not to be confused withHankel transform.

TheHankel matrix transform, or simplyHankel transform, of asequencebk{\displaystyle b_{k}} is the sequence of the determinants of the Hankel matrices formed frombk{\displaystyle b_{k}}. Given an integern>0{\displaystyle n>0}, define the corresponding(n×n){\displaystyle (n\times n)}-dimensional Hankel matrixBn{\displaystyle B_{n}} as having the matrix elements[Bn]i,j=bi+j.{\displaystyle [B_{n}]_{i,j}=b_{i+j}.} Then the sequencehn{\displaystyle h_{n}} given byhn=detBn{\displaystyle h_{n}=\det B_{n}}is the Hankel transform of the sequencebk.{\displaystyle b_{k}.} The Hankel transform is invariant under thebinomial transform of a sequence. That is, if one writescn=k=0n(nk)bk{\displaystyle c_{n}=\sum _{k=0}^{n}{n \choose k}b_{k}}as the binomial transform of the sequencebn{\displaystyle b_{n}}, then one hasdetBn=detCn.{\displaystyle \det B_{n}=\det C_{n}.}

Applications of Hankel matrices

[edit]

Hankel matrices are formed when, given a sequence of output data, a realization of an underlying state-space orhidden Markov model is desired.[3] The singular value decomposition of the Hankel matrix provides a means of computing theA,B, andC matrices which define the state-space realization.[4] The Hankel matrix formed from the signal has been found useful for decomposition of non-stationary signals and time-frequency representation.

Method of moments for polynomial distributions

[edit]

Themethod of moments applied to polynomial distributions results in a Hankel matrix that needs to beinverted in order to obtain the weight parameters of the polynomial distribution approximation.[5]

Positive Hankel matrices and the Hamburger moment problems

[edit]
Further information:Hamburger moment problem

See also

[edit]

Notes

[edit]
  1. ^Yasuda, M. (2003). "A Spectral Characterization of Hermitian Centrosymmetric and Hermitian Skew-Centrosymmetric K-Matrices".SIAM J. Matrix Anal. Appl.25 (3):601–605.doi:10.1137/S0895479802418835.
  2. ^Fuhrmann 2012, §8.3
  3. ^Aoki, Masanao (1983)."Prediction of Time Series".Notes on Economic Time Series Analysis : System Theoretic Perspectives. New York: Springer. pp. 38–47.ISBN 0-387-12696-1.
  4. ^Aoki, Masanao (1983)."Rank determination of Hankel matrices".Notes on Economic Time Series Analysis : System Theoretic Perspectives. New York: Springer. pp. 67–68.ISBN 0-387-12696-1.
  5. ^J. Munkhammar, L. Mattsson, J. Rydén (2017) "Polynomial probability distribution estimation using the method of moments". PLoS ONE 12(4): e0174573.https://doi.org/10.1371/journal.pone.0174573

References

[edit]
Matrix classes
Explicitly constrained entries
Constant
Conditions oneigenvalues or eigenvectors
Satisfying conditions onproducts orinverses
With specific applications
Used instatistics
Used ingraph theory
Used in science and engineering
Related terms
Authority control databases: NationalEdit this at Wikidata
Retrieved from "https://en.wikipedia.org/w/index.php?title=Hankel_matrix&oldid=1277636549"
Categories:
Hidden categories:

[8]ページ先頭

©2009-2025 Movatter.jp