Файл: Ганин, М. П. Прикладные методы теории марковских случайных процессов учебное пособие.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