GI/G/s待ち行列の平均待ち時間の近似式(3)

さて、今度はD/M/sの平均待ち時間CT_q(D/M/s)を表す数式を(それが近似式であっても)求めることが目標になります。これについて私はいろいろ試みたのですが、うまくいきませんでした。そこで代わりにE_2/M/sの平均待ち時間CT_q(E_2/M/sならば状態遷移図を描くことで求めることが出来るのではないか、と思い、CT_q(E_2/M/s)を求めることにしました。CT_q(E_2/M/s)からCT_q(D/M/s)を求めるには次のようにします。
また、Pageの近似式

  • CT_q(GI/G/s){\approx}c_a^2c_e^2CT_q(M/M/s)+c_a^2(1-c_e^2)CT_q(M/D/s)+c_e^2(1-c_a^2)CT_q(D/M/s)・・・・(2)

を用います。左辺をE_2/M/sに適用すると、c_a^2=1/2c_e=1なので

  • CT_q(E_2/M/s){\approx}\frac{1}{2}CT_q(M/M/s)+\frac{1}{2}CT_q(D/M/s)・・・・(12)

よって

  • 2CT_q(E_2/M/s){\approx}CT_q(M/M/s)+CT_q(D/M/s)
  • CT_q(D/M/s){\approx}2CT_q(E_2/M/s)-CT_q(M/M/s)・・・・(13)

式(13)によってCT_q(E_2/M/s)からCT_q(D/M/s)を(近似的に)求めることが出来ます。


次にCT_q(E_2/M/s)を求めることを考えます。まずs=2の場合、つまりCT_q(E_2/M/2)を考えます。E_2/M/2待ち行列の状態を、到着ジョブのフェーズと待ち行列内のジョブ数の組で表すことにします。ジョブの到着は、フェーズ0とフェーズ1の2段階を経て発生すると考えます。ジョブ到着間隔の平均は、装置が2台あることに注意すれば

  • \frac{t_e}{2u}・・・・(14)

です。そこでフェーズ0もフェーズ1も、その状態に滞在する時間が平均

  • \frac{t_e}{4u}・・・・(15)

の指数分布になっているとすれば、ジョブの到着間隔は平均

  • \frac{t_e}{2u}・・・・(14)

の2次のアーラン分布になります。このようにすると、E_2/M/2待ち行列の状態遷移は以下の図のようになります。

  • 図1

図の各楕円が状態を表し、楕円の中の2つの数字は、最初のものが到着過程のフェーズを、2番目のものが待ち行列システム内のジョブ数を表しています。また、楕円と楕円を結ぶ矢印につけられた数式は、時間間隔dtにその遷移が起こる確率密度を表しています。定常状態を求めるためには各遷移の確率密度の比だけが重要なので図1の確率密度の全てにt_eをかけて、図を下のように簡単にします。

  • 図2

ここから平衡方程式を立てると以下のようになります。なお、以下の式でp(n,m)は状態(n,m)の定常状態確率を表します。

  • 4up(0,0)=p(0,1)・・・・(16)
  • 4up(1,0)=4up(0,0)+p(1,1)・・・・(17)
  • (4u+1)p(0,1)=4up(1,0)+2p(0,2)・・・・(18)
  • (4u+1)p(1,1)=4up(0,1)+2p(1,2)・・・・(19)
  • (4u+2)p(0,2)=4up(1,1)+2p(0,3)・・・・(20)
  • (4u+2)p(1,2)=4up(0,2)+2p(1,3)・・・・(21)
  • ・・・・・・・

ところが、この連立一次方程式は簡単に解くことが出来ません。無限に一次方程式が存在するので、解くためにはシステム内のジョブ数が大きい状態の定常状態確率を無視してゼロとおいて、有限個の連立一次方程式にする必要があります。しかし、uが1に近くなると無視した状態の確率も無視出来なくなるので、この方法で定常状態確率を求めるのも考えものです。定常状態確率を求める理由は、それを用いてCT_q(E_2/M/2)を求めたいからでした。そこで作戦を変えて、全ての定常状態確率を求めることなくCT_q(E_2/M/2)を求める方法を考えてみます。