Continuous Time Markov Chain (Markov Process)

Similar documents
5. Stochastic processes (1)

Lecture 4 Notes (Little s Theorem)

IS 709/809: Computational Methods in IS Research. Queueing Theory Introduction

Vehicle Arrival Models : Headway

Solutions to Assignment 1

Chapter 3 Common Families of Distributions

Basic definitions and relations

Problem Set 9 Due December, 7

Chapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws

Introduction to Probability and Statistics Slides 4 Chapter 4

PROBLEMS FOR MATH 162 If a problem is starred, all subproblems are due. If only subproblems are starred, only those are due. SLOPES OF TANGENT LINES

Homework 4 (Stats 620, Winter 2017) Due Tuesday Feb 14, in class Questions are derived from problems in Stochastic Processes by S. Ross.

Reliability of Technical Systems

Lecture 4 Kinetics of a particle Part 3: Impulse and Momentum

Stationary Distribution. Design and Analysis of Algorithms Andrei Bulatov

20. Applications of the Genetic-Drift Model

Solutions for Assignment 2

Appendix to Creating Work Breaks From Available Idleness

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Problem Set 5. Graduate Macro II, Spring 2017 The University of Notre Dame Professor Sims

Solutions Problem Set 3 Macro II (14.452)

TMA 4265 Stochastic Processes

MANAGEMENT SCIENCE doi /mnsc ec pp. ec1 ec20

Macroeconomic Theory Ph.D. Qualifying Examination Fall 2005 ANSWER EACH PART IN A SEPARATE BLUE BOOK. PART ONE: ANSWER IN BOOK 1 WEIGHT 1/3

Transient Little s Law for the First and Second Moments of G/M/1/N Queue Measures

Transform Techniques. Moment Generating Function

u(x) = e x 2 y + 2 ) Integrate and solve for x (1 + x)y + y = cos x Answer: Divide both sides by 1 + x and solve for y. y = x y + cos x

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

Stochastic models and their distributions

SOLUTIONS TO ECE 3084

APPM 2360 Homework Solutions, Due June 10

Asymptotic Equipartition Property - Seminar 3, part 1

Math 2214 Solution Test 1A Spring 2016

Embedded Systems 4. Petri nets. Introduced in 1962 by Carl Adam Petri in his PhD thesis. Different Types of Petri nets known

Statistical Distributions

Essential Microeconomics : OPTIMAL CONTROL 1. Consider the following class of optimization problems

EXERCISES FOR SECTION 1.5

4.5 Constant Acceleration

B. Maddah INDE 504 Simulation 09/02/17

CHAPTER 2 Signals And Spectra

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

Math 115 Final Exam December 14, 2017

2 Modern Stochastic Process Methods for Multi-state System Reliability Assessment

Intermediate Macroeconomics: Mid-term exam May 30 th, 2016 Makoto Saito

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)

Lecture Notes 5: Investment

N H. be the number of living fish outside area H, and let C be the cumulative catch of fish. The behavior of N H

Math 334 Test 1 KEY Spring 2010 Section: 001. Instructor: Scott Glasgow Dates: May 10 and 11.

Math Final Exam Solutions

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

Math 10B: Mock Mid II. April 13, 2016

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

1. Consider a pure-exchange economy with stochastic endowments. The state of the economy

Solution: b All the terms must have the dimension of acceleration. We see that, indeed, each term has the units of acceleration

Math 333 Problem Set #2 Solution 14 February 2003

Our main purpose in this section is to undertake an examination of the stock

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

The Arcsine Distribution

1 birth rate γ (number of births per time interval) 2 death rate δ proportional to size of population

EP2200 Queuing theory and teletraffic systems. 3rd lecture Markov chains Birth-death process - Poisson process. Viktoria Fodor KTH EES

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

KINEMATICS IN ONE DIMENSION

Chapter 8 The Complete Response of RL and RC Circuits

Comparison between the Discrete and Continuous Time Models

Math 334 Fall 2011 Homework 11 Solutions

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates

Fault Tolerant Computing CS 530 Reliability Analysis. Yashwant K. Malaiya Colorado State University

) were both constant and we brought them from under the integral.

10. State Space Methods

E β t log (C t ) + M t M t 1. = Y t + B t 1 P t. B t 0 (3) v t = P tc t M t Question 1. Find the FOC s for an optimum in the agent s problem.

- The whole joint distribution is independent of the date at which it is measured and depends only on the lag.

Double system parts optimization: static and dynamic model

Lecture Notes 2. The Hilbert Space Approach to Time Series

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB

Radical Expressions. Terminology: A radical will have the following; a radical sign, a radicand, and an index.

CHAPTER 12 DIRECT CURRENT CIRCUITS

Stochastic Model for Cancer Cell Growth through Single Forward Mutation

8. Basic RL and RC Circuits

EE650R: Reliability Physics of Nanoelectronic Devices Lecture 9:

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.

ON DETERMINATION OF SOME CHARACTERISTICS OF SEMI-MARKOV PROCESS FOR DIFFERENT DISTRIBUTIONS OF TRANSIENT PROBABILITIES ABSTRACT

13.3 Term structure models

Online Appendix for "Customer Recognition in. Experience versus Inspection Good Markets"

= ( ) ) or a system of differential equations with continuous parametrization (T = R

non -negative cone Population dynamics motivates the study of linear models whose coefficient matrices are non-negative or positive.

Problem Set on Differential Equations

Chapter 7 Response of First-order RL and RC Circuits

Number of modes per unit volume of the cavity per unit frequency interval is given by: Mode Density, N

Anno accademico 2006/2007. Davide Migliore

SMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15.

EXPONENTIAL PROBABILITY DISTRIBUTION

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Some Basic Information about M-S-D Systems

( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively:

FORECASTS GENERATING FOR ARCH-GARCH PROCESSES USING THE MATLAB PROCEDURES

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles

HOMEWORK # 2: MATH 211, SPRING Note: This is the last solution set where I will describe the MATLAB I used to make my pictures.

ME 391 Mechanical Engineering Analysis

Transcription:

Coninuous Time Markov Chain (Markov Process) The sae sace is a se of all non-negaive inegers The sysem can change is sae a any ime ( ) denoes he sae of he sysem a ime The random rocess ( ) forms a coninuous-ime Markov chain if P( ) ( ) e, ( ) e,, ( ) e P ( ) ( ) e n 2 2 n n n n n for any n, and for any ime sequence,,, and any as hisory e, e,, e n 2 n 2 Simly saed, he fuure deends on he resen sae only The as hisory is irrelevan Sae Residence Time Suose he sysem eners sae E and leaves he sae afer an amoun of ime is referred o as he sae residence ime Suose he sysem has sayed in sae E for a duraion of ime s Wha is he robabiliy ha he sysem will leave he sae E wihin ime? By he definiion of he Markove chain, he robabiliy does no deend on he as hisory, is memoryless P s s P Therefore, for any sae E, he sae residence ime is exonenially disribued

2 Theorem Memoryless Proery: roof P s s P for any s, is an exonenially disribued random variable P s and s P s s P s Tha imlies Ps s P s F ( s) F ( s) F ( ) F ( s) Noe F () Thus we have F ( s) F () s F () F () F () s F( s) F( s) F( ) F() F ( s ) aking he limi as, F () s F () F () s mus equal o P Define R ( s) F ( s) R s R R s R Then F ( s) R ( s) and F () R () Therefore we have ( ) () ( ) Define () Then R ( s) R ( s) Solving he above equaion wih R () F (), s R ( s) e s and finally F ( s) e

3 Sae Transiion Probabiliies and Raes in a Homogeneous Markov Process Assume P s k s does no vary wih ime s, and define he sae ransiion robabiliies k ( ) P s k s for any s Define he sae ransiion raes k k k ( ) lim is he rae a which he sysem changes is sae from sae E o E k Sae robabiliies Define () P ( ) is he robabiliy ha he sysem is in sae E a ime The seady-sae robabiliies are lim ( ) We will learn how o find he seady-sae robabiliies The seady-sae robabiliies are found from Global Balance Equaions,,2, ( e) k i i k i Eq e is referred o as he global balance equaion A flow is he sum of he corresonding sae ransiion raes The flow ou of he sae E is, while he flow ino he sae E is k i i k i The global balance equaion saes ha, for any sae, he flow ou of he sae mus equal he flow ino he sae

4 We mus relace one of he equaions in Eq e by 2 roof ( ) P i i P ip i () ( ) i i Dividing he sum ino wo erms, ( ) () ( ) () ( ) i i i Rearranging erms, () i( ) i() i( ) i i ( ) () () ( ) () ( ) i i i i i Dividing boh sides by, ( ) ( ) i ( ) i ( ) () i() ( e2) In he sae sae, we require ( ) ( ) lim As and, eq e2 becomes i i i i i i i

5 Markov Processs Examle: M/M/ queue wih one waiing roomm Cusomer Arrivals Cusomers arrive in a Poisson manner wih arrival rae The iner-arrival imes beween cusomers are exonenially disribued wih df a( x) e x x An arriving cusomer eners he server if f he server is idle If he server is busy, he cusomer eners he waiing room if here is a sace in he waiing room If he waiing room is full (already occuied by someone else inn his examle), hen he cusomer leaves he sysem No reservaion No reurn The cusomer is reeced Cusomer Service The service ime of each cusomer is a random variable exonenially disribued wih mean / The service ime has he df b( x) e x x Uon comleion of service, he server immediaely begins service of a cusomer waiing in he waiing room, if any If no cusomer waiing, he server becomes idle unil he nex cusomer arrives

6 Sae Sace Le ( ) denoe he number of cusomerss in he sysem a ime Since he waiing room can accommodaee a mos one cusomer, he sae sace is,, 2 A any ime, ( ), or 2 The sae sace of he sysem is comosed of hree saes Markov Process Le ( ) denoe he number of cusomerss in he sysem a ime Then () forms a coninuous-ime Markov chain since a sae change occurs if and only if eiher an arrivall occurs or he service of a cusomer ends Boh yes of he evens occurr in a memoryless manner The sae residencee imes are memoryless s Sae Transiion Raes Consider Sae ransiion occurs from sae o sae if an arrival occurss during, and Therefore The sae ransiionn rae is is he sae ransiionn robabiliyy from sae o sae in he inerval : P an arrival during from a Poisson source wih rae P s lim ( ) s for any s s

7 Nex consider Sae ransiion occurs from sae o sae when he service iss comleedd and no arrivals occur wihin is he sae ransiionn robabiliyy from sae o sae in he inerval : P service comleion wihin b x dx e P s s for any s P no arrivals during from a Poisson source wih w rae Therefore The sae ransiionn rae is lim ( ) Likewise, we can obain 2 2 2 2

8 Global Balance Equaion For sae, Flow Ou, Flow Ino The gbe is () e For sae, Flow Ou, Flow Ino 2 The gbe is 2 (2) e For sae 2, Flow Ou, Flow Ino The gbe is 2 2 (3) e Eq e-3 are linearly deenden Relace any one wih The soluion is 2 2 2 2 Sae Residence Times The sae residencee ime iss an exonenially disribued random variable For examle, he df of is e and he df of is, e

9 Birh deah Process and Local Balance Equaions A birh-deah rocess is secial case of he homogeneous Markov rocess in which for k k E E E k- E k E k+ The global balance equaion saes k k k k k k k I can be shown ha wih a birh-deah rocess, he flow equilibrium occurs and he local balance equaions are saified: and k k k k k k k k E k- E k E k+ Flow Equilibrium

B D Process Examle: M/M/ Queue Poisson Arrivals wih rae Infinie Waiing Room for Cusomers Exonenial Service Time wih mean M/M/ queue is a birh-deah rocess, and we can use he local balance equaions 2 3 4 o o o where 2 2 2 o o n n n nn o o 2 n n o 2 o o n n n n,,2,

Performance Measuress The fracion of ime he sysem is idle is Uilizaion, he fracion of ime he sysem is busy, is The average number of cusomers in he sysem is N, The average waiing ime, in FIFO, is n / n n W N The average sysem ime, he oal ime sen in he sysem, is T W Noe T N When a cusomer leaves he sysem, here are N cusomer rs waiing behind