Файл: Ганин, М. П. Прикладные методы теории марковских случайных процессов учебное пособие.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 17.10.2024

Просмотров: 163

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Р(А) — РУ} (t, т).

Гипотеза / / г пусть означает,

что в момент

где

случайная функция равна х,, если

Х{1)— Ху; при

атом

Р{Н{) = Р ЫЬ, /'),

а

Р ( А / Н Р ц (/',

т).

Следовательно,

спра­

ведливо равенство

 

 

 

 

 

 

(t, 0

=

2

ры (t, Г) Рп (Р, т)

(t < Г < Т)

(2.13)

 

 

/=0

 

 

 

(k, j — 0, 1........ 7V)'

Данное соотношение между вероятностями перехода называется

уравнением

Колмогорова — Чепмена — Смолуховского.

что

Функция

Ру,j (4, 0 является условной вероятностью того,

в момент /

справедливо равенство X ( t ) = x j, определенной

при

условии, что X(t0) -=Ху. Если действительно в начальный момент tQ

случайная

функция

Х{() равна

xkl то

условная

вероятность

Л ] (Аь I)

совпадает с

безусловной

вероятностью Р-} (/) того, что

в момент t случайная функция X(I)

равна

х г Когда

известны ве­

роятности Ру(1о) того» что в момент /о случайная функция X(t)

равна x k (к — 0, 1, ...,

N), для безусловной вероятности Л (0 =

— Р[А'(/)— Xj] согласно

формуле полной вероятности

(2.12) полу­

чаем

 

 

Л ( 0 = 2 Л ( * о ) Л с 1 ( * , » t).

(2.14)

 

к=0

 

При фиксированном t совокупность вероятностей P\{t) п воз­ можных значений х, (j — 0, 1, ..., N) является рядом распределе­ ния случайной величины X(t). С помощью этого ряда распределения можно найти числовые характеристики случайного процесса X(i).

В частности, математическое ожидание x(t) и дисперсия Z)[2l(/)] случайной функции X(t) определяются формулами:

 

_

N

 

 

 

(2.15)

 

* ( 0 =

2 V j ( 0 ;

 

 

 

j=o

 

 

 

 

[АГ(01 = £ [x-j -

х (0 ]2Л (0 =

2

x]Pi(t) — Lx (г1)]2 .

(2.16)

j= o

 

 

j= 0

 

 

 

Если марковский случайный процесс X (/)

целочисленный,

т. е.

его возможные значения

xs — j

(/ = 0,

1,

.. .,

N), то при определе­

нии математического ожидания и дисперсии удобно пспользовать производящую функцию G(u\ t), которая зависит от вещественного аргумента и и от времени t. Данная функция определяется фор­ мулой

N

G (и; t) =

2 и1A (t),

(2.17)

так что

j= 0

 

 

(2.18)

0(1;

0 = 1 -

13


Через производные по и от производящей функции математиче­ ское ожидание и дисперсия целочисленного марковского случайного процесса X (/) выражаются равенствами:

X(t):

dG{it; t)

)

(2.19)

ди

 

11=1

 

d2G («;

t)

X(t)— \x(t)\2.

(2.20)

D\X(t) ] =

+

ди2

U -1

 

 

Кроме математического ожидания и дисперсии случайного про­ цесса X(t), в приложениях часто используется корреляционная функция Кх (/, т). Для любого вещественного случайного процесса X (/) корреляционная функция определяется формулой

K A t, т) = Ж {[AT (^) — ^ (01 [ЛГ(х) — ^ (т)]}.

(2.21)

Используя вероятности (2.1), корреляционную функцию марков­ ского случайного процесса с дискретными ординатами и непрерыв­ ным временем можно рассчитать с помощью равенства

к х (t, т )= V V

[

^ - ^ [ x j - * ( - ) ] Rki (/, А .

(2.22)

k=0 j —О

 

 

 

 

Данную формулу можно также переписать в виде

 

Кх (/, т ) = 2

2

XkXftk) (t,

х) - х (t) х (х).

(2.23)

k=0 j = 0

 

 

Дисперсия случайного процесса X(t)

связана с корреляционной

функцией соотношением

 

 

 

 

D[X(t)] = Kx(t,t).

 

(2.24)

Случайный процесс X (/)

с дискретными ординатами и непрерыв­

ным временем называется стационарным, если при любом t0> 0 и при произвольных И < k < . . . < t u где I — любое число, совпадают /-мерные ряды распределения систем случайных величин X(tx),

X(t2), ... , X(tt) и X(tx-\- /о), X(t2Jr t0), ... , X{tt +/o)-

Условие стационарности означает вероятностное равновесие слу­ чайного процесса во времени, т. е. что ряды распределения любого

порядка / для X(t)

не зависят от

начала

отсчета

аргумента /.

Дйя такого процесса математическое ожидание

x(t) ц

дис­

персия D[X(/)] постоянные, а корреляционная функция K x(t,

т) за­

висит от t и т только как от разности, т. е.

 

 

 

х (*) = const,

D [ X (*)] = const,

Кх(t,

t) = АГх(х — /).

(2.25)

14


В приложениях часто не требуется знать ряды распределения любого порядка, а достаточно иметь математическое ожидание и корреляционную функцию случайного процесса X(t). Если мате­

матическое ожидание х постоянно, а корреляционная функция (t, т) зависит от t и % только как от разности, то такой случай­ ный процесс X(t) называется стационарным в широком смысле. Любой стационарный процесс является стационарным в широком

смысле. Обратное не всегда верно, так как из равенств

(2.25) не

следует выполнение условий стационарности процесса.

 

с

Пример 2.1. Целочисленный марковский случайный.процесс X(t)

дискретными ординатами и непрерывным временем

совпадает

с

числом требований в системе массового обслуживания, которые

в систему поступают по одному независимо друг от друга. Вероят­

ность P kJ-

(t, т) того, что к моменту т >

/ в систему поступит / тре­

бований,

если к моменту t

их поступило k, выражается формулой

 

 

(*. -О =

е

~ н и т) (/ > k)

( 2 -26)

где а(/,

т) — заданная функция, которую можно представить в виде

 

 

a{t,

т) —

.

 

 

 

 

t

 

 

При t =

О'требований в системе нет.

G(и, t). Определить матема­

Составить производящую функцию

тическое ожидание, дисперсию и корреляционную функцию про­ цесса X(t). Составить двумерный ряд распределения

tfkj V, x ) = P [ X ( t ) = k , X (x )= f]

{k, j = 0, 1, ...).

Р е ш е н и е . При фиксированном t случайная величина X (t) совпадает с числом требований в системе массового обслуживания. Запишем ряд распределения этой случайной величины в виде

P l ( t ) = P [ X ( t ) = j ]

 

 

(/ = 0 , 1 , . . . ) .

 

 

 

Так как при t = 0 в

системе требований нет, то

Р о ( 0 ) = 1 ;

Pj (0) — 0 ( / =

1, 2, ...).

Воспользовавшись формулой

(2.14) при

^о = 0, N== со,

для вероятности Pj(t)

при любом t находим

 

Л ( 0 = 2 Л < ( о ) Р к](о,

t) =

poi(о,

t),

 

т. е.

к=0

 

 

 

 

 

 

 

 

 

 

 

р. (t) =

e-a(0, t)

(/ =

0, 1,

...).

 

15


Нз полученного выражения следует, что случайная величина X(t) распределена по закону Пуассона с, параметром а (0, /). Про­ изводящая функция

 

 

 

[иа(0, ^)]j

 

 

 

Л

 

. = exp \{и -

1) а (0,

t)\.

Так как

 

 

 

dG

а (0, t) G (и; t)\

д'Ч)

 

да

ди~ =

[а (0, t)YG(u\ t) ,

то согласно (2.19) и (2.20)

x{t) = D\X{t)] = a{<d, t).

При т 'У t математическое

ожидание случайного процесса X(t)

будет

 

 

 

х W = а (0, х) = Г

Г т (5)

-h

( & ) Я ,

0

о

 

t

т. е.

 

 

 

-) .

 

 

Чтобы найти корреляционную функцию ставим Х(х) в виде Х(Г) — X(t) -)-£(/, х), бований, поступивших в систему за время

Кх {t, т)

при т >

t,

пред­

где £(f,

т) — число

тре­

от t до т. Тогда

 

.

Кх (t, x) = M [ [ X ( t ) - x (*)] [X (Т) _ *

(т)]} =

= M{[X(t)- x(t)\*} + M \ { X { t ) - x ( t ) ] \ l { t ,

x ) - a ( t , т)]) .

Случайные величины X(t) и £(/, т) независимые, поэтому по­ следнее слагаемое равно нулю. Следовательно, в данном случае кор­

реляционная

функция совпадает

с дисперсией,

т.

е.

Кх (t,

т) =

= £>[*(/)] =

а (0,

t).

Так

как Кх{х,

t) -

К* (t,

х),

то Kx{t,

т) =

— а (0, т) при t >

т.

 

 

 

 

 

 

 

 

 

При j < k и т

1вероятность Rki (t,

т)

равна нулю,

а при j k

Rki(t> i) = p k(t)

РыУ,

[a(0, t)\*[a{t, T)]i-k

 

 

x) =

 

kl(j .- k)l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, искомый двумерный ряд распределения

 

 

 

 

 

0

при

/ < k\

 

 

 

P [ X ( t ) = k ,

X(x) =

j] =

И О .

t)\*[a(t,

т)р . £—а(0,т) при

/> £ .

 

 

 

 

 

6 ! (/ — &)!

 

 

 

 

16


Если t > т, то

 

 

 

 

О

при

k < j \

P \ X { t ) = k , X ( x ) = j ] :

[я (0, *)]кИ * .

t)}*-'

a(0,t) при

 

М Ь ~ У)!

 

Воспользовавшись формулой (2.23), находим

 

 

оо

со

 

 

 

 

 

 

KAt, *)= 2

2

[А- * ( * ) ] [ / - * to]як1(*>

*0 =

 

 

к — 0 j=U

 

 

 

 

 

 

: ^ -а(0.т) у

1

[£_*(*)]

[а (0’

i ;

[S-f £-Зф)]

^

 

 

к !

' 1

 

 

 

s!

j-o

 

 

 

s=0

 

 

 

Так как

 

 

 

 

 

 

 

 

 

 

 

[ a ( 0 ,

О Г

~ 0 ;

2 -

[ a t t ,

t ) | s

oa(t,‘ )

 

 

 

 

2 [ * - * ( * ) ]

k\

s!

 

 

k= 0

 

 

 

 

 

s=0

 

 

 

TO

 

 

 

 

 

 

 

 

 

tf*(*.

T) =

V

[ ^

_ 1) +

k -

(0]

HQ, t)]*

 

 

k — 0

 

 

 

 

 

k\

 

 

 

 

 

 

 

 

=

[a (0, ^)]2 +

а (0, £) — x (£) a (0,

£) =

a (0,

t),

что совпадает с полученным ранее выражением.

§ 3. ОДНОРОДНЫЙ МАРКОВСКИЙ ПРОЦЕСС

Марковский случайный процесс X (t) с дискретными ординатами и непрерывным временем называется однородным, если определяе­

мые формулой (2.3) вероятности перехода Pkj ( t, т)

зависят от t и

т как от разности, т. е. если

 

 

 

^kj (^, х) =

^kj (х —

( ~ > t )

(3.1)

 

( * , / =

0, 1, . .. , N).

 

В этом случае условные вероятности (3.1)

зависят от одного аргу­

мента, и потому вместо Ркj (t, х) при х ;> t можно

писать P ri (t)

при t > 0. Уравнение

Колмогорова — Чепмена — Смолуховского

(2.13) для однородного процесса принимает вид

 

Ло(") =

! / к i { t ) P i M - t )

( x > t )

(3.2)

 

(k, 7 =

0, 1, . .. , N).

 

2

17