Continuous Time Markov Chain

Similar documents
6. Stochastic processes (2)

6. Stochastic processes (2)

Continuous Time Markov Chains

Markov chains. Definition of a CTMC: [2, page 381] is a continuous time, discrete value random process such that for an infinitesimal

Google PageRank with Stochastic Matrix

Topics in Probability Theory and Stochastic Processes Steven R. Dunbar. Classes of States and Stationary Distributions

Applied Stochastic Processes

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

APPENDIX A Some Linear Algebra

Review: Discrete Event Random Processes. Hongwei Zhang

Random Walks on Digraphs

Introduction to Continuous-Time Markov Chains and Queueing Theory

EXPONENTIAL ERGODICITY FOR SINGLE-BIRTH PROCESSES

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

MARKOV CHAINS. 5. Recurrence and transience. and the length of the rth excursion to i by. Part IB Michaelmas 2009 YMS.

2.3 Nilpotent endomorphisms

MATH 281A: Homework #6

Perron Vectors of an Irreducible Nonnegative Interval Matrix

Lecture 3: Probability Distributions

Convexity preserving interpolation by splines of arbitrary degree

Lecture 2: Numerical Methods for Differentiations and Integrations

A note on almost sure behavior of randomly weighted sums of φ-mixing random variables with φ-mixing weights

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

DS-GA 1002 Lecture notes 5 Fall Random processes

Case Study of Markov Chains Ray-Knight Compactification

Lecture 21: Numerical methods for pricing American type derivatives

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Modelli Clamfim Equazioni differenziali 22 settembre 2016

Dynamic Systems on Graphs

Modelli Clamfim Equazioni differenziali 7 ottobre 2013

Engineering Risk Benefit Analysis

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS

Changing Topology and Communication Delays

FINITE-STATE MARKOV CHAINS

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

Chapter 7 Channel Capacity and Coding

FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

The Jacobsthal and Jacobsthal-Lucas Numbers via Square Roots of Matrices

Irreducibility. Irreducible. every state can be reached from every other state For any i,j, exist an m 0, such that. Absorbing state: p jj =1

Games of Threats. Elon Kohlberg Abraham Neyman. Working Paper

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation

Hidden Markov Models & The Multivariate Gaussian (10/26/04)

Lecture 17 : Stochastic Processes II

Entropy of Markov Information Sources and Capacity of Discrete Input Constrained Channels (from Immink, Coding Techniques for Digital Recorders)

Convergence of random processes

Introduction. - The Second Lyapunov Method. - The First Lyapunov Method

SUPER PRINCIPAL FIBER BUNDLE WITH SUPER ACTION

The Order Relation and Trace Inequalities for. Hermitian Operators

5 The Rational Canonical Form

Ballot Paths Avoiding Depth Zero Patterns

2 Finite difference basics

Affine and Riemannian Connections

TCOM 501: Networking Theory & Fundamentals. Lecture 7 February 25, 2003 Prof. Yannis A. Korilis

j) = 1 (note sigma notation) ii. Continuous random variable (e.g. Normal distribution) 1. density function: f ( x) 0 and f ( x) dx = 1

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Projective change between two Special (α, β)- Finsler Metrics

Chapter 7 Channel Capacity and Coding

AN EXTENDED CLASS OF TIME-CONTINUOUS BRANCHING PROCESSES. Rong-Rong Chen. ( University of Illinois at Urbana-Champaign)

Appendix B. Criterion of Riemann-Stieltjes Integrability

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

ECE 534: Elements of Information Theory. Solutions to Midterm Exam (Spring 2006)

THERE ARE INFINITELY MANY FIBONACCI COMPOSITES WITH PRIME SUBSCRIPTS

ISSN On error probability exponents of many hypotheses optimal testing illustrations

3 Basic boundary value problems for analytic function in the upper half plane

MEM 255 Introduction to Control Systems Review: Basics of Linear Algebra

MAT 578 Functional Analysis

CONDITIONS FOR INVARIANT SUBMANIFOLD OF A MANIFOLD WITH THE (ϕ, ξ, η, G)-STRUCTURE. Jovanka Nikić

Several generation methods of multinomial distributed random number Tian Lei 1, a,linxihe 1,b,Zhigang Zhang 1,c

2 More examples with details

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

The Feynman path integral

Modelli Clamfim Equazione del Calore Lezione ottobre 2014

MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS

Probability Theory (revisited)

ESTIMATION OF THE STATIONARY DISTRIBUTION OF A SEMI-MARKOV CHAIN

Perfect Competition and the Nash Bargaining Solution

Lecture 8 Markov Chains

Queueing Networks II Network Performance

Lectures on Stochastic Stability. Sergey FOSS. Heriot-Watt University. Lecture 5. Monotonicity and Saturation Rule

The lower and upper bounds on Perron root of nonnegative irreducible matrices

Hidden Markov Model Cheat Sheet

Problem Set 9 - Solutions Due: April 27, 2005

PHYS 705: Classical Mechanics. Canonical Transformation II

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

Three Kinds of Geometric Convergence for Markov Chains and the Spectral Gap Property

SL n (F ) Equals its Own Derived Group

APPROXIMATE PRICES OF BASKET AND ASIAN OPTIONS DUPONT OLIVIER. Premia 14

Expected Value and Variance

Numerical Heat and Mass Transfer

Tokyo Institute of Technology Periodic Sequencing Control over Multi Communication Channels with Packet Losses

Research Article Relative Smooth Topological Spaces

NUMERICAL DIFFERENTIATION

Lecture 12: Discrete Laplacian

TAIL PROBABILITIES OF RANDOMLY WEIGHTED SUMS OF RANDOM VARIABLES WITH DOMINATED VARIATION

More metrics on cartesian products

Transcription:

Contnuous Tme Markov Chan Hu Jn Department of Electroncs and Communcaton Engneerng Hanyang Unversty ERICA Campus

Contents Contnuous tme Markov Chan (CTMC) Propertes of sojourn tme Relatons Transton probablty functon Transton rate matrx Kolmogorov's forward and backward dfferental equatons Propertes of CTMC Irreducblty Recurrence Tme average Statonary dstrbuton Lmtng dstrbuton 2

Defnton Introducton A contnuous-tme stochastc process {X t, t } wth a countable state space S s sad to be a contnuous-tme Markov chan (CTMC) f t satsfes the Markov Property:,, P X j X X x u s P X j X ts s u u ts s for any, j, x S, u s, and all s, t. u The condtonal dstrbuton of any future state X t+s gven the past states and the present state X s s ndependent of the past states and depends only on the present state. Tme homogeneous (nvarant) CTMC: P{X t+s = j X s = } s ndependent of s. 3

Introducton Y T Y 2 T 2 Y 1 T 1 1 2 3 tme 4

Transton Probablty Functon Let P(t) = (p j (t)) where p j (t) = P{X t+s = j X s = }, whch s called the transton probablty functon from state to state j. The ntal dstrbuton and the matrx P(t) unquely determne the future behavor of the CTMC. For t t... t n, 1 1,, n t 1 1,, 1 1,, n n tn n tn n t 1 1 1 1,, n n tn n tn n n n1 t n1,, n1 n n1 p t t p t t p t PX n1 n 1 2 P X X X t n t n n n1 n P X X X P X X P X X P X X p t t P X X n1 n2 n1 n2 1 1 5

Notatons for a CTMC The tmes of jumps (state transtons): τ = < τ 1 < τ 2, The sojourn tmes: T n = τ n+1 - τ n, n The sequence of states vsted: Y = X, Y 1 = X τ1, Y 2 = X τ2, Y T Y 2 Y 1 T 1 T 2 1 2 3 tme We assume that the sample path s rght contnuous. 6

Propertes of the Sojourn Tme Based on the homogeneousness and the Markov property, we have, P T t s T s X P T t X Memoryless property of the sojourn tme n CTMC. Let g (t) = P{T > t X = }. Then we have PT t s T s, X PT s X condtonal probablty PT t X PT s X Markov property g t g s g t s P T t s X The only functon that satsfes the above equaton s g t e qt for some constant q. Note that g 1 The sojourn tme s exponentally dstrbuted. 7

Propertes of the State Transton For j, we have 1, 1, PY j X PT t X Markov property P Y j T t X P Y j T t X P T t X re where r j = P{Y 1 = j Y = }. j 1 qt Frst jump probablty matrx: R = (r j ) It s also the one step transton probablty matrx of the embedded DTMC {Y n } where the embedded ponts are the jump tmes. Note that the dagonal elements of R are all. 8

Instantaneous state A state for whch q =. Absorbng state A state for whch q =. A Regular CTMC A CTMC s sad to be regular (non-explosve) f, wth probablty 1, the number of transtons n any fnte length of tme s fnte. Example of a non-regular CTMC r q, 1 1, 2 Suffcent condton for a regular CTMC sup q S 9

Transton Rate Matrx Transton rate from to j : q r q j j j Defne q j rjq f j q f j Transton rate matrx: Q=(q j ) q q1 q2 q r1q r2q q1 q11 q12 r1q1 q1 r12q1 Q q2 q21 q 22 r2q2 r21q2 q 2 Note that we have Qe= where e s a column vector whose elements are all 1. 1

Propertes of q j ( j) The probablty of two or more transtons n tme t s o t. j h PY j, T h, T T h X oh p h P X j X h 1 1 qt,, P Y j T T h T t X q e dt o h 1 1, j Y P T h t T t X P Y 1 1 h q ht qt j j j h j e r q e dt o h qh j j e r q e dt o h j j h q q t r q r q e j qh q e q q q jh oh q r q h o h qh j q j o h o h qt q e dt o h 11 o h lttle o functon: lm h h

h PT h X oh qh e oh 1 q h oh Propertes of q The probablty of two or more transtons n tme t s o t. p h P X X 12

We have Propertes of q j 1 1 p h r q h o h q h o h j j j j p h q h o h q h o h q : transton rate from state to state j q j : transton rate of leavng state Therefore, we have lm h lm h p j h h p 1 h q j h j q 13

Basc Equatons for a CTMC Chapman - Kolmogorov's equaton p t s p t p s j k kj ks Kolmogorov's backward dfferental equaton d d d P t pj t qk pkj t,.e., t t. dt dt P QP ks dt Kolmogorov's forward dfferental equaton d d d P t pj t pk tqkj,.e., t t. dt dt P P Q ks dt 14

Kolmogorov's Forward Dfferental Equaton From Chapman-Kolmogorov's equaton, we have p t s p t p s j k kj ks Consequently, j j k kj j k kj jj j ks ks, k j. p t s p t p t p s p t p t p s 1 p s p t. Hence, s p 1 j t s pj t pkj s p jj lm lm pk t p s s s ks, k j s s d dt t t. P P Q ks, k j ks, k j ks p k s p pk tlm p t s kj j s s p t q p t q k kj j jj t q kj. j 1 p lm s jj t s 15

Example: The Two-State Chan A CTMC spends an exponental tme wth rate λ n sate before gong to state 1 It spends an exponental tme wth rate μ n state 1 before gong to state 1. Q Compute p (t) and p 11 (t) 1 16

Example: The Two-State Chan From Kolmogorov's forward equaton, we have p t p t p t p ' t p t p t 1 Note that 1 1 e p ' t p t e t t d t e p t e t dt t t e p t e c 17

Example: The Two-State Chan t p t ce p 1 p t e In addton, we can observe lm p t. t By symmetry, we also have p11 t e t t 18

Tme Dependent Probablty From Kolmogorov's forward dfferental equaton, the vector P(t) s gven by P(t) = P() e Qt = e Qt (Note that P() = I), where the matrx exponent functon e Q s defned by e Q = I + Q + Q 2 /(2!) + = n= Q n /(n!) 19

Irreducblty and Recurrence Irreducblty A CTMC {X t } s rreducble f any one n the followng holds: The embedded DTMC {Y n } s rreducble For any,j S, we have p j (t) > for some t >. For any,j S, we have p j (t) > for all t >. Recurrence Let nf t X, lm t X s st A state s recurrent (transent) f 1 1 P X 2

Statonary Measure = (, 1, )( ) s a statonary measure f and P(t) = for all t. Suppose that a CTMC {X t } s rreducble and recurrent. Then there exsts the unque statonary measure up to a constant multplcaton. 1., for all S 2. For any state, f we defne E I X s j ds X then ther exsts aconstant c such that π= cπ. Note that j E X j 3. j where q j j π j s the statonary measure for the embedded DTMC Yn js 21

Statonary Dstrbuton The statonary measure s, n fact, a soluton of Q =. Proof: From the Kolmogorov's forward dfferental equaton, we have πp d t π t t dt πp πp Q πq The thrd property n the prevous slde From πq, for any S we have q q q q q r j j j j j j j js js, j js q s a statonary measure of the embedded DTMC Y S n whose transton probablty matrx R= r wth r j. q q j j js, j rate out = rate n Global balance equaton 22

Postve Recurrence A recurrent state s called postve recurrent (null recurrent) f E[τ() X = ] < (= ). Therefore, an rreducble and postve recurrent CTMC has the unque statonary dstrbuton: 1 E X 23

Ergodcty For a CTMC, we don't need to consder the perodcty because all the sojourn tmes are exponentally dstrbuted. An rreducble postve recurrent CTMC s called ergodc. If an rreducble regular CTMC {X t } has a nonnegatve numbers satsfyng Q = and e = 1, then {X t } s ergodc and hence t has the unque statonary dstrbuton. If a CTMC {X t } s ergodc and s the statonary dstrbuton, then for all, j S we have lm p j t j t 24