next up previous
Postscript version of these notes

Example: tex2html_wrap_inline259 . Then

displaymath261

and the chain is otherwise specivied by tex2html_wrap_inline263 and tex2html_wrap_inline265 . The matrix tex2html_wrap_inline267 is

displaymath269

The backward equations become

align26

while the forward equations are

align40

Add tex2html_wrap_inline271 first plus tex2html_wrap_inline273 third backward equations to get

displaymath275

so

displaymath277

Put t=0 to get tex2html_wrap_inline281 . This gives

displaymath283

Plug this back in to the first equation and get

displaymath285

Multiply by tex2html_wrap_inline287 and get

displaymath289

which can be integrated to get

displaymath291

Alternative calculation:

displaymath269

can be written as

displaymath295

where

displaymath297

displaymath299

and

displaymath301

Then

align102

Now

displaymath303

so we get

align124

where

displaymath305

Notice: rows of tex2html_wrap_inline307 are a stationary initial distribution. If rows are tex2html_wrap_inline309 then

displaymath311

so

displaymath313

Moreover

displaymath315

Fact: tex2html_wrap_inline317 is long run fraction of time in state 0.

Fact:

displaymath319

Ergodic Theorem in continuous time.

Birth and Death Processes

Consider a population of X(t) individuals. Suppose in next time interval (t,t+h) probability of population increase of 1 (called a birth) is tex2html_wrap_inline325 and probability of decrease of 1 (death) is tex2html_wrap_inline327 .

Jargon: X is a birth and death process.

Special cases:

All tex2html_wrap_inline331 ; called a pure birth process.

All tex2html_wrap_inline333 (0 is absorbing): pure death process.

tex2html_wrap_inline335 and tex2html_wrap_inline337 is a linear birth and death process.

tex2html_wrap_inline339 , tex2html_wrap_inline341 : Poisson Process.

tex2html_wrap_inline343 and tex2html_wrap_inline337 is a linear birth and death process with immigration.

Queuing Theory

Ingredients of Queuing Problem:

1: Queue input process.

2: Number of servers

3: Queue discipline: first come first serve? last in first out? pre-emptive priorities?

4: Service time distribution.

Example: Imagine customers arriving at a facility at times of a Poisson Process N with rate tex2html_wrap_inline349 . This is the input process, denoted M (for Markov) in queuing literature.

Single server case:

Service distribution: exponential service times, rate tex2html_wrap_inline353 .

Queue discipline: first come first serve.

X(t) = number of customers in line at time t.

X is a Markov process called M/M/1 queue:

displaymath363

displaymath365

Example: tex2html_wrap_inline367 queue:

Customers arrive according to PP rate tex2html_wrap_inline349 . Each customer begins service immediately. X(t) is number being served at time t. X is a birth and death process with

displaymath377

and

displaymath379

Stationary Initial Distributions

We have seen that a stationary initial distribution tex2html_wrap_inline309 is a probability vector solving

displaymath315

Rewrite this as

displaymath385

Interpretation: LHS is rate at which process leaves state j; process is in state j a fraction tex2html_wrap_inline391 of time and then makes transition at rate tex2html_wrap_inline393 . RHS is total rate of arrival in state j. For each state tex2html_wrap_inline397 tex2html_wrap_inline399 is fraction of time spent in state i and then tex2html_wrap_inline403 the instantaneous rate of transition from i to j.

So equation says:

Rate of departure from j balances rate of arrival to j. This is called balance.

Application to birth and death processes:

Equation is

displaymath413

for tex2html_wrap_inline415 and

displaymath417

Notice that this permits the recursion

align197

which extends by induction to

displaymath419

Apply tex2html_wrap_inline421 to get

displaymath423

This gives the formula announced:

displaymath425

If

displaymath427

then we have defined a probability vector which solves

displaymath315

Since

displaymath431

we see that

displaymath433

so that tex2html_wrap_inline435 is constant. Put t=0 to discover that the constant is tex2html_wrap_inline309 .


next up previous



Richard Lockhart
Friday November 10 14:07:45 PST 2000