wq
-3. M/M/S
待ち行列,
アーラン
の即時式
モデル
1
金子邦彦
(待ち行列の数理)
URL:
https://www
.kkaneko.jp/cc/wq/index.html
ア
ウ
トライン
3-1
M/M/
S
待ち行列
3-
2
アー
ランの
即時式モデル
2
3-1 M/M/S
待ち行列
3
X:
到着過程
→
ポアソン分布
のとき「
M
」と書く
Y
:
処理時間分布
→
指数分布
のとき「
M
」と書く
Z:
サーバ数
個数を限定しな
いとき「
S
」
と書く
K
:
待ち行列の長さ
の制
限
制限しないときは何
も書かない
4
ケンドール記法
M/M/S
M/M/S
待ち行列モデル
サーバ数:
S
待ち行列
サーバ
ポアソン分布
で到着
処理時間は
指数分布
5
M/M/S
待ち行列モデルの解析
•
状態遷移
•
システム内のジョブ数を「状態」と見る
•
定常状態
•
定常状態での、各「状態」の確率を求める
•
システムの「処理率
」を求める
6
状態
状態0:
システム
内のジョブ数が
0
状態1:
システム内のジョブ数が
1
状態
S:
システム内のジョ
ブ数が
S
7
状態遷移
•
ジョブが到着:
状態
k
から状態
k
+
1
に遷移
•
ジョブが完了:
状態
k
から状態
kー1
に遷移
8
遷移確率に関する方程式
•
微小時間
⊿
t(限
りな
く0に近い
)につい
ての式
•
「ジョブの
到着」と「ジョブ
の完了」
は、
「同時」に
は起きない
•
「ジョ
ブの
完了」の確率
は、
⊿
t
と
μ
の式
•
「ジョ
ブの
到着」の確率
は、
⊿
t
と
λ
の式
9
•
k
≦
S
ならば(
kは状態番号)
•
サーバに空きがある
k
μ
⊿
t
(1
ー
μ
⊿
t)
≒
k
μ
⊿
t
•
k>
S
ならば(
kは状態番号)
•
サーバは全て忙しい
Sμ
⊿
t
(1
ー
μ
⊿
t)
≒
Sμ
⊿
t
k-1
S-1
10
「ジョブの完了」に関する
式
•
ジョブが到着
:
状態
k
から状態
k
+
1
に遷移
λ
⊿
t
•
ジョブが完了:
状態
k
から状態
k-1
に遷移
k
≦
S
ならば:
k
μ
⊿
t
k>
S
ならば:
Sμ
⊿
t
11
状態遷移
状態遷移図
1-
λ
状態0
状態1
λ
μ
1-
μ
-
λ
λ
2
μ
状態
S-1
状態
S
λ
Sμ
1-
Sμ
-
λ
λ
Sμ
λ
Sμ
状態
S
+
1
1-
Sμ
-
λ
1-(
S
-1)
μ
-
λ
12
定常状態方程式
0<n<
S
では
(
λ
+n
μ
)
P
n
= λPn
ー
1
+
(n+1)
μPn
+1
S
≦
n
では
(
λ
+
Sμ
)
P
n
= λPn
ー
1
+ SμPn
+1
λP
0
= μP
1
13
定常確率を
P0
で表す
0<n<
S
では
Pn =
S
≦
n
では
Pn =
λ
P
0
μ
n!
( )
n
λ
P
0
μ
S!S
( )
n
n-
S
14
システム処理能力
ρ
ρ= λ/μ
•
λ
⊿
t
:
「時間
(
t,
t
+
⊿
t
)
に到
着するジョ
ブ数」の平
均
•
μ
⊿
t
:
「サーバが
ジョブを処理中の間,
⊿
t
内に完了する処理数
」の平均
15
待ち合わせが生じる確率
•
システム内のジョブ
数が
S
以上
•
Pqueue
= ∑
Pn
= ∑
= ∑
P
0
n=S
∞
λ
P
0
μ
S!S
( )
n
n-
S
n=S
∞
n=S
∞
(
Sρ
)
1
S!
S
n
n-
S
16
3-
2
アーラン
の即時式モデル
17
X:
到着過程
→
ポアソン分布
のとき「
M
」と書く
Y
:
処理時間分布
→
指数分布
のとき「
M
」と書く
Z:
サーバ数
個数を限定しな
いとき「
S
」
と書く
K
:
待ち行列の長さ
の制
限
K = 1
18
ケンドール記法
M/M/S/1
M/M/S/1
のとき
サーバ数:
S
K=1
なので
待ち行列なし
処理時間は
指数分布
ポアソン分布
で到着
19
アーランの即時式モデル
•
待ち行列モデル:
M/M/S/S
•
入力回線数:
無限
•
待ち行列長:
L
=0
•
システム内の最大占
有サーバ数:
K
=
S
20
「即時式」の意味
K
=
1
なの
で,次の性質を持つ
•
待ち行列なし
•
待ち行列長は
0
•
システム内の占有サ
ーバ数が
S
(すべてのサー
バ
が占有されている)
のとき,到着したジョブは棄
却される
•
サーバの空きがあれ
ば,直ちに処理される
→
即時式
21
•
ジョ
ブの
到着
•
平均
λ
のポアソン分布
→
ある時刻に客が到着してから時間
t
内に次の客が到着
する確率はポアソン分布に従
う
:
1-
e
-
λ
t
•
微少時間
⊿
t
の間にジョブが到着する確率:
λ
⊿
t
•
ジョブの処
理時間
•
平均1
/μ
の指数分布
•
微少時間
⊿
t
の間に処理が終わる確率:
μ
⊿
t
システム処理能力:
ρ
=
λ/μ
22
ジョブのモデル
状態
状態0:
系内のジョブの数が
0
状態1:
系内のジョブの数が
1
状態
S:
系内のジョブの数が
S
(状態は
S
まで)
23
•
ジョブが到着:
1)
k<
S
のとき
状態
k
から状態
k
+
1
に遷移
2)
k
=S
のとき
状態
S
のまま(新しい到着は棄却される)
•
ジョブの回線の占有終わり:
状態
k
から状態
kー1
に遷移
24
状態遷移
•
ジョブが到着:
1)
k<
S
のとき
状態
k
から状態
k
+
1
に遷移:
λ
⊿
t
•
ジョブの回線の占
有終わり:
状態
k
から状態
k
ー1
に遷移
:
k
μ
⊿
t
25
状態遷移
状態遷移図
1-
λ
状態0
状態1
λ
μ
1-
μ
-
λ
λ
2
μ
状態
S-1
状態
S
λ
Sμ
μ
-
Sμ
1-(
S
-1)
μ
-
λ
λ
(
S-
1)μ
26
定常状態についての方程式
(
λ
+n
μ
)
P
n
= λPn
ー
1
+
(n+1)
μPn
+1
1
λP
0
= μP
1
λP
n-1
= SμP
n
27
•
定常状態で,系内に
n
個のジョブ(
0
≦
n
≦
s
)がある
確率を
Pn
とすると
•
棄却率:
•
システム内の占有サーバ数が
S
になる確率:
Ps
•
上の式で,
n=s
として計算する
!
2
1
!
)
(
2
s
n
t
P
s
n
n
28
棄却率