合流型待ち行列の解析

先行エントリ:大域平衡方程式と局所平衡方程式


今度は、下図のような3ステーションからなり、装置1と装置2からの流れが装置3の前で合流する待ち行列を考えます。このような待ち行列でも(処理時間が指数分布で到着間隔も指数分布ならば)積形式解になるかどうか調べます。

  • 図1

1台目の装置の平均処理時間をt_{e1}、2台目の装置の平均処理時間をt_{e2}、3台目の装置の平均処理時間をt_{e3}とします。さらに、1台目の装置の利用率u_1、2台目の装置の利用率u_2、3台目の装置の利用率u_3とします。装置1でのスループット\lambda_1、装置2でのスループット\lambda_2とします。すると

  • \lambda_1=\frac{u_1}{t_{e1}}・・・・・・(1)
  • \lambda_2=\frac{u_2}{t_{e2}}・・・・・・(2)

が成り立ちます。さらに、装置3にはこの両方の流れが合流して入って来ますので装置3でのスループット\lambda_3とすると

  • \lambda_3=\lambda_1+\lambda_2・・・・・・(3)

一方、それをt_{e3}u_3で表せば、

  • \lambda_3=\frac{u_3}{t_{e3}}・・・・・・(4)

が成り立ちます。
状態は(k,m,n)で表し、kが装置1で処理しているかあるいは待っているロットの数、mが装置2で処理しているかあるいは待っているロットの数、nが装置3で処理しているかあるいは待っているロットの数、を表すものとします。
今回は、以下の図で各状態へ入る流量(=遷移確率)と出る流量が等しいとします。(今回は合流があるので、図に示された全ての遷移の流量が等しいとは言えません。)

上図から数式を作ると、平衡状態が成り立つための式は

  • p(k,m,n)\lambda_1=p(k+1,m,n)\frac{1}{t_{e1}}・・・・・・(5)
  • p(k,m,n)\lambda_2=p(k,m+1,n)\frac{1}{t_{e2}}・・・・・・(6)
  • p(k+1,m,n)\frac{1}{t_{e1}}+p(k,m+1,n)\frac{1}{t_{e2}}=p(k,m,n+1)\frac{1}{t_{e3}}・・・・・・(7)

ここで式(7)に式(5)と式(6)を代入すると

  • p(k,m,n)\lambda_1+p(k,m,n)\lambda_2=p(k,m,n+1)\frac{1}{t_{e3}}
  • p(k,m,n)(\lambda_1+\lambda_2)=p(k,m,n+1)\frac{1}{t_{e3}}

ここで式(3)を用いて

  • p(k,m,n)\lambda_3=p(k,m,n+1)\frac{1}{t_{e3}}・・・・・・(8)

式(5)に式(1)を代入して少し変形すると

  • p(k+1,m,n)=p(k,m,n)u_1・・・・・・(9)

式(6)に式(2)を代入して少し変形すると

  • p(k,m+1,n)=p(k,m,n)u_2・・・・・・(10)

式(8)に式(4)を代入して少し変形すると

  • p(k,m,n+1)=p(k,m,n)u_3・・・・・・(11)

よって、式(9)からは

  • p(k,m,n)=p(0,m,n)u_1^k・・・・・・(12)

式(10)からは

  • p(k,m,n)=p(k,0,n)u_2^m・・・・・・(13)

式(11)からは

  • p(k,m,n)=p(k,m,0)u_3^n・・・・・・(14)

が導かれます。


式(13)からは

  • p(0,m,n)=p(0,0,n)u_2^m・・・・・・(15)

が導かれます。式(15)を式(12)に代入して

  • p(k,m,n)=p(0,0,n)u_1^ku_2^m・・・・・・(16)

となります。式(14)からは

  • p(0,0,n)=p(0,0,0)u_3^n・・・・・・(17)

が導かれます。式(17)を式(16)に代入して

  • p(k,m,n)=p(0,0,0)u_1^ku_2^mu_3^n・・・・・・(18)

となります。式(18)で分かるように、図1のような合流型の待ち行列の場合も積形式解になることが分かりました。
分岐型待ち行列の解析」に続きます。