狠狠撸

狠狠撸Share a Scribd company logo
ランダムなトポロジー
s.t.@simizut22
2021/09/23 ロマンティック数学ナイト #16
今日の内容
1 グラフと複体
2 ランダム幾何学的複体
3 漸近挙動と相転移
4 まとめ
? s.t.@simizut22 2
自己紹介
? 清水超貴
? s.t.@simizut22
? ロマ数過去登壇:
#1 : operad
#4 : ホモトピー球面のなす群
#9 : Euler 標数と Magnitude
#16 : Random Topology
? 群馬県出身
? 現在: 吉田大学 D
? s.t.@simizut22 3
1 グラフと複体
2 ランダム幾何学的複体
3 漸近挙動と相転移
4 まとめ
? s.t.@simizut22 4
グラフと複体 グラフ
(単純無向)グラフ G とは, 頂点集合 V と辺集合 E で与えらえれるもの. これ
を G = (V, E) と書くこととする.
Figure: 頂点数 10 のグラフ
Remark
? 単純とは, 多重辺と自己ループを含まないということ
? 無向とは, 辺には特別な向きがないということ
? s.t.@simizut22 5
グラフと複体 高次元への拡張
G = (V, E) をグラフとする. G の頂点 x, y, z がどの 2 点も辺を張る時面を張る
として, 面全体の集合 F を定義する. これにより定まる K = (V, E, F) を G の旗
複体(の 2-骨格)という.
(a) グラフ G (b) 旗複体 K(G)
? s.t.@simizut22 6
位相不変量
グラフ G の旗複体 K に対し,
? G の連結な極大の部分グラフ(の旗複体)をそれぞれ連結成分という.
? 連結成分の個数 β0 を 0 次 Betti 数とという.
? G の閉路を K のサイクルという.
? いくつかの面の境界の和で与えられるサイクルとバウンダリーという.
? K のバウンダリーでないサイクルで, 線形独立なものの個数 β1(K) を 1 次
Betti 数という.
Remark
β0(K), β1(K) は連続変形(ホモトピー)で不変な値であり, いわゆる位相不変量になって
いる. したがって,
旗複体 K1, K2 に対し, β0 または β1 が異なる =? K1, K2 は位相的に異なる
という分類ができる.
? s.t.@simizut22 7
位相不変量 β0
(a) 連結グラフ
β0 = 1
(b) 連結でないグラフ
β0 = 2
? s.t.@simizut22 8
位相不変量 β1
(a) β0 = 1, β1 = 0 (b) β0 = 1, β1 = 1
(c) β0 = 1, β1 = 3
? s.t.@simizut22 9
1 グラフと複体
2 ランダム幾何学的複体
3 漸近挙動と相転移
4 まとめ
? s.t.@simizut22 10
ランダム幾何学的複体 定義
Rd
の有限集合 X = {x1, x2, . . . , xn} 及び半径 r > 0 に対し, その幾何学的グラフ
G(X, r) が以下のように定義される.
? 頂点集合 V = X
? x, y ∈ X は, その距離 kx ? yk が r 以下の時に辺をなすとする. i.e.
E = {{x, y} | x, y ∈ X, kx ? yk ≤ r}
また, この旗複体を VR(X, r) と書くこととする.
? s.t.@simizut22 11
ランダム幾何学的複体 具体例
Figure: r = 3 での旗複体
? s.t.@simizut22 12
ランダム幾何学的複体 具体例
n = 10 として, r = 1 まで発展させると以下のようになる.
? s.t.@simizut22 13
ランダム幾何学的複体 Betti 数の変化
Figure: [0, 1]4
内の 1000 点から計算した Betti 数の変化
観察
? β0 は単調減少. b1 は単調ではない.
? β0 ? 1, β1 は r が大きくなると 0 になる.
? s.t.@simizut22 14
ランダム幾何学的複体 Betti 数の変化
Figure: n = 100 から n = 1000 まで動かしたグラフの変化
観察
? グラフの形は大体似ている.
? β1 のピークが左にずれていく?
? s.t.@simizut22 15
ランダム幾何学的複体 Betti 数の変化
Q
点の個数 n に応じて, 半径 r = rn も変わったとき, β0(VR(Xn, rn)), β1(VR(Xn, rn)) は n → ∞
でどのような変化をするのか?
?
A
lim
n→∞
nrd
n =
?
?
?
?
?
?
?
?
?
?
?
0, sub-critical regime,
α ∈ (0, ∞), critical regime
∞, super-critical regime
のどれになるかに応じて, 挙動が大きく異なる. 相転移現象
? s.t.@simizut22 16
1 グラフと複体
2 ランダム幾何学的複体
3 漸近挙動と相転移
4 まとめ
? s.t.@simizut22 17
漸近挙動と相転移 設定
? x1, x2, . . . , xn, . . . は Rd
に値を持つ互いに独立な確率変数列で, 共通の確率密
度関数 f : Rd
→ R を持つものとする.
? rn を半径の列とする.
? ランダム幾何学的グラフ G(Xn, rn) を Gn , その旗複体 VR(Xn, rn) を VRn と
表す.
? k = 0, 1 として VRn の k 次数 Betti 数 βk (VRn) を βk,n で表す:
β0,n = β0(VRn) = β0(Gn), and β1,n = β1(VRn)
? s.t.@simizut22 18
漸近挙動と相転移 β0
Theorem ([Penrose et al., 2003])
1. nrd
n → 0 の亞臨界域において,
E[β0,n] ≈ n (n → ∞)
2. nrd
n → λ ∈ (0, 1) の臨界域において,
? 大数の法則が成立: 0 < c < 1 が存在して,
β0,n ≈ cn (n → ∞)
? 中心極限定理が成立:
β0,n ? E[β0,n]
√
n
d
?
→ N(0, ?
σ2
)
3. nrd
n = 2d?1
dωd
(c + log n) の時,
P(β0,n = 1) → e?e?c
? s.t.@simizut22 19
漸近挙動と相転移 亞臨界域
Theorem ([Kahle and Meckes, 2013])
limn→∞ nrd
n → 0 の亞臨界域において, 次の極限が存在:(大数の法則)
E[β1,n]
n4r3d
n
→ C1.
Theorem ([Kahle and Meckes, 2013] )
1. n4
r3d
n → 0 の時,
P(β1,n = 0) → 1
2. n4
r3d
n → α ∈ (0, ∞) の時,
β1,n
d
?
→ Poisson(C1α).
3. n4
r3d
n → ∞ の時, 次の中心極限定理が成立:
β1,n ? E[β1,n]
Var(β1,n)
d
?
→ N(0, 1).
? s.t.@simizut22 20
漸近挙動と相転移 臨界域
Theorem ([Kahle and Meckes, 2013])
limn→∞ nrd
n → c の臨界域(critical regime)において,
β1,n ~ n (n → ∞).
? s.t.@simizut22 21
漸近挙動と相転移 超臨界域
Theorem ([Kahle and Meckes, 2013])
各 xi は内点を持つ凸体 K に値を持つ i.i.d. な一様確率変数とする. limn→∞ nrd
n = ∞ の超臨
界域(super-critical regime)において, K のみに依存した定数 c > 0 が存在して,
E[β1,n] ≤ MWne?cWn
となる. ただし Wn = nrd
n . 特に,
E[β1,n]
n → 0.
Theorem ([Kahle and Meckes, 2013] )
xi は先と同様とする. この時, nrd
n ≥ c log n の連結域(connected regime)の時,
P

Kn が単連結

→ 0 (n → ∞).
? s.t.@simizut22 22
1 グラフと複体
2 ランダム幾何学的複体
3 漸近挙動と相転移
4 まとめ
? s.t.@simizut22 23
拡張
? 今回扱わなかったが, 古典的なランダムトポロジーの対象として,
Erd?s-Rényi グラフ G(n, p) [Erd?s and Rényi, 1959] がある.
? 今回は 2 次元複体までを考えたが, 一般の次元の複体の一般の次元の Betti
数に対しても同様の結果が得られている.
? 今回はグラフの旗複体に制限して話をしたが, 例えば ?ech 複体がよく扱わ
れている.
? rn を rn(t) という関数として, 確率過程としての極限を考える拡張もなされ
ている.
? Persistent Betti 数への拡張もある.
? さらに, その多変数化への拡張は... あっあっあっ.
? 基本的に, critical は人間が扱うには早すぎる.
? s.t.@simizut22 24
まとめ
Critical 何もわからん...
? s.t.@simizut22 25
Ad

Recommended

topology of musical data
topology of musical data
Tatsuki SHIMIZU
?
代数トポロジー入门
代数トポロジー入门
Tatsuki SHIMIZU
?
introductino to persistent homology and topological data analysis
introductino to persistent homology and topological data analysis
Tatsuki SHIMIZU
?
情報幾何の基礎輪読会 #1
情報幾何の基礎輪読会 #1
Tatsuki SHIMIZU
?
Magnitude ~ extend the Euler Characteristics via M?bius Inversion ~
Magnitude ~ extend the Euler Characteristics via M?bius Inversion ~
Tatsuki SHIMIZU
?
とぽろじー入门(画像なし版)
とぽろじー入门(画像なし版)
Katsuya Ito
?
Introduction to Persistence Theory
Introduction to Persistence Theory
Tatsuki SHIMIZU
?
スペクトラル?クラスタリング
スペクトラル?クラスタリング
Akira Miyazawa
?
整数格子点上の劣モジュラ被覆に対する高速アルゴリズム
整数格子点上の劣モジュラ被覆に対する高速アルゴリズム
Tasuku Soma
?
Operad and Recognition Principle
Operad and Recognition Principle
Tatsuki SHIMIZU
?
PRML 10.4 - 10.6
PRML 10.4 - 10.6
Akira Miyazawa
?
アルゴリズムイントロダクション15章 動的計画法
アルゴリズムイントロダクション15章 動的計画法
nitoyon
?
动的计画法を极める!
动的计画法を极める!
HCPC: 北海道大学競技プログラミングサークル
?
PRML 第4章
PRML 第4章
Akira Miyazawa
?
ディジタル信号処理 課題解説(その1) 2014年度版
ディジタル信号処理 課題解説(その1) 2014年度版
dsp_kyoto_2014
?
ディジタル信号処理の课题解説 その3
ディジタル信号処理の课题解説 その3
noname409
?
公開鍵暗号2: NP困難性
公開鍵暗号2: NP困難性
Joe Suzuki
?
Introduction to Topological Data Analysis
Introduction to Topological Data Analysis
Tatsuki SHIMIZU
?
ディジタル信号処理 課題解説(その3) 2014年度版
ディジタル信号処理 課題解説(その3) 2014年度版
dsp_kyoto_2014
?
PRML 第14章
PRML 第14章
Akira Miyazawa
?
Introduction to Categorical Programming (Revised)
Introduction to Categorical Programming (Revised)
Masahiro Sakai
?
ディジタル信号処理の课题解説
ディジタル信号処理の课题解説
noname409
?
ディジタル信号処理 课题解説 その4
ディジタル信号処理 课题解説 その4
noname409
?
【第34回数学カフェの予习会#1】微分と代数学のつながり
【第34回数学カフェの予习会#1】微分と代数学のつながり
MathCafe
?
自动定理証明の绍介
自动定理証明の绍介
Masahiro Sakai
?
PRML_titech 8.1 - 8.2
PRML_titech 8.1 - 8.2
Takafumi Sakakibara
?
PRML 1.6 情報理論
PRML 1.6 情報理論
sleepy_yoshi
?
8.4 グラフィカルモデルによる推論
8.4 グラフィカルモデルによる推論
sleepy_yoshi
?
渡辺澄夫著「ベイズ統計の理論と方法」5.1 マルコフ連鎖モンテカルロ法
渡辺澄夫著「ベイズ統計の理論と方法」5.1 マルコフ連鎖モンテカルロ法
Kenichi Hironaka
?
几何を使った统计のはなし
几何を使った统计のはなし
Toru Imai
?

More Related Content

What's hot (20)

整数格子点上の劣モジュラ被覆に対する高速アルゴリズム
整数格子点上の劣モジュラ被覆に対する高速アルゴリズム
Tasuku Soma
?
Operad and Recognition Principle
Operad and Recognition Principle
Tatsuki SHIMIZU
?
PRML 10.4 - 10.6
PRML 10.4 - 10.6
Akira Miyazawa
?
アルゴリズムイントロダクション15章 動的計画法
アルゴリズムイントロダクション15章 動的計画法
nitoyon
?
动的计画法を极める!
动的计画法を极める!
HCPC: 北海道大学競技プログラミングサークル
?
PRML 第4章
PRML 第4章
Akira Miyazawa
?
ディジタル信号処理 課題解説(その1) 2014年度版
ディジタル信号処理 課題解説(その1) 2014年度版
dsp_kyoto_2014
?
ディジタル信号処理の课题解説 その3
ディジタル信号処理の课题解説 その3
noname409
?
公開鍵暗号2: NP困難性
公開鍵暗号2: NP困難性
Joe Suzuki
?
Introduction to Topological Data Analysis
Introduction to Topological Data Analysis
Tatsuki SHIMIZU
?
ディジタル信号処理 課題解説(その3) 2014年度版
ディジタル信号処理 課題解説(その3) 2014年度版
dsp_kyoto_2014
?
PRML 第14章
PRML 第14章
Akira Miyazawa
?
Introduction to Categorical Programming (Revised)
Introduction to Categorical Programming (Revised)
Masahiro Sakai
?
ディジタル信号処理の课题解説
ディジタル信号処理の课题解説
noname409
?
ディジタル信号処理 课题解説 その4
ディジタル信号処理 课题解説 その4
noname409
?
【第34回数学カフェの予习会#1】微分と代数学のつながり
【第34回数学カフェの予习会#1】微分と代数学のつながり
MathCafe
?
自动定理証明の绍介
自动定理証明の绍介
Masahiro Sakai
?
PRML_titech 8.1 - 8.2
PRML_titech 8.1 - 8.2
Takafumi Sakakibara
?
PRML 1.6 情報理論
PRML 1.6 情報理論
sleepy_yoshi
?
8.4 グラフィカルモデルによる推論
8.4 グラフィカルモデルによる推論
sleepy_yoshi
?
整数格子点上の劣モジュラ被覆に対する高速アルゴリズム
整数格子点上の劣モジュラ被覆に対する高速アルゴリズム
Tasuku Soma
?
Operad and Recognition Principle
Operad and Recognition Principle
Tatsuki SHIMIZU
?
アルゴリズムイントロダクション15章 動的計画法
アルゴリズムイントロダクション15章 動的計画法
nitoyon
?
ディジタル信号処理 課題解説(その1) 2014年度版
ディジタル信号処理 課題解説(その1) 2014年度版
dsp_kyoto_2014
?
ディジタル信号処理の课题解説 その3
ディジタル信号処理の课题解説 その3
noname409
?
公開鍵暗号2: NP困難性
公開鍵暗号2: NP困難性
Joe Suzuki
?
Introduction to Topological Data Analysis
Introduction to Topological Data Analysis
Tatsuki SHIMIZU
?
ディジタル信号処理 課題解説(その3) 2014年度版
ディジタル信号処理 課題解説(その3) 2014年度版
dsp_kyoto_2014
?
Introduction to Categorical Programming (Revised)
Introduction to Categorical Programming (Revised)
Masahiro Sakai
?
ディジタル信号処理の课题解説
ディジタル信号処理の课题解説
noname409
?
ディジタル信号処理 课题解説 その4
ディジタル信号処理 课题解説 その4
noname409
?
【第34回数学カフェの予习会#1】微分と代数学のつながり
【第34回数学カフェの予习会#1】微分と代数学のつながり
MathCafe
?
自动定理証明の绍介
自动定理証明の绍介
Masahiro Sakai
?
PRML 1.6 情報理論
PRML 1.6 情報理論
sleepy_yoshi
?
8.4 グラフィカルモデルによる推論
8.4 グラフィカルモデルによる推論
sleepy_yoshi
?

Similar to ロマ数16 simizut (20)

渡辺澄夫著「ベイズ統計の理論と方法」5.1 マルコフ連鎖モンテカルロ法
渡辺澄夫著「ベイズ統計の理論と方法」5.1 マルコフ連鎖モンテカルロ法
Kenichi Hironaka
?
几何を使った统计のはなし
几何を使った统计のはなし
Toru Imai
?
A summary on “On choosing and bounding probability metrics”
A summary on “On choosing and bounding probability metrics”
Kota Matsui
?
笔搁惭尝轮讲用资料10章(パターン认识と机械学习,近似推论法)
笔搁惭尝轮讲用资料10章(パターン认识と机械学习,近似推论法)
Toshiyuki Shimono
?
モンテカルロサンプリング
モンテカルロサンプリング
Kosei ABE
?
Limit as adjointment answer
Limit as adjointment answer
tetsuya_komota
?
ある反転授业の试み:正规分布の罢补测濒辞谤展开をとおして
ある反転授业の试み:正规分布の罢补测濒辞谤展开をとおして
Hideo Hirose
?
kosenconf_Tsukuba_sciences_slide
kosenconf_Tsukuba_sciences_slide
Shigeki Nakamura
?
骋贰贰(一般化推定方程式)の理论
骋贰贰(一般化推定方程式)の理论
Koichiro Gibo
?
导来代数几何入门
导来代数几何入门
Naoya Umezaki
?
基礎からのベイズ統計学 輪読会資料 第4章 メトロポリス?ヘイスティングス法
基礎からのベイズ統計学 輪読会資料 第4章 メトロポリス?ヘイスティングス法
Ken'ichi Matsui
?
Infomation geometry(overview)
Infomation geometry(overview)
Yoshitake Misaki
?
Math20160415 epsilondelta
Math20160415 epsilondelta
Atsushi Kadotani
?
マルコフ连锁モンテカルロ法入门-2
マルコフ连锁モンテカルロ法入门-2
Nagi Teramo
?
第8回Zansa 俺の人生ランダムウォーク
第8回Zansa 俺の人生ランダムウォーク
Zansa
?
外積代数で読み解く平行体 ~究極の関係式を追い求めて~
外積代数で読み解く平行体 ~究極の関係式を追い求めて~
SoshunNaito
?
自由エネルキ?ー原理入門: 正規分布を仮定した場合
自由エネルキ?ー原理入門: 正規分布を仮定した場合
Masatoshi Yoshida
?
贵惭惭の実装と导出
贵惭惭の実装と导出
Keigo Nitadori
?
渡辺澄夫著「ベイズ統計の理論と方法」5.1 マルコフ連鎖モンテカルロ法
渡辺澄夫著「ベイズ統計の理論と方法」5.1 マルコフ連鎖モンテカルロ法
Kenichi Hironaka
?
几何を使った统计のはなし
几何を使った统计のはなし
Toru Imai
?
A summary on “On choosing and bounding probability metrics”
A summary on “On choosing and bounding probability metrics”
Kota Matsui
?
笔搁惭尝轮讲用资料10章(パターン认识と机械学习,近似推论法)
笔搁惭尝轮讲用资料10章(パターン认识と机械学习,近似推论法)
Toshiyuki Shimono
?
モンテカルロサンプリング
モンテカルロサンプリング
Kosei ABE
?
Limit as adjointment answer
Limit as adjointment answer
tetsuya_komota
?
ある反転授业の试み:正规分布の罢补测濒辞谤展开をとおして
ある反転授业の试み:正规分布の罢补测濒辞谤展开をとおして
Hideo Hirose
?
kosenconf_Tsukuba_sciences_slide
kosenconf_Tsukuba_sciences_slide
Shigeki Nakamura
?
骋贰贰(一般化推定方程式)の理论
骋贰贰(一般化推定方程式)の理论
Koichiro Gibo
?
导来代数几何入门
导来代数几何入门
Naoya Umezaki
?
基礎からのベイズ統計学 輪読会資料 第4章 メトロポリス?ヘイスティングス法
基礎からのベイズ統計学 輪読会資料 第4章 メトロポリス?ヘイスティングス法
Ken'ichi Matsui
?
Infomation geometry(overview)
Infomation geometry(overview)
Yoshitake Misaki
?
マルコフ连锁モンテカルロ法入门-2
マルコフ连锁モンテカルロ法入门-2
Nagi Teramo
?
第8回Zansa 俺の人生ランダムウォーク
第8回Zansa 俺の人生ランダムウォーク
Zansa
?
外積代数で読み解く平行体 ~究極の関係式を追い求めて~
外積代数で読み解く平行体 ~究極の関係式を追い求めて~
SoshunNaito
?
自由エネルキ?ー原理入門: 正規分布を仮定した場合
自由エネルキ?ー原理入門: 正規分布を仮定した場合
Masatoshi Yoshida
?
贵惭惭の実装と导出
贵惭惭の実装と导出
Keigo Nitadori
?
Ad

More from Tatsuki SHIMIZU (11)

TDA やら Night!!
TDA やら Night!!
Tatsuki SHIMIZU
?
エキゾチック球面ナイト(浮気編)~28 日周期の彼女たち~
エキゾチック球面ナイト(浮気編)~28 日周期の彼女たち~
Tatsuki SHIMIZU
?
Practical topology
Practical topology
Tatsuki SHIMIZU
?
Euler 標数は測度ですか??
Euler 標数は測度ですか??
Tatsuki SHIMIZU
?
しかくのお勉强
しかくのお勉强
Tatsuki SHIMIZU
?
Packing
Packing
Tatsuki SHIMIZU
?
情報幾何学の基礎 第2章 4.5
情報幾何学の基礎 第2章 4.5
Tatsuki SHIMIZU
?
Jules henri poincaré
Jules henri poincaré
Tatsuki SHIMIZU
?
Effective modern-c++#9
Effective modern-c++#9
Tatsuki SHIMIZU
?
effective modern c++ chapeter36
effective modern c++ chapeter36
Tatsuki SHIMIZU
?
emc++ chapter32
emc++ chapter32
Tatsuki SHIMIZU
?
エキゾチック球面ナイト(浮気編)~28 日周期の彼女たち~
エキゾチック球面ナイト(浮気編)~28 日周期の彼女たち~
Tatsuki SHIMIZU
?
Euler 標数は測度ですか??
Euler 標数は測度ですか??
Tatsuki SHIMIZU
?
情報幾何学の基礎 第2章 4.5
情報幾何学の基礎 第2章 4.5
Tatsuki SHIMIZU
?
effective modern c++ chapeter36
effective modern c++ chapeter36
Tatsuki SHIMIZU
?
Ad

ロマ数16 simizut

  • 2. 今日の内容 1 グラフと複体 2 ランダム幾何学的複体 3 漸近挙動と相転移 4 まとめ ? s.t.@simizut22 2
  • 3. 自己紹介 ? 清水超貴 ? s.t.@simizut22 ? ロマ数過去登壇: #1 : operad #4 : ホモトピー球面のなす群 #9 : Euler 標数と Magnitude #16 : Random Topology ? 群馬県出身 ? 現在: 吉田大学 D ? s.t.@simizut22 3
  • 4. 1 グラフと複体 2 ランダム幾何学的複体 3 漸近挙動と相転移 4 まとめ ? s.t.@simizut22 4
  • 5. グラフと複体 グラフ (単純無向)グラフ G とは, 頂点集合 V と辺集合 E で与えらえれるもの. これ を G = (V, E) と書くこととする. Figure: 頂点数 10 のグラフ Remark ? 単純とは, 多重辺と自己ループを含まないということ ? 無向とは, 辺には特別な向きがないということ ? s.t.@simizut22 5
  • 6. グラフと複体 高次元への拡張 G = (V, E) をグラフとする. G の頂点 x, y, z がどの 2 点も辺を張る時面を張る として, 面全体の集合 F を定義する. これにより定まる K = (V, E, F) を G の旗 複体(の 2-骨格)という. (a) グラフ G (b) 旗複体 K(G) ? s.t.@simizut22 6
  • 7. 位相不変量 グラフ G の旗複体 K に対し, ? G の連結な極大の部分グラフ(の旗複体)をそれぞれ連結成分という. ? 連結成分の個数 β0 を 0 次 Betti 数とという. ? G の閉路を K のサイクルという. ? いくつかの面の境界の和で与えられるサイクルとバウンダリーという. ? K のバウンダリーでないサイクルで, 線形独立なものの個数 β1(K) を 1 次 Betti 数という. Remark β0(K), β1(K) は連続変形(ホモトピー)で不変な値であり, いわゆる位相不変量になって いる. したがって, 旗複体 K1, K2 に対し, β0 または β1 が異なる =? K1, K2 は位相的に異なる という分類ができる. ? s.t.@simizut22 7
  • 8. 位相不変量 β0 (a) 連結グラフ β0 = 1 (b) 連結でないグラフ β0 = 2 ? s.t.@simizut22 8
  • 9. 位相不変量 β1 (a) β0 = 1, β1 = 0 (b) β0 = 1, β1 = 1 (c) β0 = 1, β1 = 3 ? s.t.@simizut22 9
  • 10. 1 グラフと複体 2 ランダム幾何学的複体 3 漸近挙動と相転移 4 まとめ ? s.t.@simizut22 10
  • 11. ランダム幾何学的複体 定義 Rd の有限集合 X = {x1, x2, . . . , xn} 及び半径 r > 0 に対し, その幾何学的グラフ G(X, r) が以下のように定義される. ? 頂点集合 V = X ? x, y ∈ X は, その距離 kx ? yk が r 以下の時に辺をなすとする. i.e. E = {{x, y} | x, y ∈ X, kx ? yk ≤ r} また, この旗複体を VR(X, r) と書くこととする. ? s.t.@simizut22 11
  • 12. ランダム幾何学的複体 具体例 Figure: r = 3 での旗複体 ? s.t.@simizut22 12
  • 13. ランダム幾何学的複体 具体例 n = 10 として, r = 1 まで発展させると以下のようになる. ? s.t.@simizut22 13
  • 14. ランダム幾何学的複体 Betti 数の変化 Figure: [0, 1]4 内の 1000 点から計算した Betti 数の変化 観察 ? β0 は単調減少. b1 は単調ではない. ? β0 ? 1, β1 は r が大きくなると 0 になる. ? s.t.@simizut22 14
  • 15. ランダム幾何学的複体 Betti 数の変化 Figure: n = 100 から n = 1000 まで動かしたグラフの変化 観察 ? グラフの形は大体似ている. ? β1 のピークが左にずれていく? ? s.t.@simizut22 15
  • 16. ランダム幾何学的複体 Betti 数の変化 Q 点の個数 n に応じて, 半径 r = rn も変わったとき, β0(VR(Xn, rn)), β1(VR(Xn, rn)) は n → ∞ でどのような変化をするのか? ? A lim n→∞ nrd n = ? ? ? ? ? ? ? ? ? ? ? 0, sub-critical regime, α ∈ (0, ∞), critical regime ∞, super-critical regime のどれになるかに応じて, 挙動が大きく異なる. 相転移現象 ? s.t.@simizut22 16
  • 17. 1 グラフと複体 2 ランダム幾何学的複体 3 漸近挙動と相転移 4 まとめ ? s.t.@simizut22 17
  • 18. 漸近挙動と相転移 設定 ? x1, x2, . . . , xn, . . . は Rd に値を持つ互いに独立な確率変数列で, 共通の確率密 度関数 f : Rd → R を持つものとする. ? rn を半径の列とする. ? ランダム幾何学的グラフ G(Xn, rn) を Gn , その旗複体 VR(Xn, rn) を VRn と 表す. ? k = 0, 1 として VRn の k 次数 Betti 数 βk (VRn) を βk,n で表す: β0,n = β0(VRn) = β0(Gn), and β1,n = β1(VRn) ? s.t.@simizut22 18
  • 19. 漸近挙動と相転移 β0 Theorem ([Penrose et al., 2003]) 1. nrd n → 0 の亞臨界域において, E[β0,n] ≈ n (n → ∞) 2. nrd n → λ ∈ (0, 1) の臨界域において, ? 大数の法則が成立: 0 < c < 1 が存在して, β0,n ≈ cn (n → ∞) ? 中心極限定理が成立: β0,n ? E[β0,n] √ n d ? → N(0, ? σ2 ) 3. nrd n = 2d?1 dωd (c + log n) の時, P(β0,n = 1) → e?e?c ? s.t.@simizut22 19
  • 20. 漸近挙動と相転移 亞臨界域 Theorem ([Kahle and Meckes, 2013]) limn→∞ nrd n → 0 の亞臨界域において, 次の極限が存在:(大数の法則) E[β1,n] n4r3d n → C1. Theorem ([Kahle and Meckes, 2013] ) 1. n4 r3d n → 0 の時, P(β1,n = 0) → 1 2. n4 r3d n → α ∈ (0, ∞) の時, β1,n d ? → Poisson(C1α). 3. n4 r3d n → ∞ の時, 次の中心極限定理が成立: β1,n ? E[β1,n] Var(β1,n) d ? → N(0, 1). ? s.t.@simizut22 20
  • 21. 漸近挙動と相転移 臨界域 Theorem ([Kahle and Meckes, 2013]) limn→∞ nrd n → c の臨界域(critical regime)において, β1,n ~ n (n → ∞). ? s.t.@simizut22 21
  • 22. 漸近挙動と相転移 超臨界域 Theorem ([Kahle and Meckes, 2013]) 各 xi は内点を持つ凸体 K に値を持つ i.i.d. な一様確率変数とする. limn→∞ nrd n = ∞ の超臨 界域(super-critical regime)において, K のみに依存した定数 c > 0 が存在して, E[β1,n] ≤ MWne?cWn となる. ただし Wn = nrd n . 特に, E[β1,n] n → 0. Theorem ([Kahle and Meckes, 2013] ) xi は先と同様とする. この時, nrd n ≥ c log n の連結域(connected regime)の時, P Kn が単連結 → 0 (n → ∞). ? s.t.@simizut22 22
  • 23. 1 グラフと複体 2 ランダム幾何学的複体 3 漸近挙動と相転移 4 まとめ ? s.t.@simizut22 23
  • 24. 拡張 ? 今回扱わなかったが, 古典的なランダムトポロジーの対象として, Erd?s-Rényi グラフ G(n, p) [Erd?s and Rényi, 1959] がある. ? 今回は 2 次元複体までを考えたが, 一般の次元の複体の一般の次元の Betti 数に対しても同様の結果が得られている. ? 今回はグラフの旗複体に制限して話をしたが, 例えば ?ech 複体がよく扱わ れている. ? rn を rn(t) という関数として, 確率過程としての極限を考える拡張もなされ ている. ? Persistent Betti 数への拡張もある. ? さらに, その多変数化への拡張は... あっあっあっ. ? 基本的に, critical は人間が扱うには早すぎる. ? s.t.@simizut22 24