Birth-Death Processes

Birth-Death Processes: Transient Solution Poisson Process: State Distribution Poisson Process: Inter-arrival Times

Dr Conor McArdle EE414 - Birth-Death Processes 1/17 Birth-Death Processes

Birth-Death Processes are Markov chains where transitions are allowed only between neighbouring states of the process.

They may be discrete- or continuous-time. Our primary interest is in the continuous-time case.

Birth-Death processes are used to simplify the analysis of systems that see random arrivals and departures of customers (packets, requests, or calls etc.) over time.

We will develop results that are useful for performance evaluation of servers (or transmission lines) with associated input queues. The results also apply to circuit switches with limited numbers of outgoing channels.

We start in this section by seeking the transient solution to the Birth-Death system. Recall that the general transient solution π(t) of a continuous-time is given implicitly by π˙ (t) = π(t)Q where Q is the matrix of state transition-rates between each pair of states in the process. The Birth-Death assumption imposes a certain structure on Q which gives the solution π(t) a simplified form.

Dr Conor McArdle EE414 - Birth-Death Processes 2/17 Birth-Death Processes: Transition-rates

Consider the state transition-rate diagram of the general Birth-Death process below.

ë0 ë1 ëk-1 ëk

0 1 2 k-1 k k+1

ì1 ì2 ìk ìk+1

The states of the process are represented by numbered circles, the number indicating the state. For a birth-death process this number represents the number of individuals in the population, which for our purposes we interpret as the number of packets (or requests or calls) in the system.

Arrows represent the possible transitions between states and the labels on the arrows are the state transition-rates between states.

Transitions which relate to births (an increase in the population by one) are labeled with the “birth rate” (λk = birth rate when process is in state k).

Transitions which relate to deaths (a decrease in the population by one) are labeled with the “death rate” (µk = death rate when process is in state k).

Dr Conor McArdle EE414 - Birth-Death Processes 3/17 Birth-Death Processes: Transient Solution

Let us construct the matrix of transition rates Q from the Birth-Death state transition-rate diagram.

The component qij of Q is the transition rate from state i to state j (i 6= j). For j = i + 1 the rate is given by the birth-rate λi and for j = i − 1 the rate is given by the death-rate µi.

Given the stipulation that transitions are only allowed between neighbouring states, all rates for larger transition steps in the process (e.g. qij where j = i + 2) are zero. P P The rates at each state i must balance, that is ∀j qij = 0, so qii = − ∀j,i6=j qij. −qii may be interpreted as the total rate at which state i is exited to other states.

Forming the matrix Q from the above rates gives:   −λ0 λ0 0 0 ...  µ1 −(λ1 + µ1) λ1 0 ...    0 µ2 −(λ2 + µ2) λ2 ... Q =    0 0 µ3 −(λ3 + µ3) ...  . . . . .  ......

Dr Conor McArdle EE414 - Birth-Death Processes 4/17 Birth-Death Processes: Transient Solution

With this Q matrix, the transient solution equation

π˙ (t) = π(t)Q,

can be written as the set of differential equations:

dπ (t) 0 = −λ π (t) + µ π (t) k = 0 (1) dt 0 0 1 1 dπ (t) k = λ π (t) − (λ + µ )π (t) + µ π (t) k ≥ 1 (2) dt k−1 k−1 k k k k+1 k+1

The solution to these equations is not trivial and, for general sets of {λk} and {µk}, it is often the case that numerical methods are required to find a solution.

However, we will consider an important Birth-Death process, the Poisson process, for which we can easily find an exact (analytic) solution.

The Poisson process is used extensively to model the arrival times of packets (or calls) to a system.

Dr Conor McArdle EE414 - Birth-Death Processes 5/17 Birth-Death Processes

Birth-Death Processes: Transient Solution Poisson Process: State Distribution Poisson Process: Inter-arrival Times

Dr Conor McArdle EE414 - Birth-Death Processes 6/17 The Poisson Process

Consider a pure birth system where µk = 0 for all k and λk = λ for all k.

Thus our transition rate matrix Q is:

−λ λ 0 0 ...  0 −λ λ 0 ...    0 0 −λ λ . . .    0 0 0 −λ . . .  . . . . .  ......

and the system of equations π˙ (t) = π(t)Q is:

dπ (t) 0 = −λπ (t) k = 0 (3) dt 0 dπ (t) k = λπ (t) − λπ (t) k ≥ 1 (4) dt k−1 k We assume the initial condition (initial state distribution) for the system is:

 1 k = 0 π (0) = k 0 k 6= 0

Dr Conor McArdle EE414 - Birth-Death Processes 7/17 The Poisson Process

Thus (from the usual differenetial equation solution methods) the solution for π0(t) is

−λt π0(t) = e .

Using this solution in equation (4) for k = 1 gives:

dπ (t) 1 = λe−λt − λπ (t) dt 1 The solution to this differential equation is:

−λt π1(t) = λte

Continuing in this manner we have an expression for the transient solution of the pure birth process. (λt)k π (t) = e−λt k ≥ 0, t ≥ 0 (5) k k!

This pure-birth process is called the Poisson Process and equation (5) is the pmf of the state probabilities at time t (often called the Poisson distribution).

Dr Conor McArdle EE414 - Birth-Death Processes 8/17 The Poisson Process

At some fixed time τ, equation (5) gives the state distribution of the process at that time (a pmf).

For some chosen state n, equation (5) gives the probability of being in that state as time progresses from 0 to ∞ (a pdf).

¼

0231 4

Figure: The Family of Curves πk(t) as a function of k and λt

Dr Conor McArdle EE414 - Birth-Death Processes 9/17 The Poisson Process: Interpretation/Uses

The Poisson process is widely used for modelling the arrival process of packets (or calls) to a system.

If we assume arrivals of packets to a system are independent events and that each arrival occurs at the same average rate λ, then we have a Poisson process, where the state of the process corresponds to the number of arrivals since time 0.

Thus equation (5) gives a complete probabilistic description of the number of arrivals during a given time period and also a probabilistic description of how long the process stays in any one state, or alternatively, how long we expect between arrivals or until the next arrival.

We can note that a Poisson process arises from combining an infinite number of independent stationary (renewal) processes, each generating arrivals (births) at some rate. In the limit, the overall arrival process has a constant arrival rate (birth rate) λ. When we have any large population of such independent processes (e.g. a large number of people making phone calls through a switch), the Poisson process provides a very useful approximation to the overall arrival process.

We will now examine the previously mentioned probabilistic descriptions.

Dr Conor McArdle EE414 - Birth-Death Processes 10/17 The Poisson Process: Moments

Consider the random variable N(t), associated with the Poisson process, describing the number of arrivals (births) in the interval [0, t)

The pmf of N(t) is the state distribution of the Poisson process at time t, that is, the pmf of N(t) is (λt)k p (t) = π (t) = e−λt k ≥ 0, t ≥ 0 k k k!

We would expect that the average number of arrivals in this interval would be λt, given that λ is the arrival rate. So the of N(t) should be λt. We can check this: ∞ X E[N(t)] = kpk(t) k=0 ∞ ∞ ∞ X (λt)k X (λt)k X (λt)k = e−λt k = e−λt = e−λtλt k! (k − 1)! k! k=0 k=1 k=0

x x x2 Noting that the Taylor series expansion of e is: e = 1 + x + 2! + ...

⇒ E[N(t)] = λt

Dr Conor McArdle EE414 - Birth-Death Processes 11/17 The Poisson Process: Moments

We may also calculate the variance of the number of arrivals in [0, t). We first calculate: ∞ X E[N(t) · (N(t) − 1)] = k(k − 1)pk(t) k=0 ∞ X (λt)k = e−λt k(k − 1) k! k=0 ∞ X (λt)k−2 = e−λt(λt)2 (k − 2)! k=2 ∞ X (λt)k = e−λt(λt)2 k! k=0 = (λt)2

Using this result and the expression for E[N(t)] we can calculate:

Var(N(t)) = E[N 2(t)] − E2[N(t)] = E[N(t) · (N(t) − 1)] + E[N(t)] − (E[N(t)])2 = (λt)2 + λt − (λt)2 ⇒ Var(N(t)) = λt

Dr Conor McArdle EE414 - Birth-Death Processes 12/17 Birth-Death Processes

Birth-Death Processes: Transient Solution Poisson Process: State Distribution Poisson Process: Inter-arrival Times

Dr Conor McArdle EE414 - Birth-Death Processes 13/17 The Poisson Process: Inter-arrival Times

We may also calculate the distribution of the times between arrivals (births) for a Poisson process. Let τ be the random variable describing the time between adjacent arrivals (the inter-arrival time) for a Poisson process. We may set t = 0 when an arrival occurs and so, equivalently, τ is the random variable describing the time until the next arrival. ¼ Arrivals time

We know the probability of having no arrivals in [0, t), that is Pr(τ > t), from the Poisson state distribution: −λt π0(t) = e = Pr(τ > t) so the probability of the complementary event (τ ≤ t) is 1 − e−λt.

This gives the cdf of τ, the inter-arrival time: −λt Fτ (t) = Pr(τ ≤ t) = 1 − e t ≥ 0.

The pdf of τ is given by the derivative of Fτ (t): −λt fτ (t) = λe t ≥ 0

Dr Conor McArdle EE414 - Birth-Death Processes 14/17 The Poisson Process: Inter-arrival Times

We may calculate the expected value of τ as: Z ∞ Z ∞ −λt 1 E[τ] = tfτ (t) dt = tλe dt = 0 0 λ This is the average time between arrivals, which is as expected since the arrival rate is equal to λ.

To summarise, a pure-birth process with constant birth rate λ is described by a Poisson process whose mean state is λt. The times between births (arrivals) are exponentially distributed with mean 1/λ. The Poisson process approximates the combined arrival process of a large group of independent arrival processes.

Dr Conor McArdle EE414 - Birth-Death Processes 15/17 The Poisson Process : The Memoryless Property

Consider the Poisson arrival process where the inter-arrival time is the random variable τ. Suppose that an arrival has occurred at time zero. Let time t0 pass with no arrivals. What is the probability that the next arrival will occur within t seconds from now? This is the conditional probability:

Pr [τ ≤ t + t0 | τ > t0] Pr [t < τ ≤ t + t ] = 0 0 Pr [τ > t0] Pr [τ ≤ t + t ] − Pr [τ ≤ t ] = 0 0 Pr [τ > t0] 1 − e−λ(t+t0) − (1 − e−λt0 ) = 1 − (1 − eλt0 ) = 1 − e−λt

So the distribution of remaining time until the next arrival, given that t0 seconds has elapsed since the last arrival, is identically equal to the unconditional distribution of the inter-arrival time. This fact is called the memoryless property of the . The exponential distribution is the only continuous distribution with this property.

Dr Conor McArdle EE414 - Birth-Death Processes 16/17 Birth-Death Process: Conclusion

We have found the form of the general transient solution of the Birth-Death process and found a closed-form solution for the pure birth (Poisson) process.

The Poisson process is useful for describing arrival processes.

We will next turn to Birth-Death processes which describe queueing systems. The processes we will study have stationary state distributions and we are most interested in performance measures relating to mean values of their stationary distributions.

Dr Conor McArdle EE414 - Birth-Death Processes 17/17