网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

基于指标比对串联排队系统平均排队时间的近似方法  PDF

  • 吴登磊
  • 赵宁
  • 刘文奇
昆明理工大学理学院, 昆明, 650500

中图分类号: O226

最近更新:2020-08-25

DOI:10.16356/j.1005-2615.2020.04.017

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
目录contents

摘要

串联排队系统是构成排队网络的基本结构,但是除了满足马尔可夫性或服务时间为常数的串联排队系统外,一般的串联排队系统的平均排队时间难以精确计算。为了刻画串联排队系统各个站之间的关联性,本文提出指标比的概念,基于指标比对系统的平均排队时间进行研究;通过分析指标比的数值特性,得到指标比的拟合表达式,进而对系统下游工作站的平均排队时间提出近似方法。数值实验结果显示,本文提出的近似方法对串联排队系统平均排队时间的估计效果较好。

排队论是用来研究随机服务系统的一门学科,通常用于分析现实生活具有可变性或不规则流量的服务系统的性能,例如,交[

1]、生[2]等各项资源共享的服务系统。研究排队系统的目的是正确设计和有效运行各个服务系统,使之发挥最佳效益。串联排队系统是排队网络的一种基本结构,各服务站之间相互关联,上游工作站的输出过程是下游工作站的输入过程。除了满足马尔可夫性或服务时间为常数的串联排队系统外,其他串联排队系统的平均排队时间的解析表达式很难求解。目前,关于服务时间服从一般分布的串联排队系统的研究均采用近似分析的方法。

对于G/G/1排队系统,Kingman[

3]提出了平均排队时间的上限。Gelenbe[4],Kimura[5], Kobayashi[6]基于到达时间和服务时间的分布提出了平均排队时间的二阶矩近似方法。Wu和Shen[7]对GI/G/1排队系统的平均排队时间提出三阶矩近似方法。虽然近似分析方法不同,但大多数近似公式中都包含M/M/1排队系统的平均排队时间的表达式。

在Jackson网络中,外部到达过程为泊松过程,服务时间服从指数分布,每个站可以看成是独立的M/M/1排队系统,其排队时间存在解析表达[

8]。当串联排队系统不满足马尔可夫性时,各个站的离去过程通常是非更新过程,这导致串联排队系统的性能很难用解析的方法来分析。串联排队系统平均排队时间的近似方法得到广泛研[9⁃15]。Whitt[9]在20世纪80年代提出了排队网络分析方法(Queueing network analysis, QNA),Dai[10]以及Harrison[11]提出了排队网络方法(Queueing network, QNET),Wu[12]利用固有比率的近线性的特征提出了串联排队系统的近似方法。Wu和Zhao[13]研究了串联排队系统各个站之间的依赖性。受文献[13]的启发,本文提出采用指标比(Performance ratio,PR)的方法分析串联排队系统上游工作站对下游工作站的影响。

1 问题描述

本文研究的串联排队系统如图1所示,系统由2个站串联而成,每站有一个服务器。作业的到达过程为泊松过程,平均速率为λ,到达时间间隔的平方变异系数(定义1)等于1(Ca2=1),作业依次由两个站完成服务。2个站的平均服务时间分别为S1S2(单位:min),第1站和第2站的服务时间的平方变异系数分别为CS12CS22,每个服务器前的缓冲区为无限大。系统中第i个站的利用率为ρi=λ/μi,其中μi=1/Si,利用率高的工作站为瓶颈站,第i个站的平均排队时间为TQii=1,2。本文所研究的系统假设第2站为瓶颈站,作业遵循先到先服务(First come first serve, FCFS)的服务规则,且CSi21i=1,2

图 1 两个站的串联排队系统

Fig.1 Tandem queuing systems with two stations

定义1   假设随机变量X的期望为EX,方差为DX,则随机变量X的平方变异系数为C2=DXE2X

图1所示系统中,由于第2站的缓冲区为无穷大,第1站的输出过程不受第2站的影响,所以第1站的排队时间可以通过单一服务器的排队系统分析。根据文献[

14],第1站的平均排队时间等于Ca2+CS122ρ11-ρ11μ1。如果串联排队系统第1站的服务时间满足S1=0,则系统中第1站对第2站没有任何影响,第2站的输入过程为第1站的到达过程,其平均排队时间等于Ca2+CS222ρ21-ρ21μ2。然而,如果第1站的服务时间S10,且服从一般分布,第2站的输入过程一般是非更新过程,则难以用解析的方法求解第2站的平均排队时间。

下面使用指标比的方法刻画第1站对第2站的影响,PR的定义如下。

定义2   串联排队系统第2站的PR

PR=TQ2Ca2+CS222ρ21-ρ21μ2 (1)

式中TQ2表示串联排队系统第2站的平均排队时间。

串联排队系统第2站的指标比为串联排队系统在S10情况下与S1=0的情况下第2站的平均排队时间之比。由(1)式知,TQ2=PR×Ca2+CS222ρ21-ρ21μ2。因此如果PR已知,则可以通过PR计算系统第2站的平均排队时间。下一节重点分析PR的值。

2 指标比的值

2.1 服务时间服从指数分布

假设系统中2个服务器的服务时间均服从指数分布(即CS12=CS22=1),外部到达过程为泊松过程,由文献[

15]可知,第1站的离去过程为泊松过程,第2站的输入过程与系统的到达过程独立同分布,因此,TQ2=Ca2+CS222ρ21-ρ21μ2。由式(1)得,当服务时间服从指数分布,第2站的PR

PR=1 (2)

2.2 服务时间为常数

当系统中2个服务器的服务时间都是常数(即CS12=CS22=0)时,由于系统的平均排队时间只受瓶颈站影[

16],故

TQ1+TQ2=Ca22ρ21-ρ21μ2 (3)

式中,TQ1=Ca22ρ11-ρ11μ1ρ1=ρ2μ2μ1。由式(3)TQ2=Ca22ρ21-ρ21μ2-Ca22ρ11-ρ11μ1。将TQ2μi=1/Si代入式(1),得服务时间为常数时第2站的PR

PR=TQ2Ca22ρ21-ρ21μ2=1-S12S221-ρ21-ρ2S1S2 (4)

S15,10,15,20,25,29ρ2{0.1,, 0.95}S2=30。当系统的服务时间为常数时,PR与系统各参数的关系如图2所示。

图2 PR与参数之间的关系

Fig.2 Relationship between PR and system parameters

图2可知:当系统的服务时间为常数时,PRρ2的增大而增大,随S1的增大而减小,PR的取值均小于1。

2.3 服务时间服从一般分布

当系统中2个服务器的服务时间服从一般分布,很难求解指标比的解析表达式。为了分析指标比与各个参数之间的关系,模拟了4 860种不同参数组合下的串联排队系统,这4 860个系统的服务时间均服从Gamma分布,参数取值如下:Ca2=1,S15,10,15,20,25,29S2=30CS12{0.1,0.2, ,0.9}CS220.1,0.2,,0.9ρ2{0.1,0.2,, 0.9,0.95}

对于每组给定的参数,模拟运行30个样本,每个样本取第400 001个至第600 000个作业的平均排队时间,平均排队时间模拟值的置信水平大于95%,保证了模拟数据的可靠性。根据模拟得到的平均排队时间TQ2计算PR的值,数值结果如图3所示。图3显示PRρ2的单调增函数,且PR<1PR的取值受如下系统参数的影响:Ca2,S1,CS12,S2,CS22,ρ2

图 3 不同参数下的指标比

Fig.3 Performance ratio with different parameters

为了便于分析,将服务时间服从一般分布的串联排队系统的第2站的指标比记为PR=f(Ca2,S1,CS12,S2,CS22,ρ2)。下面采用控制变量法,分别讨论如下4种情况下指标比与各参数的关系:f(1,S1,0.5,30,0.9,0.5)f(1,15,CS12,30,0.2,0.1)f(1,15,0.3,30,CS22,0.3)f(1,10,0.1,30,0.5,ρ2)。4种情况下指标比的模拟数值结果如图4所示,PR具有以下特性:(1)PRS1的增大而减小;(2)PR关于CS12近似线性递增,其中0CS121;(3)PRCS22的增大而增大,且具有一定的上凸性;(4)PRρ2的增大而增大,并且具有一定的下凸性。

图4 PR与系统参数的关系

Fig.4 Relationship between PR and system parameters

由于PR关于CS12近似线性递增(0CS121),且线段上任意一点可以由2个端点线性组合,因此如果0<CS12<1,其PR值可以表示为系统服务时间分别为常数(CS12=0)和指数分布(CS12=1)情况下PR的线性组合,即式(2)式(4)的线性组合。记PR的近似值为PR',则

PR'k1×1-S12S221-ρ21-ρ2S1S2+k2 (5)

根据模拟得到的4 860个PR数据进行拟合,得

k1=(1-CS12)S2S1(1+CS22)2S1S2,k2=1-(1-CS12)S2S1(1+CS22)2S1S2

因此,PR的近似表达式为

PR'1-S12S221-ρ21-ρ2S1S2(1-CS12)S2S1(1+CS22)2S1S2 (6)

注:(1)如果串联排队系统的服务时间服从指数分布或为常数,则PR'=PR。(2)虽然式(6)是根据服务时间服从Gamma分布的串联排队系统的数值结果拟合而成,但也适用于其他平方变异系数小于1的任意分布(见3.2节)。

3 误差分析

3.1 指标比的误差分析

第2.3节给出了指标比PR的拟合表达式,即式(6),下面通过模拟验证其有效性。令S1=20,CS12=0.9,S2=30,CS22=0.5ρ2{0.1,0.2,,0.9}表1分别比较了不同ρ2下第二个站的平均排队时间TQ2、指标比PR、拟合指标比PR',及PRPR'的相对误差。

表1 PR的相对误差
Table 1 Relative error of PR
ρ2TQ2PRPR'PR'-PRPR/%
0.1 2.5 0.981 0.992 1.076
0.2 5.5 0.982 0.992 1.090
0.3 9.5 0.981 0.993 1.201
0.4 14.7 0.983 0.993 1.051
0.5 22.1 0.982 0.994 1.191
0.6 33.2 0.984 0.995 1.103
0.7 51.6 0.984 0.995 1.172
0.8 88.8 0.987 0.996 0.947
0.9 200.5 0.990 0.998 0.774

表1所示,第2站的平均排队时间TQ2ρ2的增大而增大,拟合指标比的相对误差均小于2%,拟合效果较好。此外,针对2.3节中模拟得到的4 860个PR与拟合值PR'比较,发现相对误差小于5%的拟合指标比达到98.1%。

3.2 平均排队时间的误差分析

研究指标比的目的是求解第2站的平均排队时间。由式(1)和(6)得TQ2的近似值TQ2'如下:

TQ2'=1-S12S221-ρ21-ρ2S1S2(1-CS12)S2S1(1+CS22)2S1S2Ca2+CS222·ρ21-ρ2S2 (7)

下面通过模拟服务时间服从Gamma分布、均匀分布和三角分布的串联排队系统,验证TQ2'的准确性。

3.2.1 系统的服务时间服从Gamma分布

假设串联排队系统中2个站的服务时间均服从Gamma分布,选取34组不同参数的串联排队系统(见表2),分别比较不同参数下第2站的平均排队时间TQ2、近似值TQ2'及两者的相对误差,如表2所示。表2中的34组数据结果显示近似值TQ2'的相对误差均小于3%,运用指标比PR的拟合表达式可以较准确地计算平均排队时间TQ2

表2 TQ2的相对误差
Table 2 Relative error of TQ2
类别ρ2S1CS12S2CS22TQ2TQ2'TQ2'-TQ2TQ2/%

ρ2

0.10 20 0.5 30 0.5 2.255 2.279 1.07
0.20 20 0.5 30 0.5 5.090 5.150 1.17
0.30 20 0.5 30 0.5 8.726 8.871 1.66
0.40 20 0.5 30 0.5 13.606 13.877 1.99
0.50 20 0.5 30 0.5 20.494 20.956 2.25
0.60 20 0.5 30 0.5 30.917 31.691 2.50
0.70 20 0.5 30 0.5 48.566 49.797 2.54
0.80 20 0.5 30 0.5 84.368 86.470 2.49
0.90 20 0.5 30 0.5 194.792 197.867 1.58
0.95 20 0.5 30 0.5 416.534 422.165 1.35

S1

0.80 5 0.5 30 0.5 89.283 89.992 0.79
0.80 10 0.5 30 0.5 89.740 88.941 0.90
0.80 15 0.5 30 0.5 88.181 88.750 0.65
0.80 20 0.5 30 0.5 86.470 84.368 2.49
0.80 25 0.5 30 0.5 81.696 79.732 2.46
0.80 29 0.5 30 0.5 73.458 73.213 0.33

CS12

0.80 15 0.1 30 0.5 85.950 85.978 0.03
0.80 15 0.2 30 0.5 86.800 86.779 0.02
0.80 15 0.3 30 0.5 87.550 87.739 0.22
0.80 15 0.4 30 0.5 88.200 87.591 0.70
0.80 15 0.5 30 0.5 88.750 88.181 0.65
0.80 15 0.6 30 0.5 89.200 87.950 1.42
0.80 15 0.7 30 0.5 89.550 87.931 1.84
0.80 15 0.8 30 0.5 89.800 88.345 1.65
0.80 15 0.9 30 0.5 89.950 88.090 2.11

CS22

0.80 5 0.9 30 0.1 66.000 65.912 0.13
0.80 5 0.9 30 0.2 72.000 72.098 0.14
0.80 5 0.9 30 0.3 78.000 78.253 0.32
0.80 5 0.9 30 0.4 84.000 84.005 0.01
0.80 5 0.9 30 0.5 90.000 89.962 0.04
0.80 5 0.9 30 0.6 96.000 95.560 0.46
0.80 5 0.9 30 0.7 102.000 101.519 0.47
0.80 5 0.9 30 0.8 108.000 107.462 0.50
0.80 5 0.9 30 0.9 114.000 113.971 0.03

3.2.2 系统的服务时间服从均匀分布

假设串联排队系统中第1站和第2站的服务时间分别服从均匀分布U0.77,29.23U1.54,58.46,则S1=15,S2=30,CS12=0.3,CS22=0.3

ρ2=0.1,0.2,,0.9,0.95,对串联排队系统进行仿真模拟,从而得到第2站的平均排队时间的模拟值。将模拟值与式(7)计算得到的近似平均排队时间进行比较,数值结果如图5所示。图5中平均排队时间的模拟值与拟合值几乎重合,误差较小。

图 5 服务时间为均匀分布情况下的误差分析

Fig.5 Error analysis for case with uniform distribution

3.2.3 系统的服务时间服从三角分布

假设串联排队系统中第1站和第2站的服务时间分别服从参数为ai,bi,ci的三角分布(i=1,2),其概率密度为

f(xai,bi,ci)=2(x-ai)(bi-ai)(ci-ai)ai<xci2(bi-x)(bi-ai)(bi-ci)ci<xbi                0       

a1=1.83b1=38.17c1=5a2=3.67b2=76.33c2=10,则S1=15S2=30CS12=0.3CS22=0.3

ρ2=0.1,0.2,,0.9,0.95,对服务时间服从三角分布串联排队系统进行仿真模拟,得到第2站的平均排队时间的模拟值。图6对模拟值与式(7)计算得到的近似平均排队时间进行比较。图6显示通过仿真模拟得到的平均排队时间与拟合得到的平均排队时间误差较小,从而验证了通过指标比近似计算串联排队系统平均排队时间方法的可靠性。

图 6 服务时间为三角分布情况下的误差分析

Fig.6 Error analysis the case with triangular distribution

4 结 论

本文首次提出用指标比来刻画串联排队系统上游工作站对下游工作站的影响,并基于指标比提出串联排队系统平均排队时间的近似方法。通过分析指标比与系统各参数之间的关系,拟合得到指标比的近似表达式,利用指标比可以较准确地估计串联排队系统瓶颈站的平均排队时间。本文运用指标比的方法分析了2个站的串联排队系统,其中第2站为瓶颈站,未来可以运用该方法研究其他排队网络,例如:瓶颈站在前的串联排队系统、缓冲区有限的串联排队系统、复杂网络等。

参考文献

1

冯霞孟金双. 基于排队论的航班滑出时间预测[J]. 南京航空航天大学学报, 2016,48(5): 772-780. [百度学术

FENG Xia, MENG Jinshuang. Aircraft taxi-out time prediction based on queuing theory[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2016,48(5): 772-780. [百度学术

2

李蕊赵宁刘文奇, 等待时间有限的串行生产系统的缓冲区优化[J]. 控制与决策2018, 33(2): 345-350. [百度学术

LI Rui, ZHAO Ning, LIU Wenqi. Optimization of buffer for series production systems with waiting time constraints[J]. Control and Decision, 2018, 33(2): 345-350. [百度学术

3

KINGMAN J. On queues in heavy traffic[J]. Journal of the Royal Statistical Society, Series B (Methodological), 1962,24(2): 383-392. [百度学术

4

GELENBE E. On approximate computer system models[J]. Journal of the ACM, 1975, 22(2): 261-269. [百度学术

5

KIMURA T. Refining diffusion approximations for GI/G/1 queues: A tight discretization method[C]// Proceedings of Teletraffic Issues in an Advanced Information Society. Netherlands: Elsevier Science Publishers B. V., 1985: 317-323. [百度学术

6

KOBAYASHI H. Application of the diffusion approximation to queueing networks I: Equilibrium queue distributions[J]. Journal of the ACM, 1974, 21(2): 316-328. [百度学术

7

WU Kan, SHEN Yichi. Three-moment approximation for the mean queue time of a GI/G/1 queue[J]. IISE Transactions, 2018, 50(2): 63-73. [百度学术

8

JACKSON J R. Networks of waiting lines[J]. Operations Research, 1957, 5(4): 518-521. [百度学术

9

WHITT W. The queueing network analyzer[J]. Bell System Technical Journal, 1983, 62(9): 2779-2815. [百度学术

10

DAI J, HARRISON J M. Reflected Brownian motion in an orthant: Numerical methods for steady-state analysis[J]. The Annals of Applied Probability, 1992, 2(1): 65-86. [百度学术

11

HARRISON J M, Nguyen V. The QNET method for two-moment analysis of open queueing networks[J]. Queueing Systems, 1990, 6(1): 1-32. [百度学术

12

WU Kan, MCGINNIS L. Interpolation approximations for queues in series[J]. IIE Transactions, 2013,45(3): 273-290. [百度学术

13

WU Kan, ZHAO Ning. Dependence among single stations in series and its applications in productivity improvement[J]. European Journal of Operational Research, 2015, 247(1): 245-258. [百度学术

14

KINGMAN J. Some inequalities for the queue GI/G/1[J]. Biometrika, 1962, 49(3/4): 315-324. [百度学术

15

BURKE P J. The output of a queueing system[J]. Operations Research, 1965, 4(6): 699-704. [百度学术

16

FRIEDMAN H D. Reduction methods for tandem queueing systems[J]. Operations Research, 1965, 13(1): 121-131. [百度学术

您是第位访问者
网站版权 © 南京航空航天大学学报
技术支持:北京勤云科技发展有限公司
请使用 Firefox、Chrome、IE10、IE11、360极速模式、搜狗极速模式、QQ极速模式等浏览器,其他浏览器不建议使用!