Файл: Казаков В.А. Введение в теорию марковских процессов и некоторые радиотехнические задачи.pdf

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

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

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

Добавлен: 28.06.2024

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

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

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

з а т е л ь н о м у

( э к с п о н е н ц и а л ь н о м у )

закону

с параметром

X. Показательное распределение

интерва­

лов времени между событиями обладает одним замеча­ тельным свойством: если после очередного события про­

шло уже некоторое

время

Х\,

то

закон

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

оставшейся

части промежутка

0

= 7"ті

остается неиз­

менным при любом

Т і > 0 ,

т. е.

 

 

 

 

 

и,(Ѳ)=; яе-м .

 

(2.17)

Докажем

это свойство

(см., например, [4, 5]). Совме­

стим начало отсчета с моментом появления некоторого

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

(0, %\) и

(fi,

Ті + Ѳ).

Для

этих интервалов справедливо

соотно-.

шение (2.10):

 

 

 

 

 

 

 

 

 

 

Р,1 + е(0) =

*Ч(°)М°)-

(2-1 8 )

 

Кроме того, согласно (2.9)

имеем

 

Из

(2.18),

(2.19)

находим

 

 

 

 

 

 

 

Р0(О) = ^

^

=

е - х \

(2.20)

что

полностью

совпадает с формулой (2.14), на

основа­

нии которой был получен закон распределения

(2.16).

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

выражение

(2.17)

верно.

 

Совпадение

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

w{x)

и ИУ{Ѳ) подтвержда­

ет свойство марковости пуассоновского-потока. Дейст­ вительно, закон распределения ш(Ѳ) оставшейся части промежутка («будущее») не зависит от того, сколь долго

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

(0, ті) с момента

осу­

ществления последнего события

(«прошлое»), но

при

этом плотность ш(Ѳ) зависит от Ті («настоящее»), так как величина Ѳ отсчитывается с м о м е н т а ті.

Из всех видов распределений интервалов между со­ бытиями только экспоненциальное распределение обла­ дает указанным свойством, следовательно, пуассоновский поток (в общем случае нестационарный) является единственным, который может определять моменты вре­ мени переходов в математической модели, описанной в '§ 2.1. Поскольку каждый пуассоновский поток пол­ ностью характеризуется интенсивностью \(t), то для полного задания разрывного марковского процесса с ди-

43


скретными состояниями необходимо располагать N зна­

чениями

величины Х(ХІ, t)

( i = l , 2, ..., N).

При контину­

альном

множестве

возможных состояний

марковский

разрывной

процесс

должен

определяться

функцией

Х(х, t).

 

 

 

 

 

 

 

Марковские процессы,

у

которых

интенсивность

является

функцией

состояний х, играют большую роль

в теории

массового

обслуживания и связанных с ней

прикладных

областях [5—8].

 

 

 

2.3. Дифференциальные уравнения для разрывных марковских процессов с дискретными состояниями

Итак, для описания системы с конечным числом со­ стояний N следует задать в общем случае N. пуассоновских потоков, которые характеризуют случайные мо-

^менты переходов системы из состояния в состояние. Однако для определения вероятностей перехода в систе­

ме знания одних потоков, очевидно, недостаточно. Необ­

ходимо еще

ввести в

рассмотрение о т н о с и т е л ь н ы е

вероятности

переходов

ац(і), которые являются услов­

ными по отношению к факту возникновения скачка. Ве­

роятность

ga(t) есть

вероятность

того,

что

система,

находившаяся до момента t в состоянии ХІ,

перейдет в со­

стояние Xj

п р и у с л о в и и ,

ч т о

в м о м е н т

t п р о ­

и з о ш е л

с к а ч о к .

Таким

образом,

если скачок не

происходит, то система

остается в прежнем состоянии х\.

Это обстоятельство делает естественным

предположение

о том, что элементы главной диагонали матрицы отно­

сительных вероятностей перехода Q(t)

следует положить

равными нулю, т. е.

 

 

 

 

qu(it)=0

для всех

г.

(2 . 21)

Поскольку

матрица Q{t)

является

стохастической, то

для каждой

ее строки справедливо условие

нормировки

Запишем теперь собственно

вероятности переходов

в системе за малый промежуток

(t, t+àt), предположив,

что поглощающие состояния отсутствуют. Выберем не­ которое состояние ХІ, которым управляет пуассоновский 44


поток с интенсивностью li(t).

В силу малости

из (2.7)

имеем

 

 

 

 

 

 

 

 

 

 

 

Ai(t,bt)=

 

t+àt

%i{t')dt'

=Xt{t)àt.

 

 

 

j

 

 

 

 

 

t

 

 

 

 

интервале (t,

Тогда вероятность

отсутствия скачка

на

£+ At)"согласно

(2.6) равна

 

R

 

 

 

 

Pt.

At (°) = е _

Ч ( ° " * 1 -

k (0

 

(2-2 3 )

С учетом

(2.2),

(2.23)

вероятность

того, что на

рассма­

триваемом

интервале

произойдет, скачок,

определяется

соотношением

 

 

 

 

 

 

 

 

 

 

 

 

PtiU(l)*b®M.

 

 

 

(2.24)

Выражения (2.21) — (2.23)

позволяют

записать

вероят­

ности перехода.из состояния ХІ В виде

 

 

 

 

 

pii{t,t+M)

 

= \—%i{t)U,

 

 

(2.25)

 

 

Pij{t,

t+ät)=U(t)btqijit).

 

 

 

(2.26)

Безусловная

вероятность pj(t+At)

 

того, что система

в .момент

t+At

окажется в состоянии Xj равна

 

 

^ (*+• дд

 

 

[ 1 - я , (О Д*] +

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

Pi{t)li{t)Atqu{t).

 

 

 

(2.27)

t=i

Перенесем pj(t) в левую часть, разделим обе части урав­ нения (2,27) на At и перейдем к пределу при Arf—>0. В результате получим

 

 

N

^

= -li(f)Pi(t)+^Pi(t)lt(f)gm

0 = 1 . 2 , ...,#).

(2.28)

Система уравнений (2.28) характеризует изменение без­ условных вероятностей с течением времени. Впервые она была получена А. Н. Колмогоровым в работе [9].

Положим теперь

- Я , ( / ) = а и ( 0 ; M O t t H O ^ ^ j P ) . іфІ-

(2.29)

45


Коэффициенты üij(t) *) позволяют заіписать систему уравнений (2.28) в компактном виде

 

 

dpi (0_<

N

 

 

 

 

 

 

(2.30)

 

 

dt

=1

 

 

 

 

і

 

 

или

в матричной

форме

 

 

 

 

 

£ £ L =

P(*)-A(Q,

(2.31)

где

Р(і)матрица-строка безусловных

вероятностей;

А(і)

—матрица,

составленная

из элементов ua(t).

Рассмотрим вывод системы дифференциальных урав­

нений (2.30), основывающийся

на ином,

нежели (2.25),

(2.26), представлении вероятностей перехода. Предполо­

жим, что при малом

любую вероятность перехода

можно записать в

виде

 

 

 

Pi){t,

t+At)

=

ôij+%ij(t)M,

(2.32)

где •ôfj — символ Кронекера;

Xij(t)—интенсивность

пу-

ассоновского потока, переводящего систему из состояния

ХІ в состояние Xj. Матрица

вероятностей

перехода с эле­

ментами- (2.32) имеет вид

 

 

 

 

 

 

 

 

 

P(f,f

+

Af)

=

 

 

 

"1 +

\ u (0 At

А1 2 (/) At

 

.

. .

Kw

(t) At

 

Х2 І

(/) At

1 + Х(t)At

.

. .

X2N

(t) At

(2.33),

 

 

 

 

 

 

 

 

 

X,„ (t) At

Лот(0

à

t

 

 

1 + Х Ш ( 0 At

 

L "Ni

 

 

 

 

 

 

 

Форма записи вероятностей перехода (2.32) учиты­ вает ряд обстоятельств. Во-первых, при At—И) конечное состояние Xj совпадает с начальным х,-, т. е.

 

Pü(t, '0 =6ij-,

(2.34)

во-вторых, при

\фі

 

 

 

 

 

(2.35)

Соотношение (2.35) означает, что вероятность пере­

хода из ХІ в Xj

( і ф і )

выражается

через интенсивность

*> Величины a.ij(t)At

называют инфинитезимальнымп (локаль­

ными) переходными вероятностями.

 

46


kij(t)

пуассоновского

потока,

переводящего

систему

из

г-го состояния в /-е. Сопоставление

выражений

\2.26)

и

(2.35)

позволяет уяснить разницу

в

описании

системы.

В первом случае для каждого

состояния задается о д и н

пуассоновский поток с интенсивностью Xi(t)

и набор от­

носительных

'вероятностей

перехода

qij{t).

Во

втором

случае

для

каждого

состояния задается

н е с к о л ь к о

(в общем случае N1)

пуассоновских потоков,

перево­

дящих систему из t'-ro состояния во все остальные.

 

Что касается «перехода» ХІ*-Х{,

 

то, исходя

из усло­

вия нормировки

 

 

 

 

 

 

 

 

 

 

I ! Ра (t, t + M) =

S

[8„ +

Ut (t) Ы] = 1,

 

 

 

/=і

 

 

/=і

 

 

 

 

 

 

 

откуда следует

N

 

 

S

а«(*)д* = о.

 

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

 

 

Xu (0 =

- 2 Ui (*) < 0-

(2.36)

Соотношение (2.36) противоречит определению интен­ сивности потока, приведенному в § 2.2, и поэтому функ­ цию А,ІІ(./) можно лишь условно назвать интенсивностью. Несмотря на это, выражение (2.36) хорошо согласуется с описанием протекающего в системе процесса при по­ мощи графив (см. ниже).

Подставим вероятность перехода (2.32) в уравнение

Колмогорова—Чепмена

(2.1):

 

 

Pii(ta, t + Д*) •= S

pa (fet) [8ц +

Хц {t) At] ,=

= . р « ( ' о . 0 + Д ' І !

h){t)Pü{ta,t).

 

г=і

 

 

Отсюда обычным образом получим

систему уравнений

d-M^n =

j^h{t)Pil(t0lt),

 

(2.37)

решение которой показывает изменение вероятностей пе­ рехода во времени. Начальными условиями для системы (2.37) являются соотношения (2.34).