数学建模社区-数学中国

标题: 排队论模型(二):生灭过程 、 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 y2 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 M8 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 b2.3 忙期和闲期  d! }2 `: N- b9 w1 c+ W; j

; T' U1 d( \: a' p4 N7 u0 c$ m( U0 Y' X9 y

- Z, ^; Y8 v) G个顾客在系统内的平均逗留时间应等于服务员平均连续忙的时间。
$ i+ D7 ]- ]$ t, a, t# \
3 V! Z. c: ?) X, ]6 J3 与排队论模型有关的 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) |) cp0=1-Pwait;$ o" L2 Y# o4 y. K. Y* t( E1 ?
Pt_gt_10=@exp(-1);
! l$ Q4 V+ K5 wend + 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' u4 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# z7 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: P3 K2 n# o; J( J, h4 T  h

: P  Y6 F+ ?& t; B7 e% I5 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% k6 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! lmodel: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 KW_s=L_s/lamda;
+ P4 [) y  ]: M1 R6 gend4 }! 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 a9 \) 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