微分方程式と力学系

エルゴード理論

力学系の長期的な平均挙動を研究し、測度を保つ変換と結びつける分野。

直観同じ半端な角度でずっと回転し続けたら、やがてすべての場所を訪れるだろうか?

円形の目盛り盤に点を一つ記し、完全な一回転の無理数倍にあたる一定の角度——例えばヒマワリの種や松ぼっくりの鱗片が育つのと同じ、約 $137.5^\circ$ の黄金角——で何度も何度も回転させる。角度が無理数であるため、その点は決して出発点に正確には戻らない。そして驚くべきことに、やがて円周上のあらゆる点に任意に近づき、それぞれの小さな弧を、その弧の長さに比例した頻度で訪れる。ここにはどこにもランダム性は関与していない——規則はただ一つの硬い回転が永遠に繰り返されるだけである——それでも長期的な統計は、まるで点が一様ランダムに選ばれたかのように見える。この均等分布現象、そして決定論的な規則が無限の時間にわたってどのような平均的性質を生み出すかという問いこそが、エルゴード理論の主題である。

thetaスライダーを黄金角

137.5^\circ

だけ繰り返しドラッグしてみよう。印を付けた点は決して正確には繰り返さないが、十分な回数の後には円周上に均等に広がった稠密な点集合を掃き出している——これが無理数回転

T(x)=x+\alpha \pmod 1

の均等分布性とエルゴード性の背後にある幾何学的な描像である。

大学測度を保つ変換

定義: 測度を保つ変換

$(X,\mathcal F,\mu)$ を確率空間(集合 $X$ 、可測部分集合の $\sigma$ -代数 $\mathcal F$ 、 $\mu(X)=1$ を満たす測度 $\mu$ )とする。写像 $T:X\to X$ が測度を保つとは、任意の $A\in\mathcal F$ に対して $\mu(T^{-1}A) = \mu(A)$ が成り立つことをいう—— $A$ に入ることになる点の集合の測度が、 $A$ 自身の測度に等しいということである。直感的には、 $T$ を適用しても確率質量が生成されたり消滅したりすることはなく、ただその位置が並び替えられるだけである。

\mu(T^{-1}A) = \mu(A) \quad \forall\, A \in \mathcal{F}

この理論を支える三つの例がある。 $\alpha$ が無理数である円 $[0,1)$ 上の無理数回転 $T(x)=x+\alpha \pmod 1$ は通常の(ルベーグ)長さを保つ。弧を回転させてもその長さは変わらないからである。倍加写像 $T(x)=2x\bmod 1$ もルベーグ測度を保つ——これはちょうど2対1の写像であり、ある区間の二つの逆像の枝はそれぞれ因子 $2$ で圧縮されるので、それらの合計の長さは元の長さを正確に再構成する。そしてベルヌーイシフト(独立なコイン投げを一歩ずつずらす)は、無限のコイン投げの列の空間上の自然な直積確率測度を保つ。独立な投げの列を一歩ずらしても、依然として独立で同一分布の投げが残るからである。

T^{-1}A=A \;\Longrightarrow\; \mu(A)\in\{0,1\}

測度を保つ $T$ がエルゴード的であるとは、あらゆる $T$ -不変集合が本質的に自明であること、すなわち任意の可測な $A$ に対して $T^{-1}A=A \;\Longrightarrow\; \mu(A)\in\{0,1\}$ が成り立つことをいう。同値な言い方をすれば、系は $T$ が決して混ぜ合わせない正の測度を持つ二つの部分に分割できない——力学系を非自明な形で分解する方法は存在しない。無理数回転も倍加写像もルベーグ測度に関してエルゴード的であるが、その理由はまったく異なり、それを以下の表で正確に示す。

三つの典型例におけるエルゴード性と混合性
変換	エルゴード的か?	混合的か?
無理数回転 $T(x)=x+\alpha \pmod 1$	はい——あらゆる軌道が均等分布する	いいえ——剛体回転は二つの弧を混ぜ合わせることはない
倍加写像 $T(x)=2x\bmod 1$	はい	はい——遠く離れた時刻間の相関はゼロへ減衰する
恒等写像 $T(x)=x$	いいえ( $X$ が一点でない限り)——あらゆる集合が不変である	いいえ

大学バーコフのエルゴード定理とポアンカレの再帰定理

定理: バーコフの各点エルゴード定理

$T$ を確率空間 $(X,\mathcal F,\mu)$ 上の測度を保つ変換とし、 $f\in L^1(\mu)$ とする。このとき時間平均 $\frac{S_nf(x)}{n}$ は $\mu$ -ほとんど全ての $x$ に対して、 $T$ -不変な極限 $\bar f(x)$ に収束し、 $\int_X \bar f\,d\mu = \int_X f\,d\mu$ を満たす。さらに $T$ がエルゴード的であれば、 $\bar f(x)=\int_X f\,d\mu\ \text{a.e.}$ となる——一つの軌道に沿った時間平均が空間平均に等しくなる。

なぜ正しいのか?

これは、物理学者たちが証明なしに何十年も用いてきた直感(統計力学におけるボルツマンのエルゴード仮説)を厳密にしたものである。あるものの長期平均を計算するには、一つの系を非常に長い時間発展させて観察するか、ある瞬間における系全体の集団にわたって平均を取るかのいずれかでよく、エルゴード系についてはこれら二つのまったく異なって聞こえる計算がまさに同じ答えを与える。

証明

ステップ1(不変な limsup と liminf)。 $f^*(x)=\limsup_{n\to\infty}\frac{S_nf(x)}{n}$ と $f_*(x)=\liminf_{n\to\infty}\frac{S_nf(x)}{n}$ を定義する。 $S_{n+1}f(x) = f(x) + S_nf(Tx)$ であるから、 $n+1$ で割って $n\to\infty$ とすると $f^*(Tx)=f^*(x)$ および $f_*(Tx)=f_*(x)$ が示され、両者とも $T$ -不変な関数である。

ステップ2(最大エルゴード定理)。 $\lambda\in\mathbb R$ に対して、 $A_\lambda=\{x: \sup_n S_nf(x)/n>\lambda\}$ を、時間平均の途中経過が $\lambda$ を超えたことがある点の集合とする。鍵となる技術的補題( $g=f-\lambda$ と部分和の最大値 $M_n(x)=\max(0,S_1g(x),\dots,S_ng(x))$ を考え、 $M_n(Tx)\ge S_kg(Tx)$ を各項ごとに用いて $M_n>0$ となる集合上で積分することで証明される)により $\int_{A_\lambda} f\,d\mu \ge \lambda\,\mu(A_\lambda)$ が得られる。

ステップ3(f と f_ を挟み込む)。背理法として、正の測度を持つ集合上で $f^*>f_*$ であると仮定する。すると有理数 $\alpha<\beta$ が存在して、 $E=\{x: f_*(x)<\alpha<\beta<f^*(x)\}$ は正の測度を持つ。 $E$ は( $f^*,f_*$ が不変であるため) $T$ -不変であるから、 $E$ 上だけに注目してよい。 $E$ 上で $f-\beta$ に最大不等式を適用すると $\int_E f\,d\mu\ge\beta\mu(E)$ が強制され、 $\alpha-f$ に同様に適用すると $\int_E f\,d\mu\le\alpha\mu(E)$ が強制される。 $\alpha<\beta$ かつ $\mu(E)>0$ であるからこれらは矛盾する。したがって $f^*=f_*$ がほとんど至るところで成り立ち、極限 $\bar f(x)=\lim_n S_nf(x)/n$ はほとんど至るところで存在し $T$ -不変である。

ステップ4(積分を一致させる、およびエルゴードの場合)。優収束論法( $f$ を切り詰め、最大不等式を再び用いて裾を制御する)により $\int_X \bar f\,d\mu = \int_X f\,d\mu$ が示される。最後に、 $T$ がエルゴード的であれば、 $T$ -不変な関数 $\bar f$ はほとんど至るところで定数でなければならない(そのレベル集合 $\{\bar f\le c\}$ のそれぞれが $T$ -不変であり、したがって測度が $0$ か $1$ であるという、エルゴード性の定義そのものによる)。積分が一致することと合わせると、その定数は $\int_X f\,d\mu$ でなければならず、 $\bar f(x)=\int_X f\,d\mu\ \text{a.e.}$ が得られる。

定理: ポアンカレの再帰定理

$T$ を確率空間(より一般には有限測度空間) $(X,\mathcal F,\mu)$ 上の測度を保つ変換とし、 $A\in\mathcal F$ が $\mu(A)>0$ を満たすとする。このとき $A$ のほとんど全ての点は無限回 $A$ に戻る:ほとんど全ての $x\in A$ について、無限個の $n\ge1$ に対して $T^nx\in A$ が成り立つ。

なぜ正しいのか?

空間が有限であり、何も破壊されない(測度が保たれる)ならば、ある領域が永遠に、まったく新しい、これまで訪れたことのない領域へと点を送り続けることはできない——やがて系は、単に戻ってくる測度を置く新しい場所がもう残っていないという理由だけで、すでに訪れた場所を再訪し始めざるを得なくなる。

証明

ステップ1(決して戻らない点)。 $A$ の点のうち決して $A$ に戻らないものの集合として $A_0=\{x\in A: T^nx\notin A\ \forall n\ge1\}$ とおく。集合 $A_0, T^{-1}A_0, T^{-2}A_0,\dots$ は互いに素である。 $x\in T^{-i}A_0\cap T^{-j}A_0$ ( $i<j$ )とすると、 $T^ix\in A_0$ であるが同時に $j-i\ge1$ に対し $T^j x = T^{j-i}(T^ix) \in A$ となり、 $T^ix\in A_0$ が決して $A$ に戻らないことに矛盾する。したがって $T^{-i}A_0 \cap T^{-j}A_0=\varnothing\ (i\ne j)$ 。

ステップ2(決して戻らない集合の測度はゼロ)。 $T$ は測度を保つので、任意の $i$ について $\mu(A_0)=\mu(T^{-i}A_0)$ 。もし $\mu(A_0)>0$ ならば、可算個の互いに素な集合 $T^{-i}A_0$ ( $i=0,1,2,\dots$ )はすべてこの同じ正の測度を持つことになり、それらの和集合は無限の測度を持つことになる—— $\mu(X)<\infty$ (実際には $\mu(X)=1$ )であるから不可能である。したがって $\mu(A_0)=0$ 。

ステップ3(有限回しか戻らない点の集合も測度ゼロ)。 $B=\{x\in A: x\ \text{returns to}\ A\ \text{only finitely often}\}$ とおく。 $B$ を、シフトされた力学系のもとでの $T^kA$ の(本質的には)決して戻らない集合にわたる $k$ についての可算和として $B=\bigcup_{k\ge0} T^{-k}\{x\in T^kA: x\ \text{never returns to}\ T^kA\}$ と書けば、各項は( $\mu(T^kA)=\mu(A)$ を用いて) $A$ の代わりに $T^kA$ に対してステップ1–2の議論をそのまま適用することで測度ゼロとなる。可算劣加法性により $\mu(B)=0$ である。

ステップ4(結論)。 $\mu(B)=0$ であるから $A$ において全測度を持つ任意の $x\in A\setminus B$ は、 $B$ の定義により無限回 $A$ に戻る。これがまさに定理の主張である。

発展混合:エルゴード性より強い性質

定義: 強混合

測度を保つ $T$ が(強)混合的であるとは、任意の可測な $A,B$ に対して $\lim_{n\to\infty}\mu(T^{-n}A\cap B)=\mu(A)\mu(B)$ が成り立つことをいう: $n$ ステップ後に $B$ のうち $A$ に戻ってくる割合が、 $A$ と $B$ が統計的に独立であった場合の値に収束する。混合はエルゴード性を含意するが( $T^{-1}A=A$ を満たす $B=A$ を取ると $\mu(A)=\mu(A\cap A)\to\mu(A)^2$ となり $\mu(A)\in\{0,1\}$ が強制される)、逆は成り立たない。

無理数回転はなぜエルゴード的でありながら混合的ではないのか。 $A=B$ を小さな弧とする:それを $n\alpha$ だけ回転させても、同じ小さな弧が円周上を剛体的に移動するだけなので、 $T^{-n}A\cap A$ は空であるか、あるいは(回転が均等分布するため無限個の $n$ について)再び $A$ のほとんど全体に非常に近くなるかのいずれかである——重なりの割合は独立な場合の値 $\mu(A)^2$ に決して落ち着かず、 $\mu(A)$ 付近まで何度も揺り戻される。対照的に倍加写像は、あらゆる小さな区間を空間全体にわたって指数関数的な速さで引き伸ばし折りたたみ、真に撹拌する——これが混合性の原因であり、最終的には $r=4$ のロジスティック力学のようなカオス的写像を統計的な目的で「ランダムと同じくらい良い」ものとして扱うことにつながる。

大学実世界での応用と具体例

エルゴード理論は「無限時間にわたる平均的挙動」を計算可能で証明可能な主張に変える。これはまさに統計力学が、単一の物理系の長期平均を集団平均で置き換えることを正当化するために必要とするもの(ボルツマンのエルゴード仮説)であり、データ圧縮が、ある情報源をどれだけ圧縮できるかを制限するエントロピー率を定義するために必要とするもの(対応するシフトのコルモゴロフ・シナイ・エントロピーを通じて)であり、GoogleのPageRankのような検索エンジンが、ランダムにウェブを閲覧する者に対して一意な長期訪問頻度を保証するために必要とするもの(マルコフ連鎖の不変測度)であり、そしてカオス的写像から構築される暗号論的擬似乱数生成器が、その出力を統計的にランダムとみなすことを正当化するために必要とするものである。いずれの場合も、エルゴード性と混合性こそが、単一の決定論的軌道をあたかも系全体についての統計的情報を担っているかのように扱うことを許すのである。

例: 無理数回転は与えられた弧をどのくらいの頻度で訪れるか?

$\alpha$ を無理数として、円 $[0,1)$ 上にルベーグ測度とともに作用する $T(x)=x+\alpha\pmod 1$ を考え、 $A=[0,0.3)$ とする。典型的な出発点 $x_0$ について、最初の $n$ 個の反復 $x_0,x_1,\dots,x_{n-1}$ のうち $A$ に入るものの割合は、 $n\to\infty$ のときどうなるか?

解答

これはまさに $f=\mathbf 1_A$ ( $A$ の指示関数)によるバーコフのエルゴード定理の計算である。時間平均 $\frac1n\sum_{k=0}^{n-1}\mathbf 1_A(x_k)$ は、まさに最初の $n$ 個の反復のうち $A$ に入るものの割合である。

無理数回転はルベーグ測度に関するエルゴード変換の古典的な例である(これは実質的にワイルの均等分布定理である。任意の $T$ -不変集合をフーリエ級数に展開すると、回転がそれらを $e^{2\pi i k\alpha}\ne1$ 倍することにより、ゼロでないフーリエ係数はすべて消えざるを得ず、定数関数だけが残る)。

$T$ がエルゴード的であるから、バーコフの定理は最も強い形で適用される。すなわち、時間平均は、出発点にわたる平均だけでなく、あらゆる典型的な出発点について空間平均に等しい。空間平均は $\int_X \mathbf 1_A\,d\mu = \mu(A) = 0.3-0 = 0.3$ である。

したがって、軌道が $A$ で過ごす時間の長期的な割合は、どの $x_0$ から出発しても(測度ゼロの例外集合を除いて)正確に $0.3$ である——「決定論的な回転は、統計的には、毎回 $A$ 内の一様ランダムな点を選ぶのとまったく同じように振る舞う」ということの厳密な言い換えである。

例: 偏りのあるデータ源はどこまで圧縮できるか?ベルヌーイシフトのコルモゴロフ・シナイ・エントロピー

あるデータ源が独立なビットを出力し、それぞれ確率 $p=0.3$ で $1$ 、確率 $0.7$ で $0$ であるとする——これは積測度 $(0.3,0.7)$ を持つ列の上のベルヌーイシフトとしてエルゴード的にモデル化される。そのエントロピー率(シフトのコルモゴロフ・シナイ・エントロピーであり、ここでは単一の記号のシャノンエントロピーに等しい)を計算せよ。これはシャノンの情報源符号化定理により、この情報源を可逆圧縮するために平均して必要な記号あたりの最小ビット数である。

解答

独立同分布(ベルヌーイ)情報源では、シフト写像のコルモゴロフ・シナイ・エントロピーは単一の記号の通常のシャノンエントロピーに帰着する: $h=-\sum_i p_i\log_2 p_i$ 。

$p_1=0.3$ 、 $p_0=0.7$ を代入すると: $h=-0.3\log_2 0.3-0.7\log_2 0.7$ 。

各項を計算する: $-0.3\log_2 0.3 = 0.3\times1.737=0.521$ ビット、 $-0.7\log_2 0.7=0.7\times0.515=0.361$ ビット( $\log_2 0.3\approx-1.737$ 、 $\log_2 0.7\approx-0.515$ を用いた)。

これらを足し合わせると $h\approx0.881\ \text{bits}$ となる。したがって平均して、どのような可逆符号もこの情報源を記号あたり約 $0.881$ ビットより小さく圧縮することはできない(公平なコインに必要な記号あたり $1$ ビットよりも明らかに少ない。これはまさに偏りによって情報源がより予測しやすくなり、したがってより圧縮しやすくなるためである)——そしてシャノンの定理は、この限界が達成可能であることを保証する。

よくある誤り. エルゴード性と混合性はしばしば混同されるが、実際には異なる概念である。エルゴード性は、系が二つの不変な部分に分割できないこと(同値に、時間平均が空間平均に等しいこと)しか述べないが、混合性はさらに、事象間の相関が時間とともにゼロへ減衰することを述べる。あらゆる混合系はエルゴード的であるが、上述の無理数回転は、エルゴード性が混合性を含意しないことを示す標準的な反例である。もう一つのよくある誤りは、ポアンカレの再帰定理を過剰に解釈することである。それは

A

のほとんど全ての点が無限回

A

に戻ることを保証するだけで、いつ戻るかについては何も述べていない——期待される再帰時間は天文学的に長くなりうる(カックの補題により、

A

への平均再帰時間は

1/\mu(A)

であり、小さな集合

A

に対しては非常に大きい)。またこの定理はまったく同じ点に戻ることについては何の情報も与えず、集合

A

に戻ることのみを保証する。

歴史的ノート

この分野の動機となった問いは19世紀の統計力学に由来する。ルートヴィヒ・ボルツマンのエルゴード仮説(1870年代)は、気体の相空間内の軌道が、そのエネルギーと矛盾しないあらゆる状態に最終的に任意に近づくと提唱し、時間平均をはるかに計算しやすい集団平均で置き換えることを正当化した。アンリ・ポアンカレは1890年、三体問題に関する懸賞論文の中の技術的な一段階として再帰定理を証明し、有界な力学系が最終的に以前の配置に任意に近く戻らなければならないことを示した——これは当初ほとんど逆説的だと見なされた結果であった(系が戻ってくるならエントロピーはどうして増加しうるのか?)。これは、天文学的に長い再帰時間が正しく理解されて初めて解消された。厳密なエルゴード理論という数学分野は1931年から1932年にかけて始まった。ジョージ・バーコフとジョン・フォン・ノイマンは互いに数か月のうちに独立に研究し、それぞれエルゴード定理の一つの版を証明した——フォン・ノイマンの $L^2$ 平均版が先に得られたがやや遅れて発表され、一方バーコフのより強い各点(ほとんど至るところ)版は1931年12月に米国科学アカデミー紀要に掲載され、提唱されてから半世紀以上を経て、ついにボルツマンの仮説を厳密な数学的基盤の上に置いた。

アンリ・ポアンカレジョン・フォン・ノイマン

研究の最前線 2026年時点

エルゴード理論は、その起源である統計力学をはるかに超えた不可欠な道具となっている。ピーター・サルナックのメビウス非相関予想(2010年頃に提唱)は、数論のメビウス関数が、位相的エントロピーがゼロである任意の決定論的系の軌道と漸近的に無相関であることを予測する——「素数は十分にカオス的でないあらゆる時計にとってランダムに見える」ということの正確なエルゴード理論的定式化であり、多くの特定のクラスの系について検証されているが、2026年現在一般には未解決である。算術級数に関するセメレディの定理のヒレル・ファーステンバーグによるエルゴード理論的な再証明(1977年)は、エルゴード・ラムゼー理論という分野を切り開き、その多重再帰定理は、より一般的な群作用や配置へとさらに拡張され続けている。より応用的な側面では、分子動力学シミュレーションから、ノイズを伴う力学系とみなされたニューラルネットワークの学習ダイナミクスに至るまで、現実的な高次元系についてエントロピーと混合率を厳密に計算あるいは評価することは大部分未解決のままであり、量子多体系がどれほど速く「熱平衡化」する(固有状態熱化仮説によって定式化される、エルゴード的な挙動に近づく)かを理解することは、エルゴード理論を物性物理学や量子情報とつなぐ活発な研究分野である。

変換 $T:X\to X$ が確率測度 $\mu$ を保つ、すなわち $\mu(T^{-1}A) = \mu(A)$ であるとする。この等式はどのような集合 $A$ について成り立たなければならないか?

バーコフのエルゴード定理により、エルゴード的な無理数回転 $T(x)=x+\alpha \pmod 1$ と $A=[0,0.3)$ について、典型的な軌道が $A$ で過ごす長期的な時間の割合はどれだけか?

ポアンカレの再帰定理によれば、有限測度空間上の測度を保つ系で $\mu(A)>0$ のとき、 $A$ のほとんど全ての点は……

あるソースが $P(1)=0.3$ の独立同分布ビットを発する。そのエントロピー率(対応するベルヌーイシフトのコルモゴロフ・シナイ・エントロピー)は、ビット単位で小数点以下2桁に丸めると、どれに最も近いか?

参考文献

Peter Walters (1982). An Introduction to Ergodic Theory
George D. Birkhoff (1931). Proof of the Ergodic Theorem
John von Neumann (1932). Proof of the Quasi-Ergodic Hypothesis
Hillel Furstenberg (1981). Recurrence in Ergodic Theory and Combinatorial Number Theory