M/M/1における待ち時間の式の導出(2)
ここでは、M/M/m待ち行列での考察をなぞる方法を紹介します。
到着時間の間隔が指数分布であり、処理時間が指数分布であることから、任意の時刻からの間にロットが到着する確率や、装置の処理が終了する確率は常に一定です(記憶なし特性)。平均処理時間はなので、の間に装置の処理が完了する確率は
- ・・・・・・(1)
となります。一方、装置の利用率が100%(つまり1)の時の装置にやってくるロットの到着間隔の平均は、平均処理時間と等しくなりますからになります。さらに、ロットの到着間隔は装置の利用率と反比例しますから
となります。よって、の間にロットが到着する確率は
- ・・・・・・(2)
次に、システム(待ち行列と装置を合わせた全体)内にあるロットの数が個である確率をとします。またシステム内にロットが個ある状態を状態と呼ぶことにします。定常状態が存在すると仮定して、定常状態のの値をこれから求めます。
まず、の確率を考えます。任意の時刻からの間に他の状態から状態0に入る遷移は、状態1から、ロットの処理が終了して状態0になる遷移しかありません。また、任意の時刻からの間に状態0から他の状態に出ていく遷移は、ロットが到着して状態1に出ていく遷移しかありません。それぞれの遷移の確率は、状態1→状態0が、式(1)から
- ・・・・・・(3)
状態0→状態1が、式(2)から
- ・・・・・・(4)
定常状態ではの増減がないので、(3)と(4)は等しくなければなりません。
よって、
- ・・・・・・(5)
となります。
今度は、状態1と状態2の関係を考えてみます。同様に式(1)と(2)を用いると下図のような結果になります。
ところが、
であることが判明していますので、が一定であるためには今度は
- 状態2→状態1の遷移確率 と
- 状態1→状態2の遷移確率
が等しくなければなりません。よって
- ・・・・・・(6)
以下、同様に考えれば、
- ・・・・・・(7)
になります。式(7)から
- ・・・・・・(8)
が導かれます。さて、ここでは、装置が空いている確率ですからになります。よって
- ・・・・・・(9)
これで、定常状態での各状態の確率が求まりました。ここから、待ち行列に並んでいるロットの数の平均、すなわち行列の長さの平均を求めることが出来ます。待ち行列に1個ロットが並んでいる状態は、装置がロットを1個処理中であって、待ち行列に1個並んでいる状態ですので、その確率はになります。同様に考えて、待ち行列に個ロットが並んでいる状態の発生確率はになります。
よって、行列の平均の長さは
- ・・・・・・(10)
で与えられることになります。式(10)に式(9)を代入して
- ・・・・・・(11)
ここで
- ・・・・・・(12)
とおくと
- ・・・・・・(13)
よって式(12)の両辺から式(13)の両辺を引くと
よって
これを式(11)に代入して
- ・・・・・・(14)
これで行列の平均の長さが求まりました。次に平均待ち時間を求めるために、リトルの法則
を利用します。(平均WIP)にあたるのがです。求めたいのは平均待ち時間ですが、これは(平均サイクルタイム)にあたります。(平均スループット)は到着間隔の逆数であることを考慮すれば、
となります。リトルの法則から
となりますので、
ここに式(14)を代入すれば、
- ・・・・・・(15)
これでM/M/1待ち行列の待ち時間の式を導出することが出来ました。