NSM-029

行列指数関数の数学的基盤

作成日: 2026-06-03 / 正方行列・行列ノルム・バナッハ代数・マクローリン展開・コンパクト集合 / Mathematics in Neuroscience プレゼン準備

このノートの位置づけ:
神経科学のマルコフ連鎖モデルでは $e^{Qt}$ という行列指数関数が登場する。このノートは「なぜ $e^{Qt}$ が定義できるのか・なぜ収束するのか・なぜ微分できるのか」を理解するための数学的基盤を、関数解析の予備知識なしで段階的に解説する。
対象読者: 大学学部レベルの線形代数は知っているが、関数解析は未学習の神経科学研究者。
【本ノートでの用語規約】
ノルム(norm): ベクトル空間の各元に「大きさ」を割り当てる関数 $\|\cdot\|$。距離の概念を一般化する。行列への拡張が「行列ノルム」。
完備性(completeness): ノルム空間において「コーシー列が必ず収束する」性質。実数列が「$\varepsilon$-差が 0 に近づけば収束する」のと同様の概念。
劣乗法性(submultiplicativity): $\|AB\| \leq \|A\|\,\|B\|$ という行列ノルムの性質。「積のノルムは各ノルムの積以下」。「積との互換性」とも呼ぶ。

1. 正方行列の集合 $M_n(\mathbb{C})$

教科書確認済み

$n \times n$ 行列($n$ 行 $n$ 列で行数と列数が等しい行列)を 正方行列 という。$n \times n$ の複素正方行列全体の集合を $M_n(\mathbb{C})$ と書く。

線形空間としての構造

$M_n(\mathbb{C})$ には次の 2 つの演算が定義される:

これにより $M_n(\mathbb{C})$ は $n^2$ 次元の線形空間(ベクトル空間)になる。$n = 3$ なら $3^2 = 9$ 次元の空間だ。

行列の積(バナッハ「代数」になる鍵)

さらに 行列の積 $AB$ が定義される($(AB)_{ij} = \sum_k A_{ik} B_{kj}$)。積が定義される線形空間を 代数 という。この「代数」という構造が後で登場する「バナッハ代数」の核心となる。

なぜ正方行列が必要か: $e^{Qt}$ を定義するには $Q^2, Q^3, \ldots$ という積を繰り返す必要がある。$A \in M_n(\mathbb{C})$ ならば $A \cdot A$ が同じ空間 $M_n(\mathbb{C})$ に属する(「閉じている」)ことが保証される。正方でない行列では積が閉じない。

具体例: 2状態イオンチャネルの Q 行列

Closed ⇌ Open の 2 状態モデル(開口速度 $\alpha$、閉口速度 $\beta$)の Q 行列は $M_2(\mathbb{R}) \subset M_2(\mathbb{C})$ の元:

$$Q = \begin{pmatrix} -\alpha & \alpha \\ \beta & -\beta \end{pmatrix}$$

行和が 0($-\alpha + \alpha = 0$, $\beta + (-\beta) = 0$)になることは Q 行列の基本性質。$Q \in M_2(\mathbb{R})$ は 2×2 の正方行列なので、$Q^2, Q^3, \ldots$ が全て $M_2(\mathbb{R})$ に属することが保証される。

M_n(C) の構造 — 3 つの演算が同時に定義される 和・スカラー倍 A + B, αA → 線形空間 (n² 次元) 行列の積 AB, A², A³, … → 代数の構造 (べき乗が閉じる) ノルム (距離) ‖A‖ ∈ ℝ≥0 → 収束を測る (次節で定義) 3 つ全てを備えた構造 → バナッハ代数(§3 で解説)
図 1. $M_n(\mathbb{C})$ が持つ 3 つの構造。この 3 つが揃うことでバナッハ代数が構成される。

2. 行列ノルム

教科書確認済み — Rudin, Functional Analysis, Chapter 1; Horn & Johnson, Matrix Analysis, §5.6

ノルムとは何か

ノルムはベクトルや行列の「大きさ」を測る関数 $\|\cdot\|$ であり、次の 3 条件を満たす:

  1. $\|A\| \geq 0$、かつ $\|A\| = 0 \Leftrightarrow A = O$(零行列のみ大きさが 0)
  2. $\|\alpha A\| = |\alpha|\,\|A\|$(スカラー倍に比例)
  3. $\|A + B\| \leq \|A\| + \|B\|$(三角不等式)

行列ノルムの定義(誘導ノルム)

ベクトルノルム $\|v\|$(例: ユークリッドノルム $\|v\| = \sqrt{v_1^2 + \cdots + v_n^2}$)が与えられたとき、行列 $A$ の 誘導ノルム を次で定義する:

$$\|A\| = \sup_{v \neq 0} \frac{\|Av\|}{\|v\|}$$

これは「$A$ が単位ベクトルを最大でどれだけ引き伸ばすか」を測る。$\|A\| = 2$ であれば、$A$ は任意のベクトルを長さの 2 倍以内に引き伸ばすことを意味する。

フロベニウスノルム(もう 1 つの選択肢):
$$\|A\|_F = \sqrt{\sum_{i,j} |A_{ij}|^2}$$ 成分を全て 2 乗して和の平方根を取る。計算が簡単で、誘導ノルムと同等の性質を持つ。本ノートの議論はどちらのノルムでも成立する。

最重要性質: 劣乗法性

命題(行列ノルムの劣乗法性)
$$\|AB\| \leq \|A\|\,\|B\|$$

積のノルムは、各ノルムの積以下である。

直感: $A$ が「最大 $\|A\|$ 倍の引き伸ばし」、$B$ が「最大 $\|B\|$ 倍の引き伸ばし」なら、$AB$ は「最大 $\|A\|\cdot\|B\|$ 倍の引き伸ばし」を超えない。

なぜこの性質が必要か

$e^{At}$ のマクローリン展開には $\frac{(At)^k}{k!}$ の和が登場する。各項のノルムを抑えるために:

$$\|(At)^k\| = \|At \cdot At \cdots At\| \leq \|At\|^k = \bigl(\|A\|\,|t|\bigr)^k$$

劣乗法性を $k$ 回適用することでこの不等式が得られる。これが「無限和が収束する」ことの鍵となる(§4 で詳述)。

劣乗法性 ‖AB‖ ≤ ‖A‖‖B‖ の直感 v ‖v‖ = 1 B 最大 ‖B‖ 倍 Bv ‖Bv‖ ≤ ‖B‖ A 最大 ‖A‖ 倍 ABv ‖ABv‖ ≤ ‖A‖‖B‖ ∴ ‖AB‖ ≤ ‖A‖‖B‖ ‖ABv‖ ≤ ‖A‖‖Bv‖ ≤ ‖A‖‖B‖‖v‖ = ‖A‖‖B‖(‖v‖=1 の場合)→ sup を取って ‖AB‖ ≤ ‖A‖‖B‖
図 2. 劣乗法性の直感。B が「最大 ‖B‖ 倍」、A が「最大 ‖A‖ 倍」伸ばすなら、合成 AB は「最大 ‖A‖‖B‖ 倍」を超えない。

3. バナッハ代数

教科書確認済み — Rudin, Functional Analysis, Chapter 10; Conway, A Course in Functional Analysis, §VIII.1

バナッハ空間: ノルム空間 + 完備性

定義(バナッハ空間)

ノルム空間とはノルム $\|\cdot\|$ を備えた線形空間。さらに 完備(complete)であるとは、全ての コーシー列 が収束することをいう。完備なノルム空間を バナッハ空間 という。

コーシー列: $m, n \to \infty$ のとき $\|A_m - A_n\| \to 0$ となる列 $\{A_n\}$。

完備性の直感: 「どんどん近づいていく点列には、必ず近づいていく先がある」という性質。実数 $\mathbb{R}$ は完備だが、有理数 $\mathbb{Q}$ は完備でない($\sqrt{2}$ に収束するコーシー列の極限が $\mathbb{Q}$ に存在しない)。

バナッハ代数: バナッハ空間 + 積との互換

定義(バナッハ代数)

バナッハ空間 $(X, \|\cdot\|)$ に積演算 $\cdot : X \times X \to X$ が定義され、次の条件を満たすとき バナッハ代数 という:

  1. 積は結合的: $(AB)C = A(BC)$
  2. 積は分配的: $A(B+C) = AB + AC$
  3. 劣乗法性: $\|AB\| \leq \|A\|\,\|B\|$

$M_n(\mathbb{C})$ がバナッハ代数である理由

必要な性質 $M_n(\mathbb{C})$ での根拠
ノルム空間 誘導ノルム(またはフロベニウスノルム)が定義できる
完備性 有限次元ノルム空間は常に完備($\mathbb{R}^{n^2}$ と同型なので)
積演算 行列の積 $AB$ が定義される
劣乗法性 誘導ノルムは定義から $\|AB\| \leq \|A\|\,\|B\|$ を満たす
核心: $M_n(\mathbb{C})$ は有限次元なので完備性は「自動的に」成立する。完備性の証明は必要ない。この「自動性」が有限次元行列の扱いを大きく簡単にしている。

4. マクローリン展開による定義

教科書確認済み — Horn & Johnson, Matrix Analysis, §2.10; Higham, Functions of Matrices, Chapter 10

スカラーの場合

実数 $x$ に対するテイラー展開(マクローリン展開):

$$e^x = \sum_{k=0}^{\infty} \frac{x^k}{k!} = 1 + x + \frac{x^2}{2} + \frac{x^3}{6} + \cdots$$

この級数は全ての $x \in \mathbb{R}$ で絶対収束する。

行列への拡張

$A \in M_n(\mathbb{C})$、$t \in \mathbb{R}$ に対して:

$$e^{At} = \sum_{k=0}^{\infty} \frac{(At)^k}{k!} = I + At + \frac{(At)^2}{2!} + \frac{(At)^3}{3!} + \cdots$$

ここで $(At)^0 = I$(単位行列)。この形がスカラーと完全に同じ形をしているのは偶然ではなく、バナッハ代数では「絶対収束する冪級数が定義できる」からだ。

なぜ収束するか

部分和 $S_N = \sum_{k=0}^{N} \frac{(At)^k}{k!}$ がコーシー列になることを示す。劣乗法性より $\|(At)^k\| \leq \|At\|^k$ なので:

$$\|S_N - S_M\| \leq \sum_{k=M+1}^{N} \frac{\|(At)^k\|}{k!} \leq \sum_{k=M+1}^{N} \frac{\|At\|^k}{k!}$$

右辺は $e^{\|At\|}$(スカラーの指数関数)の部分和の差なので $N, M \to \infty$ で $0$ に収束する。よって $\{S_N\}$ はコーシー列。$M_n(\mathbb{C})$ は完備(バナッハ空間)なので、コーシー列は必ず収束する。$\square$

バナッハ代数の役割の整理: この 2 点がなければ、$e^{At}$ という行列が「存在する」ことすら言えない。

部分和の具体例(2×2 の場合)

$Q = \begin{pmatrix} -\alpha & \alpha \\ \beta & -\beta \end{pmatrix}$ で $t = 0.1$、$\alpha = 3$、$\beta = 7$ とする。

第0項 (k=0):  I = ((1, 0), (0, 1))

第1項 (k=1):  Qt = ((-0.3, 0.3), (0.7, -0.7))

第2項 (k=2):  (Qt)²/2! : (Qt)² = ((-0.3)²+0.3×0.7, …) を計算
              ≈ (( 0.03, -0.03), (-0.07,  0.07)) / ... (実際は行列計算)

第3項以降:    急速に小さくなる(‖Qt‖ = 1.0 程度なので各項は 1/k! で減衰)

数値計算では第 10〜15 項程度でほぼ収束する。

5. コンパクト集合と一様絶対収束

教科書確認済み — Rudin, Principles of Mathematical Analysis, Chapter 7 (一様収束); Rudin, Real and Complex Analysis, Theorem 10.28 (項別微分)

コンパクト集合とは

$\mathbb{R}$ のコンパクト集合の最も身近な例は 有界閉区間 $[0, T]$ である。時間 $t$ を考えるときに $t \in [0, T]$(今から $T$ 秒後まで)と限定するのが典型的な状況。

コンパクト集合(直感的定義)

$\mathbb{R}^n$ のコンパクト集合とは 有界かつ閉じた集合(Heine-Borel の定理)。直感的には「全体が有限の大きさに収まっており、端点も含む」集合。

各点収束と一様収束の違い

無限級数 $e^{At} = \sum_{k=0}^{\infty} \frac{(At)^k}{k!}$ が $t$ の関数として収束するとき、2 種類の収束がある:

各点収束 vs 一様収束(区間 [0, T] 上) 各点収束(点ごとに別々に収束) 各 t を固定すると S_N(t) → e^{At} ただし「何項で収束するか」は t により異なる t ‖e‖ t が大きいほど収束が遅い(誤差幅が変わる) 一様収束(区間全体で同時に収束) 全ての t ∈ [0, T] で同じ N 項あれば十分 誤差 = sup_{t ∈ [0,T]} ‖S_N(t) − e^{At}‖ → 0 t ‖e‖ 一定幅 ε どの t でも N 項以降は誤差が ε 未満(N は t によらない)
図 3. 各点収束(左)と一様収束(右)の違い。一様収束では「全ての $t \in [0,T]$ で同じ項数 $N$ があれば収束が保証される」。

なぜコンパクト集合上で一様絶対収束するか

$t \in [0, T]$(コンパクト集合)に限定すると:

$$\sum_{k=0}^{\infty} \frac{\|(At)^k\|}{k!} \leq \sum_{k=0}^{\infty} \frac{(\|A\|\,|t|)^k}{k!} \leq \sum_{k=0}^{\infty} \frac{(\|A\| T)^k}{k!} = e^{\|A\| T} < \infty$$

右辺 $e^{\|A\|T}$ は $T$ に依存するが、$t$ には依存しない定数である。Weierstrass M-test(各項のノルムの和が収束する定数で抑えられるなら、級数は一様絶対収束する)を適用することで、$e^{At}$ は $[0, T]$ 上で一様絶対収束することが示される。

なぜ「コンパクト集合」が必要か:
$t \in [0, \infty)$(非コンパクト)では $e^{\|A\|T}$ の $T \to \infty$ 極限が発散するため、全区間での一様収束は保証されない。$[0, T]$ という有界集合に限定することで、上界 $e^{\|A\|T}$ が有限に保たれる。これが「コンパクト集合上の一様収束」という表現の実質的な意味。

一様収束がなぜ重要か: 項別微分の正当化

一様絶対収束が確立されると、次の項別微分が正当化される(Rudin, Principles, Theorem 7.17):

$$\frac{d}{dt} e^{At} = \frac{d}{dt} \sum_{k=0}^{\infty} \frac{(At)^k}{k!} = \sum_{k=1}^{\infty} \frac{k A (At)^{k-1}}{k!} = A \sum_{k=0}^{\infty} \frac{(At)^k}{k!} = A e^{At}$$

これにより前進方程式 $\frac{dP}{dt} = QP$ の解が $P(t) = e^{Qt}$ であることが数学的に確立される。

6. 全体のつながり — 論理の流れ

行列指数関数が定義・微分できるまでの論理連鎖 M_n(ℂ) に行列ノルムを入れる ‖A‖ = sup{ ‖Av‖/‖v‖ : v ≠ 0 } 有限次元 → 完備性自動 バナッハ代数が構成される 完備ノルム空間 + 積演算 + ‖AB‖ ≤ ‖A‖‖B‖ 劣乗法性 → 各項を抑える マクローリン展開が絶対収束 e^{At} = Σ (At)^k/k! は M_n(ℂ) の元として収束する t ∈ [0,T] → 上界 e^{‖A‖T} が有限 コンパクト集合上で一様絶対収束 sup_{t∈[0,T]} Σ ‖(At)^k‖/k! ≤ e^{‖A‖T} < ∞ Weierstrass M-test + 項別微分定理 項別微分が正当化される d/dt e^{At} = A e^{At} → 前進方程式 dP/dt = QP が解ける 各ステップが次のステップを可能にする — 飛ばせるステップは 1 つもない
図 4. 行列指数関数が「定義でき・微分できる」ための論理連鎖。5 段階が順番に積み重なる。

7. まとめ

各概念の役割を 1 行で整理する:

概念 $e^{Qt}$ における役割 なければどうなるか
$M_n(\mathbb{C})$ $e^{Qt}$ の「住処」。積が閉じた線形空間 $Q^k$ が同じ空間に属することを言えない
行列ノルム 収束を測る道具。$\|(At)^k\| \leq \|At\|^k$ を与える 各項の大きさを比較できない
バナッハ代数 「コーシー列には収束先がある」の保証 $e^{Qt}$ という行列が存在するか不明
マクローリン展開 $e^{Qt}$ の定義そのもの (これが出発点)
コンパクト集合 一様絶対収束と項別微分の前提 「$t$ によらず同じ精度」が保証されない
ひとことまとめ:
$M_n(\mathbb{C})$ はバナッハ代数なので、スカラーの $e^x$ のマクローリン展開を行列版 $e^{At}$ にそのまま持ち込める。コンパクト集合 $[0,T]$ 上に限定することで一様絶対収束が確立し、項別微分が正当化される。これにより $\frac{d}{dt}e^{Qt} = Qe^{Qt}$ が成り立ち、前進方程式 $\frac{dP}{dt} = QP$ の解 $P(t) = e^{Qt}$ が数学的に確立される。

参考文献

関連項目


NSM-029 / 作成: 2026-06-03 / Mathematics in Neuroscience プレゼン準備