Probability concept
連続 時間マルコフ連鎖 ( CTMC )は、各状態において、 指数確率変数 に従って状態が変化し、その後、 確率行列 の確率によって指定された別の状態に移行する連続 確率過程 です。同等の定式化では、この過程は、遷移可能な状態ごとに1つずつ、一連の指数確率変数の最小値に従って状態が変化すると記述され、パラメータは現在の状態によって決定されます
3つの状態を持つCTMCの例は次のとおりです。プロセスは 、ホールド時間 (指数確率変数) で指定された時間後に遷移します。 ここで、 i は現在の状態です。各確率変数は独立しており 、、、 となります。遷移が発生する場合、プロセスは ジャンプ チェーン (確率行列を持つ 離散時間マルコフ連鎖 )に従って移動します { 0 , 1 , 2 } {\displaystyle \{0,1,2\}} E i {\displaystyle E_{i}} E 0 ∼ Exp ( 6 ) {\displaystyle E_{0}\sim {\text{Exp}}(6)} E 1 ∼ Exp ( 12 ) {\displaystyle E_{1}\sim {\text{Exp}}(12)} E 2 ∼ Exp ( 18 ) {\displaystyle E_{2}\sim {\text{Exp}}(18)}
[ 0 1 2 1 2 1 3 0 2 3 5 6 1 6 0 ] . {\displaystyle {\begin{bmatrix}0&{\frac {1}{2}}&{\frac {1}{2}}\\{\frac {1}{3}}&0&{\frac {2}{3}}\\{\frac {5}{6}}&{\frac {1}{6}}&0\end{bmatrix}}.} 同様に、競合指数関数 の性質により 、このCTMCは 、独立で、かつ、 パラメータが Q行列 によって与えられる2つの確率変数の最小値に従って状態 iから状態を変化させます。 E i , j ∼ Exp ( q i , j ) {\displaystyle E_{i,j}\sim {\text{Exp}}(q_{i,j})} i ≠ j {\displaystyle i\neq j} Q = ( q i , j ) {\displaystyle Q=(q_{i,j})}
[ − 6 3 3 4 − 12 8 15 3 − 18 ] . {\displaystyle {\begin{bmatrix}-6&3&3\\4&-12&8\\15&3&-18\end{bmatrix}}.} 非対角成分はそれぞれ、ジャンプ連鎖が状態 iから状態 j に 移動する確率を、状態 i の期待保持時間で割ったものとして計算できます 。対角成分は、各行の合計が0になるように選択されます。 q i , j {\displaystyle q_{i,j}}
CTMCは、指数分布と離散時間マルコフ連鎖の無記憶性により、その動作は現在の状態のみ 依存し、過去の動作には依存しないというマルコフ 性を満たします。
定義 を確率空間とし、を可算な空でない集合とし 、 (「時間」)とします 。 離散 計量 をとすることで、 関数 の 正しい連続性 を理解できます 。連続時間マルコフ連鎖は次のように定義されます。 [1] ( Ω , A , Pr ) {\displaystyle (\Omega ,{\cal {A}},\Pr )} S {\displaystyle S} T = R ≥ 0 {\displaystyle T=\mathbb {R} _{\geq 0}} T {\displaystyle T} S {\displaystyle S} R ≥ 0 → S {\displaystyle \mathbb {R} _{\geq 0}\to S}
上の 確率ベクトル (以下では 、マルコフ連鎖の 初期分布 として解釈します)、および λ {\displaystyle \lambda } S {\displaystyle S} 上の 速度行列 、つまり、 となる 関数 Q {\displaystyle Q} S {\displaystyle S} Q : S 2 → R {\displaystyle Q:S^{2}\to \mathbb {R} } すべての異なる に対して 、 i , j ∈ S , 0 ≤ q i , j {\displaystyle i,j\in S,0\leq q_{i,j}} すべての に対して となります ( が無限大であっても、この和は 事前に 明確に定義されます( に等しい可能性があります )。なぜなら、和に現れる各項は非負だからです。 事後的には 、が に等しいと仮定し、 が実数値である と仮定しているので 、和も有限である( に等しくない) ことがわかります。 一部の著者は、 という修正された規定を除いて逐語的に同じ定義を使用し 、が 安定して いる、 または 完全に安定している と言って、すべてのエントリが実数値であることを 意味します。) [2] i ∈ S , {\displaystyle i\in S,} ∑ j ∈ S : j ≠ i q i , j = − q i , i . {\displaystyle \sum _{j\in S:j\neq i}q_{i,j}=-q_{i,i}.} S {\displaystyle S} + ∞ {\displaystyle +\infty } + ∞ {\displaystyle +\infty } − q i , i {\displaystyle -q_{i,i}} Q {\displaystyle Q} Q : S 2 → R ∪ { − ∞ } {\displaystyle Q:S^{2}\to \mathbb {R} \cup \{-\infty \}} Q {\displaystyle Q} range Q ⊆ R {\displaystyle \operatorname {range} Q\subseteq \mathbb {R} } の行和は 0: 、より簡潔に言えば であることに注意してください。この状況は 、遷移行列のすべての行和が1に等しい 離散時間マルコフ連鎖 の状況とは対照的です Q {\displaystyle Q} ∀ i ∈ S , ∑ j ∈ S q i , j = 0 , {\displaystyle \forall i\in S,~\sum _{j\in S}q_{i,j}=0,} Q ⋅ 1 = 0 {\displaystyle Q\cdot 1=0}
ここで、 が-可測で あるとします。 初期分布 と速度行列 を用いてマルコフであることを定義するには 、遷移確率を用いるか、ジャンプ連鎖と保持時間を用いるかの 3つの同等な方法があります。 [5] X : T → S Ω {\displaystyle X:T\to S^{\Omega }} ∀ t ∈ T X ( t ) {\displaystyle \forall t\in T~X(t)} ( A , P ( S ) ) {\displaystyle ({\cal {A}},{\cal {P}}(S))} X {\displaystyle X} λ {\displaystyle \lambda } Q {\displaystyle Q}
遷移確率定義の前提として、まず 正規 速度行列 の定義を説明します。遷移速度行列を用いて 、次の定理により ( ) 上の 遷移行列 の集合を生成することで、マルコフ連鎖のダイナミクスを指定します。 Q {\displaystyle Q} P ( t ) {\displaystyle P(t)} S {\displaystyle S} t ∈ R ≥ 0 {\displaystyle t\in \mathbb {R} _{\geq 0}}
コルモゴロフ後退方程式 ( [6] ) の解の存在 — すべての に対して、 要素 が 微分可能であり、 コルモゴロフ後退方程式を 満たす が 存在する 。 P ∈ ( [ 0 , 1 ] S × S ) T {\displaystyle P\in ([0,1]^{S\times S})^{T}} i , j ∈ S {\displaystyle i,j\in S} ( P ( t ) i , j ) t ∈ T {\displaystyle (P(t)_{i,j})_{t\in T}} P {\displaystyle P}
P ( 0 ) = ( [ i = j ] ) i , j ∈ S , ∀ t ∈ T ∀ i , j ∈ S ( P ( t ) i , j ) ′ = ∑ k ∈ S q i , k P ( t ) k , j . {\displaystyle P(0)=([i=j])_{i,j\in S},~\forall t\in T~\forall i,j\in S~~(P(t)_{i,j})'=\sum _{k\in S}q_{i,k}P(t)_{k,j}.} 0
が正則で ある と言うのは 、上記の系が一意であること、すなわち解が1つだけ存在することを意味します。 [7] 不規則 である と言うのは、が正則ではないこと を意味します 。 が有限である場合、解は1つだけ存在し、 したがって 正則です。そうでない場合、 は無限であり、上に不規則な遷移率行列が存在します 。 [a] が正則である場合 、唯一の解 に対して 、各 に対して 、は 確率行列 になります 。 [6] 次のサブセクションの最初からこのセクションの最後まで、 が正則であると 仮定しますが、この仮定を含めないのが慣例です 。(専門家への注意:したがって、一般的な連続時間マルコフ連鎖を定義するのではなく、 非爆発的な 連続時間マルコフ連鎖のみを定義します。) Q {\displaystyle Q} Q {\displaystyle Q} Q {\displaystyle Q} S {\displaystyle S} P = ( e t Q ) t ∈ T , {\displaystyle P=(e^{tQ})_{t\in T},} Q {\displaystyle Q} S {\displaystyle S} S {\displaystyle S} Q {\displaystyle Q} P {\displaystyle P} t ∈ T {\displaystyle t\in T} P ( t ) {\displaystyle P(t)} Q {\displaystyle Q}
遷移確率の定義 をシステム( 0 )の(唯一の)解とします 。(唯一のことは、が正則であるという仮定によって保証されます 。) 初期分布 と速度行列 を持つマルコフとは、任意の非負整数に対して 、 すべてのに対して となること を意味します P {\displaystyle P} Q {\displaystyle Q} X {\displaystyle X} λ {\displaystyle \lambda } Q {\displaystyle Q} n ≥ 0 {\displaystyle n\geq 0} t 0 , … , t n + 1 ∈ T {\displaystyle t_{0},\dots ,t_{n+1}\in T} t 0 < ⋯ < t n + 1 , {\displaystyle t_{0}<\dots <t_{n+1},} i 0 , … , i n + 1 ∈ I , {\displaystyle i_{0},\dots ,i_{n+1}\in I,}
Pr ( X 0 = i 0 , … , X t n + 1 = i n + 1 ) = λ i 0 ∏ k ∈ Z : 0 ≤ k ≤ n P ( t k + 1 − t k ) i k , i k + 1 {\displaystyle \Pr(X_{0}=i_{0},\dots ,X_{t_{n+1}}=i_{n+1})=\lambda _{i_{0}}\prod _{k\in \mathbb {Z} :0\leq k\leq n}P(t_{k+1}-t_{k})_{i_{k},i_{k+1}}} 。 1
帰納法と、 ( 1 )を含む上記の命題と次の命題 との同値性を示すことができるという事実を用いて、すべての に対して 、任意の非負整数 に対して 、すべての に対して、 すべての に対して ( したがって となる )、 ∀ A , B ∈ A Pr ( B ) ≠ 0 → Pr ( A ∩ B ) = Pr ( A ∣ B ) Pr ( B ) , {\displaystyle \forall A,B\in {\cal {A}}~~\Pr(B)\neq 0\rightarrow \Pr(A\cap B)=\Pr(A\mid B)\Pr(B),} i ∈ I , Pr ( X 0 = i ) = λ i {\displaystyle i\in I,~\Pr(X_{0}=i)=\lambda _{i}} n ≥ 0 {\displaystyle n\geq 0} t 0 , … , t n + 1 ∈ T {\displaystyle t_{0},\dots ,t_{n+1}\in T} t 0 < ⋯ < t n + 1 , {\displaystyle t_{0}<\dots <t_{n+1},} i 0 , … , i n + 1 ∈ I {\displaystyle i_{0},\dots ,i_{n+1}\in I} 0 < Pr ( X 0 = i 0 , … , X t n = i n ) {\displaystyle 0<\Pr(X_{0}=i_{0},\dots ,X_{t_{n}}=i_{n})} 0 < Pr ( X t n = i n ) {\displaystyle 0<\Pr(X_{t_{n}}=i_{n})}
Pr ( X t n + 1 = i n + 1 ∣ X t n = i n , … , X t 0 = i 0 ) = P ( t n + 1 − t n ) i n , i n + 1 . {\displaystyle \Pr(X_{t_{n+1}}=i_{n+1}\mid X_{t_{n}}=i_{n},\dots ,X_{t_{0}}=i_{0})=P(t_{n+1}-t_{n})_{i_{n},i_{n+1}}.} 2
関数 ( ) の連続性から、軌道は( 上の 離散計量 に関して) ほぼ確実に 右連続で あることが分かる。つまり、 となる -零集合が 存在する 。 ( P ( t ) i , j ) t ∈ T {\displaystyle (P(t)_{i,j})_{t\in T}} i , j ∈ S {\displaystyle i,j\in S} ( X t ( ω ) ) t ∈ T {\displaystyle (X_{t}(\omega ))_{t\in T}} S {\displaystyle S} Pr {\displaystyle \Pr } N {\displaystyle N} { ω ∈ Ω : ( X t ( ω ) ) t ∈ T is right continuous } ⊆ N {\displaystyle \{\omega \in \Omega :(X_{t}(\omega ))_{t\in T}{\text{ is right continuous}}\}\subseteq N}
ジャンプ連鎖/保持時間の定義
右連続関数に関連付けられたシーケンス を右連続と仮定する(離散計量 を と仮定した場合 ) 。 定義 する f : T → S {\displaystyle f:T\to S} S {\displaystyle S}
h = h ( f ) = ( inf { u ∈ ( 0 , + ∞ ) : f ( t + u ) ≠ f ( t ) } ) t ∈ T ) ∪ { + ∞ , 0 } , {\displaystyle h=h(f)=(\inf\{u\in (0,+\infty ):f(t+u)\neq f(t)\})_{t\in T})\cup \{+\infty ,0\},} とする
H = H ( f ) = ( h ∘ n 0 ) n ∈ Z ≥ 0 {\displaystyle H=H(f)=(h^{\circ n}0)_{n\in \mathbb {Z} _{\geq 0}}} を に関連付けられた 保持時間シーケンス とし 、 を選択し 、 とする f {\displaystyle f} s ∈ S , {\displaystyle s\in S,}
y = y ( f ) = ( { f ( ∑ k ∈ n H k ) if ∑ k ∈ n H k < + ∞ , s else ) n ∈ ω {\displaystyle y=y(f)=\left({\begin{cases}f(\sum _{k\in n}H_{k})&{\text{ if }}\sum _{k\in n}H_{k}<+\infty ,\\s&{\text{ else}}\end{cases}}\right)_{n\in \omega }} を に関連付けられた「 状態シーケンス 」とする 。 f {\displaystyle f}
ジャンプ行列 Π の定義 ジャンプ 行列 は 、 への依存性を強調したい場合は、 と表記され、 が 関数の 零点集合 で ある 行列である Π {\displaystyle \Pi } Π ( Q ) {\displaystyle \Pi (Q)} Q {\displaystyle Q} Π = ( [ i = j ] ) i ∈ Z , j ∈ S ∪ ⋃ i ∈ S ∖ Z ( { ( ( i , j ) , − q i , j / q i , i ) : j ∈ S ∖ { i } } ∪ { ( ( i , i ) , 0 ) } ) , {\displaystyle \Pi =([i=j])_{i\in Z,j\in S}\cup \bigcup _{i\in S\setminus Z}(\{((i,j),-q_{i,j}/q_{i,i}):j\in S\setminus \{i\}\}\cup \{((i,i),0)\}),} Z = Z ( Q ) = { k ∈ S : q k , k = 0 } {\displaystyle Z=Z(Q)=\{k\in S:q_{k,k}=0\}} ( q k , k ) k ∈ S . {\displaystyle (q_{k,k})_{k\in S}.}
ジャンプ連鎖/保持時間特性 は初期分布 と速度行列 を持つマルコフ連鎖である とは 、次のことを意味します。の軌跡は ほぼ確実に右連続であり、 を(どこでも)右連続の軌跡を持つように 修正すると、 ほぼ確実に(専門家への注意:この条件は 非爆発的であると言います)、状態系列 は初期分布 (ジャンプ連鎖特性)と遷移行列 と (保持時間特性)を持つ離散時間マルコフ連鎖です。 X {\displaystyle X} λ {\displaystyle \lambda } Q {\displaystyle Q} X {\displaystyle X} f {\displaystyle f} X {\displaystyle X} ∑ n ∈ Z ≥ 0 H ( f ( ω ) ) n = + ∞ {\displaystyle \sum _{n\in \mathbb {Z} _{\geq 0}}H(f(\omega ))_{n}=+\infty } X {\displaystyle X} y ( f ( ω ) ) {\displaystyle y(f(\omega ))} λ {\displaystyle \lambda } Π ( Q ) , {\displaystyle \Pi (Q),} ∀ n ∈ Z ≥ 0 ∀ B ∈ B ( R ≥ 0 ) Pr ( H n ( f ) ∈ B ) = Exp ( − q Y n , Y n ) ( B ) {\displaystyle \forall n\in \mathbb {Z} _{\geq 0}~\forall B\in {\cal {B}}(\mathbb {R} _{\geq 0})~\Pr(H_{n}(f)\in B)=\operatorname {Exp} (-q_{Y_{n},Y_{n}})(B)}
無限小定義 連続時間マルコフ連鎖は、遷移速度、つまり状態iとjの間の遷移確率の時間に関する微分によって特徴付けられます。 は 初期分布 と速度行列 を持つマルコフ連鎖であるとは 、すべての およびすべてのに対して 、すべての に対して 、そして の小さな正の値に対して 、すべての に対して以下が成り立つ ことを意味します 。 X {\displaystyle X} λ {\displaystyle \lambda } Q {\displaystyle Q} i ∈ S , {\displaystyle i\in S,} Pr ( X ( 0 ) = i ) = λ i {\displaystyle \Pr(X(0)=i)=\lambda _{i}} i , j {\displaystyle i,j} t {\displaystyle t} h {\displaystyle h} t ∈ T {\displaystyle t\in T} 0 < Pr ( X ( t ) = i ) {\displaystyle 0<\Pr(X(t)=i)}
Pr ( X ( t + h ) = j ∣ X ( t ) = i ) = [ i = j ] + q i , j h + o ( h ) {\displaystyle \Pr(X(t+h)=j\mid X(t)=i)=[i=j]+q_{i,j}h+o(h)} 、 ここで、項 は の場合 、それ以外の 場合は であり、 小文字のo 項は に特定の方法で依存します 。 [ i = j ] {\displaystyle [i=j]} 1 {\displaystyle 1} i = j {\displaystyle i=j} 0 {\displaystyle 0} o ( h ) {\displaystyle o(h)} i , j , h {\displaystyle i,j,h}
上記の式は、が に対して から への遷移がどれだけ速く起こり 、 に対して から への遷移がどれだけ速く起こるかを測るものと見ることができる こと を示しています 。 q i , j {\displaystyle q_{i,j}} i {\displaystyle i} j {\displaystyle j} i ≠ j {\displaystyle i\neq j} i {\displaystyle i} i = j {\displaystyle i=j}
特性
伝達クラス 伝達クラス、過渡性、再帰性、正再帰性および零再帰性は、 離散時間マルコフ連鎖 の場合と同じように定義されます。
過渡的動作 要素p ij = P( X t = j | X 0 = i )を持つ行列を P( t ) と書きます。すると、行列 P( t ) は順方向方程式、つまり 1階微分方程式を満たします。
P ′ ( t ) = P ( t ) Q {\displaystyle P'(t)=P(t)Q} 、 ここで、プライムはt に関する微分を表します 。この方程式の解は、 指数行列で与えられます。
P ( t ) = e t Q {\displaystyle P(t)=e^{tQ}} 。 状態空間 {1,2} 上の CTMC のような単純なケースでは、このようなプロセスの一般的な Q行列は、 α 、 β > 0 の次の 2 × 2 行列です。
Q = ( − α α β − β ) . {\displaystyle Q={\begin{pmatrix}-\alpha &\alpha \\\beta &-\beta \end{pmatrix}}.} この場合、順方向行列に関する上記の関係は、次のように明示的に解くことができます
P ( t ) = ( β α + β + α α + β e − ( α + β ) t α α + β − α α + β e − ( α + β ) t β α + β − β α + β e − ( α + β ) t α α + β + β α + β e − ( α + β ) t ) {\displaystyle P(t)={\begin{pmatrix}{\frac {\beta }{\alpha +\beta }}+{\frac {\alpha }{\alpha +\beta }}e^{-(\alpha +\beta )t}&{\frac {\alpha }{\alpha +\beta }}-{\frac {\alpha }{\alpha +\beta }}e^{-(\alpha +\beta )t}\\{\frac {\beta }{\alpha +\beta }}-{\frac {\beta }{\alpha +\beta }}e^{-(\alpha +\beta )t}&{\frac {\alpha }{\alpha +\beta }}+{\frac {\beta }{\alpha +\beta }}e^{-(\alpha +\beta )t}\end{pmatrix}}} 。 大きな行列では、直接解の計算は複雑です。Q が 行列の 半群 の生成元であるという事実が
P ( t + s ) = e ( t + s ) Q = e t Q e s Q = P ( t ) P ( s ) {\displaystyle P(t+s)=e^{(t+s)Q}=e^{tQ}e^{sQ}=P(t)P(s)} 使用されます。
定常分布 既約回帰CTMCの定常分布は、 t の大きな値に対して過程が収束する確率分布です 。先に検討した2状態過程において、P( t )が次のように与えられる
ことに注目してください。
P ( t ) = ( β α + β + α α + β e − ( α + β ) t α α + β − α α + β e − ( α + β ) t β α + β − β α + β e − ( α + β ) t α α + β + β α + β e − ( α + β ) t ) {\displaystyle P(t)={\begin{pmatrix}{\frac {\beta }{\alpha +\beta }}+{\frac {\alpha }{\alpha +\beta }}e^{-(\alpha +\beta )t}&{\frac {\alpha }{\alpha +\beta }}-{\frac {\alpha }{\alpha +\beta }}e^{-(\alpha +\beta )t}\\{\frac {\beta }{\alpha +\beta }}-{\frac {\beta }{\alpha +\beta }}e^{-(\alpha +\beta )t}&{\frac {\alpha }{\alpha +\beta }}+{\frac {\beta }{\alpha +\beta }}e^{-(\alpha +\beta )t}\end{pmatrix}}} 、 t → ∞ になると、 分布は
P π = ( β α + β α α + β β α + β α α + β ) {\displaystyle P_{\pi }={\begin{pmatrix}{\frac {\beta }{\alpha +\beta }}&{\frac {\alpha }{\alpha +\beta }}\\{\frac {\beta }{\alpha +\beta }}&{\frac {\alpha }{\alpha +\beta }}\end{pmatrix}}} 。 各行は同じ分布を持つことに注目してください。これは開始状態に依存しないためです。行ベクトル π は、
π Q = 0 {\displaystyle \pi Q=0} 制約付きで
∑ i ∈ S π i = 1 {\displaystyle \sum _{i\in S}\pi _{i}=1} 。
例1 金融市場の状態を記述する連続時間マルコフ連鎖の有向グラフ表現(注:数字は架空のものです。) 右の図は、状態空間{強気市場、弱気市場、停滞市場}と 遷移率行列を持つ連続時間マルコフ連鎖を表しています。
Q = ( − 0.025 0.02 0.005 0.3 − 0.5 0.2 0.02 0.4 − 0.42 ) . {\displaystyle Q={\begin{pmatrix}-0.025&0.02&0.005\\0.3&-0.5&0.2\\0.02&0.4&-0.42\end{pmatrix}}.} この連鎖の定常分布は、要素の合計が1になるという制約の下で、を 解くことで見つけることができます π Q = 0 {\displaystyle \pi Q=0}
π = ( 0.885 0.071 0.044 ) . {\displaystyle \pi ={\begin{pmatrix}0.885&0.071&0.044\end{pmatrix}}.}
例2 遷移確率付き遷移グラフ。状態1、5、6、8の例を示します。状態2と8の間には双方向の秘密通路があります 右の図は、状態空間{1,2,3,4,5,6,7,8,9}を持つ パックマンを モデル化した離散時間マルコフ連鎖を表しています。プレイヤーはパックマンを操作して迷路を進み、パックドットを食べます。その間、ゴーストに追われます。便宜上、迷路は3x3の小さなグリッドとし、ゴーストは水平方向と垂直方向にランダムに移動します。状態2と状態8の間には、両方向に秘密の通路が設けられています。以下の遷移率行列では、確率0の要素は削除されます。
Q = ( − 1 1 2 1 2 1 4 − 1 1 4 1 4 1 4 1 2 − 1 1 2 1 3 − 1 1 3 1 3 1 4 1 4 − 1 1 4 1 4 1 3 1 3 − 1 1 3 1 2 − 1 1 2 1 4 1 4 1 4 − 1 1 4 1 2 1 2 − 1 ) {\displaystyle Q={\begin{pmatrix}-1&{\frac {1}{2}}&&{\frac {1}{2}}\\{\frac {1}{4}}&-1&{\frac {1}{4}}&&{\frac {1}{4}}&&&{\frac {1}{4}}\\&{\frac {1}{2}}&-1&&&{\frac {1}{2}}\\{\frac {1}{3}}&&&-1&{\frac {1}{3}}&&{\frac {1}{3}}\\&{\frac {1}{4}}&&{\frac {1}{4}}&-1&{\frac {1}{4}}&&{\frac {1}{4}}\\&&{\frac {1}{3}}&&{\frac {1}{3}}&-1&&&{\frac {1}{3}}\\&&&{\frac {1}{2}}&&&-1&{\frac {1}{2}}\\&{\frac {1}{4}}&&&{\frac {1}{4}}&&{\frac {1}{4}}&-1&{\frac {1}{4}}\\&&&&&{\frac {1}{2}}&&{\frac {1}{2}}&-1\end{pmatrix}}}
このマルコフ連鎖は既約ではありません。なぜなら、ゴーストは有限の時間内にあらゆる状態からあらゆる状態へ移動できるからです。秘密の通路があるため、マルコフ連鎖は非周期的でもあります。なぜなら、ゴーストは任意の状態から任意の状態へ、偶数回と奇数回の両方の状態遷移で移動できるからです。したがって、唯一の定常分布が存在し、 要素の合計が1になるという制約の下で を解くことで見つけることができます。制約の下でのこの線形方程式の解は です。 中央の状態と隣接する秘密の通路の境界状態2と8が最も多く訪れられ、コーナーの状態は最も少なく訪れられます。 π Q = 0 {\displaystyle \pi Q=0} π = ( 7.7 , 15.4 , 7.7 , 11.5 , 15.4 , 11.5 , 7.7 , 15.4 , 7.7 ) % . {\displaystyle \pi =(7.7,15.4,7.7,11.5,15.4,11.5,7.7,15.4,7.7)\%.}
時間反転 CTMC X t の場合、時間反転プロセスは と定義されます 。 ケリーの補題 により、このプロセスは順方向プロセスと同じ定常分布を持ちます X ^ t = X T − t {\displaystyle {\hat {X}}_{t}=X_{T-t}}
連鎖は、逆の過程が順方向の過程と同じである場合、可逆的であると言われます。 コルモゴロフの基準 によれば、過程が可逆的であるための必要十分条件は、閉ループの周りの遷移率の積が両方向で同じでなければならないということです。
埋め込みマルコフ連鎖 エルゴード連続時間マルコフ連鎖 Q の 定常 確率分布 π を 見つける 1つの方法は 、まず 埋め込みマルコフ連鎖(EMC) を見つけることです。厳密に言えば、EMCは正則離散時間マルコフ連鎖です。EMCの1ステップ遷移確率行列 Sの各要素は s ij で表され 、 状態 iから状態 j への遷移の 条件付き確率 を表します。これらの条件付き確率は次のように見つけることができます。
s i j = { q i j ∑ k ≠ i q i k if i ≠ j , 0 otherwise . {\displaystyle s_{ij}={\begin{cases}{\frac {q_{ij}}{\sum _{k\neq i}q_{ik}}}&{\text{if }}i\neq j,\\0&{\text{otherwise}}.\end{cases}}} このことから、 S は 次のように表すことができます
S = I − ( diag ( Q ) ) − 1 Q {\displaystyle S=I-\left(\operatorname {diag} (Q)\right)^{-1}Q} ここで、 Iは 単位行列 であり 、diag( Q ) は 行列 Qから 主対角要素 を選択し 、他のすべての要素をゼロに設定すること
によって形成される 対角行列です。
定常確率分布ベクトルを求めるには、次に、 となるようなもの
を見つけなければなりません。 φ {\displaystyle \varphi }
φ S = φ , {\displaystyle \varphi S=\varphi ,} は行ベクトルで あり、すべての要素が 0より大きく、 = 1となるようなもの。これから、 πは 次のように求められます
。 φ {\displaystyle \varphi } φ {\displaystyle \varphi } ‖ φ ‖ 1 {\displaystyle \|\varphi \|_{1}}
π = − φ ( diag ( Q ) ) − 1 ‖ φ ( diag ( Q ) ) − 1 ‖ 1 . {\displaystyle \pi ={-\varphi (\operatorname {diag} (Q))^{-1} \over \left\|\varphi (\operatorname {diag} (Q))^{-1}\right\|_{1}}.} ( Sは 周期的であっても、 Qは 周期的でなくても構いません。πが求められたら 、 単位ベクトル に正規化する必要があります 。)
連続時間マルコフ連鎖から導出できるもう1つの離散時間プロセスは、δスケルトンです。これは、 X ( t )をδ単位の時間間隔で観測することによって形成される(離散時間)マルコフ連鎖です。確率変数 X (0), X (δ), X (2δ), ...は、δスケルトンが訪れる状態の系列を与えます。
参照
注釈 ^ Ross, SM (2010). Introduction to Probability Models (10 ed.). Elsevier. ISBN 978-0-12-375686-2 。 ^ Anderson 1991, 64ページの定義を参照。 ^ Norris 1997、定理2.8.4および定理2.8.2(b) ^ ab Anderson 1991、定理2.2.2(1)、70ページ ^ Anderson 1991、定義81ページ ^ Anderson 1991、20ページ
参考文献 Anderson, William J. (1991). Continuous-time Markov chains: an applications-oriented approach . Springer Leo Breiman (1992) [1968] 確率論 。初版はAddison-Wesley社発行、 応用数学協会により再版。ISBN 0-89871-296-3 (第7章参照) Chen, Mu-Fa (2004). マルコフ連鎖から非平衡粒子系へ (第2版). World Scientific. Chen, Mu-Fa; Mao, Yong-Hua (2021). 確率過程入門 . World Scientific. JL Doob (1953) . 確率過程 . ニューヨーク:John Wiley and Sons ISBN 0-471-52369-0 。 AA Markov (1971). 「確率論の極限定理の連鎖的変数の和への拡張」。R. Howard著『 動的確率システム』第1巻:マルコフ連鎖 の付録Bに転載。John Wiley and Sons. Markov, AA (2006). 「連鎖的サンプルの連結に関するテキスト『エフゲニー・オネーギン』の統計的調査の例」。Science in Context . 19 (4). Link, Davidによる翻訳: 591–600 . doi :10.1017/s0269889706001074. S2CID 144854176 SP MeynとRL Tweedie (1993) 『マルコフ連鎖と確率的安定性 』ロンドン:Springer-Verlag ISBN 0-387-19832-6 オンライン:MCSS。第2版、ケンブリッジ大学出版局、2009年。 Kemeny, John G.、Hazleton Mirkil、J. Laurie Snell、Gerald L. Thompson (1959) 『有限数学構造 (第1版)』、ニュージャージー州エングルウッドクリフス:Prentice-Hall, Inc. 米国議会図書館カード目録番号 59-12841。 古典文献。第6章「有限マルコフ連鎖」 384ページ 以降を参照。 John G. Kemeny & J. Laurie Snell (1960) 『有限マルコフ連鎖』 、D. van Nostrand Company ISBN 0-442-04328-7 E. ヌメリン著『 既約マルコフ連鎖と非負値作用素 』ケンブリッジ大学出版局、1984年、2004年。ISBN 0-521-60494-X ノリス、JR (1997) 『マルコフ連鎖 』 doi :10.1017/CBO9780511810633.005。ISBN 9780511810633 。 セネタ、E. 『非負値行列とマルコフ連鎖 』第2版、1981年、XVI、288ページ、ソフトカバー、Springer Series in Statistics。(初版:Allen & Unwin Ltd.、ロンドン、1973年) ISBN 978-0-387-29765-1 ユーリ・スホフ、マーク・ケルバート(2008年) 『マルコフ連鎖:ランダム過程とその応用の入門 』ケンブリッジ大学出版局 ^ 例えば、 となるような(唯一の)遷移率行列 と を考えてみましょう 。 ( この 場合、 の残りの要素は すべてゼロになります。 誕生過程 を参照してください。)この場合 、 は不規則行列です。そして、一般の無限 の場合 、 非負整数でインデックスを付ける と、上記の行列を適切に修正したバージョンは 不規則行列になります。 [9] S = Z ≥ 0 {\displaystyle S=\mathbb {Z} _{\geq 0}} Q {\displaystyle Q} S {\displaystyle S} ∀ i ∈ Z ≥ 0 Q i , i + 1 = i 2 , Q i , i = − i 2 {\displaystyle \forall i\in \mathbb {Z} _{\geq 0}~~Q_{i,i+1}=i^{2},~Q_{i,i}=-i^{2}} Q {\displaystyle Q} Q {\displaystyle Q} S {\displaystyle S} S {\displaystyle S} Z ≥ 0 {\displaystyle \mathbb {Z} _{\geq 0}} Q {\displaystyle Q}
単一待ち行列ノード 到着過程 待ち行列ネットワーク サービスポリシー 主要概念 極限定理 拡張 情報システム