ランダムウォークの原点への復帰(return of a random walk to the origin)とは、原点から出発し、確率 $\frac12$ ずつで独立に $+1$ か $-1$ 動く数直線上の点(1次元のランダムウォーク)が、再び原点に戻る現象である。$2n$ 回目に原点にいる確率は $u_{2n}=\binom{2n}{n}\frac{1}{4^n}$ で、$2n$ 回目まで一度も原点に戻らない確率も $u_{2n}$ に等しい。$\frac{1}{2\sqrt n}\le u_{2n}\le\frac{1}{\sqrt{2n+1}}$ から、$2n$ 回目までに原点に戻る確率 $1-u_{2n}$ は $n\to\infty$ で $1$ に近づく。$2n$ 回目までに戻る回数の期待値は $(2n+1)u_{2n}-1$ で、限りなく大きくなる。
前提知識: 反復試行の確率, 最短経路の数え上げと鏡像原理, 確率変数の期待値と分散
数直線の原点に点を置き、硬貨を投げる。表が出たら右へ $1$、裏が出たら左へ $1$ だけ点を動かす。これを何回もくり返すと、点は右へ行ったり左へ行ったりしながら、ときどき原点に戻ってくる。図 1 は、$20$ 回ずつ投げた 3 回の試しの記録で、横に回数、縦に点の位置をとり、原点に戻った回に丸を付けた。
硬貨を 20 回投げて、表なら +1、裏なら −1 動いた点の位置の記録を 3 回分並べた図。丸は原点に戻った回で、1 本目は 5 回、2 本目は 3 回、3 本目は 5 回戻っている
2 本目の記録では、$13$ 回目まで一度も原点に戻っていない。表と裏は同じ確率で出るのに、点が片側に長くとどまることがある。まず、硬貨を 4 回投げる場合を全部数えてみる。
表を $+$、裏を $-$ と書く。4 回の結果の並びは $2^4=16$ 通りで、どれも確率 $\dfrac1{16}$ で起こる。
(1) 4 回目に原点にいるのは、$+$ と $-$ がちょうど 2 個ずつのときで、並び方は $\dbinom42=6$ 通りある:
$$
{+}{+}{-}{-},\quad {+}{-}{+}{-},\quad {+}{-}{-}{+},\quad {-}{+}{+}{-},\quad {-}{+}{-}{+},\quad {-}{-}{+}{+}
$$
よって、4 回目に原点にいる確率は $\dfrac6{16}=\dfrac38$ である。
(2) 1 回目から 4 回目まで、一度も原点に戻らない並びを数える。1 回目が $+$ のとき、位置は $1$ で、2 回目に $-$ が出ると原点に戻るので 2 回目は $+$ で位置 $2$ である。3 回目は $+$(位置 $3$)でも $-$(位置 $1$)でもよい。位置 $3$ からは 4 回目に何が出ても原点に戻らず、位置 $1$ からは 4 回目に $+$ でなければならない。よって
$$
{+}{+}{+}{+},\quad {+}{+}{+}{-},\quad {+}{+}{-}{+}
$$
の 3 通りである。1 回目が $-$ のときも、符号を入れ替えた 3 通りがある。合わせて 6 通りで、確率は $\dfrac6{16}=\dfrac38$ である。
(1) と (2) の確率は、どちらも $\dfrac38$ で等しい(図 2)。
上の段は 4 回目に原点にいる 6 本、下の段は 4 回目まで一度も原点に戻らない 6 本の記録。どちらも 6 本ずつある
「4 回目に原点にいる」と「4 回目まで一度も原点に戻らない」は、まったく違う出来事である。それなのに確率が同じになった。この記事では、これが偶然ではなく、何回投げても成り立つことを示し、それを使って次の問いに答える。
| 高校の計算 | この記事の言葉 | 大学の言葉 |
|---|---|---|
| 表と裏の並べ方を数える | 長さ $n$ の歩み | 確率変数の列 $X_1,X_2,\ldots$ |
| 碁盤の目の最短経路 | 歩みを格子の道とみる | ランダムウォークの経路 |
| 鏡像原理・投票定理 | 原点に戻らない歩みの数 | 鏡像原理 |
| 中央の二項係数 $\binom{2n}n$ の大きさ | $u_{2n}$ の大きさ | Stirling の公式による漸近 |
次の仮定のもとで、数直線の上を動く点を考える。
(R1) $k=1,2,3,\ldots$ について、$k$ 回目の動き $X_k$ は $+1$ か $-1$ で、それぞれ確率 $\dfrac12$ をとる。$X_1,X_2,\ldots$ は互いに独立である(硬貨を独立に投げる)。
(R2) 点ははじめ原点にあり、$n$ 回動いた後の位置を $S_n=X_1+X_2+\cdots+X_n$ とする($S_0=0$)。
この点の動きを 1 次元のランダムウォーク という。$n$ 回目までの動きの並び $(X_1,\ldots,X_n)$ を 長さ $n$ の歩み という。長さ $n$ の歩みは $2^n$ 通りあり、条件 (i) から、どれも確率 $\left(\dfrac12\right)^n$ で起こる。$k\ge1$ について $S_k=0$ となるとき、点は $k$ 回目に原点に戻る という(出発の $S_0=0$ は戻ったと数えない)。
確率がすべて等しいので、$n$ 回目までの出来事の確率は「その出来事が起こる長さ $n$ の歩みの数」を $2^n$ で割れば求まる。この記事の計算はすべてこの形で、反復試行の確率 と同じく、並び方を数える問題になる。
長さ $6$ の歩み $(+1,-1,-1,+1,+1,+1)$ では、位置は
$$
S_1=1,\quad S_2=0,\quad S_3=-1,\quad S_4=0,\quad S_5=1,\quad S_6=2
$$
と動き、2 回目と 4 回目に原点に戻る。長さ $6$ の歩み $(+1,+1,-1,+1,-1,-1)$ では、位置は $1,2,1,2,1,0$ と動き、6 回目に初めて原点に戻る。
原点に戻れるのは偶数回目だけである。
$n$ 回動いた後、$+1$ の動きが $a$ 回、$-1$ の動きが $b$ 回あったとすると、$S_n=a-b$、$n=a+b$ である。特に、$S_n=0$ となるのは $n$ が偶数で $a=b=\dfrac n2$ のときに限る。したがって、奇数回目に原点にいる確率は $0$ である。
$S_n$ は $+1$ を $a$ 個、$-1$ を $b$ 個足したものなので $S_n=a-b$ であり、動きの個数を数えて $n=a+b$ である。$S_n=0$ なら $a=b$ で、$n=a+b=2a$ は偶数、$a=b=\dfrac n2$ である。$n$ が奇数なら $a=b$ となることはないので、$S_n=0$ となる歩みは 1 つもなく、確率は $0$ である。$\square$
そこで、偶数回目に原点にいる確率に名前を付ける。
$n=1,2,3,\ldots$ について、$2n$ 回目に原点にいる確率を
$$
u_{2n}=P(S_{2n}=0)
$$
と書く。また $u_0=1$ と決める(出発の時点では必ず原点にいる)。
歩みは、最短経路の数え上げと鏡像原理 の格子の道で表せる。$+1$ の動きを右へ 1(R)、$-1$ の動きを上へ 1(U)に対応させると、長さ $n$ の歩みは、原点から出発する R と U の $n$ 歩の道になる。$k$ 歩目の後の点を $(a_k,b_k)$ とすると、$a_k$ は $+1$ の回数、$b_k$ は $-1$ の回数なので
$$
S_k=a_k-b_k
$$
である。したがって
ex-rwr-start では、4 回目まで一度も原点に戻らない確率と、4 回目に原点にいる確率が等しかった。6 回でも数えてみる。
長さ $6$ の歩みで、$S_1,\ldots,S_6$ がすべて正のものを、最後の位置 $S_6$ で分けて数える。$S_6$ は偶数で(prop-rwr-parity と同じ理由で $S_6=a-b$、$a+b=6$)、正なので $2,4,6$ のどれかである。
(1) $S_6=6$:$+1$ が 6 回の 1 通り。
(2) $S_6=4$:$+1$ が 5 回、$-1$ が 1 回。$-1$ を 1 回目に置くと $S_1=-1$、2 回目に置くと $S_2=0$ になるので、$-1$ は 3〜6 回目のどこかで、4 通り。この 4 通りでは、$-1$ の直前の位置が $2$ 以上なので、位置は $1$ 以上のままである。
(3) $S_6=2$:$+1$ が 4 回、$-1$ が 2 回の並び $\dbinom62=15$ 通りのうち、位置がずっと正のものは
$$
{+}{+}{+}{+}{-}{-},\ \ {+}{+}{+}{-}{+}{-},\ \ {+}{+}{+}{-}{-}{+},\ \ {+}{+}{-}{+}{+}{-},\ \ {+}{+}{-}{+}{-}{+}
$$
の 5 通りである。
合わせて $1+4+5=10$ 通りある。すべて負のものも、符号を入れ替えて $10$ 通りあるので、6 回目まで一度も原点に戻らない歩みは $20$ 通りで、$\dbinom63=20$ と一致する。確率は $\dfrac{20}{64}=\dfrac5{16}$ で、$u_6=\dbinom63\dfrac1{4^3}=\dfrac{20}{64}=\dfrac5{16}$ に等しい。
ex-rwr-n3 の (3) の $5$ 通りは、最短経路の数え上げと鏡像原理 の投票定理の場面そのものである。$+1$ を候補者 A の票、$-1$ を候補者 B の票とみると、「位置がずっと正」は「開票の間、A の票数がつねに B より多い」ことになる。投票定理によれば、A が $a$ 票、B が $b$ 票($a>b$)のとき、その確率は $\dfrac{a-b}{a+b}$ である。$a=4$、$b=2$ なら $\dfrac26\cdot\dbinom62=\dfrac13\cdot15=5$ 通りで、数え上げと一致する。一般の場合の証明は、この見方で最後の位置ごとに数えて足す。
$n$ を正の整数とする。
(1) $2n$ 回目に原点にいる確率は
$$
u_{2n}=\binom{2n}n\frac1{4^n}
$$
である。
(2) 長さ $2n$ の歩みのうち、$S_1,S_2,\ldots,S_{2n}$ がどれも $0$ でないもの($2n$ 回目まで一度も原点に戻らないもの)の数は $\dbinom{2n}n$ である。したがって
$$
P(S_1\ne0,\ S_2\ne0,\ \ldots,\ S_{2n}\ne0)=u_{2n}
$$
である。
方針:(1) は並べ方を数えるだけである。(2) は、まず「原点に戻らない歩みは、ずっと正か、ずっと負か」に分ける。ずっと正の歩みを最後の位置 $S_{2n}=2j$ ごとに投票定理で数え、その和が望遠鏡のように消えて $\dfrac12\dbinom{2n}n$ になることを示す。
段 1((1) の証明)。prop-rwr-parity により、$S_{2n}=0$ となるのは、$2n$ 回の動きのうち $+1$ がちょうど $n$ 回、$-1$ が $n$ 回のときである。$+1$ を置く $n$ か所の選び方は $\dbinom{2n}n$ 通りある(格子の道で言えば、$(0,0)$ から $(n,n)$ への最短経路の数である)。長さ $2n$ の歩みは $2^{2n}=4^n$ 通りで、どれも同じ確率なので、$u_{2n}=\dbinom{2n}n\Big/4^n$ である。
段 2(ずっと正か、ずっと負か)。$S_1,\ldots,S_{2n}$ がどれも $0$ でない歩みを考える。$S_1$ は $+1$ か $-1$ である。$S_1=1$ なのに、ある $k$ で $S_k<0$ となったとする。位置は整数で、1 回に $1$ ずつしか変わらないので、$1$ 以上の値から $-1$ 以下の値に移る途中で、必ず $0$ を通る。これは「どれも $0$ でない」に反する。よって $S_1=1$ なら $S_1,\ldots,S_{2n}$ はすべて正である。同じように、$S_1=-1$ ならすべて負である。
また、歩みのすべての動きの符号を入れ替える($X_k$ を $-X_k$ にする)と、位置もすべて符号が入れ替わる。この入れ替えは、ずっと正の歩みとずっと負の歩みを 1 対 1 に対応させる(2 回入れ替えると元に戻る)。したがって
$$
(\text{原点に戻らない歩みの数})=2\times(\text{ずっと正の歩みの数})
$$
である。
段 3(ずっと正の歩みを、最後の位置ごとに数える)。ずっと正の歩みでは $S_{2n}>0$ で、prop-rwr-parity と同じく $S_{2n}=a-b$、$a+b=2n$ なので $S_{2n}$ は偶数である。$S_{2n}=2j$($j=1,2,\ldots,n$)とおくと、$+1$ の回数は $a=n+j$、$-1$ の回数は $b=n-j$ である。$+1$ を A の票、$-1$ を B の票とみると、「$S_1,\ldots,S_{2n}$ がすべて正」は「最初の 1 票から最後まで、A の票数がつねに B の票数より多い」ことと同じである。$a>b\ge0$ なので、最短経路の数え上げと鏡像原理 の定理「投票定理」が使えて、A と B の並べ方 $\dbinom{2n}{n+j}$ 通りのうち、その割合は $\dfrac{a-b}{a+b}=\dfrac{2j}{2n}=\dfrac jn$ である。よって、$S_{2n}=2j$ となるずっと正の歩みの数は
$$
\frac jn\binom{2n}{n+j}
$$
である。
段 4(差の形に書きかえる)。$m=n+j$ とおく($n+1\le m\le2n$)。階乗で書くと
$$
\binom{2n-1}{m-1}=\frac{(2n-1)!}{(m-1)!\,(2n-m)!}=\frac m{2n}\cdot\frac{(2n)!}{m!\,(2n-m)!}=\frac m{2n}\binom{2n}m
$$
である(分母と分子に $2n$ と $m$ を掛けて、$(2n)!=2n\cdot(2n-1)!$、$m!=m\cdot(m-1)!$ を使った)。同じように、$m\le2n-1$ なら
$$
\binom{2n-1}{m}=\frac{(2n-1)!}{m!\,(2n-1-m)!}=\frac{2n-m}{2n}\cdot\frac{(2n)!}{m!\,(2n-m)!}=\frac{2n-m}{2n}\binom{2n}m
$$
であり、$m=2n$ のときは両辺とも $0$ である($\dbinom{2n-1}{2n}=0$ と約束する)。2 つの式を引くと
$$
\binom{2n-1}{m-1}-\binom{2n-1}{m}=\frac{m-(2n-m)}{2n}\binom{2n}m=\frac{2j}{2n}\binom{2n}{n+j}=\frac jn\binom{2n}{n+j}
$$
となる($m-(2n-m)=2m-2n=2j$)。つまり、段 3 の数は
$$
\frac jn\binom{2n}{n+j}=\binom{2n-1}{n+j-1}-\binom{2n-1}{n+j}
$$
と差の形に書ける。
段 5(足すと消える)。段 4 の式を $j=1,2,\ldots,n$ について足すと、
$$
\sum_{j=1}^n\left\{\binom{2n-1}{n+j-1}-\binom{2n-1}{n+j}\right\}=\left\{\binom{2n-1}{n}-\binom{2n-1}{n+1}\right\}+\left\{\binom{2n-1}{n+1}-\binom{2n-1}{n+2}\right\}+\cdots+\left\{\binom{2n-1}{2n-1}-\binom{2n-1}{2n}\right\}
$$
となり、隣り合う中かっこの間で項が打ち消し合って、最初の項と最後の項だけが残る:
$$
\binom{2n-1}{n}-\binom{2n-1}{2n}=\binom{2n-1}{n}
$$
段 4 の 2 つ目の式で $m=n$ とおくと $\dbinom{2n-1}n=\dfrac{2n-n}{2n}\dbinom{2n}n=\dfrac12\dbinom{2n}n$ なので、ずっと正の歩みの数は $\dfrac12\dbinom{2n}n$ である。
段 6(結論)。段 2 より、原点に戻らない歩みの数は $2\times\dfrac12\dbinom{2n}n=\dbinom{2n}n$ である。これを $4^n$ で割ると、確率は $\dbinom{2n}n\Big/4^n=u_{2n}$ で、段 1 の値に等しい。$\square$
$n=3$ では、段 3 の数は $j=1,2,3$ について
$$
\frac13\binom64=\frac{15}3=5,\qquad\frac23\binom65=\frac{2\cdot6}3=4,\qquad\frac33\binom66=1
$$
で、ex-rwr-n3 の (3)・(2)・(1) の数と一致する。段 4 の差の形では
$$
5=\binom53-\binom54=10-5,\qquad4=\binom54-\binom55=5-1,\qquad1=\binom55-\binom56=1-0
$$
で、足すと $10-5+5-1+1-0=10=\binom53=\dfrac12\binom63$ となる。
$u_{2n}$ の値は、次の漸化式で順に計算できる。
$n\ge1$ について
$$
u_{2n}=\frac{2n-1}{2n}\,u_{2n-2}
$$
である。したがって
$$
u_{2n}=\frac12\cdot\frac34\cdot\frac56\cdots\frac{2n-1}{2n}
$$
である。
thm-rwr-main の (1) から
$$
\frac{u_{2n}}{u_{2n-2}}=\frac{\binom{2n}n}{\binom{2n-2}{n-1}}\cdot\frac{4^{n-1}}{4^n}=\frac{(2n)!\,\{(n-1)!\}^2}{(n!)^2\,(2n-2)!}\cdot\frac14=\frac{2n(2n-1)}{n^2}\cdot\frac14=\frac{2n-1}{2n}
$$
である($(2n)!=2n(2n-1)\cdot(2n-2)!$、$n!=n\cdot(n-1)!$ で約分した)。$n=1$ でも $u_2=\dfrac12=\dfrac12u_0$ で成り立つ。この式を $n,n-1,\ldots,1$ と順に使い、$u_0=1$ とすると積の形を得る。$\square$
| $n$ | $1$ | $2$ | $3$ | $4$ | $5$ |
|---|---|---|---|---|---|
| $\binom{2n}n$ | $2$ | $6$ | $20$ | $70$ | $252$ |
| $4^n$ | $4$ | $16$ | $64$ | $256$ | $1024$ |
| $u_{2n}$ | $\frac12$ | $\frac38$ | $\frac5{16}$ | $\frac{35}{128}$ | $\frac{63}{256}$ |
| 小数 | $0.500$ | $0.375$ | $0.313$ | $0.273$ | $0.246$ |
漸化式でも、$u_4=\dfrac34\cdot\dfrac12=\dfrac38$、$u_6=\dfrac56\cdot\dfrac38=\dfrac5{16}$、$u_8=\dfrac78\cdot\dfrac5{16}=\dfrac{35}{128}$ と確かめられる。$u_{2n}$ は $n$ とともに減っていくが、減り方はゆっくりで、$u_{100}$($n=50$)でも約 $0.0796$ ある。
thm-rwr-main の (2) により、「$2n$ 回目までに少なくとも 1 回原点に戻る」確率は、余事象の確率を $1$ から引いて $1-u_{2n}$ である。$n=2$ なら $1-\dfrac38=\dfrac58$、$n=5$ なら $1-\dfrac{63}{256}=\dfrac{193}{256}\approx0.754$ である。$n$ を大きくするとこの確率が $1$ に近づくかどうかは、$u_{2n}$ が $0$ に近づくかどうかで決まる。それを、prop-rwr-recur の積の形を上下から評価して確かめる。
$n\ge1$ について
$$
\frac1{2\sqrt n}\le u_{2n}\le\frac1{\sqrt{2n+1}}
$$
である。
方針:積 $u_{2n}=\dfrac12\cdot\dfrac34\cdots\dfrac{2n-1}{2n}$ の 2 乗を考え、各因子の 2 乗を、掛けると約分で消えていく分数で上下からはさむ。
段 1(上からの評価)。$k\ge1$ について、$(2k-1)(2k+1)=4k^2-1<4k^2$ なので、両辺を正の数 $2k(2k+1)$ で割ると
$$
\frac{2k-1}{2k}<\frac{2k}{2k+1}
$$
である。両辺に正の数 $\dfrac{2k-1}{2k}$ を掛けると
$$
\left(\frac{2k-1}{2k}\right)^2<\frac{2k-1}{2k}\cdot\frac{2k}{2k+1}=\frac{2k-1}{2k+1}
$$
である。これを $k=1,2,\ldots,n$ について掛け合わせると、右辺の積は
$$
\frac13\cdot\frac35\cdot\frac57\cdots\frac{2n-1}{2n+1}=\frac1{2n+1}
$$
と、隣り合う分数の分子と分母が約分で消える。よって $u_{2n}^2<\dfrac1{2n+1}$ で、両辺の正の平方根をとって $u_{2n}<\dfrac1{\sqrt{2n+1}}$ である。
段 2(下からの評価)。$k\ge2$ について、$(2k-1)^2=4k^2-4k+1>4k^2-4k=2k(2k-2)$ なので、両辺を正の数 $2k(2k-1)$ で割ると
$$
\frac{2k-1}{2k}>\frac{2k-2}{2k-1}
$$
である。両辺に $\dfrac{2k-1}{2k}$ を掛けると
$$
\left(\frac{2k-1}{2k}\right)^2>\frac{2k-1}{2k}\cdot\frac{2k-2}{2k-1}=\frac{2k-2}{2k}=\frac{k-1}k
$$
である。$k=1$ の因子は $\left(\dfrac12\right)^2=\dfrac14$ のままにして、$k=2,\ldots,n$ の不等式を掛け合わせると
$$
u_{2n}^2\ge\frac14\cdot\frac12\cdot\frac23\cdots\frac{n-1}n=\frac14\cdot\frac1n=\frac1{4n}
$$
である($n=1$ では右辺の積は $\dfrac14$ だけで、等号が成り立つ)。正の平方根をとって $u_{2n}\ge\dfrac1{2\sqrt n}$ である。$\square$
$n\ge1$ について、$2n$ 回目までに少なくとも 1 回原点に戻る確率は $1-u_{2n}$ であり、
$$
\lim_{n\to\infty}(1-u_{2n})=1
$$
である。
方針:前半は thm-rwr-main の (2) の余事象である。後半は lem-rwr-bounds の上の評価とはさみうちの原理による。
段 1(前半)。「$2n$ 回目までに少なくとも 1 回原点に戻る」の余事象は「$S_1,\ldots,S_{2n}$ がどれも $0$ でない」である。thm-rwr-main の (2) より余事象の確率は $u_{2n}$ なので、求める確率は $1-u_{2n}$ である。
段 2(後半)。lem-rwr-bounds より $0< u_{2n}\le\dfrac1{\sqrt{2n+1}}$ である。$n\to\infty$ で $\dfrac1{\sqrt{2n+1}}\to0$ なので、はさみうちの原理により $u_{2n}\to0$ であり、$1-u_{2n}\to1$ である。$\square$
「いつかは原点に戻る確率」は、無限に続く歩みについての確率で、長さが有限の歩みを数えるこの記事の方法では直接は扱えない。そこで、GS06 の Example 12.1 と同じく、「$2n$ 回目までに戻る確率」の $n\to\infty$ での極限を「いつかは戻る確率」と定める。thm-rwr-eventual は、この意味で 原点に戻る確率は $1$ であることを示している。
$1-u_{2n}$ は $1$ に近づくが、その近づき方は遅い(図 3)。ex-rwr-bounds の $n=50$ でも、$100$ 回投げて一度も原点に戻らない確率が約 $8\%$ ある。
2n 回目までに原点に戻る確率 1 − u_2n の n による変化。1 に近づくが遅い。帯は上下の評価からわかる範囲
$u_{2n}$ の大きさは、もっと正確に言える。Wallisの公式(高校数学) の系「Wallis の積と中央の二項係数」は、$\dbinom{2n}n\sim\dfrac{4^n}{\sqrt{\pi n}}$、すなわち
$$
\lim_{n\to\infty}\frac{\binom{2n}n}{4^n}\sqrt{\pi n}=1
$$
を示している。これを $u_{2n}$ の言葉で書くと
$$
\lim_{n\to\infty}u_{2n}\sqrt{\pi n}=1
$$
である。これは「$u_{2n}$ と $\dfrac1{\sqrt{\pi n}}$ の 比 が $1$ に近づく」という主張である(2 つの差も $0$ に近づくが、どちらも $0$ に近づくので、差が $0$ に近づくことだけでは大きさの見積もりにならない)。比 $u_{2n}\sqrt{\pi n}$ は $n=5$ で約 $0.975$、$n=50$ で約 $0.9975$ である(有限の $n$ の値は目安で、極限の証明は上の系による)。GS06 の第 12.1 節の練習問題 4 は、これを Stirling の公式から導く問題である。lem-rwr-bounds を書きかえると $\dfrac{\sqrt\pi}2\le u_{2n}\sqrt{\pi n}\le\sqrt{\dfrac{\pi n}{2n+1}}<\sqrt{\dfrac\pi2}$ で、比を $0.886$ 以上 $1.254$ 未満にはさむ評価になっている。極限の値 $1$ はこの間にある。
原点に戻る確率が $1$ なら、1 回戻った後も同じことがくり返され、何回も戻ってきそうである。$2n$ 回目までに原点に戻る回数の期待値を求めて、これを確かめる。
$2n$ 回目までに原点に戻る回数($S_k=0$ となる $k$、$1\le k\le2n$ の個数)を $N_{2n}$ とすると
$$
E[N_{2n}]=u_2+u_4+\cdots+u_{2n}=(2n+1)u_{2n}-1
$$
である。さらに $E[N_{2n}]>\sqrt n-1$ で、$n\to\infty$ のとき $E[N_{2n}]$ は限りなく大きくなる。
方針:$N_{2n}$ を「$k$ 回目に原点にいれば $1$、いなければ $0$」という確率変数の和に分け、期待値の線形性を使う。和の公式は $n$ についての帰納法で示す。
段 1(和に分ける)。$k=1,\ldots,2n$ について、$S_k=0$ なら $1$、そうでなければ $0$ をとる確率変数を $I_k$ とする。$N_{2n}=I_1+I_2+\cdots+I_{2n}$ である。$I_k$ の期待値は $1\cdot P(S_k=0)+0\cdot P(S_k\ne0)=P(S_k=0)$ である。確率変数の期待値と分散 の命題「期待値の線形性」により
$$
E[N_{2n}]=\sum_{k=1}^{2n}P(S_k=0)
$$
である。prop-rwr-parity より奇数の $k$ の項は $0$ で、偶数の $k=2j$ の項は $u_{2j}$ なので、$E[N_{2n}]=u_2+u_4+\cdots+u_{2n}$ である。
段 2(和の公式)。$u_2+\cdots+u_{2n}=(2n+1)u_{2n}-1$ を $n$ についての帰納法で示す。$n=1$ では、左辺 $u_2=\dfrac12$、右辺 $3\cdot\dfrac12-1=\dfrac12$ で成り立つ。$n$ で成り立つとすると、$n+1$ では
$$
u_2+\cdots+u_{2n}+u_{2n+2}=(2n+1)u_{2n}-1+u_{2n+2}
$$
である。prop-rwr-recur より $u_{2n+2}=\dfrac{2n+1}{2n+2}u_{2n}$、つまり $(2n+1)u_{2n}=(2n+2)u_{2n+2}$ なので、右辺は
$$
(2n+2)u_{2n+2}-1+u_{2n+2}=(2n+3)u_{2n+2}-1=\bigl(2(n+1)+1\bigr)u_{2(n+1)}-1
$$
となり、$n+1$ でも成り立つ。
段 3(限りなく大きくなる)。lem-rwr-bounds の下の評価 $u_{2n}\ge\dfrac1{2\sqrt n}$ から
$$
E[N_{2n}]=(2n+1)u_{2n}-1\ge\frac{2n+1}{2\sqrt n}-1>\frac{2n}{2\sqrt n}-1=\sqrt n-1
$$
である。$n\to\infty$ で $\sqrt n-1$ は限りなく大きくなるので、$E[N_{2n}]$ も限りなく大きくなる。$\square$
| $n$ | $u_{2n}$ | $1-u_{2n}$ | $E[N_{2n}]$ | $\sqrt n-1$ |
|---|---|---|---|---|
| $5$ | $0.246$ | $0.754$ | $1.707$ | $1.236$ |
| $10$ | $0.176$ | $0.824$ | $2.700$ | $2.162$ |
| $50$ | $0.0796$ | $0.920$ | $7.039$ | $6.071$ |
| $100$ | $0.0564$ | $0.944$ | $10.326$ | $9.000$ |
Wallisの公式(高校数学) の系をもう一度使うと、$\dfrac{E[N_{2n}]+1}{2\sqrt{n/\pi}}=\dfrac{2n+1}{2n}\cdot u_{2n}\sqrt{\pi n}$ は $1$ に近づく($\dfrac{2n+1}{2n}\to1$ と、上の極限による)。つまり $E[N_{2n}]+1$ と $2\sqrt{n/\pi}$ の比が $1$ に近づき、戻る回数の期待値は $\sqrt n$ の程度で増える。GS06 の Example 12.3 は、同じ期待値を母関数を使って求めている。
Banachのマッチ箱問題 の定理「残りの本数の期待値」は、マッチが $N$ 本ずつ入った箱で、空に気づいたときにもう一方に残っている本数 $R$ の期待値を $E[R]=\dfrac{(2N+1)\binom{2N}N}{4^N}-1$ と求めている。これは cor-rwr-count の $(2n+1)u_{2n}-1$ で $n=N$ としたものと同じ式である。$N=n=1,\ldots,8$ について全部の歩みを数えて比べると、期待値だけでなく、「$2n$ 回目までに原点に戻る回数が $r$ である確率」と「残りの本数が $r$ である確率」$\dbinom{2n-r}n\left(\dfrac12\right)^{2n-r}$ が、すべての $r$ で一致した。一般の $n$ での分布の一致は、この記事では扱わない(証明しない)。
thm-rwr-main の (2) の等式「$2n$ 回目まで一度も原点に戻らない確率 $=$ $2n$ 回目に原点にいる確率」は、def-rwr-walk の条件と、回数が偶数であることに支えられている。
| 外す条件 | 反例 | 成り立たなくなること |
|---|---|---|
| 条件 (i) の確率 $\frac12$(左右の確率を偏らせる。独立性は保つ) | ex-rwr-cx-biased:表の確率 $p=\frac23$、$n=1$ | 戻らない確率 $\frac59$ と原点にいる確率 $\frac49$ が等しくない。戻る回数の期待値も $8$ を超えない |
| 条件 (i) の「$+1$ か $-1$」(その場にとどまる動きを許す。各動きの確率も $\frac13$ に変わる) | ex-rwr-cx-lazy:$-1,0,+1$ を $\frac13$ ずつ、$n=1$ | 戻らない確率 $\frac49$ と原点にいる確率 $\frac13$ が等しくない |
| 回数が偶数 $2n$ であること | ex-rwr-cx-odd:回数 $3$ | 3 回目まで戻らない確率は $\frac12$、3 回目に原点にいる確率は $0$ |
条件 (i) を変えて、$X_k=+1$ となる確率を $p=\dfrac23$、$X_k=-1$ となる確率を $q=\dfrac13$ とする(独立性はそのまま)。
(1) $n=1$(2 回)で比べる。2 回目まで一度も原点に戻らないのは $++$ と $--$ のときで、確率は $p^2+q^2=\dfrac49+\dfrac19=\dfrac59$ である。2 回目に原点にいるのは $+-$ と $-+$ のときで、確率は $2pq=\dfrac49$ である。$\dfrac59\ne\dfrac49$ なので、thm-rwr-main の (2) の等式は成り立たない。証明の段 3 では「どの歩みも同じ確率」を使って、確率を歩みの数に置きかえていた。$p\ne q$ では $+1$ の多い歩みほど起こりやすく、数えることと確率を求めることが一致しない。
(2) $2n$ 回目に原点にいる確率は、$+1$ と $-1$ が $n$ 回ずつの歩みがどれも確率 $p^nq^n$ で起こるので、$\dbinom{2n}n(pq)^n$ である。$\dbinom{2n}n$ は二項定理の $(1+1)^{2n}=4^n$ を展開した項の 1 つなので $4^n$ 以下で、
$$
\binom{2n}n(pq)^n\le(4pq)^n=\left(\frac89\right)^n
$$
である。cor-rwr-count の段 1 と同じ計算で、$2n$ 回目までに原点に戻る回数の期待値は
$$
\sum_{k=1}^n\binom{2k}k(pq)^k\le\sum_{k=1}^n\left(\frac89\right)^k<\frac{8/9}{1-8/9}=8
$$
で、$n$ をどれほど大きくしても $8$ を超えない。$p=\dfrac12$ のときに限りなく大きくなったこと(cor-rwr-count)と対照的である。$p\ne\dfrac12$ では、いつか原点に戻る確率も $1$ より小さくなることが知られているが、この記事では証明しない。
条件 (i) を変えて、1 回の動きが $-1$、$0$、$+1$ のどれかで、それぞれ確率 $\dfrac13$ ずつとする(独立性はそのまま)。この点は奇数回目にも原点にいることがある(1 回目に $0$ が出れば $S_1=0$)ので、prop-rwr-parity も成り立たない。
2 回目までで比べる。2 回目に原点にいるのは、動きが $(0,0)$、$(+1,-1)$、$(-1,+1)$ の 3 通りで、確率は $\dfrac39=\dfrac13$ である。2 回目まで一度も原点に戻らないのは、1 回目が $\pm1$(2 通り)で、2 回目に原点に戻らない動き(位置 $1$ からなら $0$ か $+1$、位置 $-1$ からなら $0$ か $-1$ の 2 通り)のときで、$2\cdot2=4$ 通り、確率は $\dfrac49$ である。$\dfrac49\ne\dfrac13$ なので、等式は成り立たない。
def-rwr-walk のランダムウォークのまま、回数を奇数 $3$ にする。prop-rwr-parity より $S_3\ne0$ がつねに成り立つので、「3 回目まで一度も原点に戻らない」は「2 回目まで一度も原点に戻らない」と同じ出来事で、確率は thm-rwr-main の (2) により $u_2=\dfrac12$ である。一方、3 回目に原点にいる確率は $0$ である。等式を奇数の回数 $2n+1$ で使うときは、「$2n+1$ 回目まで一度も原点に戻らない確率は $u_{2n}$」と読みかえる。
GS06 の第 12.1 節は、$2m$ 回目に初めて原点に戻る確率 $f_{2m}$ と $u_{2m}$ の関係(Theorem 12.2)を母関数の式に直して、$f_{2m}$ を求め(Theorem 12.3)、その母関数に $x=1$ を代入して、原点に戻る確率が $1$ であることを示している(Example 12.1)。この記事は母関数を使わず、鏡像原理(投票定理)と、積の上下の評価だけで同じ結論を得た。$f_{2m}$ は 初めて原点に戻る時刻 で、母関数を使わずに差をとる方法で求める。
平面の格子の上を、上下左右に確率 $\frac14$ ずつ動く点や、空間の格子の上を 6 方向に確率 $\frac16$ ずつ動く点についても、同じ問いを考えられる。GS06 の Example 12.2 は、平面では原点に戻る確率が $1$ であり、空間では $1$ より小さいことを示している(GS06 はこれらの結果を Pólya によるものとしている)。この記事では証明しない。数直線と平面では「いつか必ず戻る」のに、空間では戻らないことがある、という違いは、ランダムウォーク で扱う話題である。
長さ $8$ の歩みで、$S_1,\ldots,S_8$ がすべて正のものを、最後の位置 $S_8=2j$($j=1,2,3,4$)ごとに thm-rwr-main の証明の段 3 の式で数え、原点に戻らない歩みの数が $\dbinom84$ になることを確かめよ。
段 3 の式 $\frac j4\binom8{4+j}$ で、$j=1$:$\frac14\cdot56=14$、$j=2$:$\frac24\cdot28=14$、$j=3$:$\frac34\cdot8=6$、$j=4$:$\frac44\cdot1=1$。合計 $14+14+6+1=35$ で、ずっと負のものも $35$ 通りなので、原点に戻らない歩みは $70$ 通り。$\binom84=70$ と一致し、確率は $\frac{70}{256}=\frac{35}{128}=u_8$ である。
$2n$ 回目までに原点に戻る確率が $0.9$ 以上になる最小の $n$ を求めたい。(1) lem-rwr-bounds を使って、$n\ge50$ なら確率が $0.9$ 以上であること、$n\le24$ なら $0.9$ 未満であることを示せ。(2) $u_{2n}$ を計算して、最小の $n$ を求めよ。
(1) 確率は $1-u_{2n}$ なので、$u_{2n}\le0.1$ かどうかを調べればよい。$n\ge50$ なら $2n+1\ge101>100$ で、$u_{2n}\le\frac1{\sqrt{2n+1}}<\frac1{10}$。$n\le24$ なら $u_{2n}\ge\frac1{2\sqrt n}\ge\frac1{2\sqrt{24}}>\frac1{2\cdot5}=0.1$($\sqrt{24}<5$)。
(2) prop-rwr-recur の漸化式で順に計算すると、$u_{62}\approx0.1009$($n=31$)、$u_{64}\approx0.0993$($n=32$)。よって最小の $n$ は $32$ で、$64$ 回投げれば、原点に戻る確率がはじめて $0.9$ を超える。
Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する