acpass

モデリング用語

ひとことで言うと

ある状態に「戻ってこられる歩数」が何歩おきかを表すのが周期です。戻れる歩数の最大公約数で定義し、それが1なら非周期、2以上ならその間隔でしか戻れない規則的な連鎖だ、ということになります。非周期であることが、分布が一定値に収束する(極限分布が存在する)ための鍵になります。

状態0に自己ループP_00=0.7>0を持つ2状態連鎖。戻れるステップ数の集合が{1,2,3,...}となり最大公約数d(0)=1で非周期。これに対し自己ループのない交互推移ではd=2の周期をもつ。周期は戻りうるステップ数のgcdで判定する非周期例 P=[[0.7,0.3],[0.4,0.6]]:自己ループ有りで d=1010.30.40.70.6P_00=0.7>0 なので戻れる歩数は{1,2,...}、d(0)=gcd=1(非周期)対照:自己ループ無しの交互推移は戻れる歩数が偶数のみで d=2

非周期な例 P=\begin{psmallmatrix}0.7&0.3\\0.4&0.6\end{psmallmatrix}。状態0に自己ループ P00=0.7>0P_{00}=0.7>0 があるので、戻れる歩数の集合は {1,2,3,}\{1,2,3,\dots\}、最大公約数 d(0)=1d(0)=1 で非周期。自己ループの無い交互推移なら戻れる歩数は偶数のみで周期 d=2d=2 になる。

数式で表すと

d(i)=gcd{n:Piin>0}d(i)=\gcd\{n: P^n_{ii}>0\}

状態に戻りうるステップ数の最大公約数。1なら非周期で極限分布が存在。

状態 ii の周期は、その状態に正の確率で戻ってこられるステップ数すべての最大公約数 d(i)=gcd{n1: Piin>0}d(i)=\gcd\{n\ge1:\ P^n_{ii}>0\} で定義されます。d(i)=1d(i)=1 なら非周期、d(i)2d(i)\ge2 ならその値が周期で、d(i)d(i) の倍数の歩数でしか元に戻れない規則性を意味します。周期はconcept: 極限分布・concept: エルゴード性を理解するうえでの判定道具で、「非周期かどうか」が極限分布の存在を左右します。 計算の手順を、2つの例で具体的に追います。まず周期的な例 P=(0110)P=\begin{pmatrix}0&1\\1&0\end{pmatrix}(状態0と1を確率1で交互に推移)。状態0に戻るには必ず偶数歩かかります(P00=0P_{00}=0 なので1歩では戻れない、P002=1P^2_{00}=1 なので2歩では戻れる)。戻れる歩数の集合は {2,4,6,}\{2,4,6,\dots\}、最大公約数は d(0)=gcd{2,4,6,}=2d(0)=\gcd\{2,4,6,\dots\}=2。よってこの連鎖は周期2で、非周期ではありません。これがconcept: 極限分布で PnP^nIIPP を往復して収束しなかった原因そのものです。 次に非周期な例 P=(0.70.30.40.6)P=\begin{pmatrix}0.7&0.3\\0.4&0.6\end{pmatrix}。状態0には自己ループ P00=0.7>0P_{00}=0.7>0 があるので、1歩で戻れます。さらに2歩でも3歩でも戻れるので、戻れる歩数の集合は {1,2,3,}\{1,2,3,\dots\}、最大公約数は d(0)=gcd{1,2,3,}=1d(0)=\gcd\{1,2,3,\dots\}=1 で非周期です。判定の実用的なコツは、自己ループ Pii>0P_{ii}>0 が1つでもあれば即座に d(i)=1d(i)=1(非周期)が確定する、という点です。一般に、既約な連鎖では全状態の周期が一致するので、代表として1状態を調べれば足ります。

試験に出る性質

定義

d(i)=gcd{n1:Piin>0}d(i)=\gcd\{n\ge1:P^n_{ii}>0\}。状態 ii に戻れる歩数すべての最大公約数。

非周期の判定

d(i)=1d(i)=1 なら非周期、d(i)2d(i)\ge2 ならその値が周期で倍数歩でしか戻れない。

自己ループは即非周期

Pii>0P_{ii}>0 が1つでもあれば1歩で戻れるので d(i)=1d(i)=1 が確定する。実用的な近道。

既約なら周期は共通

既約な連鎖では全状態の周期が一致するので、1状態を調べれば全体が分かる。

極限分布との関係

非周期(d=1d=1)なら極限分布が存在し定常分布に収束。周期 2\ge2 なら振動して極限分布なし(concept: 極限分布)。

例で見る

周期例 P=(0110)P=\begin{pmatrix}0&1\\1&0\end{pmatrix}:状態0に戻れるのは偶数歩のみ(P002=1,P004=1,P^2_{00}=1,P^4_{00}=1,\dots、奇数歩は0)。 d(0)=gcd{2,4,6,}=2d(0)=\gcd\{2,4,6,\dots\}=2 で周期2(非周期ではない)。 非周期例 P=(0.70.30.40.6)P=\begin{pmatrix}0.7&0.3\\0.4&0.6\end{pmatrix}:自己ループ P00=0.7>0P_{00}=0.7>0 なので1歩で戻れる。 戻れる歩数は {1,2,3,}\{1,2,3,\dots\}d(0)=gcd{1,2,3,}=1d(0)=\gcd\{1,2,3,\dots\}=1(非周期)。d=1d=1 なので極限分布が存在する。

つまずきポイント

  • 周期を「戻れる最小歩数」と勘違いする(正しくは戻れる歩数すべての最大公約数。最小ではない)
  • 自己ループがあるのに周期2以上と判定する(Pii>0P_{ii}>0 なら1歩で戻れて必ず d(i)=1d(i)=1=非周期)
  • 周期と極限分布の存在を切り離して考える(周期1でないと極限分布は存在しない。判定の目的を見失わない)

定着クイズ

状態 ii の周期 d(i)d(i) の定義は?

自己ループ Pii>0P_{ii}>0 がある状態の周期は?

P=[[0,1],[1,0]]P=[[0,1],[1,0]] の状態0の周期 d(0)d(0) は?

この用語を扱う問題(1