M/M/2の出発過程はポアソン過程

M/M/2の出発過程は何?(2)」の続きです。
前回の結論は下の図

と、下の一連の微分方程式

  • \frac{p({\ge}2}2,t)}{dt}=\frac{2u}{t_e}p(1,t)-\frac{2}{t_e}p({\ge}2,t)・・・・・・(8)
  • \frac{p(1,t)}{dt}=\frac{2u}{t_e}p(0,t)-\frac{1}{t_e}p(1,t)-\frac{2u}{t_e}p(1,t)・・・・・・(9)
  • \frac{p(0,t)}{dt}=-\frac{2u}{t_e}p(0,t)・・・・・・(10)

と初期条件

  • p({\ge}2,0)=\frac{2u^2}{1+u}・・・・・・(5)
  • p(1,0)=\frac{2u(1-u)}{1+u}・・・・・・(6)
  • p(0,0)=\frac{1-u}{1+u}・・・・・・(7)

でした。(8)〜(10)のうちすぐ解けるのは(10)で、(10)からは

  • p(0,t)=A\exp\left(-\frac{2u}{t_e}t\right)

が導かれます。ここにA積分定数です。(7)から

  • A=\frac{1-u}{1+u}

であることが分かり、

  • p(0,t)=\frac{1-u}{1+u}\exp\left(-\frac{2u}{t_e}t\right)・・・・・・(11)

となります。


次に(11)を(9)に代入して

  • \frac{p(1,t)}{dt}=\frac{2u}{t_e}\frac{1-u}{1+u}\exp\left(-\frac{2u}{t_e}t\right)-\frac{1}{t_e}p(1,t)-\frac{2u}{t_e}p(1,t)・・・・・・(12)

これをどうやって解けばよろしいでしょうか? (6)を考慮して

  • p(1,t)=\frac{2u(1-u)}{1+u}\exp(-at)・・・・・・(13)

であると仮定してみましょう。ここでaはまだ値の分からない定数であるとします。(13)を(12)に代入すると

  • -a\frac{2u(1-u)}{1+u}\exp(-at)=\frac{2u}{t_e}\frac{1-u}{1+u}\exp\left(-\frac{2u}{t_e}t\right)-\frac{1}{t_e}\frac{2u(1-u)}{1+u}\exp(-at)-\frac{2u}{t_e}\frac{2u(1-u)}{1+u}\exp(-at)
  • \frac{1}{t_e}\frac{2u(1-u)}{1+u}\exp(-at)+\frac{2u}{t_e}\frac{2u(1-u)}{1+u}\exp(-at)-a\frac{2u(1-u)}{1+u}\exp(-at)=\frac{2u}{t_e}\frac{1-u}{1+u}\exp\left(-\frac{2u}{t_e}t\right)
  • \left[\frac{1}{t_e}+\frac{2u}{t_e}-a\right]\frac{2u(1-u)}{1+u}\exp(-at)=\frac{2u}{t_e}\frac{1-u}{1+u}\exp\left(-\frac{2u}{t_e}t\right)
  • \left[\frac{1}{t_e}+\frac{2u}{t_e}-a\right]\exp(-at)=\frac{1}{t_e}\exp\left(-\frac{2u}{t_e}t\right)・・・・・・(14)

式(14)が成り立つためには

  • \frac{1}{t_e}+\frac{2u}{t_e}-a=\frac{1}{t_e}・・・・・・(15)

  • a=\frac{2u}{t_e}・・・・・・(16)

が成り立つ必要があります。(15)を変形すると

  • \frac{2u}{t_e}-a=0
  • a=\frac{2u}{t_e}

となり、(16)と同じ結果になるので、(16)が成り立てば(15)も成り立つことが分かり、よって(14)も成り立ちます。よって(13)に(16)を代入して

  • p(1,t)=\frac{2u(1-u)}{1+u}\exp\left\(-\frac{2u}{t_e}t\right)・・・・・・(17)

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


次に(17)を(8)に代入します。

  • \frac{p({\ge}2}2,t)}{dt}=\frac{2u}{t_e}\frac{2u(1-u)}{1+u}\exp\left\(-\frac{2u}{t_e}t\right)-\frac{2}{t_e}p({\ge}2,t)・・・・・・(18)

これも、先ほどと同じように

  • p({\ge}2,t)=\frac{2u^2}{1+u}\exp(-bt)・・・・・・(19)

と仮定して解きます。(19)の形を決めるには(5)を考慮しました。(19)を(18)に代入して

  • -b\frac{2u^2}{1+u}\exp(-bt)=\frac{2u}{t_e}\frac{2u(1-u)}{1+u}\exp\left\(-\frac{2u}{t_e}t\right)-\frac{2}{t_e}\frac{2u^2}{1+u}\exp(-bt)
  • \frac{2}{t_e}\frac{2u^2}{1+u}\exp(-bt)-b\frac{2u^2}{1+u}\exp(-bt)=\frac{2u}{t_e}\frac{2u(1-u)}{1+u}\exp\left\(-\frac{2u}{t_e}t\right)
  • \left[\frac{2}{t_e}-b\right]\frac{2u^2}{1+u}\exp(-bt)=\frac{2u}{t_e}\frac{2u(1-u)}{1+u}\exp\left\(-\frac{2u}{t_e}t\right)
  • \left[\frac{2}{t_e}-b\right]\exp(-bt)=\frac{2(1-u)}{t_e}\exp\left\(-\frac{2u}{t_e}t\right)・・・・・・(20)

式(20)が成り立つためには

  • \frac{2}{t_e}-b=\frac{2(1-u)}{t_e}・・・・・・(21)

  • b=\frac{2u}{t_e}・・・・・・(22)

が成り立つ必要があります。(21)を変形すると

  • b=\frac{2}{t_e}-\frac{2(1-u)}{t_e}=\frac{2u}{t_e}

となって(22)に一致するので、(22)が成り立てば(21)も成り立ち、従って(20)も成り立つことが分かります。よって(18)の解は、(19)に(22)を代入した

  • p({\ge}2,t)=\frac{2u^2}{1+u}\exp\left(-\frac{2u}{t_e}t\right)・・・・・・(23)

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


では、ジョブが出発する確率密度はtの経過とともにどのように変化するでしょうか? ジョブが出発することに対応する遷移は状態「2」から状態「出発後:1か2」への遷移と、状態「1」から状態「出発後:0」への遷移の2つがあります。状態「2」から状態「出発後:1か2」への遷移の確率密度は

  • \frac{2}{t_e}p({\ge}2,t)=\frac{4u^2}{(1+u)t_e}\exp\left(-\frac{2u}{t_e}t\right)・・・・・・(24)

状態「1」から状態「出発後:0」への遷移の確率密度は

  • \frac{1}{t_e}p(1,t)=\frac{2u(1-u)}{(1+u)t_e}\exp\left(-\frac{2u}{t_e}t\right)・・・・・・(25)

式(24)と(25)の右辺を足して

  • \frac{4u^2+2u(1-u)}{(1+u)t_e}\exp\left(-\frac{2u}{t_e}t\right)=\frac{4u^2+2u-2u^2}{(1+u)t_e}\exp\left(-\frac{2u}{t_e}t\right)
    • =\frac{2u(u+1)}{(1+u)t_e}\exp\left(-\frac{2u}{t_e}t\right)=\frac{2u}{t_e}\exp\left(-\frac{2u}{t_e}t\right)

結局、ジョブが出発する確率密度は

  • \frac{2u}{t_e}\exp\left(-\frac{2u}{t_e}t\right)

で表されます。これは平均間隔

  • \frac{t_e}{2u}

ポアソン分布の式にほかなりません。よって、M/M/2ステーションからのジョブの出発過程はポアソン過程になります。


M/M/mの出発過程はポアソン過程(1)」に続きます。