カーレマン行列

数学におけるカーレマン行列(カーレマンぎょうれつ、: Carleman matrix)は、函数の合成行列の積に読み替えるために用いられる行列である。反復理論において、パターン認識のみでは反復の出来ない連続な反復函数を見つけるために用いられる。その他、確率母函数マルコフ連鎖の理論で用いられる。

定義

函数 f ( x ) {\displaystyle f(x)} カーレマン行列は次のように定義される。

M [ f ] j k = 1 k ! [ d k d x k ( f ( x ) ) j ] x = 0 , {\displaystyle M[f]_{jk}={\frac {1}{k!}}\left[{\frac {d^{k}}{dx^{k}}}(f(x))^{j}\right]_{x=0},}

したがって、次の(テイラー級数の)方程式が満たされる。

( f ( x ) ) j = k = 0 M [ f ] j k x k . {\displaystyle (f(x))^{j}=\sum _{k=0}^{\infty }M[f]_{jk}x^{k}.}

例えば、 f ( x ) {\displaystyle f(x)} を計算すると

f ( x ) = k = 0 M [ f ] 1 , k x k {\displaystyle f(x)=\sum _{k=0}^{\infty }M[f]_{1,k}x^{k}}

のように、 M [ f ] {\displaystyle M[f]} の第 1 行と列ベクトル [ 1 , x , x 2 , x 3 , ] {\displaystyle [1,x,x^{2},x^{3},\ldots ]^{\top }} (⊤ は転置)のドット積で与えられる。

M [ f ] {\displaystyle M[f]} の次の行の成分は、以下のような f ( x ) {\displaystyle f(x)} の 2 次のベキを与える。

f ( x ) 2 = k = 0 M [ f ] 2 , k x k . {\displaystyle f(x)^{2}=\sum _{k=0}^{\infty }M[f]_{2,k}x^{k}.}

そして、 f ( x ) {\displaystyle f(x)} のゼロ次のベキを M [ f ] {\displaystyle M[f]} に含めるように、行 0 は第一成分を除いてすべてゼロであるようにする。すなわち

f ( x ) 0 = 1 = k = 0 M [ f ] 0 , k x k = 1 + k = 1 0 x k {\displaystyle f(x)^{0}=1=\sum _{k=0}^{\infty }M[f]_{0,k}x^{k}=1+\sum _{k=1}^{\infty }0\cdot x^{k}}

が成り立つ。

したがって、 M [ f ] {\displaystyle M[f]} と列ベクトル [ 1 , x , x 2 , ] {\displaystyle [1,x,x^{2},\ldots ]^{\top }} のドット積は、列ベクトル [ 1 , f ( x ) , f ( x ) 2 , ] {\displaystyle [1,f(x),f(x)^{2},\ldots ]^{\top }} を導く。すなわち、

M [ f ] [ 1 , x , x 2 , x 3 , ] = [ 1 , f ( x ) , ( f ( x ) ) 2 , ( f ( x ) ) 3 , ] . {\displaystyle M[f]\cdot [1,x,x^{2},x^{3},\ldots ]^{\top }=[1,f(x),(f(x))^{2},(f(x))^{3},\ldots ]^{\top }.}

ベル行列

函数 f ( x ) {\displaystyle f(x)} ベル行列(Bell matrix)は次のように定義される。

B [ f ] j k = 1 j ! [ d j d x j ( f ( x ) ) k ] x = 0 . {\displaystyle B[f]_{jk}={\frac {1}{j!}}\left[{\frac {d^{j}}{dx^{j}}}(f(x))^{k}\right]_{x=0}.}

したがって次の方程式が満たされる。

( f ( x ) ) k = j = 0 B [ f ] j k x j . {\displaystyle (f(x))^{k}=\sum _{j=0}^{\infty }B[f]_{jk}x^{j}.}

これは上述のカーレマン行列の転置である。

ジャボチンスキー行列

エリ・ジャボチンスキー(Eri Jabotinsky)は 1947 年、多項式の畳み込みを表現する目的で行列の概念を開発した。このため何人かの研究者はベル行列のことをジャボチンスキー行列と呼んでおり、今後この名がより正式なものになる可能性もある。

一般化

函数のカーレマン行列の一般化は、任意の点の周りで次のように定義される。

M [ f ] x 0 = M x [ x x 0 ] M [ f ] M x [ x + x 0 ] {\displaystyle M[f]_{x_{0}}=M_{x}[x-x_{0}]M[f]M_{x}[x+x_{0}]}

あるいは M [ f ] x 0 = M [ g ] {\displaystyle M[f]_{x_{0}}=M[g]} where g ( x ) = f ( x + x 0 ) x 0 {\displaystyle g(x)=f(x+x_{0})-x_{0}} 。この定義より、行列のベキを次のように表現できる。

( M [ f ] x 0 ) n = M x [ x x 0 ] M [ f ] n M x [ x + x 0 ] {\displaystyle (M[f]_{x_{0}})^{n}=M_{x}[x-x_{0}]M[f]^{n}M_{x}[x+x_{0}]}

行列の性質

これまでに紹介した行列は、次の基本的な関係式を満たす。

  • M [ f g ] = M [ f ] M [ g ] , {\displaystyle M[f\circ g]=M[f]M[g],}
  • B [ f g ] = B [ g ] B [ f ] . {\displaystyle B[f\circ g]=B[g]B[f].}

このことよりカーレマン行列 M f ( x ) {\displaystyle f(x)} の(順)表現であり、ベル行列 B f ( x ) {\displaystyle f(x)} の逆表現(anti-representation)である。ここで項 f g {\displaystyle f\circ g} は函数の合成 f ( g ( x ) ) {\displaystyle f(g(x))} を意味する。

その他の性質には、次のようなものがある。

  • M [ f n ] = M [ f ] n {\displaystyle M[f^{n}]=M[f]^{n}} 。但し f n {\displaystyle f^{n}} 反復函数
  • M [ f 1 ] = M [ f ] 1 {\displaystyle M[f^{-1}]=M[f]^{-1}} 。但し f 1 {\displaystyle f^{-1}} は(カーレマン行列が可逆であるなら)逆函数

  • 定数のカーレマン行列は次で与えられる。
    M [ a ] = ( 1 0 0 a 0 0 a 2 0 0 ) {\displaystyle M[a]={\begin{pmatrix}1&0&0&\cdots \\a&0&0&\cdots \\a^{2}&0&0&\cdots \\\vdots &\vdots &\vdots &\ddots \end{pmatrix}}}
  • 恒等函数のカーレマン行列は次で与えられる。
    M x [ x ] = ( 1 0 0 0 1 0 0 0 1 ) {\displaystyle M_{x}[x]={\begin{pmatrix}1&0&0&\cdots \\0&1&0&\cdots \\0&0&1&\cdots \\\vdots &\vdots &\vdots &\ddots \end{pmatrix}}}
  • 定数の和に関するカーレマン行列は次で与えられる。
    M x [ a + x ] = ( 1 0 0 a 1 0 a 2 2 a 1 ) {\displaystyle M_{x}[a+x]={\begin{pmatrix}1&0&0&\cdots \\a&1&0&\cdots \\a^{2}&2a&1&\cdots \\\vdots &\vdots &\vdots &\ddots \end{pmatrix}}}
  • 定数倍に関するカーレマン行列は次で与えられる。
    M x [ c x ] = ( 1 0 0 0 c 0 0 0 c 2 ) {\displaystyle M_{x}[cx]={\begin{pmatrix}1&0&0&\cdots \\0&c&0&\cdots \\0&0&c^{2}&\cdots \\\vdots &\vdots &\vdots &\ddots \end{pmatrix}}}
  • 一次函数のカーレマン行列は次で与えられる。
    M x [ a + c x ] = ( 1 0 0 a c 0 a 2 2 a c c 2 ) {\displaystyle M_{x}[a+cx]={\begin{pmatrix}1&0&0&\cdots \\a&c&0&\cdots \\a^{2}&2ac&c^{2}&\cdots \\\vdots &\vdots &\vdots &\ddots \end{pmatrix}}}
  • 函数 f ( x ) = k = 1 f k x k {\displaystyle f(x)=\sum _{k=1}^{\infty }f_{k}x^{k}} のカーレマン行列は次で与えられる。
    M [ f ] = ( 1 0 0 0 f 1 f 2 0 0 f 1 2 ) {\displaystyle M[f]={\begin{pmatrix}1&0&0&\cdots \\0&f_{1}&f_{2}&\cdots \\0&0&f_{1}^{2}&\cdots \\\vdots &\vdots &\vdots &\ddots \end{pmatrix}}}
  • 函数 f ( x ) = k = 0 f k x k {\displaystyle f(x)=\sum _{k=0}^{\infty }f_{k}x^{k}} のカーレマン行列は次で与えられる。
    M [ f ] = ( 1 0 0 f 0 f 1 f 2 f 0 2 2 f 0 f 1 f 1 2 + 2 f 0 f 2 ) {\displaystyle M[f]={\begin{pmatrix}1&0&0&\cdots \\f_{0}&f_{1}&f_{2}&\cdots \\f_{0}^{2}&2f_{0}f_{1}&f_{1}^{2}+2f_{0}f_{2}&\cdots \\\vdots &\vdots &\vdots &\ddots \end{pmatrix}}}

関連項目

参考文献

  • R Aldrovandi, Special Matrices of Mathematical Physics: Stochastic, Circulant and Bell Matrices, World Scientific, 2001. (preview)
  • R. Aldrovandi, L. P. Freitas, Continuous Iteration of Dynamical Maps, online preprint, 1997.
  • P. Gralewicz, K. Kowalski, Continuous time evolution from iterated maps and Carleman linearization, online preprint, 2000.
  • K Kowalski and W-H Steeb, Nonlinear Dynamical Systems and Carleman Linearization, World Scientific, 1991. (preview)
  • D. Knuth, Convolution Polynomials arXiv online print, 1992
  • Jabotinsky, Eri: Representation of Functions by Matrices. Application to Faber Polynomials in: Proceedings of the American Mathematical Society, Vol. 4, No. 4 (Aug., 1953), pp. 546- 553 Stable jstor-URL

en:Carleman matrix