Section 12 Continuous Time Markov Chains
Section 12 Continuous Time Markov Chains
Markov Chains
Professor Izhak Rubin
Electrical Engineering Department
UCLA
[email protected]
P X t y X s x , x,y S , t s 0
P x, y 1, t 0, x S
yS
3. Pt s x, y Pt x, z Ps z , y , x, y S , s,t 0.
zS
Chapman-Kolmogorov Equation
and Poisson Counting Process
Proof of Chapman Kolmogorov Equation:
Pt s x, y P X t s y X 0 x
P X t s y, X t z X 0 x
zS
P X t s y X t z, X 0 x P X t z X 0 x
zS
P X t s y X t z P X t z X 0 x Pt x, z Ps z , y
zS
zS
s,t 0, x,y S.
Example: Let N N t , t 0 be the Poisson Counting Process with intensity ,
over the state space of non-negative integers S = {0,1,2,...}.
Then, we readily show that:
1. N satisfies the Markov Property
2. Pt x, y e
t I y x 0 ,
y x !
yx
t 0, x,y S.
Probabilistic Structure
We consider a CTMC X with a standard TPF for which
lim Pt i, j ij , i, j S
t 0
Xt
Y3
Y1
Y0
A0 A1
Y4
Y2
A2
A3
A4
t
6
R i, j 1.
jS
Example :
For Poisson counting process N N t , t 0 with intensity ,
P An 1 An t | Yn i e t , t 0, qi 0, each i S.
For the embedded state sequence Y , the TPF is given as
R i, j j ,i 1 , i,j S; where j ,i 1 is the Kronecker Delta function, so that
R i, j qi
0, i j, q 0
i
Transition
Intensities
We define:
qij
d
1
Pt i, j lim [ Ph i, j P0 i, j ]; P0 i, j ij
t 0 dt
h 0 h
lim
so that
Ph (i, j ) qij h o(h); for i j; where lim [o(h)/h] = 0
h 0
and that
q
jS
ij
0.
Thus :
1
1
For i j: qij lim [ Ph i, j ] lim [ Avg.# transitions i j over (0, h)]
h 0 h
h 0 h
transition rate from state i to state j;
1
1
qi lim [1 Ph i, i ] lim [ Avg.# transitions out of state i over (0, h)]
h0 h
h0 h
transition rate out of state i.
Prof. Izhak Rubin
1
1
Pt h i, j Pt i, j Pt i, k Ph k , j kj
h
h
k
1
d
lim Pt h i, j Pt i, j Pt i, j Pt i, k qkj
h 0 h
dt
k
Prof. Izhak Rubin
10
jS
Using KFE,
d
lim Pt i, j 0
t dt
P k qkj 0, j S
kS
jS
(2.1)
k j
P j 1.
(2.2)
jS
11
CT Birth-and-Death Process
A continuous time process X X t , t 0 , S = {0,1,2,...}, is defined to be a
Birth and Death process if it is a CTMC with the following generator:
j i 1, i 0
i ,
,
j i 1, i 1
i
qij
i i , j i, i 0
0,
o.w.
where qi i i , 0 0; i 0, i 0; i 0, i 1.
The TPF for the embedded DTMC Y is thus given by:
i
, j i 1
qij i i
R i, j
qi i
, j i 1
i i
for qi 0 and i 1; R 0,1 1.
Prof. Izhak Rubin
12
P 0 0 P 1 1.
P j j j P j 1 j 1 P j 1 j 1.
P 0 0 P 1 1 0
P j j P j 1 j 1 P j 1 j 1 P j j 0
0
P
1
P 0
j 1
i
P j P 0
i 0 i 1
j 1
Let a j
i 0
i
for j 1, and set a0 1 ; then :
i 1
P j P 0 a j , j 0 .
Prof. Izhak Rubin
13
Case 1:
a
j 0
P 0 1
a
j 0
, P i ai
a
j 0
, i 0.
lim Pi X t j 0,
t
j 0.
14
qij i
i i , j i, N i 0
0,
o.w.
where qi i i , 0 0, N 0, and i 0, i 0, otherwise.
The TPF for the embedded DTMC Y is thus given by:
i
, j i 1
qij i i
R i, j
qi i
, j i 1
i i
for qi 0 and i 1, j N ; R 0,1 1, R N , N 1 1.
Prof. Izhak Rubin
15
Since
a
j 0
(assuming positive birth and rate transition intensities, except on the boundaries)
and is given by
P i P 0 ai ai
noting that
P 0 1
a
j 0
, i0
a
j 0
.
Prof. Izhak Rubin
16