Ｍ／Ｄ／１の定常状態分布の求め方（１） - 工場統計力学（建設中！）

Ｍ／Ｄ／１待ち行列が定常状態の時のシステム内のジョブ数の分布の求め方が分かりました。これは計算方法が分かったというだけであって、きれいな式としてお見せすることは出来ません。
その計算方法をご紹介します。

Ｍ／Ｄ／１待ち行列システム内にジョブが $k$ 個存在する確率は、「Ｍ／Ｇ／ｍの定常状態のジョブ数分布について」で示したように時間平均でみても処理終了時点平均で見ても等しいことになります。そこで、ここでは処理終了直後のジョブ数を考察することにします。定常状態になっていれば、ある処理終了直後とその次の処理終了直後でシステム内にジョブが $k$ 個存在する確率は等しいはずです。このことを用いて定常状態確率を求めていきます。システム内にジョブが $k$ 個存在する定常状態確率を $p(k)$ で表すことにします。

装置の利用率を $u$ で表すと、明らかに

$p(0)=1-u$ ・・・・・（１）

であることが分かります。
次に、１回の装置処理時間 $t_e$ の間にジョブがいくつ到着するかを検討しまする。到着はポアソン過程（Ｍ）なので、ポアソン分布の式（「ポアソン分布」参照）

$p(k,t)=\frac{({\lambda}t)^k}{k!}\exp(-{\lambda}t)$ ・・・・・（２）

の $t$ に $t_e$ を代入すれば、処理時間 $t_e$ の間にジョブが $k$ 個到着する確率 $A(k)$ は、

$A(k)=\frac{({\lambda}t_e)^k}{k!}\exp(-{\lambda}t_e)$ ・・・・・（３）

となることが分かります。ここで

${\lambda}t_e=u$

なので式（３）は

$A(k)=\frac{u^k}{k!}\exp(-u)$ ・・・・・（４）

と書き表すことが出来ます。

２個以上一度に処理終了することが出来ないから、ある処理終了直後のジョブ数が１つ前の処理終了直後のジョブ数より２個以上少ないことはあり得ません。ですので、ある処理終了直後にシステムにジョブが０個であったならば、その１つ前の処理終了直後にはジョブが１個だった場合と、０個だった場合が考えられます。
１つ前の処理終了直後１個ある状態が、処理中にジョブが０個到着して次の処理終了直後に０個になる確率は、