logo

パーシステント・モジュール 📂位相データ分析

パーシステント・モジュール

定義 1

20220710_132230.png

RRとし、拡張された整数集合 Z\mathbb{Z}Z:=Z{±}\overline{\mathbb{Z}} := \mathbb{Z} \cup \left\{ \pm \infty \right\}としよう。

  1. 次のように鎖複体群Ci\mathsf{C}_{\ast}^{i}間のチェーンマップchain Mapfi:CiCi+1f^{i} : \mathsf{C}_{\ast}^{i} \to \mathsf{C}_{\ast}^{i+1}が与えられた場合、C:={(Ci,i,fi):i0}\mathcal{C} := \left\{ \left( \mathsf{C}_{\ast}^{i} , \partial_{i} , f^{i} \right) : i \ge 0 \right\}持続複体群persistent Complexと言う。 C0f0C1f1C2f2 \mathsf{C}_{\ast}^{0} \overset{f^{0}}{\longrightarrow} \mathsf{C}_{\ast}^{1} \overset{f^{1}}{\longrightarrow} \mathsf{C}_{\ast}^{2} \overset{f^{2}}{\longrightarrow} \cdots
  2. RR-モジュール MiM^{i}間のホモモルフィズム φi:MiMi+1\varphi^{i} : M^{i} \to M^{i+1}が与えられた場合、M:={(Mi,φi):i0}\mathcal{M} := \left\{ \left( M^{i} , \varphi^{i} \right) : i \ge 0 \right\}持続モジュールpersistent Moduleと言う。
  3. 持続複体群C\mathcal{C}や持続モジュールM\mathcal{M}の各コンポーネントCi\mathsf{C}_{\ast}^{i}MiM^{i}が、RR-モジュール有限生成であり、あるmZm \in \mathbb{Z}に対して全てのfif^{i}φi\varphi^{i}imi \ge m同型である場合、C\mathcal{C}M\mathcal{M}有限型と言う。
  4. 0ijZ0 \le i \le j \in \overline{Z}を満たす区間(i,j)(i,j)P\mathcal{P}-インターバルp-intervalと言う。

説明

これは何の馬鹿げた定義だ?

これらの定義はゾモロディアンzomorodianの論文から引用したもので、見ての通り、三つの定義は関係がなさそうだ。彼の論文を読んでいると、この部分が一番問題だった。こういった定義が次々に紹介されるので、何か関係があるはずなのに、どうしてもただの定義同士に関係が見えない。持続複体群は「持続」と全く関係がなく、P\mathcal{P}-インターバルにはP\mathcal{P}が出てこない。さらに言えば、持続モジュールは「モジュール」ではなく「モジュールのファミリー」なのに、単にモジュールと呼ばれている。

20220710_143709.png

これら自体の意味にこだわると、純度100%の精神病になる。誰かが「実際には全く関係がない」と一言教えてくれれば、もっと時間を節約して、快適に学習できただろう。定義はそのまま受け入れて、実際に重要な議論に移ろう。

トポロジカルデータ解析

基本的にトポロジカルデータ解析topological Data analysisでは、フィルター付き複体群でどのようにいくつかのトポロジカルな性質が持続性persistencyを持つかに興味を持つ。例えば、ビエトリス-リプス複体群で半径ε>0\varepsilon > 0を増やしながらさまざまな複体群を列挙してフィルター付き複体群を作れば、その結果の変化がまさにfif^{i}で表されることになる。p\partial_{p}はフィルター付き複体群に関係なく、各複体群に伴うべき境界マップなので、この文脈では気にする必要がない。これは、TDAの直感的なプロセスを数学的表現に移したものであり、持続モジュールは特に取り扱いやすい条件を満たせるようにした持続複体群の中でもCi=Mi\mathsf{C}_{\ast}^{i} = M^{i}でありfi=φif^{i} = \varphi^{i}である場合を想定しているのだ。この文脈を理解すれば、そういった定義の紹介もそんなに変ではない。問題は、そういった良い条件を与えられた持続モジュールM\mathcal{M}もまだ扱いづらいということだが、ここでグレードモジュールが彗星のように現れる。

M\mathcal{M}RR上での持続モジュールであるとき、R[t]R[t]の文脈で α(M):=i=0Mi \alpha \left( \mathcal{M} \right) := \bigoplus_{i=0}^{\infty} M^{i} 群作用t(m0,m1,m2,)=(0,φ0(m0),φ1(m1),φ2(m2),) t \cdot \left( m^{0} , m^{1} , m^{2} , \cdots \right) = \left( 0, \varphi^{0} \left( m^{0} \right) , \varphi^{1} \left( m^{1} \right) , \varphi^{2} \left( m^{2} \right) , \cdots \right) 与えることができる。つまり、ttが作用するとか、かけられるというのは、MiM^{i}からMi+1M^{i+1}に要素を昇格させるということだ。この対応関係α\alphaは、次の定理に従って、持続モジュールの問題をグレードモジュールの問題に移す。

可換代数学のアーティン-リース理論: 対応関係α\alphaは、RR上の有限型持続ホモロジーのカテゴリとR[t]R[t]上の有限生成標準グレードモジュールのカテゴリとの間のカテゴリ同等性を定義する。

グレードモジュールは、言葉が少し複雑だが、我々が学部時代によく扱った多項式環に過ぎないかもしれず、それよりもさらにシンプルな形であるかもしれない。これらの構築を受け入れることができれば、もはや抽象的な代数トポロジーを超えて、ホモロジーというものを計算するアルゴリズムの領域に進む準備ができたということだ。 Q(i,j):=iF[t]/(tji) Q (i,j) := \sum^{i} F[t] / \left( t^{j-i} \right) 我々は、FFに対してグレード環でありPIDであるF[t]F[t]P\mathcal{P}-インターバル群S\mathcal{S}と次のような全射QQで関連付けることになる。直感的な例を想像してみると、我々のデータがε5=0.6\varepsilon_{5} = 0.6からε12=1.5\varepsilon_{12} = 1.5までの11番目のベッチ数β1=2\beta_{1} = 2だった場合、次のような関係を見つけたことになる。 Q(5,12)=iF[t]/(t125) Q (5,12) = \sum^{i} F[t] / \left( t^{12-5} \right) そして、ベッチ数で説明されるトポロジカルな性質が持続された期間は1.50.6=0.91.5 - 0.6 = 0.9である。明示的ではないが、この話を聞いてみれば、持続複体群などになぜ「持続」persistentという表現がつけられたのか、少しは共感できるようになるだろう。


  1. Zomorodian. (2005). Computing Persistent Homology: ch3 ↩︎