M/M/1における待ち時間の式の導出(3)
ここでは、M/G/1待ち行列での考察をなぞる方法を紹介します。
ロットが到着した時に装置が他のロットを処理中かどうかということを考えます。到着間隔が指数分布、つまりデタラメに到着する、ということから、ロットが到着した時に装置が処理中である確率はそもそも利用率になることが分かります。
- このあたりの議論は厳密にはPASTAという概念の議論になります。PASTAとはPoisson Arrival See Time Averageの略でポアソン到着(到着間隔が指数分布の到着)は(いろいろな現象の)時間平均を見るだろう、ということです。PASTAについてはここでは詳しく述べません。
さらに、到着した時に待ち行列で待っているロットの個数の平均値はこれもPASTAにより待ち行列の(時間)平均の長さになります。また、装置が処理中の時の処理の残り時間はもちろん確率変数ですが、処理時間が指数分布(M)の確率変数であることから残り時間の平均値は処理時間の平均値に等しいことが分かります。よって、このロットが待ち行列で待つ時間の平均値は、今、装置が処理中であるならば、
- (装置が処理中である場合の待ち行列の平均の長さ)×+
であり、今、装置が空いているならば
- ゼロ
です。装置が処理中である確率はですから、ロットが待ち行列で待つ時間の平均は
- ×(装置が処理中である場合の待ち行列の平均の長さ)×+
となります。
- ×(装置が処理中である場合の待ち行列の平均の長さ)
は
- (装置が処理中である確率)×(装置が処理中である場合の待ち行列の平均の長さ)+(装置が処理中でない確率)×0
ここで装置が処理中でない場合は、待ち行列の長さはゼロであることを考慮すれば、上の式はさらに
と書けることになり、よってこれは待ち行列の平均の長さにほかならないことになります。よって、ロットが待ち行列で待つ時間の平均は
- ・・・・・・(1)
と書けます。ここでに対してリトルの法則を適用します。ロットの到着するスループットをとすると
になります。スループットはロットの平均到着間隔の逆数であり、ロットの平均到着間隔は
ですので
よって
これを式(1)に代入すると
よって
- ・・・・・・(2)
これでM/M/1待ち行列の待ち時間の式を導出することが出来ました。