数学建模之排队论模型

1第五讲排队论模型本讲主要内容1.排队论的基本概念2.单服务台的排队模型3.多服务台的排队模型4.排队系统的最优化问题5.数学建模实例：校园网的设计和调节收费问题5.1排队论的基本概念5.1.1什么是排队系统排队论也称随机服务系统理论，它是20世纪初由丹麦数学家Erlang应用数学方法在研究电话话务理论过程中而发展起来的一门学科，在实际中有广泛的应用。它涉及的是建立一些数学模型，藉以对随机发生的需求提供服务的系统预测其行为。现实世界中排队的现象比比皆是，如到商店购货、轮船进港、病人就诊、机器等待修理等等。排队的内容虽然不同，但有如下共同特征：（1）有请求服务的人或物，如候诊的病人、请求着陆的飞机等，我们将此称为“顾客”。（2）有为顾客提供服务的人或物，如医生、飞机跑道等，我们称此为“服务员”。由顾客和服务员就组成服务系统。（3）顾客随机地一个一个（或者一批一批）来到服务系统，每位顾客需要服务的时间不一定是确定的，服务过程的这种随机性造成某个阶段顾客排长队，而某些时候服务员又空闲无事。为了叙述一个给定的排队系统，必须规定系统的下列组成部分：1.输入过程即顾客来到服务台的概率分布。排队问题首先要根据原始资料，由顾客到达的规律、作出经验分布，然后按照统计学的方法（如卡方检验法）确定服从哪种理论分布，并估计它的参数值。我们主要讨论顾客来到服务台的概率分布服从泊松分布，且顾客的达到是相互独立的、平稳的输入过程。所谓“平稳”是指分布的期望值和方差参数都不受时间的影响。2.排队规则即顾客排队和等待的规则。排队规则一般有即时制和等待制两种。所谓即时制就是服务台被占用时顾客便随即离去；等待制就是服务台被占用时，顾客便排队等候服务。等待制服务的次序规则有先到先服务、随机服务、有优先权的先服务等，我们主要讨论先到先服务的系统。3.服务机构服务机构可以是没有服务员的，也可以是一个或多个服务员的；可以对单【修理工录用问题】工厂平均每天有一台机器发生故障而需要修理，机器的故障数服从泊松分布。修理一台机器平均花费20元。现有技术水平不同的修理工人A和B，A种修理工平均每天能修理1.2台机器，每天工资3元；B种修理工平均每天能修理1.5台机器，每天工资5元，两种修理工修理机器的时间为负指数分布。问工厂录用哪种工人较合算？PDF文件使用pdfFactoryPro试用版本创建独顾客进行服务，也可以对成批顾客进行服务。和输入过程一样，多数的服务时间都是随机的，且我们总是假定服务时间的分布是平稳的。若以ξn表示服务员为第n个顾客提供服务所需的时间，则服务时间所构成的序列{ξn}，n=1，2，…所服从的概率分布表达了排队系统的服务机制，一般假定，相继的服务时间ξ1，ξ2，…是独立同分布的，并且任意两个顾客到来的时间间隔序列{Tn}也是独立的。如果按服务系统的以上三个特征的各种可能情形来对服务系统进行分类，那么分类就太多了。因此，现在已被广泛采用的是按顾客相继到达时间间隔的分布、服务时间的分布和服务台的个数进行分类。排队论主要是对服务系统建立数学模型，研究如下内容：（1）排队系统的概率分布问题，主要是研究队长分布、等待时间分布和忙期分布等；（2）最优化问题：分为静态最优化和动态最优化，即为系统的最优设计和系统的最优运行问题；（3）排队系统的统计推断：判断一个给定的排队系统符合哪种模型，以便于根据排队理论进行分析研究。5.1.2排队模型的标准形式排队模型的标准形式为X/Y/Z/A/B/C，其中：X表示顾客来到时间间隔的分布类型；Y表示服务时间的分布类型；Z表示服务员个数；A系统容量;B顾客源个数;C服务规则.例如先来先服务的等待排队模型主要由三参数法即X/Y/Z，“M/M/1/k/∞/FCFS”表示顾客到达间隔时间和服务时间均服从负指数分布，一个服务台，系统至多容纳k个顾客潜在的顾客数不限，先来先服务的排队系统。“M/M/c”即Poisson输入，负指数服务时间分布，c个服务台的等待制排队模型。“M/G/1”即Poisson输入，一般服务时间分布，单个服务台的等待制排队模型。5.1.3排队系统的运行指标研究排队问题的目的，是研究排队系统的运行效率，估计服务质量，确定系统参数的最优值，以决定系统的结构是否合理，设计改进措施等。所以，必须确定用来判断系统运行优劣的基本数量指标，这些数量指标通常是：（1）队长指排队系统中的顾客数，它的期望值记为Ls；排队长，指在排队系统中排队等待服务的顾客数，其期望值记为Lq。系统中的顾客数=等待服务的顾客数+正被服务的顾客数所以Lq（或Ls）越大，说明服务效率越低。（2）逗留时间指一个顾客在排队系统中的停留时间，即顾客从进入服务系统到服务完毕的整个时间。其期望值记为Ws。等待时间，指一个顾客在排队系统中等待服务的时间，其期望值记为Wq。逗留时间=等待时间+服务时间（3）忙期指从顾客到达空闲服务机构起到服务机构再次为空闲这段时间长度，即服务机构连续工作的时间长度。它关系到服务员的工作长度，即服务机构连续工作的时间长度。PDF文件使用pdfFactoryPro试用版本创建它关系到服务员的工作强度、忙期的长度和一个忙期中平均完成服务的顾客数，这些都是衡量服务效率的指标。要计算以上这些指标必须知道系统状态的概率，所谓系统状态即时刻t时排队系统中的顾客数。如果时刻t时排队系统中有n个顾客，就说系统的状态是n，其概率一般用Pn(t)表示。求Pn(t)的方法，首先要建立含Pn(t)的关系式，因t为连续变量而n只取非负整数，所以建立的Pn(t)的关系式一般是微分差分方程，这时要求方程的解是不容易的，有时即使求出也很难利用。因此，往往只求稳态解Pn，求Pn并不一定求t→∞时的Pn(t)极限，而只需由)('tPn=0，用Pn代替Pn(t)即可。5.2单服务台的排队模型设系统的输入过程服从泊松分布，服务时间服从负指数分布，单服务台的排队系统有以下三种情形：（1）标准型：M/M/1(M/M/1/∞/∞);（2）系统容量有限制：M/M/1/N/∞;（3）顾客源为有限的：M/M/1/∞/m.5.2.1标准型：M/M/1M/M/1模型是指顾客源为无限，顾客到达相互独立，到达过程是平稳的，到达率数服从参数为l的泊松分布：单服务台、队长无限、先到先服务；各顾客的服务时间服从参数为μ的负指数分布，且相互独立。首先求出排队系统在任意时刻t的、状态为n的概率Pn(t)，已知顾客到达率服从参数为λ的泊松分布，服务时间服从参数为μ的负指数分布，由此决定了[t，t+△t]时间间隔内：（1）有1个顾客到达的概率为λ△t+o(△t)，没有顾客到达的概率是1-λ△t+o(△t)。（2）当有顾客在接受服务时，1个顾客被服务完了的概率是μ△t+o(△t)，没有服务完的概率是1-μ△t+o(△t)。（3）多于一个顾客到达或服务完的概率为o(△t)，均可忽略。注1：因为单位时间内顾客到达数X～P（λ），所以Δt时间间隔内顾客到达数Y～P（λΔt），因而在Δt时间间隔内有1个顾客到达的概率为：P{Y=1}=λΔt·e-λΔt=λΔt+o(Δt)，没有顾客到达的概率为P{Y=0}=e-λΔt=1-λΔt+o(Δt)。注2：由于服务时间T～E（μ），故在有顾客接受服务时，1个顾客被服务完的概率为P{T≤Δt}=1-e-μΔt=μΔt+o(Δt)，没有被服务完的概率为1-μΔt+o(Δt)。在t+△t时刻，系统中有n个顾客的状态由t时刻的以下状态转化而来：①t时刻系统中有n个顾客，没有顾客到达且没有顾客服务完毕，其概率为：[1-λ△t+o(△t)][1-μ△t+o(△t)]=(1-λ△t-μ△t)+o(△t)；②t时刻系统中有n+1个顾客，没有顾客到达且有1个顾客服务完毕，其概率为：[1-λ△t+o(△t)][μ△t+o(△t)]=μ△t+o(△t)；③t时刻系统中有n-1个顾客，有1个顾客到达且没有顾客服务完毕，其概率为：[λ△t+o(△t)][1-μ△t+o(△t)]=λ△t+o(△t)；④其他状态的概率为o(△t)。因此，在t+△t时刻，系统中有n个顾客的概率Pn(t+△t)满足：)()()()1)(()(11tttPttPtttPttPnnnnΔ+Δ+Δ+Δ-Δ-=Δ+-+olmml.移项整理，两边同除以△t，得PDF文件使用pdfFactoryPro试用版本创建ΔΔ++-+=Δ-Δ++-)()()()()()()(11omlml.令△t→0，得L2,1)()()()()(11=+-+=+-ntPtPtPdttdPnnnnmlml.当n=0时，因为)()1)(()1)(()(100ttttPttPttPΔ+ΔΔ-+Δ-=Δ+omll所以有)()()(100tPtPdttdPml+-=.对于稳态情形，与t无关，其导数为零。因此，得到⎩⎨⎧=+-=+-++-01,0)(1011PPnPPPnnnmlmlml这是关于Pn的差分方程，也反映出了系统状态的转移关系，即每一状态都是平衡的，求解得()01/PPml=，递推可得())1(/0≥=nPPnnml.由概率的性质知∑∞==01nnP，将上式代入λ/μ1时可得到ml/10-=P()()nnPmlml//1-=.因为顾客到达规律服从参数为λ的泊松分布，服务时间服从参数为μ的负指数分布，其期望值就分别为λ，1/μ。所以λ表示单位时间内平均到达的顾客数，μ表示单位时间内能服务完的顾客数。如果令ρ=λ/μ，这时ρ就表示相同时间内顾客到达的平均数与能被服务的平均数之比，它是刻画服务效率和服务机构利用程度的重要标志，称ρ为服务强度。上面在ρ1的条件下得到了稳定状态下的概率Pn，n=0，1，2，…。其实，如果ρ1，可以证明排队长度将是无限增加的，即使ρ=1的情况下，P0(t)也是随时间而变化的，系统达不到稳定状态.因此，这里只讨论ρ1时情况，从上面的推导知Pn=(1-ρ)ρnn=0,1,2,….下面计算出系统的运行指标.(1)队长（平均顾客数）：由于系统的状态为n时即系统中有n个顾客，由期望的定义得)1/()1(10rrrr-=-==∑∑∞=∞=nnnnsnnpL)./(lml-=(2)排队长：（等待的平均顾客数）PDF文件使用pdfFactoryPro试用版本创建∑∑∞=∞=--=-=11)1()1()1(nnnnqnpnLrr)1/(2rr-=)./(lmrl-=可以证明，顾客在系统中逗留时间服从参数为μ-l的负指数分布。因此，有(3)系统中顾客的平均逗留时间：)/(1lm-=sW.(4)系统中顾客的平均等待时间：)/(1lmrm-=-=sqWW由以上结论可以看出，各指标之间有如下关系：ssWLl=；；qqWLl=；；mlm/,/1+=+=qsqsLLWW，这些关系式称为Little公式.在指标的计算过程中，一般只要计算其中一个，其它的指标便可随之导出。5.2.2系统容量有限制：M/M/1/N/∞因为是单服务台，排队系统的容量为N，即是排队等待的顾客最多为N-1，在某时刻一顾客到达时，如系统中已有N个顾客，那么这个顾客就被拒绝进入系统。假设顾客平均到达率为l，平均服务率为μ，在研究系统中有n个顾客的概率Pn(t)时，和标准型M/M/1研究方法相同，当n=N时有)()()(1'tPtPtPNNNml-=-在稳态情形下，令mlr=，得⎜⎜⎜⎝⎛=-=+=+=--+111011,,2,1,)1(NNnnnPPNnPPPPPrrrrL在条件∑==NiiP01下解上式得到⎪⎪⎩⎪⎪⎨⎧≤≤--=≠--=++.1111,11110NnPPnNnNrrrrrr注：PDF文件使用pdfFactoryPro试用版本创建（1）如果ρ=1（即l=μ），由1110+====NPPPNL，即到达率和服务率相等，在稳态情况下系统不会出现排队等待现象；（2）这里因为