行列値関数

行列値関数（ぎょうれつちかんすう、テンプレート:Lang-en-short）とは行列を変数に持つ特殊関数の総称である^[1]^[2]。

定義

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

初等関数

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

f (z) = c_{0} + c_{1} z + \dots + c_{m} z^{m}

このときは当然ながら、

f (A) = c_{0} I + c_{1} A + \dots + c_{m} A^{m}

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

f (z) : = \sum_{k = 0}^{\infty} c_{k} z^{k}

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

f (A) : = \sum_{k = 0}^{\infty} c_{k} A^{k}

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

\exp A : = \sum_{k = 0}^{\infty} \frac{1}{k!} A^{k},

\sin A : = \sum_{k = 0}^{\infty} \frac{(- 1)^{k}}{(2 k + 1)!} A^{2 k + 1}, \cos A : = \sum_{k = 0}^{\infty} \frac{(- 1)^{k}}{(2 k)!} A^{2 k} .

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

特殊関数

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

などの関数^[8]^[9]、もしくはそのq類似についても行列バージョンを考えることができる^{[Q 1]}^{[Q 2]}^{[Q 3]}^{[Q 4]}^{[Q 5]}^{[Q 6]}。例えば、行列からなる無限乗積の収束を適切に定義したうえで^[10]^[11]、qポッホハマー記号の行列バージョンは次のように定義できる^{[Q 1]}^{[Q 2]}。

(A; q)_{n} : = \prod_{k = 0}^{n - 1} (I - A q^{k}), (A; q)_{\infty} : = \lim_{n \to \infty} (A; q)_{n}, | q | < 1, A \in ℂ^{n \times n}

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

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

工学的重要性

行列指数関数はen:exponential integratorなどの常微分方程式の数値解法において必要である他^{[A 1]}^{[A 2]}^{[A 3]}^{[A 4]}^{[A 5]}^{[A 6]}^{[A 7]}^{[A 8]}、統計学などにおいて重要視されている^[1]^[12]^[13]。このような背景の下で、数値線形代数の研究者たちは行列値関数の高精度計算^[1]^[14]^{[B 1]}・精度保証付き数値計算^[15]^[16]^[17]^[18]の研究に積極的に取り組んでいる。具体的には、以下の関数が取り組まれている。

行列指数関数^{[B 2]}^{[B 3]}^{[B 4]}^{[B 5]}^{[B 6]}
行列の平方根^{[B 7]}^{[B 8]}^{[B 9]}^{[B 10]}
行列の三角関数^{[B 11]}
行列の実数乗^{[B 12]}^{[B 13]}
行列の対数^{[B 13]}^{[B 14]}
行列のガンマ関数^{[B 15]}
行列の超幾何級数^{[B 16]}^{[B 17]}

参考文献

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.

外部リンク

Highamによる著作

行列指数関数

テンプレート:Math-stub テンプレート:Linear-algebra-stub テンプレート:Linear algebra

↑ ^1.0 ^1.1 ^1.2 ^1.3 ^1.4 ^1.5 Higham, N. J. (2008). Functions of matrices: theory and computation. SIAM.
↑ ^2.0 ^2.1 ^2.2 ^2.3 千葉克裕; 行列の関数とジョルダン標準形, 2010. サイエンティスト社
↑ Jodar, L., & Cortés, J. C. (1998). Some properties of Gamma and Beta matrix functions. Applied Mathematics Letters, 11(1), 89-93.
↑ Kontsevich, M. (1992). Intersection theory on the moduli space of curves and the matrix Airy function. en:Communications in Mathematical Physics, 147(1), 1-23.
↑ Herz, C. S. (1955). Bessel functions of matrix argument. Annals of Mathematics(2), 61, 474-523.
↑ Sinap, A., & Van Assche, W. (1996). Orthogonal matrix polynomials and applications. en:Journal of Computational and Applied Mathematics, 66(1-2), 27-52.
↑ Koev, P., & Edelman, A. (2008). Hypergeometric function of a matrix argument. Department of Mathematics, Massachusetts Institute of Technology (April 11 2008). URL https://math.mit.edu/~plamen/software/mhgref.html.
↑ Gross, K. I., & Richards, D. S. P. (1987). Special functions of matrix argument. I. Algebraic induction, zonal polynomials, and hypergeometric functions. Transactions of the American Mathematical Society, 301(2), 781-811.
↑ Abdalla, M. (2018). Special matrix functions: Characteristics, achievements and future directions. Linear and Multilinear Algebra, 1-28.
↑ Trench, W. F. (1999). Invertibly convergent infinite products of matrices. en:Journal of computational and applied mathematics, 101(1-2), 255-263.
↑ Trench, W. F. (1995). Invertibly convergent infinite products of matrices, with applications to difference equations. Computers & Mathematics with Applications, 30(11), 39-46.
↑ Leach, B. G. (1969). Bessel functions of matrix argument with statistical applications (Doctoral dissertation).
↑ James, A. T. (1975). Special functions of matrix and single argument in statistics. In Theory and Application of Special Functions (pp. 497-520). en:Academic Press.
↑ 数値線形代数の数理とHPC, 櫻井鉄也, 松尾宇泰, 片桐孝洋編（シリーズ応用数理 / 日本応用数理学会監修, 第6巻）共立出版, 2018.8
↑ Miyajima, S. (2019). Verified computation of the matrix exponential. Advances in Computational Mathematics, 45(1), 137-152.
↑ Miyajima, S. (2019). Verified computation for the matrix principal logarithm. Linear Algebra and its Applications, 569, 38-61.
↑ Miyajima, S. (2018). Fast verified computation for the matrix principal pth root. en:Journal of Computational and Applied Mathematics, 330, 276-288.
↑ Shinya Miyajima, Verified computation for the matrix Lambert W function, Applied Mathematics and Computation, Volume 362, Pages 1-15, December 2019.

引用エラー: 「Q」という名前のグループの <ref> タグがありますが、対応する <references group="Q"/> タグが見つかりません
引用エラー: 「A」という名前のグループの <ref> タグがありますが、対応する <references group="A"/> タグが見つかりません
引用エラー: 「B」という名前のグループの <ref> タグがありますが、対応する <references group="B"/> タグが見つかりません

[higham-1] 1.0 ^1.1 ^1.2 ^1.3 ^1.4 ^1.5 Higham, N. J. (2008). Functions of matrices: theory and computation. SIAM.

[chiba-2] 2.0 ^2.1 ^2.2 ^2.3 千葉克裕; 行列の関数とジョルダン標準形, 2010. サイエンティスト社

[3] Jodar, L., & Cortés, J. C. (1998). Some properties of Gamma and Beta matrix functions. Applied Mathematics Letters, 11(1), 89-93.

[4] Kontsevich, M. (1992). Intersection theory on the moduli space of curves and the matrix Airy function. en:Communications in Mathematical Physics, 147(1), 1-23.

[5] Herz, C. S. (1955). Bessel functions of matrix argument. Annals of Mathematics(2), 61, 474-523.

[6] Sinap, A., & Van Assche, W. (1996). Orthogonal matrix polynomials and applications. en:Journal of Computational and Applied Mathematics, 66(1-2), 27-52.

[7] Koev, P., & Edelman, A. (2008). Hypergeometric function of a matrix argument. Department of Mathematics, Massachusetts Institute of Technology (April 11 2008). URL https://math.mit.edu/~plamen/software/mhgref.html.

[8] Gross, K. I., & Richards, D. S. P. (1987). Special functions of matrix argument. I. Algebraic induction, zonal polynomials, and hypergeometric functions. Transactions of the American Mathematical Society, 301(2), 781-811.

[9] Abdalla, M. (2018). Special matrix functions: Characteristics, achievements and future directions. Linear and Multilinear Algebra, 1-28.

[16] Trench, W. F. (1999). Invertibly convergent infinite products of matrices. en:Journal of computational and applied mathematics, 101(1-2), 255-263.

[17] Trench, W. F. (1995). Invertibly convergent infinite products of matrices, with applications to difference equations. Computers & Mathematics with Applications, 30(11), 39-46.

[26] Leach, B. G. (1969). Bessel functions of matrix argument with statistical applications (Doctoral dissertation).

[27] James, A. T. (1975). Special functions of matrix and single argument in statistics. In Theory and Application of Special Functions (pp. 497-520). en:Academic Press.

[hpc-28] 数値線形代数の数理とHPC, 櫻井鉄也, 松尾宇泰, 片桐孝洋編（シリーズ応用数理 / 日本応用数理学会監修, 第6巻）共立出版, 2018.8

[30] Miyajima, S. (2019). Verified computation of the matrix exponential. Advances in Computational Mathematics, 45(1), 137-152.

[31] Miyajima, S. (2019). Verified computation for the matrix principal logarithm. Linear Algebra and its Applications, 569, 38-61.

[32] Miyajima, S. (2018). Fast verified computation for the matrix principal pth root. en:Journal of Computational and Applied Mathematics, 330, 276-288.

[33] Shinya Miyajima, Verified computation for the matrix Lambert W function, Applied Mathematics and Computation, Volume 362, Pages 1-15, December 2019.

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

[9]

[Q 1]

[Q 2]

[Q 3]

[Q 4]

[Q 5]

[Q 6]

[10]

[11]

[A 1]

[A 2]

[A 3]

[A 4]

[A 5]

[A 6]

[A 7]

[A 8]

[12]

[13]

[14]

[B 1]

[15]

[16]

[17]

[18]

[B 2]

[B 3]

[B 4]

[B 5]

[B 6]

[B 7]

[B 8]

[B 9]

[B 10]

[B 11]

[B 12]

[B 13]

[B 14]

[B 15]

[B 16]

[B 17]

行列値関数

目次

定義

初等関数

特殊関数

工学的重要性

関連項目

関連分野

研究者

主な行列値関数

出典

$q -$ 類似

ODEの数値計算