Pollaczek-Khinchin formula for the
M / G / I
queue in
discrete time with vacations
W.C.Chan R.-J.Chen T.-C. LU
Indexing term.s: Pollaclrk-Khinchiii formula, IM/G/I queue in discrete time
Abstract: The continuous-time MIGI1 queue with
vacations has been studied by many researchers. In the paper the authors report on an investigation of the discrete-time M/G/l queue using Little's formula and conditional expectation. This direct approach can also be adopted to study the continuous-time case.
1 Introduction
Consider a discrete-time M/G/l queue in which the server begins a vacation of random length each time that the system becomes empty. If the server returns from a vacation to find one or more customers waiting, the server works until the system empties, and then begins another vacation. If the server returns from a vacation to find no customers waiting, the server begins another vacation immediately. We assume that the lengths of vacations are independent and identically distributed (i.i.d.) random variables and are independ- ent of the arrival process as well as the service times of customers.
The continuous-time model has been analysed in a number of papers [l-61. A key result from these analyses is that the number of customers present in the system at a random point in time in equilibrium is the sum of two independent random variables: the number of Poisson arrivals during a time interval of the residual vacation time, and the number of customers present at a random point in time in equilibrium in the corresponding standard MIGI1 queue.
The purpose of this paper is to present a discrete- time model using a simple and direct method. The Pollaczek-Khinchin (P-K) mean-value formula, either in continuous-time or discrete-time models, can be derived using Little's formula and conditional expectation.
0 IEE, 1997
IEE Proceedings online no. 1997 1225
Paper first received 13th June 1996 and in revised form 21st February 1997
W.C. Chan is with the Department of Electrical and Computer Engineer- ing, The University of Calgary, Alberta, Canada T2N IN4
R.J. Chen and T.C. Lu are with the Department of Computer Science and Inlormation Engineering, National Chiao Tung University, Hsinchu. Taiwan
2 Discrete-time M / G / I model
A discrete-time queueing system, such as an ATM network, is characterised by time-slotted and synchronous services [7]. The time axis is divided into equal intervals, called slots. Arriving customers that find the server busy wait in a queue. Servicing of customers is synchronised to start only at slot boundaries.
We will use the following definition for the discrete- time models. Departures take place only at slot bound- aries. Without loss of generality, we normalise the length of a slot to unit time. Slots are sequentially numbered in nonnegative integers so that thejth slot is located in the time interval [i - 1, j), where j = 1, 2,
_ _ _
. The two time points immediately before and after time j are denoted j - and j'. A customer completing service in slot j will leave the system at timej- and a customer starting service in slot (j +1) will begin the service at time j.The service times of customers are i.i.d. random vari- ables with probability density function
p k = P{service time equals k slots} k
2
INow
we define two discrete-time models for two differ- ent memoryless arrival processes.Dejhition
I:
Model A (bulk arrivals at slot boundaries) Bulk arrivals occur at slot boundaries with a Poisson distributed size of meanA
customers, that is,f ( k ) = P { k customers arrive at a slot boundary}
X"-X
k ! k = 0 , 1 , 2 , . . . __
-
Definition 2: Model B (exponential interarrival time) Interarrival times are identical and exponentially dis- tributed with mean l / A slots. The probability density function of the interarrival time is
f(x) = Xe-'" 2 2 0
Equivalently the number of arrivals in a slot has the same Poisson distribution as that of model A, How- ever, customers in model B have to wait, on average, an additional half-slot more than that of customers in model A.
Proposition 1: The mean service time of customers in model B is a half-slot more than that of customers in model A.
3
Consider the test customer in model A who arrives at a discrete-time MiGil queueing system and finds the server busy. The time interval between the arrival epoch of this test customer and the completion of the service of the customer being served is known as the residual service time R as shown in Fig. 1.
Residual service time in discrete-time model
I , , , , ,
R #
I , , ,
From Fig. 1 we see that the test customer must wait at least an amount of time equal to the residual service time R. We shall determine the probability distribu- tions of both the special service time Y and the residual service time R in terms of the distribution pk of the general service time X and of its first and second moments (m, and m2, respectively).
Definition 3: Define the following random variables: X = the general service time with probability distribu-
tion { p k , k = 0 , I , ...}, the first moment m l and second moment m2;
Y = the special service time within which the test cus- tomer arrives; and
R
= the residual service time.By definition, P k = P { X = k } is the probability that X equals exactly k slots in length. Now let
4 k = P{Y = k } k = 0,1,.
.
.be the probability that Y equals exactly k slots in length.
Proposition 2: Given { p k ) , the probability distribution of Y can be expressed as
where m1 is the first moment or mean of X ,
Prooj Since the arrival process is Poisson, which is independent of the service process, the probability of the event { Y = k } is proportional to k and p k . Thus qk
can be expressed as the product of kpk multiplied by a proportional constant c,
q k =cl%pl, k = O , l ,
. . .
where the right-hand side expresses the linear weighting with respect to the length of Y and includes a constant
c which must be evaluated so as to properly normalise this probability distribution (see p.171 of [SI). It follows that CO k=O k=O Thus 1 c = - ml
where ml =
XEo
kpk, is the mean of XProposition 3: Let r, denote the probability that the residual service time R equals exactly k slots in length.
IEE Proc -Comput D g i t Tech Vol 144 N o 4 July 1997
Then
r j = P { R = j }
Prooj! Suppose that the special service time Y is exactly k slots in length and that the arrival process is Poisson. The arrival epoch of the test customer occurs at any one of the k time points of Y with equal probability 1/ k. In other words, given that Y = k the conditional probability of R being exactlyj slots in length is
It follows that the joint probability
P { R = j , Y = k } = P { R = jlY = k } P { Y = k } 4 k k - P k ml r3 =
P { R
=j }
- - - __ - Therefore CO = P { R = j , Y = k } k = j + l CO= E "
k=3+1 m1Applying the results of proposition 2 and proposition 3
we obtain the mean value of the residual time as fol- lows:
Theorem 1: The mean residual service time R A for model A is given by
where ml, m2 are the first and second moments of X , respectively.
Proof: By definition, we write
00 R A =
Cjr3
j=1 1 = - C j P { X>
j } proposition 2 ml j=o - m2 -ml - 2mlApplying the results of proposition 1 and theorem 1,
we obtain the mean residual service time for model B as follows:
Theorem 2: The mean residual service time
RB
for the discrete-time model B is given bym2 2ml -
RE = ~
where m l , m2 are the first and second moments of X ,
respectively.
4 Pollaczek-Khinchin mean-value formula for the M / G / I queue in discrete time
We use the conditional expectation and Little's formula to derive the P-K formula. Model B is adopted here because its mean residual service time has the same form as the continuous-time model. The method devel- oped in this Section is also useful for deriving the P-K formula in continuous time.
Depending on the state of the server, the waiting time of an arriving customer can have two different values. Let us define the following events:
B = the arriving customer finds the server busy and N, customers waiting in the queue;
I the arriving customer finds the server idle
Let W denote the waiting time. The average waiting time of a customer can be expressed as
where p = A E [ d is the probability that the server is busy. E [ W ] = p E [ W / B ]
+
(1 - p ) E [ W / I ] (1) test customer o r r i v i n g time time~y-x,,q
W.%j
Fig.2 Waiting time for. the test customer
For the event B, we have the waiting time as shown in Fig. 2 and E'[WIB] =
E[R
+
x1
+
x2+ . .
'+
X v q / B ] = E [ R / B ]+
E [ X 1+
X 2+
. .
+
Xhr,iB] = E [ R / B ]+
E [ X / B ] E [ N , / B ] = E [ R / B ]+
E [ W ] = E [ R / B ]+
E[X]E[1Vq]/p Thus we find E[WIB] = E [ R / B ]+
E [ W ] where X N q = 0 if Nq = 0 and E [ X N , p ] = E [ X ] = 7x1 E[NIBl = E[N,lp = W W I / phave been used. Note that
X
and B are independent random variables.For the event I, the waiting time becomes zero because the arriving customer is served immediately. Thus
E[WII] = 0
E[WI = P(E[RIBl
+
E [ Wand hence the second term on the right-hand side of eqn. 1 vanishes. Therefore, we obtain from eqn. 1 or
Substituting the average residual service time R B for E[RIB] in this expression yields the desired Pollaczek- Khinchin formula.
To calculate the mean number of customers in the sys- tem, we have to obtain the total system time T , which is the sum of the waiting time and the service time. It follows that
E [ T ]
=E [ W ]
+
E [ X ]
Let N be the mean number of customers in the system. By Little's formula, we have
= XE[T]
(4)
4.7 Mean waiting time in the M/G/I queue with vacations
Generally a queueing system will have busy periods with at least one customer present and idle periods with no customer present. A busy period is a time interval that begins when an arriving customer finds the system empty. An idle period is the period between two successive busy periods. Clearly, busy and idle periods occur alternatively and form a cycle.
Suppose that at the end of each busy period, the server goes on vacation for a random interval of time with first moment vi and second moment v2. For com- puter communication networks, vacations correspond to transmissions of various kinds of control and record-keeping packets when there is little traffic or
when the transmit queue is empty.
We shall derive an expression for the average waiting time in the MIGIl queue with vacations. In this case the second term on the right-hand side of eqn. 1 is no longer zero and is simply the delay caused by the resid- ual vacation interval equal to (1 - p)v2/2v,. Now from eqns. 1 and 2 and theorem 2 we obtain
Solving for E[W yields the mean waiting time for model B:
(5)
This expression may be regarded as an extension of the Pollaczek-Khinchin formula (eqn. 3) for the mean waiting time in the M/G/l queue in discrete time to the case with vacations.
5 Application
Consider a polling system for transmission of cells from several cell-based streams into a statistical multi- plexing system where the cell size for each stream is constant. This situation arises often in multiaccess channels.
As a typical example, we consider a communication IEE PIOC -Camput Dzgit Tech, Val 144, No 4 July 1997
channel that can be accessed by m spatially separated stations. These m stations are connected by cables in a unidirectional loop. Each station transmits the back- logged packets when it is polled. The interpolling times can be regarded as vacations. A station transmits pack- ets when it receives the poll.
We will consider m homogeneous traffic streams. Each consists of three Poisson arrival streams with rates AV, A, and
Ad
which correspond to the arrival rate of video, audio and data packet streams, respectively. For simplicity, we assume the lengths of these packets are l,, 1, and ld cells for each type of stream (Fig. 3). The total packet arrival rate of the combined stream is and the total cell arrival rate of the combined stream isA s = Aa f
+
Ad=
Asia
+
A v l v f A d l dThe total packet arrival rate of all the m streams is then m 5 .
A" A d A 0
dota flow
J-Lv- L i d - L t a -
time
Fig.3 Input packets for a single stream
Consider the network as an MiGI1 queue. The mean residual service time is given by
where
and
Now we calculate the mean vacation time v1 which is the sum of the other m - 1 mean busy periods and m token releasing time z. Initially, there are A5q packets in the MiGI1 queue. Consider the successive arrivals for a last-in first-out (LIFO) queueing discipline [8]. The mean busy period is equal to AJvlBM,G,l, where BM,GII is the mean busy period in the MiGI1 queue. From [[8], p. 2331 we have
E[XI
1 - P
B M / G / l = -
Then the mean vacation time can be expressed as vI = ( m - l)X,ulBM/G/l
+
m r or - m7 711 = 1 - ( m ~ 1 ) X s B ~ / ~ / iIEE Proc.-Comput. Digit. Tedi., Vol. 144, No. 4, J~rly 1997
mr -
-
1 - ( m - 1)J- 1 - P - m-ru - P ) - 1 - m pwhere p is the utilisation factor of the M/G/1 queue with arrival rate A,s. For exponentially distributed vaca- tion times with mean vl, we have v2 = 2vI2. From eqn. 5 we obtain the mean waiting time
6 Discussion
It is interesting to note that if model
A
is adopted, then using similar arguments as in Section 4 would result in the following expressions:The Pollaczek-Khinchin formula in discrete time becomes P -E [ W ] =
-
RA 1 - P - X ( m 2 - m l ) - 2(1 - P )The mean number of customers in the system is
and the mean waiting time in the MIGII queue with vacations now becomes
The results of eqns. 6, 7 and 8 for model A correspond to eqns. 3, 4 and 5 for model B, respectively.
7
ConclusionsIn the past, most queueing analysis has been based on queueing phenomena in continuous time. Recently in the telecommunication industries, B-ISDN (Broadband Integrated Services Digital Network) has received con- siderable attention. B-ISDN can provide a common interface for future communication including the trans- mission of video, data and speech. Since information in B-ISDN is transported by means of discrete units of 53-octet ATM (Asynchronous Transfer Mode) cells, it appears that the analysis of a queueing system in dis- crete time is more natural.
In this paper, we derive the Pollaczek-Khinchin for- mula for the M/Gi1 queue in discrete time with vaca- tions. We simply use Little's formula and conditional expectation. This mean value analysis provides an alternative method of deriving performance measures for either continuous-time or discrete-time MiGI1 queues.
8 Acknowledgment 3 4
The authors wish to thank the referees for their careful reviews and constructive comments which improved the
quality of the paper. 5
226
6 References
7 COOPER, R.B.: ‘Queues served in cyclic order: Waiting times’, Bell Syst. Tech. J., 1970, 49, pp. 399413
LEVY, Y., and YECHIALI, U,: ‘Utilization of idle time in an M i Gil queueing system’, Manage. Sei., 1975, 22, pp. 202-211 8
HEYMAN. D.P.: ‘The T-oolicv for the MIGI1 aueue’. Munape I ;
sei., 1971, i3, pp. 775-718
Y
SCHOLL, M., and KLEINROCK, L.: ‘On the MIGIl queue with rest period and certain service-independent queueing dkiplines’,
Oper Res., 1983, 31, pp. 705-719
FUHRMANN, S.W.: ‘A note on the MIGI1 queue with server vacations’, Oper Res., 1984, 32, (6), pp. 1368-1373
FUHRMANN, S.W., and COOPER, R.B.: ‘Stochastic decompo- sitions in the M/G/l queues with generalized vacations’, Oper. Res., 1985, 33, (S), pp. 117-1129
BRUNEEL, H., and KIM, B.: ‘Discrete-time models for commu- nication systems including ATM’ (Kluwer Academic, Boston,
1993)
KLEINROCK, L.: ‘Queueing systems, vol. 1: theory’ (John Wiley & Sons, New York, 1975)