GI/G/s待ち行列の到着時刻状態分布と定常状態分布の関係を求める試み(失敗)

以前、自分が書いた「GI/M/s待ち行列の到着時刻状態分布の近似式(1)」を見直していて、実はGI/M/s待ち行列で用いた論理の展開がGI/G/sの場合にも適用出来るのではないか、という気がしてきました。そこで、このことを以下に試してみます。


GI/G/s待ち行列において、時間平均で、システム内にジョブがk個ある確率(定常状態分布)をp(k)で表し、ジョブ到着時にシステム内にジョブが(その到着するジョブを除いて)k個ある確率(到着時刻状態分布)を\pi(k)で表すことにします。そして、この両者の関係を求めてみます。


到着直前にシステム内のジョブがk個だった場合、ジョブ到着によってk+1個になります。このk{\rightar}k+1の変化が単位時間あたり平均何回起きるか考えてみます。装置の稼働率u、装置の平均処理時間をt_eとすると、単位時間内の平均ジョブ到着数は

  • \frac{su}{t_e}・・・・(1)

となります。到着直前にシステム内のジョブ数がk個である確率は\pi(k)になりますので、単位時間あたりk{\rightar}k+1の変化が起きる平均回数は

  • \frac{su}{t_e}\pi(k)・・・・(2)

になります。次に、単位時間あたりk+1{\rightar}kの変化が起きる平均回数を考えます。任意の時刻にシステム内にジョブがk+1個あったとしてそれがk個になるというのは、ジョブの処理が終了するということです。k+1{\ge}sの場合は、全ての装置が処理中ですから、このジョブ終了は(処理時間分布が指数分布でなくて一般の分布であったとしても)、単位時間あたり

  • \frac{s}{t_e}・・・・(3)

回になるような気がします。気がしますが、本当でしょうか? もしこれが本当ならば、任意の時刻にシステム内にジョブがk+1個ある確率はp(k+1)なので、k+1{\rightar}kの変化は単位時間あたり平均

  • \frac{s}{t_e}p(k+1)・・・・(4)

回になります。今は定常状態を考えていますから、式(2)と式(4)は等しいはずです。よって

  • k+1{\ge}sの場合
    • \frac{su}{t_e}\pi(k)=\frac{s}{t_e}p(k+1)

よって

  • k+1{\ge}sの場合
    • p(k+1)=u\pi(k)・・・・(5)

となります。


では、式(3)がGI/G/sで成り立つというのは本当でしょうか?
最初、私は本当のような気がしたのですが、調べてみると反例が出てきました。問題を扱いやすくするためにs=1とします。そしてM/D/1待ち行列を考えてみます。式(3)が成り立つとすれば、式(5)からk{\ge}0の場合、つまり任意のkについて

  • p(k+1)=u\pi(k)・・・・(6)

となります。さらにM/D/1では到着過程がM(ポアソン過程)なので[PASTA]を適用することが出来て、

  • \pi(k)=p(k)・・・・(7)

となります。式(6)と(7)から

  • p(k+1)=up(k)・・・・(8)

よって

  • p(k)=u^kp(0)・・・・(9)

となります。明らかに

  • p(0)=1-u・・・・(10)

なので

  • p(k)=u^k(1-u)・・・・(11)

となります。しかし、これはM/M/1の定常状態分布の式であって、M/D/1ではこの式は成り立ちません。


ということは上の

k+1{\ge}sの場合は、全ての装置が処理中ですから、このジョブ終了は(処理時間分布が指数分布でなくて一般の分布であったとしても)、単位時間あたり

  • \frac{s}{t_e}・・・・(3)

回になるような気がします。

というのが私の間違いだったということになります。しかし、どこを直せばよいのか、まだ、私には分かりません。今後の宿題です。