メビウスの帯と微分方程式

← back メビウスの帯と微分方程式 2026年8月8日 代数的トポロジーは、クラインの壺やメビウスの帯のような、しばしば面白いとされる図形を研究する数学の一分野です。しかし、この紹介の仕方は、現代数学のほぼ全ての分野にとって代数的トポロジーがいかに不可欠であるかを伝えきれていないかもしれません。読者にその重要性を少しでも味わってもらうために、今日の投稿ではまずメビウスの帯の面白い性質を探求し、次に同じ面白い現象が特定の微分方程式を解こうとする際にどのように現れるかを示します。 では、メビウスの帯とは何で、どのような面白い性質があるのでしょうか？ メビウスの帯 紙のシートを用意してください。2つの対辺をテープで貼り合わせると、円筒ができます。代わりに、一方の辺をねじってから貼り合わせると、メビウスの帯が得られます。 メビウスの帯 drag to rotate メビウスの帯は、最初は円筒とそれほど違わないように見えますが、非常に奇妙な性質が1つあります。それは「片面性」です！しかし、それは正確にはどういう意味でしょうか？メビウスの帯の面を指定することは、紙のどちらの面が「上」でどちらが「下」かを決定することと同じです。これらの決定は局所的に一貫している必要があります。つまり、近くの点はどちらが上か下かについて合意する必要があります。言い換えれば、私がメビウスの帯の片面にいるアリだとしたら、数歩歩いても同じ面にいるはずです。これは、任意の曲面に対して2つの面を見つけるのは非常に簡単であることを示唆しています。開始点を選び、一方の面を任意に「上」、もう一方を「下」とラベル付けします。次に、上の意味をコピーしながら移動します。この様子を短いアニメーションで示します。1点での上向きの概念（上向きベクトルで示される）から始め、ベクトルを近くの点に移動させて、一貫した上の概念を与えます。 Local Orientability ▶ Play Press play to animate 局所的には、法線は常に同じ方向を向いて戻ってくる 残念ながら、この手順には小さな問題があります。全体をループしてメビウスの帯を一周する長い旅をすると、出発点に戻ったときに「上」という概念が反転してしまいます！これを観察するための短いアニメーションを作成しました。アニメーションが再生されたら、メビウスの帯をドラッグして回転させ、完全な経路をよりよく見ることができます。 Non-Orientability ▶ Play Press play to animate 1回の完全なループの後、ベクトルが反転する — 同じ表面、反対の向き これが、メビウスの帯に「上」や「下」という概念を定義できない理由です。小さな領域では意味がありますが、メビウスの帯全体を一周すると、「上」と「下」の役割が逆転してしまいます！紙で作られたメビウスの帯があれば、現実世界でこれを観察できます。指でメビウスの帯をなぞってみてください。もし指が帯の「上」から始まったとしたら、それは「下」に終わるでしょう。 微分方程式における問題 メビウスの帯を紹介したので、次に微分の問題に移りましょう。偉大な数学者リーマンは、複素数上で微分方程式を解くことに興味を持っていました。例として、次の方程式を考えてみましょう。 df/dz = (1/2z)f(z) リーマンは、この微分方程式の複素数解を見つけることに興味を持っていました。f(z)のテイラー級数を見つけることによって、このタイプの微分方程式を解くのは簡単です。その方法を説明しましょう。まず、方程式の右辺には2zによる除算が含まれていることに注意してください。0による除算は危険なので、z=1の近くで解を探しましょう。その場合、f(z)のテイラー級数を次のように書くことができます。 f(z) = a0 + a1(z-1) + a2(z-1)^2 + a3(z-1)^3 + ... ここで、a0, a1, a2, ... は解くべき係数です。これは1階微分方程式なので、自由度が1つあることに注意してください。したがって、a0を任意の値を設定できます。以下の計算を簡単にするために、a0 = 1と設定します。つまり、微分方程式にf(1) = 1という初期条件を与えます。 df/dz = (1/2z)f(z) は次と同等です。 2z(df/dz) = f(z) テイラー級数を項ごとに微分すると、次のようになります。 df/dz = a1 + 2a2(z-1) + 3a3(z-1)^2 + ... 上記の式に2zを掛けるのは少しトリッキーです。なぜなら、2zは変数zのべき級数ですが、df/dzはz=1を中心に展開されたべき級数だからです。したがって、乗算を行う前に、次のように書き換えます。 2z = 2 + 2(z-1) したがって、 2z(df/dz) = (2 + 2(z-1)) * (a1 + 2a2(z-1) + 3a3(z-1)^2 + ...) これを展開すると、 2z(df/dz) = 2a1 + (2a1 + 4a2)(z-1) + (4a2 + 6a3)(z-1)^2 + ... したがって、方程式 2z(df/dz) = f(z) は次のように書き換えられます。 2a1 + (2a1 + 4a2)(z-1) + (4a2 + 6a3)(z-1)^2 + ... = 1 + a1(z-1) + a2(z-1)^2 + ... 両辺の(z-1)^nの項の係数を等しくすると、次のようになります。 2a1 = 1 2a1 + 4a2 = a1 4a2 + 6a3 = a2 そして、それ以降も同様です。最初の式2a1 = 1から、a1 = 1/2が得られます。しかし、a1 = 1/2であることがわかったので、この値を2番目の式2a1 + 4a2 = a1に代入して、次のように書き換えることができます。 1 + 4a2 = 1/2 これにより、a2 = -1/8が得られます。同様に、a2 = -1/8であることがわかったので、これを3番目の式4a2 + 6a3 = a2に代入すると、次のようになります。 -4/8 + 6a3 = -1/8 これにより、a3 = 1/16が得られます。 方程式の連鎖をたどるだけで、f(z)のテイラー係数を再帰的に解くことができます。したがって、次の方程式の解f(z)のテイラー係数を好きなだけ求めることができます。 df/dz = (1/2z)f(z) この方法はz=1の解を見つけるのに優れていますが、テイラー級数の欠点は、常に収束するとは限らないことです。この特定のテイラー級数は収束半径が1であることが判明しました。これは、z=1から1単位未満の距離にある任意の複素数zに対して、先ほど見つけた形式的べき級数解が収束することを意味します。ここでは、べき級数のいくつかの値を評価して試すことができます。以下の円盤内の点をクリックすると、そこで級数が評価されます。 z — plane 残念ながら、強調表示された青い円の外側のzの値に対しては、先ほど見つけた形式的べき級数解が常に収束するとは限りません。例として、z = -1を級数に代入すると、次のようになります。 f(-1) = 1 - 1 - 1/2 - 1/2 - 5/8 - 7/8 - 21/16 - ... 最初の1-1の後、項がどんどん負になっていくことに注目してください。これは続き、f(-1) = -∞となることが判明しました。なぜなら、より大きな数を引き続けるからです。これは困ります！複素平面全体で微分方程式を解きたいのですが、今のところ小さな円盤でしか解けません。しかし、このべき級数はz=1を中心に項ごとに展開することによって見つけたことに注意してください。別の中心点周りで項ごとに展開することに理由がないわけではありません。例えば、青い円の端に近い点α = 0.6 + 0.8iを取りましょう。最初のべき級数解を使用して、f(α) ≈ 0.8944 + 0.4472iが見つかります。これを初期条件として使用すると、以前の再帰方法を使用して、αを中心に展開された形式的べき級数解を見つけることができます。 df/dz = (1/2z)f(z) 以下のウィジェットでは、異なる中心点周りに展開された2つの異なる形式的べき級数解を表示します。 expansion around 1 (same one as above) expansion around a = 0.6 + 0.8i z — plane disk around a = 1 disk around a; = 0.6+0.8i ウィジェットでは、2つの異なる形式的べき級数解を、それぞれ異なる点を中心に展開して表示します。fi