pij pik u pkj u k E
для любых u 0, .
Дополнительно мы будем предполагать, что выполняется следующее свойство:
5. (Непрерывность переходной функции в нуле)
lim p |
p |
0 для всех i, j E. |
0 ij |
ij |
|
Для цепей Маркова с непрерывным временем можно ввести матричную
переходную функцию (t) |
|
pij (t) |
|
|
, |
имеющую тот же смысл, что |
для |
|||||
|
|
|
|
|
|
N,N |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
||||
дискретных цепей Маркова имела матрица переходных вероятностей. |
|
|||||||||||
Интенсивности переходов |
|
|||||||||||
Теорема (без доказательства). |
|
|
Пусть t;t 0; - однородная |
цепь |
||||||||
Маркова с непрерывным временем, |
i, j E, |
i j . Тогда существует (но |
||||||||||
может быть бесконечным) предел lim |
|
pii t 1 |
, |
который мы обозначим , и |
||||||||
|
|
|||||||||||
|
|
|
|
t 0 |
t |
|
|
ii |
||||
существует конечный предел lim |
pij |
|
t |
, |
который мы обозначим . |
|
||||||
|
|
|
|
|||||||||
t 0 t |
|
|
ij |
|
||||||||
|
|
|
|
|||||||||
Определение. Величину ij назовем интенсивностью перехода из i-ого
состояния в |
j-ое; а величину ii |
назовем интенсивностью перехода из i-ого |
||||||||||||
состояния в i-ое. |
|
|
|
|
|
|||||||||
Матрица |
|
|
|
|
ij |
|
|
|
N,N |
называется инфинитезимальной матрицей, или |
||||
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|||||
матрицей интенсивностей. |
|
|
|
|
||||||||||
Утверждение. Пусть |
t;t 0; |
- |
однородная цепь |
Маркова с |
||||||||||
непрерывным |
временем. |
Тогда |
для |
ее |
интенсивностей |
выполняется |
||||||||
неравенство ij ii .
j i
Доказательство. Это неравенство имеет место не только в случае конечной цепи Маркова с непрерывным временем. Рассмотрим сумму
pij (t) 1,
j
pij(t) 1 pii (t).
j i
Для любого M такого, что M i, имеет место неравенство
M
pij (t) 1 pii (t)
|
j 1, j i |
1 pii t |
|
||
1 |
M |
|
|||
|
|
pij t |
|
. |
|
t |
t |
||||
j i |
|
||||
36
Устремим t 0 , тогда
M
ij ii . j i
M
Это неравенство для ij верно при любом M , следовательно, верно и
j i
для ij .
j i
Утверждение доказано. Следствие. В условиях утверждения для любого i E выполняется
неравенство
ij 0. j E
Определение. Состояние i E из фазового пространства однородной цепи Маркова с непрерывным временем называется мгновенным, если ii .
Состояние i E называется задерживающим, если ii 0. Состояние i E называется поглощающим, если ii 0.
Состояние i E называется регулярным, если оно задерживающее и
ij ii . j i
Определение. Цепь Маркова, все состояния которой регулярны, называется
консервативной.
Теорема (прямая система уравнений Колмогорова). Переходные вероятности консервативной цепи Маркова с конечным множеством состояний E, | E | N , удовлетворяют системе дифференциальных уравнений
dp (t) |
N |
1, |
i j |
|
ij |
pik (t) kj , по всем i, j E, с начальными условиями |
pij 0 |
, |
|
dt |
||||
k 1 |
0, |
i j |
i, j E.
Вматричном виде эта система записывается как
' t t ,
где ' t || pij '(t)|| - матрица, состоящая из производных переходных функций с начальными условиями 0 EN N - единичная матрица порядка N .
Доказательство. Зафиксируем h 0. Тогда для любых |
i, j E верно |
равенство |
|
N |
|
pij t h pij t pik t pkj h pij t |
|
k 1 |
|
pik t pkj h pij t pjj h pij t . |
|
k j |
|
Разделим левую и правую часть равенства на h. Тогда |
|
pij t h pij t |
pik t |
pkj h |
pij t |
pjj h 1 |
. |
|
h |
h |
|
||||
k j |
|
|
h |
|||
|
|
37 |
|
|
|
|
При h 0 предел правой части существует, значит, существует и предел
левой части, который мы обозначим через p '( ) t (производная |
p |
t |
справа). |
|||||||||||
|
|
|
|
|
ij |
|
|
|
|
ij |
|
|
|
|
|
pij '( ) t pik t kj pij t jj |
N |
|
|
|
|
||||||||
|
pik t kj . |
|
|
|
|
|||||||||
|
|
k j |
|
|
|
k 1 |
|
|
|
|
||||
Теперь выведем аналогичное уравнение для левой производной. |
|
|
|
|||||||||||
Для любых i, j E верно равенство |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
N |
|
|
|
|
|
|
|
|
|
|
pij t h pij t pij t h pik t h pkj (h) |
|
|
|
||||||||||
|
|
|
|
|
k 1 |
|
|
|
|
|
|
|
|
|
|
pij t h pij t h pjj h pik t h pkj h . |
|
|
|
||||||||||
|
|
|
|
|
k j |
|
|
|
|
|
|
|
|
|
Разделим левую и правую часть равенства на h. Тогда верны равенства |
||||||||||||||
|
pij t h pij t |
pij t h |
1 pjj h |
|
pik t h |
pkj h |
|
|
||||||
|
|
|
|
|
|
|||||||||
|
h |
|
h |
k j |
h |
|
|
|
||||||
|
pij t h |
pjj h 1 |
pik t h |
pkj h |
. |
|
|
|
|
|||||
|
|
|
|
|
|
|
||||||||
|
|
|
h |
|
k j |
|
h |
|
|
|
|
|||
Предел правой части существует при h 0 , |
значит, существует и предел |
|||||||||||||
левой части, который мы обозначим через p '( ) t (производная |
p |
t |
слева). |
|||||||||||
|
|
|
|
|
ij |
|
|
|
|
ij |
|
|
|
|
N
pij '( ) t pij t jj pik t kj pik t kj . k j k 1
Мы доказали, что
|
p '( ) t p '( ) t . |
||
|
ij |
ij |
|
|
dpij (t) |
N |
|
Следовательно, производная |
существует и равна pik t kj . |
||
dt |
|||
|
k 1 |
||
Теорема доказана.
Теорема (обратная система уравнений Колмогорова). Переходные вероятности консервативной цепи Маркова с конечным множеством состояний E, | E | N , удовлетворяют системе дифференциальных уравнений
dp (t) |
N |
1, |
i j |
|
ij |
ik pkj(t), по всем i, j E, с начальными условиями |
pij 0 |
, |
|
dt |
||||
k 1 |
0, |
i j |
i, j E.
Вматричном виде эта система записывается как
' t t
сначальными условиями 0 EN N .
Докажите эту теорему самостоятельно. Пусть нам дана некоторая матрица

ij 
N,N
со свойствами:
38
1. ij 0 и ii 0 при всех i j из E.
N
2. ij 0.
j 1
Рассмотрим вопрос, существует ли цепь Маркова, для которой произвольная матрица , удовлетворяющая свойствам 1 и 2, является матрицей интенсивностей.
Теорема (без доказательства). Если дана 
ij 
N,N , удовлетворяющая
сформулированным выше свойствам 1 и 2, то прямая и обратная системы уравнений Колмогорова имеют единственное решение
t 
pij t 
N N
со свойствами:
1. 0 pij t .
N
2. pij t 1.
j1
3.s t s t .
Отметим, что в некоторых пособиях по теории случайных процессов под системами дифференциальных уравнений Колмогорова понимают уравнения связывающие вероятности pk t P t k того, что консервативная цепь
Маркова в момент времени t находится в состоянии k с их производными по t. В матричном виде такие системы имеют вид:
p' t p t ,
где |
p t |
pk t ,k |
|
|
- распределение |
цепи |
Маркова с непрерывным |
|||||||||
1,N |
||||||||||||||||
|
|
|
|
|
|
|
|
|
dp |
t |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
временем в |
момент |
времени t, p' t |
|
k |
|
pk t ',k 1,N . Систему |
||||||||||
|
dt |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
N |
|
|
|
|
|
|
|||
решают при |
условии, |
что |
pk t 1. Подобные |
системы позволяют легко |
||||||||||||
k 1
находить стационарные распределения для цепей Маркова с непрерывным временем.
Процессы размножения и гибели
Определение. Консервативная цепь Маркова t;t 0; называется
процессом размножения и гибели, если для ее переходных функций выполняются следующие условия:
1). pi,i 1 t i t o t , i . 2). pi,i 1 t i t o t , i .
3). pi,i t 1 i i t o t , i .
39
1, |
i j |
4). pi,j 0 |
. |
0, |
i j |
5). 0 0, 0 0; для всех i 1 величины i и i больше нуля. Это пример процесса с бесконечным множеством состояний. Легко видеть, что верно равенство
pi,i 1 t pi,i 1 t pi,i t 1 o t ,
а матрица интенсивностей (бесконечного размера) будет иметь вид
|
|
0 |
0 |
0 |
0 |
|
|
|
1 |
( 1 1) |
1 |
0 |
|
|
|
0 |
2 |
( 2 2) |
2 |
|
|
|
|
|
|
|
|
Покажем, как будут выглядеть прямые и обратные системы уравнений
Колмогорова для такого процесса. |
|
' t t , как легко видеть, для всех |
|||||||||||||
Уравнения из прямой системы |
|||||||||||||||
i, j E будут иметь вид |
|
|
|
|
|
|
|
|
|
|
|
||||
|
dpi,j t |
pi,j 1 t j 1 pi,j t j |
j pi,j 1 t j 1. |
||||||||||||
|
|
|
|||||||||||||
|
|
dt |
|
0. |
|
|
|
|
|
|
|
|
|||
Мы будем считать, что |
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
Начальные условия стандартны: |
pi,j 0 |
1, |
i j |
|
|
||||||||||
|
|
|
. |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
0, |
i j |
|
|
||
Уравнения из обратной системы ' t t |
для всех i, j E будут иметь |
||||||||||||||
вид |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
dpi, j t |
p |
|
t |
|
p |
t p |
i 1,j |
t |
|||||
|
|
|
|
||||||||||||
|
|
dt |
i |
i 1, j |
|
i |
i |
|
i, j |
i |
|
||||
|
|
|
|
|
|
|
0,i j |
|
|
|
|||||
с теми же начальными условиями pi,j (0) |
. |
|
|
||||||||||||
|
|
|
|
|
|||||||||||
Обозначим pn t P t |
n , |
|
|
|
1,i j |
|
pn 0 P 0 n по |
||||||||
где |
n 0 . Тогда набор |
||||||||||||||
всем n 0 даст начальное распределение процесса размножения и гибели. Очевидно равенство
pn t pi 0 pi,n(t).
|
|
i 0 |
|
|
|
|
|
|
|
Определение. Вектор |
|
qo,q1,... , где |
qi 0, |
qi 1, называется |
q |
||||
|
|
|
|
i 0 |
стационарным распределением процесса размножения и гибели t;t 0; ,
|
|
если выполняется равенство qn qi pin (t) для всех n 0 |
и для любого t 0. |
i 0 |
|
40