<< Chapter < Page Chapter >> Page >

The following important conditions are intuitive and may be established rigorously:

  • i j implies i is recurrent iff j is recurrent
  • i j and i recurrent implies i j
  • i j and i recurrent implies j recurrent

Limit theorems for finite state space sequences

The following propositions may be established for Markov sequences with finite state space:

  • There are no null states, and not all states are transient.
  • If a class of states is irreducible (i.e.,has no proper closed subsets), then
    • All states are recurrent
    • All states are aperiodic or all are periodic with the same period.
    • If a class C is closed, irreducible, and i is a transient state (necessarily not in C ),
      • then F ( i , j ) = F ( i , k ) for all j , k C .

A limit theorem

If the states in a Markov chain are ergodic (i.e., positive, recurrent, aperiodic), then

lim n p n ( i , j ) = π j > 0 j = 1 M π j = 1 π j = i = 1 M π i p ( i , j )

If, as above, we let

π ( n ) = [ p 1 ( n ) p 1 ( n ) p M ( n ) ] so that π ( n ) = π ( 0 ) P n

the result above may be written

π ( n ) = π ( 0 ) P n π ( 0 ) P 0

where

P 0 = π 1 π 2 π m π 1 π 2 π m π 1 π 2 π m

Each row of P 0 = lim n P n is the long run distribution π = lim n π ( n ) .

Definition . A distribution is stationary iff

π = π P

The result above may be stated by saying that the long-run distribution is the stationary distribution. A generating function analysisshows the convergence is exponential in the following sense

| P n - P 0 | a | λ | n

where | λ | is the largest absolute value of the eigenvalues for P other than λ = 1 .

The long run distribution for the inventory example

We use MATLAB to check the eigenvalues for the transition probability P and to obtain increasing powers of P . The convergence process is readily evident.

P = 0.0803 0.1839 0.3679 0.36790.6321 0.3679 0 0 0.2642 0.3679 0.3679 00.0803 0.1839 0.3679 0.3679 E = abs(eig(P))E = 1.00000.2602 0.26020.0000 format longN = E(2).^[4 8 12] N = 0.00458242348096 0.00002099860496 0.00000009622450>>P4 = P^4 P4 =0.28958568915950 0.28593792666752 0.26059678211310 0.16387960205989 0.28156644866011 0.28479107531968 0.26746979455342 0.166172681466790.28385952806702 0.28250048636032 0.26288737107246 0.17075261450021 0.28958568915950 0.28593792666752 0.26059678211310 0.16387960205989>>P8 = P^8 P8 =0.28580046500309 0.28471421248816 0.26315895715219 0.16632636535655 0.28577030590344 0.28469190218618 0.26316681807503 0.166370973835350.28581491438224 0.28471028095839 0.26314057837998 0.16633422627939 0.28580046500309 0.28471421248816 0.26315895715219 0.16632636535655>>P12 = P^12 P12 =0.28579560683438 0.28470680858266 0.26315641543927 0.16634116914369 0.28579574073314 0.28470680714781 0.26315628010643 0.166341172012610.28579574360207 0.28470687626748 0.26315634631961 0.16634103381085 0.28579560683438 0.28470680858266 0.26315641543927 0.16634116914369>>error4 = max(max(abs(P^16 - P4))) % Use P^16 for P_0 error4 = 0.00441148012334% Compare with 0.0045824...>>error8 = max(max(abs(P^16 - P8))) error8 = 2.984007206519035e-05 % Compare with 0.00002099>>error12 = max(max(abs(P^16 - P12))) error12 = 1.005660185959822e-07 % Compare with 0.00000009622450

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Applied probability. OpenStax CNX. Aug 31, 2009 Download for free at http://cnx.org/content/col10708/1.6
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Applied probability' conversation and receive update notifications?

Ask