排队论及相关程序.docx
《排队论及相关程序.docx》由会员分享，可在线阅读，更多相关《排队论及相关程序.docx（14页珍藏版）》请在冰豆网上搜索。
排队论及相关程序
实用排队论
排队论又称随机服务系统，它应用于一切服务系统，包括生产管理系统、通信系统、交通系统、计算机存储系统。
它通过建立一些数学模型，以对随机发生的需求提供服务的系统预测。
现实生活中如排队买票、病人排队就诊、轮船进港、高速路上汽车通过收费站、机器等待修理等等。
一、排队论的基本构成。
（1）输入过程。
输入过程是描述顾客是按照怎样的规律到达排队系统的。
包括①顾客总体：
顾客的来源是有限的还是无限的。
②到达的类型：
顾客到达是单个到达还是成批到达。
③相继顾客到达的时间间隔：
通常假定是相互独立同分布，有的是等间隔到达，有的是服从负指数分布，有的是服从k阶Erlang分布。
（2）排队规则排队规则指顾客按怎样的规定的次序接受服务。
常见的有等待制，损失制，混合制，闭
合制。
当一个顾客到达时所有服务台都不空闲，则此顾客排队等待直到得到服务后离开，称为等待制。
在等待制中，可以采用先到先服务，如排队买票；也有后到先服务，如天气预报；也有随机服务，如电话服务；也有有优先权的服务，如危重病人可优先看病。
当一个顾客到来时，所有服务台都不空闲，则该顾客立即离开不等待，称为损失制。
顾客排队等候的人数是有限长的，称为混合制度。
当顾客对象和服务对象相同且固定时是闭合制。
如几名维修工人固定维修某个工厂的机器就属于闭合制。
（3）服务机构
服务机构主要包括：
服务台的数量；服务时间服从的分布。
常见的有定长分布、负指数分布、几何分布等。
二、排队系统的数量指标
（1）队长与等待队长
队长（通常记为Ls）是指系统中的平均顾客数（包括正在接受服务的顾客）。
等待队
长（通常记为Lq）指系统中处于等待的顾客的数量。
显然，队长等于等待队长加上正在服务的顾客数。
（2）等待时间
等待时间包括顾客的平均逗留时间（通常记为Ws）和平均等待时间（通常记为Wq）。
顾客的平均逗留时间是指顾客进入系统到离开系统这段时间，包括等待时间和接受服务的时间。
顾客的平均等待时间是指顾客进入系统到接受服务这段时间。
（3）忙期
从顾客到达空闲的系统，服务立即开始，直到再次变为空闲，这段时间是系统连续繁忙的时期，称之为系统的忙期。
它反映了系统中服务机构工作强度，是衡量服务系统利用效率的指标，即
服务强度=忙期/服务总时间=1—闲期/服务总时间闲期与忙期对应的系统的空闲时间，也就是系统连续保持空闲的时间长度。
三、排队论中的符号表示
排队论中的记号是20世纪50年代初由D.G.Kendall引入的，通常由3~5个字母组成，形式为：
A/B/C/n
其中A表示输入过程，B代表服务时间，C代表服务台数量，n表示系统空间数。
如
（1）M/M/S/a表示输入过程是Poisson流，服务时间服从负指数分布，系统有S个服务
台平行服务，系统容量为无穷大的等待制排队系统。
（2）M/G/S/a表示输入过程是Poisson流，服务时间服从一般概率分布，系统有S个服务台平行服务，系统容量为无穷大的等待制排队系统。
（3）D/M/S/K表示顾客相继到达时间间隔独立、服从定长分布，服务时间服从负指数分布，系统有S个服务台平行服务，系统容量为K个的混合制系统。
（4）M/M/S/S表示输入过程是Poisson流，服务时间服从负指数分布，系统有S个服务
台平行服务，顾客到达后不等待的损失制系统。
（5）M/M/S/K/K表示输入过程是Poisson流，服务时间服从负指数分布，系统有S个
服务台平行服务，系统容量和顾客容量都为K个的闭合制系统。
四、排队论中四种重要的模型
的Poisson分布，在［0,t］时间内到达的顾客数
1•等待制模型M/M/S/a
该模型中顾客到达规律服从参数为
X（t）服从的的分布为:
其单位时间到达的顾客平均数为，［0,t］时间内到达的顾客平均数为t。
顾客接受服务的时间服从负指数分布，单位时间服务的顾客平均数为，服务时间的
分布为：
f（t）
（2）
1
每个顾客接受服务的平均时间为-。
下面分别给出S=1和S>1的一些主要结果。
1.1只有一个服务台的S=1情形
Pn
（1）nn0,1,2,3L
（3）
可以计算出稳定状态下系统有n个顾客的概率:
其中一称为系统的服务强度。
则系统没有顾客的概率为:
Po11-
系统中顾客的平均队长为:
Lsn.Pn
（1）n.
n0n01
（4）
系统中顾客的平均等待队长为：
2
Lq（n1）.Pn
（1）（n1）.n
n1n11
2
（5）
（）
系统中顾客的平均逗留时间为：
1
Ws
（6）
系统中顾客的平均等待时间为:
11
\/\/
（7）
VVq
q（）
从（4）~（6）式可以看出：
LsWs，LqWq
（8）
LsLq
或Ws-，Wq
（9）
该公式称为Little公式。
在其它排队论模型中依然适用。
Little公式的直观意义：
LsWs表明排队系统的队长等于一个顾客平均逗留时间内到达的顾客数。
LqWq表明排队系统的等待队长等于一个顾客平均等待时间内到达的顾客数。
1.2系统有多个服务台S>1情形
当系统中有s个服务台，系统服务能力为S，服务强度为-
s
系统中顾客的平均队长为:
系统中顾客的逗留时间为：
WsL（11）
系统中顾客的平均等待时间为：
1
WqWs—（12）
系统中顾客的平均等待队长为：
LqWq（13）
1.3LINGO中的相关函数及相关参数计算公式
（1）顾客等待概率的公式：
（14）
Pvait=@peb（load,S）
其中S为服务台服务台个数，load为系统到达的载荷，即load-o
（2）
T
Sload
（15）
顾客的平均等待时间公式:
其中T为顾客接受服务的平均时间，有
当load>s时无意义，表示当系统负荷超过服务台个数时，排队系统达到不稳定状态,队伍将越排越长。
一1
（3）
（16）
系统中顾客的平均逗留时间眼Wq—
（4）系统中顾客的的平均队长LsWs（17）
（5）系统中顾客的的平均等待队长LqWq（18）
例1某机关接待室只有1名对外接待人员，每天工作10小时，来访人员和接待时间都是随机的。
设来访人员按照Poisson流到达，到达速率为8人/小时，接待人员的服务速
率为9人/小时，接待时间服从负指数分布。
（1）计算来访人员的平均等待时间，等候的平均人数。
（2）若到达速率增大为20人/小时，每个接待人员的服务速率不变，为使来访问人员平均等待时间不超过半小时，最少应该配置几名接待人员。
解答：
（1）该问题属于M/M/1/g排队模型。
S=1,8,9需要计算来访人员的平均等待时间Wq,等候的平均人数Lq。
LINGO程序为：
model:
lp=8;
u=9;
T=1/u;
load=lp/u;
S=1;
Pwait=@PEB（load,S）;!
等待概率;
W_q=Pwait*T/（S-load）;!
平均等待时间;
L_q=lp*W_q;!
顾客的平均等待队长;end
计算结果：
来访人员的平均等待时间Wq0.89小时=53分钟，等候的平均人数Lq7.1人。
（2）该问题属于M/M/S/a排队模型。
求最小的s使，来访人员的平均等待时间wq0.5。
LINGO程序为：
model:
min=S;
lp=20;
u=9;!
服务率;
T=1/u;
load=lp/u;
Pwait=@PEB（load,S）;!
接待人员的等待概率;
W_q=Pwait*T/（S-load）;!
平均等待时间;
W_q<=0.5;
L_q=lp*W_q;!
顾客的平均等待队长;
TT=W_q*60;
@gin（S）;
end
计算结果为：
最少需要接待人员S=3人。
来访人员等待概率为0.55，排队等待平均时间为4.7分钟，队长平均长度为1.58人。
2.损失制模型M/M/S/S
M/M/S/S模型表示顾客到达人数服从Poisson分布，单位时间到达率为，服务台服
务时间服从负指数分布，单位时间服务平均人数为。
当S个服务台被占用后，顾客自动
离开，不再等待。
这里我们给出LINGO中的有关函数及相关参数的计算公式
（1）系统损失概率
19）
Plost=@pel（load,S）
其中S为服务台服务台个数，load为系统到达的载荷，即load-o
损失概率表示损失的顾客所占的比率。
（20）
（2）单位时间内进入系统的平均顾客数
e（1Rost）
（3）系统中顾客的平均队长（系统在单位时间内占用服务台的均值）
Ls
（4）系统中顾客的平均逗留时间（服务时间）
1
WsT（22）
（5）系统服务台的效率
Ls
s
在损失制排队模型中，顾客平均等待时间Wq0，平均等待队长Lq0，因为没有顾
客等待。
例2某单位电话交换台有一部300门内线电话的总机，已知上班时间有30%的内线分
机平均每30分钟要一次外线电话，70%的分机每隔70分钟时要一次外线电话。
又知从外单
位打来的电话的呼唤率平均30秒一次，设与外线的平均通话时间为3分钟，以上时间都服
从负指数分布。
如果要求外线电话接通率为95%以上，问电话交换台应设置多少外线？
Post5%
LINGO程序为:
model:
min=S;
lp=480;!
每小时平均到达电话数;
u=20;!
服务率;
load=lp/u;
Plost=@PEL（load,S）;!
损失率;
Plost<=0.05;
lpe=lp*（1-Plost）;
L_s=lpe/u;!
顾客的平均队长;
eta=L_s/S;!
系统服务台的效率;
@gin（S）;
end计算结果为：
最小的电话交换台为S=30。
电话损失率为Plost=0.04，实际进入系统的电话平均为e460.7，平均队长Ls23，
系统服务台的效率0.768。
3.混合制模型M/M/S/K
混合制模型M/M/S/K，表示顾客到达人数服从Poisson分布，单位时间到达率为，服务台服务时间服从负指数分布，单位时间服务平均人数为，系统有S个服务台，系统对顾客的容量为K。
当K个位置被顾客占用时，新到的顾客自动离去。
当系统中有空位置时，新到的顾客进入系统排队等候。
对混合制模型，LINGO没有相关函数计算参数。
需要自己编程计算。
（1）混合制模型基本公式：
设稳定状态下系统有i（i0,1,2丄,K）个顾客的概率为p（i0,1,2,l,K）。
p0表示系统空闲的概率。
因此：
K
pi1pi0,i1,2,L,K（24）
i0
设出0,1,2,L,K）表示系统中有i个顾客时的输入强度，i（i0,1,2,L,K）表示
系统中有i个顾客时的服务强度。
在稳定状态下，可建立平衡方程：
1p10p0
0p02p2（11）p1
（25）
i1pi1i1pi1（ii）pi（i1,2,L,K1）
K1pK1K.pK
对于混合制系统M
升级会员 