ある規則で点を動かし続けたとき、1 つの点の動きを長く平均したものと、空間全体の平均とは、いつ同じになるでしょうか。この記事は、その問いに答える 2 つの定理(時間平均が収束する、その極限が空間平均になる)を、例を確かめながら証明します。読み終えると、サイコロを何度も振った平均が期待値に近づく大数の強法則が、この一般論の特別な場合として出てくることが分かります。以下の要点は全体の見取り図なので、用語が分からなくても先へ進んでかまいません。
確率空間 ( X , B , μ ) (X,\mathcal{B},\mu) ( X , B , μ ) 上の可測写像 T T T が μ ( T − 1 A ) = μ ( A ) \mu(T^{-1}A)=\mu(A) μ ( T − 1 A ) = μ ( A ) を満たすとき測度保存系 といいます。T T T で動かしても消えない集合が自明なものだけのときエルゴード的 といい、これは「T T T 不変な関数は定数だけ」と同値です。
円の無理数回転はエルゴード的ですが混合的ではありません。2 倍写像と Bernoulli シフトは混合的です。どれも Fourier 展開か「柱集合による近似」(座標を有限個だけ指定して決まる集合で、全体を近似すること)で示せます。
von Neumann の平均エルゴード定理 : Hilbert 空間 の縮小作用素(ノルムを増やさない線形作用素)V V V について、1 n ∑ k < n V k x \frac1n\sum_{k<n}V^kx n 1 ∑ k < n V k x は不動点空間への直交射影 P x Px P x にノルム収束します。証明は H = ker ( I − V ) ⊕ R a n ( I − V ) ‾ \mathcal{H}=\ker(I-V)\oplus\overline{\mathrm{Ran}(I-V)} H = ker ( I − V ) ⊕ Ran ( I − V ) という分解だけです。
Birkhoff の各点エルゴード定理 : f ∈ L 1 ( μ ) f\in L^1(\mu) f ∈ L 1 ( μ ) の時間平均 1 n ∑ k < n f ( T k x ) \frac1n\sum_{k<n}f(T^kx) n 1 ∑ k < n f ( T k x ) はほとんどすべての x x x で収束し、極限は不変 σ-加法族への条件付き期待値 E [ f ∣ I ] E[f\mid\mathcal{I}] E [ f ∣ I ] です。エルゴード的なら極限は空間平均 ∫ f d μ \int f\,d\mu ∫ f d μ です。鍵は最大エルゴード補題 (Garsia の 3 行の証明)です。
応用として、独立同分布の確率変数列の大数の強法則が Bernoulli シフトのエルゴード性から従います。次章では、この定理を劣加法的な列に拡張して Lyapunov 指数を扱います。
この記事で使う事前知識
Hilbert 空間 : 内積をもち、そこから決まる距離で完備な線形空間。→ Hilbert 空間
直交射影 : 閉部分空間への、最も近い点を対応させる写像。→ 直交射影
正規直交基底 : 長さ 1 で互いに直交するベクトルの組で、すべてのベクトルをその線形結合の極限で表せるもの。→ 正規直交基底
条件付き期待値 : 部分 σ-加法族に属する集合の上では、元の関数と同じ積分をもつ、部分 σ-加法族に関して可測な関数。→ 条件付き期待値
直積 σ-加法族 : 座標写像がすべて可測になる最小の σ-加法族。→ 直積 σ-加法族
概収束 : ほとんどすべての点で値が収束すること。→ 概収束
気体の分子を 1 つ追いかけ、ある量(たとえば運動エネルギー)を長い時間にわたって平均したもの(時間平均 )と、ある瞬間にすべての分子について平均したもの(空間平均 、相平均)は一致するでしょうか。Boltzmann は統計力学の基礎づけにこの一致を仮定しました(エルゴード仮説)。これを数学の問題にしたのが 1931〜32 年の von Neumann と Birkhoff の定理で、「軌道がほとんどすべての点で空間全体をまんべんなく訪れる」ための条件としてエルゴード性 を取り出しました。
確率論から見ると、エルゴード定理は大数の強法則の一般化です。独立同分布の列 X 1 , X 2 , … X_1,X_2,\dots X 1 , X 2 , … は、列全体の空間 R N \mathbb{R}^{\mathbb{N}} R N の上で「1 つずらす」写像(シフト)の軌道とみなせます。独立性は使わず、ずらしても分布が変わらない(定常性)ことと、不変な事象が自明であること(エルゴード性)だけで、標本平均の概収束が言えます。
この章の結果は、準周期的な Schrödinger 作用素(Harper 模型とほぼ Mathieu 作用素 )で、位相 θ \theta θ を θ + α \theta+\alpha θ + α に動かす無理数回転の上の「転送行列の積」を調べるときの土台になります。
以下 ( X , B , μ ) (X,\mathcal{B},\mu) ( X , B , μ ) は確率空間(μ ( X ) = 1 \mu(X)=1 μ ( X ) = 1 )です。
定義 2.1 (測度保存系 ) ↗
可測写像 T : X → X T:X\to X T : X → X がすべての A ∈ B A\in\mathcal{B} A ∈ B について μ ( T − 1 A ) = μ ( A ) \mu(T^{-1}A)=\mu(A) μ ( T − 1 A ) = μ ( A ) を満たすとき、T T T は μ \mu μ を保存する といい、( X , B , μ , T ) (X,\mathcal{B},\mu,T) ( X , B , μ , T ) を測度保存系 という。T 0 = i d T^0=\mathrm{id} T 0 = id 、T n = T ∘ T n − 1 T^{n}=T\circ T^{n-1} T n = T ∘ T n − 1 とし、可測関数 f f f について
S n f = ∑ k = 0 n − 1 f ∘ T k , A n f = 1 n S n f ( n ≥ 1 ) , S 0 f = 0 S_nf=\sum_{k=0}^{n-1}f\circ T^k,\qquad A_nf=\frac1nS_nf\quad(n\ge1),\qquad S_0f=0 S n f = k = 0 ∑ n − 1 f ∘ T k , A n f = n 1 S n f ( n ≥ 1 ) , S 0 f = 0 を Birkhoff 和 、Birkhoff 平均 (時間平均)という。
T T T が可逆である必要はありません。条件は逆像 T − 1 A T^{-1}A T − 1 A について書かれていて、これは「T T T で 1 歩進んだ後に A A A にいる点の集合」です。
補題 2.2 (積分の不変性と Koopman 作用素 ) ✓
( X , B , μ , T ) (X,\mathcal{B},\mu,T) ( X , B , μ , T ) を測度保存系とする。
可測関数 f : X → [ 0 , ∞ ] f:X\to[0,\infty] f : X → [ 0 , ∞ ] について ∫ f ∘ T d μ = ∫ f d μ \int f\circ T\,d\mu=\int f\,d\mu ∫ f ∘ T d μ = ∫ f d μ 。
f ∈ L 1 ( μ ) f\in L^1(\mu) f ∈ L 1 ( μ ) なら f ∘ T ∈ L 1 ( μ ) f\circ T\in L^1(\mu) f ∘ T ∈ L 1 ( μ ) で ∫ f ∘ T d μ = ∫ f d μ \int f\circ T\,d\mu=\int f\,d\mu ∫ f ∘ T d μ = ∫ f d μ 。
1 ≤ p < ∞ 1\le p<\infty 1 ≤ p < ∞ について U T f = f ∘ T U_Tf=f\circ T U T f = f ∘ T は L p ( μ ) L^p(\mu) L p ( μ ) 上の well-defined な線形等長作用素である。特に ∥ A n f ∥ p ≤ ∥ f ∥ p \|A_nf\|_p\le\|f\|_p ∥ A n f ∥ p ≤ ∥ f ∥ p 。
証明(補題 2.2) 1。f = 1 A f=\mathbf{1}_A f = 1 A なら 1 A ∘ T = 1 T − 1 A \mathbf{1}_A\circ T=\mathbf{1}_{T^{-1}A} 1 A ∘ T = 1 T − 1 A で、両辺は μ ( T − 1 A ) = μ ( A ) \mu(T^{-1}A)=\mu(A) μ ( T − 1 A ) = μ ( A ) です。非負単関数については線形性から従います。一般の f f f には単関数近似定理(定理 3.7[ルベーグ積分の定義と収束定理] )で 0 ≤ s k ↑ f 0\le s_k\uparrow f 0 ≤ s k ↑ f を取ると s k ∘ T ↑ f ∘ T s_k\circ T\uparrow f\circ T s k ∘ T ↑ f ∘ T なので、単調収束定理(定理 5.2[ルベーグ積分の定義と収束定理] )で両辺の極限を取れます。
2。f = f + − f − f=f^+-f^- f = f + − f − に分けて 1 を使います。
3。f = g f=g f = g a.e. なら { f ∘ T ≠ g ∘ T } = T − 1 { f ≠ g } \{f\circ T\ne g\circ T\}=T^{-1}\{f\ne g\} { f ∘ T = g ∘ T } = T − 1 { f = g } は測度 0 0 0 なので、U T U_T U T は同値類の上で定まります。1 を ∣ f ∣ p |f|^p ∣ f ∣ p に使うと ∥ f ∘ T ∥ p = ∥ f ∥ p \|f\circ T\|_p=\|f\|_p ∥ f ∘ T ∥ p = ∥ f ∥ p です。A n f A_nf A n f については Minkowski の不等式(定理 4.5[L^p 空間と関数解析への導入] )から ∥ A n f ∥ p ≤ 1 n ∑ k < n ∥ f ∘ T k ∥ p = ∥ f ∥ p \|A_nf\|_p\le\frac1n\sum_{k<n}\|f\circ T^k\|_p=\|f\|_p ∥ A n f ∥ p ≤ n 1 ∑ k < n ∥ f ∘ T k ∥ p = ∥ f ∥ p です。
∎
U T U_T U T を Koopman 作用素 といいます。非線形な写像 T T T を、関数空間の上の線形作用素に置き換えるのがエルゴード理論の基本的な発想です。
例 2.3 (円の回転 )
T = R / Z \mathbb{T}=\mathbb{R}/\mathbb{Z} T = R / Z を [ 0 , 1 ) [0,1) [ 0 , 1 ) と同一視し、Lebesgue 測度 m m m を入れます。α ∈ R \alpha\in\mathbb{R} α ∈ R について R α x = x + α ( m o d 1 ) R_\alpha x=x+\alpha\pmod 1 R α x = x + α ( mod 1 ) とすると、R α − 1 A R_\alpha^{-1}A R α − 1 A は A A A を − α -\alpha − α だけ平行移動して [ 0 , 1 ) [0,1) [ 0 , 1 ) に折り返した集合で、Lebesgue 測度の平行移動不変性から m ( R α − 1 A ) = m ( A ) m(R_\alpha^{-1}A)=m(A) m ( R α − 1 A ) = m ( A ) です。
例 2.4 (2 倍写像 )
D x = 2 x ( m o d 1 ) D x=2x\pmod1 D x = 2 x ( mod 1 ) を [ 0 , 1 ) [0,1) [ 0 , 1 ) 上で考えます。区間について D − 1 [ a , b ) = [ a 2 , b 2 ) ∪ [ a + 1 2 , b + 1 2 ) D^{-1}[a,b)=[\tfrac a2,\tfrac b2)\cup[\tfrac{a+1}2,\tfrac{b+1}2) D − 1 [ a , b ) = [ 2 a , 2 b ) ∪ [ 2 a + 1 , 2 b + 1 ) なので m ( D − 1 [ a , b ) ) = b − a m(D^{-1}[a,b))=b-a m ( D − 1 [ a , b )) = b − a です。区間 [ a , b ) [a,b) [ a , b ) (と空集合)は π-系で Borel 集合族を生成するので、測度の一意性(定理 4.5[測度空間と Carathéodory の拡張定理] )から m ∘ D − 1 = m m\circ D^{-1}=m m ∘ D − 1 = m です。D D D は 2 対 1 で可逆ではありません。x x x の 2 進展開 x = 0. x 1 x 2 x 3 … x=0.x_1x_2x_3\dots x = 0. x 1 x 2 x 3 … に対して D x = 0. x 2 x 3 … Dx=0.x_2x_3\dots D x = 0. x 2 x 3 … なので、D D D は「2 進展開の桁を 1 つずらす」写像です。
次の補題は、コインを無限回投げる試行のように、無限個の座標をもつ空間に測度を入れる話です。座標を最初の m m m 個だけ指定した集合(柱集合)の測度を ν ⊗ m \nu^{\otimes m} ν ⊗ m で決めておき、それを全体へ延ばすのが Kolmogorov の拡張定理の役目です。
補題 2.5 (Bernoulli シフトの構成 ) ✓
ν \nu ν を R \mathbb{R} R 上の確率測度、Ω = R N \Omega=\mathbb{R}^{\mathbb{N}} Ω = R N 、F \mathcal{F} F を座標写像 x ↦ x k x\mapsto x_k x ↦ x k で生成される直積 σ-加法族とする。F \mathcal{F} F 上の確率測度 P P P で、任意の m m m と Borel 集合 B ⊂ R m B\subset\mathbb{R}^m B ⊂ R m について P ( { ( x 1 , … , x m ) ∈ B } ) = ν ⊗ m ( B ) P(\{(x_1,\dots,x_m)\in B\})=\nu^{\otimes m}(B) P ({( x 1 , … , x m ) ∈ B }) = ν ⊗ m ( B ) を満たすもの P = ν ⊗ N P=\nu^{\otimes\mathbb{N}} P = ν ⊗ N がただ 1 つ存在する。さらにシフト σ ( x ) k = x k + 1 \sigma(x)_k=x_{k+1} σ ( x ) k = x k + 1 は可測で P ∘ σ − 1 = P P\circ\sigma^{-1}=P P ∘ σ − 1 = P を満たす。
証明(補題 2.5) 有限集合 J ⊂ N J\subset\mathbb{N} J ⊂ N ごとの積測度 μ J = ν ⊗ J \mu_J=\nu^{\otimes J} μ J = ν ⊗ J (定理 3.2[直積測度と Fubini–Tonelli の定理] )は、J ⊂ J ′ J\subset J' J ⊂ J ′ のとき μ J ′ \mu_{J'} μ J ′ の R J \mathbb{R}^J R J への周辺分布が μ J \mu_J μ J になる(余分な座標で全空間 R \mathbb{R} R を取れば ν ( R ) = 1 \nu(\mathbb{R})=1 ν ( R ) = 1 が掛かるだけ)ので整合的です。Kolmogorov の拡張定理(定理 4.1[Kolmogorov の拡張定理と連続性定理] )から、柱集合の上で μ J \mu_J μ J と一致する確率測度 P P P が F \mathcal{F} F 上にただ 1 つあります。
σ \sigma σ の可測性は、x k ∘ σ = x k + 1 x_k\circ\sigma=x_{k+1} x k ∘ σ = x k + 1 が各 k k k で可測であることから従います。柱集合 C = { ( x 1 , … , x m ) ∈ B } C=\{(x_1,\dots,x_m)\in B\} C = {( x 1 , … , x m ) ∈ B } について σ − 1 C = { ( x 2 , … , x m + 1 ) ∈ B } \sigma^{-1}C=\{(x_2,\dots,x_{m+1})\in B\} σ − 1 C = {( x 2 , … , x m + 1 ) ∈ B } で、P P P の下での座標 ( x 2 , … , x m + 1 ) (x_2,\dots,x_{m+1}) ( x 2 , … , x m + 1 ) の分布は μ { 2 , … , m + 1 } = ν ⊗ m \mu_{\{2,\dots,m+1\}}=\nu^{\otimes m} μ { 2 , … , m + 1 } = ν ⊗ m なので P ( σ − 1 C ) = ν ⊗ m ( B ) = P ( C ) P(\sigma^{-1}C)=\nu^{\otimes m}(B)=P(C) P ( σ − 1 C ) = ν ⊗ m ( B ) = P ( C ) です。柱集合の全体は代数(特に π-系)で F \mathcal{F} F を生成するので、測度の一意性(定理 4.5[測度空間と Carathéodory の拡張定理] )から P ∘ σ − 1 = P P\circ\sigma^{-1}=P P ∘ σ − 1 = P です。
∎
例 2.6 (Bernoulli シフト )
補題 2.5 の ( R N , F , ν ⊗ N , σ ) (\mathbb{R}^{\mathbb{N}},\mathcal{F},\nu^{\otimes\mathbb{N}},\sigma) ( R N , F , ν ⊗ N , σ ) を Bernoulli シフト とよびます。座標 x 1 , x 2 , … x_1,x_2,\dots x 1 , x 2 , … は P P P の下で独立で、各々の分布は ν \nu ν です。シフトは「1 回分の試行を忘れて次から数え直す」操作で、測度保存系の最も基本的な例です。ν \nu ν が 2 点 { 0 , 1 } \{0,1\} { 0 , 1 } に重み 1 2 \frac12 2 1 ずつ乗る場合は、2 倍写像と(零集合を除いて)同じものです。
測度保存性だけから、すでに次の定性的な結論が出ます。
定理 2.7 (Poincaré の回帰定理 ) ✓
( X , B , μ , T ) (X,\mathcal{B},\mu,T) ( X , B , μ , T ) を測度保存系、A ∈ B A\in\mathcal{B} A ∈ B とする。ほとんどすべての x ∈ A x\in A x ∈ A について、T n x ∈ A T^nx\in A T n x ∈ A となる n ≥ 1 n\ge1 n ≥ 1 が無限個ある。
証明(定理 2.7) B = A ∖ ⋃ n ≥ 1 T − n A B=A\setminus\bigcup_{n\ge1}T^{-n}A B = A ∖ ⋃ n ≥ 1 T − n A (A A A から出て二度と戻らない点)とおきます。0 ≤ n < m 0\le n<m 0 ≤ n < m で x ∈ T − n B ∩ T − m B x\in T^{-n}B\cap T^{-m}B x ∈ T − n B ∩ T − m B なら、y = T n x ∈ B y=T^nx\in B y = T n x ∈ B かつ T m − n y ∈ B ⊂ A T^{m-n}y\in B\subset A T m − n y ∈ B ⊂ A で、y ∈ B y\in B y ∈ B に反します。よって B , T − 1 B , T − 2 B , … B,T^{-1}B,T^{-2}B,\dots B , T − 1 B , T − 2 B , … は互いに素で、どれも測度 μ ( B ) \mu(B) μ ( B ) です。和の測度は 1 1 1 以下なので μ ( B ) = 0 \mu(B)=0 μ ( B ) = 0 です。
x ∈ A x\in A x ∈ A が有限回しか戻らないとし、最後に A A A にいた時刻を m ≥ 0 m\ge0 m ≥ 0 とすると、T m x ∈ B T^mx\in B T m x ∈ B です。したがってそのような x x x は ⋃ m ≥ 0 T − m B \bigcup_{m\ge0}T^{-m}B ⋃ m ≥ 0 T − m B に入り、これは零集合の可算和で零集合です。
∎
有限測度が本質的です。R \mathbb{R} R 上の平行移動 x ↦ x + 1 x\mapsto x+1 x ↦ x + 1 は Lebesgue 測度を保ちますが、[ 0 , 1 ) [0,1) [ 0 , 1 ) の点は二度と戻りません。
定義 3.1 (不変集合とエルゴード性 ) ↗
( X , B , μ , T ) (X,\mathcal{B},\mu,T) ( X , B , μ , T ) を測度保存系とする。T − 1 A = A T^{-1}A=A T − 1 A = A を満たす A ∈ B A\in\mathcal{B} A ∈ B を不変集合 といい、その全体を I \mathcal{I} I と書く(逆像は集合演算と可換なので I \mathcal{I} I は σ-加法族である)。すべての A ∈ I A\in\mathcal{I} A ∈ I が μ ( A ) ∈ { 0 , 1 } \mu(A)\in\{0,1\} μ ( A ) ∈ { 0 , 1 } を満たすとき、T T T (または μ \mu μ )はエルゴード的 であるという。
不変集合が測度 0 0 0 と 1 1 1 の中間にあれば、X X X は互いに行き来しない 2 つの部分に分かれ、その両方で別々に時間平均を取ることになります。エルゴード性はそれが起きないという条件です。実際に使うときは「ほとんど不変」な集合や関数で十分なことを、次の補題で保証しておきます。
この補題のねらいは、測度 0 0 0 の違いを無視してよい「ほとんど不変」な集合を、本当に不変な集合に取り替えられると示すことです。そうすると、エルゴード性を確かめるときに、厳密な不変性まで要求しなくてすみます。
補題 3.2 (ほとんど不変な集合 ) ✓
A ∈ B A\in\mathcal{B} A ∈ B が μ ( T − 1 A △ A ) = 0 \mu(T^{-1}A\,\triangle\,A)=0 μ ( T − 1 A △ A ) = 0 を満たすなら、A ∞ = ⋂ N ≥ 1 ⋃ n ≥ N T − n A A_\infty=\bigcap_{N\ge1}\bigcup_{n\ge N}T^{-n}A A ∞ = ⋂ N ≥ 1 ⋃ n ≥ N T − n A は不変集合で、μ ( A △ A ∞ ) = 0 \mu(A\,\triangle\,A_\infty)=0 μ ( A △ A ∞ ) = 0 である。
証明(補題 3.2) T − k T^{-k} T − k は対称差と可換なので、測度保存性から μ ( T − ( k + 1 ) A △ T − k A ) = μ ( T − k ( T − 1 A △ A ) ) = 0 \mu(T^{-(k+1)}A\,\triangle\,T^{-k}A)=\mu(T^{-k}(T^{-1}A\,\triangle\,A))=0 μ ( T − ( k + 1 ) A △ T − k A ) = μ ( T − k ( T − 1 A △ A )) = 0 です。三角不等式 C △ E ⊂ ( C △ D ) ∪ ( D △ E ) C\,\triangle\,E\subset(C\,\triangle\,D)\cup(D\,\triangle\,E) C △ E ⊂ ( C △ D ) ∪ ( D △ E ) を繰り返して μ ( T − n A △ A ) ≤ ∑ k < n μ ( T − ( k + 1 ) A △ T − k A ) = 0 \mu(T^{-n}A\,\triangle\,A)\le\sum_{k<n}\mu(T^{-(k+1)}A\,\triangle\,T^{-k}A)=0 μ ( T − n A △ A ) ≤ ∑ k < n μ ( T − ( k + 1 ) A △ T − k A ) = 0 です。x ∈ A ∞ ∖ A x\in A_\infty\setminus A x ∈ A ∞ ∖ A なら x ∈ T − n A ∖ A x\in T^{-n}A\setminus A x ∈ T − n A ∖ A となる n n n があり、x ∈ A ∖ A ∞ x\in A\setminus A_\infty x ∈ A ∖ A ∞ なら十分大きいすべての n n n で x ∈ A ∖ T − n A x\in A\setminus T^{-n}A x ∈ A ∖ T − n A です。どちらも ⋃ n ( T − n A △ A ) \bigcup_n(T^{-n}A\,\triangle\,A) ⋃ n ( T − n A △ A ) に入るので、μ ( A △ A ∞ ) = 0 \mu(A\,\triangle\,A_\infty)=0 μ ( A △ A ∞ ) = 0 です。最後に T − 1 A ∞ = ⋂ N ⋃ n ≥ N T − ( n + 1 ) A = A ∞ T^{-1}A_\infty=\bigcap_N\bigcup_{n\ge N}T^{-(n+1)}A=A_\infty T − 1 A ∞ = ⋂ N ⋃ n ≥ N T − ( n + 1 ) A = A ∞ です。
∎
命題 3.3 (エルゴード性の言い換え ) ✓
測度保存系について次は同値である。
T T T はエルゴード的である。
μ ( T − 1 A △ A ) = 0 \mu(T^{-1}A\,\triangle\,A)=0 μ ( T − 1 A △ A ) = 0 なる A ∈ B A\in\mathcal{B} A ∈ B はすべて μ ( A ) ∈ { 0 , 1 } \mu(A)\in\{0,1\} μ ( A ) ∈ { 0 , 1 } を満たす。
可測関数 f : X → R f:X\to\mathbb{R} f : X → R が f ∘ T = f f\circ T=f f ∘ T = f a.e. を満たせば、f f f は a.e. で定数である。
f ∈ L 2 ( μ ) f\in L^2(\mu) f ∈ L 2 ( μ ) (複素数値)が U T f = f U_Tf=f U T f = f を満たせば、f f f は a.e. で定数である。
証明(命題 3.3) 1 ⇒ 2:補題 3.2 の A ∞ A_\infty A ∞ は不変集合なので μ ( A ) = μ ( A ∞ ) ∈ { 0 , 1 } \mu(A)=\mu(A_\infty)\in\{0,1\} μ ( A ) = μ ( A ∞ ) ∈ { 0 , 1 } です。
2 ⇒ 3:c ∈ R c\in\mathbb{R} c ∈ R について T − 1 { f > c } = { f ∘ T > c } T^{-1}\{f>c\}=\{f\circ T>c\} T − 1 { f > c } = { f ∘ T > c } は { f > c } \{f>c\} { f > c } と零集合しか違わないので、μ ( f > c ) ∈ { 0 , 1 } \mu(f>c)\in\{0,1\} μ ( f > c ) ∈ { 0 , 1 } です。φ ( c ) = μ ( f > c ) \varphi(c)=\mu(f>c) φ ( c ) = μ ( f > c ) は単調減少で、f f f が実数値なので測度の連続性から c → − ∞ c\to-\infty c → − ∞ で 1 1 1 、c → + ∞ c\to+\infty c → + ∞ で 0 0 0 に近づきます。c 0 = sup { c : φ ( c ) = 1 } c_0=\sup\{c:\varphi(c)=1\} c 0 = sup { c : φ ( c ) = 1 } は有限で、各 k k k について μ ( f > c 0 + 1 k ) = 0 \mu(f>c_0+\frac1k)=0 μ ( f > c 0 + k 1 ) = 0 、μ ( f > c 0 − 1 k ) = 1 \mu(f>c_0-\frac1k)=1 μ ( f > c 0 − k 1 ) = 1 です。k k k について可算個の和・共通部分を取って μ ( f ≠ c 0 ) = 0 \mu(f\ne c_0)=0 μ ( f = c 0 ) = 0 です。
3 ⇒ 4:U T f = f U_Tf=f U T f = f は f ∘ T = f f\circ T=f f ∘ T = f a.e. のことで、実部と虚部が別々に同じ等式を満たします。
4 ⇒ 1:A ∈ I A\in\mathcal{I} A ∈ I なら 1 A ∘ T = 1 T − 1 A = 1 A \mathbf{1}_A\circ T=\mathbf{1}_{T^{-1}A}=\mathbf{1}_A 1 A ∘ T = 1 T − 1 A = 1 A なので 1 A \mathbf{1}_A 1 A は a.e. 定数、すなわち μ ( A ) ∈ { 0 , 1 } \mu(A)\in\{0,1\} μ ( A ) ∈ { 0 , 1 } です。
∎
定義 3.4 (混合 ) ↗
測度保存系が混合的 (強混合)であるとは、すべての A , B ∈ B A,B\in\mathcal{B} A , B ∈ B について
lim n → ∞ μ ( T − n A ∩ B ) = μ ( A ) μ ( B ) \lim_{n\to\infty}\mu(T^{-n}A\cap B)=\mu(A)\mu(B) n → ∞ lim μ ( T − n A ∩ B ) = μ ( A ) μ ( B ) が成り立つことをいう。
T − n A T^{-n}A T − n A は「n n n 歩後に A A A にいる」事象で、混合性は「十分時間がたつと、現在の事象 B B B と未来の事象が独立になる」ことを意味します。
証明(命題 3.5) A ∈ I A\in\mathcal{I} A ∈ I なら T − n A = A T^{-n}A=A T − n A = A なので、A = B A=B A = B として定義の式は μ ( A ) = μ ( A ) 2 \mu(A)=\mu(A)^2 μ ( A ) = μ ( A ) 2 、すなわち μ ( A ) ∈ { 0 , 1 } \mu(A)\in\{0,1\} μ ( A ) ∈ { 0 , 1 } です。
∎
混合性を確かめるには、関数の言葉に直して稠密な部分集合で調べるのが便利です。以下 L 2 ( μ ) L^2(\mu) L 2 ( μ ) の内積を ⟨ g , f ⟩ = ∫ g ˉ f d μ \langle g,f\rangle=\int\bar gf\,d\mu ⟨ g , f ⟩ = ∫ g ˉ f d μ (第 2 変数について線形)とします。
補題 3.6 (稠密な集合での混合 ) ✓
D ⊂ L 2 ( μ ) \mathcal{D}\subset L^2(\mu) D ⊂ L 2 ( μ ) を稠密な部分集合とし、すべての f , g ∈ D f,g\in\mathcal{D} f , g ∈ D について ⟨ g , U T n f ⟩ → ⟨ g , 1 ⟩ ⟨ 1 , f ⟩ \langle g,U_T^nf\rangle\to\langle g,1\rangle\langle1,f\rangle ⟨ g , U T n f ⟩ → ⟨ g , 1 ⟩ ⟨ 1 , f ⟩ (n → ∞ n\to\infty n → ∞ )とする。このとき同じことがすべての f , g ∈ L 2 ( μ ) f,g\in L^2(\mu) f , g ∈ L 2 ( μ ) で成り立ち、特に T T T は混合的である。
証明(補題 3.6) f , g ∈ L 2 f,g\in L^2 f , g ∈ L 2 と ε > 0 \varepsilon>0 ε > 0 について f ′ , g ′ ∈ D f',g'\in\mathcal{D} f ′ , g ′ ∈ D を ∥ f − f ′ ∥ 2 , ∥ g − g ′ ∥ 2 < ε \|f-f'\|_2,\|g-g'\|_2<\varepsilon ∥ f − f ′ ∥ 2 , ∥ g − g ′ ∥ 2 < ε と取ります。U T n U_T^n U T n は等長(補題 2.2 の 3)なので Cauchy–Schwarz の不等式から
∣ ⟨ g , U T n f ⟩ − ⟨ g ′ , U T n f ′ ⟩ ∣ ≤ ∥ 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 , U T n f ⟩ − ⟨ g ′ , U T n f ′ ⟩ ∣ ≤ ∥ g − g ′ ∥ 2 ∥ f ∥ 2 + ∥ g ′ ∥ 2 ∥ f − f ′ ∥ 2 ≤ ε ( ∥ f ∥ 2 + ∥ g ∥ 2 + ε ) , 右辺の ⟨ g , 1 ⟩ ⟨ 1 , f ⟩ \langle g,1\rangle\langle1,f\rangle ⟨ g , 1 ⟩ ⟨ 1 , f ⟩ も同じ程度しか変わりません。したがって lim sup n ∣ ⟨ g , U T n f ⟩ − ⟨ g , 1 ⟩ ⟨ 1 , f ⟩ ∣ \limsup_n|\langle g,U_T^nf\rangle-\langle g,1\rangle\langle1,f\rangle| lim sup n ∣ ⟨ g , U T n f ⟩ − ⟨ g , 1 ⟩ ⟨ 1 , f ⟩ ∣ は ε \varepsilon ε の定数倍以下で、ε → 0 \varepsilon\to0 ε → 0 とすれば 0 0 0 です。f = 1 A f=\mathbf{1}_A f = 1 A 、g = 1 B g=\mathbf{1}_B g = 1 B とすると ⟨ g , U T n f ⟩ = μ ( B ∩ T − n A ) \langle g,U_T^nf\rangle=\mu(B\cap T^{-n}A) ⟨ g , U T n f ⟩ = μ ( B ∩ T − n A ) 、⟨ g , 1 ⟩ ⟨ 1 , f ⟩ = μ ( B ) μ ( A ) \langle g,1\rangle\langle1,f\rangle=\mu(B)\mu(A) ⟨ g , 1 ⟩ ⟨ 1 , f ⟩ = μ ( B ) μ ( A ) です。
∎
T = [ 0 , 1 ) \mathbb{T}=[0,1) T = [ 0 , 1 ) 上で e k ( x ) = e 2 π i k x e_k(x)=e^{2\pi ikx} e k ( x ) = e 2 π ik x (k ∈ Z k\in\mathbb{Z} k ∈ Z )とおきます。Parseval の等式(定理 5.1[フーリエ級数] )を変数変換 θ = 2 π x \theta=2\pi x θ = 2 π x で移すと、{ e k } \{e_k\} { e k } は L 2 ( T , m ) L^2(\mathbb{T},m) L 2 ( T , m ) の正規直交基底で、f f f は Fourier 係数 f ^ ( k ) = ⟨ e k , f ⟩ \hat f(k)=\langle e_k,f\rangle f ^ ( k ) = ⟨ e k , f ⟩ で決まり、三角多項式(有限個の e k e_k e k の線形結合)は L 2 L^2 L 2 で稠密です。
1 と 2 は「無理数のときだけ、不変な関数が定数に限られる」ことを Fourier 係数で確かめます。3 の証明の方針は、重なりの測度 μ ( R α − n A ∩ A ) \mu(R_\alpha^{-n}A\cap A) μ ( R α − n A ∩ A ) が、ある n n n では μ ( A ) 2 = 1 4 \mu(A)^2=\frac14 μ ( A ) 2 = 4 1 から遠いままであることを示すことです。そのために、n α n\alpha n α が整数にとても近くなる n n n を鳩の巣原理で見つけます。
定理 3.7 (円の回転のエルゴード性 ) ✓
R α R_\alpha R α を 例 2.3 の回転とする。
α ∉ Q \alpha\notin\mathbb{Q} α ∈ / Q なら R α R_\alpha R α はエルゴード的である。
α ∈ Q \alpha\in\mathbb{Q} α ∈ Q なら R α R_\alpha R α はエルゴード的でない。
どの α \alpha α についても R α R_\alpha R α は混合的でない。
証明(定理 3.7) 1。命題 3.3 の 4 を確かめます。f ∈ L 2 f\in L^2 f ∈ L 2 が f ( x + α ) = f ( x ) f(x+\alpha)=f(x) f ( x + α ) = f ( x ) a.e. を満たすとすると、平行移動 y = x + α y=x+\alpha y = x + α で積分変数を変えて(Lebesgue 測度の平行移動不変性)
f ^ ( k ) = ∫ 0 1 f ( x + α ) e − 2 π i k x d x = ∫ 0 1 f ( y ) e − 2 π i k ( y − α ) d y = e 2 π i k α 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) f ^ ( k ) = ∫ 0 1 f ( x + α ) e − 2 π ik x d x = ∫ 0 1 f ( y ) e − 2 π ik ( y − α ) d y = e 2 π ik α f ^ ( k ) です。k ≠ 0 k\ne0 k = 0 なら k α ∉ Z k\alpha\notin\mathbb{Z} k α ∈ / Z なので e 2 π i k α ≠ 1 e^{2\pi ik\alpha}\ne1 e 2 π ik α = 1 、よって f ^ ( k ) = 0 \hat f(k)=0 f ^ ( k ) = 0 です。Parseval の等式から ∥ f − f ^ ( 0 ) ∥ 2 2 = ∑ k ≠ 0 ∣ f ^ ( k ) ∣ 2 = 0 \|f-\hat f(0)\|_2^2=\sum_{k\ne0}|\hat f(k)|^2=0 ∥ f − f ^ ( 0 ) ∥ 2 2 = ∑ k = 0 ∣ f ^ ( k ) ∣ 2 = 0 で、f f f は a.e. で定数です。
2。α = p / q \alpha=p/q α = p / q (q ≥ 1 q\ge1 q ≥ 1 )なら e q ( x + α ) = e 2 π i p e q ( x ) = e q ( x ) e_q(x+\alpha)=e^{2\pi ip}e_q(x)=e_q(x) e q ( x + α ) = e 2 π i p e q ( x ) = e q ( x ) で、e q e_q e q は定数でない不変関数です。
3。A = [ 0 , 1 2 ) A=[0,\frac12) A = [ 0 , 2 1 ) とすると、R α − n A = A − n α R_\alpha^{-n}A=A-n\alpha R α − n A = A − n α (mod 1)と A A A の共通部分の測度は、t = n α t=n\alpha t = n α から最も近い整数までの距離を ∥ t ∥ \|t\| ∥ t ∥ として 1 2 − ∥ n α ∥ \frac12-\|n\alpha\| 2 1 − ∥ n α ∥ です(長さ 1 2 \frac12 2 1 の 2 つの弧を ∥ t ∥ \|t\| ∥ t ∥ だけずらした重なり)。α = p / q \alpha=p/q α = p / q なら n = q , 2 q , … n=q,2q,\dots n = q , 2 q , … でこれは 1 2 \frac12 2 1 です。α \alpha α が無理数でも、N N N を固定して 0 , α , 2 α , … , N α 0,\alpha,2\alpha,\dots,N\alpha 0 , α , 2 α , … , N α を mod 1 で考えると、[ 0 , 1 ) [0,1) [ 0 , 1 ) を長さ 1 N \frac1N N 1 の N N N 個の区間に分けたとき、鳩の巣原理で 2 つが同じ区間に入り、その差から 1 ≤ q ≤ N 1\le q\le N 1 ≤ q ≤ N で ∥ q α ∥ ≤ 1 N \|q\alpha\|\le\frac1N ∥ q α ∥ ≤ N 1 となる q q q が得られます。N → ∞ N\to\infty N → ∞ で ∥ q j α ∥ → 0 \|q_j\alpha\|\to0 ∥ q j α ∥ → 0 、かつ α \alpha α が無理数なので q j → ∞ q_j\to\infty q j → ∞ と取れます(有界な q q q では ∥ q α ∥ \|q\alpha\| ∥ q α ∥ は正の最小値をもつ)。この部分列で重なりは 1 2 \frac12 2 1 に近づき、μ ( A ) 2 = 1 4 \mu(A)^2=\frac14 μ ( A ) 2 = 4 1 に収束しません。
∎
Harper 模型とほぼ Mathieu 作用素 の 補題 7.1[Harper 模型とほぼ Mathieu 作用素] は、すべての点で不変な有界関数についての同じ主張で、1 の特別な場合です。
0 1 2 3 0, 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 倍写像 D D D は Lebesgue 測度について混合的であり、特にエルゴード的である。
証明(命題 3.8) e k ∘ D ( x ) = e 2 π i k ( 2 x ) = e 2 k ( x ) e_k\circ D(x)=e^{2\pi ik(2x)}=e_{2k}(x) e k ∘ D ( x ) = e 2 π ik ( 2 x ) = e 2 k ( x ) (2 x 2x 2 x と 2 x m o d 1 2x\bmod1 2 x mod 1 は整数だけ違う)なので U D n e k = e 2 n k U_D^ne_k=e_{2^nk} U D n e k = e 2 n k です。k ≠ 0 k\ne0 k = 0 と任意の l l l について、2 n ∣ k ∣ > ∣ l ∣ 2^n|k|>|l| 2 n ∣ k ∣ > ∣ l ∣ となる n n n 以降は ⟨ e l , U D n e k ⟩ = ⟨ e l , e 2 n k ⟩ = 0 \langle e_l,U_D^ne_k\rangle=\langle e_l,e_{2^nk}\rangle=0 ⟨ e l , U D n e k ⟩ = ⟨ e l , e 2 n k ⟩ = 0 で、これは ⟨ e l , 1 ⟩ ⟨ 1 , e k ⟩ = 0 \langle e_l,1\rangle\langle1,e_k\rangle=0 ⟨ e l , 1 ⟩ ⟨ 1 , e k ⟩ = 0 に一致します。k = 0 k=0 k = 0 なら両辺とも δ l 0 \delta_{l0} δ l 0 です。双線形性から、三角多項式 f , g f,g f , g について十分大きい n n n で ⟨ g , U D n f ⟩ = ⟨ g , 1 ⟩ ⟨ 1 , f ⟩ \langle g,U_D^nf\rangle=\langle g,1\rangle\langle1,f\rangle ⟨ g , U D n f ⟩ = ⟨ g , 1 ⟩ ⟨ 1 , f ⟩ が成り立ちます。三角多項式は稠密なので 補題 3.6 から混合性が従い、命題 3.5 からエルゴード性が従います。
∎
Bernoulli シフトでは Fourier 展開の代わりに、柱集合で一般の事象を近似します。
補題 3.9 (代数による近似 ) ✓
A \mathcal{A} A を集合の代数、μ \mu μ を σ ( A ) \sigma(\mathcal{A}) σ ( A ) 上の確率測度とする。任意の B ∈ σ ( A ) B\in\sigma(\mathcal{A}) B ∈ σ ( A ) と ε > 0 \varepsilon>0 ε > 0 について、μ ( A △ B ) < ε \mu(A\,\triangle\,B)<\varepsilon μ ( A △ B ) < ε となる A ∈ A A\in\mathcal{A} A ∈ A が存在する。
証明(補題 3.9) 主張が成り立つ B B B の全体を G \mathcal{G} G とします。A ⊂ G \mathcal{A}\subset\mathcal{G} A ⊂ G です。A c △ B c = A △ B A^c\,\triangle\,B^c=A\,\triangle\,B A c △ B c = A △ B なので G \mathcal{G} G は補集合で閉じています。B 1 , B 2 , ⋯ ∈ G B_1,B_2,\dots\in\mathcal{G} B 1 , B 2 , ⋯ ∈ G 、B = ⋃ i B i B=\bigcup_iB_i B = ⋃ i B i とします。測度の連続性から μ ( B ∖ ⋃ i ≤ N B i ) < ε / 2 \mu(B\setminus\bigcup_{i\le N}B_i)<\varepsilon/2 μ ( B ∖ ⋃ i ≤ N B i ) < ε /2 となる N N N があり、A i ∈ A A_i\in\mathcal{A} A i ∈ A を μ ( A i △ B i ) < ε / 2 i + 1 \mu(A_i\,\triangle\,B_i)<\varepsilon/2^{i+1} μ ( A i △ B i ) < ε / 2 i + 1 と取って A = ⋃ i ≤ N A i ∈ A A=\bigcup_{i\le N}A_i\in\mathcal{A} A = ⋃ i ≤ N A i ∈ A とすると
A △ B ⊂ ( B ∖ ⋃ i ≤ N B i ) ∪ ⋃ i ≤ N ( A i △ B i ) 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 ⊂ ( B ∖ i ≤ N ⋃ B i ) ∪ i ≤ N ⋃ ( A i △ B i ) なので μ ( A △ B ) < ε \mu(A\,\triangle\,B)<\varepsilon μ ( A △ B ) < ε です。よって G \mathcal{G} G は A \mathcal{A} A を含む σ-加法族で、σ ( A ) ⊂ G \sigma(\mathcal{A})\subset\mathcal{G} σ ( A ) ⊂ G です。
∎
命題 3.10 (Bernoulli シフトは混合的 ) ✓
補題 2.5 の Bernoulli シフト ( R N , F , ν ⊗ N , σ ) (\mathbb{R}^{\mathbb{N}},\mathcal{F},\nu^{\otimes\mathbb{N}},\sigma) ( R N , F , ν ⊗ N , σ ) は混合的であり、特にエルゴード的である。
証明(命題 3.10) A , B A,B A , B が座標 1 , … , m 1,\dots,m 1 , … , m で決まる柱集合、A = { ( x 1 , … , x m ) ∈ A ′ } A=\{(x_1,\dots,x_m)\in A'\} A = {( x 1 , … , x m ) ∈ A ′ } 、B = { ( x 1 , … , x m ) ∈ B ′ } B=\{(x_1,\dots,x_m)\in B'\} B = {( x 1 , … , x m ) ∈ B ′ } のとき、n ≥ m n\ge m n ≥ m なら σ − n A = { ( x n + 1 , … , x n + m ) ∈ A ′ } \sigma^{-n}A=\{(x_{n+1},\dots,x_{n+m})\in A'\} σ − n A = {( x n + 1 , … , x n + m ) ∈ A ′ } は B B B と重ならない座標で決まるので、座標 1 , … , n + m 1,\dots,n+m 1 , … , n + m の分布が ν ⊗ ( n + m ) \nu^{\otimes(n+m)} ν ⊗ ( n + m ) であること(補題 2.5 )から P ( σ − n A ∩ 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) P ( σ − n A ∩ B ) = ν ⊗ m ( A ′ ) ν ⊗ m ( B ′ ) = P ( A ) P ( B ) です。一般の A , B ∈ F A,B\in\mathcal{F} A , B ∈ F と ε > 0 \varepsilon>0 ε > 0 について、補題 3.9 で柱集合 A 0 , B 0 A_0,B_0 A 0 , B 0 を P ( A △ A 0 ) , P ( B △ B 0 ) < ε P(A\,\triangle\,A_0),P(B\,\triangle\,B_0)<\varepsilon P ( A △ A 0 ) , P ( B △ B 0 ) < ε と取ります(座標の個数は大きい方にそろえる)。P ( σ − n ( A △ A 0 ) ) = P ( A △ A 0 ) P(\sigma^{-n}(A\,\triangle\,A_0))=P(A\,\triangle\,A_0) P ( σ − n ( A △ A 0 )) = P ( A △ A 0 ) なので
∣ P ( σ − n A ∩ B ) − P ( σ − n A 0 ∩ B 0 ) ∣ ≤ 2 ε , ∣ P ( A ) P ( B ) − P ( A 0 ) P ( B 0 ) ∣ ≤ 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 ∣ P ( σ − n A ∩ B ) − P ( σ − n A 0 ∩ B 0 ) ∣ ≤ 2 ε , ∣ P ( A ) P ( B ) − P ( A 0 ) P ( B 0 ) ∣ ≤ 2 ε で、lim sup n ∣ P ( σ − n A ∩ B ) − P ( A ) P ( B ) ∣ ≤ 4 ε \limsup_n|P(\sigma^{-n}A\cap B)-P(A)P(B)|\le4\varepsilon lim sup n ∣ P ( σ − n A ∩ B ) − P ( A ) P ( B ) ∣ ≤ 4 ε です。
∎
まず L 2 L^2 L 2 での収束を、Hilbert 空間の一般論として示します。V V V を Hilbert 空間 H \mathcal{H} H 上の有界作用素で ∥ V ∥ ≤ 1 \|V\|\le1 ∥ V ∥ ≤ 1 を満たすもの(縮小作用素 )とします。Koopman 作用素 U T U_T U T は等長なので縮小作用素です。
この補題は、次の定理の証明の準備です。空間を「動かない部分」と「x − V x x-Vx x − V x の形のベクトルで張られる部分」に直交分解し、前者では平均が変わらず、後者では平均が打ち消し合って 0 0 0 に近づく、という見通しを与えます。
補題 4.1 (縮小作用素の不動点 ) ✓
∥ V ∥ ≤ 1 \|V\|\le1 ∥ V ∥ ≤ 1 なら、V x = x Vx=x V x = x と V ∗ x = x V^*x=x V ∗ x = x は同値である。さらに
H = ker ( I − V ) ⊕ R a n ( I − V ) ‾ (直交直和) \mathcal{H}=\ker(I-V)\oplus\overline{\mathrm{Ran}(I-V)}\qquad\text{(直交直和)} H = ker ( I − V ) ⊕ Ran ( I − V ) (直交直和) が成り立つ。
証明(補題 4.1) V x = x Vx=x V x = x とします。∥ V ∗ ∥ = ∥ V ∥ ≤ 1 \|V^*\|=\|V\|\le1 ∥ V ∗ ∥ = ∥ V ∥ ≤ 1 (命題 3.3[有界作用素と共役作用素] )と ⟨ x , V ∗ x ⟩ = ⟨ V x , x ⟩ = ∥ x ∥ 2 \langle x,V^*x\rangle=\langle Vx,x\rangle=\|x\|^2 ⟨ x , V ∗ x ⟩ = ⟨ V x , x ⟩ = ∥ x ∥ 2 から
∥ V ∗ x − x ∥ 2 = ∥ V ∗ x ∥ 2 − 2 R e ⟨ 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 ∗ x − x ∥ 2 = ∥ V ∗ x ∥ 2 − 2 Re ⟨ x , V ∗ x ⟩ + ∥ x ∥ 2 ≤ ∥ x ∥ 2 − 2∥ x ∥ 2 + ∥ x ∥ 2 = 0 です。V ∗ ∗ = V V^{**}=V V ∗∗ = V なので逆も同じです。次に、y ⊥ R a n ( I − V ) y\perp\mathrm{Ran}(I-V) y ⊥ Ran ( I − V ) はすべての x x x で ⟨ ( I − V ∗ ) y , x ⟩ = ⟨ y , ( I − V ) x ⟩ = 0 \langle(I-V^*)y,x\rangle=\langle y,(I-V)x\rangle=0 ⟨( I − V ∗ ) y , x ⟩ = ⟨ y , ( I − V ) x ⟩ = 0 となること、すなわち y ∈ ker ( I − V ∗ ) = ker ( I − V ) y\in\ker(I-V^*)=\ker(I-V) y ∈ ker ( I − V ∗ ) = ker ( I − V ) と同値です。M = R a n ( I − V ) ‾ M=\overline{\mathrm{Ran}(I-V)} M = Ran ( I − V ) は閉部分空間で M ⊥ = R a n ( I − V ) ⊥ = ker ( I − V ) M^\perp=\mathrm{Ran}(I-V)^\perp=\ker(I-V) M ⊥ = Ran ( I − V ) ⊥ = ker ( I − V ) なので、射影定理(定理 3.4[Hilbert 空間の基礎] の 2)から主張の分解が得られます。
∎
定理 4.2 (von Neumann の平均エルゴード定理 ) ✓
V V V を Hilbert 空間 H \mathcal{H} H 上の縮小作用素、P P P を閉部分空間 ker ( I − V ) \ker(I-V) ker ( I − V ) への直交射影(定義 3.5[Hilbert 空間の基礎] )とする。すべての x ∈ H x\in\mathcal{H} x ∈ H について
∥ 1 n ∑ k = 0 n − 1 V k x − P x ∥ → 0 ( n → ∞ ) . \Bigl\|\frac1n\sum_{k=0}^{n-1}V^kx-Px\Bigr\|\to0\qquad(n\to\infty). n 1 k = 0 ∑ n − 1 V k x − P x → 0 ( n → ∞ ) .
証明(定理 4.2) M n = 1 n ∑ k < n V k M_n=\frac1n\sum_{k<n}V^k M n = n 1 ∑ k < n V k とおくと ∥ M n ∥ ≤ 1 \|M_n\|\le1 ∥ M n ∥ ≤ 1 です。補題 4.1 の分解で x = P x + ( x − P x ) x=Px+(x-Px) x = P x + ( x − P x ) 、x − P x ∈ R a n ( I − V ) ‾ x-Px\in\overline{\mathrm{Ran}(I-V)} x − P x ∈ Ran ( I − V ) と書きます。P x Px P x は V V V で動かないので M n P x = P x M_nPx=Px M n P x = P x です。z = y − V y ∈ R a n ( I − V ) z=y-Vy\in\mathrm{Ran}(I-V) z = y − V y ∈ Ran ( I − V ) なら和が打ち消し合って
M n z = 1 n ( y − V n y ) , ∥ M n z ∥ ≤ 2 ∥ y ∥ n → 0 M_nz=\frac1n(y-V^ny),\qquad\|M_nz\|\le\frac{2\|y\|}n\to0 M n z = n 1 ( y − V n y ) , ∥ M n z ∥ ≤ n 2∥ y ∥ → 0 です。w ∈ R a n ( I − V ) ‾ w\in\overline{\mathrm{Ran}(I-V)} w ∈ Ran ( I − V ) と ε > 0 \varepsilon>0 ε > 0 について ∥ w − z ∥ < ε \|w-z\|<\varepsilon ∥ w − z ∥ < ε となる z ∈ R a n ( I − V ) z\in\mathrm{Ran}(I-V) z ∈ Ran ( I − V ) を取れば、∥ M n w ∥ ≤ ∥ M n ( w − z ) ∥ + ∥ M n z ∥ ≤ ε + ∥ M n z ∥ \|M_nw\|\le\|M_n(w-z)\|+\|M_nz\|\le\varepsilon+\|M_nz\| ∥ M n w ∥ ≤ ∥ M n ( w − z ) ∥ + ∥ M n z ∥ ≤ ε + ∥ M n z ∥ で、lim sup n ∥ M n w ∥ ≤ ε \limsup_n\|M_nw\|\le\varepsilon lim sup n ∥ M n w ∥ ≤ ε です。w = x − P x w=x-Px w = x − P x として M n x − P x = M n ( x − P x ) → 0 M_nx-Px=M_n(x-Px)\to0 M n x − P x = M n ( x − P x ) → 0 を得ます。
∎
系 4.3 (測度保存系での平均エルゴード定理 ) ✓
測度保存系と f ∈ L 2 ( μ ) f\in L^2(\mu) f ∈ L 2 ( μ ) について、A n f A_nf A n f は L 2 ( μ ) L^2(\mu) L 2 ( μ ) で U T U_T U T の不動点空間への射影 P f Pf P f に収束する。T T T がエルゴード的なら P f = ∫ f d μ Pf=\int f\,d\mu P f = ∫ f d μ (定数関数)である。
証明(系 4.3) A n f = 1 n ∑ k < n U T k f A_nf=\frac1n\sum_{k<n}U_T^kf A n f = n 1 ∑ k < n U T k f なので、前半は 定理 4.2 を V = U T V=U_T V = U T に使ったものです。エルゴード的なら 命題 3.3 の 4 から ker ( I − U T ) \ker(I-U_T) ker ( I − U T ) は定数関数の空間 C 1 \mathbb{C}1 C 1 で、∥ 1 ∥ 2 = 1 \|1\|_2=1 ∥1 ∥ 2 = 1 なので射影は P f = ⟨ 1 , f ⟩ 1 = ∫ f d μ Pf=\langle1,f\rangle1=\int f\,d\mu P f = ⟨ 1 , f ⟩ 1 = ∫ f d μ です。
∎
例 4.4 (回転での平均の計算 )
α ∉ Q \alpha\notin\mathbb{Q} α ∈ / Q 、k ≠ 0 k\ne0 k = 0 とすると、U R α e k = e 2 π i k α e k U_{R_\alpha}e_k=e^{2\pi ik\alpha}e_k U R α e k = e 2 π ik α e k なので、等比級数の和から
A n e k = 1 n ∑ j = 0 n − 1 e 2 π i j k α e k = 1 − e 2 π i n k α n ( 1 − e 2 π i k α ) e k , ∥ A n e k ∥ ∞ ≤ 2 n ∣ 1 − e 2 π i k α ∣ 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}|} A n e k = n 1 j = 0 ∑ n − 1 e 2 π ij k α e k = n ( 1 − e 2 π ik α ) 1 − e 2 π ink α e k , ∥ A n e k ∥ ∞ ≤ n ∣1 − e 2 π ik α ∣ 2 です。収束は一様ですが、その速さ 1 / ( n ∣ 1 − e 2 π i k α ∣ ) 1/(n|1-e^{2\pi ik\alpha}|) 1/ ( n ∣1 − e 2 π ik α ∣ ) は k α k\alpha k α が整数にどれだけ近いかで決まり、k k k について一様ではありません。平均エルゴード定理は収束を保証しますが、一般には速さについて何も言いません。
平均エルゴード定理は L 2 L^2 L 2 の収束で、個々の点 x x x での時間平均については何も言いません。各点の収束には、時間平均の「最大値」を制御する不等式が要ります。以下 f f f は実数値とします(複素数値なら実部と虚部に分けます)。
補題 5.1 (最大エルゴード補題 ) ✓
( X , B , μ , T ) (X,\mathcal{B},\mu,T) ( X , B , μ , T ) を測度保存系、f ∈ L 1 ( μ ) f\in L^1(\mu) f ∈ L 1 ( μ ) を実数値とし、N ≥ 1 N\ge1 N ≥ 1 について M N = max 0 ≤ n ≤ N S n f M_N=\max_{0\le n\le N}S_nf M N = max 0 ≤ n ≤ N S n f (S 0 f = 0 S_0f=0 S 0 f = 0 を含むので M N ≥ 0 M_N\ge0 M N ≥ 0 )とおく。このとき
∫ { M N > 0 } f d μ ≥ 0 , ∫ { sup n ≥ 1 S n f > 0 } f d μ ≥ 0. \int_{\{M_N>0\}}f\,d\mu\ge0,\qquad\int_{\{\sup_{n\ge1}S_nf>0\}}f\,d\mu\ge0 . ∫ { M N > 0 } f d μ ≥ 0 , ∫ { s u p n ≥ 1 S n f > 0 } f d μ ≥ 0.
証明(補題 5.1) (Garsia の証明)0 ≤ n ≤ N 0\le n\le N 0 ≤ n ≤ N について M N ∘ T ≥ S n f ∘ T M_N\circ T\ge S_nf\circ T M N ∘ T ≥ S n f ∘ T なので、f + M N ∘ T ≥ f + S n f ∘ T = S n + 1 f f+M_N\circ T\ge f+S_nf\circ T=S_{n+1}f f + M N ∘ T ≥ f + S n f ∘ T = S n + 1 f です。n n n について最大を取ると
f + M N ∘ T ≥ max 1 ≤ n ≤ N + 1 S n f ≥ max 1 ≤ n ≤ N S n f . f+M_N\circ T\ge\max_{1\le n\le N+1}S_nf\ge\max_{1\le n\le N}S_nf . f + M N ∘ T ≥ 1 ≤ n ≤ N + 1 max S n f ≥ 1 ≤ n ≤ N max S n f . { M N > 0 } \{M_N>0\} { M N > 0 } の上では最大値は n = 0 n=0 n = 0 の項 0 0 0 ではないので M N = max 1 ≤ n ≤ N S n f M_N=\max_{1\le n\le N}S_nf M N = max 1 ≤ n ≤ N S n f で、したがって f ≥ M N − M N ∘ T f\ge M_N-M_N\circ T f ≥ M N − M N ∘ T です。∣ M N ∣ ≤ ∑ k < N ∣ f ∣ ∘ T k |M_N|\le\sum_{k<N}|f|\circ T^k ∣ M N ∣ ≤ ∑ k < N ∣ f ∣ ∘ T k は可積分(補題 2.2 の 2)で、{ M N > 0 } \{M_N>0\} { M N > 0 } の外では M N = 0 M_N=0 M N = 0 、また M N ∘ T ≥ 0 M_N\circ T\ge0 M N ∘ T ≥ 0 なので
∫ { M N > 0 } f d μ ≥ ∫ { M N > 0 } ( M N − M N ∘ T ) d μ ≥ ∫ X M N d μ − ∫ X M N ∘ 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 ∫ { M N > 0 } f d μ ≥ ∫ { M N > 0 } ( M N − M N ∘ T ) d μ ≥ ∫ X M N d μ − ∫ X M N ∘ T d μ = 0 です(最後は 補題 2.2 の 2)。N → ∞ N\to\infty N → ∞ で { M N > 0 } \{M_N>0\} { M N > 0 } は { sup n ≥ 1 S n f > 0 } \{\sup_{n\ge1}S_nf>0\} { sup n ≥ 1 S n f > 0 } に増加し、∣ f 1 { M N > 0 } ∣ ≤ ∣ f ∣ |f\mathbf{1}_{\{M_N>0\}}|\le|f| ∣ f 1 { M N > 0 } ∣ ≤ ∣ f ∣ なので、優収束定理(定理 7.3[ルベーグ積分の定義と収束定理] )で 2 つ目の不等式を得ます。
∎
この補題の気持ちは「時間平均がどこかで正になる点の集合の上では、f f f の平均は非負」ということです。次の定理ではこれを f − b f-b f − b や a − f a-f a − f に使い、時間平均の上極限と下極限が食い違う集合の測度を 0 0 0 に押さえ込みます。
証明の方針は次のとおりです。時間平均の上極限 f ∗ f^* f ∗ と下極限 f ∗ f_* f ∗ が食い違う点の集合は不変集合になります。その上で f f f から定数 b b b を引く(または a a a から f f f を引く)と、その集合の上では時間平均が正になる点がある状況をつくれるので、最大エルゴード補題から積分の符号が決まります。2 つを足すと、食い違う集合の測度が 0 0 0 だと分かります。
定理 5.2 (Birkhoff の各点エルゴード定理 ) ✓
( X , B , μ , T ) (X,\mathcal{B},\mu,T) ( X , B , μ , T ) を測度保存系、I \mathcal{I} I を不変 σ-加法族、f ∈ L 1 ( μ ) f\in L^1(\mu) f ∈ L 1 ( μ ) とする。I \mathcal{I} I -可測な f ˉ ∈ L 1 ( μ ) \bar f\in L^1(\mu) f ˉ ∈ L 1 ( μ ) が存在して次が成り立つ。
ほとんどすべての x x x で A n f ( x ) = 1 n ∑ k = 0 n − 1 f ( T k x ) → f ˉ ( x ) A_nf(x)=\frac1n\sum_{k=0}^{n-1}f(T^kx)\to\bar f(x) A n f ( x ) = n 1 ∑ k = 0 n − 1 f ( T k x ) → f ˉ ( x ) 。
∥ A n f − f ˉ ∥ 1 → 0 \|A_nf-\bar f\|_1\to0 ∥ A n f − f ˉ ∥ 1 → 0 、∥ f ˉ ∥ 1 ≤ ∥ f ∥ 1 \|\bar f\|_1\le\|f\|_1 ∥ f ˉ ∥ 1 ≤ ∥ f ∥ 1 。
すべての B ∈ I B\in\mathcal{I} B ∈ I で ∫ B f ˉ d μ = ∫ B f d μ \int_B\bar f\,d\mu=\int_Bf\,d\mu ∫ B f ˉ d μ = ∫ B f d μ 。すなわち f ˉ = E [ f ∣ I ] \bar f=E[f\mid\mathcal{I}] f ˉ = E [ f ∣ I ] である。
T T T がエルゴード的なら、f ˉ = ∫ f d μ \bar f=\int f\,d\mu f ˉ = ∫ f d μ a.e. である。
証明(定理 5.2) f f f は実数値とし、代表元をすべての点で有限な値をとるように選びます。
不変性。 f ∗ = lim sup n A n f f^*=\limsup_nA_nf f ∗ = lim sup n A n f 、f ∗ = lim inf n A n f f_*=\liminf_nA_nf f ∗ = lim inf n A n f (値は [ − ∞ , ∞ ] [-\infty,\infty] [ − ∞ , ∞ ] )とおきます。( n + 1 ) A n + 1 f ( x ) = f ( x ) + n A n f ( T x ) (n+1)A_{n+1}f(x)=f(x)+nA_nf(Tx) ( n + 1 ) A n + 1 f ( x ) = f ( x ) + n A n f ( T x ) から
A n + 1 f ( x ) = f ( x ) n + 1 + n n + 1 A n f ( T x ) A_{n+1}f(x)=\frac{f(x)}{n+1}+\frac n{n+1}A_nf(Tx) A n + 1 f ( x ) = n + 1 f ( x ) + n + 1 n A n f ( T x ) で、f ( x ) f(x) f ( x ) は有限なので、すべての x x x で f ∗ ( T x ) = f ∗ ( x ) f^*(Tx)=f^*(x) f ∗ ( T x ) = f ∗ ( x ) 、f ∗ ( T x ) = f ∗ ( x ) f_*(Tx)=f_*(x) f ∗ ( T x ) = f ∗ ( x ) です。したがって { f ∗ > c } \{f^*>c\} { f ∗ > c } などは不変集合で、f ∗ , f ∗ f^*,f_* f ∗ , f ∗ は I \mathcal{I} I -可測です。
上極限と下極限の一致。 有理数 a < b a<b a < b について E = { f ∗ < a } ∩ { f ∗ > b } ∈ I E=\{f_*<a\}\cap\{f^*>b\}\in\mathcal{I} E = { f ∗ < a } ∩ { f ∗ > b } ∈ I とおきます。g = ( f − b ) 1 E g=(f-b)\mathbf{1}_E g = ( f − b ) 1 E とすると、1 E ∘ T k = 1 T − k E = 1 E \mathbf{1}_E\circ T^k=\mathbf{1}_{T^{-k}E}=\mathbf{1}_E 1 E ∘ T k = 1 T − k E = 1 E なので S n g = 1 E ( S n f − n b ) = n 1 E ( A n f − b ) S_ng=\mathbf{1}_E\,(S_nf-nb)=n\mathbf{1}_E(A_nf-b) S n g = 1 E ( S n f − nb ) = n 1 E ( A n f − b ) です。E E E の上では lim sup A n f > b \limsup A_nf>b lim sup A n f > b なので sup n ≥ 1 S n g > 0 \sup_{n\ge1}S_ng>0 sup n ≥ 1 S n g > 0 、E E E の外では S n g = 0 S_ng=0 S n g = 0 です。よって { sup n ≥ 1 S n g > 0 } = E \{\sup_{n\ge1}S_ng>0\}=E { sup n ≥ 1 S n g > 0 } = E で、補題 5.1 から ∫ E ( f − b ) d μ ≥ 0 \int_E(f-b)\,d\mu\ge0 ∫ E ( f − b ) d μ ≥ 0 です。同様に h = ( a − f ) 1 E h=(a-f)\mathbf{1}_E h = ( a − f ) 1 E について、E E E の上では lim inf A n f < a \liminf A_nf<a lim inf A n f < a なので { sup n ≥ 1 S n h > 0 } = E \{\sup_{n\ge1}S_nh>0\}=E { sup n ≥ 1 S n h > 0 } = E で、∫ E ( a − f ) d μ ≥ 0 \int_E(a-f)\,d\mu\ge0 ∫ E ( a − f ) d μ ≥ 0 です。足すと ( a − b ) μ ( E ) ≥ 0 (a-b)\mu(E)\ge0 ( a − b ) μ ( E ) ≥ 0 、すなわち μ ( E ) = 0 \mu(E)=0 μ ( E ) = 0 です。a < b a<b a < b を有理数全体で動かした可算和も零集合なので、ほとんどすべての点で f ∗ = f ∗ f_*=f^* f ∗ = f ∗ 、すなわち A n f A_nf A n f は [ − ∞ , ∞ ] [-\infty,\infty] [ − ∞ , ∞ ] の中で収束します。
可積分性。 ∣ A n f ∣ ≤ A n ∣ f ∣ |A_nf|\le A_n|f| ∣ A n f ∣ ≤ A n ∣ f ∣ と 補題 2.2 から ∫ ∣ A n f ∣ d μ ≤ ∥ f ∥ 1 \int|A_nf|\,d\mu\le\|f\|_1 ∫ ∣ A n f ∣ d μ ≤ ∥ f ∥ 1 です。Fatou の補題(補題 6.1[ルベーグ積分の定義と収束定理] )で ∫ ∣ f ∗ ∣ d μ ≤ lim inf n ∫ ∣ A n f ∣ d μ ≤ ∥ f ∥ 1 \int|f^*|\,d\mu\le\liminf_n\int|A_nf|\,d\mu\le\|f\|_1 ∫ ∣ f ∗ ∣ d μ ≤ lim inf n ∫ ∣ A n f ∣ d μ ≤ ∥ f ∥ 1 なので、f ∗ f^* f ∗ は a.e. で有限です。f ˉ = f ∗ 1 { ∣ f ∗ ∣ < ∞ } \bar f=f^*\mathbf{1}_{\{|f^*|<\infty\}} f ˉ = f ∗ 1 { ∣ f ∗ ∣ < ∞ } とおくと({ ∣ f ∗ ∣ < ∞ } ∈ I \{|f^*|<\infty\}\in\mathcal{I} { ∣ f ∗ ∣ < ∞ } ∈ I )、f ˉ \bar f f ˉ は I \mathcal{I} I -可測、可積分で、1 と ∥ f ˉ ∥ 1 ≤ ∥ f ∥ 1 \|\bar f\|_1\le\|f\|_1 ∥ f ˉ ∥ 1 ≤ ∥ f ∥ 1 が成り立ちます。
L 1 L^1 L 1 収束。 ∣ f ∣ ≤ K |f|\le K ∣ f ∣ ≤ K なら ∣ A n f ∣ ≤ K |A_nf|\le K ∣ A n f ∣ ≤ K なので、優収束定理(定理 7.3[ルベーグ積分の定義と収束定理] )で ∥ A n f − f ˉ ∥ 1 → 0 \|A_nf-\bar f\|_1\to0 ∥ A n f − f ˉ ∥ 1 → 0 です。一般の f f f では f K = max ( − K , min ( f , K ) ) f_K=\max(-K,\min(f,K)) f K = max ( − K , min ( f , K )) とおくと、優収束定理で ∥ f − f K ∥ 1 → 0 \|f-f_K\|_1\to0 ∥ f − f K ∥ 1 → 0 (K → ∞ K\to\infty K → ∞ )です。A n A_n A n は線形なので a.e. で f ˉ − f K ‾ = lim n A n ( f − f K ) \bar f-\overline{f_K}=\lim_nA_n(f-f_K) f ˉ − f K = lim n A n ( f − f K ) で、上の可積分性の議論を f − f K f-f_K f − f K に使えば ∥ f ˉ − f K ‾ ∥ 1 ≤ ∥ f − f K ∥ 1 \|\bar f-\overline{f_K}\|_1\le\|f-f_K\|_1 ∥ f ˉ − f K ∥ 1 ≤ ∥ f − f K ∥ 1 です。したがって
∥ A n f − f ˉ ∥ 1 ≤ ∥ A n ( f − f K ) ∥ 1 + ∥ A n f K − f K ‾ ∥ 1 + ∥ f K ‾ − f ˉ ∥ 1 ≤ 2 ∥ f − f K ∥ 1 + ∥ A n f K − f K ‾ ∥ 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 ∥ A n f − f ˉ ∥ 1 ≤ ∥ A n ( f − f K ) ∥ 1 + ∥ A n f K − f K ∥ 1 + ∥ f K − f ˉ ∥ 1 ≤ 2∥ f − f K ∥ 1 + ∥ A n f K − f K ∥ 1 で、K K K を大きく取ってから n → ∞ n\to\infty n → ∞ とすれば 0 0 0 に近づきます。
条件付き期待値。 B ∈ I B\in\mathcal{I} B ∈ I なら 1 B = 1 B ∘ T k \mathbf{1}_B=\mathbf{1}_B\circ T^k 1 B = 1 B ∘ T k なので、補題 2.2 から ∫ B f ∘ T k d μ = ∫ ( 1 B f ) ∘ T k d μ = ∫ B f d μ \int_Bf\circ T^k\,d\mu=\int(\mathbf{1}_Bf)\circ T^k\,d\mu=\int_Bf\,d\mu ∫ B f ∘ T k d μ = ∫ ( 1 B f ) ∘ T k d μ = ∫ B f d μ 、よって ∫ B A n f d μ = ∫ B f d μ \int_BA_nf\,d\mu=\int_Bf\,d\mu ∫ B A n f d μ = ∫ B f d μ です。L 1 L^1 L 1 収束から ∫ B f ˉ d μ = ∫ B f d μ \int_B\bar f\,d\mu=\int_Bf\,d\mu ∫ B f ˉ d μ = ∫ B f d μ です。f ˉ \bar f f ˉ は I \mathcal{I} I -可測で可積分なので、これは条件付き期待値の定義(定義 3.1[条件付き期待値] )の 3 条件で、一意性(定理 3.2[条件付き期待値] )から f ˉ = E [ f ∣ I ] \bar f=E[f\mid\mathcal{I}] f ˉ = E [ f ∣ I ] です。
エルゴード的な場合。 f ˉ \bar f f ˉ は I \mathcal{I} I -可測なので f ˉ ∘ T = f ˉ \bar f\circ T=\bar f f ˉ ∘ T = f ˉ で、命題 3.3 の 3 から a.e. で定数 c c c です。3 で B = X B=X B = X として c = ∫ f d μ c=\int f\,d\mu c = ∫ f d μ です。
∎
f ∈ L 2 f\in L^2 f ∈ L 2 なら、系 4.3 の L 2 L^2 L 2 極限 P f Pf P f と f ˉ \bar f f ˉ は一致します。確率空間では ∥ ⋅ ∥ 1 ≤ ∥ ⋅ ∥ 2 \|\cdot\|_1\le\|\cdot\|_2 ∥ ⋅ ∥ 1 ≤ ∥ ⋅ ∥ 2 (Cauchy–Schwarz)なので A n f A_nf A n f は L 1 L^1 L 1 で P f Pf P f にも f ˉ \bar f f ˉ にも収束し、L 1 L^1 L 1 の極限は一意だからです。特に P P P は L 2 L^2 L 2 上の条件付き期待値 E [ ⋅ ∣ I ] E[\cdot\mid\mathcal{I}] E [ ⋅ ∣ I ] で、これは 定理 5.1[条件付き期待値] と整合します。
系 5.4 (大数の強法則(エルゴード定理による証明) ) ✓
X 1 , X 2 , … X_1,X_2,\dots X 1 , X 2 , … を確率空間 ( Ω ′ , F ′ , P ′ ) (\Omega',\mathcal{F}',P') ( Ω ′ , F ′ , P ′ ) 上の独立同分布な確率変数列で E [ ∣ X 1 ∣ ] < ∞ E[|X_1|]<\infty E [ ∣ X 1 ∣ ] < ∞ とする。このとき P ′ P' P ′ -a.s. で 1 n ∑ k = 1 n X k → E [ X 1 ] \frac1n\sum_{k=1}^nX_k\to E[X_1] n 1 ∑ k = 1 n X k → E [ X 1 ] である。
証明(系 5.4) ν \nu ν を X 1 X_1 X 1 の分布、( R N , F , P , σ ) (\mathbb{R}^{\mathbb{N}},\mathcal{F},P,\sigma) ( R N , F , P , σ ) を 補題 2.5 の Bernoulli シフトとします。Φ ( ω ) = ( X 1 ( ω ) , X 2 ( ω ) , … ) \Phi(\omega)=(X_1(\omega),X_2(\omega),\dots) Φ ( ω ) = ( X 1 ( ω ) , X 2 ( ω ) , … ) は各座標が可測なので、座標写像で生成される直積 σ-加法族 F \mathcal{F} F について可測です。独立性から ( X 1 , … , X m ) (X_1,\dots,X_m) ( X 1 , … , X m ) の分布は直方体の上で ν ⊗ m \nu^{\otimes m} ν ⊗ m と一致し、直方体は π-系なので測度の一意性(定理 4.5[測度空間と Carathéodory の拡張定理] )から分布は ν ⊗ m \nu^{\otimes m} ν ⊗ m です。よって Q = P ′ ∘ Φ − 1 Q=P'\circ\Phi^{-1} Q = P ′ ∘ Φ − 1 は柱集合の上で P P P と一致し、同じ一意性から Q = P Q=P Q = P です。
f ( x ) = x 1 f(x)=x_1 f ( x ) = x 1 は ∫ ∣ f ∣ d P = E [ ∣ X 1 ∣ ] < ∞ \int|f|\,dP=E[|X_1|]<\infty ∫ ∣ f ∣ d P = E [ ∣ X 1 ∣ ] < ∞ を満たし、f ∘ σ k ( x ) = x k + 1 f\circ\sigma^k(x)=x_{k+1} f ∘ σ k ( x ) = x k + 1 なので A n f ( x ) = 1 n ∑ k = 1 n x k A_nf(x)=\frac1n\sum_{k=1}^nx_k A n f ( x ) = n 1 ∑ k = 1 n x k です。命題 3.10 からシフトはエルゴード的なので、定理 5.2 の 4 により、集合 G = { x : 1 n ∑ k = 1 n x k → ∫ x 1 d P } G=\{x:\frac1n\sum_{k=1}^nx_k\to\int x_1\,dP\} G = { x : n 1 ∑ k = 1 n x k → ∫ x 1 d P } (各点の上極限・下極限が可測なので G ∈ F G\in\mathcal{F} G ∈ F )は P ( G ) = 1 P(G)=1 P ( G ) = 1 です。∫ x 1 d P = E [ X 1 ] \int x_1\,dP=E[X_1] ∫ x 1 d P = E [ X 1 ] で、P ′ ( Φ − 1 G ) = Q ( G ) = P ( G ) = 1 P'(\Phi^{-1}G)=Q(G)=P(G)=1 P ′ ( Φ − 1 G ) = Q ( G ) = P ( G ) = 1 です。
∎
大数の法則と中心極限定理 の 定理 4.5[大数の法則と中心極限定理] の証明は独立性を分散の計算に使いますが、上の証明は「シフト不変な事象が自明である」ことしか使っていません。そのため、独立でなくても定常でエルゴード的な列(たとえば既約で非周期的な有限 Markov 連鎖の定常列)に同じ結論が成り立ちます。
例 5.5 (Borel の正規数定理(2 進) )
x ∈ [ 0 , 1 ) x\in[0,1) x ∈ [ 0 , 1 ) の 2 進展開 x = 0. x 1 x 2 … x=0.x_1x_2\ldots x = 0. x 1 x 2 … (1 が無限に続く展開は使わない)について、ほとんどすべての x x x で
lim n → ∞ # { k ≤ n : x k = 1 } n = 1 2 \lim_{n\to\infty}\frac{\#\{k\le n:x_k=1\}}n=\frac12 n → ∞ lim n # { k ≤ n : x k = 1 } = 2 1 です。実際、例 2.4 の 2 倍写像について D k − 1 x D^{k-1}x D k − 1 x の第 1 桁が x k x_k x k で、x k = 1 [ 1 2 , 1 ) ( D k − 1 x ) x_k=\mathbf{1}_{[\frac12,1)}(D^{k-1}x) x k = 1 [ 2 1 , 1 ) ( D k − 1 x ) です。命題 3.8 から D D D はエルゴード的なので、定理 5.2 の 4 から時間平均は m ( [ 1 2 , 1 ) ) = 1 2 m([\frac12,1))=\frac12 m ([ 2 1 , 1 )) = 2 1 に収束します。同じ議論を 1 [ j / 2 ℓ , ( j + 1 ) / 2 ℓ ) \mathbf{1}_{[j/2^\ell,(j+1)/2^\ell)} 1 [ j / 2 ℓ , ( j + 1 ) / 2 ℓ ) に使えば、長さ ℓ \ell ℓ のどの 0-1 列も頻度 2 − ℓ 2^{-\ell} 2 − ℓ で現れることが分かります。
演習 6.1 易
2 倍写像 D D D と f = 1 [ 1 2 , 1 ) f=\mathbf{1}_{[\frac12,1)} f = 1 [ 2 1 , 1 ) について、x = 1 3 x=\frac13 x = 3 1 と x = 1 7 x=\frac17 x = 7 1 での時間平均 lim n A n f ( x ) \lim_nA_nf(x) lim n A n f ( x ) を求めてください。結果が 例 5.5 と矛盾しない理由を述べてください。
解答 1 3 = 0.010101 … ( 2 ) \frac13=0.010101\ldots_{(2)} 3 1 = 0.010101 … ( 2 ) で、D ( 1 3 ) = 2 3 D(\frac13)=\frac23 D ( 3 1 ) = 3 2 、D ( 2 3 ) = 1 3 D(\frac23)=\frac13 D ( 3 2 ) = 3 1 と周期 2 で回ります。f ( 1 3 ) = 0 f(\frac13)=0 f ( 3 1 ) = 0 、f ( 2 3 ) = 1 f(\frac23)=1 f ( 3 2 ) = 1 なので A n f ( 1 3 ) → 1 2 A_nf(\frac13)\to\frac12 A n f ( 3 1 ) → 2 1 です。1 7 = 0.001001 … ( 2 ) \frac17=0.001001\ldots_{(2)} 7 1 = 0.001001 … ( 2 ) で、軌道は 1 7 → 2 7 → 4 7 → 1 7 \frac17\to\frac27\to\frac47\to\frac17 7 1 → 7 2 → 7 4 → 7 1 、f f f の値は 0 , 0 , 1 0,0,1 0 , 0 , 1 なので A n f ( 1 7 ) → 1 3 A_nf(\frac17)\to\frac13 A n f ( 7 1 ) → 3 1 です。有理数(さらに一般に周期点や最終的に周期的な点)は可算個で Lebesgue 測度 0 0 0 なので、「ほとんどすべての x x x で 1 2 \frac12 2 1 」という主張とは矛盾しません。エルゴード定理の例外集合は空でないことがあり、しかも稠密になりえます。
演習 6.2 標準
測度保存系がエルゴード的であるための必要十分条件は、すべての A , B ∈ B A,B\in\mathcal{B} A , B ∈ B について
lim n → ∞ 1 n ∑ k = 0 n − 1 μ ( T − k A ∩ B ) = μ ( A ) μ ( B ) \lim_{n\to\infty}\frac1n\sum_{k=0}^{n-1}\mu(T^{-k}A\cap B)=\mu(A)\mu(B) n → ∞ lim n 1 k = 0 ∑ n − 1 μ ( T − k A ∩ B ) = μ ( A ) μ ( B ) が成り立つことであることを示してください。混合性(定義 3.4 )との違いも述べてください。
解答 (必要性)1 n ∑ k < n μ ( T − k A ∩ B ) = ⟨ 1 B , A n 1 A ⟩ \frac1n\sum_{k<n}\mu(T^{-k}A\cap B)=\langle\mathbf{1}_B,A_n\mathbf{1}_A\rangle n 1 ∑ k < n μ ( T − k A ∩ B ) = ⟨ 1 B , A n 1 A ⟩ です。系 4.3 からエルゴード的なら A n 1 A → μ ( A ) A_n\mathbf{1}_A\to\mu(A) A n 1 A → μ ( A ) が L 2 L^2 L 2 で成り立ち、内積の連続性から左辺は ⟨ 1 B , μ ( A ) 1 ⟩ = μ ( A ) μ ( B ) \langle\mathbf{1}_B,\mu(A)1\rangle=\mu(A)\mu(B) ⟨ 1 B , μ ( A ) 1 ⟩ = μ ( A ) μ ( B ) に収束します。
(十分性)A ∈ I A\in\mathcal{I} A ∈ I を取り B = A B=A B = A とすると T − k A = A T^{-k}A=A T − k A = A なので左辺は μ ( A ) \mu(A) μ ( A ) 、右辺は μ ( A ) 2 \mu(A)^2 μ ( A ) 2 で、μ ( A ) ∈ { 0 , 1 } \mu(A)\in\{0,1\} μ ( A ) ∈ { 0 , 1 } です。
混合性は数列 μ ( T − k A ∩ B ) \mu(T^{-k}A\cap B) μ ( T − k A ∩ B ) そのものの収束、エルゴード性はその Cesàro 平均の収束です。無理数回転(定理 3.7 )は後者だけを満たす例です。
演習 6.3 難
α ∉ Q \alpha\notin\mathbb{Q} α ∈ / Q とします。T \mathbb{T} T 上のすべての連続関数 f f f とすべての x ∈ T x\in\mathbb{T} x ∈ T について
1 n ∑ k = 0 n − 1 f ( x + k α ) → ∫ 0 1 f d m \frac1n\sum_{k=0}^{n-1}f(x+k\alpha)\to\int_0^1f\,dm n 1 k = 0 ∑ n − 1 f ( x + k α ) → ∫ 0 1 f d m が成り立つこと(Weyl の一様分布定理)を示してください。不連続な f f f では「すべての x x x 」が成り立たない例も挙げてください。
解答 f = e k f=e_k f = e k (k ≠ 0 k\ne0 k = 0 )なら 例 4.4 の計算で sup x ∣ A n e k ( x ) ∣ ≤ 2 / ( n ∣ 1 − e 2 π i k α ∣ ) → 0 = ∫ e k d m \sup_x|A_ne_k(x)|\le2/(n|1-e^{2\pi ik\alpha}|)\to0=\int e_k\,dm sup x ∣ A n e k ( x ) ∣ ≤ 2/ ( n ∣1 − e 2 π ik α ∣ ) → 0 = ∫ e k d m 、f = e 0 f=e_0 f = e 0 なら両辺 1 1 1 です。線形性から三角多項式 p p p について A n p → ∫ p d m A_np\to\int p\,dm A n p → ∫ p d m が x x x について一様に成り立ちます。連続な f f f と ε > 0 \varepsilon>0 ε > 0 について、Fejér の定理(定理 4.2[フーリエ級数] の 1 を θ = 2 π x \theta=2\pi x θ = 2 π x で移したもの)から sup ∣ f − p ∣ < ε \sup|f-p|<\varepsilon sup ∣ f − p ∣ < ε となる三角多項式 p p p があります。∣ A n f − A n p ∣ ≤ ε |A_nf-A_np|\le\varepsilon ∣ A n f − A n p ∣ ≤ ε 、∣ ∫ f − ∫ p ∣ ≤ ε |\int f-\int p|\le\varepsilon ∣ ∫ f − ∫ p ∣ ≤ ε なので lim sup n sup x ∣ A n f ( x ) − ∫ f ∣ ≤ 2 ε \limsup_n\sup_x|A_nf(x)-\int f|\le2\varepsilon lim sup n sup x ∣ A n f ( x ) − ∫ f ∣ ≤ 2 ε です。
不連続な例:f = 1 { k α m o d 1 : k ≥ 0 } f=\mathbf{1}_{\{k\alpha\bmod1:k\ge0\}} f = 1 { k α mod 1 : k ≥ 0 } (0 0 0 の軌道の指示関数)は可算集合の指示関数なので ∫ f = 0 \int f=0 ∫ f = 0 ですが、x = 0 x=0 x = 0 ではすべての項が 1 1 1 で時間平均は 1 1 1 です。
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.
執筆: Claude Code(Opus 5.5)、Claude Code(Sonnet 5.5)・監査 3 回(Claude Code)・最終 2026-10-12