更多资料获取在文章末尾!!!!!!!!!!!!!!!!
第六章 排队论模型 |
排队论起源于 1909 年丹麦电话工程师 A. K.爱尔朗的工作,他对电话通话拥挤问 题进行了研究。1917 年,爱尔朗发表了他的著名的文章—“自动电话交换中的概率理 论的几个问题的解决”。排队论已广泛应用于解决军事、运输、维修、生产、服务、库 存、医疗卫生、教育、水利灌溉之类的排队系统的问题,显示了强大的生命力。 排队是在日常生活中经常遇到的现象,如顾客到商店购买物品、病人到医院看病常 常要排队。此时要求服务的数量超过服务机构(服务台、服务员等)的容量。也就是说, 到达的顾客不能立即得到服务,因而出现了排队现象。这种现象不仅在个人日常生活中 出现,电话局的占线问题,车站、码头等交通枢纽的车船堵塞和疏导,故障机器的停机 待修,水库的存贮调节等都是有形或无形的排队现象。由于顾客到达和服务时间的随机 性。可以说排队现象几乎是不可避免的。 |
排队论(Queuing Theory)也称随机服务系统理论,就是为解决上述问题而发展 的一门学科。它研究的内容有下列三部分: |
(i)性态问题,即研究各种排队系统的概率规律性,主要是研究队长分布、等待 时间分布和忙期分布等,包括了瞬态和稳态两种情形。 |
(ii)最优化问题,又分静态最优和动态最优,前者指最优设计。后者指现有排队 系统的最优运营。 |
(iii)排队系统的统计推断,即判断一个给定的排队系统符合于哪种模型,以便 根据排队理论进行分析研究。 这里将介绍排队论的一些基本知识,分析几个常见的排队模型。 |
§1 基本概念 |
1.1 排队过程的一般表示 下图是排队论的一般模型。 |
图1 排队模型 |
图中虚线所包含的部分为排队系统。各个顾客从顾客源出发,随机地来到服务机构,按 一定的排队规则等待服务,直到按一定的服务规则接受完服务后离开排队系统。 凡要求服务的对象统称为顾客,为顾客服务的人或物称为服务员,由顾客和服务员 组成服务系统。对于一个服务系统来说,如果服务机构过小,以致不能满足要求服务的 众多顾客的需要,那么就会产生拥挤现象而使服务质量降低。 因此,顾客总希望服务 机构越大越好,但是,如果服务机构过大,人力和物力方面的开支也就相应增加,从而 会造成浪费,因此研究排队模型的目的就是要在顾客需要和服务机构的规模之间进行权 衡决策,使其达到合理的平衡。 |
1.2 排队系统的组成和特征 一般的排队过程都由输入过程、排队规则、服务过程三部分组成,现分述如下: 1.2.1 输入过程 |
输入过程是指顾客到来时间的规律性,可能有下列不同情况: (i)顾客的组成可能是有限的,也可能是无限的。 |
-118- |
(ii)顾客到达的方式可能是一个—个的,也可能是成批的。 (iii)顾客到达可以是相互独立的,即以前的到达情况对以后的到达没有影响; |
否则是相关的。 |
(iv)输入过程可以是平稳的,即相继到达的间隔时间分布及其数学期望、方差等 数字特征都与时间无关,否则是非平稳的。 1.2.2 排队规则 |
排队规则指到达排队系统的顾客按怎样的规则排队等待,可分为损失制,等待制和 混合制三种。 |
(i)损失制(消失制)。当顾客到达时,所有的服务台均被占用,顾客随即离去。 (ii)等待制。当顾客到达时,所有的服务台均被占用,顾客就排队等待,直到接 受完服务才离去。例如出故障的机器排队等待维修就是这种情况。 (iii)混合制。介于损失制和等待制之间的是混合制,即既有等待又有损失。有 队列长度有限和排队等待时间有限两种情况,在限度以内就排队等待,超过一定限度就 离去。 |
排队方式还分为单列、多列和循环队列。 1.2.3 服务过程 |
(i)服务机构。主要有以下几种类型:单服务台;多服务台并联(每个服务台同 时为不同顾客服务);多服务台串联(多服务台依次为同一顾客服务);混合型。 (ii)服务规则。按为顾客服务的次序采用以下几种规则: ①先到先服务,这是通常的情形。 |
②后到先服务,如情报系统中,最后到的情报信息往往最有价值,因而常被优先处 理。 |
③随机服务,服务台从等待的顾客中随机地取其一进行服务,而不管到达的先后。 ④优先服务,如医疗系统对病情严重的病人给予优先治疗。 1.3 排队模型的符号表示 |
排队模型用六个符号表示,在符号之间用斜线隔开,即X /Y / Z / A/ B /C 。第一 个符号X 表示顾客到达流或顾客到达间隔时间的分布;第二个符号Y 表示服务时间的 分布;第三个符号Z 表示服务台数目;第四个符号A 是系统容量限制;第五个符号B 是 顾客源数目;第六个符号C 是服务规则,如先到先服务 FCFS,后到先服务 LCFS 等。并 约定,如略去后三项,即指X /Y / Z / ∞ / ∞ / FCFS的情形。我们只讨论先到先服务 FCFS 的情形,所以略去第六项。 |
表示顾客到达间隔时间和服务时间的分布的约定符号为: M —指数分布(M 是 Markov 的字头,因为指数分布具有无记忆性,即 Markov 性); |
D —确定型(Deterministic); Ek —k 阶爱尔朗(Erlang)分布; G —一般(general)服务时间的分布; |
GI —一般相互独立(General Independent)的时间间隔的分布。 例如,M / M /1表示相继到达间隔时间为指数分布、服务时间为指数分布、单服 务台、等待制系统。D / M / c 表示确定的到达时间、服务时间为指数分布、c 个平行 服务台(但顾客是一队)的模型。 |
1.4 排队系统的运行指标 为了研究排队系统运行的效率,估计其服务质量,确定系统的最优参数,评价系统 的结构是否合理并研究其改进的措施,必须确定用以判断系统运行优劣的基本数量指 |
-119- |
标,这些数量指标通常是: (i)平均队长:指系统内顾客数(包括正被服务的顾客与排队等待服务的顾客)的 数学期望,记作Ls 。 |
(ii)平均排队长:指系统内等待服务的顾客数的数学期望,记作Lq 。 |
(iii)平均逗留时间:顾客在系统内逗留时间(包括排队等待的时间和接受服务的 时间)的数学期望,记作Ws 。 (iv)平均等待时间:指一个顾客在排队系统中排队等待时间的数学期望,记作 |
Wq 。 |
(v)平均忙期:指服务机构连续繁忙时间(顾客到达空闲服务机构起,到服务机 构再次空闲止的时间)长度的数学期望,记为Tb 。 |
还有由于顾客被拒绝而使企业受到损失的损失率以及以后经常遇到的服务强度等, 这些都是很重要的指标。 |
计算这些指标的基础是表达系统状态的概率。所谓系统的状态即指系统中顾客数, 如果系统中有n 个顾客就说系统的状态是n ,它的可能值是 (i)队长没有限制时,n = 0,1,2,L, |
(ii)队长有限制,最大数为N 时,n = 0,1,L,N , (iii)损失制,服务台个数是c 时,n = 0,1,L,c 。 |
这些状态的概率一般是随时刻t 而变化,所以在时刻t 、系统状态为n 的概率用 P (t) 表示。稳态时系统状态为n 的概率用P 表示。 |
n |
n |
§2 输入过程与服务时间的分布 |
排队系统中的事件流包括顾客到达流和服务时间流。由于顾客到达的间隔时间和服 务时间不可能是负值,因此,它的分布是非负随机变量的分布。最常用的分布有泊松分 布、确定型分布,指数分布和爱尔朗分布。 |
2.1 泊松流与指数分布 设N(t) 表示在时间区间[0,t) 内到达的顾客数(t > 0),令P (t ,t )表示在时间区 |
正如在概率论中所学过的,我们说随机变量{N(t) = N(s + t) − N(s)}服从泊松分 |
布。它的数学期望和方差分别是 |
E[N(t)] = λt ;Var[N(t)] = λt 。 当输入过程是泊松流时,那么顾客相继到达的时间间隔T 必服从指数分布。这是 |
由于 |
P{T > t} = P{[0,t) 内呼叫次数为零} = P (t) = e−λ |
t |
0 |
那么,以F(t) 表示T 的分布函数,则有 |
⎧ − −λt |
≥ 0 |
1 e , t |
P{T ≤ t} = F(t) = |
⎨ |
0, |
t < 0 |
⎩ |
而分布密度函数为 |
f (t) = e , t > 0 . |
λ |
−λt |
-121 |
1 |
对于泊松流,λ 表示单位时间平均到达的顾客数,所以 就表示相继顾客到达平均 |
λ |
间隔时间,而这正和ET 的意义相符。 |
对一顾客的服务时间也就是在忙期相继离开系统的两顾客的间隔时间,有时也服从 指数分布。这时设它的分布函数和密度函数分别是 |
( ) = 1− −μ ,g t G t |
e |
t |
( ) = μe−μ |
t |
我们得到 |
P{T ≤ t + Δt |T > t} Δt |
P{t < T ≤ t + Δt} ΔtP{T > t} |
lim |
= lim |
= μ |
Δt→0 |
Δt→0 |
这表明,在任何小的时间间隔[t,t + Δt) 内一个顾客被服务完了(离去)的概率是 |
1 μ |
μΔt + o(Δt) 。μ 表示单位时间能被服务完成的顾客数,称为平均服务率,而 表示 |
一个顾客的平均服务时间。 2.2 常用的几种概率分布及其产生 2.2.1 常用的连续型概率分布 |
我们只给出这些分布的参数、记号和通常的应用范围,更详细的内容参看专门的概 率论书籍。 (i)均匀分布 区间(a,b) 内的均匀分布记作U(a,b) 。服从U(0,1) 分布的随机变量又称为随机 数,它是产生其它随机变量的基础。如若X 为U(0,1) 分布,则Y = a + (b − a)X 服从 |
U(a,b) 。 |
(ii)正态分布 |
以μ 为期望,σ 2 为方差的正态分布记作N(μ,σ ) 。正态分布的应用十分广泛。 |
2 |
正态分布还可以作为二项分布一定条件下的近似。 (iii)指数分布 |
指数分布是单参数λ 的非对称分布,记作Exp(λ),概率密度函数为: |
⎧λ |
e−λt , t |
≥ 0 |
f (t) = |
⎨ |
⎩ 0, |
t < 0 |
1 λ |
1 λ 2 |
它的数学期望为 ,方差为 。指数分布是唯一具有无记忆性的连续型随机变量,即 |
有P(X > t + s | X > t) = P(X > s) ,在排队论、可靠性分析中有广泛应用。 |
(iv)Gamma 分布 Gamma 分布是双参数α,β 的非对称分布,记作G(α,β ) ,期望是αβ 。α = 1时蜕 化为指数分布。n 个相互独立、同分布(参数λ )的指数分布之和是 Gamma 分布 (α = n,β = λ) 。Gamma 分布可用于服务时间,零件寿命等。 |
Gamma 分布又称爱尔朗分布。 (v)Weibull 分布 Weibull 分布是双参数α,β 的非对称分布,记作W(α,β ) 。α = 1时蜕化为指数分 |
布。作为设备、零件的寿命分布在可靠性分析中有着非常广泛的应用。 (vi)Beta 分布 |
-122 |
Beta 分布是区间(0,1) 内的双参数、非均匀分布,记作B(α, β ) 。 |
2.2.2 常用的离散型概率分布 (i)离散均匀分布 |
(ii)Bernoulli 分布(两点分布) Bernoulli 分布是 x = 1,0 处取值的概率分别是 p 和1− p 的两点分布,记作 Bern( p) 。用于基本的离散模型。 |
(iii)泊松(Poisson)分布 |
泊松分布与指数分布有密切的关系。当顾客平均到达率为常数λ 的到达间隔服从 指数分布时,单位时间内到达的顾客数K 服从泊松分布,即单位时间内到达k 位顾客 的概率为 |
记作B(n, p) 。二项分布是n 个独立的 Bernoulli 分布之和。它在产品检验、保险、生 |
§3 生灭过程 一类非常重要且广泛存在的排队系统是生灭过程排队系统。生灭过程是一类特殊的 随机过程,在生物学、物理学、运筹学中有广泛的应用。在排队论中,如果N(t) 表示 |
时刻t 系统中的顾客数,则{N(t),t ≥ 0}就构成了一个随机过程。如果用“生”表示顾 |
客的到达,“灭”表示顾客的离去,则对许多排队过程来说,{N(t),t ≥ 0}就是一类特 殊的随机过程-生灭过程。下面结合排队论的术语给出生灭过程的定义。 定义1 设{N(t),t ≥ 0}为一个随机过程。若N(t) 的概率分布具有以下性质: (1)假设N(t) = n,则从时刻t 起到下一个顾客到达时刻止的时间服从参数为λn 的负指数分布,n = 0,1,2,L。 (2)假设N(t) = n,则从时刻t 起到下一个顾客离去时刻止的时间服从参数为μn 的负指数分别,n = 0,1,2,L。 |
(3)同一时刻只有一个顾客到达或离去。 则称{N(t),t ≥ 0}为一个生灭过程。 一般来说,得到N(t) 的分布 pn (t) = P{N(t) = n}(n = 0,1,2,L)是比较困难的, 因此通常是求当系统到达平衡后的状态分布,记为 pn ,n = 0,1,2,L。 为求平稳分布,考虑系统可能处的任一状态n 。假设记录了一段时间内系统进入状 态n 和离开状态n 的次数,则因为“进入”和“离开”是交替发生的,所以这两个数要 |
-123- |