曲面のホモロジー群

[size=150][b]このワークシートは[url=https://www.geogebra.org/m/twxxx3yq]Math by Code[/url]の一部です。[br][br][/b][size=100]前回は、図形を単純な単体からできたものと考えたり、[br]単体を基底にしてつないだ鎖としての自由なたし算の群、鎖群を考えたり、[br]鎖群の境界を求めたりしたね。[br][br]図形のつながり方が同じときに同相といった。[br][/size][size=100]同相ならば、[color=#0000ff]オイラー指標が一致する[/color]ことは単体にして数えたり、規則性から推理できる。[br]同相ならば、変わらない量を不変量と呼ぶ。[br][b]同相ならば、同型になるものがある。それが、[color=#0000ff]ホモロジー群[/color]だ。[br][/b][color=#ff0000]※これは同相の条件ではなく、性質なので注意したい。[br][/color][/size][b][br]<ホモロジー群>[br][/b]鎖の複体では、どの段階でも鎖K=C[sub]n+1[/sub]の境界の境界はゼロ、dn(dn+a(K))=0となる。[br][/size]言い換えると、前回B1がZ1の部分で、Z1がC1の部分だった。[br][br][color=#0000ff][b][size=150]Z1の部分群B1にとよる商群(同値類群)を[br]1次のホモロジー群といい、H1とかく。H1=Z1/B1[/size][/b][/color][br]一般に、n+1複体では、BnがZnの部分でZnがCnの部分となる。[br]n次のホモロジー群はHnとかき、[br][color=#0000ff][b][size=200]Hn=Zn/Bn[br][/size][/b][/color]さらに、Hnを行列にしたときのランク、線形独立の階数をn次のベッチ数βnという。[br]HnがZと同型ならばベッチ数は1[br]HnがZ×Zと同型ならばベッチ数は2[br]Hnがゼロ写像と同型ならばベッチ数は0。[br][br]ホモロジー群がベクトル空間であるとみなしたときの次元数です。[br][b][size=150][size=200]βn=rank(Hn)[/size][/size][/b][br][br]H0​は、点の連結したかたまり数。[br]H1​は、面上の独立な[b]ループ数[/b]。[br]H2​は、面に囲われた閉じた[b]空洞の数[/b]。
1.長方形のホモロジー群H1
[b]<穴のない長方形>[/b][br]穴のない長方形K=ABCDの対角線ACをつけた図形のH1を求めてみよう。[br]K={A,B,C,D,AB,BC,CD,DA,AC,ABC,ACD}[br]2鎖C2(K)=sABC+tACD,[br]1鎖C1(K)=sAB+tBC+uCD+vDA+wAC[br]・輪体Z1=m(AB+BC+CA)+n(AC+CD+DA)=[b]<AB+BC+CA, AC+CD+DA>[/b][br]・境界B1=d2(C2(K))=d2(sABC+tACD)=s(BC-AC+AB)+t(CD-AD+AC)[br] =s(AB+BC+CA)+t(AC+CD+DA)=[b]<AB+BC+CA, AC+CD+DA>[br][/b]・Z1とB1で基底は一致しているから、[b][color=#0000ff]H1=Z1/B1={0}[/color][/b][br][b][br]基点に回収できないループはない。H1の自由度は0次元だから、ベッチ数は0だね。[/b]
穴のない図形ABCD
[b]<穴のある長方形>[/b][br]長方形ABCDに穴ABCをあけた図形のH1を求めてみよう。[br]K={A,B,C,D,AB,BC,CD,DA,AC,ACD}[br]2鎖C2(K)=tACD[br]1鎖C1(K)=sAB+tBC+uCD+vDA+wAC[br]・輪体Z1=m(AB+BC+CA)+n(AC+CD+DA)=<[b]AB+BC+CA, AC+CD+DA[/b]>[br]・境界B1=d2(C2(K))=d2(tACD)=t(CD-AD+AC)=<[b]AC+CD+DA[/b]>[br]・Z1はB1に対して、[b]m(AB+BC+CA)の自由度を多く持っている[/b]。[br] mは整数だから、[b][color=#0000ff]H1=Z1/B1=Z[/color][/b]となる。[br][br][b]三角形ABCに穴のあいたループが1つある。H1の自由度は1次元だから、ベッチ数は1だね。[/b]
穴のある図形ABCD
2.曲面の展開図のホモロジー群
長方形の分析と同じ要領で、[br]ユークリッド3次元空間内の曲面での穴、ループについての情報を得たい。[br]曲面の展開図を使ってみよう。[br][br]計算の量を減らすために、2次元単体としては3角形ではなく2角形もありにしよう。[br]面を細分したり、その逆に統合しても同相であるかどうかは変わらないということを前提とする。[br][br][b][size=150]<球面の展開図>[br][/size][/b]Aが北極、Bが南極とする。2本の経線をC、Dとする。日本の東半球をE、西半球をWとする。[br]球の展開図のH1を求めてみよう。[br]K={A,B,C,D,E,W}[br]2鎖C2(K)=tE+sW[br]1鎖C1(K)=sC+tD[br]・d1(C1(K)=d1(sC+tD)=s(B-A)+t(B-A)=(s+t)(B-A)=0となるのはs=-tとして、[br] 輪体Z1=sC-sD=s(C-D)=<C-D>[br]・境界B1=d2(C2(K))=d2(tE+sW)=t(C-D)+s(C-D)=<[b]C-D[/b]>[br]・Z1とB1で基底は一致しているから、[b][color=#0000ff]H1=Z1/B1={0}[/color][/b][br][b][br]基点に回収できないループはない。H1の自由度は0次元だから、ベッチ数は0だね。[/b]
球の展開図
[size=150][b]<トーラスの展開図>[br][/b][size=100]トーラス面をループvで切ってループuが直線になるように変形すると、[br]トーラス面は円柱の側面の形になるね。Aは2つの円周にあり、uの両端にある。[br]次にループuに沿って切って広げると長方形になる。[br]このとき、Aは2つの円周ループvが直線AAになり、ループuも2直線AAになる。[br]ということは、展開図の長方形はAAAAとなるね。[br][/size][/size]トーラスの展開図AAAAのH1を求めてみよう。[br]K={A,u,v,T}[br]2鎖C2(K)=tT[br]1鎖C1(K)=su+tv[br]・d1(C1(K)=d1(su+tv)=s(A-A)+t(A-A)=(s+t)(A-A)=(s+t)O=0となるのはs,tは任意。[br] 輪体Z1=su+tv=Z<u, v>[br]・境界B1=d2(C2(K))=d2(tT)=t(u-v-u+v)=(t-t)(u-v)={0}[br][b][color=#ff0000] {0}は要素数としては1だから、/B1は/0ではないので、意味がある。[br][/color][/b]・[b][color=#0000ff]H1=Z1/B1=Z<u,v>[/color][/b][br][b][br]1次元のループは2種類あり、H1の自由度はuとvの2次元だから、ベッチ数は2だね。[br][/b]剰余類としては、Z1のu,vのあるなしで2[sup]2[/sup]=4がH1のコセット数になるね。[br][br]このように、見た目で穴数がすぐわかるような図形のホモロジー群を計算するのは[br]それほど意味を感じられないかもしれない。[br]しかし、複雑な曲面や図形でも、根本のロジックは単純だ。[br]それが[b]行列、ベクトルを使った代数計算でできる[/b]というところに可能性を感じられるね。[br][br][color=#9900ff][u][b][size=150]質問:ホモロジー群のベッチ数を求めるコードはどうかけばよいでしょうか。[br][/size][/b][/u][/color][br]行列の階数計算のために、sympyを「pip install sympy」などを実行して環境に入れましょう。[br]・基底ベクトルの名前ではなく、0次から2次の鎖の基底の次元を入れます。[br]・境界写像d3は0->C3への写像だから、いつも0写像です。[br] 境界写像d0はC0->0の写像だから、いつも0写像です。だから、ともにランクは0ですね。[br] 残りの境界写像d1,d2は曲面と用意した点、辺、面の数で変わるので読み取らせましょう。[br]・Znの自由度はC1の自由度からd1、つまりBn-1の自由度を引いた残りです。[br] Bnの自由度は、境界写像dn+1の自由度を入れます。[br] Hn=Zn-Bnから、残った自由度がn次のホモロジー群のベッチ数になりますね。[br][color=#0000ff][In]Python[br]#=============================================[br][/color][color=#0000ff]#複体の点、辺、面の単体次元dimC0,dimC1,dimC2[br]#境界写像d1:C1->C0, d2:C2->C1を引数にして、複体の0次、1次、2次のベッチ数を[br]#ホモロジー群のランクから求める。[br][/color]import numpy as np[br]from sympy import Matrix[br][b][color=#0000ff]def homology_rank(dimC0,dimC1,dimC2,d1,d2):[br] #d3,d0はゼロ写像だから、ランクは0[br] Z0 = dimC0[br] B0 = d1.rank()[br] [br] Z1 = dimC1 - B0[br] B1 = d2.rank()[br][br] Z2 = dimC2 - B1[br] B2 = 0[br] print(f"H0,H1,H2 = (Z0 - B0, Z1 - B1, Z2 - B2 )={Z0 - B0, Z1 - B1, Z2 - B2}")[br][/color][/b][br]# トーラス=====================================================[br][b]# K={A,u,v,T}[br]# 基底の次元 # {A},{u,v},{T}[br][b]# 境界写像d1: C1 -> C0 (d1(u)=d1(v) = A - A = 0だからゼロ写像です)[br][/b]# 境界写像d2: C2 -> C1 (Tの境界は u-v-u+v=0 だからゼロ写像です)[br][/b]homology_rank(1,2,1,Matrix([[0, 0]]),Matrix([[0],[0]]))[br]#==============================================[br][OUT][br][b][color=#0000ff]H0,H1,H2 = (Z0 - B0, Z1 - B1, Z2 - B2 )=(1, 2, 1)[br][/color][/b][br]上の図をみて、球面のホモロジー群のランクを求めてみよう。[br]# 球面=====================================================[br][b]# K={A,B,C,D,E,W}[br]# 基底の次元 # {A,B},{C,D},{E,W}[br][b]# 境界写像d1: C1 -> C0 (d1(C)=d1(D) = B - A)[br][/b]# 境界写像d2: C2 -> C1 (d2(E)=d2(W) = C - D)[br][/b]homology_rank(2,2,2,Matrix([[1, -1],[1, -1]]),Matrix([[1, -1],[1, -1]]))[br]#==============================================[br][OUT][br][b][color=#0000ff]H0,H1,H2 = (Z0 - B0, Z1 - B1, Z2 - B2 )=(1, 0, 1)[/color][/b]
ホモロジー群のランクは前後のBの残差の自由度
3.メビウスの帯など
メビウスの帯の展開図をかいて、[br]ホモロジー群のランクがどうなるかをさぐってみよう。[br]その前に円柱の側面(アニュラス)で練習しておこう。[br]トーラスをループvで切ってのばしたのがアニュラスだ。[br]開くことで、Aの点がAとBに分かれるね。それに伴って辺も1本wが増えるね。[br][br]# アニュラス(円柱の側面)=====================================================[br]# G={A,B,u,v,w,T}[br]# 基底の次元 # {A,B},{u,v,w},{T}[br]# 境界写像d1: C1 -> C0 (d1(u)= B - A ,d1(v) = A - A = 0, d1(w)= B - B = 0)[br]# 行列のu列だけB行とA行で係数1,-1が入り、v列、w列は0にする。[br]# 境界写像d2: C2 -> C1 (Gの境界は u-w-u+v = v-w )[br] # 行列のu行、v行、w行の順にT列の値に係数をかける。[br][b]homology_rank(2,3,1,Matrix([[1, 0, 0],[-1, 0, 0]]),Matrix([[0],[1],[-1]]))[br][/b]#===============================================[br][OUT][br][b]H0,H1,H2 = (Z0 - B0, Z1 - B1, Z2 - B2 )=(1, 1, 0)[br][/b][br]メビウスの帯は同じ辺uどうしを逆向きにして貼り合わせる。[br]# メビウスの帯=====================================================[br]# G={A,B,u,v,w,T}[br]# 基底の次元 # {A,B},{u,v,w},{T}[br]# 境界写像d1: C1 -> C0 (d1(u)= B - A ,d1(v) = A - B , d1(w)= B - A) [br]# 行列のu列、v列、w列の順にB係数行は1,-1,1、A係数行は-1,1,-1[br]# 境界写像d2: C2 -> C1 (Gの境界は u-w+u+v =2u+v-w )[br]# 行列のu行、v行、w行の順にT列の値に係数をかける。[br][b]homology_rank(2,3,1,Matrix([[1, -1, 1],[-1, 1, -1]]),Matrix([[2],[1],[-1]]))[br][/b]#===============================================[br][b]H0,H1,H2 = (Z0 - B0, Z1 - B1, Z2 - B2 )=(1, 1, 0)[br][/b][br]1回ひねっているのに円柱の側面と同じランクになってしまった。[br][br]H0=1なのは、1つながりだから。[br]H1=1なのは、くだにまきつけたループが1種類だけ可能で、あとは0に回収できる。[br]H2=0なのは、曲面の空洞はトーラスや球面のように囲われた空洞がないから。[br][br][b][color=#9900ff]同相な図形はホモロジー群が同型になるという性質が間違っているわけではない。[br]ホモロジー群が同型であることが、同相な図形の条件にはならないということだね。[br][br][/color][/b]面の向きや、ループの向きの情報までは、ホモロジー群でとらえることは難しそうだね。
円柱の側面
メビウスの帯

Information: 曲面のホモロジー群