Ｄ／Ｍ／１待ち行列のジョブのサイクルタイムの変動係数

待ち行列理論

今度はＤ／Ｍ／１待ち行列におけるジョブのサイクルタイムの変動係数を求めてみます。


「Ｄ／Ｍ／１待ち行列の到着時刻状態分布（２）」の式（１４）（ここでは番号を振り直して式（１）とします）が示すようにジョブが待ち行列システムに到着した時、このシステムにジョブが（処理中のものと待っているものの両方を含めて）k個ある確率は、

になります。そうすると以下「Ｍ／Ｍ／１待ち行列のジョブのサイクルタイムの変動係数」での検討内容をそのまま適用することが出来ます。まず、装置の処理時間の平均をt_eとします。次に、ジョブが到着した時、システム内にジョブがk個ある場合の、この到着したジョブのサイクルタイムは確率的に変動するので確率変数CT_kで表します。そうするとその平均は

  • E(CT_k)=(k+1)t_e・・・・・（３）

で表され、その分散は

  • V(CT_k)=(k+1)t_e^2・・・・・（４）

で表されます。
さて、ジョブが到着した時、システム内にジョブがk個ある確率は式（１）で与えられるので、全ての場合に渡って平均をとった、ジョブのサイクルタイムの平均E(CT)は、式（１）と（３）から

  • E(CT)=\Bigsum_{k=0}^{\infty}p(k)E(CT_k)=\Bigsum_{k=0}^{\infty}(1-b)b^k(k+1)t_e=(1-b)t_e\Bigsum_{k=0}^{\infty}(k+1)b^k
    • =(1-b)t_e\Bigsum_{k=1}^{\infty}kb^{k-1}

ここで「補足」の式（２）を用いれば、

  • (1-b)t_e\Bigsum_{k=1}^{\infty}kb^{k-1}=(1-b)t_e\frac{1}{(1-b)^2}=\frac{t_e}{1-b}

よって

  • E(CT)=\frac{t_e}{1-b}・・・・・（５）


同様にして全ての場合に渡ってのサイクルタイムの分散を求めます。そのためにまず、全ての場合に渡っての２乗平均を求めます。ジョブが到着した時、システム内にジョブがk個ある場合のCT_kの２乗平均E(CT_k^2は

  • E(CT_k^2)=V(CT_k)+E(CT_k)^2=(k+1)t_e^2+(k+1)^2t_e^2=(k+2)(k+1)t_e^2・・・・・（６）

で求めることが出来ます。この値を全ての場合に渡って平均します。

  • E(CT^2)=\Bigsum_{k=0}^{\infty}p(k)E(CT_k^2)=\Bigsum_{k=0}^{\infty}(1-b)b^k(k+2)(k+1)t_e^2=(1-b)t_e^2\Bigsum_{k=0}^{\infty}(k+2)(k+1)b^k
    • =(1-b)t_e^2\Bigsum_{k=2}^{\infty}k(k-1)b^k-2

ここで「補足」の式（３）を用いれば、

  • (1-b)t_e^2\Bigsum_{k=2}^{\infty}k(k-1)b^k-2=(1-b)t_e^2\frac{2}{(1-b)^3}=\frac{2t_e^2}{(1-b)^2}

よって

  • E(CT^2)=\frac{2t_e^2}{(1-b)^2}・・・・・（７）

サイクルタイムの分散V(CT)は式（５）と（７）から求めることが出来ます。

  • V(CT)=E(CT^2)-E(CT)^2=\frac{2t_e^2}{(1-b)^2}-\frac{t_e^2}{(1-b)^2}=\frac{t_e^2}{(1-b)^2}・・・・・（８）

よってCTの標準偏差STD(CT)は

  • STD(CT)=\frac{t_e}{1-b}・・・・・（９）

式（５）と（９）から変動係数cは、

  • c=\frac{STD(CT)}{E(CT)}=\frac{t_e}{1-b}\frac{1-b}{t_e}=1

つまり

  • c=1・・・・・（１０）

となります。つまり、Ｄ／Ｍ／１でもＭ／Ｍ／１でもサイクルタイムの変動係数は同じ値になるという結果になりました。おそらくＧＩ／Ｍ／１ならばサイクルタイムの変動係数は１になるのでしょう。