摘要
串联排队系统是构成排队网络的基本结构,但是除了满足马尔可夫性或服务时间为常数的串联排队系统外,一般的串联排队系统的平均排队时间难以精确计算。为了刻画串联排队系统各个站之间的关联性,本文提出指标比的概念,基于指标比对系统的平均排队时间进行研究;通过分析指标比的数值特性,得到指标比的拟合表达式,进而对系统下游工作站的平均排队时间提出近似方法。数值实验结果显示,本文提出的近似方法对串联排队系统平均排队时间的估计效果较好。
排队论是用来研究随机服务系统的一门学科,通常用于分析现实生活具有可变性或不规则流量的服务系统的性能,例如,交通[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,)的方法分析串联排队系统上游工作站对下游工作站的影响。
本文研究的串联排队系统如图1所示,系统由2个站串联而成,每站有一个服务器。作业的到达过程为泊松过程,平均速率为,到达时间间隔的平方变异系数(定义1)等于1(),作业依次由两个站完成服务。2个站的平均服务时间分别为和(单位:min),第1站和第2站的服务时间的平方变异系数分别为和,每个服务器前的缓冲区为无限大。系统中第个站的利用率为,其中,利用率高的工作站为瓶颈站,第个站的平均排队时间为,。本文所研究的系统假设第2站为瓶颈站,作业遵循先到先服务(First come first serve, FCFS)的服务规则,且,。
Fig.1 Tandem queuing systems with two stations
定义1 假设随机变量的期望为,方差为,则随机变量的平方变异系数为。
在图1所示系统中,由于第2站的缓冲区为无穷大,第1站的输出过程不受第2站的影响,所以第1站的排队时间可以通过单一服务器的排队系统分析。根据文献[14],第1站的平均排队时间等于。如果串联排队系统第1站的服务时间满足,则系统中第1站对第2站没有任何影响,第2站的输入过程为第1站的到达过程,其平均排队时间等于。然而,如果第1站的服务时间,且服从一般分布,第2站的输入过程一般是非更新过程,则难以用解析的方法求解第2站的平均排队时间。
下面使用指标比的方法刻画第1站对第2站的影响,的定义如下。
定义2 串联排队系统第2站的为
式中表示串联排队系统第2站的平均排队时间。
串联排队系统第2站的指标比为串联排队系统在情况下与的情况下第2站的平均排队时间之比。由(1)式知,。因此如果已知,则可以通过计算系统第2站的平均排队时间。下一节重点分析的值。
假设系统中2个服务器的服务时间均服从指数分布(即),外部到达过程为泊松过程,由文献[15]可知,第1站的离去过程为泊松过程,第2站的输入过程与系统的到达过程独立同分布,因此,。由式(1)得,当服务时间服从指数分布,第2站的为
当系统中2个服务器的服务时间都是常数(即)时,由于系统的平均排队时间只受瓶颈站影响[16],故
式中,,。由式(3)得。将及代入式(1),得服务时间为常数时第2站的为
令, ,。当系统的服务时间为常数时,与系统各参数的关系如图2所示。
Fig.2 Relationship between and system parameters
由图2可知:当系统的服务时间为常数时,随的增大而增大,随的增大而减小,的取值均小于1。
当系统中2个服务器的服务时间服从一般分布,很难求解指标比的解析表达式。为了分析指标比与各个参数之间的关系,模拟了4 860种不同参数组合下的串联排队系统,这4 860个系统的服务时间均服从Gamma分布,参数取值如下:,,, ,, 。
对于每组给定的参数,模拟运行30个样本,每个样本取第400 001个至第600 000个作业的平均排队时间,平均排队时间模拟值的置信水平大于95%,保证了模拟数据的可靠性。根据模拟得到的平均排队时间计算的值,数值结果如图3所示。图3显示是的单调增函数,且,的取值受如下系统参数的影响:。
Fig.3 Performance ratio with different parameters
为了便于分析,将服务时间服从一般分布的串联排队系统的第2站的指标比记为。下面采用控制变量法,分别讨论如下4种情况下指标比与各参数的关系:,,,。4种情况下指标比的模拟数值结果如图4所示,具有以下特性:(1)随的增大而减小;(2)关于近似线性递增,其中;(3)随的增大而增大,且具有一定的上凸性;(4)随的增大而增大,并且具有一定的下凸性。
Fig.4 Relationship between and system parameters
由于关于近似线性递增(),且线段上任意一点可以由2个端点线性组合,因此如果,其值可以表示为系统服务时间分别为常数(=0)和指数分布(=1)情况下的线性组合,即式(2)和式(4)的线性组合。记的近似值为,则
根据模拟得到的4 860个数据进行拟合,得
因此,的近似表达式为
注:(1)如果串联排队系统的服务时间服从指数分布或为常数,则。(2)虽然式(6)是根据服务时间服从Gamma分布的串联排队系统的数值结果拟合而成,但也适用于其他平方变异系数小于1的任意分布(见3.2节)。
第2.3节给出了指标比的拟合表达式,即式(6),下面通过模拟验证其有效性。令,,,表1分别比较了不同下第二个站的平均排队时间、指标比、拟合指标比,及和的相对误差。
表1 PR的相对误差
Table 1 Relative error of PR
| | | | |
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站的平均排队时间随的增大而增大,拟合指标比的相对误差均小于2%,拟合效果较好。此外,针对2.3节中模拟得到的4 860个与拟合值比较,发现相对误差小于5%的拟合指标比达到98.1%。
研究指标比的目的是求解第2站的平均排队时间。由式(1)和(6)得的近似值如下:
下面通过模拟服务时间服从Gamma分布、均匀分布和三角分布的串联排队系统,验证的准确性。
假设串联排队系统中2个站的服务时间均服从Gamma分布,选取34组不同参数的串联排队系统(见表2),分别比较不同参数下第2站的平均排队时间、近似值及两者的相对误差,如表2所示。表2中的34组数据结果显示近似值的相对误差均小于3%,运用指标比的拟合表达式可以较准确地计算平均排队时间。
表2 TQ2的相对误差
Table 2 Relative error of TQ2
类别 | | | | | | | | |
变 化 |
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 |
变 化 |
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 |
变 化 |
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 |
变 化 |
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 |
假设串联排队系统中第1站和第2站的服务时间分别服从均匀分布和,则,。
令,对串联排队系统进行仿真模拟,从而得到第2站的平均排队时间的模拟值。将模拟值与式(7)计算得到的近似平均排队时间进行比较,数值结果如图5所示。图5中平均排队时间的模拟值与拟合值几乎重合,误差较小。
Fig.5 Error analysis for case with uniform distribution
假设串联排队系统中第1站和第2站的服务时间分别服从参数为的三角分布(),其概率密度为
令,,,,,,则,,,。
令,对服务时间服从三角分布串联排队系统进行仿真模拟,得到第2站的平均排队时间的模拟值。图6对模拟值与式(7)计算得到的近似平均排队时间进行比较。图6显示通过仿真模拟得到的平均排队时间与拟合得到的平均排队时间误差较小,从而验证了通过指标比近似计算串联排队系统平均排队时间方法的可靠性。
Fig.6 Error analysis the case with triangular distribution
本文首次提出用指标比来刻画串联排队系统上游工作站对下游工作站的影响,并基于指标比提出串联排队系统平均排队时间的近似方法。通过分析指标比与系统各参数之间的关系,拟合得到指标比的近似表达式,利用指标比可以较准确地估计串联排队系统瓶颈站的平均排队时间。本文运用指标比的方法分析了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.
[百度学术]