GI/G/s待ち行列の到着時刻状態分布と定常状態分布の関係を求める試み(失敗)
以前、自分が書いた「GI/M/s待ち行列の到着時刻状態分布の近似式(1)」を見直していて、実はGI/M/s待ち行列で用いた論理の展開がGI/G/sの場合にも適用出来るのではないか、という気がしてきました。そこで、このことを以下に試してみます。
GI/G/s待ち行列において、時間平均で、システム内にジョブが個ある確率(定常状態分布)をで表し、ジョブ到着時にシステム内にジョブが(その到着するジョブを除いて)個ある確率(到着時刻状態分布)をで表すことにします。そして、この両者の関係を求めてみます。
到着直前にシステム内のジョブが個だった場合、ジョブ到着によって個になります。このの変化が単位時間あたり平均何回起きるか考えてみます。装置の稼働率を、装置の平均処理時間をとすると、単位時間内の平均ジョブ到着数は
- ・・・・(1)
となります。到着直前にシステム内のジョブ数が個である確率はになりますので、単位時間あたりの変化が起きる平均回数は
- ・・・・(2)
になります。次に、単位時間あたりの変化が起きる平均回数を考えます。任意の時刻にシステム内にジョブが個あったとしてそれが個になるというのは、ジョブの処理が終了するということです。の場合は、全ての装置が処理中ですから、このジョブ終了は(処理時間分布が指数分布でなくて一般の分布であったとしても)、単位時間あたり
- ・・・・(3)
回になるような気がします。気がしますが、本当でしょうか? もしこれが本当ならば、任意の時刻にシステム内にジョブが個ある確率はなので、の変化は単位時間あたり平均
- ・・・・(4)
回になります。今は定常状態を考えていますから、式(2)と式(4)は等しいはずです。よって
- の場合
よって
- の場合
- ・・・・(5)
となります。
では、式(3)がGI/G/sで成り立つというのは本当でしょうか?
最初、私は本当のような気がしたのですが、調べてみると反例が出てきました。問題を扱いやすくするためにとします。そしてM/D/1待ち行列を考えてみます。式(3)が成り立つとすれば、式(5)からの場合、つまり任意のについて
- ・・・・(6)
となります。さらにM/D/1では到着過程がM(ポアソン過程)なので[PASTA]を適用することが出来て、
- ・・・・(7)
となります。式(6)と(7)から
- ・・・・(8)
よって
- ・・・・(9)
となります。明らかに
- ・・・・(10)
なので
- ・・・・(11)
となります。しかし、これはM/M/1の定常状態分布の式であって、M/D/1ではこの式は成り立ちません。
ということは上の
の場合は、全ての装置が処理中ですから、このジョブ終了は(処理時間分布が指数分布でなくて一般の分布であったとしても)、単位時間あたり
- ・・・・(3)
回になるような気がします。
というのが私の間違いだったということになります。しかし、どこを直せばよいのか、まだ、私には分かりません。今後の宿題です。