【管理运筹学】第 10 章 | 排队论(3,标准的 M/M/1 排队系统)
对标准的 M/M/1 排队系统进行了介绍,包括模型特征及分析以及各指标的计算公式。
引言
前两篇文章,分别对基本的排队论概念如 Kendall 记号和三个常见的分布如普阿松分布、负指数分布等作了介绍。有了这些基础后,我们便可以排队系统进行分析,我们首先讨论标准的 M / M / 1 M/M/1 M/M/1 模型即( M / M / 1 / ∞ / ∞ M/M/1/\infty/\infty M/M/1/∞/∞ )。
一、模型特征及分析
标准的 M / M / 1 M/M/1 M/M/1 模型是指下列条件下的排队系统:
- 输入过程 —— 顾客源无限,到达过程普阿松分布。
- 排队规则 —— 一队,且队长无限制,先到先服务。
- 服务机构 —— 一个服务台,各个顾客的服务时间是相互独立的,且服从相同的负指数分布。
此外,假定到达间隔时间和服务时间是相互独立的。
在分析标准的 M / M / 1 M/M/1 M/M/1 模型时,首先要求出系统在任意时刻 t t t 的状态为 n n n 的概率 P n ( t ) P_n(t) Pn(t) 其表示长为 t t t 的时间内,系统内有 n n n 个顾客的概率,它决定了系统运行的特征。
已知顾客到达时间服从参数为 λ \lambda λ 的普阿松过程,服务时间服从参数为 λ \lambda λ 的负指数分布,所以在区间 [ t , t + Δ t ] [t,t+\Delta t] [t,t+Δt] 内有以下三种情况:
- 有一个顾客的概率为 λ Δ t + ο ( Δ t ) \lambda\Delta t+\omicron(\Delta t) λΔt+ο(Δt) ;没有顾客到达的概率为 1 − λ Δ t + ο ( Δ t ) 1-\lambda\Delta t+\omicron(\Delta t) 1−λΔt+ο(Δt) 。
- 当有顾客在接受服务时,1 个顾客被服务完离去的概率是 μ Δ t + ο ( Δ t ) \mu\Delta t+\omicron(\Delta t) μΔt+ο(Δt) ,没有离去(仍在服务)的概率是 1 − μ Δ t + ο ( Δ t ) 1-\mu\Delta t+\omicron(\Delta t) 1−μΔt+ο(Δt) 。
- 多于一个到达或离去的概率为 ο ( Δ t ) \omicron(\Delta t) ο(Δt) ,可以忽略。
在时刻
[
t
,
t
+
Δ
t
]
[t,t+\Delta t]
[t,t+Δt] ,系统中有
n
n
n 个顾客(
n
>
0
n>0
n>0),存在如下表所示的四种可能情况(到达或离去两个以上的不计,√ 表示发生一个,× 表示没发生)。
以上四种情况是互不相容的,所以 P n ( t + Δ t ) P_n(t+\Delta t) Pn(t+Δt) 是各情况概率之和,只研究稳态的情形,可以得到: { λ P 0 = μ P 1 λ P n − 1 + μ P n + 1 = ( λ + μ ) P n , n ≥ 1 (1) \begin{cases} \lambda P_0=\mu P_1 \\ \lambda P_{n-1}+\mu P_{n+1}=(\lambda+\mu)P_n,n\geq1 \end{cases}\tag{1} {λP0=μP1λPn−1+μPn+1=(λ+μ)Pn,n≥1(1) 其中 λ , μ \lambda,\mu λ,μ 分别为到达率和服务率(单位时间内到达或服务的顾客数)。式 (1) 是关于 P n P_n Pn 的差分方程。它表明了各状态之间的转移关系,可以用下图(来源:知乎-运筹 OR 帷幄)表示。
在标准的
M
/
M
/
1
M/M/1
M/M/1 模型中,到达率和服务率是不随时间变化的,故上图中的
λ
1
,
μ
1
,
⋯
\lambda_1,\mu_1,\cdots
λ1,μ1,⋯ 均为对应的常数
λ
,
μ
\lambda,\mu
λ,μ 。
由式 (1) 的第一个式子,可知 P 1 = ( λ / μ ) P 0 P_1=(\lambda/\mu)P_0 P1=(λ/μ)P0 ,将其代入第二个式子,可以得到: μ P 2 = ( λ + μ ) ( λ / μ ) P 0 − λ P 0 ⟹ P 2 = ( λ / μ ) 2 P 0 \mu P_2=(\lambda+\mu)(\lambda/\mu)P_0-\lambda P_0\Longrightarrow P_2=(\lambda/\mu)^2P_0 μP2=(λ+μ)(λ/μ)P0−λP0⟹P2=(λ/μ)2P0 同理,可推得 P n = ( λ / μ ) n P 0 P_n=(\lambda/\mu)^nP_0 Pn=(λ/μ)nP0 。设 ρ = λ / μ < 1 \rho=\lambda/\mu<1 ρ=λ/μ<1 ,保证队列不会一直排下去,由概率的性质有 ∑ n = 1 ∞ P n = 1 (2) \sum_{n=1}^\infty P_n=1 \tag{2} n=1∑∞Pn=1(2) 将得到的 P n P_n Pn 的关系代入,有 ∑ n = 0 ∞ P n = P 0 ∑ n = 0 ∞ ρ n = P 0 1 1 − ρ = 1 \sum_{n=0}^\infty P_n=P_0\sum_{n=0}^\infty \rho^n=P_0\frac{1}{1-\rho}=1 n=0∑∞Pn=P0n=0∑∞ρn=P01−ρ1=1 于是我们可以得到 { P 0 = 1 − ρ P n = ( 1 − ρ ) ρ n ( n ≥ 1 , ρ < 1 ) (3) \begin{cases} P_0=1-\rho \\ P_n=(1-\rho)\rho^n \end{cases}(n\geq1,\rho<1) \tag{3} {P0=1−ρPn=(1−ρ)ρn(n≥1,ρ<1)(3) 这就是系统状态为 n n n 的概率。
上述构造的 ρ \rho ρ 具有实际意义,其根据表达式不同,可以有不同的解释。当表示为 ρ = λ / μ \rho=\lambda/\mu ρ=λ/μ 时,它是平均到达率与平均服务率之比;若表示为 ( 1 / μ ) / ( 1 / λ ) (1/\mu)/(1/\lambda) (1/μ)/(1/λ) ,它是一个顾客的服务时间与到达间隔时间之比,称其为服务强度(Traffic intensity),或称为话务强度; ρ \rho ρ 还可写为 1 − P 0 1-P_0 1−P0 ,表示服务机构的繁忙程度(1 减去顾客数为 0 的概率即有顾客的概率),所以又称为服务机构的利用率。
二、系统指标
利用式 (3) ,我们可以算出排队系统的一些指标。
1. 在系统中的平均顾客数(队长的期望)
用 L s L_s Ls 表示排队系统中的平均顾客数,由概率论的知识,我们有 L s = ∑ n = 0 ∞ n P n = ∑ n = 0 ∞ n ⋅ ( 1 − ρ ) ρ n = ρ 1 − ρ = λ μ − λ . L_s=\sum_{n=0}^\infty nP_n=\sum_{n=0}^\infty n\cdot(1-\rho)\rho^n=\frac{\rho}{1-\rho}=\frac{\lambda}{\mu-\lambda}. Ls=n=0∑∞nPn=n=0∑∞n⋅(1−ρ)ρn=1−ρρ=μ−λλ.
2. 在队列中的平均顾客数(队列长的期望)
用 L q L_q Lq 表示队列中等待的平均顾客数(队列长的期望)。 L q = ∑ n = 1 ∞ ( n − 1 ) P n = ρ 2 1 − ρ = L s − ρ = ρ λ μ − λ . L_q=\sum_{n=1}^\infty(n-1)P_n=\frac{\rho^2}{1-\rho}=L_s-\rho=\frac{\rho\lambda}{\mu-\lambda}. Lq=n=1∑∞(n−1)Pn=1−ρρ2=Ls−ρ=μ−λρλ.
3. 在系统中顾客平均逗留时间
用 W s W_s Ws 表示顾客在系统中平均逗留时间。在 M / M / 1 M/M/1 M/M/1 情形下,顾客逗留时间服从参数为 μ − λ \mu-\lambda μ−λ 的负指数分布,于是 W s = 1 μ − λ . W_s=\frac{1}{\mu-\lambda}. Ws=μ−λ1.
4. 在队列中顾客的平均等待时间
用 W q W_q Wq 表示在队列中顾客的平均等待时间,其值为平均逗留时间减去一个平均服务时间,即 W q = W s − 1 μ = ρ μ − λ . W_q=W_s-\frac{1}{\mu}=\frac{\rho}{\mu-\lambda}. Wq=Ws−μ1=μ−λρ.
以上各式归纳如下: ( 1 ) L s = λ μ − λ , ( 2 ) L q = ρ λ μ − λ ; ( 3 ) W s = 1 μ − λ , ( 4 ) W q = ρ μ − λ . (1)L_s=\frac{\lambda}{\mu-\lambda},(2)L_q=\frac{\rho\lambda}{\mu-\lambda};(3)W_s=\frac{1}{\mu-\lambda},(4)W_q=\frac{\rho}{\mu-\lambda}. (1)Ls=μ−λλ,(2)Lq=μ−λρλ;(3)Ws=μ−λ1,(4)Wq=μ−λρ. 它们的相互关系如下: ( 1 ) L s = λ W s , ( 2 ) L q = λ W q ; ( 3 ) W s = W q + 1 μ , ( 4 ) L s = L q + λ μ . (1)L_s=\lambda W_s,(2)L_q=\lambda W_q;(3)W_s=W_q+\frac{1}{\mu},(4)L_s=L_q+\frac{\lambda}{\mu}. (1)Ls=λWs,(2)Lq=λWq;(3)Ws=Wq+μ1,(4)Ls=Lq+μλ. 上式称为 Little 公式。
写在最后
标准的 M / M / 1 M/M/1 M/M/1 模型还是相对简单些的,用一些基础的概率论和高数即可对各结论进行数学推导。
开放原子开发者工作坊旨在鼓励更多人参与开源活动,与志同道合的开发者们相互交流开发经验、分享开发心得、获取前沿技术趋势。工作坊有多种形式的开发者活动,如meetup、训练营等,主打技术交流,干货满满,真诚地邀请各位开发者共同参与!
更多推荐
所有评论(0)