Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Ганин, М. П. Прикладные методы теории марковских случайных процессов учебное пособие

.pdf
Скачиваний:
36
Добавлен:
22.10.2023
Размер:
12 Мб
Скачать

Имеем

 

 

 

/>к+1,1 — — [(^i +

>-2 +

ЫА’М — ^lPk-l.l — ^гРк.о]

(26-26)

Г1

 

 

 

(k =

i, 2, . . .) .

 

Определив р\,\ по формуле

(26.22) и вероятность ро,\ с

помощью

равенства

 

 

 

Рол— —

[O-i~Ь^2)Ро.о— PiPi,o],

(26.27)

tx2

 

 

рекуррентное соотношение (26.26) можно использовать для вычис­

ления вероятностей pk+i,i при

1.

Аналогично находятся вероят­

ности р k.j при любых значениях k и

 

Согласно (26.20) и (26.6)

имеем

 

 

О (и,

1 ) = 2

Л « Ь =

1 — «1

(26.28)

 

к=0

 

1 — Я ,И

 

 

 

Так как при а1 \а\ <

1 справедливо

разложение

(1 — я^г) - 1 = 2 (aiK)k , k=0

то для вероятности рк нахождения в системе k требований из пер­ вого потока справедливо равенство

 

Л = ( 1 - * . Х

(k— O, 1, ...).

(26.29)

Математическое ожидание числа Xi находящихся в системе тре­

бований из первого потока

 

 

 

 

сЮ (и, 1)

 

1

(26.30)

 

du

u=i

 

 

1 —' -ч

 

Математическое ожидание числа Y\ обслуживаемых требований из

первого потока

 

 

 

 

I

ео

 

 

(26.31)

 

У1 = 1 - 2 p i + s =

1 ~

Ро:

 

s=0

 

 

 

Математическое ожидание числа требований из первого потока, ожидающих начала обслуживания,

 

 

«2

 

Z, = X1 — у,

1

л

(26.32)

— GCj

 

 

Математическое ожидание времени пребывания требования в си­ стеме

*С1 =

1

(26.33)

th (1 -

 

ai)

200

М а т е м а т и ч е с к о е о ж и д а н и е вр е м е н и о ж и д а н и я н а ч а л а о б с л у ж и в а ­

н и я л ю б о го т р е б о в а н и я и з

п е р в о го п о ток а

 

 

('ОЖ1 '

 

 

х х

1

__

а,

(26.34)

 

 

xi

I1!

~ Pi

(1 — ai)

 

 

 

 

Согласно (26.30) и (26.7) имеем

 

 

о ( 1, ® ) = 2

/?*®i =

T

1 — gt — a2

(26-35)

*iu2(v) — a2v

 

 

 

j=o

 

 

 

 

 

Вероятность p* того,

что

в системе

имеется / требований из вто­

рого потока, находится с помощью равенства

 

 

Р,

1

d>G(1, v)

 

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

(26.36)

]'■

dvj

v=0

 

 

 

 

Математическое ожидание числа Х2 находящихся в системе требо­ ваний из второго потока можно определить с помощью равенства

dG ( 1, v)

 

1

CCr) Ct2

 

du2 (v)

 

X, =

dv

1

Ct i

 

dv

V=1

Так как

 

 

 

 

 

 

 

du2(v)

2Xi

X2 "T

X2 (X1+

Iх!)*

?M ’

dv

v-l

Pi — xi

p, _

 

TO

 

 

 

 

 

 

 

 

1

— ar

 

P2al

 

(26.37)

 

 

(1

"

ai)

 

 

 

Математическое ожидание общего числа требований в системе

 

 

 

х Х\ -(- х2.

 

 

(26.38)

Математическое ожидание времени нахождения требования из вто­ рого потока в системе

I

l^ i

(26.39)

iс2

Pi (1 — ai)

1а2 (1 — *1 — “ г)

 

Так как математическое ожидание времени обслуживания требова-

1

ния из второго потока равно---- , то математическое ожидание вре-

мени ожидания начала обслуживания

(последнего) можно найти

с помощью равенства

 

tож2 tc2 ~~

(26.40)

201

М а т е м а т и ч е с к о е о ж и д а н и е ч и с л а

н а х о д я щ и х с я в о чер ед и

тр е б о в а ­

н и й и з п ер во го п о то к а

 

 

22=^2^0ж2 1

(26.41)

Математическое ожидание общего числа требований в очереди

5 =

г, + *2.

(26.42)

Математическое ожидание числа обслуживаемых требований

~y =

x — z.

(26.43)

Математическое ожидание числа обслуживаемых требований из второго потока

У2 = У-У1 = У — аи

(26.44)

Чтобы найти математические ожидания tc и t0K времени пре­ бывания требования в системе и времени ожидания начала обслу­ живания (независимо от номера потока), воспользуемся формулой

М(Т) — Р(Н{)М (71///,) Р(Н2)М(Т/Н2),

(26.45)

где гинйтеза # j означает,

что

требование

принадлежит

/-му по-

 

 

j

 

то

 

току. Так как Р(НЛ — - — —— ( / = 1 , 2),

 

 

Ki

т" h

 

 

 

 

).

| >

^с1 4 ^2^с2)’

(26.46)

А>ж

^

у

(^1^ож1 4" ^ожг) •

(26.47)

Пример 26.1. Одна линия связи используется для передачи

срочных и обычных

сообщений. Срочные

сообщения не

теряются

и передаются так, будто линия связи существует только для их передачи. Обычные сообщения передаются при освобождении ли­ нии, причем если передача обычного сообщения была прервана, то это сообщение передается вновь сразу после освобождения линии связи. Потоки срочных и обычных сообщений независимые с ин­ тенсивностями А.1 = 0,04 l/мин и %2 = 0,09 1/мин соответственно. Среднее время передачи срочного сообщения 5 мин, а обычного —

3 мин 20 сек.

Определить

показатели эффективности установившегося ре­

жима функционирования линии связи.

0,09 1/мин;

щ =

Р е ш е н и е .

По условию А,] = 0,04 1/мин-, %2=

= - i - = 0 , 2

1/мин; р2=::^ - = 0,3 1/лш«. Тогда cti = — —0,2;

а2 =

t\h

 

Сч

Iх!

 

202

Вероятность

того,

что линия связи

свободна,

ро.о =

1 — ai —

— 02 =

0,5. Вероятность того, что не передается срочное сообщение,

Ро — 1 —ai =

0,8. Вероятность наличия очереди из

k срочных сооб­

щений

/?к+1 = (1 — <*i) а^+! = 0,8 (0,2)k+1

(к = 0, 1 , ...).

 

 

Математическое ожидание числа срочных сообщений: в системе

а,

0,25;

в процессе

передачи

0,2; в

очереди

Х\= --------- =

у\ = а, =

_

0,05.

Математическое ожидание времени от момента

z,=X i — г/1 =

поступления срочного сообщения до окончания его

передачи

^ci =

—-^ -= 6 ,2 5

мин,

а математическое ожидание времени ожидания

начала передачи toxl =

— =

1,25 мин.

 

 

 

 

 

 

 

 

 

 

ку

 

 

 

 

 

 

Математическое ожидание числа обычных сообщений в системе

х 2

 

 

 

1 +

р,а,

6,6

=

0,825 .

 

 

 

 

В-1 (1 — ai)

 

 

 

 

 

 

 

8

 

 

 

 

Тогда х = Ху -j-x2— 1,075.

времени

от момента

поступления

Математическое

ожидание

 

 

 

 

 

 

 

-j

 

х 2

0,82о

обычного сообщения до окончания его передачи гс2= - т — — п лп ~

= 9 мин 10 сек,

а

математическое ожидание времени до начала

передачи (последней) £ож2 = tc2— 3 мин 20 сек = 5 мин 50 сек.

Ма­

тематическое ожидание числа обычных сообщений

в

очереди

z2 =

=-■ h tn

0,09

35

0,525.

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

Общее число сообщений, ожидающих передачи, z — zt -j- z2 —

=0,575. Математическое ожидание числа передаваемых сообщений

у= х z = 0,5, в том числе обычных сообщений у2 = у У\~ 0,3.

Математическое ожидание времени от момента поступления со­ общения до его полной передачи

^ = - ^ - ( 2 5 + 82,5] : 8,27 мин,

а до момента начала передачи

7

52,5) = 4,42

мин

^ож

^ W “1

ГЛАВА 3

ПРИКЛАДНЫЕ МЕТОДЫ ТЕОРИИ НЕПРЕРЫВНЫХ МАРКОВСКИХ ПРОЦЕССОВ

Непрерывные марковские случайные процессы, т. е. марковские процессы с непрерывными ординатами и непрерывным временем, широко используются в приложениях к различным областям науки и техники. Эти процессы тесно связаны со стохастическими диф­ ференциальными уравнениями, с помощью которых описывается функционирование различных физических систем. Основная слож­ ность применения теории непрерывных марковских процессов со­ стоит в возможности получения решения соответствующего диф­ ференциального уравнения в частных производных при заданных начальном и граничных условиях.

В данной главе излагаются прикладные методы теории непре­ рывных марковских случайных процессов. Дается вывод уравнений Колмогорова, которым удовлетворяет плотность распределения про­ цесса, и устанавливается связь этих уравнений со стохастическими дифференциальными уравнениями. Рассматриваются начальные и граничные условия, которым удовлетворяет искомая плотность рас­ пределения. Природятся некоторые точные и приближенные методы решения уравнений Колмогорова, исследуются выбросы одномер­ ных марковских процессов. Теория одномерных процессов обоб­ щается для многомерных марковских случайных процессов. Рас­ сматривается возможность замены реального случайного процесса непрерывным марковским процессом.

§27. ТИПЫ СЛУЧАЙНЫХ ФУНКЦИЙ

СНЕПРЕРЫВНЫМИ ОРДИНАТАМИ

При любом фиксированном значении неслучайного аргумента t

случайная функция X(t) с непрерывными

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

непрерывной случайной величиной. Если t\,

t2, .. ., tt — произволь­

ные значения аргумента t, который будем называть временем, при­

чем ^1 <

t2< . . . <

tb ТО совокупность X(t\), X(t2) ,

. . . , X (t[)

образует

систему

I непрерывных случайных величин.

Полной ха-

204

рактеристикой этой системы является /-мерная функция распреде­ ления

 

и

■ • •- Щ tx, t,,

. . . , i,) =

 

=

P\X{ty) < x x,

X (t2) < x 2, .

.

. , * ( * , ) < * , ] ,

(27.1)

которая зависит от / аргументов хи х2,

 

... , хх и от / параметров

t\, t2,

th Функция

распределения

(27.1) непрерывная и диф­

ференцируемая по всем аргументам. Вместо этой функции часто удобнее использовать /-мерную плотность распределения, которая определяется формулой

 

fl(.Xj, Х2, . • ■,

Xf,

tx,

t2t • .

•,

/;)

 

__

(^i, X2, . .

. ,

Xx,

tx, t2, .

.

. ,

tj)

(07 0\

 

dxxdx2 .

.

.. dxt

 

 

'

' '

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

только математическое ожидание x(t), дисперсия D[X(/)] и корре­ ляционная функция Для их определения достаточно знать законы распределения первого и второго порядка, так как

 

 

x ( t) =

j

xfj(x;

t)dx;

 

 

 

(27.3)

 

 

 

—• oo

 

 

 

 

 

 

 

 

 

 

 

 

00

 

 

 

 

 

 

D \X(/)] =

Kx(t, t ) =

j

\x — x(t)}2ft(x; t)dx\

(27.4)

Kx (?, t”) =

M {\X (f) -

 

x {t’)\[X (/") -

x (/")]! =

 

 

00

 

 

 

 

 

 

 

 

 

 

 

= JJ \х?-х(Г)\ \x” - * ( / " ) ] f, (x't x"\ t\

t") dx' dx",

(27.5)

причем

 

 

J h{x, x'-, (, t')dx'.

 

 

 

 

f\(x;

0 =

 

 

(27.6)

Случайная функция

X(t)

называется стационарной,

если

для

нее законы распределения любого порядка /

(/ =

1 , 2, ...)

не зави­

сят от начала отсчета параметра /, т. е. если

при

любом /о

tx

справедливы равенства

 

 

 

 

 

 

 

 

 

 

fi(Xi, х 2, .

. . ,

Хх,

/],

t2, . .

. ,

/;) =

 

 

 

— fi(xi> x 2> • • ’ > x i\

tx

 

t0,

t2

^oi •

• - i

h

to)

(27.7)

 

 

(* = 1 , 2, , . . ) ,

 

 

 

 

 

205

В ч а с т н о с т и , есл и п р и н я т ь io = t\, то

f\(xi;

= (xi;

0);

(27.8)

f2(-^ь x2\ tu t2)=

f2(x1, *ai

0, t2 — tx),

(27.9)

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

пределения системы [X(/t), X(t2)]

от параметров tx и t2 зависит как

от разности t2t\. Из (27.3) — (27.5)

следует, что в

этом случае

математическое ожидание x(t) и

дисперсия Z)[Z(/)]

постоянные,

а корреляционная функция Кх ((',

I")

зависит от ¥ и t"

как от раз­

ности ¥' ¥, т. е.

 

 

 

x(/) = const; D\X(/)]■=

(0) = const;

(27.10)

KAt\ n = KAt” - t ' ) .

(27.11)

Из условий стационарности (27.7)

случайной функции X (/) кроме

(27.10) и (27.11) можно получить также другие свойства различ­ ных вероятностных характеристик. Ясно, что при выполнении (27.10) и (27.11) условия стационарности (27.7) случайной функ­ ции X(t) могут не выполняться, а потому при этом X(t) может быть и нестационарной случайной функцией.

В приложениях наиболее часто используется так называемая корреляционная теория случайный функций, когда необходимо

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

тематическое ожидание x(t) постоянно, а корреляционная функ­ ция Кт. {¥, ¥') от ¥ и t" зависит как от разности, т. е. выполняются условия (27.10) и (27.11). Любая стационарная функция является стационарной и в широком смысле. Если же случайная функция X(t) стационарна в широком смысле, то она не обязательно яв­ ляется стационарной.

Случайная функция X{t) называется нормальной, если при лю­ бом I система случайных величин X(ti), X(t2), .. ., X{tt) нормаль­ ная. Известно, что /-мерная плотность распределения нормальной

системы

(Хи Х ^ . .. , Xt) полностью определяется матричной стро­

кой (Xi,

х2, . . . , Х[), составленной из математических ожиданий ис­

ходных случайных величин, и элементами

k}S(/,

s = 1 , 2, ..., /)

корреляционной матрицы

/С =

|^JSII•

Так

как

Х^ —- X (/j) ( / =

= 1) 2, ..., /), то х j = х (/j)

(/ =

1 , 2, ... , /);

 

 

kis= М [(A'j — XjK^s— *з)] = Kx(tp

tt) (/,

s =

1, 2, . . . , /).

206

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

ческое ожидание x(t) этой случайной функции в произвольный мо­ мент времени / и через корреляционную функцию Кх(/', t") для любых f и t".

Если нормальная случайная функция X(t) стационарная в ши­ роком смысле, то х(^) = const, a Kx(t', t") = / ( х(т), где %= ?' t'. В этом случае /-мерная плотность распределения полностью опре­

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

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

1

f1 (*; t) —

У & Щ х Щ

| Х - X ( t ) ] »

2D[X(t)]

(27.12)

Если нормальная случайная функция X(t) стационарная, то (27.12) принимает вид

 

 

 

 

(х - х)«

 

 

fl (х) =

 

1

е

( 2 7 . 1 3 )

 

ах1/

 

 

 

 

2тс

 

ГД6

sl = D [ X ( t ) } = K A 0).

 

Двумерная плотность распределения стационарной нормальной слу­

чайной функции X(t)

определяется формулой

 

 

х 2; ')

 

 

 

ехр

-

1

-[(Xj — X)2

 

 

 

2а2[1

(х)]

2ттох /

1

- Г

2Х(Т)

 

 

— 2rx(x)(xi

-

х)(х2 — х ) + ( х 2 — х )2]| ,

(27.14)

где

 

 

 

 

 

 

 

 

 

X

 

* =

 

.

(27.15)

 

 

 

 

 

 

 

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

 

i

 

 

 

fi (*i,

. . . ) Х[, ti, /2, . . . , /;) — Г

"

/j)

(27.16)

 

]«>i

 

 

 

 

(/= = 1, 2, . . . ) .

 

 

 

207

При этом ординаты X (tl) l X (t2),

... , X (/,) случайной

функции

X{t)

взаимно независимы, а потому и

не коррелированы, т. е.

Knit',

t") = 0 при f ф t". Математическое ожидание x(t)

и диспер­

сия D[X(/)] случайной функции

X (/)

находятся

по

формулам

(27.3)

и (27.4). Из (27.16) следует, что

процесс с

независимыми

ординатами полностью определяется одномерной плотностью рас­ пределения / i(x; t).

Стационарный чисто случайный процесс X(t) с бесконечной дисперсией П[Х(()] называется белым шумом. Корреляционная

функция белого шума определяется формулой

 

Кп(х)=с8(т),

(27.17)

где с — положительная постоянная, называемая интенсивностью бе­ лого шума, а 6 (т) — дельта-функция Дирака, т. е.

(

0

при

т ф 0 ;

{

со

при

(27-18)

т = 0.

Если математическое ожидание случайной функции X{t) по­ стоянно, а корреляционная функция выражается формулой (27.17), то X[t) называется стационарным в широком смысле белым шу­ мом. Ординаты этого случайного процесса некоррелированные, но в отличие от чисто случайного процесса могут быть зависимыми. Спектральная плотность белого шума (стационарного п стационар­ ного в широком смысле)

5 х(ш) =

е - 1°«КЛт)<Ь = 4 -

(27.19)

 

 

 

Название «белый шум» для стационарной случайной функции с по­ стоянной спектральной плотностью объясняется аналогией с белым светом, который представляет собой сумму (суперпозицию) всех спектральных составляющих с одинаковыми интенсивностями. Бе­ лый шум можно представить в виде суммы гармонических колеба­ ний всех частот с одной и той же дисперсией амплитуды.

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

ментов ti < t2< ... < t,

случайная величина X(t[)

и разности

X (£j) — X (^j—0 0 — 2, 3,

. . . , / ) взаимно независимы.

Чтобы найти

/-мерную характеристическую функцию для X(t), воспользуемся равенством

2

BjA'(^j) = ®1A ’U 1) + 2 ® . I A , ( M - A '( ^ _ 1)],

(27-20)

j—1

s=2

 

где

i

 

 

 

 

-ys= 2 « j .

(27.21)

 

j-s

 

208

Используя определение характеристической функции для системы [^ (/1), X{t2), ... , X{tt)\ случайных величин, получаем /-мерную ха­ рактеристическую функцию случайной функции X{t) в виде

Ех {щ, и2, .

.

. , щ\ tu /2, . .

. , tt) = М

j=i

 

 

 

 

: М [eiV‘X(U] f ] М {^iVstX(ts)-X (ts_ 1)l j

(27.22)

 

 

s=2

 

 

 

Из этого выражения следует, что

 

 

 

Ei {lli,

м25 • • •i Щ> / ц ^2i • • •»

/;) =

 

 

 

I

 

 

 

=

Ei(Vu ^ ) П ^ г ( —

*V. *.-i,

tt).

(27.23)

 

 

s=2

 

 

 

Так как одномерная характеристическая функция £i(o; t) по­ лучается из двумерной характеристической функции E2(v, v t , t') с помощью равенства

Ex{v, t) = E2(v, 0; t, Г),

(27.24)

то /-мерная характеристическая функция случайного процесса с не­ зависимыми приращениями при любом / определяется через дву­ мерную характеристическую функцию. В свою очередь, характери­ стическая функция E2(v, v'\ t, t') и двумерная плотность распре­ деления f2(x, х'; /, Е) связаны равенствами:

 

со

 

 

 

 

М * . ■*'; t,

f) = ^ y r J j* « - lvxiv'x'£ '2 ('y. v'-,

t,

t')dvdv'\

(27.25)

 

— oo

 

 

 

 

J

v'\ t, t')— J j

 

 

 

(27.26)

E2(v,

f2{x, x';

t,

t')dxdx' .

 

«—oo

 

 

 

 

Следовательно, случайный процесс с независимыми приращениями полностью определяется двумерной характеристической функцией или двумерной плотностью распределения.

Для любой случайной функции /-мерную плотность распределе­ ния можно представить в виде

f i { x 1 » -^2» ■ •

• »

Xi,

ti, t2, . •

. , ^ l ) f i {Xi\

t\)f{x2,

t2jXi,

/j) X

Х/С^з»

X2,

/j,

/ 2). . . / (Xi,

tijxj) АГ2) •. ■;

•£;—liA;

/г> ■••;

i)-

 

 

 

 

 

 

 

(27.27)

14

209

Соседние файлы в папке книги из ГПНТБ