M/M/sの平均待ち時間とΠの関係式

Πを使ってM/G/mの平均待ち行列長を表す」の式(4)を求める別のやり方を見つけましたので、書いておきます。これは「M/G/1における待ち時間の式の導出(1)」で行った方法をM/M/sに適用したものです。到着過程がポアソン過程なのでPASTAを利用出来ます。ジョブが到着した時にそのジョブが見る待ち行列長はL_qになるはずです。これらのジョブ、つまり到着したジョブより先に待っているジョブのうち、先頭のジョブが処理を開始してから、この到着したジョブが処理開始するまでに要する時間の平均値は、装置の平均処理時間がt_eで装置がs台なので

  • L_q\frac{t_e}{s}・・・・(1)

になります。さらに今処理中のジョブがあります。もし、このジョブ到着時に装置が1台も空いていないとしますと、それからどれかの装置が空くまでの時間の平均値は、処理時間分布が指数分布であることを考慮すれば、

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

となります。これがすなわち「到着したジョブより先に待っているジョブのうち、先頭のジョブが処理を開始」するまでに待つ時間です。しかし、ジョブが到着した時に、いつも装置が1台も空いていないとは限りません。装置が全て処理中である確率は\Pi(M/M/s)で表されるのでした。装置が全て処理中である場合は式(2)だけ待ちますが、1台でも空いている装置があると待ち時間はゼロです。よって、先頭のジョブが待つ平均時間は

  • \Pi(M/M/s)\frac{t_e}{s}・・・・(3)

になります。以上のことから、到着したジョブが処理開始を待つ平均時間CT_q

  • CT_q=L_q\frac{t_e}{s}+\Pi(M/M/s)\frac{t_e}{s}・・・・(4)

となることが言えます。ここで待ち行列L_qをWIPとみなしてリトルの法則を適用してみます。この場合サイクルタイムはCT_qスループット\frac{su}{t_e}になります。よって

  • L_q=CT_q\frac{su}{t_e}・・・・(5)

をいうことが出来ます。式(5)を式(4)に代入すると

  • CT_q=CT_q\frac{su}{t_e}\times\frac{t_e}{s}+\Pi(M/M/s)\frac{t_e}{s}

よって

    • CT_q=CT_qu+\Pi(M/M/s)\frac{t_e}{s}
    • (1-u)CT_q=\Pi(M/M/s)\frac{t_e}{s}
  • CT_q=\frac{\Pi(M/M/s)}{s(1-u)}t_e・・・・(6)

これは「Πを使ってM/G/mの平均待ち行列長を表す」の式(4)にほかなりません。このようにして式(6)を導くことが出来ました。