行列値関数(ぎょうれつちかんすう、英: matrix-valued functions)とは行列を変数に持つ特殊関数の総称である。

定義

ガンマ関数などを除けば、通常の特殊関数は多くの場合に常微分方程式の解 (可積分系の厳密解) として定義される。しかし行列値関数の場合は異なる。

初等関数

A C n × n , f : C C {\displaystyle \forall A\in \mathbb {C} ^{n\times n},\quad f:\mathbb {C} \to \mathbb {C} } が与えられたとする。このとき f ( z ) {\displaystyle f(z)} がどのような関数であれば f ( A ) {\displaystyle f(A)} に行列としての意味を持たせられるかを考える。自然に思いつくのは多項式の場合:

f ( z ) = c 0 c 1 z c m z m {\displaystyle f(z)=c_{0} c_{1}z \cdots c_{m}z^{m}}

このときは当然ながら、

f ( A ) = c 0 I c 1 A c m A m {\displaystyle f(A)=c_{0}I c_{1}A \cdots c_{m}A^{m}}

と定義するのが合理的である。この考えを発展させることで

f ( z ) := k = 0 c k z k {\displaystyle f(z):=\sum _{k=0}^{\infty }c_{k}z^{k}}

と定義されているときには

f ( A ) := k = 0 c k A k {\displaystyle f(A):=\sum _{k=0}^{\infty }c_{k}A^{k}}

と定義すればよいということが言える (もちろん行列からなる無限列の収束を適切に定義することも必要不可欠である)。例えば行列指数関数などの初等関数は次のように定められる:

exp A := k = 0 1 k ! A k , {\displaystyle \exp A:=\sum _{k=0}^{\infty }{\frac {1}{k!}}A^{k},}
sin A := k = 0 ( 1 ) k ( 2 k 1 ) ! A 2 k 1 , cos A := k = 0 ( 1 ) k ( 2 k ) ! A 2 k . {\displaystyle \sin A:=\sum _{k=0}^{\infty }{\frac {(-1)^{k}}{(2k 1)!}}A^{2k 1},\quad \cos A:=\sum _{k=0}^{\infty }{\frac {(-1)^{k}}{(2k)!}}A^{2k}.}

もしも f ( z ) {\displaystyle f(z)} がベキ級数表示を持たない場合はLagrange-Sylvester 多項式という道具を使って f ( A ) {\displaystyle f(A)} を定めることができる。

特殊関数

代表的な特殊関数、具体的には

  • ガンマ関数
  • エアリー関数
  • ベッセル関数
  • 直交多項式
  • 超幾何級数

などの関数、もしくはそのq類似についても行列バージョンを考えることができる。例えば、行列からなる無限乗積の収束を適切に定義したうえで、qポッホハマー記号の行列バージョンは次のように定義できる。

( A ; q ) n := k = 0 n 1 ( I A q k ) , ( A ; q ) := lim n ( A ; q ) n , | q | < 1 , A C n × n {\displaystyle (A;q)_{n}\!:=\!\prod _{k=0}^{n-1}(I-Aq^{k}),\quad (A;q)_{\infty }\!:=\!\lim _{n\to \infty }(A;q)_{n},\quad |q|<1,\quad A\in \mathbb {C} ^{n\times n}}

これを使って、en:q-gamma functionの行列バージョンも導入できる。

Γ q ( A ) := ( q ; q ) ( q A ; q ) 1 ( 1 q ) I A , | q | < 1 {\displaystyle \Gamma _{q}(A):=(q;q)_{\infty }(q^{A};q)_{\infty }^{-1}(1-q)^{I-A},\quad |q|<1}

工学的重要性

行列指数関数はen:exponential integratorなどの常微分方程式の数値解法において必要である他、統計学などにおいて重要視されている。このような背景の下で、数値線形代数の研究者たちは行列値関数の高精度計算・精度保証付き数値計算の研究に積極的に取り組んでいる。具体的には、以下の関数が取り組まれている。

  • 行列指数関数
  • 行列の平方根
  • 行列の三角関数
  • 行列の実数乗
  • 行列の対数
  • 行列のガンマ関数
  • 行列の超幾何級数

関連項目

関連分野

  • 行列解析 (en) (線型代数学の中で行列値関数の性質を調べる分科)
  • 数値線形代数

研究者

  • クリーブ・モラー (en)
  • ニコラス・ハイアム (en) (行列値関数の研究で多くの業績がある)

主な行列値関数

  • 行列指数関数
  • 行列の対数
  • 行列の平方根

出典

q {\displaystyle q-} 類似

ODEの数値計算

行列値関数の高精度計算

参考文献

  • Higham, N. J. (2006). Functions of matrices. Manchester Institute for Mathematical Sciences, School of Mathematics, The University of Manchester.
  • Higham, N. J. (2002). The matrix computation toolbox.
  • A Survey of the Matrix Exponential Formulae with Some Applications (2016), Baoying Zheng, Lin Zhang, Minhyung Cho, and Junde Wu. J. Math. Study Vol. 49, No. 4, pp. 393-428.

外部リンク

  • 行列の指数関数とその性質
  • 行列関数の数値計算アルゴリズム開発
  • 行列引数超幾何関数の数値計算と holonomic D加群、researchmapより
  • ウィシャート分布に現れる行列変数の超幾何関数に対するホロノミック勾配法 (PDF)
  • Computation of matrix gamma function (PDF)

Highamによる著作

  • Matrix Functions and their Sensitivity (PDF)
  • The Matrix Logarithm: from Theory to Computation (PDF)
  • Functions of a Matrix: Theory, Applications and Computation (PDF)
  • Challenges in Multivalued Matrix Functions (PDF)

行列指数関数

  • expmv - GitHub
  • How and How Not to Compute the Exponential of a Matrix (PDF)
  • Charlie Van Loan and the Matrix Exponential (PDF)

行列の指数関数 リレーエッセイ 静岡産業大学|経営・スポーツ・静岡県の大学

行列の指数関数 リレーエッセイ 静岡産業大学|経営・スポーツ・静岡県の大学

第5回 配列. ppt download

INDEX関数の構文 指定した行列の値を取得しよう 配列形式編 DL可 データ分析ドットコム

約数関数を成分に持つ行列の行列式 Mathlog