「《ポアソン過程と出生死滅過程》」の版間の差分

2007年7月9日 (月) 11:21時点における版

【ぽあそんかていとしゅっせいしめつかてい (Poisson process and birth and death process) 】

　ポアソン過程 (Poisson process) は, ランダムに生起する事象を表す基本的な確率過程で, 客の到着や故障の発生, 個体の出生など様々な現象のモデル化に使われる. 一方, 出生死滅過程は個体の出生だけでなくランダムな死滅も考慮した確率過程で, 待ち行列理論をはじめ広く利用されている.

ポアソン過程　事象の生起時点列を構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle 0 \le T_1 \le T_2 \le ...\, } とし, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle N(t)\, } を区間構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle [0, t]\, } における事象の生起数, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle N(u,v) = N(v) - N(u)\, } を区間構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (u, v]\, } での生起数とする. このような確率過程構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{N(t), t\ge 0\}\, } は一般に計数過程と呼ばれる. 計数過程構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{N(t)\}\, } がポアソン過程であるとは, 正の実数構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda\, } が存在して任意の構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle t\ge 0\, } および構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle h>0\, } に対して

構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{array}{lll} &\mathrm{P}(N(t,t+h) = 1 \ | \ T_1,...,T_{N(t)}) = \lambda h + o(h), & \qquad (1)\\ \\ &\mathrm{P}(N(t,t+h) \geq 2 \ | \ T_1,...,T_{N(t)}) = o(h). & \qquad (2) \end{array} }

が成り立つことである.

　(1), (2) はランダムな事象の生起を3つの点で特徴付けている. 第1は, 微小区間構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (t, t+h]\, } に事象が生起する確率は時刻構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle t\, } 以前の挙動に独立であるという点, 第2は, 微小区間に2つ以上の事象が生起する確率は無視できるという点, 第3は, 微小区間に事象の生起する確率が時刻によらない点である. 式 (1) の構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda\, } を強度 (intensity) または生起率と呼ぶ. これは単位時間あたりの平均生起数を表す. 強度を時間の関数構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda(t)\, } に拡張したものは非定常ポアソン過程と呼ばれる. 以下はポアソン過程の性質であり, それぞれがポアソン過程の同値な定義でもある.

性質1　ポアソン過程構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{N(t)\}\, } において,事象の生起間隔の列構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle U_i =T_{i+1} - T_i\, } は互いに独立で平均構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle 1/\lambda\, } の指数分布に従う. \medskip

性質2　ポアソン過程構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{N(t)\}\, } は独立増分過程で, 任意の構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle s<t\, } に対して構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle N(s,t)\, } は平均構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda (t-s)\, } のポアソン分布に従う.

　性質1は指数分布の無記憶性から自然に導かれる. また, 性質2より複数の独立なポアソン過程の重ね合わせは, それぞれの強度の和を強度に持つポアソン過程となることが分かる. また, 次の定理は確率変数の和に対する少数の法則の確率過程版である.

定理1　各構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle k\, } に対して構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \ell_k\, } 個の計数過程構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{N_{k1}(t)\}, \cdots, \{N_{k\ell_k}(t)\}\, } を考え, その重ね合わせを構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle N_k(t) =N_{k1}(t)+ \cdots +N_{k\ell_k}(t)\, } とする. 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \textstyle \lim_{k\to\infty} \ell_k=\infty\, } で, かつ (a) 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{N_{ki}(t)\}, \, i=1, \ldots , \ell_k\, } は互いに独立, (b) 任意の構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle u<v\, } に対して構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \textstyle \lim_{k\to\infty} \sup_{1\le i \le \ell_k} \mathrm{P}(N_{ki}(u,v) \ge 1) = 0\, } が成り立つとすると, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \textstyle k\to\infty\, } のとき構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{N_k(t)\}\, } が平均測度構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{\Lambda(t)\}\, } の (非定常) ポアソン過程に収束するための必要十分条件は, 任意の構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle u<v\, } に対して, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \textstyle \lim_{k\to\infty} \sum_{i=1}^{\ell_k} \mathrm{P}(N_{ki}(u,v)=1) =\Lambda(v) - \Lambda(u)\, } および $\textstyle \lim _{k\to \infty }\sum _{i=1}^{\ell _{k}}\mathrm {P} (N_{ki}(u,v)>1)=0\,$ が成り立つことである. なお, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \Lambda(t)\, } が微分可能ならば強度は構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda(t) = \mbox{d}\Lambda(t)/\mbox{d}t\, } となる.

　定理1は, 実際に起こる様々な現象をポアソン過程を用いて表わすことの妥当性を示唆している. 例えば, 電話網のある回線群への接続要求 (呼) は非常に多くの電話機からかかってくる呼の重ね合わせとみなせる. この場合, 各電話機は独立に使われており (仮定 (a)), その頻度は十分小さい (仮定 (b)) と考えられるため, この回線群への呼の発生はポアソン過程としてモデル化できるであろう. この他にも, マルチンゲールによるポアソン過程の特徴付けや, 事象平均と時間平均の同等性を示すPASTA (Poisson arrivals see time averages) など, ポアソン過程には興味深い性質が多い.

ポアソン過程の一般化　ポアソン過程を特徴付ける3つの条件のうち第2の条件を緩め, 事象の生起時点列はポアソン過程であるが, 各生起時点で同時に発生する事象の数は独立で同一の分布に従う確率変数である場合, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle N(t)\, } は複合ポアソン過程と呼ばれる. また, 非定常ポアソン過程の強度 $\lambda (t)\,$ を確率過程に拡張したものは2重確率ポアソン過程 (doubly stochastic Poisson process) と呼ばれる. 例えば, マルコフ変調ポアソン過程は構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda(t)\, } が連続時間マルコフ連鎖に従う例である.

出生過程　性質1より, ポアソン過程は状態空間構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{0, 1, ...\}\, } 上の連続時間マルコフ連鎖であることがわかる. 推移速度行列を構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \boldsymbol{Q} =(q_{ij})\, } とすると, 性質1から構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle q_{i,i+1} = -q_{ii} = \lambda, \, i\ge 0\, } でその他の構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \boldsymbol{Q}\, } の要素は全て0となる. これを一般化して, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle i\, } から構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle i+1\, } への推移速度が構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle i\, } に依存して $\lambda _{i}\,$ で定まるマルコフ連鎖を出生過程 (birth process)と呼ぶ. 出生過程の推移速度行列は構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle q_{i,i+1} = -q_{ii} = \lambda_i, \, i\ge 0\, } で, その他の要素は0である.

出生死滅過程　出生過程では, 状態は構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle i\, } から構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle i+1\, } というように1ずつ進んでいくが, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle i\, } から構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle i-1\, } へ戻ることも許すように一般化すると, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle q_{i,i+1} = \lambda_i, \, q_{i+1,i} = \mu_{i+1}, \, i\ge 0\, } かつ構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle q_{00} =-\lambda_0, \, q_{ii} = -(\lambda_i + \mu_i), \, i\ge 1\, } で, その他の要素は0の推移速度行列が得られる. このような3重対角の推移速度行列に従う連続時間マルコフ連鎖を出生死滅過程 (birth and death process) という. また, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda_i\, } , 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \mu_i\, } はそれぞれ状態構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle i\, } での出生率, 死滅率と呼ばれる. 出生死滅過程では, 状態構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle i\; (\ge 1)\, } に滞在する時間の長さはパラメータ $\lambda _{i}+\mu _{i}\,$ の指数分布に従い, 滞在時間を終えると確率構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda_i/(\lambda_i+\mu_i)\, } で状態構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle i+1\, } へ, 確率 $\mu _{i}/(\lambda _{i}+\mu _{i})\,$ で状態構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle i-1\, } へ推移する.

　出生死滅過程は隣り合う状態間でのみ推移が起きるという特徴を持つため, 定常分布などの特性量が陽な形で得られる. 例えば, 応用上重要な構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda_i=\lambda\, } , 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \mu_i=\mu\, } の出生死滅過程は, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda < \mu\, } のとき正再帰的で, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \rho=\lambda/\mu\, } とすると状態構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle j\, } にいる定常確率は $\pi _{j}=(1-\rho )\rho ^{j},\;j=0,1,\ldots \,$ という幾何分布となる. なお, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda = \mu\, } のときは零再帰的, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda > \mu\, } のときは一時的となり定常分布は存在しない. この例はM/M/1 待ち行列モデルに相当する出生死滅過程であるが, 出生死滅過程はより一般的なM/M/c 待ち行列モデル (M/M/構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle c\, } 待ち行列モデル) などのマルコフ型の待ち行列モデルや, 機械修理モデルを解析する上でも重要な確率過程となっている.

参考文献

[1] P. Brémaud, Point Processes and Queues, Springer-Verlag, 1981.

[2] D. R. Cox and V. Isham, Point Processes, Chapman and Hall, 1980.

[3] R. W. Wolff, Stochastic Modeling and the Theory of Queues, Prentice-Hall, 1989.

[4] 宮沢政清, 『確率と確率過程』, 近代科学社, 1993.

@@ 3行目: / 3行目: @@
 　[[ポアソン過程]] (Poisson process) は, ランダムに生起する事象を表す基本的な[[確率過程]]で, 客の到着や故障の発生, 個体の出生など様々な現象のモデル化に使われる. 一方, [[出生死滅過程]]は個体の出生だけでなくランダムな死滅も考慮した確率過程で, [[待ち行列理論]]をはじめ広く利用されている.
-'''ポアソン過程'''　事象の生起時点列を $<math>0 \le T_1 \le T_2 \le ...\, </math>$ とし, $<math>N(t)\, </math>$ を区間 $<math>[0, t]\, </math>$ における事象の生起数, $<math>N(u,v) = N(v) - N(u)\, </math>$ を区間 $<math>(u, v]\, </math>$ での生起数とする. このような確率過程$<math>\{N(t), t\ge 0\}\, </math>$ は一般に計数過程と呼ばれる. 計数過程 $<math>\{N(t)\}\, </math>$ がポアソン過程であるとは, 正の実数 $<math>\lambda\, </math>$ が存在して任意の $<math>t\ge 0\, </math>$ および $<math>h>0\, </math>$ に対して
+'''ポアソン過程'''　事象の生起時点列を <math>0 \le T_1 \le T_2 \le ...\, </math> とし, <math>N(t)\, </math> を区間 <math>[0, t]\, </math> における事象の生起数, <math>N(u,v) = N(v) - N(u)\, </math> を区間 <math>(u, v]\, </math> での生起数とする. このような確率過程<math>\{N(t), t\ge 0\}\, </math> は一般に計数過程と呼ばれる. 計数過程 <math>\{N(t)\}\, </math> がポアソン過程であるとは, 正の実数 <math>\lambda\, </math> が存在して任意の <math>t\ge 0\, </math> および <math>h>0\, </math> に対して
-\begin{eqnarray}
+:<math>
-\mathrm{P}(N(t,t+h) = 1 \, | \, T_1,...,T_{N(t)}) &=& \lambda h + o(h),
+\begin{array}{lll}
-\label{B-D-04+def11} \\
+&\mathrm{P}(N(t,t+h) = 1 \ | \ T_1,...,T_{N(t)}) = \lambda h + o(h), & \qquad (1)\\
-\mathrm{P}(N(t,t+h) \geq 2 \, | \, T_1,...,T_{N(t)}) &=& o(h). \label{B-D-04+def12}
+\\
-\end{eqnarray}
+&\mathrm{P}(N(t,t+h) \geq 2 \ | \ T_1,...,T_{N(t)}) = o(h). & \qquad (2)
+\end{array}
+</math>
 が成り立つことである.
-　(1), (2) はランダムな事象の生起を3つの点で特徴付けている. 第1は, 微小区間 $<math>(t, t+h]\, </math>$ に事象が生起する確率は時刻 $<math>t\, </math>$ 以前の挙動に独立であるという点, 第2は, 微小区間に2つ以上の事象が生起する確率は無視できるという点, 第3は, 微小区間に事象の生起する確率が時刻によらない点である. 式 (1) の $<math>\lambda\, </math>$ を強度 (intensity)  または生起率と呼ぶ. これは単位時間あたりの平均生起数を表す. 強度を時間の関数 $\<math>lambda(t)\, </math>$ に拡張したものは[[非定常ポアソン過程]]と呼ばれる. 以下はポアソン過程の性質であり, それぞれがポアソン過程の同値な定義でもある.
+　(1), (2) はランダムな事象の生起を3つの点で特徴付けている. 第1は, 微小区間 <math>(t, t+h]\, </math> に事象が生起する確率は時刻 <math>t\, </math> 以前の挙動に独立であるという点, 第2は, 微小区間に2つ以上の事象が生起する確率は無視できるという点, 第3は, 微小区間に事象の生起する確率が時刻によらない点である. 式 (1) の <math>\lambda\, </math> を強度 (intensity)  または生起率と呼ぶ. これは単位時間あたりの平均生起数を表す. 強度を時間の関数 <math>\lambda(t)\, </math> に拡張したものは[[非定常ポアソン過程]]と呼ばれる. 以下はポアソン過程の性質であり, それぞれがポアソン過程の同値な定義でもある.
-'''性質1'''　ポアソン過程 $<math>\{N(t)\}\, </math>$ において,事象の生起間隔の列 $<math>U_i =T_{i+1} - T_i\, </math>$ は互いに独立で平均 $<math>1/\lambda\, </math>$ の
+'''性質1'''　ポアソン過程 <math>\{N(t)\}\, </math> において,事象の生起間隔の列 <math>U_i =T_{i+1} - T_i\, </math> は互いに独立で平均 <math>1/\lambda\, </math> の
 [[指数分布]]に従う.
 \medskip
-'''性質2'''　ポアソン過程 $<math>\{N(t)\}\, </math>$ は[[独立増分過程]]で, 任意の $<math>s<t\, </math>$ に対して $<math>N(s,t)\, </math>$ は平均 $<math>\lambda (t-s)\, </math>$ の[[ポアソン分布]]に従う.
+'''性質2'''　ポアソン過程 <math>\{N(t)\}\, </math> は[[独立増分過程]]で, 任意の <math>s<t\, </math> に対して <math>N(s,t)\, </math> は平均 <math>\lambda (t-s)\, </math> の[[ポアソン分布]]に従う.
 　性質1は[[指数分布の無記憶性 (指数分布の)|指数分布の無記憶性]]から自然に導かれる. また, 性質2より複数の独立なポアソン過程の重ね合わせは, それぞれの強度の和を強度に持つポアソン過程となることが分かる. また, 次の定理は確率変数の和に対する[[少数の法則]]の確率過程版である.
-'''定理1'''　各 $<math>k\, </math>$ に対して $<math>\ell_k\, </math>$ 個の計数過程 $<math>\{N_{k1}(t)\}, \cdots, \{N_{k\ell_k}(t)\}\, </math>$ を考え, その重ね合わせを $<math>N_k(t) =N_{k1}(t)+ \cdots +N_{k\ell_k}(t)\, </math>$ とする. $<math>\lim_{k\to\infty} \ell_k=\infty\, </math>$ で, かつ (a) $<math>\{N_{ki}(t)\}, \, i=1, \ldots , \ell_k\, </math>$ は互いに独立, (b) 任意の $<math>u<v\, </math>$ に対して $<math>\lim_{k\to\infty} \sup_{1\le i \le \ell_k} \mathrm{P}(N_{ki}(u,v) \ge 1) = 0\, </math>$ が成り立つとすると, $<math>k\to\infty\, </math>$ のとき $<math>\{N_k(t)\}\, </math>$ が[[平均測度]] $<math>\{\Lambda(t)\}\, </math>$ の (非定常) ポアソン過程に収束するための必要十分条件は, 任意の $<math>u<v\, </math>$ に対して, $<math>\lim_{k\to\infty} \sum_{i=1}^{\ell_k} \mathrm{P}(N_{ki}(u,v)=1) =\Lambda(v) - \Lambda(u)\, </math>$ および $<math>\lim_{k\to\infty} \sum_{i=1}^{\ell_k} \mathrm{P}(N_{ki}(u,v)>1) = 0\, </math>$が成り立つことである. なお, $<math>\Lambda(t)\, </math>$ が微分可能ならば強度は $<math>\lambda(t) = \mbox{d}\Lambda(t)/\mbox{d}t\, </math>$ となる.
+'''定理1'''　各 <math>k\, </math> に対して <math>\ell_k\, </math> 個の計数過程 <math>\{N_{k1}(t)\}, \cdots, \{N_{k\ell_k}(t)\}\, </math> を考え, その重ね合わせを <math>N_k(t) =N_{k1}(t)+ \cdots +N_{k\ell_k}(t)\, </math> とする. <math>\textstyle \lim_{k\to\infty} \ell_k=\infty\, </math> で, かつ (a) <math>\{N_{ki}(t)\}, \, i=1, \ldots , \ell_k\, </math> は互いに独立, (b) 任意の <math>u<v\, </math> に対して <math>\textstyle \lim_{k\to\infty} \sup_{1\le i \le \ell_k} \mathrm{P}(N_{ki}(u,v) \ge 1) = 0\, </math> が成り立つとすると, <math>\textstyle k\to\infty\, </math> のとき <math>\{N_k(t)\}\, </math> が[[平均測度]] <math>\{\Lambda(t)\}\, </math> の (非定常) ポアソン過程に収束するための必要十分条件は, 任意の <math>u<v\, </math> に対して, <math>\textstyle \lim_{k\to\infty} \sum_{i=1}^{\ell_k} \mathrm{P}(N_{ki}(u,v)=1) =\Lambda(v) - \Lambda(u)\, </math> および <math>\textstyle \lim_{k\to\infty} \sum_{i=1}^{\ell_k} \mathrm{P}(N_{ki}(u,v)>1) = 0\, </math>が成り立つことである. なお, <math>\Lambda(t)\, </math> が微分可能ならば強度は <math>\lambda(t) = \mbox{d}\Lambda(t)/\mbox{d}t\, </math> となる.
 　定理1は, 実際に起こる様々な現象をポアソン過程を用いて表わすことの妥当性を示唆している. 例えば, 電話網のある回線群への接続要求 (呼) は非常に多くの電話機からかかってくる呼の重ね合わせとみなせる. この場合, 各電話機は独立に使われており (仮定 (a)), その頻度は十分小さい (仮定 (b)) と考えられるため, この回線群への呼の発生はポアソン過程としてモデル化できるであろう. この他にも, [[マルチンゲール]]によるポアソン過程の特徴付けや, 事象平均と時間平均の同等性を示す[[PASTA]] (Poisson arrivals see time averages) など, ポアソン過程には興味深い性質が多い.
-'''ポアソン過程の一般化'''　ポアソン過程を特徴付ける3つの条件のうち第2の条件を緩め, 事象の生起時点列はポアソン過程であるが, 各生起時点で同時に発生する事象の数は独立で同一の分布に従う確率変数である場合, $<math>N(t)\, </math>$ は複合ポアソン過程と呼ばれる. また, 非定常ポアソン過程の強度 $<math>\lambda(t)\, </math>$ を確率過程に拡張したものは2重確率ポアソン過程 (doubly stochastic Poisson process) と呼ばれる. 例えば, [[マルコフ変調ポアソン過程]]は $<math>\lambda(t)\, </math>$ が連続時間マルコフ連鎖に従う例である.
+'''ポアソン過程の一般化'''　ポアソン過程を特徴付ける3つの条件のうち第2の条件を緩め, 事象の生起時点列はポアソン過程であるが, 各生起時点で同時に発生する事象の数は独立で同一の分布に従う確率変数である場合, <math>N(t)\, </math> は複合ポアソン過程と呼ばれる. また, 非定常ポアソン過程の強度 <math>\lambda(t)\, </math> を確率過程に拡張したものは2重確率ポアソン過程 (doubly stochastic Poisson process) と呼ばれる. 例えば, [[マルコフ変調ポアソン過程]]は <math>\lambda(t)\, </math> が連続時間マルコフ連鎖に従う例である.
-'''出生過程'''　性質1より, ポアソン過程は[[状態空間]] $<math>\{0, 1, ...\}\, </math>$ 上の[[連続時間マルコフ連鎖]]であることがわかる. [[推移速度行列]]を $<math>\mbox{\boldmath$Q$} =(q_{ij})\, </math>$ とすると, 性質1から $<math>q_{i,i+1} = -q_{ii} = \lambda, \, i\ge 0\, </math>$ でその他の $<math>\mbox{\boldmath$Q$}\, </math>$ の要素は全て0となる. これを一般化して, $<math>i\, </math>$ から $<math>i+1\, </math>$ への推移速度が $<math>i\, </math>$ に依存して $<math>\lambda_i\, </math>$ で定まるマルコフ連鎖を[[出生過程]] (birth process)と呼ぶ. 出生過程の推移速度行列は$<math>q_{i,i+1} = -q_{ii} = \lambda_i, \, i\ge 0\, </math>$ で, その他の要素は0である.
+'''出生過程'''　性質1より, ポアソン過程は[[状態空間]] <math>\{0, 1, ...\}\, </math> 上の[[連続時間マルコフ連鎖]]であることがわかる. [[推移速度行列]]を <math>\boldsymbol{Q} =(q_{ij})\, </math> とすると, 性質1から <math>q_{i,i+1} = -q_{ii} = \lambda, \, i\ge 0\, </math> でその他の <math>\boldsymbol{Q}\, </math> の要素は全て0となる. これを一般化して, <math>i\, </math> から <math>i+1\, </math> への推移速度が <math>i\, </math> に依存して <math>\lambda_i\, </math> で定まるマルコフ連鎖を[[出生過程]] (birth process)と呼ぶ. 出生過程の推移速度行列は<math>q_{i,i+1} = -q_{ii} = \lambda_i, \, i\ge 0\, </math> で, その他の要素は0である.
-'''出生死滅過程'''　出生過程では, 状態は $<math>i\, </math>$ から <math>$i+1\, </math>$ というように1ずつ進んでいくが, $<math>i\, </math>$ から $<math>i-1\, </math>$ へ戻ることも許すように一般化すると, $<math>q_{i,i+1} = \lambda_i, \, q_{i+1,i} = \mu_{i+1}, \, i\ge 0\, </math>$ かつ $<math>q_{00} =-\lambda_0, \, q_{ii} = -(\lambda_i + \mu_i), \, i\ge 1\, </math>$ で, その他の要素は0の推移速度行列が得られる. このような3重対角の推移速度行列に従う連続時間マルコフ連鎖を[[出生死滅過程]] (birth and death process) という. また, $<math>\lambda_i\, </math>$, $<math>\mu_i\, </math>$ はそれぞれ状態 $<math>i\, </math>$ での出生率, 死滅率と呼ばれる. 出生死滅過程では, 状態 $<math>i\; (\ge 1)\, </math>$ に滞在する時間の長さはパラメータ $<math>\lambda_i+\mu_i\, </math>$ の指数分布に従い, 滞在時間を終えると確率 $<math>\lambda_i/(\lambda_i+\mu_i)\, </math>$ で状態 $<math>i+1\, </math>$ へ, 確率 $<math>\mu_i/(\lambda_i+\mu_i)\, </math>$ で状態 $<math>i-1\, </math>$ へ推移する.
+'''出生死滅過程'''　出生過程では, 状態は <math>i\, </math> から <math>i+1\, </math> というように1ずつ進んでいくが, <math>i\, </math> から <math>i-1\, </math> へ戻ることも許すように一般化すると, <math>q_{i,i+1} = \lambda_i, \, q_{i+1,i} = \mu_{i+1}, \, i\ge 0\, </math> かつ <math>q_{00} =-\lambda_0, \, q_{ii} = -(\lambda_i + \mu_i), \, i\ge 1\, </math> で, その他の要素は0の推移速度行列が得られる. このような3重対角の推移速度行列に従う連続時間マルコフ連鎖を[[出生死滅過程]] (birth and death process) という. また, <math>\lambda_i\, </math>, <math>\mu_i\, </math> はそれぞれ状態 <math>i\, </math> での出生率, 死滅率と呼ばれる. 出生死滅過程では, 状態 <math>i\; (\ge 1)\, </math> に滞在する時間の長さはパラメータ <math>\lambda_i+\mu_i\, </math> の指数分布に従い, 滞在時間を終えると確率 <math>\lambda_i/(\lambda_i+\mu_i)\, </math> で状態 <math>i+1\, </math> へ, 確率 <math>\mu_i/(\lambda_i+\mu_i)\, </math> で状態 <math>i-1\, </math> へ推移する.
-　出生死滅過程は隣り合う状態間でのみ[[推移]]が起きるという特徴を持つため, [[定常分布]]などの特性量が陽な形で得られる. 例えば, 応用上重要な $<math>\lambda_i=\lambda\, </math>$, $<math>\mu_i=\mu\, </math>$ の出生死滅過程は, $<math>\lambda < \mu\, </math>$ のとき[[正再帰的]]で, $<math>\rho=\lambda/\mu\, </math>$ とすると状態 $<math>j$\, </math> にいる定常確率は $<math>\pi_j = (1 - \rho)\rho^j, \; j=0,1,\ldots\, </math>$ という[[幾何分布]]となる. なお, $<math>\lambda = \mu\, </math>$ のときは零再帰的, $<math>\lambda > \mu$\, </math> のときは一時的となり定常分布は存在しない. この例は[[M/M/1 待ち行列モデル]]に相当する出生死滅過程であるが, 出生死滅過程はより一般的な[[M/M/c 待ち行列モデル]] (M/M/<math>$c$\, </math> 待ち行列モデル) などのマルコフ型の待ち行列モデルや, [[機械修理モデル]]を解析する上でも重要な確率過程となっている.
+　出生死滅過程は隣り合う状態間でのみ[[推移]]が起きるという特徴を持つため, [[定常分布]]などの特性量が陽な形で得られる. 例えば, 応用上重要な <math>\lambda_i=\lambda\, </math>, <math>\mu_i=\mu\, </math> の出生死滅過程は, <math>\lambda < \mu\, </math> のとき[[正再帰的]]で, <math>\rho=\lambda/\mu\, </math> とすると状態 <math>j\, </math> にいる定常確率は <math>\pi_j = (1 - \rho)\rho^j, \; j=0,1,\ldots\, </math> という[[幾何分布]]となる. なお, <math>\lambda = \mu\, </math> のときは零再帰的, <math>\lambda > \mu\, </math> のときは一時的となり定常分布は存在しない. この例は[[M/M/1 待ち行列モデル]]に相当する出生死滅過程であるが, 出生死滅過程はより一般的な[[M/M/c 待ち行列モデル]] (M/M/<math>c\, </math> 待ち行列モデル) などのマルコフ型の待ち行列モデルや, [[機械修理モデル]]を解析する上でも重要な確率過程となっている.

「《ポアソン過程と出生死滅過程》」の版間の差分

2007年7月9日 (月) 11:21時点における版

案内メニュー

検索