B.3. Markovian Systems in Equilibrium

B.3.1. Discrete-time Markov Chains

A Markov chain discrete-time Markov chain - DTMC is a random process that undergoes transitions from one state to another on a state space. It must possess the “memorylessness” property , namely, that the probability distribution of the next state depends only on the current state and not on the sequence of events that preceded it. The material presented in this subsection will be useful for the study of queuing systems treated in the later subsections of Appendix B.3 and Chapters 4 and 5.
More formally, a stochastic process {Xn,n=0,1,2,}image that takes a finite or countable number of possible values is said to be in state iimage at time nimage if Xn=iimage. We assume that there is a fixed probability Pijimage that if the process is in state iimage it will next be in state jimage,

Pr{Xn+1=jXn=i,Xn1=in1,,X1=i1,X0=i0}=Pr{Xn+1=jXn=i}=Pij

image (B.2)

for all states i0,i1,,in1,i,jimage and all n0image. Such a stochastic process is known as a (discrete) Markov chain. The Markovian property, expressed by Eq. (B.2), is that the conditional distribution of any future state Xn+1image given the past states X0,X1,,Xn1image and the present state Xnimage is independent of the past states and depends only on the present state. Thus, a Markov chain is a sequence of random variables X1,X2,X3,image with the Markov property. Since the process must make a transition into some state, j=0Pij=1,i=0,1,,Pij0,i,j0image.
A Markov chain is said to be homogeneous in time if the transition probability between two distinct states depends only on the time step difference, i.e. Pr{Xm+n=jXm=i}=Pijnimage. Here, we focus on homogeneous Markov processes.
Eq. (B.2) defines the so-called one-step transition probabilities , namely probabilities that describe the transition of the chain in one time step. These transitions can be depicted via a state diagram, where a state diagram is a directed graph used to picture the state transitions. The states that the chain transitions to will be called neighboring states. Similarly, the nimage-step probabilities that the process will be in state jimage starting at iimage after nimage successive transitions (nimage successive time steps) are

Pijn=Pr{Xn+m=jXm=i},n,i,j0.

image (B.3)

In order to compute the nimage-step transition probabilities, the Chapman-Kolmogorov equations shown in the following can be employed:

Pijn+m=k=0PiknPkjm,for alln,m0,alli,j.

image (B.4)

If fijnimage is the probability that starting in state iimage, the first transition to jimage occurs at time nimage, then state jimage is said to be recurrent if fjj=1image, which means that starting at state jimage the process will return to state jimage with probability 1 for some time nimage (fjj=n=1fijnimage). A state is transient if fjj<1image. If state jimage is recurrent and μj=n=1nfjjnimage denotes the expected number of transitions needed to return to start jimage (mean recurrence time - expected return time in state jimage), then state jimage is positive recurrent if μj<image and null recurrent if μj=image. A state jimage is said to be accessible from a state iimage (written ijimage) if a system that started in state iimage has a nonzero probability of transitioning into state jimage at some point (Pr{Xn=jX0=i}=Pijn>0image for some n>0image). A state iimage has period kimage if any return to state iimage must occur in multiples of kimage time steps. If k=1image, then the state is said to be aperiodic, namely, returns to state iimage can occur at irregular times. A Markov chain is aperiodic if every state is aperiodic and an irreducible Markov chain only needs one aperiodic state to imply all states are aperiodic. A state iimage is called absorbing if it is impossible to leave this state.
A state iimage is said to be ergodic if it is aperiodic and positive recurrent, namely, if it is recurrent, has a period of 1, and it has finite mean recurrence time. If all states in an irreducible Markov chain are ergodic, then the chain is said to be ergodic. It can be shown that a finite state irreducible Markov chain is ergodic if it has an aperiodic state.
If an irreducible aperiodic Markov chain is positive recurrent, then we can define as πj=limnPijnimage, where Pijnimage is obtained from (B.3) for n=1image. Let Simage denote the set of possible states that the process can make a transition to. If there is a probability distribution over states πimage such that

πj=iSπiPr{Xn+1=jXn=i}

image (B.5)

for every state jimage and every time nimage then πimage is a stationary distribution of the Markov chain. An irreducible chain has a stationary distribution if and only if all of its states are positive recurrent. In that case, πimage is unique and related to the expected return time Mjimage by πj=CMjimage, where Cimage is a normalizing constant independent of the state. The chain converges to the stationary distribution regardless of the initial state distribution. Such πimage is called the equilibrium distribution of the chain. For a finite state space, a stationary distribution πimage is a (row) vector, whose entries are non-negative and sum to 1, it is unchanged by the operation of transition matrix Pimage and it is defined by πP=πimage. In addition, 0πj1image and jSπj=1image. Then the transition probability distribution can be represented by the transition matrixPimage, with the (i,j)imageth element of Pimage equal to πj=Pr{Xn+1=jXn=i}image.

B.3.2. Continuous-time Markov Processes

A continuous-time Markov process may be in one of the finite or infinite many (however countably infinite many) states at a time. The probability that the process transitions from state iimage to jimage at time t0+timage is defined as Pr{X(t0+t)=jX(t0)=j}image and it is independent of the behavior of X(t)image prior to t0image. If X(t)image is a homogeneous Markov process, the above probability depends only on the elapsed time required for the transition and not the initial time t0image. Thus in this case, Pij(t)=Pr{X(t0+t)=jX(t0)=i}image or Pij(t)=Pr{X(t)=jX(0)=i}image, denoted by state transition rates as opposed to the corresponding state transition probabilities in the discrete case (Pijnimage). The state transition rates are obtained as derivatives of the state transition probabilities at t=0image. Markov processes, in general, are characterized by the basic property that the past has no influence on the future if the present is specified and can be expressed in the continuous form as

Pr{X(n)=xnX(n1),.X(1)}=Pr{X(n)=xnX(n1)}.

image (B.6)

All Markov processes share the interesting property that the time it takes for a change of state, called sojourn time is an exponentially distributed random variable [174].
The initial probability distribution of state iimage is given by Pi(0)=P{X(0)}=iimage, defined analogously for the rest of states, and the unconditional probability that the process is in state jimage at time timage is defined as Pj(t)=iPi(0)Pij(t)image. For arbitrary timage and simage, the continuous version of the Chapman-Kolmogorov equations is given by

Pij(t+s)=P{X(t+s)=jX(0)=i}=kPik(t)Pkj(s).

image (B.7)

B.3.3. Birth-and-Death Processes

The birth-death process is a special case of continuous-time Markov chain (CTMC), where the state transitions are of two types only, namely, “births,” which increase the state variable by one and “deaths,” which decrease the state by one [174]. A typical example is that of some banks which have a security door that allows only one person at a time to enter/leave the bank. In this case, irrespective of the system state, at each time instant the next event is either the arrival of a customer in the bank, or the departure of a customer that completed the required service and leaves the building. In these and similar cases, when a birth occurs, the process goes from state nimage to n+1image. When a death occurs, the process goes from state nimage to state n1image. The state diagram evolution for such a process is depicted in Fig. B.3. Such processes are specified by their birth rates {λi}i=0image and death rates {μi}i=1image. A pure birth process is a birth-death process where μi=0for alli0image, while a pure death process is a birth-death process where λi=0for alli0image. A (homogeneous) Poisson process is a pure birth process where λi=λfor alli0image.
Fig. B.3
FIGURE B.3 State diagram for the birth-death process.
The global balance equations (also known as full balance equations) are a set of equations that in principle can be solved to give the equilibrium distribution of a Markov process (when such a distribution exists). For a continuous-time Markov process with state space Simage, transition rate from state iimage to jimage given by qijimage and equilibrium distribution given by πimage, the global balance equations are given for every state iimage in Simage by

jS{i}πiqij=jS{i}πjqji,

image (B.8)

where πiqijimage represents the probability flux for state iimage to state jimage.
For a CTMC with transition rate matrix Q=[qij]image, if πiimage can be found such that for every pair of states iimage and jimage, πiqij=πjqjiimage holds, then the global balance equations are satisfied and πimage is the stationary distribution of the process. If such a solution can be found, the resulting equations are usually much easier than directly solving the global balance equations. Furthermore, a CTMC is reversible if and only if the detailed balance conditions are satisfied for every pair of states iimage and jimage.
In queuing theory, the birth-death process is the most fundamental example of a queuing model. In the following, we provide a brief overview of simple and characteristic birth-death queuing systems that have particular interest for the analysis of queuing networks and some of the methodologies presented in this book, especially in Chapter 4. The rest of the systems presented in this subsection with an infinite queue are all characteristic examples of birth-death systems defined in this subsection.

B.3.4. The MM1image Queuing System

The MM1image queue is a system having a single server, where arrivals are determined by a Poisson process and job service times have an exponential distribution. Arrivals are assumed to occur at rate λimage and service times have an exponential distribution with mean 1μimage. The single server operates according to a FIFO discipline. When the service of a customer is complete, the customer at the head of the queue leaves the storage space entering the server, the customer that completed service leaves the systems, and the number of customers in the system reduces by one. The buffer (queue) of the system is of infinite size.
The model is considered stable only if λ<μimage. If, on average, arrivals happen faster than service completions the queue will grow indefinitely long and the system will not have a stationary distribution. Assuming ρ=λμimage for the utilization of the buffer and require ρ<1image for the queue to be stable, ρimage represents the average proportion of time for which the server is occupied. In order to see this more intuitively, consider the version of Little’s law with respect to the server, which gives the average number of customers in the server Nwimage with respect to the arrival rate λimage and the mean service time 1μimage, Nw=λμ=ρimage. For the MM1image system to be stable, at most one customer at a time should be in service, as expected by the definition of its operation.
The corresponding state diagram is of birth-death type, similar to Fig. B.3. Useful closed-form expressions may be obtained for various quantities of interest. More details on the analysis methodology of the MM1image queue can be found in [25,142], etc. The most useful results obtained include the following:
• The steady-state probability distribution is given by pn=ρn(1ρ),n=0,1,image
• The average number of customers in the system in the steady state is N=λμλimage.
• The average system delay per customer (waiting time plus service time) is T=1μλimage.
• The average waiting time in queue is W=ρμλimage.
• The average number of customers in the queue is NQ=ρ21ρimage.
Additional quantities of interest can be obtained via analysis as shown in [25] and other references targeted explicitly to queuing theory analysis [33,36,90,142,209,211,229].

B.3.5. The MMmimage System and Other Multiserver Queuing Systems

The MM1image queue is a very simple system that emerges frequently in various occasions. However, oftentimes, multiserver systems are employed. In this subsection, we overview the most fundamental of them. As with the MM1image system, for all these multiserver systems, it is possible to derive a set of equations from the corresponding state transition diagram and solve for the steady-state (equilibrium) probabilities. Using Little’s law allows computing the average delay per customer and other quantities of interest.

The MMmimage queuing system

The MMmimage system is similar to the MM1image with respect to the arrival process (Poisson), the independence of service times, and the exponentially distributed and independent interarrival times. However, this system has mimage parallel servers, meaning that a customer at the head of the queue (ready to receive service) can be routed to any of the server currently available, once a customer completes service in one of them.
Applying the analysis methodology of balance equations, the corresponding outcomes for the quantities of interest obtained in the MM1image case are as follows:
• The probability that an arriving customer will find all servers buys and will be forced to queue is PQ=p0(mρ)mm!(1ρ)image, where p0image is given by the expression

p0=[1+n=1m1(mρ)nn!+n=m(mρ)nm!1mnm]1.

image (B.9)

• The average number of customers in the system in the steady state is N=mρ+ρPQ1ρimage.
• The average system delay per customer (waiting time plus service time) is T=1μ+PQmμλimage.
• The average waiting time in queue is W=ρPQλ(1ρ)image.
• The average number of customers in the queue is NQ=PQρ1ρimage.
Now the utilization of the system is defined as ρ=λmμ<1image. In order to understand this stability condition for the system, note that now mimage servers operate in parallel each with an exponentially distributed service rate of μimage. This means that the combined server will be operating as a single server with service rate which will be the minimum of mimage exponential distributions. This has been proven to be again exponentially distributed with rate the sum of partial rates mμimage. Thus, for the mean service time expression in the similar application of Little’s law that we performed for the MM1image system, now for the case for the MMmimage would yield λmμimage.

The MMimage queuing system

In several cases, the available service has such a capacity that practically, any customer entering the system receives immediate service. Typical examples include web platforms receiving requests and serving them immediately (up to a certain point of course). In practice, the number of available servers is always finite. However, whenever the number of available servers is such that no queuing at all is observed, the system might well be considered as of an MMimage type.
The MMimage system can be considered as the limit of the MMmimage system when mimage. Consequently, in the limiting case where the number of servers is high (theoretically infinite number of servers), m=image, one obtains a MMimage system for which the quantities of interest become the following:
• The average number of customers in the system in the steady state N=λμimage.
• The average system delay per customer (waiting time plus service time) T=1μimage.
Clearly, for such a system where a customer always finds an available empty server, there is no queuing and thus NQ=0,W=0image. The response time for each arriving job is a single exponential distribution with parameter μimage as already explained and it thus coincides with the average system delay per customer given above. One should also notice that the stationary distribution of a more general queue with practical applications, the MGimage queue where the service discipline can be arbitrary, but arrivals are Poisson and there exists an infinite number of servers, is the same as that of the MMimage queue. This implies that an overprovisioned system, i.e. one with practically infinite number of servers, behaves the same irrespective of the service policy and the behavior is rather simple, since each customer is essentially served immediately. Thus, the mean number of customers in the system will be N=λμimage in any case.
The steady-state distribution of the system is given by the expression

pn=(λμ)neλμn!,n=0,1,

image (B.10)

with p0=eλμimage. From the above expression, it can be seen that the system is always stable, as no special constraint on the latter distribution is required. This is also intuitively verified by the fact that no matter how heavy the input load of customers, there are always available servers to satisfy it (infinite number). Thus, the number of customers waiting in the systems never grows arbitrarily and the system will be stable. In fact, it can be shown that the number of customers in the system is Poisson distributed even if the service time distribution is not exponential, e.g. in MGimage[25]. The MMimage system is also known as the Erlang-C model, as opposed to the Erlang-B that will be described in the following.

The MMmmimage queuing system

This system is identical to an MMmimage system with the difference that there is no queuing buffer space. The system consists of a facility where arrivals form a single queue and are governed by a Poisson process, there are mimage servers and job service times are exponentially distributed. Thus, the total capacity of the system is defined by the number of customers the system can serve instantaneously, which also means that when an arriving customer finds all mimage servers of the system busy, it is blocked from entering the system and thus discarded. For this reason, the model is called a loss system. The lost customers are said to be blocked. Practical examples of such systems are telephone exchanges, cellular phones, and other more general flow-service based systems.
In general, it is tough to obtain analytic expressions for the quantities of interest mentioned above for the other systems presented. However, it is possible to obtain the probability of iimage customers in the system and the blocking probability pmimage, namely, the probability that an arriving customer will find all servers busy and will be discarded from the system

pm=(λμ)mm!n=0m(λμ)nn!.

image (B.11)

The latter is known as the Erlang B formula. This is also called the Erlang-B loss function and tabular solutions are usually provided in the related references. More details can be obtained in the corresponding references [25,33,36,90,142,209,211,229]. It should be noted that as mimage, the above distribution approaches the Poisson distribution of the MMimage system, described above.
..................Content has been hidden....................

You can't read the all page of ebook, please click here login for view all page.
Reset