• 検索結果がありません。

D Law of iterated logarithm .

ドキュメント内 数理物理学 - Random walk と self-avoiding walk (ページ 106-112)

Mean square displacementの指数E[Wnp]∼n がpathの「ぎざぎざ」の度合いを示す最初の指標だと すると,その次の指標として到達範囲の評価 定理100のアナロジーが成り立つ.(なお,その具体形はν (と 空間のフラクタル次元等)で決まるので,mean square displacementより詳細な性質なのに,本質は(適当 な仮定の下で) mean square displacementと同時に定まりそうに見える.厳密な結果の詳細は略す([36]の 文献を参照)が,§A.11にそのことを支持する直感的議論を掲げた.)

定理100タイプの評価の,pathの形状の性質における一つの帰結としてlaw of iterated logarithmを

(一般化を意識しながら)まず1次元simple random walkについて説明する.Law of iterated logarithm は,ブラウン運動と SRWに共通な pathの性質の一つで古くから知られている[10,§VIII.5], [8,§7.9]65

LILは[10]に言及されているように,広くマルコフ過程全体に対して成り立つことが古くから知られて いるが,当時の議論はlower boundに関してはマルコフ性が本質だった.

最後にくりこみ群の視点からの別証明の可能性に言及しておく.

D.1 Upper bound from decay estimate.

LIL上からの評価はhitting timeの上からの評価があればどんな確率連鎖でも成り立つ.

定理 95 Wn,n∈Z+,を実確率変数列(1次元確率連鎖)とする.(SRWである必要はない.)もし,C+,i>0, i= 1,2,3,u+> C+,2(1ν) と0< ν <1 が存在して,

P[ max

0kn|Wk|x]C+,1eC+,2(xnν)1/(1ν), 0< xu+nν(log logn)1ν, nC+,3, (192) が成り立てば,

nlim→∞

|Wn|

ψ+(n) C+,2(1ν) (193)

が成り立つ.ここで ψ+(n) =nν(log logn)1ν,n∈N,とおいた.

36 (i) 結論のψ+(n)の nν の係数に log logという関数があることからiterated logarithm の名前が ついている.

証明からわかるように,仮定の評価 (192)の指数部の形からこの関数形が決まってしまうので,(192) の指数部の形を表現した定理と言ってよいだろう.

(ii) Mean square displacementの観点からはxnν の程度である.仮定のxの範囲はこれからごくわ ずか遠いところまで届いたpath についても一様な上からの評価があれば十分であることを言ってい る.うまく行くとき(例えば SRW)x=O(nν+,)で仮定と同様の評価を得ることがやさしい(こ のとき,もちろん仮定の uは自由にとれる).

通常の中心極限定理はy=xnν を固定してn→ ∞としたときに仮定のような評価が(等号で)成 り立つことを主張する.必要な仮定は中心極限定理の結果と似ているが,変数の取りうる範囲がわず かに広い.

SRW(2項分布)の中心極限定理は最初から適用範囲の広い状態で知られていた[10,§VII.2,5]

(iii) (192)の評価はx=o(nν(log logn)1ν) ではいらないが,C+,1 を十分大きく取れば x=O(nν)で は必ず成り立ってしまうので,排除できる範囲はわずかだから省略した.

定理95 の証明.

定理95の証明にはBorel–Cantelli第1定理 定理56を使う.

λ >1 を任意にとる.仮定からλ > aν > λ 1

uC+,21ν なるa >1,C+,2(1ν) < u < u+ が存在する.ψ+ が 増加関数であることに注意して,仮定をn= [aj+1],x=λψ+(aj)で使うことにより,十分大きいj0 に対

65但し,SRWでは歩数無限の方向の漸近形,ブラウン運動では(self-similarityがあるのでどっち向きとも言えるが,H¨older 続性との関連でpathの性質として特に)短時間側の性質として定式化されるのが普通である点で,みかけが違う.

して

jj0

P[ max

ajk<aj+1

|Wk|

ψ+(k) C+,2(1ν)λ]

jj0

P[ max

k<aj+1|Wk|λC+,2(1ν)ψ+(aj) ]

jj0

C+,1(jloga)(λa−ν)1/(1−ν)

<∞.

よってBorel–Cantelli第1定理 定理56から,殆ど全てのω∈Ωに対して,あるj1(ω)以降のjについて max

ajk<aj+1

|Wk(ω)|

ψ+(k) < λC+,2(1ν). よって lim

k→∞

|Wk(ω)|

ψ+(k) < λC+,2(1ν).これが任意のλ >1に対して成り立つから(λ= 1 +p1,p∈N, で共

通部分をとるいつもの論法で)主張を得る.

D.2 Upper bound (Simple random walk).

Wn が 1次元SRWのとき 定理95の仮定が成り立つことを示す.

まず,命題3 の例題を復習しておく.

96 (命題3の系) {Wn}を1次元simple random walkとし,m∈Zのhitting timeTmと書くとき,

P[ max

0kn

Wkm] = P[Tmn] = 2P[Wn+1> m], m, n∈N.

証明. Hitting time の定義と,命題3 (で x方向にm ずらして時間を反転させた主張)より直ちに結論を

得る.

命題 97 {Wn} が1次元SRWのとき,任意の( >0と C >0 に対して,C >0が存在して

P[Wnx]ex2/(2n), 0< xCn1,, nC. (194) 注 37 (i) 対称性から, P[|Wn|> x√

n] についても定数を除いて同様の評価が成り立つ.従って特に,

1次元 SRWに対して 定理95の仮定が成り立つ.

(ii) もっと精密な結果が容易に得られる (定理100,101)が,定理95 の仮定の証明にはこれで十分で ある.その上証明がやたら簡単である(定理100ではStirling公式とテーラー定理の剰余項の一様評 価が不可避と思う)

なお,定理95で必要なx,nの関係に比べてやたらに広い範囲で成り立っているが,xのべきを落と した非精密な内容(定理100 を参照)なのと,下からの評価が狭い範囲でしか成り立たない.

証明. チェビシェフの不等式から,

P[Wnx]eλxE[eλX1 ]n=eλx(coshλ)n が任意のx∈Rとλ0 とn∈Z+ に対して成り立つので,λ=x/nとすると,

P[Wn x]ex2/n(coshx

n)n, x >0, n >0. C>0を cosh C

C, 3となるようにとると,0< xCn1,, nC のとき,

coshx

n 1 + x2 2n2 + x4

24n4coshx

n 1 + x2 2n2 + x4

8n4 ex2/(2n2)

となるのでP[Wn x]ex2/(2n) を得る. 定理95,系96,命題97から1次元SRWについてlaw of iterated logarithmの上からの評価を得る.

98 1次元simple random walk では(193)が成り立つ.

D.3 Lower bound (Simple random walk).

Random walk (独立確率変数列の和)の LILのlower boundは伝統的に(192)に類似の形の下からの評価 と,random walkの(時間方向の)Markov性を用いてきた.(前者については,simple random walkではド モアブル−ラプラスの中心極限定理によって可能.)

本当は細かい構造の(統計的な)繰り返し(スケール不変性)が本質であるから,この証明は不満.Zを 含むある種のフラクタル上の RWや self-repelling walk では「細かさの方向」の branching process (ス ケール方向のMarkov性)によるLILの証明が可能で,そのほうが良い証明だと思うが,技術的に難しいの で伝統的な証明をまず紹介する.

定理 99 Wn, n∈Z+,を時間的に一様な1次元random walk (実確率変数列でi.i.d.の和になっているも の)とする.もし,C±,i>0,i= 1,2,3,K >0 とu±> C±,2(1ν) と0< ν <1 が存在して,(192)

C,1eC−,2(xn−ν)1/(1−ν) P[|Wn|x], Knν< xu nν(log logn)1ν, nC3,, (195) が成り立てば,

nlim→∞

|Wn|

ψ(n) C2,1 (196)

が成り立つ.ここで ψ(n) =nν(log logn)1ν,n∈N,とおいた.

38 (i) C2, =C2,+ とできればψ+=ψ =ψ とおけるので,古典的な low of iterated logarithm

nlim→∞

|Wn|

ψ(n)= 1 (197)

を得る(1次元SRWではC2, =C2,+ でこのことが成り立っている定理100)

(ii)

(iii) (192)のときと違って,(195) の適用範囲の xの下限を明示的に限った.これは x=o(nν) のとき

PC1 を主張してしまうので,不要に強い仮定が入るのを嫌ったため(LILが成り立つ状況でこれが 成り立たないとは考えにくいが).また,x=o(nν(log logn)1ν)ではいらないので下限はそれだけ 上がるが,式を面倒にするのを避けた.

(iv) Upper bound, lower bound ともに連続極限(例えばブラウン運動)でも同様の性質が知られているが,

このときは x=O(tν) で評価があれば十分なようである.この違いがどこから来るのか詳細を確認 していないが,連続極限では自己相似性が正確に成り立つ(確率連鎖では,1歩より細かい構造はな い)ことから来るのだろうと思われる.

今日ではまずブラウン運動について LILを証明し,そのあとで SRWはブラウン運動に近い(almost sure invariance principle)ことを言ってSRWLILを証明するのが常道のようである[8,§7.9].こ れは上記の視点であろう.極めて重要な視点だが,連続極限そのものへの言及は確率過程論の講義に 期待して,この講義では意識的に避ける.

定理99 の証明.

定理99の証明にはBorel–Cantelli 第2定理 定理57を用いる.

0< λ <1を任意に取る.仮定からC,2(1ν)< u < u がとれ,また,a0>1がとれて,a > a0 なら ばλ

a a−1

ν

<1 とできる.

Dj = W[aj] −W[aj−1] とおくと時間的に一様なRWという仮定から{Dj} は独立で Dj の分布は W[aj][aj−1] に等しい.ψ が増加関数であることに注意して,仮定をn= [aj][aj1],x=λψ(aj)C2,1 で使うことにより,十分大きいj0に対して

jj0

P[ |Dj|

ψ(aj)λC2,1 ]

jj0

C,1(jloga)(λ(a+0a−1)ν)1/(1−ν)

=∞.

a+ 0と書いたのはaより大きい数で,仮定の範囲でうまく取る.)

よってBorel–Cantelli第2定理 定理57から(定理95の証明の最後の部分と同様の論理で),

lim

j→∞

|W[aj]−W[aj−1]|

ψ(aj) C2,1, a.s..

他方,仮定の下で 定理95 が成り立つから,殆ど全てのω と任意の( >0 に対してj1 がとれて,j > j0

ならば

W[aj]

ψ(aj)C2,1 1−(−(1 +()ψ(aj1) ψ+(aj)

C2,

C2,+

1ν

1−(−(1 +() 1 aν

C2,

C2,+

1ν

. よって,

nlim→∞

|Wn|

ψ(n)C2,1 lim

j→∞

|W[aj]|

ψ(aj)C2,1 1 1 aν

C2,

C2,+

1ν

. これが任意のa > a0 で成り立つから,主張を得る.

次の性質は中心極限定理の結果と似ているが,変数の取りうる範囲が 命題1より強い([10,§VII.2]).今 まで通りAn/Bn1 となることをAn∼Bn 等と書く.

定理 100 (ドモアブル−ラプラスの極限定理) n∈Z+,y∈Rに対して hn = 2/√

n, xy = (y−n2)hとお く.α, β∈Rが hnx3α0,hx3β 0,を満たしながらn→ ∞とするとき,

2n

αkβ

n k

xβ−1/2 xα−1/2

ey2/2 dy

2π.

22 以下に[10,§VII.2] に書かれている定理100の証明の骨子を引用するが,いい書き方とは思えない.

例えば,原典ではx3k/√

n→0のとき n

k

についての漸近形が成り立つことから,両者の比の1 からの ずれがA/n+B|xk+1|3/√

nで押さえられる,と書いてあるが,これは間違いである.例えば

|xk+1|3/√ n で押さえられていれば十分.はっきり言うと,原典の書き方では o(1)ということしか分からない.

n k

kについての和の漸近的評価が必要なので,o(1)評価では不十分ではないかという疑念が生 じる.よって,以下の証明は書き方が不十分である.これを補足して証明を完成せよ.(Statement の書き

方も(195) に合った記述の方が望ましい.)

定理100の証明. テーラーの定理から,各 k 毎に

xk+1/2 xk−1/2

ey2/2dy =heξ2/2 を満たす xk −hn 2 < ξ <

xk+hn

2 が存在するので,

hnex2k/2=e(ξ2x2k)/2

xk+1/2 xk−1/2

ey2/2dy.

nを十分大きく取ればk に無関係に2−x2k|を小さくとれるので,

(∀( >0)∃n0; (∀n > n0)e,

xβ+1/2 xα−1/2

ey2/2dyhn

αkβ

ex2k/2e,

xβ+1/2 xα−1/2

ey2/2dy が任意の(nとともに変化してもかまわない)α,β に対して成り立つ.ここまでは問題ないだろう.

一方,Stirlingの公式(§A.1)と(1 + 1n)n∼eから,

2n n

k

∼hnex2k/2, hnx3k0.

これより(ここが私は心配),

2n

αkβ

n k

∼hn

αkβ

ex2k/2, hnx3α0, hnx3β0.

[10,§VII.2]では以上をもって良しとしているように見える.

1次元SRWではP[Wn =k] =

n (n+k)/2

だから,誤差関数の漸近形 命題55を 定理100に適 用すると次を得る.

101 P[|Wn|x] n

2π 1

xex2/(2n),x3/n20,n→ ∞. 今までの諸結果とその注を合わせると次の最終結果を得る.

102 (1次元 SRWLIL) 1次元simple random walkでは (197)が成り立つ.

D.4 Hitting timeshort time estimate から lower bound への十分条件.

§D.5でLILのlower bound の別証明を与えるための準備を行う.

命題 103 Zに値を取る確率変数列Wk,k∈Z+,に対して,0 を除く G˜n = 2nZの hitting timeTn と 書く:

Tn= inf{k∈N|Wk∈G˜n\ {0}}

このとき,λ >2(即ち0< ν = log 2

logλ <1)に対して

(logn)(1ν)λnTnc , i.o., (198) となる66 c >0とがあれば,

lim

k→∞

|Wk| ψ(k) cν が成り立つ.ここで ψ(k) =kν(log logk)1ν

39 SRWでなくてよい!また,LILの下限は遠くに行くという主張だから,1歩毎に隣に移るwalk

ある必要もない.

証明. (198)の不等式が成り立つ w∈Ωとn∈Z+ をとって k=Tn(w)とおくと,

|Wk(w)|2n, kcλn(logn)(1ν). よって

2ncνkν(logn)1ν. (199)

先ずn3で

nlogk−logc

logλ + (1−ν)log logn

logλ logk−logc logλ を得て,これからさらに,

lognlog logk+ log(1 logc

logk)log logλ となるので,任意の ( >0 に対して,k0 がとれて

logn(1−() log logk, kk0.

これを(199)に代入すると

|Wk(w)|2ncνkν(1−()1ν(log logk)1ν を得る.仮定から,これが成り立つkが無数にあるから

nlim→∞

|Wk(w)|

ψ(k) (1−()1νcν が成り立つ.( >0 は任意だったから,

nlim→∞

|Wk(w)| ψ(k) cν だが,このようなwが確率1 で存在するのだから主張を得る.

66i.o.infinitely often.与式の成り立つnが無限個あること.もちろん,そうなる確率が1ということ.

ドキュメント内 数理物理学 - Random walk と self-avoiding walk (ページ 106-112)