左から右への各上昇斜対角線が一定である正方行列
線型代数学において、ハンケル行列(カタレクティカント行列)は、ヘルマン・ハンケルにちなんで名付けられ、左から右への各上昇する歪対角成分が一定である行列である。例えば、
より一般的には、ハンケル行列とは、以下の形式の
行列である。

成分の観点から見ると、の元が で表され、 と仮定すると、すべての に対して




プロパティ
- 任意のハンケル行列は対称です。
- を交換行列とします。 がハンケル行列である場合、 はテプリッツ行列です。





- が実対称ならば、は符号までと同じ固有値を持つ。 [1]



- ヒルベルト行列はハンケル行列の例です。
- ハンケル行列の行列式はカタレクティカントと呼ばれます。
ハンケル演算子
形式ローラン級数
が与えられた場合、
対応するハンケル演算子は次のように定義されます[2]これは多項式
を取り、それを積 に渡しますが、非負の指数を持つの累乗をすべて破棄して、 の要素、つまり厳密に負の指数を持つ形式冪級数を与えます。写像は自然に -線型であり、その行列は要素とに関してハンケル行列です
。任意の
ハンケル行列はこのようにして生じます。クロネッカーによる定理によれば、この行列の階数はが有理関数、つまり2つの多項式の分数で
ある場合に限り有限となります。
![{\displaystyle g\in \mathbf {C} [z]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3df7d024d0b972a46da42c84aeff389ff5458f1f)


![{\displaystyle z^{-1}\mathbf {C} [[z^{-1}]]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5a964b9b781a86975a982d8911e5bef310129193)

![{\displaystyle \mathbf {C} [z]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4df688b72293d8ac85f15267559929137e35ff1c)
![{\displaystyle 1,z,z^{2},\dots \in \mathbf {C} [z]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2814851a432ce68c94aab843fac3ba4882bbd3a6)
![{\displaystyle z^{-1},z^{-2},\dots \in z^{-1}\mathbf {C} [[z^{-1}]]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bbf789a2c0db73be54af3b1919900b73cf98bad9)


近似値
ハンケル演算子の近似、おそらくは低次の演算子による近似に関心を持つことがよくあります。演算子の出力を近似するには、スペクトル ノルム (演算子 2 ノルム) を使用して近似のエラーを測定することができます。これは、演算子のアクションを近似する可能性のある手法として
特異値分解を示唆しています。
行列は有限である必要はないことに注意してください。行列が無限である場合、個々の特異ベクトルを計算する従来の方法は直接機能しません。また、近似はハンケル行列である必要があり、これは AAK 理論で示すことができます。

数列のハンケル行列変換、または単にハンケル変換は、から形成されるハンケル行列の行列式の数列です。整数 が与えられた場合、対応する 次元ハンケル行列を行列要素を持つものとして定義します。すると、によって与えられる
数列は
数列のハンケル変換です。ハンケル変換は数列の二項変換に対して不変です。つまり、 を
数列の二項変換として書くと
、 




![{\displaystyle [B_{n}]_{i,j}=b_{i+j}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/eec3523edf9c2a5ccedd03a4c69779106d89269b)





ハンケル行列の応用
ハンケル行列は、出力データのシーケンスが与えられたときに、基礎となる状態空間モデルまたは隠れマルコフモデルの実現が望まれる場合に形成されます。[3]ハンケル行列の特異値分解は、状態空間実現を定義するA、B、およびC行列 を計算する手段を提供します。 [4]信号から形成されたハンケル行列は、非定常信号の分解と時間周波数表現に有用であることがわかっています。
多項式分布のモーメント法
多項式分布にモーメント法を適用するとハンケル行列が得られ、これを逆行列にすることで多項式分布近似の重みパラメータを得ることができる。[5]
正のハンケル行列とハンバーガーモーメント問題
参照
注記
- ^ Yasuda, M. (2003). 「エルミート中心対称およびエルミート歪中心対称 K 行列のスペクトル特性」SIAM J. Matrix Anal. Appl . 25 (3): 601–605. doi :10.1137/S0895479802418835.
- ^ フーアマン 2012、§8.3
- ^ 青木正直(1983) 「時系列の予測」経済時系列分析に関するノート:システム理論的視点ニューヨーク: シュプリンガー pp. 38–47. ISBN 0-387-12696-1。
- ^ 青木正尚 (1983) 「ハンケル行列のランク決定」経済時系列分析に関するノート:システム理論的視点ニューヨーク: シュプリンガー pp. 67–68. ISBN 0-387-12696-1。
- ^ J. Munkhammar、L. Mattsson、J. Rydén (2017)「モーメント法を用いた多項式確率分布推定」PLoS ONE 12(4): e0174573. https://doi.org/10.1371/journal.pone.0174573
参考文献
- Brent RP (1999)、「構造化線形システムのための高速アルゴリズムの安定性」、構造付き行列のための高速で信頼性の高いアルゴリズム(編集者—T. Kailath、AH Sayed)、第 4 章 ( SIAM )。
- Fuhrmann, Paul A. (2012).線形代数への多項式アプローチ。Universitext (第2版)。ニューヨーク、NY: Springer。doi : 10.1007 /978-1-4614-0338-8。ISBN 978-1-4614-0337-1.ZBL1239.15001 。