Файл: Ганин, М. П. Прикладные методы теории марковских случайных процессов учебное пособие.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 17.10.2024
Просмотров: 172
Скачиваний: 0
Умножив уравнение для вероятности Pkj(/0, /) из (5.6) на и) и просуммировав результат по всем возможным значениям /, при ходим к равенству
^ P k](t0, t) |
|
|
O |
00 |
|
= |
- T (t) |
2 uipk, (t0, |
j=k+1 |
о |
|
J'-k |
|
|
_j=k |
|
|
С учетом (5.8) |
это выражение переписывается в виде |
|
|||
dGk (u; |
t0, |
t) |
l )Gk(a; |
t0, t). |
(5.10) |
dt |
|
= |
|||
|
|
|
|
|
В данном уравнении переменная и является параметром. Поэтому равенство (5.10) представляет собой обыкновенное дифференциаль ное уравнение первого порядка с разделяющимися переменными.
Интегрируя это уравнение, |
находим |
|
|
|
|
|||
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
(u — 1) |
( 1 (Е) dE |
|
|
Gk(u; t0, |
t) = |
Gk(ir, tm t0)e |
. |
(5.11) |
|||
Положим |
|
|
|
|
|
|
|
|
|
a(t0, |
t ) = |
|
|
|
|
(5.12) |
|
|
|
|
|
to |
|
|
|
|
Тогда с учетом (5.9) из (5.11) получаем |
|
|
|
|||||
Имеем |
Ок(и; t0, 0 |
= « k«(u_I),(t,,t)- |
(5.13) |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
eua |
2 |
(ua)s |
|
|
|
|
|
|
= |
|
|
|
|
|
|
|
|
|
s=0 s\ |
|
|
||
Поэтому |
разложение |
функции |
(5.13) |
в ряд по степеням |
пара |
|||
метра и будет следующим: |
|
|
|
|
|
|
||
|
Gk (и; t0, t) = |
|
|
|
|
= |
|
|
|
|
а(1о, 1) у |
и> \ajto, |
Q ] j - |
k |
(5.14) |
||
Сравнивая коэффициенты при одинаковых степенях параметра и |
||||||||
в (5.8) и |
(5.14), для искомых вероятностей перехода Р^Л^о, |
0 по |
||||||
лучаем |
|
|
[a (t0, |
t)V~k |
|
|
||
|
Лд (to, |
|
|
(5.15) |
||||
|
t) — |
|
fe)l |
|
|
|||
|
|
|
(/' ~ |
|
|
|
||
|
(i = |
k, k + l , |
. . . ; |
k = |
0, 1, |
...). |
|
31
Так как /0 и t произвольные, то вероятность перехода при т ~>t определяется формулой
|
Ш , *)Р~к |
kj(r’ |
Ь —k)\ |
J / - * ! * |
|
|
{К j = 0, 1, . .. ; / > A). . |
Ру (Л. т) I
(5.16)
Обратная система дифференциальных уравнений Колмогорова (4.9) для пуассоновского процесса имеет вид:
dPkk(t, -■)
dt |
T(*)^kk(*. |
т); |
|
|
|
|
|
|
|
|
|
|
— Т V) |
\Pki(t, |
r ) - |
Pk+U(t, т)1 |
(5.17) |
(/ = * + 1 , k + |
2, . .. ; |
k = |
0, 1, ...). |
|
О помощью непосредственной подстановки нетрудно убедиться, что решение этой системы также определяется формулой (5.16).
Система уравнений (4.6) для вероятностей Рj (t) при пуассонов ском процессе записывается в виде:
' W = - 7 |
( * ) P 0(*); |
j |
Pj (0 = - |
Т(*) [Pj (0 — Pj-i(01 |
(5.18) |
|
(/ = 1 , 2 , . . . ) . |
I |
При t = 0 физическая система, соответствующая рассматривае мому процессу, находится в состоянии С0. Поэтому начальные усло вия для искомых функций Рj (t) следующие:
Р0(0) = 1; Pj (0) = 0 (У=1,2,...). |
(5.19) |
Решение системы (5.18) может быть получено по аналогии с ре шением системы (5.6). При этом Pj (t) совпадает с Ру(Аь 0 при k = 0 и to= 0. Следовательно, вероятность того, что в момент t фи зическая система находится в состоянии Cj, определяется формулой
= |
e-»(0,t) 0 = 0 , 1 , . . . ) , |
(5.20) |
где |
t |
|
|
|
|
fl(0, |
*) = J 7 (5) |
(5.21) |
|
о |
|
Математическое ожидание x(t), дисперсия D[X(t)] и корреля ционная функция Kx(t, %) при %> t для пуассоновского процесса
равны а(0, t), т. е. (см. пример 2.1) |
|
x { t ) = D [X (0] = Кх {t, т) = а (0, t). |
(5.22) |
32
Обозначим через P(k-\-s, x/k, t) вероятность того, что за про межуток времени от t до т физическая система из состояния Ск пе рейдет в Ck+S. Эта вероятность совпадает с Pk,k+s (t, т), п потому согласно (5.16)
Р [k + s, x/k, t) = ^ |
е - а (*• |
(5.23) |
(s = 0, 1, ...)•
Из выражения (5.23) следует, что за любой промежуток вре мени от t до т случайное число смен состояний системы при пуассонс/вском процессе подчиняется распределению Пуассона ■с пара метром a{t, т). Данный параметр равен математическому ожиданию (и дисперсии) числа смен состояний физической системы за ука занный промежуток времени. Математическое ожидание числа смен состояний в единицу времени в интервале от t до т будет
т
X {t, х) = |
т (Ц) dl |
(5.24) |
t
Мгновенная интенсивность (плотность) смены состояний физиче ской системы при пуассоновском процессе получается равной
Х(£) = НтН£, т ) = |
Пт |
Т $ ) d l |
T~t |
--t |
|
Раскрывая неопределенность по правилу Лопиталя, находим
М 0 = Г (0 - |
(5.25) |
Таким образом, интенсивность (плотность) смены состояний X(t) физической системы при пуассоновском случайном процессе в лю бой момент времени t совпадает с временной плотностью вероятно сти у (0 перехода системы в следующее состояние.
При однородном пуассоновском процессе интенсивность смены состояний К, совпадающая с у, является постоянной величиной. Данная постоянная равна математическому ожиданию числа смен состояний физической системы в единицу времени. Заменяя в (5.5) Pkj {to, t) на Pkj(0 и у на Я,, прямую систему дифференциальных уравнений Колмогорова при однородном пуассоновском процессе получаем в виде:
Pkk (^) — ~~ ^Лек if) !
P k j ( 0 = - X [ P k k W - P k , , - i W ] |
' (5.26) |
{i = k + 1, k + 2, . . . ; k = 0, 1, ...).
3 |
33 |
Согласно (5.15) при a.(to, t) — Kt решение этой системы
|
|
W |
(иу~* |
-xt |
(5.27) |
|
|
(j — k)\ |
|||
|
|
|
|
|
|
(/ = |
k + 1, k -J- 2, |
. . . ; k — Q, 1, |
...). |
||
Обратная система |
уравнений |
Колмогорова |
(5.17) при однород |
||
ном процессе будет |
|
|
|
|
|
а д |
= |
- х/экк (г); |
|
|
|
Рк) (0 ~ |
^kj (О + |
^Ли-I.j (О |
|
(5.28) |
(i = k + 1) k + 2, . . . ; & = 0, 1, ...).
Система уравнений (5.18) для вероятностей Рj (t) при однород ном процессе записывается так же, как при неоднородном процессе, т. е.
P'0(t) = - l P 0(t);
PJ (t) = - l P ]{t) + lP^(t) |
(5.29) |
(/= 1, 2, ...). |
|
Уравнения (5.29) относительно просто можно получить без ис пользования общих систе!м дифференциальных уравнений Колмого рова. Для этого только нужно с помощью формулы полной вероят ности найти вероятность Pj (t -f- Дt) того, что в момент t -)- At фи зическая система находится в состоянии Cj (/ = 0, 1, ...), когда гипотезой Hs является факт нахождения системы в момент t в со стоянии Cs (s — 0, 1, ...). За малый промежуток времени от t до t -j- At физическая система не может изменить более одного состоя ния. Поэтому данная система может оказаться в состоянии Ci в мо мент t -f- At только в двух случаях, когда за время At исходное со стояние Cj не изменяется и когда за это время происходит переход из состояния Сj_i в Cj. Следовательно,
Р<> “Ь — Ро (t)Лю {С t+ &.t);
Pi{t+U) = Pj (0 Pjj {t, t+ \t) + P j_1(t) Pj-i.j (t, t + \t) (5.30)
(/ = 1, 2, ...).
Введем функцию 0(Д/), которая по отношению к At имеет поря док малости выше первого, а потому
|
lim |
0Щ ) |
= 0. |
(5.31) |
|
it -.0 |
м |
|
|
Тогда |
14" Д£) — ^kt -f- 0 (&t); |
| |
||
Pj— |
||||
Pjj (t, t + Д*) = 1 - |
Ш + 0 (\t). |
(5.32) |
||
j |
34
Подставляя эти выражения в (5.30), получаем |
|
|
||||||
P0(t + M) = |
P0(t)(l - Ш ) + 0 ( Д * ) ; |
|
|
(5.33) |
||||
Pi(t + |
b t ) = ; P i ( t ) ( l - ХМ) + |
(t) Ш + |
О (М) |
|||||
или |
|
(7 = |
1 , 2 , . . . ) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
P0(t + |
M ) - P 0(t) |
|
, п m , |
0 (ДО. |
|
|
||
' |
It |
~ |
“ |
aF° {t) + |
~м~ ’ |
|
|
|
|
|
( 0 _ = |
_ |
ХР. {t) + |
хр._1(0 + |
9 W ) |
(5.34) |
|
|
|
(/ = 1 ,2 , . . . ) . |
|
|
|
|
||
Переходя в (5.34) к пределу при At -> 0, |
получим систему диф |
|||||||
ференциальных уравнений (5.29). |
|
|
|
|
||||
Решение системы (5.29) |
может быть найдено, как и выше, с по |
|||||||
мощью производящей функции. При этом |
Pj (t) ~ Poj (t). |
Чтобы |
||||||
произвести другое решение этой системы, положим |
|
|||||||
|
|
Pi(t) = |
e~u Qi (t). |
|
|
(5.35) |
||
Так как |
|
|
|
|
|
|
|
|
|
|
р\(0 = [QJ W —xQj(01 e~xt, |
|
|
то система (5.29) преобразуется к виду
|
Qi(f) = |
0, |
Q;.(0 = |
^Qj-i(0 |
(/ = |
1 , 2 , . . . ) , |
(5.36) |
|
причем |
|
|
|
|
|
|
|
|
|
Qo(0) = |
1, |
Qj (0) = |
0 (/ = |
1, 2, |
. |
. . ). |
(5.37) |
Из (5.36) последовательно находим: |
Q0(7) = l; |
Ql(t)=%t-J |
||||||
0.2{t) — |
• Полагая |
|
|
|
|
|
||
|
|
|
|
(W )i-1 |
|
|
|
|
|
|
|
Ql-i(0 = |
|
|
|
|
|
|
|
|
|
( / - i ) i ’ |
|
|
|
|
получаем |
|
|
|
(Ity |
|
|
|
|
|
|
|
Qi(t) |
|
|
(5.38) |
||
|
|
|
}• |
|
|
|
||
|
|
|
|
|
|
|
|
|
Следовательно, |
решение системы (5.29) |
записывается в виде |
||||||
|
PiV) = м |
/! |
(/ — 0, |
1, . |
. |
. ). |
(5.39) |
|
|
|
|
|
|
|
|
|
;35