数学建模社区-数学中国
标题:
排队论模型(二):生灭过程 、 M / M /s 等待制排队模型、多服务台模型
[打印本页]
作者:
浅夏110
时间:
2020-6-12 09:59
标题:
排队论模型(二):生灭过程 、 M / M /s 等待制排队模型、多服务台模型
1 生灭过程
9 Z5 Y/ |) D1 B
一类非常重要且广泛存在的排队系统是生灭过程排队系统。生灭过程是一类特殊的随机过程,在生物学、物理学、运筹学中有广泛的应用。在排队论中,如果 N(t) 表示 时刻t 系统中的顾客数,则{N(t),t ≥ 0}就构成了一个随机过程。如果用“生”表示顾 客的到达,“灭”表示顾客的离去,则对许多排队过程来说,{N(t),t ≥ 0}就是一类特殊的随机过程-生灭过程。
6 z/ v2 `6 U' J
2 ~: l; {" x$ {) F4 V
下面结合排队论的术语给出生灭过程的定义。
5 P2 ?- O1 D, }. F# W: X, a
3 s6 o3 }' H0 ^: t% D4 N
# {) T5 U% M" T2 H& d! I; v2 N- M
; ]! o/ P) G `2 K
为求平稳分布,考虑系统可能处的任一状态 n 。假设记录了一段时间内系统进入状 态n 和离开状态 n 的次数,则因为“进入”和“离开”是交替发生的,所以这两个数要么相等,要么相差为 1。但就这两种事件的平均发生率来说,可以认为是相等的。即当 系统运行相当时间而到达平衡状态后,对任一状态 n 来说,单位时间内进入该状态的平 均次数和单位时间内离开该状态的平均次数应该相等,这就是系统在统计平衡下的“流 入=流出”原理。根据这一原理,可得到任一状态下的平衡方程如下:
% f/ w$ h8 [$ K$ n" M3 y
2 Q+ u3 t/ d9 k# E
( }; }( a! j2 S( F$ x
+ n3 S- O( w3 c3 P4 m9 t2 Y3 O' r
) r. I/ m) z( F I# s7 [. u
* Q1 ?" A* Q X( ]
述公式得到平稳状态的概率分布。
* t7 A$ F+ g# H$ m9 }0 f' G0 M
8 p% v: l+ l# n, y. f! i
2 M / M /s 等待制排队模型
8 z2 V* A, e1 v' b: G0 z9 U
2.1 单服务台模型
% |7 G! a- _5 N" b5 w' D
单服务台等待制模型 M / M /1/ ∞ 是指:顾客的相继到达时间服从参数为λ 的负指 数分布,服务台个数为 1,服务时间V 服从参数为 μ 的负指数分布,系统空间无限, 允许无限排队,这是一类最简单的排队系统。
8 Y2 C1 U% P' | d- I# {( X# q
8 R* `0 g9 e. x/ {
2.1 队长的分布
{* k; x% X: K" w8 o) T4 G
+ h4 {- v7 l( P9 d% X& \1 v
4 h h" ^! v' F6 e" ?
+ Y: Y, a* p/ w: z8 `. f- r
2.2 几个主要数量指标
. ~2 h" k% x) ^
对单服务台等待制排队系统,由已得到的平稳状态下队长的分布,可以得到平均队 长
3 X: ^: ~6 F5 v( m- W d
9 r0 b1 r& {; L" i E; i3 H8 F
9 p* A9 |2 E: i- g \
8 C& o. l( n8 }4 ~
# }9 J9 j# p$ ?' W5 ?8 C& P! [
5 {" `; c3 o# d' U
式(14)和式(15)通常称为 Little 公式,是排队论中一个非常重要的公式。
% i$ k2 v* M. x
. R' N$ d( _' R! M5 p0 K7 h3 b
2.3 忙期和闲期
d! }2 `: N- b9 w1 c+ W; j
; T' U1 d( \: a' p4 N
7 u0 c$ m( U0 Y' X9 y
- Z, ^; Y8 v) G
个顾客在系统内的平均逗留时间应等于服务员平均连续忙的时间。
$ i+ D7 ]- ]$ t, a, t# \
3 V! Z. c: ?) X, ]6 J
3 与排队论模型有关的 LINGO 函数
1 m) k: `& T8 ~8 B( d: X. o, l
(1)@peb(load,S) 该函数的返回值是当到达负荷为 load,服务系统中有 S 个服务台且允许排队时系 统繁忙的概率,也就是顾客等待的概率。
' m5 c5 X. p$ v
) p, A2 t8 ]% k" X
(2)@pel(load,S) 该函数的返回值是当到达负荷为 load,服务系统中有 S 个服务台且不允许排队时 系统损失概率,也就是顾客得不到服务离开的概率。
- G* O; m) c6 W5 B7 w. v
' X$ A3 Z7 u9 H! H
(3)@pfs(load,S,K) 该函数的返回值是当到达负荷为 load,顾客数为 K,平行服务台数量为 S 时,有限 源的 Poisson 服务系统等待或返修顾客数的期望值。
6 L$ _$ O3 G0 e! l2 M2 ^* M. U
0 u- ?7 I6 Q- S) p
例 1 某修理店只有一个修理工,来修理的顾客到达过程为 Poisson 流,平均 4 人 /h;修理时间服从负指数分布,平均需要 6min。试求:(1)修理店空闲的概率;(2) 店内恰有 3 个顾客的概率;(3)店内至少有 1 个顾客的概率;(4)在店内的平均顾客数; (5)每位顾客在店内的平均逗留时间;(6)等待服务的平均顾客数;(7)每位顾客平 均等待服务时间;(8)顾客在店内等待时间超过 10min 的概率。
! g* A% m* ^6 v' Q9 M" n
! K0 x+ w" e+ d6 b
7 x8 V" c( S) w2 M' L% i& q/ b
) C* W9 [: ~( x7 ?/ Z
编写 LINGO 程序如下:
1 [% a- D8 n6 @ u9 i/ L( `) x- Z+ g
2 z- g: I% J5 a7 V
model:
4 W( V6 C3 a2 b
s=1;lamda=4;mu=10;rho=lamda/mu;
( j6 |; v! J; w6 b
Pwait=@peb(rho,s);
8 n B# y& Y) |) c
p0=1-Pwait;
$ o" L2 Y# o4 y. K. Y* t( E1 ?
Pt_gt_10=@exp(-1);
! l$ Q4 V+ K5 w
end
+ D0 `( G/ p1 @& G( Y
4 多服务台模型( M / M /s/ ∞ )
# r& m/ ]+ i7 `( _4 V
设顾客单个到达,相继到达时间间隔服从参数为λ 的负指数分布,系统中共有 s 个 服务台,每个服务台的服务时间相互独立,且服从参数为 μ 的负指数分布。当顾客到 达时,若有空闲的服务台则马上接受服务,否则便排成一个队列等待,等待时间为无限。
) u2 X# s9 o; U7 ?" Z* x5 u
V2 n% C7 m; c9 R8 f- e2 j& M
L9 V( k: t' u
4 n: T; y. [! K9 Z+ V3 H% ]! I7 ]
公式(19)和式(20)给出了在平衡条件下系统中顾客数为 n 的概率,当 n ≥ s 时,即 系统中顾客数大于或等于服务台个数,这时再来的顾客必须等待,因此记
+ B H' A' ~6 K: d
- B4 Q" E& m4 p. D
! _% u: x3 F6 K+ Z# z
7 u- I* s+ z, _' _
式(21)称为 Erlang 等待公式,它给出了顾客到达系统时需要等待的概率。 对多服务台等待制排队系统,由已得到的平稳分布可得平均排队长 Lq 为:
2 D% U6 @6 S) G9 \
0 C# P, z' v; j- \$ F
" k- r) t, f* I1 o: P
3 K2 n# o; J( J, h4 T h
: P Y6 F+ ?& t; B7 e% I
5 r( n6 h5 h( F
对多服务台系统,Little 公式依然成立,即有
. I6 E" l0 o/ `7 O% ]5 n5 \
7 B% D ?1 L" q/ G' | W: Q9 j
0 m- h, w: g6 u% Q) u9 y, F, R( ^' X! Y
( c8 O$ C* r; Y5 S* S2 y% _- ^! n: B
例 2 某售票处有 3 个窗口,顾客的到达为 Poisson 流,平均到达率为 λ = 0.9人/ min ;服务(售票)时间服从负指数分布,平均服务率 μ = 0.4人/ min 。 现设顾客到达后排成一个队列,依次向空闲的窗口购票,这一排队系统可看成是一个
& {& H# g6 y' o3 N6 `
" _4 ^$ a' v7 }
M / M / s/ ∞ 系统,其中
/ D2 k+ ~" q/ X% k
6 r+ v" M; u$ f' j p
) O) Q; w" p4 X5 N p2 Z
9 w/ O9 P3 _4 t. H6 x& r
- z0 U4 ]$ D0 ~) x T& n. [' h1 L
* `4 I0 l) p! q# O r9 v4 G. b& Q
求解的 LINGO 程序如下:
6 J/ |% S4 j7 ?' T9 E
& e$ M% v0 T5 D* y! l
model:
3 p: O$ L+ _: {7 A% O) E
s=3;lamda=0.9;mu=0.4;rho=lamda/mu;rho_s=rho/s;
x# a) f- b2 a( R$ f
P_wait=@peb(rho,s);
) ^* S* I2 o* D$ Z
p0=6*(1-rho_s)/rho^3*P_wait;
& ~9 t9 x3 T0 a( G; g; O
L_q=P_wait*rho_s/(1-rho_s);
: e3 ]* u9 p h6 f6 s- y
L_s=L_q+rho;
% R% o- j6 h' O8 ]
W_q=L_q/lamda;
- J7 M6 c- C8 K
W_s=L_s/lamda;
+ P4 [) y ]: M1 R6 g
end
4 }! z9 Q' r+ O
l) b+ @" ^9 c$ D
————————————————
" i' a M0 J! Y5 @: _7 c/ e4 y- \" r
版权声明:本文为CSDN博主「wamg潇潇」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
; W3 Y. R" v% Y. t" D- S
原文链接:https://blog.csdn.net/qq_29831163/java/article/details/89735349
4 k; J: R8 d" L: t" h; F5 a
9 \) t: {9 ^, a, \( g. n
$ y d4 F5 \- {# ?2 Z7 ~, r2 C. |
作者:
SHINee0525
时间:
2020-8-17 14:46
感觉很好 很实用
% r- N6 Q3 B8 B4 u) V
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5