第六章 排队系统分析
(Queuing Systems Analysis)
6.1 排队的基本概念
6.2 到达与服务的规律
6.3 M/M/1排队模型
6.4 M/M/C排队模型
6.5 M/G/1排队模型
6.6 排队系统优化第六章 排队系统分析
6.1 排队的基本概念一、排队系统的组成队列服务机构顾客源到达 离去第六章 排队系统分析到达的顾客 服务内容 服务台收费站排队的车辆病人到达机场上空的飞机到达港口的货船进入我方阵地的敌机需要加油的车辆电话呼唤出故障的机器收费看病降落装货(卸货)
我方的防空火力加油通话修理收费站医生跑道装卸码头我方的防空导弹加油站的加油机交换台修理技工现实系统中形形色色的排队系统第六章 排队系统分析
1、输入过程
(1)顾客源:分为? 无限 (如电话呼唤)
有限 m(如车间里待修理的机器)

(2)到达规律:指到达间隔时间 T 的分布分为? 定长 D
负指数 M
k阶爱尔朗 E
k
第六章 排队系统分析
2,排队规则
(1)损失制,指顾客到达时若所有服务设施均被占用,则顾客自动离去。
(2)等待制:指顾客到达时若所有服务实施均被占用,则留下来等待,直至被服务完离去。
等待的服务规则又可分为? 先到先服务 (FCFS)
后到先服务 (LCFS)
带优先权服务(PS)
(3)混合制,分为? 系统容量有限制
等待时间有限制第六章 排队系统分析
3.服务机构
(2)服务规律:指服务时间 v 的分布分为? 定长 D
负指数 M
k阶爱尔朗 E
k
一般分布 G
(1)服务台个数 C
>
=
1
1
(并列多台)
第六章 排队系统分析二、排队模型的表示用记号(X/Y/Z/A/B/C)表示,其中
X:顾客到达时间间隔的分布
Y:服务时间的分布
Z:服务台个数
A:系统容量 N
B:顾客源数量 m
C:服务规则第六章 排队系统分析
X—表示顾客相继到达间隔时间分布,常用下列符号:
M——表示到达过程为泊松过程或负指数分布;
D——表示定长输入;
Ek——表示k阶爱尔朗分布;
G——表示一般相互独立的随机分布。
第六章 排队系统分析
Y—表示服务时间分布,所用符号与表示顾客到达间隔时间分布相同。
M——表示服务过程为泊松过程或负指数分布;
D——表示定长分布;
E
k
——表示k阶爱尔朗分布;
G——表示一般相互独立的随机分布。
第六章 排队系统分析例1 ( M / M / 1 / FCFS)表示:
到达间隔为负指数分布,服务时间也为负指数分布,1个服务台,顾客源无限,系统容量也无限,先到先服务。
//∞∞
第六章 排队系统分析三、排队问题的求解主要是计算描述系统运行状态的指标:
1,队长和排队长队长,系统中的顾客数;其概率分布称状态概率,记为 P
n
,表示系统中有 n个顾客的概率;队长的平均值记为 L
s

排队长:系统中正在排队等待的顾客数,记其均值为 L
q

第六章 排队系统分析
2,逗留时间和等待时间逗留时间:一个顾客在系统中的停留时间,记为 W,
其均值记为 W
s

等待时间:一个顾客在系统中排队等待的时间,记其均值为 W
q

第六章 排队系统分析
6.2 到达与服务规律一、到达的规律描述顾客到达规律可从两方面
到达数(人数)
到达间隔(时间)
现实中有许多服务系统,其顾客的到达具有下述特征:
( 1)无后效性:任一时段的到达数不受前一时段的影响;
( 2)平稳性:顾客到达是均匀的;( 3)稀有性:瞬时内只可能有 1个顾客到达。
第六章 排队系统分析称具有上述特征的输入为泊松流,其在 t 时段内到达 n个顾客的概率为
null,1,0,
!
)(
)( ==
ne
n
t
tP
t
n
n
λ
λ
即参数为 的泊松分布。

第六章 排队系统分析由概率论知识可知,泊松分布的参数即其均值。因此,
的含义是单位时间到达系统的平均顾客数,即到达率。
λ
下面考察,当顾客按泊松流到达时,其到达的间隔时间T 是服从什么分布呢?
tt
etFetTP
λλ
=?=≤ 1)(1)(,即而这正是负指数分布的分布函数,说明 T 服从负指数分布,且参数同为 。 λ
第六章 排队系统分析可证反之也成立。于是得到关于到达规律的重要性质:
到达数为泊松流 到达间隔服从负指数分布(同参数)。
由概率论知识可知,负指数分布的表达式(密度函数)为参数 即其均值的倒数。因此,的含义是平均间隔时间,
这与 为单位时间到达系统的平均顾客数的含义一致。
<

=
0,0
0,
)(
t
te
tf
t
T
λ
λ
λ
λ
1
λ
第六章 排队系统分析负指数分布有一个有趣的性质:无记忆性,即
00
()()P Tt tTt PTt>+ > = >
事实上,
0
0
00
00
0
()
0
0
()()
()
()
()
()
()
tt
t
t
P Tt t Tt
PT t tT t
PT t
PT t t e
ePTt
PT t e
λ
λ
λ
+
>+ >
>+ > =
>
>+
====>
>

直观上看,在已知 T >t
0
的条件下估计 T >t 的概率,与无条件时估计 T >t 的概率相同,把以前的 t
0
时间给忘了。
第六章 排队系统分析假若T表示某种电子元件的寿命,则当元件已使用了 t
0
时间后估计它还能再使用 t 时间的概率,
与刚开始用时的概率一样。说明这种元件是高度耐磨损的。
第六章 排队系统分析二、服务的规律主要讨论服务时间 v 服从负指数分布的情形,参数为,即
μ
由于 v 的均值为,即平均对每位顾客的服务时间为,可得
μ
1
μ
1
注:负指数分布的一般化 ——爱尔朗分布,可用于描述由 k 道程序组成的 1个服务台的服务时间的分布。
参数 的含义:服务率,即单位时间平均服务完 人。
μ μ
<

=
0 0,
0,
)(
t
te
tf
t
v
μ
μ
第六章 排队系统分析
6.3 M/M/1排队模型一、标准的M/M/1模型( M/M/1/ )
∞∞/
1、问题的一般提法设:泊松输入/ 负指服务/ 单服务台/ 系统无限制/ 顾客源无限制求:(1)系统状态概率 P
n;
(2)系统运行指标 L
s
,L
q
,W
s
,W
q

第六章 排队系统分析
2,系统状态概率
( 1)利用状态转移图列出平衡方程状态转移图是处理稳态M/M/C系统的一种工具,设到达与服务率分别为,则
μλ和
0
n
1 2 3
λλλλλ
μμμμμ
第六章 排队系统分析
≥+=+
=
+?
1,)(
11
10
nPPP
PP
nnn
μλμλ
μλ
( 2)由此列出平衡方程:
可解得状态概率:
≥?=
=
1),1()(
1
0
nP
P
n
n
μ
λ
μ
λ
μ
λ
第六章 排队系统分析记,称为服务强度,规定,则
ρ
μ
λ
=
1<ρ
=
=
0
0
1
PP
P
n
n
ρ
ρ
Why?
第六章 排队系统分析
3,系统运行指标
( 1) L
s
与 L
q
λμ
λ
ρ
ρ
ρ
ρρ
ρρ
ρρ
ρ
ρ
ρρ
ρ
ρ
ρρ
ρρρρρ
=
=
=
=
=?=
=?==∴
∑∑
∑∑∑

=

=

=

=

=
1
)(1
1
)(1
1
1
d
d
)(1
d
d
)(1
d
d
)(1
)(1)1(
2
01n
n
1
1
00
n
n
n
n
nn
n
ns
s
nnnpL
L 数,由期望定义,表示系统中的平均顾客∵
第六章 排队系统分析
ρ?=
=?=?=
∑∑∑

=

=

=
s
s
nnn
nnnq
L
PLPnPPnL
)1()1(
0
101
。其中 1<=
μ
λ
ρ
——因为是均值。
)呢?而不是问题:为什么 1(1 =<=? ρ
qs
LL
第六章 排队系统分析
( 2) Ws与 Wq
1
1
s
qs
W
W
WW
μλ
μλ
μ
=
=?
首先可证,逗留时间 服从参数为 的负指数分布,
而负指数分布的均值等于其参数的倒数,故平均逗留时间平均等待时间等于平均逗留时间减去平均服务时间,即第六章 排队系统分析
( 3)上述 4个指标之间的关系——Little公式
1
μμ
λ
λλ =?=?==
qsqsqqss
WWLLWLWL
。统容量无限制,故系统率。本模型中因系际进入,称有效到达率,即实应为一般的里特公式中
λλ
λλ
=
e
e
第六章 排队系统分析例2 某修理店只有一个修理工人,来修理的顾客到达数服从泊松分布,平均每小时4人;修理时间服从负指数分布,平均需6分钟。求:(1)修理店空闲的概率;(2)店内有3个顾客的概率;(3)店内至少有1个顾客的概率;(4)店内顾客的平均数;(5)
顾客在店内的平均逗留时间;(6)等待服务的顾客平均数;(7)平均等待修理时间;(8)必须在店内消耗15分钟以上的概率。
第六章 排队系统分析

人)小时小时)人人)小时小时)人
0.223)
4
1
(1)
4
1
(1)
4
1
( (8);/(
15
1
10
1
6
11
(7);/(
15
4
5
2
3
2
(6);/(
6
11
(5);/(
3
2
6
4
(4);
5
2
1 (3)
0.0384;)
5
3
()
5
2
()(1 (2);
5
3
1 (1)
1.5
4
1
4)(10
0
33
3
0
===?=<?=≥
=?=?=
=?=?=
==
==
=
=?
==?=
=?=

eeFWPWP
WW
LL
LW
L
P
P
P
sq
sq
ss
s
μ
ρ
λ
λμ
λ
ρρ
ρ

小时,人分钟人小时,人模型,解:此为标准的
5
2
/10/
6
1
/4M/M/1
==
===
μ
λ
ρ
μλ
第六章 排队系统分析二、系统容量有限的M/M/1 模型(M/M/1/ )
∞/N
1、与( M/M/1/ )的区别∞∞/
)1(
NNNe
PPP
Nn
Nnλ
Nn
=+?=

<
=
λλλ 0)(1
0,
当,
(2);10 (1)
故平均到达率

当进入系统的速率
,,,系统状态null

故速率应等于离去速率,达到统计平衡,即进入注:由于系统稳态时应
)P-(1)(1
0
μλ =?
N
P
第六章 排队系统分析
2、状态概率
=
=+=+
=
+?
NN-
nnn
PP
NnPPP
PP
μλ
μλμλ
μλ
1
1-,1,,)(
11
10
null
由此列出平衡方程:
0
n-1
1 2
λλλλ
μμμμ
N-1
λ
μ
N
n
第六章 排队系统分析
=
=
=
=+++=
=
+
=
+

0
1
0
0
0
1
0000
0
1
1
1
1
1
PP
P
PPPPPP
PP
N
N
N
N
n
N
N
n
n
n
ρ
ρ
ρ
ρ
ρ
ρρ
ρ

,可解得再由
,先解得
null
第六章 排队系统分析
3,系统运行指标




e
q
q
e
s
s
sq
N
N
N
n
ns
L
W
L
W
PLL
N
nPL
λ
λ
ρ
ρ
ρ
ρ
=
=
=
+
==
+
+
=

)(1
1
1)(
1
0
1
1
0
为有效到达率。其中 )1(
Ne
P?=λλ
第六章 排队系统分析例3 某修理站只有1个修理工,且站内最多只能停放3台待修理的机器。设待修理的机器按泊松流到达,平均每小时到达1台;修理时间服从负指数分布,平均每1.25小时可修理1台。试求:(1)站内空闲率;(2)顾客损失率;(3)有效到达率;(4)
站内平均队长;(5)机器为修理而需等待的平均时间。
第六章 排队系统分析
。,,排队系统,解:此为 25.1
8.0
1
8.01)/4/1//( ====∞ ρμλMM
。小时;台
)2.23(
0.702
0.122)(12.44)(1
(5)
)2.44(
1.251
1.255
1.251
1.25
1
1)(4
1
(4)
0.702;0.298)(11)(1 (3)
0.298;0.1221.25 (2)
0.122;
1.251
1.251
1
1
(1)
0
5
5
5
5
4
4
0
4
4
514
0
=

=

==
=
×
=
+
=
=?×=?=
=×==
=
=
=
+
e
s
e
q
q
s
e
PL
L
W
L
P
PP
P
λλ
ρ
ρ
ρ
ρ
λλ
ρ
ρ
ρ
第六章 排队系统分析例4:为开办一个小型汽车冲洗站,必须决定提供等待汽车使用的场地大小。设要冲洗的汽车到达服从泊松分布,平均每4分钟1辆,冲洗的时间服从负指数分布,平均每3分钟洗1辆。试计算当所提供的场地仅能容纳(a)1辆;(b)3辆;(c)5辆(包括正在被冲洗的1辆)时,由于等待场地不足而转向其它冲洗站的汽车的比例。
第六章 排队系统分析
。分钟,辆分钟,辆,解:
4
3
/
3
1
/
4
1
4
1
==== ρμλ
λ
428.0
7
3
,
7
4
16
9
1
4
3
1
1
1
,1 )(
01
2
0
====
=
== PPPNa ρ
ρ
ρ
154.0
175
27
,
175
64
256
81
1
4
3
1
1
1
,3 )(
0
3
3
4
0
====
=
== PPPNb ρ
ρ
ρ
072.0304.0237.0
,304.0
178.01
4
3
1
1
1
,5 )(
0
5
5
6
0
=×==
=
=
==
PP
PNc
ρ
ρ
ρ
第六章 排队系统分析三、顾客源有限的M/M/1 模型(M/M/1/ )
/m∞
1、与( M/M/1/ )的区别∞∞/
)/(;
:/

∞∞
=
)。次(每人率同即单位时间每人到率,表每个顾客的平均到达人(每人率不同)即单位时间平均到率,表全体顾客的平均到达)(
的含义
,,,系统状态
λλ
λ
λ
,
(2);10 (1)
m
mn null
第六章 排队系统分析
=
=?=
=
=
<?

∞∞
,)(
)()()(
])[(
,,0
,,
/
/
)3(
s
e
Lm
nEmnEmE
nmE
mn
mnnm
m
λ
λλλλ
λλ
λ
λ
故平均到达率
(为什么?))(
):(
(与状态无关);):(
实际进入率

第六章 排队系统分析说明(进入率与状态有关):
如m=5,n=3,如下图所示
null
null
null
丙乙甲
λ3
nullnull
进入的或甲或乙或丙,故 λ3
第六章 排队系统分析
2、状态概率
0 m-11 2
mλ(m-1)λ
λ
μμμ
m
(m-n+1) λ
n-1
μ
n
=
=+=++?
=
+?
1
11
10
1-,1,,])([)1(
mm
nnn
PP
mnPm-nPPnm
PPm
μλ
μλμλ
μλ
null
由此列出平衡方程:
第六章 排队系统分析
=
=
=

=
mnP
im
m
P
im
m
P
n
n
m
i
i
,,1,)(
)!(
!
)(
)!(
!
1
0
0
0
null
μ
λ
μ
λ
解得第六章 排队系统分析
3、系统运行指标
).1(
,1;
)(
),(
0
0
PmL
PLL
Lm
LL
Lm
s
qs
se
qs
se
=
=?
==?∴
=
λ
μ
μ
λ
μ
λ
λλ
解得另一方面,
由里特公式,

问题,的直观意义为何?
s
L
第六章 排队系统分析例5,某车间有5台机器,每台机器的连续运转时间服从负指数分布,平均连续运转时间为15分钟。有1
个修理工,每次修理时间服从负指数分布,平均每次需12分钟。
求(1)修理工空闲的概率;(2)5台机器都出故障的概率;(3)出故障机器的平均台数;(4)等待修理机器的平均台数;(5)每台机器的平均停工时间;(6)每台机器的平均等待修理时间。
第六章 排队系统分析;台 )3.76()0073.01(
8.0
1
5 (3)
0.287;(0.8)
!0
!5
(2)
0.0073](0.8)
!0
!5
(0.8)
!1
!5
(0.8)
!2
!5
(0.8)
!3
!5
(0.8)
!4
!5
(0.8)
!5
!5
[ (1)
0
5
5
1-543210
0
==
==
=+++++=
s
L
PP
P
=
=
=

=
mnP
im
m
P
im
m
P
n
n
m
i
i
,,1,)(
)!(
!
)(
)!(
!
1
0
0
0
null
μ
λ
μ
λ
。,,排队系统,解:此为8.0
15
12
12
1
15
1
)5//1//( ====∞ ρμλMM
第六章 排队系统分析
。分钟分钟;台
)(431246 (6)
)(4615
)0073.01(
12
1
5
)5(
)(77.2)0073.01(76.3 )4(
=?=
=?
=
==
q
s
q
W
W
L
由此可对该排队系统做何分析?
——机器停工时间过长,修理工几乎没有空闲时间应当提高服务率或增加修理工,或购置高效机器减少需修理率。
第六章 排队系统分析一、前提,单队、并列C 台
6.4 M/M/C排队模型
2
……
1
C
μ
μ
μ
null
(/ /):
( / / )
(/ /)
G
G
G
N
m
∞∞?


标准的模型仍可分为我们仅讨论标准 的M/M/C
第六章 排队系统分析二、(M/M/C),系统
//()G∞ ∞
服务率与服务强度:
(/ /1),
1
(/ /):
C
=
C
MM
nnC
MMC
nC
μ
λ
ρ
μ
μ
μ
λ
ρ
μ
<
<

服务率与系统状态无关,皆为,服务强度
=

服务率与系统状态有关,为

,表示每台单位时间内的平均负荷
1、与 (M/M/1/ )的区别//G∞∞
第六章 排队系统分析
2、状态概率由此列出平衡方程:
0
n
1 2
n-1
λλλλ
μ



n+1
CC-1
λλ
cμcμ
C+1
λλ
cμcμ
λλ
第六章 排队系统分析
01
11
11
(1) ( )
()
nn n
n- n n
PP
P nP nPnC
PCP CPnC
λμ
λμ
λμλμ
+
+
=
+ +=+ <
+=+ ≥


先解得,
102 13 2
23
PPP PP P
λλλ
μμμ
== = nullnull,,,
10n1
1
()
!
C
CC n
P PPPPnC
CC C
λλ λ
μμ μ

== = >,()
第六章 排队系统分析
0
0
1
1
0
0
0
0
1,
()1 ()
1
!1 !
1
( )
!
1
()
!
n
n
Cn
C
n
n
n
n
nC
PP
CC
P
C
PnC
n
P
PnC
CC
ρρλ
ρ
ρμ
λ
μ
λ
μ

=
=
=

=+ =<?


<

=



再由 解出 得:



第六章 排队系统分析
3.运行指标
00 1 1
0
0
2
1
()
( )
() ()
!(1)
C
snn n n
n n nC nC
C
qnq
n
C
qn
nC
q
s
sq
L nP nP n C P CP
LC CnPL
P
LnCPC
C
L
L
WW
λ
μ
ρ
ρ
ρ
λ λ
∞∞∞
= = =+ =+
=

=+
==+?+
=+ =+
=?=
==
∑∑∑ ∑



第六章 排队系统分析
[]
[]
λ
λμ
μ
μμ
λ
=
∑∑

CC
nn
n=0 n=0
n=0 n=0
C
n
n=0
(1)解释C- (C-n)P 的直观意义:此式即 = C- (C-n)P,
其中 (C-n)P为平均闲着的台数,C- (C-n)P为平均忙着的台
数,为每台的服务率,C- (C-n)P 为系统的平均服务率,
由统计平衡,它应等于平均到达率 。
//1
//1
q
ss
q
s
sq
LL MM
LL
L
L
WW MM
λ
μ
λ
μ
λλ λ μ
=
+
s
(2),这里与 结果相同,而与C无关,从而
1
- = - = = 也与 相同注:
第六章 排队系统分析
0.1111 0.0101 0.0014 0.0002
0.2500 0.0417 0.0103 0.0030
0.4286 0.0989 0.0333 0.0132
C
q
W μ
C
λ
μ
服务台数C
C=1 C=2 C=3 C=4
0.3
0.2
0.1
…… ……………… ……
(3) M/M/C指标有表可查:
第六章 排队系统分析
(4) 单队C台与C个单队单台系统比较
λ μ设C=2,=4,=5
------显然,单队C 台效率高 !
0.8
0.8
()
q
W
λ
ρ
μ
λ
μμ λ
==
==

()a
4λ=
4λ =
5μ =
5μ =
……
……
5μ =
()b
8λ=
5μ =
1
0
2
0.8
2
0.35
!(1 )
C
q
C
P
W
C
λ
ρ
μ
λρ
μρ
==
==

……
第六章 排队系统分析
6.5 M/G/1排队模型以上讨论了M/M/1和M/M/C系统,其前提均为泊松输入和负指数服务处理,这类系统的工具是生灭工程状态转移图。在实际中,有时到达仍为泊松过程,但服务时间并不服从负指数分布,即
M/G/1系统 这时不能用生灭过程处理,而主要依据布拉切克-钦辛公式(P-K公式)。
第六章 排队系统分析一、(M/G/1)
:
//1
E
EMM
υυσυ
ρλ υ
2
服务时间 服从任意分布,( )与 ( )存在并已
知,服务强度 = ( )<1。其他条件同 。
:系统运设求行指标。
第六章 排队系统分析
222
()
2(1 )
S
S
SqS qq
L
L
WWWELW
ρλσυ
ρ
ρ
υ λ
λ
+
=+
==?=
:由布拉切克-钦辛 (P-K)公式:
由里特公式:
,(),
解第六章 排队系统分析二、(M/D/1)系统(定长服务时间)
0Eυυ συ≡ =
2
这时 ( ),( )
2
2(1 )
s
L
ρ
ρ
ρ
∴ =+
1
:()E υ
μ
=若设
2
2( )
s
L
λλ
μ μμ λ
=+

2
,
2( ) 2( )
//1
q
qs q
L
LL W
MM
λλ λ
μ μμ λ λ μμ λ
∴ =?= = =

均为 相应指标的一半。
可见,内部越有规律越省时间第六章 排队系统分析三、(M/ /1 ) (k阶爱尔郎服务时间)
1
k
ii
i
υυ υ
=
=

设,每个 服从同参数的负指数分布
k
E
2
2
11
() ()EE
k
λ
υσ ρλυ
μ μμ
==于是,,令 = ( )=
2
2
22
11
()
(1)
2(1 )
2(1 )
s
k
k
L
k
λ
λ ρ
μμ
ρ
λ
μ ρ
μ
+
+
=+ =+

(1)
2(1 )
q
s
qs
L
Lk
LWW
k
ρ
ρ λλ
+
===
,,
由 P-K公式,
可见,k=1时即(M/M/1),k 时即(M/D/1) →∞
由里特公式,
注,对于到达与服务均为任意分布的情况,可采用随机模拟的方法求近似解。
第六章 排队系统分析
6.6 排队系统最优化
μ?
系统设计优化(静态):如何设计一个系统(如何定,C,服务
规则)使费用最经济(未必最小)。
系统控制优化(动态):一个给定的系统如何运营使目标最优。
排队系统优化
我们主要研究静态优化,
目标,费用( 损失) 最小。
*
C
μ
*
服务率待定指标服务台数第六章 排队系统分析一、标准的M/M/1系统的最优服务率
:,
sw
CC为对每个顾客的单位时间服务费 为每个顾客在
系统停留单位时间的损失费,z为总费用。

,( )
sws s
Min z C C L L
λ
μ
μλ
=+ =
单位时间费用最小,,目标
2
*
,0
()
sw
w
s
dz
CC
d
C
C
λ
μμ
μλ λ
=?=
=+
解令得
*
μ
*
μ μ
z
s
C μ
ws
CL
z
0
:,(1)
s
Max z P G C
G
μμ=若目标为服务利润最大 表达式为
其中 为单位时间对每位顾客服务注的收入。
第六章 排队系统分析二.M/M/C系统的最佳服务台数
:
s
C
w
为每台单位时间服务成本,C 为每个顾客在
系统停留单位时间的损失费,z为总费用。

(),
sws
Min z C C C L=+ 单位时间费用最小,目标
*
**
**
** * *
:()
() ( 1)
() ( 1)
,( ) ( 1) ( 1) ( )
s
w
zzC C
C
zC zC
zC zC
C
LC LC LC LC
C
=
≤?

≤+
+≤≤

解 是 的离散函数,不能求导,采用边际分析法。
是极小点
应满足化简得
*
C