コンテンツにスキップ

エルゴード定理:測度保存系・混合・von Neumann と Birkhoff の定理

対象読者: 本テーマの条件付き期待値とコルモゴロフの拡張定理の章、およびヒルベルト空間の射影定理を学んだ人。ルベーグ積分とフーリエ級数の初歩は既知とし、力学系の知識は仮定しません。

前提となる記事:条件付き期待値:σ-加法族で条件づけ、ラドン・ニコディムで存在を保証する、Hilbert 空間の基礎:射影定理・リースの表現定理・正規直交基底、Kolmogorov の拡張定理と連続性定理:有限次元分布から Wiener 測度へ

生 Markdown

ある規則で点を動かし続けたとき、1 つの点の動きを長く平均したものと、空間全体の平均とは、いつ同じになるでしょうか。この記事は、その問いに答える 2 つの定理(時間平均が収束する、その極限が空間平均になる)を、例を確かめながら証明します。読み終えると、サイコロを何度も振った平均が期待値に近づく大数の強法則が、この一般論の特別な場合として出てくることが分かります。以下の要点は全体の見取り図なので、用語が分からなくても先へ進んでかまいません。

  • 確率空間 (X,B,μ)(X,\mathcal{B},\mu) 上の可測写像 TT が μ(T−1A)=μ(A)\mu(T^{-1}A)=\mu(A) を満たすとき測度保存系といいます。TT で動かしても消えない集合が自明なものだけのときエルゴード的といい、これは「TT 不変な関数は定数だけ」と同値です。
  • 円の無理数回転はエルゴード的ですが混合的ではありません。2 倍写像と Bernoulli シフトは混合的です。どれも Fourier 展開か「柱集合による近似」(座標を有限個だけ指定して決まる集合で、全体を近似すること)で示せます。
  • von Neumann の平均エルゴード定理: Hilbert 空間の縮小作用素(ノルムを増やさない線形作用素)VV について、1n∑k<nVkx\frac1n\sum_{k<n}V^kx は不動点空間への直交射影 PxPx にノルム収束します。証明は H=ker⁡(I−V)⊕Ran(I−V)‾\mathcal{H}=\ker(I-V)\oplus\overline{\mathrm{Ran}(I-V)} という分解だけです。
  • Birkhoff の各点エルゴード定理: f∈L1(μ)f\in L^1(\mu) の時間平均 1n∑k<nf(Tkx)\frac1n\sum_{k<n}f(T^kx) はほとんどすべての xx で収束し、極限は不変 σ-加法族への条件付き期待値 E[f∣I]E[f\mid\mathcal{I}] です。エルゴード的なら極限は空間平均 ∫f dμ\int f\,d\mu です。鍵は最大エルゴード補題(Garsia の 3 行の証明)です。
  • 応用として、独立同分布の確率変数列の大数の強法則が Bernoulli シフトのエルゴード性から従います。次章では、この定理を劣加法的な列に拡張して Lyapunov 指数を扱います。

気体の分子を 1 つ追いかけ、ある量(たとえば運動エネルギー)を長い時間にわたって平均したもの(時間平均)と、ある瞬間にすべての分子について平均したもの(空間平均、相平均)は一致するでしょうか。Boltzmann は統計力学の基礎づけにこの一致を仮定しました(エルゴード仮説)。これを数学の問題にしたのが 1931〜32 年の von Neumann と Birkhoff の定理で、「軌道がほとんどすべての点で空間全体をまんべんなく訪れる」ための条件としてエルゴード性を取り出しました。

確率論から見ると、エルゴード定理は大数の強法則の一般化です。独立同分布の列 X1,X2,…X_1,X_2,\dots は、列全体の空間 RN\mathbb{R}^{\mathbb{N}} の上で「1 つずらす」写像(シフト)の軌道とみなせます。独立性は使わず、ずらしても分布が変わらない(定常性)ことと、不変な事象が自明であること(エルゴード性)だけで、標本平均の概収束が言えます。

この章の結果は、準周期的な Schrödinger 作用素(Harper 模型とほぼ Mathieu 作用素)で、位相 θ\theta を θ+α\theta+\alpha に動かす無理数回転の上の「転送行列の積」を調べるときの土台になります。

以下 (X,B,μ)(X,\mathcal{B},\mu) は確率空間(μ(X)=1\mu(X)=1)です。

定義 2.1(測度保存系)↗

可測写像 T:X→XT:X\to X がすべての A∈BA\in\mathcal{B} について μ(T−1A)=μ(A)\mu(T^{-1}A)=\mu(A) を満たすとき、TT は μ\mu を保存するといい、(X,B,μ,T)(X,\mathcal{B},\mu,T) を測度保存系という。T0=idT^0=\mathrm{id}、Tn=T∘Tn−1T^{n}=T\circ T^{n-1} とし、可測関数 ff について

Snf=∑k=0n−1f∘Tk,Anf=1nSnf(n≥1),S0f=0S_nf=\sum_{k=0}^{n-1}f\circ T^k,\qquad A_nf=\frac1nS_nf\quad(n\ge1),\qquad S_0f=0

を Birkhoff 和、Birkhoff 平均(時間平均)という。

TT が可逆である必要はありません。条件は逆像 T−1AT^{-1}A について書かれていて、これは「TT で 1 歩進んだ後に AA にいる点の集合」です。

補題 2.2(積分の不変性と Koopman 作用素)✓

(X,B,μ,T)(X,\mathcal{B},\mu,T) を測度保存系とする。

  1. 可測関数 f:X→[0,∞]f:X\to[0,\infty] について ∫f∘T dμ=∫f dμ\int f\circ T\,d\mu=\int f\,d\mu。
  2. f∈L1(μ)f\in L^1(\mu) なら f∘T∈L1(μ)f\circ T\in L^1(\mu) で ∫f∘T dμ=∫f dμ\int f\circ T\,d\mu=\int f\,d\mu。
  3. 1≤p<∞1\le p<\infty について UTf=f∘TU_Tf=f\circ T は Lp(μ)L^p(\mu) 上の well-defined な線形等長作用素である。特に ∥Anf∥p≤∥f∥p\|A_nf\|_p\le\|f\|_p。
証明(補題 2.2)

1。f=1Af=\mathbf{1}_A なら 1A∘T=1T−1A\mathbf{1}_A\circ T=\mathbf{1}_{T^{-1}A} で、両辺は μ(T−1A)=μ(A)\mu(T^{-1}A)=\mu(A) です。非負単関数については線形性から従います。一般の ff には単関数近似定理(定理 3.7[ルベーグ積分の定義と収束定理])で 0≤sk↑f0\le s_k\uparrow f を取ると sk∘T↑f∘Ts_k\circ T\uparrow f\circ T なので、単調収束定理(定理 5.2[ルベーグ積分の定義と収束定理])で両辺の極限を取れます。

2。f=f+−f−f=f^+-f^- に分けて 1 を使います。

3。f=gf=g a.e. なら {f∘T≠g∘T}=T−1{f≠g}\{f\circ T\ne g\circ T\}=T^{-1}\{f\ne g\} は測度 00 なので、UTU_T は同値類の上で定まります。1 を ∣f∣p|f|^p に使うと ∥f∘T∥p=∥f∥p\|f\circ T\|_p=\|f\|_p です。AnfA_nf については Minkowski の不等式(定理 4.5[L^p 空間と関数解析への導入])から ∥Anf∥p≤1n∑k<n∥f∘Tk∥p=∥f∥p\|A_nf\|_p\le\frac1n\sum_{k<n}\|f\circ T^k\|_p=\|f\|_p です。

UTU_T を Koopman 作用素といいます。非線形な写像 TT を、関数空間の上の線形作用素に置き換えるのがエルゴード理論の基本的な発想です。

例 2.3(円の回転)

T=R/Z\mathbb{T}=\mathbb{R}/\mathbb{Z} を [0,1)[0,1) と同一視し、Lebesgue 測度 mm を入れます。α∈R\alpha\in\mathbb{R} について Rαx=x+α(mod1)R_\alpha x=x+\alpha\pmod 1 とすると、Rα−1AR_\alpha^{-1}A は AA を −α-\alpha だけ平行移動して [0,1)[0,1) に折り返した集合で、Lebesgue 測度の平行移動不変性から m(Rα−1A)=m(A)m(R_\alpha^{-1}A)=m(A) です。

例 2.4(2 倍写像)

Dx=2x(mod1)D x=2x\pmod1 を [0,1)[0,1) 上で考えます。区間について D−1[a,b)=[a2,b2)∪[a+12,b+12)D^{-1}[a,b)=[\tfrac a2,\tfrac b2)\cup[\tfrac{a+1}2,\tfrac{b+1}2) なので m(D−1[a,b))=b−am(D^{-1}[a,b))=b-a です。区間 [a,b)[a,b)(と空集合)は π-系で Borel 集合族を生成するので、測度の一意性(定理 4.5[測度空間と Carathéodory の拡張定理])から m∘D−1=mm\circ D^{-1}=m です。DD は 2 対 1 で可逆ではありません。xx の 2 進展開 x=0.x1x2x3…x=0.x_1x_2x_3\dots に対して Dx=0.x2x3…Dx=0.x_2x_3\dots なので、DD は「2 進展開の桁を 1 つずらす」写像です。

次の補題は、コインを無限回投げる試行のように、無限個の座標をもつ空間に測度を入れる話です。座標を最初の mm 個だけ指定した集合(柱集合)の測度を ν⊗m\nu^{\otimes m} で決めておき、それを全体へ延ばすのが Kolmogorov の拡張定理の役目です。

補題 2.5(Bernoulli シフトの構成)✓

ν\nu を R\mathbb{R} 上の確率測度、Ω=RN\Omega=\mathbb{R}^{\mathbb{N}}、F\mathcal{F} を座標写像 x↦xkx\mapsto x_k で生成される直積 σ-加法族とする。F\mathcal{F} 上の確率測度 PP で、任意の mm と Borel 集合 B⊂RmB\subset\mathbb{R}^m について P({(x1,…,xm)∈B})=ν⊗m(B)P(\{(x_1,\dots,x_m)\in B\})=\nu^{\otimes m}(B) を満たすもの P=ν⊗NP=\nu^{\otimes\mathbb{N}} がただ 1 つ存在する。さらにシフト σ(x)k=xk+1\sigma(x)_k=x_{k+1} は可測で P∘σ−1=PP\circ\sigma^{-1}=P を満たす。

証明(補題 2.5)

有限集合 J⊂NJ\subset\mathbb{N} ごとの積測度 μJ=ν⊗J\mu_J=\nu^{\otimes J}(定理 3.2[直積測度と Fubini–Tonelli の定理])は、J⊂J′J\subset J' のとき μJ′\mu_{J'} の RJ\mathbb{R}^J への周辺分布が μJ\mu_J になる(余分な座標で全空間 R\mathbb{R} を取れば ν(R)=1\nu(\mathbb{R})=1 が掛かるだけ)ので整合的です。Kolmogorov の拡張定理(定理 4.1[Kolmogorov の拡張定理と連続性定理])から、柱集合の上で μJ\mu_J と一致する確率測度 PP が F\mathcal{F} 上にただ 1 つあります。

σ\sigma の可測性は、xk∘σ=xk+1x_k\circ\sigma=x_{k+1} が各 kk で可測であることから従います。柱集合 C={(x1,…,xm)∈B}C=\{(x_1,\dots,x_m)\in B\} について σ−1C={(x2,…,xm+1)∈B}\sigma^{-1}C=\{(x_2,\dots,x_{m+1})\in B\} で、PP の下での座標 (x2,…,xm+1)(x_2,\dots,x_{m+1}) の分布は μ{2,…,m+1}=ν⊗m\mu_{\{2,\dots,m+1\}}=\nu^{\otimes m} なので P(σ−1C)=ν⊗m(B)=P(C)P(\sigma^{-1}C)=\nu^{\otimes m}(B)=P(C) です。柱集合の全体は代数(特に π-系)で F\mathcal{F} を生成するので、測度の一意性(定理 4.5[測度空間と Carathéodory の拡張定理])から P∘σ−1=PP\circ\sigma^{-1}=P です。

例 2.6(Bernoulli シフト)

補題 2.5 の (RN,F,ν⊗N,σ)(\mathbb{R}^{\mathbb{N}},\mathcal{F},\nu^{\otimes\mathbb{N}},\sigma) を Bernoulli シフトとよびます。座標 x1,x2,…x_1,x_2,\dots は PP の下で独立で、各々の分布は ν\nu です。シフトは「1 回分の試行を忘れて次から数え直す」操作で、測度保存系の最も基本的な例です。ν\nu が 2 点 {0,1}\{0,1\} に重み 12\frac12 ずつ乗る場合は、2 倍写像と(零集合を除いて)同じものです。

測度保存性だけから、すでに次の定性的な結論が出ます。

定理 2.7(Poincaré の回帰定理)✓

(X,B,μ,T)(X,\mathcal{B},\mu,T) を測度保存系、A∈BA\in\mathcal{B} とする。ほとんどすべての x∈Ax\in A について、Tnx∈AT^nx\in A となる n≥1n\ge1 が無限個ある。

証明(定理 2.7)

B=A∖⋃n≥1T−nAB=A\setminus\bigcup_{n\ge1}T^{-n}A(AA から出て二度と戻らない点)とおきます。0≤n<m0\le n<m で x∈T−nB∩T−mBx\in T^{-n}B\cap T^{-m}B なら、y=Tnx∈By=T^nx\in B かつ Tm−ny∈B⊂AT^{m-n}y\in B\subset A で、y∈By\in B に反します。よって B,T−1B,T−2B,…B,T^{-1}B,T^{-2}B,\dots は互いに素で、どれも測度 μ(B)\mu(B) です。和の測度は 11 以下なので μ(B)=0\mu(B)=0 です。

x∈Ax\in A が有限回しか戻らないとし、最後に AA にいた時刻を m≥0m\ge0 とすると、Tmx∈BT^mx\in B です。したがってそのような xx は ⋃m≥0T−mB\bigcup_{m\ge0}T^{-m}B に入り、これは零集合の可算和で零集合です。

有限測度が本質的です。R\mathbb{R} 上の平行移動 x↦x+1x\mapsto x+1 は Lebesgue 測度を保ちますが、[0,1)[0,1) の点は二度と戻りません。

定義 3.1(不変集合とエルゴード性)↗

(X,B,μ,T)(X,\mathcal{B},\mu,T) を測度保存系とする。T−1A=AT^{-1}A=A を満たす A∈BA\in\mathcal{B} を不変集合といい、その全体を I\mathcal{I} と書く(逆像は集合演算と可換なので I\mathcal{I} は σ-加法族である)。すべての A∈IA\in\mathcal{I} が μ(A)∈{0,1}\mu(A)\in\{0,1\} を満たすとき、TT(または μ\mu)はエルゴード的であるという。

不変集合が測度 00 と 11 の中間にあれば、XX は互いに行き来しない 2 つの部分に分かれ、その両方で別々に時間平均を取ることになります。エルゴード性はそれが起きないという条件です。実際に使うときは「ほとんど不変」な集合や関数で十分なことを、次の補題で保証しておきます。

この補題のねらいは、測度 00 の違いを無視してよい「ほとんど不変」な集合を、本当に不変な集合に取り替えられると示すことです。そうすると、エルゴード性を確かめるときに、厳密な不変性まで要求しなくてすみます。

補題 3.2(ほとんど不変な集合)✓

A∈BA\in\mathcal{B} が μ(T−1A △ A)=0\mu(T^{-1}A\,\triangle\,A)=0 を満たすなら、A∞=⋂N≥1⋃n≥NT−nAA_\infty=\bigcap_{N\ge1}\bigcup_{n\ge N}T^{-n}A は不変集合で、μ(A △ A∞)=0\mu(A\,\triangle\,A_\infty)=0 である。

証明(補題 3.2)

T−kT^{-k} は対称差と可換なので、測度保存性から μ(T−(k+1)A △ T−kA)=μ(T−k(T−1A △ A))=0\mu(T^{-(k+1)}A\,\triangle\,T^{-k}A)=\mu(T^{-k}(T^{-1}A\,\triangle\,A))=0 です。三角不等式 C △ E⊂(C △ D)∪(D △ E)C\,\triangle\,E\subset(C\,\triangle\,D)\cup(D\,\triangle\,E) を繰り返して μ(T−nA △ A)≤∑k<nμ(T−(k+1)A △ T−kA)=0\mu(T^{-n}A\,\triangle\,A)\le\sum_{k<n}\mu(T^{-(k+1)}A\,\triangle\,T^{-k}A)=0 です。x∈A∞∖Ax\in A_\infty\setminus A なら x∈T−nA∖Ax\in T^{-n}A\setminus A となる nn があり、x∈A∖A∞x\in A\setminus A_\infty なら十分大きいすべての nn で x∈A∖T−nAx\in A\setminus T^{-n}A です。どちらも ⋃n(T−nA △ A)\bigcup_n(T^{-n}A\,\triangle\,A) に入るので、μ(A △ A∞)=0\mu(A\,\triangle\,A_\infty)=0 です。最後に T−1A∞=⋂N⋃n≥NT−(n+1)A=A∞T^{-1}A_\infty=\bigcap_N\bigcup_{n\ge N}T^{-(n+1)}A=A_\infty です。

命題 3.3(エルゴード性の言い換え)✓

測度保存系について次は同値である。

  1. TT はエルゴード的である。
  2. μ(T−1A △ A)=0\mu(T^{-1}A\,\triangle\,A)=0 なる A∈BA\in\mathcal{B} はすべて μ(A)∈{0,1}\mu(A)\in\{0,1\} を満たす。
  3. 可測関数 f:X→Rf:X\to\mathbb{R} が f∘T=ff\circ T=f a.e. を満たせば、ff は a.e. で定数である。
  4. f∈L2(μ)f\in L^2(\mu)(複素数値)が UTf=fU_Tf=f を満たせば、ff は a.e. で定数である。
証明(命題 3.3)

1 ⇒ 2:補題 3.2 の A∞A_\infty は不変集合なので μ(A)=μ(A∞)∈{0,1}\mu(A)=\mu(A_\infty)\in\{0,1\} です。

2 ⇒ 3:c∈Rc\in\mathbb{R} について T−1{f>c}={f∘T>c}T^{-1}\{f>c\}=\{f\circ T>c\} は {f>c}\{f>c\} と零集合しか違わないので、μ(f>c)∈{0,1}\mu(f>c)\in\{0,1\} です。φ(c)=μ(f>c)\varphi(c)=\mu(f>c) は単調減少で、ff が実数値なので測度の連続性から c→−∞c\to-\infty で 11、c→+∞c\to+\infty で 00 に近づきます。c0=sup⁡{c:φ(c)=1}c_0=\sup\{c:\varphi(c)=1\} は有限で、各 kk について μ(f>c0+1k)=0\mu(f>c_0+\frac1k)=0、μ(f>c0−1k)=1\mu(f>c_0-\frac1k)=1 です。kk について可算個の和・共通部分を取って μ(f≠c0)=0\mu(f\ne c_0)=0 です。

3 ⇒ 4:UTf=fU_Tf=f は f∘T=ff\circ T=f a.e. のことで、実部と虚部が別々に同じ等式を満たします。

4 ⇒ 1:A∈IA\in\mathcal{I} なら 1A∘T=1T−1A=1A\mathbf{1}_A\circ T=\mathbf{1}_{T^{-1}A}=\mathbf{1}_A なので 1A\mathbf{1}_A は a.e. 定数、すなわち μ(A)∈{0,1}\mu(A)\in\{0,1\} です。

定義 3.4(混合)↗

測度保存系が混合的(強混合)であるとは、すべての A,B∈BA,B\in\mathcal{B} について

lim⁡n→∞μ(T−nA∩B)=μ(A)μ(B)\lim_{n\to\infty}\mu(T^{-n}A\cap B)=\mu(A)\mu(B)

が成り立つことをいう。

T−nAT^{-n}A は「nn 歩後に AA にいる」事象で、混合性は「十分時間がたつと、現在の事象 BB と未来の事象が独立になる」ことを意味します。

命題 3.5(混合ならエルゴード的)✓

混合的な測度保存系はエルゴード的である。

証明(命題 3.5)

A∈IA\in\mathcal{I} なら T−nA=AT^{-n}A=A なので、A=BA=B として定義の式は μ(A)=μ(A)2\mu(A)=\mu(A)^2、すなわち μ(A)∈{0,1}\mu(A)\in\{0,1\} です。

混合性を確かめるには、関数の言葉に直して稠密な部分集合で調べるのが便利です。以下 L2(μ)L^2(\mu) の内積を ⟨g,f⟩=∫gˉf dμ\langle g,f\rangle=\int\bar gf\,d\mu(第 2 変数について線形)とします。

補題 3.6(稠密な集合での混合)✓

D⊂L2(μ)\mathcal{D}\subset L^2(\mu) を稠密な部分集合とし、すべての f,g∈Df,g\in\mathcal{D} について ⟨g,UTnf⟩→⟨g,1⟩⟨1,f⟩\langle g,U_T^nf\rangle\to\langle g,1\rangle\langle1,f\rangle(n→∞n\to\infty)とする。このとき同じことがすべての f,g∈L2(μ)f,g\in L^2(\mu) で成り立ち、特に TT は混合的である。

証明(補題 3.6)

f,g∈L2f,g\in L^2 と ε>0\varepsilon>0 について f′,g′∈Df',g'\in\mathcal{D} を ∥f−f′∥2,∥g−g′∥2<ε\|f-f'\|_2,\|g-g'\|_2<\varepsilon と取ります。UTnU_T^n は等長(補題 2.2 の 3)なので Cauchy–Schwarz の不等式から

∣⟨g,UTnf⟩−⟨g′,UTnf′⟩∣≤∥g−g′∥2∥f∥2+∥g′∥2∥f−f′∥2≤ε(∥f∥2+∥g∥2+ε),|\langle g,U_T^nf\rangle-\langle g',U_T^nf'\rangle|\le\|g-g'\|_2\|f\|_2+\|g'\|_2\|f-f'\|_2\le\varepsilon(\|f\|_2+\|g\|_2+\varepsilon),

右辺の ⟨g,1⟩⟨1,f⟩\langle g,1\rangle\langle1,f\rangle も同じ程度しか変わりません。したがって lim sup⁡n∣⟨g,UTnf⟩−⟨g,1⟩⟨1,f⟩∣\limsup_n|\langle g,U_T^nf\rangle-\langle g,1\rangle\langle1,f\rangle| は ε\varepsilon の定数倍以下で、ε→0\varepsilon\to0 とすれば 00 です。f=1Af=\mathbf{1}_A、g=1Bg=\mathbf{1}_B とすると ⟨g,UTnf⟩=μ(B∩T−nA)\langle g,U_T^nf\rangle=\mu(B\cap T^{-n}A)、⟨g,1⟩⟨1,f⟩=μ(B)μ(A)\langle g,1\rangle\langle1,f\rangle=\mu(B)\mu(A) です。

3.1. 例:回転・2 倍写像・Bernoulli シフト

Section titled “3.1. 例:回転・2 倍写像・Bernoulli シフト”

T=[0,1)\mathbb{T}=[0,1) 上で ek(x)=e2πikxe_k(x)=e^{2\pi ikx}(k∈Zk\in\mathbb{Z})とおきます。Parseval の等式(定理 5.1[フーリエ級数])を変数変換 θ=2πx\theta=2\pi x で移すと、{ek}\{e_k\} は L2(T,m)L^2(\mathbb{T},m) の正規直交基底で、ff は Fourier 係数 f^(k)=⟨ek,f⟩\hat f(k)=\langle e_k,f\rangle で決まり、三角多項式(有限個の eke_k の線形結合)は L2L^2 で稠密です。

1 と 2 は「無理数のときだけ、不変な関数が定数に限られる」ことを Fourier 係数で確かめます。3 の証明の方針は、重なりの測度 μ(Rα−nA∩A)\mu(R_\alpha^{-n}A\cap A) が、ある nn では μ(A)2=14\mu(A)^2=\frac14 から遠いままであることを示すことです。そのために、nαn\alpha が整数にとても近くなる nn を鳩の巣原理で見つけます。

定理 3.7(円の回転のエルゴード性)✓

RαR_\alpha を 例 2.3 の回転とする。

  1. α∉Q\alpha\notin\mathbb{Q} なら RαR_\alpha はエルゴード的である。
  2. α∈Q\alpha\in\mathbb{Q} なら RαR_\alpha はエルゴード的でない。
  3. どの α\alpha についても RαR_\alpha は混合的でない。
証明(定理 3.7)

1。命題 3.3 の 4 を確かめます。f∈L2f\in L^2 が f(x+α)=f(x)f(x+\alpha)=f(x) a.e. を満たすとすると、平行移動 y=x+αy=x+\alpha で積分変数を変えて(Lebesgue 測度の平行移動不変性)

f^(k)=∫01f(x+α)e−2πikx dx=∫01f(y)e−2πik(y−α) dy=e2πikαf^(k)\hat f(k)=\int_0^1f(x+\alpha)e^{-2\pi ikx}\,dx=\int_0^1f(y)e^{-2\pi ik(y-\alpha)}\,dy=e^{2\pi ik\alpha}\hat f(k)

です。k≠0k\ne0 なら kα∉Zk\alpha\notin\mathbb{Z} なので e2πikα≠1e^{2\pi ik\alpha}\ne1、よって f^(k)=0\hat f(k)=0 です。Parseval の等式から ∥f−f^(0)∥22=∑k≠0∣f^(k)∣2=0\|f-\hat f(0)\|_2^2=\sum_{k\ne0}|\hat f(k)|^2=0 で、ff は a.e. で定数です。

2。α=p/q\alpha=p/q(q≥1q\ge1)なら eq(x+α)=e2πipeq(x)=eq(x)e_q(x+\alpha)=e^{2\pi ip}e_q(x)=e_q(x) で、eqe_q は定数でない不変関数です。

3。A=[0,12)A=[0,\frac12) とすると、Rα−nA=A−nαR_\alpha^{-n}A=A-n\alpha(mod 1)と AA の共通部分の測度は、t=nαt=n\alpha から最も近い整数までの距離を ∥t∥\|t\| として 12−∥nα∥\frac12-\|n\alpha\| です(長さ 12\frac12 の 2 つの弧を ∥t∥\|t\| だけずらした重なり)。α=p/q\alpha=p/q なら n=q,2q,…n=q,2q,\dots でこれは 12\frac12 です。α\alpha が無理数でも、NN を固定して 0,α,2α,…,Nα0,\alpha,2\alpha,\dots,N\alpha を mod 1 で考えると、[0,1)[0,1) を長さ 1N\frac1N の NN 個の区間に分けたとき、鳩の巣原理で 2 つが同じ区間に入り、その差から 1≤q≤N1\le q\le N で ∥qα∥≤1N\|q\alpha\|\le\frac1N となる qq が得られます。N→∞N\to\infty で ∥qjα∥→0\|q_j\alpha\|\to0、かつ α\alpha が無理数なので qj→∞q_j\to\infty と取れます(有界な qq では ∥qα∥\|q\alpha\| は正の最小値をもつ)。この部分列で重なりは 12\frac12 に近づき、μ(A)2=14\mu(A)^2=\frac14 に収束しません。

Harper 模型とほぼ Mathieu 作用素 の 補題 7.1[Harper 模型とほぼ Mathieu 作用素] は、すべての点で不変な有界関数についての同じ主張で、1 の特別な場合です。

01230, 5, 10, …α = (√5−1)/2α = 2/5
円の回転の軌道。左:α = (√5−1)/2(無理数)の 0〜29 歩。最初の 8 点を大きく描いた。点は円周上に一様に散らばっていく。右:α = 2/5 では 5 点を周期的にめぐるだけで、たとえば 5 点の小さな近傍の和は測度が 0 と 1 の中間の不変集合になる。

命題 3.8(2 倍写像は混合的)✓

例 2.4 の 2 倍写像 DD は Lebesgue 測度について混合的であり、特にエルゴード的である。

証明(命題 3.8)

ek∘D(x)=e2πik(2x)=e2k(x)e_k\circ D(x)=e^{2\pi ik(2x)}=e_{2k}(x)(2x2x と 2x mod 12x\bmod1 は整数だけ違う)なので UDnek=e2nkU_D^ne_k=e_{2^nk} です。k≠0k\ne0 と任意の ll について、2n∣k∣>∣l∣2^n|k|>|l| となる nn 以降は ⟨el,UDnek⟩=⟨el,e2nk⟩=0\langle e_l,U_D^ne_k\rangle=\langle e_l,e_{2^nk}\rangle=0 で、これは ⟨el,1⟩⟨1,ek⟩=0\langle e_l,1\rangle\langle1,e_k\rangle=0 に一致します。k=0k=0 なら両辺とも δl0\delta_{l0} です。双線形性から、三角多項式 f,gf,g について十分大きい nn で ⟨g,UDnf⟩=⟨g,1⟩⟨1,f⟩\langle g,U_D^nf\rangle=\langle g,1\rangle\langle1,f\rangle が成り立ちます。三角多項式は稠密なので 補題 3.6 から混合性が従い、命題 3.5 からエルゴード性が従います。

Bernoulli シフトでは Fourier 展開の代わりに、柱集合で一般の事象を近似します。

補題 3.9(代数による近似)✓

A\mathcal{A} を集合の代数、μ\mu を σ(A)\sigma(\mathcal{A}) 上の確率測度とする。任意の B∈σ(A)B\in\sigma(\mathcal{A}) と ε>0\varepsilon>0 について、μ(A △ B)<ε\mu(A\,\triangle\,B)<\varepsilon となる A∈AA\in\mathcal{A} が存在する。

証明(補題 3.9)

主張が成り立つ BB の全体を G\mathcal{G} とします。A⊂G\mathcal{A}\subset\mathcal{G} です。Ac △ Bc=A △ BA^c\,\triangle\,B^c=A\,\triangle\,B なので G\mathcal{G} は補集合で閉じています。B1,B2,⋯∈GB_1,B_2,\dots\in\mathcal{G}、B=⋃iBiB=\bigcup_iB_i とします。測度の連続性から μ(B∖⋃i≤NBi)<ε/2\mu(B\setminus\bigcup_{i\le N}B_i)<\varepsilon/2 となる NN があり、Ai∈AA_i\in\mathcal{A} を μ(Ai △ Bi)<ε/2i+1\mu(A_i\,\triangle\,B_i)<\varepsilon/2^{i+1} と取って A=⋃i≤NAi∈AA=\bigcup_{i\le N}A_i\in\mathcal{A} とすると

A △ B⊂(B∖⋃i≤NBi)∪⋃i≤N(Ai △ Bi)A\,\triangle\,B\subset\Bigl(B\setminus\bigcup_{i\le N}B_i\Bigr)\cup\bigcup_{i\le N}(A_i\,\triangle\,B_i)

なので μ(A △ B)<ε\mu(A\,\triangle\,B)<\varepsilon です。よって G\mathcal{G} は A\mathcal{A} を含む σ-加法族で、σ(A)⊂G\sigma(\mathcal{A})\subset\mathcal{G} です。

命題 3.10(Bernoulli シフトは混合的)✓

補題 2.5 の Bernoulli シフト (RN,F,ν⊗N,σ)(\mathbb{R}^{\mathbb{N}},\mathcal{F},\nu^{\otimes\mathbb{N}},\sigma) は混合的であり、特にエルゴード的である。

証明(命題 3.10)

A,BA,B が座標 1,…,m1,\dots,m で決まる柱集合、A={(x1,…,xm)∈A′}A=\{(x_1,\dots,x_m)\in A'\}、B={(x1,…,xm)∈B′}B=\{(x_1,\dots,x_m)\in B'\} のとき、n≥mn\ge m なら σ−nA={(xn+1,…,xn+m)∈A′}\sigma^{-n}A=\{(x_{n+1},\dots,x_{n+m})\in A'\} は BB と重ならない座標で決まるので、座標 1,…,n+m1,\dots,n+m の分布が ν⊗(n+m)\nu^{\otimes(n+m)} であること(補題 2.5)から P(σ−nA∩B)=ν⊗m(A′)ν⊗m(B′)=P(A)P(B)P(\sigma^{-n}A\cap B)=\nu^{\otimes m}(A')\nu^{\otimes m}(B')=P(A)P(B) です。一般の A,B∈FA,B\in\mathcal{F} と ε>0\varepsilon>0 について、補題 3.9 で柱集合 A0,B0A_0,B_0 を P(A △ A0),P(B △ B0)<εP(A\,\triangle\,A_0),P(B\,\triangle\,B_0)<\varepsilon と取ります(座標の個数は大きい方にそろえる)。P(σ−n(A △ A0))=P(A △ A0)P(\sigma^{-n}(A\,\triangle\,A_0))=P(A\,\triangle\,A_0) なので

∣P(σ−nA∩B)−P(σ−nA0∩B0)∣≤2ε,∣P(A)P(B)−P(A0)P(B0)∣≤2ε|P(\sigma^{-n}A\cap B)-P(\sigma^{-n}A_0\cap B_0)|\le2\varepsilon,\qquad|P(A)P(B)-P(A_0)P(B_0)|\le2\varepsilon

で、lim sup⁡n∣P(σ−nA∩B)−P(A)P(B)∣≤4ε\limsup_n|P(\sigma^{-n}A\cap B)-P(A)P(B)|\le4\varepsilon です。

4. von Neumann の平均エルゴード定理

Section titled “4. von Neumann の平均エルゴード定理”

まず L2L^2 での収束を、Hilbert 空間の一般論として示します。VV を Hilbert 空間 H\mathcal{H} 上の有界作用素で ∥V∥≤1\|V\|\le1 を満たすもの(縮小作用素)とします。Koopman 作用素 UTU_T は等長なので縮小作用素です。

この補題は、次の定理の証明の準備です。空間を「動かない部分」と「x−Vxx-Vx の形のベクトルで張られる部分」に直交分解し、前者では平均が変わらず、後者では平均が打ち消し合って 00 に近づく、という見通しを与えます。

補題 4.1(縮小作用素の不動点)✓

∥V∥≤1\|V\|\le1 なら、Vx=xVx=x と V∗x=xV^*x=x は同値である。さらに

H=ker⁡(I−V)⊕Ran(I−V)‾(直交直和)\mathcal{H}=\ker(I-V)\oplus\overline{\mathrm{Ran}(I-V)}\qquad\text{(直交直和)}

が成り立つ。

証明(補題 4.1)

Vx=xVx=x とします。∥V∗∥=∥V∥≤1\|V^*\|=\|V\|\le1(命題 3.3[有界作用素と共役作用素])と ⟨x,V∗x⟩=⟨Vx,x⟩=∥x∥2\langle x,V^*x\rangle=\langle Vx,x\rangle=\|x\|^2 から

∥V∗x−x∥2=∥V∗x∥2−2 Re⟨x,V∗x⟩+∥x∥2≤∥x∥2−2∥x∥2+∥x∥2=0\|V^*x-x\|^2=\|V^*x\|^2-2\,\mathrm{Re}\langle x,V^*x\rangle+\|x\|^2\le\|x\|^2-2\|x\|^2+\|x\|^2=0

です。V∗∗=VV^{**}=V なので逆も同じです。次に、y⊥Ran(I−V)y\perp\mathrm{Ran}(I-V) はすべての xx で ⟨(I−V∗)y,x⟩=⟨y,(I−V)x⟩=0\langle(I-V^*)y,x\rangle=\langle y,(I-V)x\rangle=0 となること、すなわち y∈ker⁡(I−V∗)=ker⁡(I−V)y\in\ker(I-V^*)=\ker(I-V) と同値です。M=Ran(I−V)‾M=\overline{\mathrm{Ran}(I-V)} は閉部分空間で M⊥=Ran(I−V)⊥=ker⁡(I−V)M^\perp=\mathrm{Ran}(I-V)^\perp=\ker(I-V) なので、射影定理(定理 3.4[Hilbert 空間の基礎] の 2)から主張の分解が得られます。

定理 4.2(von Neumann の平均エルゴード定理)✓

VV を Hilbert 空間 H\mathcal{H} 上の縮小作用素、PP を閉部分空間 ker⁡(I−V)\ker(I-V) への直交射影(定義 3.5[Hilbert 空間の基礎])とする。すべての x∈Hx\in\mathcal{H} について

∥1n∑k=0n−1Vkx−Px∥→0(n→∞).\Bigl\|\frac1n\sum_{k=0}^{n-1}V^kx-Px\Bigr\|\to0\qquad(n\to\infty).
証明(定理 4.2)

Mn=1n∑k<nVkM_n=\frac1n\sum_{k<n}V^k とおくと ∥Mn∥≤1\|M_n\|\le1 です。補題 4.1 の分解で x=Px+(x−Px)x=Px+(x-Px)、x−Px∈Ran(I−V)‾x-Px\in\overline{\mathrm{Ran}(I-V)} と書きます。PxPx は VV で動かないので MnPx=PxM_nPx=Px です。z=y−Vy∈Ran(I−V)z=y-Vy\in\mathrm{Ran}(I-V) なら和が打ち消し合って

Mnz=1n(y−Vny),∥Mnz∥≤2∥y∥n→0M_nz=\frac1n(y-V^ny),\qquad\|M_nz\|\le\frac{2\|y\|}n\to0

です。w∈Ran(I−V)‾w\in\overline{\mathrm{Ran}(I-V)} と ε>0\varepsilon>0 について ∥w−z∥<ε\|w-z\|<\varepsilon となる z∈Ran(I−V)z\in\mathrm{Ran}(I-V) を取れば、∥Mnw∥≤∥Mn(w−z)∥+∥Mnz∥≤ε+∥Mnz∥\|M_nw\|\le\|M_n(w-z)\|+\|M_nz\|\le\varepsilon+\|M_nz\| で、lim sup⁡n∥Mnw∥≤ε\limsup_n\|M_nw\|\le\varepsilon です。w=x−Pxw=x-Px として Mnx−Px=Mn(x−Px)→0M_nx-Px=M_n(x-Px)\to0 を得ます。

系 4.3(測度保存系での平均エルゴード定理)✓

測度保存系と f∈L2(μ)f\in L^2(\mu) について、AnfA_nf は L2(μ)L^2(\mu) で UTU_T の不動点空間への射影 PfPf に収束する。TT がエルゴード的なら Pf=∫f dμPf=\int f\,d\mu(定数関数)である。

証明(系 4.3)

Anf=1n∑k<nUTkfA_nf=\frac1n\sum_{k<n}U_T^kf なので、前半は 定理 4.2 を V=UTV=U_T に使ったものです。エルゴード的なら 命題 3.3 の 4 から ker⁡(I−UT)\ker(I-U_T) は定数関数の空間 C1\mathbb{C}1 で、∥1∥2=1\|1\|_2=1 なので射影は Pf=⟨1,f⟩1=∫f dμPf=\langle1,f\rangle1=\int f\,d\mu です。

例 4.4(回転での平均の計算)

α∉Q\alpha\notin\mathbb{Q}、k≠0k\ne0 とすると、URαek=e2πikαekU_{R_\alpha}e_k=e^{2\pi ik\alpha}e_k なので、等比級数の和から

Anek=1n∑j=0n−1e2πijkαek=1−e2πinkαn(1−e2πikα)ek,∥Anek∥∞≤2n∣1−e2πikα∣A_ne_k=\frac1n\sum_{j=0}^{n-1}e^{2\pi ijk\alpha}e_k=\frac{1-e^{2\pi ink\alpha}}{n(1-e^{2\pi ik\alpha})}e_k,\qquad\|A_ne_k\|_\infty\le\frac{2}{n|1-e^{2\pi ik\alpha}|}

です。収束は一様ですが、その速さ 1/(n∣1−e2πikα∣)1/(n|1-e^{2\pi ik\alpha}|) は kαk\alpha が整数にどれだけ近いかで決まり、kk について一様ではありません。平均エルゴード定理は収束を保証しますが、一般には速さについて何も言いません。

5. Birkhoff の各点エルゴード定理

Section titled “5. Birkhoff の各点エルゴード定理”

平均エルゴード定理は L2L^2 の収束で、個々の点 xx での時間平均については何も言いません。各点の収束には、時間平均の「最大値」を制御する不等式が要ります。以下 ff は実数値とします(複素数値なら実部と虚部に分けます)。

補題 5.1(最大エルゴード補題)✓

(X,B,μ,T)(X,\mathcal{B},\mu,T) を測度保存系、f∈L1(μ)f\in L^1(\mu) を実数値とし、N≥1N\ge1 について MN=max⁡0≤n≤NSnfM_N=\max_{0\le n\le N}S_nf(S0f=0S_0f=0 を含むので MN≥0M_N\ge0)とおく。このとき

∫{MN>0}f dμ≥0,∫{sup⁡n≥1Snf>0}f dμ≥0.\int_{\{M_N>0\}}f\,d\mu\ge0,\qquad\int_{\{\sup_{n\ge1}S_nf>0\}}f\,d\mu\ge0 .
証明(補題 5.1)

(Garsia の証明)0≤n≤N0\le n\le N について MN∘T≥Snf∘TM_N\circ T\ge S_nf\circ T なので、f+MN∘T≥f+Snf∘T=Sn+1ff+M_N\circ T\ge f+S_nf\circ T=S_{n+1}f です。nn について最大を取ると

f+MN∘T≥max⁡1≤n≤N+1Snf≥max⁡1≤n≤NSnf.f+M_N\circ T\ge\max_{1\le n\le N+1}S_nf\ge\max_{1\le n\le N}S_nf .

{MN>0}\{M_N>0\} の上では最大値は n=0n=0 の項 00 ではないので MN=max⁡1≤n≤NSnfM_N=\max_{1\le n\le N}S_nf で、したがって f≥MN−MN∘Tf\ge M_N-M_N\circ T です。∣MN∣≤∑k<N∣f∣∘Tk|M_N|\le\sum_{k<N}|f|\circ T^k は可積分(補題 2.2 の 2)で、{MN>0}\{M_N>0\} の外では MN=0M_N=0、また MN∘T≥0M_N\circ T\ge0 なので

∫{MN>0}f dμ≥∫{MN>0}(MN−MN∘T) dμ≥∫XMN dμ−∫XMN∘T dμ=0\int_{\{M_N>0\}}f\,d\mu\ge\int_{\{M_N>0\}}(M_N-M_N\circ T)\,d\mu\ge\int_XM_N\,d\mu-\int_XM_N\circ T\,d\mu=0

です(最後は 補題 2.2 の 2)。N→∞N\to\infty で {MN>0}\{M_N>0\} は {sup⁡n≥1Snf>0}\{\sup_{n\ge1}S_nf>0\} に増加し、∣f1{MN>0}∣≤∣f∣|f\mathbf{1}_{\{M_N>0\}}|\le|f| なので、優収束定理(定理 7.3[ルベーグ積分の定義と収束定理])で 2 つ目の不等式を得ます。

この補題の気持ちは「時間平均がどこかで正になる点の集合の上では、ff の平均は非負」ということです。次の定理ではこれを f−bf-b や a−fa-f に使い、時間平均の上極限と下極限が食い違う集合の測度を 00 に押さえ込みます。

証明の方針は次のとおりです。時間平均の上極限 f∗f^* と下極限 f∗f_* が食い違う点の集合は不変集合になります。その上で ff から定数 bb を引く(または aa から ff を引く)と、その集合の上では時間平均が正になる点がある状況をつくれるので、最大エルゴード補題から積分の符号が決まります。2 つを足すと、食い違う集合の測度が 00 だと分かります。

定理 5.2(Birkhoff の各点エルゴード定理)✓

(X,B,μ,T)(X,\mathcal{B},\mu,T) を測度保存系、I\mathcal{I} を不変 σ-加法族、f∈L1(μ)f\in L^1(\mu) とする。I\mathcal{I}-可測な fˉ∈L1(μ)\bar f\in L^1(\mu) が存在して次が成り立つ。

  1. ほとんどすべての xx で Anf(x)=1n∑k=0n−1f(Tkx)→fˉ(x)A_nf(x)=\frac1n\sum_{k=0}^{n-1}f(T^kx)\to\bar f(x)。
  2. ∥Anf−fˉ∥1→0\|A_nf-\bar f\|_1\to0、∥fˉ∥1≤∥f∥1\|\bar f\|_1\le\|f\|_1。
  3. すべての B∈IB\in\mathcal{I} で ∫Bfˉ dμ=∫Bf dμ\int_B\bar f\,d\mu=\int_Bf\,d\mu。すなわち fˉ=E[f∣I]\bar f=E[f\mid\mathcal{I}] である。
  4. TT がエルゴード的なら、fˉ=∫f dμ\bar f=\int f\,d\mu a.e. である。
証明(定理 5.2)

ff は実数値とし、代表元をすべての点で有限な値をとるように選びます。

不変性。 f∗=lim sup⁡nAnff^*=\limsup_nA_nf、f∗=lim inf⁡nAnff_*=\liminf_nA_nf(値は [−∞,∞][-\infty,\infty])とおきます。(n+1)An+1f(x)=f(x)+nAnf(Tx)(n+1)A_{n+1}f(x)=f(x)+nA_nf(Tx) から

An+1f(x)=f(x)n+1+nn+1Anf(Tx)A_{n+1}f(x)=\frac{f(x)}{n+1}+\frac n{n+1}A_nf(Tx)

で、f(x)f(x) は有限なので、すべての xx で f∗(Tx)=f∗(x)f^*(Tx)=f^*(x)、f∗(Tx)=f∗(x)f_*(Tx)=f_*(x) です。したがって {f∗>c}\{f^*>c\} などは不変集合で、f∗,f∗f^*,f_* は I\mathcal{I}-可測です。

上極限と下極限の一致。 有理数 a<ba<b について E={f∗<a}∩{f∗>b}∈IE=\{f_*<a\}\cap\{f^*>b\}\in\mathcal{I} とおきます。g=(f−b)1Eg=(f-b)\mathbf{1}_E とすると、1E∘Tk=1T−kE=1E\mathbf{1}_E\circ T^k=\mathbf{1}_{T^{-k}E}=\mathbf{1}_E なので Sng=1E (Snf−nb)=n1E(Anf−b)S_ng=\mathbf{1}_E\,(S_nf-nb)=n\mathbf{1}_E(A_nf-b) です。EE の上では lim sup⁡Anf>b\limsup A_nf>b なので sup⁡n≥1Sng>0\sup_{n\ge1}S_ng>0、EE の外では Sng=0S_ng=0 です。よって {sup⁡n≥1Sng>0}=E\{\sup_{n\ge1}S_ng>0\}=E で、補題 5.1 から ∫E(f−b) dμ≥0\int_E(f-b)\,d\mu\ge0 です。同様に h=(a−f)1Eh=(a-f)\mathbf{1}_E について、EE の上では lim inf⁡Anf<a\liminf A_nf<a なので {sup⁡n≥1Snh>0}=E\{\sup_{n\ge1}S_nh>0\}=E で、∫E(a−f) dμ≥0\int_E(a-f)\,d\mu\ge0 です。足すと (a−b)μ(E)≥0(a-b)\mu(E)\ge0、すなわち μ(E)=0\mu(E)=0 です。a<ba<b を有理数全体で動かした可算和も零集合なので、ほとんどすべての点で f∗=f∗f_*=f^*、すなわち AnfA_nf は [−∞,∞][-\infty,\infty] の中で収束します。

可積分性。 ∣Anf∣≤An∣f∣|A_nf|\le A_n|f| と 補題 2.2 から ∫∣Anf∣ dμ≤∥f∥1\int|A_nf|\,d\mu\le\|f\|_1 です。Fatou の補題(補題 6.1[ルベーグ積分の定義と収束定理])で ∫∣f∗∣ dμ≤lim inf⁡n∫∣Anf∣ dμ≤∥f∥1\int|f^*|\,d\mu\le\liminf_n\int|A_nf|\,d\mu\le\|f\|_1 なので、f∗f^* は a.e. で有限です。fˉ=f∗1{∣f∗∣<∞}\bar f=f^*\mathbf{1}_{\{|f^*|<\infty\}} とおくと({∣f∗∣<∞}∈I\{|f^*|<\infty\}\in\mathcal{I})、fˉ\bar f は I\mathcal{I}-可測、可積分で、1 と ∥fˉ∥1≤∥f∥1\|\bar f\|_1\le\|f\|_1 が成り立ちます。

L1L^1 収束。 ∣f∣≤K|f|\le K なら ∣Anf∣≤K|A_nf|\le K なので、優収束定理(定理 7.3[ルベーグ積分の定義と収束定理])で ∥Anf−fˉ∥1→0\|A_nf-\bar f\|_1\to0 です。一般の ff では fK=max⁡(−K,min⁡(f,K))f_K=\max(-K,\min(f,K)) とおくと、優収束定理で ∥f−fK∥1→0\|f-f_K\|_1\to0(K→∞K\to\infty)です。AnA_n は線形なので a.e. で fˉ−fK‾=lim⁡nAn(f−fK)\bar f-\overline{f_K}=\lim_nA_n(f-f_K) で、上の可積分性の議論を f−fKf-f_K に使えば ∥fˉ−fK‾∥1≤∥f−fK∥1\|\bar f-\overline{f_K}\|_1\le\|f-f_K\|_1 です。したがって

∥Anf−fˉ∥1≤∥An(f−fK)∥1+∥AnfK−fK‾∥1+∥fK‾−fˉ∥1≤2∥f−fK∥1+∥AnfK−fK‾∥1\|A_nf-\bar f\|_1\le\|A_n(f-f_K)\|_1+\|A_nf_K-\overline{f_K}\|_1+\|\overline{f_K}-\bar f\|_1\le2\|f-f_K\|_1+\|A_nf_K-\overline{f_K}\|_1

で、KK を大きく取ってから n→∞n\to\infty とすれば 00 に近づきます。

条件付き期待値。 B∈IB\in\mathcal{I} なら 1B=1B∘Tk\mathbf{1}_B=\mathbf{1}_B\circ T^k なので、補題 2.2 から ∫Bf∘Tk dμ=∫(1Bf)∘Tk dμ=∫Bf dμ\int_Bf\circ T^k\,d\mu=\int(\mathbf{1}_Bf)\circ T^k\,d\mu=\int_Bf\,d\mu、よって ∫BAnf dμ=∫Bf dμ\int_BA_nf\,d\mu=\int_Bf\,d\mu です。L1L^1 収束から ∫Bfˉ dμ=∫Bf dμ\int_B\bar f\,d\mu=\int_Bf\,d\mu です。fˉ\bar f は I\mathcal{I}-可測で可積分なので、これは条件付き期待値の定義(定義 3.1[条件付き期待値])の 3 条件で、一意性(定理 3.2[条件付き期待値])から fˉ=E[f∣I]\bar f=E[f\mid\mathcal{I}] です。

エルゴード的な場合。 fˉ\bar f は I\mathcal{I}-可測なので fˉ∘T=fˉ\bar f\circ T=\bar f で、命題 3.3 の 3 から a.e. で定数 cc です。3 で B=XB=X として c=∫f dμc=\int f\,d\mu です。

f∈L2f\in L^2 なら、系 4.3 の L2L^2 極限 PfPf と fˉ\bar f は一致します。確率空間では ∥⋅∥1≤∥⋅∥2\|\cdot\|_1\le\|\cdot\|_2(Cauchy–Schwarz)なので AnfA_nf は L1L^1 で PfPf にも fˉ\bar f にも収束し、L1L^1 の極限は一意だからです。特に PP は L2L^2 上の条件付き期待値 E[⋅∣I]E[\cdot\mid\mathcal{I}] で、これは 定理 5.1[条件付き期待値] と整合します。

注意 5.3(仮定について)

  • 測度の有限性は落とせません。R\mathbb{R} 上の x↦x+1x\mapsto x+1 と Lebesgue 測度、f=1[0,1)f=\mathbf{1}_{[0,1)} では、すべての xx で Anf(x)→0A_nf(x)\to0 ですが、空間平均は意味をもちません(エルゴード定理の 4 にあたる主張が無い)。σ\sigma-有限な無限測度でも 1 の a.e. 収束そのものは成り立ちますが、この記事では扱いません。
  • 「ほとんどすべての xx」は「すべての xx」にできません。2 倍写像の周期点 x=13=0.0101…(2)x=\frac13=0.0101\ldots_{(2)} では 1[12,1)\mathbf{1}_{[\frac12,1)} の時間平均は 12\frac12 ですが、x=17=0.001001…(2)x=\frac17=0.001001\ldots_{(2)} では 13\frac13 です(演習 6.1)。

5.1. 応用:大数の強法則と正規数

Section titled “5.1. 応用:大数の強法則と正規数”

系 5.4(大数の強法則(エルゴード定理による証明))✓

X1,X2,…X_1,X_2,\dots を確率空間 (Ω′,F′,P′)(\Omega',\mathcal{F}',P') 上の独立同分布な確率変数列で E[∣X1∣]<∞E[|X_1|]<\infty とする。このとき P′P'-a.s. で 1n∑k=1nXk→E[X1]\frac1n\sum_{k=1}^nX_k\to E[X_1] である。

証明(系 5.4)

ν\nu を X1X_1 の分布、(RN,F,P,σ)(\mathbb{R}^{\mathbb{N}},\mathcal{F},P,\sigma) を 補題 2.5 の Bernoulli シフトとします。Φ(ω)=(X1(ω),X2(ω),… )\Phi(\omega)=(X_1(\omega),X_2(\omega),\dots) は各座標が可測なので、座標写像で生成される直積 σ-加法族 F\mathcal{F} について可測です。独立性から (X1,…,Xm)(X_1,\dots,X_m) の分布は直方体の上で ν⊗m\nu^{\otimes m} と一致し、直方体は π-系なので測度の一意性(定理 4.5[測度空間と Carathéodory の拡張定理])から分布は ν⊗m\nu^{\otimes m} です。よって Q=P′∘Φ−1Q=P'\circ\Phi^{-1} は柱集合の上で PP と一致し、同じ一意性から Q=PQ=P です。

f(x)=x1f(x)=x_1 は ∫∣f∣ dP=E[∣X1∣]<∞\int|f|\,dP=E[|X_1|]<\infty を満たし、f∘σk(x)=xk+1f\circ\sigma^k(x)=x_{k+1} なので Anf(x)=1n∑k=1nxkA_nf(x)=\frac1n\sum_{k=1}^nx_k です。命題 3.10 からシフトはエルゴード的なので、定理 5.2 の 4 により、集合 G={x:1n∑k=1nxk→∫x1 dP}G=\{x:\frac1n\sum_{k=1}^nx_k\to\int x_1\,dP\}(各点の上極限・下極限が可測なので G∈FG\in\mathcal{F})は P(G)=1P(G)=1 です。∫x1 dP=E[X1]\int x_1\,dP=E[X_1] で、P′(Φ−1G)=Q(G)=P(G)=1P'(\Phi^{-1}G)=Q(G)=P(G)=1 です。

大数の法則と中心極限定理 の 定理 4.5[大数の法則と中心極限定理] の証明は独立性を分散の計算に使いますが、上の証明は「シフト不変な事象が自明である」ことしか使っていません。そのため、独立でなくても定常でエルゴード的な列(たとえば既約で非周期的な有限 Markov 連鎖の定常列)に同じ結論が成り立ちます。

例 5.5(Borel の正規数定理(2 進))

x∈[0,1)x\in[0,1) の 2 進展開 x=0.x1x2…x=0.x_1x_2\ldots(1 が無限に続く展開は使わない)について、ほとんどすべての xx で

lim⁡n→∞#{k≤n:xk=1}n=12\lim_{n\to\infty}\frac{\#\{k\le n:x_k=1\}}n=\frac12

です。実際、例 2.4 の 2 倍写像について Dk−1xD^{k-1}x の第 1 桁が xkx_k で、xk=1[12,1)(Dk−1x)x_k=\mathbf{1}_{[\frac12,1)}(D^{k-1}x) です。命題 3.8 から DD はエルゴード的なので、定理 5.2 の 4 から時間平均は m([12,1))=12m([\frac12,1))=\frac12 に収束します。同じ議論を 1[j/2ℓ,(j+1)/2ℓ)\mathbf{1}_{[j/2^\ell,(j+1)/2^\ell)} に使えば、長さ ℓ\ell のどの 0-1 列も頻度 2−ℓ2^{-\ell} で現れることが分かります。

注意 5.6(エルゴード分解(主張のみ))

エルゴード的でない系は、エルゴード的な系の「重ね合わせ」に分解できます。XX がコンパクト距離空間で TT が同相写像のとき、TT 不変な確率測度の全体は凸集合で、その端点はちょうどエルゴード的な測度であり、各不変測度はエルゴード測度の上の確率測度の重心としてただ 1 通りに書けます(状態の凸分解 の 命題 7.6[状態の凸分解])。より一般に XX が標準 Borel 空間なら、μ=∫μx dμ(x)\mu=\int\mu_x\,d\mu(x) と、不変 σ-加法族 I\mathcal{I} に関する条件付き分布 μx\mu_x がほとんどすべての xx でエルゴード的になるように分解できます。Birkhoff の極限 fˉ(x)=E[f∣I](x)\bar f(x)=E[f\mid\mathcal{I}](x) は、xx が属する「エルゴード成分」μx\mu_x での空間平均 ∫f dμx\int f\,d\mu_x です。証明には条件付き分布の存在が要るので、この記事では扱いません(Einsiedler–Ward 第 6 章)。

演習 6.1易

2 倍写像 DD と f=1[12,1)f=\mathbf{1}_{[\frac12,1)} について、x=13x=\frac13 と x=17x=\frac17 での時間平均 lim⁡nAnf(x)\lim_nA_nf(x) を求めてください。結果が 例 5.5 と矛盾しない理由を述べてください。

解答

13=0.010101…(2)\frac13=0.010101\ldots_{(2)} で、D(13)=23D(\frac13)=\frac23、D(23)=13D(\frac23)=\frac13 と周期 2 で回ります。f(13)=0f(\frac13)=0、f(23)=1f(\frac23)=1 なので Anf(13)→12A_nf(\frac13)\to\frac12 です。17=0.001001…(2)\frac17=0.001001\ldots_{(2)} で、軌道は 17→27→47→17\frac17\to\frac27\to\frac47\to\frac17、ff の値は 0,0,10,0,1 なので Anf(17)→13A_nf(\frac17)\to\frac13 です。有理数(さらに一般に周期点や最終的に周期的な点)は可算個で Lebesgue 測度 00 なので、「ほとんどすべての xx で 12\frac12」という主張とは矛盾しません。エルゴード定理の例外集合は空でないことがあり、しかも稠密になりえます。

演習 6.2標準

測度保存系がエルゴード的であるための必要十分条件は、すべての A,B∈BA,B\in\mathcal{B} について

lim⁡n→∞1n∑k=0n−1μ(T−kA∩B)=μ(A)μ(B)\lim_{n\to\infty}\frac1n\sum_{k=0}^{n-1}\mu(T^{-k}A\cap B)=\mu(A)\mu(B)

が成り立つことであることを示してください。混合性(定義 3.4)との違いも述べてください。

解答

(必要性)1n∑k<nμ(T−kA∩B)=⟨1B,An1A⟩\frac1n\sum_{k<n}\mu(T^{-k}A\cap B)=\langle\mathbf{1}_B,A_n\mathbf{1}_A\rangle です。系 4.3 からエルゴード的なら An1A→μ(A)A_n\mathbf{1}_A\to\mu(A) が L2L^2 で成り立ち、内積の連続性から左辺は ⟨1B,μ(A)1⟩=μ(A)μ(B)\langle\mathbf{1}_B,\mu(A)1\rangle=\mu(A)\mu(B) に収束します。

(十分性)A∈IA\in\mathcal{I} を取り B=AB=A とすると T−kA=AT^{-k}A=A なので左辺は μ(A)\mu(A)、右辺は μ(A)2\mu(A)^2 で、μ(A)∈{0,1}\mu(A)\in\{0,1\} です。

混合性は数列 μ(T−kA∩B)\mu(T^{-k}A\cap B) そのものの収束、エルゴード性はその Cesàro 平均の収束です。無理数回転(定理 3.7)は後者だけを満たす例です。

演習 6.3難

α∉Q\alpha\notin\mathbb{Q} とします。T\mathbb{T} 上のすべての連続関数 ff とすべての x∈Tx\in\mathbb{T} について

1n∑k=0n−1f(x+kα)→∫01f dm\frac1n\sum_{k=0}^{n-1}f(x+k\alpha)\to\int_0^1f\,dm

が成り立つこと(Weyl の一様分布定理)を示してください。不連続な ff では「すべての xx」が成り立たない例も挙げてください。

解答

f=ekf=e_k(k≠0k\ne0)なら 例 4.4 の計算で sup⁡x∣Anek(x)∣≤2/(n∣1−e2πikα∣)→0=∫ek dm\sup_x|A_ne_k(x)|\le2/(n|1-e^{2\pi ik\alpha}|)\to0=\int e_k\,dm、f=e0f=e_0 なら両辺 11 です。線形性から三角多項式 pp について Anp→∫p dmA_np\to\int p\,dm が xx について一様に成り立ちます。連続な ff と ε>0\varepsilon>0 について、Fejér の定理(定理 4.2[フーリエ級数] の 1 を θ=2πx\theta=2\pi x で移したもの)から sup⁡∣f−p∣<ε\sup|f-p|<\varepsilon となる三角多項式 pp があります。∣Anf−Anp∣≤ε|A_nf-A_np|\le\varepsilon、∣∫f−∫p∣≤ε|\int f-\int p|\le\varepsilon なので lim sup⁡nsup⁡x∣Anf(x)−∫f∣≤2ε\limsup_n\sup_x|A_nf(x)-\int f|\le2\varepsilon です。

不連続な例:f=1{kα mod 1:k≥0}f=\mathbf{1}_{\{k\alpha\bmod1:k\ge0\}}(00 の軌道の指示関数)は可算集合の指示関数なので ∫f=0\int f=0 ですが、x=0x=0 ではすべての項が 11 で時間平均は 11 です。

  • P. Walters, An Introduction to Ergodic Theory, Springer, 1982(Graduate Texts in Mathematics 79)— 第 1 章(測度保存系、エルゴード性、混合)。
  • M. Einsiedler, T. Ward, Ergodic Theory with a view towards Number Theory, Springer, 2011(Graduate Texts in Mathematics 259)— 第 2 章(平均・各点エルゴード定理)、第 6 章(エルゴード分解)。
  • K. Petersen, Ergodic Theory, Cambridge University Press, 1983 — 第 2 章。
  • A. M. Garsia, “A simple proof of E. Hopf’s maximal ergodic theorem”, J. Math. Mech. 14 (1965), 381–382.
  • G. D. Birkhoff, “Proof of the ergodic theorem”, Proc. Natl. Acad. Sci. USA 17 (1931), 656–660.
  • J. von Neumann, “Proof of the quasi-ergodic hypothesis”, Proc. Natl. Acad. Sci. USA 18 (1932), 70–82.

このテーマの 100 本ノック(100 問)

執筆: Claude Code(Opus 5.5)、Claude Code(Sonnet 5.5)・監査 3 回(Claude Code)・最終 2026-10-12

この記事の誤りを報告する ・運営: 夢現技研合同会社 ・料金 ・利用条件 ・特定商取引法に基づく表記

© 2026 夢現技研合同会社 ・コード例は MIT ライセンスです。