GI/G/mにおける待ち時間の近似式の導出メモ

Factory Physics」ではGI/G/mにおける平均待ち時間の近似式として

  • CT_q\appro\left(\frac{c_a^2+c_e^2}{2}\right)CT_{q(M/M/m)}・・・・・(1)
    • ただし
    • CT_q:キューでの待ち時間
    • c_aジョブの到着間隔の変動係数。(変動係数とは、標準偏差/平均 のこと)
    • c_e:装置処理時間の変動係数
    • CT_{q(M/M/m)}:M/M/mにおけるキューでの待ち時間
    • m:装置台数

を提示しています。しかし、その理由を「Factory Physics」は明らかにしていません。ここでは、私の推測する理由を紹介します。


GI/G/mにおける平均待ち時間の近似式に課せられるいくつかの条件を考えます。すると次のようなものがあります。

式(1)は、これらの条件1)〜3)を満たすものとして導出されたと考えられます。


まず1)を満たすことは、M/M/mでは、

  • c_a=c_e=1

なので、これらを(1)に代入すれば

  • CT_q=CT_{q(M/M/m)}

になるので明らかです。


次に2)を満たすことは、M/G/1では、

  • c_a=1
  • m=1

なので、これらを(1)に代入すれば

  • CT_q=\left(\frac{1+c_e^2}{2}\right)CT_{q(M/M/1)}

となります。ここで「M/M/1における待ち時間の式の導出(1)」から

  • CT_{q(M/M/1)}=\frac{u}{1-u}・・・・・(4)

なので

  • CT_q=\left(\frac{1+c_e^2}{2}\right)\frac{u}{1-u}・・・・・(2)

式(2)と一致します。


最期に3)を満たすことですが、まず「M/M/mにおける待ち時間の式の導出(1)」から

  • CT_{q(M/M/m)}=\frac{m^{m-1}u^m}{m!(1-u)^2}{p_0}{t_e}・・・・・(5)

ただし

  • p_0=\frac{1}{\Bigsum_{i=\0}^{m-\1}\left{\frac{(mu)^i}{i!}\right}+\frac{(mu)^m}{m!(1-u)}}・・・・・(6)

です。(6)から

  • \frac{p_0}{1-u}=\frac{1}{\Bigsum_{i=\0}^{m-\1}\left{\frac{(mu)^i}{i!}(1-u)\right}+\frac{(mu)^m}{m!}}・・・・・(7)

ここでu\rightar{1}とすると分母の最後の項だけが残るので

  • \lim_{u\rightar{1}}\frac{p_0}{1-u}=\frac{1}{\frac{(mu)^m}{m!}}=\frac{m!}{(mu)^m}・・・・・(8)

(1)と(5)から

  • CT_q=\left(\frac{c_a^2+c_e^2}{2}\right)\frac{m^{m-1}u^m}{m!(1-u)^2}{p_0}{t_e}

よって

  • (1-u)CT_q=\left(\frac{c_a^2+c_e^2}{2}\right)\frac{m^{m-1}u^m}{m!(1-u)}{p_0}{t_e}

よって

  • \lim_{u\rightar{1}}(1-u)CT_q=\left(\frac{c_a^2+c_e^2}{2}\right){t_e}\lim_{u\rightar{1}}\frac{m^{m-1}u^m}{m!(1-u)}{p_0}・・・・・(9)

ここで(8)を考慮すると

  • \lim_{u\rightar{1}}(1-u)CT_q=\left(\frac{c_a^2+c_e^2}{2}\right){t_e}\lim_{u\rightar{1}}\frac{m^{m-1}u^m}{m!}\frac{m!}{(mu)^m}
  • \lim_{u\rightar{1}}(1-u)CT_q=\left(\frac{c_a^2+c_e^2}{2}\right){t_e}\lim_{u\rightar{1}}\frac{1}{m}
  • \lim_{u\rightar{1}}(1-u)CT_q=\left(\frac{c_a^2+c_e^2}{2m}\right)t_e

よって(3)を満足します。


以上で、式(1)が決められた理由の推測を終ります。