工場統計力学(建設中!) このページをアンテナに追加 RSSフィード

いまだ輝かざる暁の数は、げに多かり    (リグ・ヴェーダ)

2011-12-23

[]「待ち行列システムGI/G/1における待ちについての近似公式」の内容検討(4) 「待ち行列システムGI/G/1における待ちについての近似公式」の内容検討(4)を含むブックマーク 「待ち行列システムGI/G/1における待ちについての近似公式」の内容検討(4)のブックマークコメント

「待ち行列システムGI/G/1における待ちについての近似公式」の内容検討(3)』の続きです。

H2/M/1の場合のbの値を何とか求めることは出来ないでしょうか? ここでH2というのはバランスの取れた超指数分布を表しています。超指数分布については「バランスのとれた平均を持つ超指数分布」を参照して下さい。バランスの取れた超指数分布は1より大きい変動係数を持っています。そしてその変動係数はバランスの取れた超指数分布のパラメータを調整することでいろいろな値を取らせることが出来ます。例えば変動係数の2乗が2になるような超指数分布を考えることが出来ます。


さて、H2/M/1待ち行列はGI/M/1待ち行列の一種です。ということはGI/M/1待ち行列について述べられた

  • b=¥Bigint_0^{¥infty}¥exp¥left(-¥frac{(1-b)t}{t_e}¥right)g(t)dt・・・・(28)

が成り立ちます。(「GI/M/1待ち行列の到着時状態分布(2)」の式(16)参照。ここでは番号を振り直して式(28)とします) ここにg(t)は到着間隔の分布関数です。H2/M/1の場合、このg(t)c_a^2=2のバランスのとれた超指数分布になるわけです。そこで、まずg(t)を求め、それを式(28)に代入して整理し、そして得た式からbの値を数値計算してみようと思います。そしてその結果が私の近似式

  • b=¥frac{(c_a^2+c_e^2¥[c_a^2(1-u)+u¥])u}{(1+c_e^2)-(1-c_a^2-c_e^2¥[c_a^2(1-u)+u-1¥])u}・・・・(8)

による計算結果と、Wolfgang Kraemer氏とManfred Lagenbach-Belz氏の近似式

  • b=u+(c_a^2-1)u(1-u)h(u,c_a^2,c_e^2)・・・・(11)
  • ただし
    • h(u,c_a^2,c_e^2)=¥frac{4u}{c_a^2+u^2(4c_a^2+c_e^2)} (c_a^2{¥ge}1)・・・・(12-2)

による計算結果のどちらにより近いか比較してみようと思います。


まず、c_a^2=2のバランスの取れた超指数分布g(t)を求めます。超指数分布は以下のように表されます。

そしてバランスのとれたという言葉に対応してバランス条件

があります。この分布の変動係数c_apの間には

という関係があります。今、欲しいのはc_a^2であるような超指数分布なので、式(31)にc_a^2=2を代入して

  • p=¥frac{1}{2}¥left¥[1+¥sqrt{¥frac{2-1}{2+1}}¥right¥]=¥frac{1}{2}¥left¥[1+¥frac{1}{¥sqrt{3}}¥right¥]

よって

  • p=¥frac{1}{2}+¥frac{1}{2¥sqrt{3}}・・・・(32)

となります。

次に¥gamma_1¥gamma_2を求めます。¥gamma_1¥gamma_2については

という関係が成り立ちます。ここでE(T)は確率分布g(t)で変動する確率変数T平均値を表しています。今の場合は到着間隔の平均値です。これは少し考えれば

  • ¥frac{t_e}{u}

であることが分かります。よって

  • E(T)=¥frac{t_e}{u}・・・・(35)

となり、これを式(33)(34)に代入し、さらに式(32)を代入すると

  • ¥gamma_1=¥frac{¥left(1+¥frac{1}{¥sqrt{3}}¥right)u}{t_e}・・・・(36)
  • ¥gamma_2=¥frac{¥left(1-¥frac{1}{¥sqrt{3}}¥right)u}{t_e}・・・・(37)

となります。式(32)(36)(37)を式(29)に代入すれば欲しいg(t)の式を得るのですが、これはあまりに複雑な形になってしまいます。あまり複雑になると今後の積分の計算に支障をきたすおそれがあります。そこでpはそのままにしておくことにします。式(33)(34)に(35)を代入して

  • ¥gamma_1=¥frac{2pu}{t_e}・・・・(38)
  • ¥gamma_2=¥frac{2(1-p)u}{t_e}・・・・(39)

これらを式(29)に代入して

  • g(t)=¥frac{2p^2u}{t_e}¥exp¥left(-¥frac{2pu}{t_e}t¥right)+¥frac{2(1-p)^2u}{t_e}¥exp¥left(-¥frac{2(1-p)u}{t_e}t¥right)・・・・(40)

となります。