Marshallの公式に向けて(1)

待ち行列における近似モデル 逆瀬川浩孝氏」を読んでいると、こんな一節がありました。

5.平均値の近似

 混雑の指標として、何かの特性量の期待値だけを知りたい(求めたい)という場合がある。システム全体の解析を経ることなく、簡単な計算によって、そのような指標が求められるならば、いろいろなパラメータを動かした時の変化の様子を調べて最適なパラメータを選ぶという作業が容易になるであろう。たとえば単一窓口モデルGI/G/1の待ち時間を考えてみよう。n番目に到着した客の待ち時間をW_n、サービス時間をB_nn番目とn+1番目の客の到着間隔をA_nとすれば、

  • W_{n+1}=\max(0,W_n+B_n-A_n){\equiv}(W_n+B_n-A_b)^+

という関係が成り立つ。適当な条件のもとでW_nnを大きくすると平衡状態での待ち時間Wに収束し、平均待ち時間は上の式を使って、

  • EW=\frac{Var(U)}{2E(-U)}-\frac{Var(Y)}{2E(-U)}
    • (ただし、U=B-AY=-\min(0,W+U)

と表すことが出来る。

まず、GI/G/1でなぜ式(1)

  • W_{n+1}=\max(0,W_n+B_n-A_n){\equiv}(W_n+B_n-A_b)^+・・・・・(1)

が成り立つか、ですが、これは図を書けば直感的に分かります。
まずW_n+B_n{\ge}A_nの場合ですが、これは下図のようになります。

  • 図1

この時、

  • W_{n+1}=W_n+B_n-A_n・・・・・(2)

となっているのが分かります。
次に[tex:W_n+B_n

  • 図2

このときには

  • W_{n+1}=0・・・・・(3)

になっていることが分かります。(2)と(3)を一つの式にまとめると(1)になります。


私は以上は理解したのですが、次になぜ(1)ならば

  • EW=\frac{Var(U)}{2E(-U)}-\frac{Var(Y)}{2E(-U)}・・・・・(4)

が成り立つのか、まだ分かりません。

しかし(4)が成り立つとすると、Whittの論文「The Queueing Network Analyzer」で私が分からなかった式の一つである、GI/G/1待ち行列における出発間隔時間の2乗係数c_d^2についてのMarshallの公式

  • c_d^2=c_a^2+2\rho^2c_s^2-2\rho(1-\rho){\mu}EW・・・・・(5)

を導くことが出来ます。ところで、ここでは私の記述法に従って(5)を以下のように書き換えます。ただしEWは私の記述法ではCT_qですが、これはEWのままにしておきます。(そのほうが他の式を直さなくてすむので)

  • c_d^2=c_a^2+2u^2c_e^2-2u(1-u)\frac{EW}{t_e}・・・・・(6)

では、式(4)から(6)を導き出してみましょう。


まず(4)から

  • 2E(-U)EW=Var(U)-Var(Y)
  • Var(Y)=Var(U)-2E(-U)EW・・・・・(7)

となります。


Marshallの公式に向けて(2)」に続きます。