Nikolaos Limnios

Queueing Theory 1


Скачать книгу

service times will be general, we will be using the discrete phase type distribution to represent them based on the fact any discrete general distributions, with finite support, can be represented by a Markov chain (Alfa 2004).

      Consider the simple Geo/Geo/1 system where the probability of arrival is a and of service completion is b. Let Xn(a, b), n = 0,1,2, ∙∙∙ be the number in the system at time n for a specified pair (a, b), with

for the case of a stable system. Based on the results in Alfa (2016), the associated discrete time Markov chain (DTMC) has the transition matrix P written as

      For a stable system with

, we have

      and

      where

is the expected value of
.

      Next we consider the cases of a = f(b). We require that the following condition applies all the time in this section:

      [1.1]

      We consider two main cases: one when f(∙) is a decreasing function and the second one when it is an increasing function.

      When we consider the case where a = f (b), we have

      and

      where

      Next, we study the behavior of this system as a function of b.

      1.2.1.1. Consider x0(b)

      By differentiating x0(b) with respect to b, we have

      PROPOSITION 1.1.– By studying the derivative of x0(b) w.r.t. b, we observe

       1) if f’(b) is negative, then an increase in b leads to an increase in the probability that the system is empty;

       2) if f’(b) is positive, then the probability that the system is empty could be increasing or decreasing, depending on whether respectively. Hence, an increase in the arrival rate as a function of the service rate does not necessarily mean a decrease in the emptiness probability.

      We obtain additional information by considering the rate of increase in the probability of empty system as a rate of increase in the service completion probability, i.e.

      PROPOSITION 1.2.– If f(b) is convex and f’(b) < 0, then x0(b) is definitely concave. if f(b) is decreasing fast, then x0(b) is increasing slowly, i.e. f(b) decreasing faster than increases. If, however, f(b) is convex and

then x0(b) is convex.

      Next let us consider E(X(b)). By differentiating E(X(b)) w.r.t. b, wehave

      Define

      For a stable system, we know that bf(b) and 1–f(b) are always positive, hence both A(b) and B(b) are always positive if

.

      We thus have that the behavior of E(X(b)) depends entirely on f’(b), A(b) and B(b) and whether

or not. E(X(b)) depends on the sign and values of the first three parameters, and the value of the fourth one.

      PROPOSITION 1.3.– For a stable system,

       1) if then E(X(b)) will be decreasing as b increases;

       2) if then E(X(b)) will be increasing as b increases.

      The first part of proposition 1.3 is straightforward in that an increase in the service rate reduces the arrival, thereby reducing the traffic intensity and decreasing the expected number in the system. However, an increase to the service rate could also lead to an increase in the mean number in the system as observed in the second part of the proposition, which can be explained by the rate at which the arrival probability increases compared to the service completion probability (b), especially if that rate is high.

      1.2.1.3. Special cases

      CASE 1: Consider the case where f(b) is a decreasing function of b, say for example a = f(b) = 1 – b2, with the strict condition that f(b) = 1 – b2 < b. For this, we have

      REMARK 1.1.– For the case where a = f(b) = 1 – b2, the arrival probability is a decreasing concave function of the service completion probability and the mean number in the system is decreasing in a convex form.

      CASE 2: Next consider the case where f(b) is an increasing function of b, say for example a = f(b) = b2, with the strict condition that f(b) = b2 < b. In this case, we have

      In this case, as b increases we see that E(X(b)) increases because the arrival probability increases faster than service completion increases, i.e. f′(b) = 2 b.

      REMARK 1.2.– For the case where a = f(b) = b2, with arrival probability increasing