Файл: Казаков В.А. Введение в теорию марковских процессов и некоторые радиотехнические задачи.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 + Х2І (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). |
Во |
втором |
||||||
случае |
для |
каждого |
состояния задается |
н е с к о л ь к о |
|||||||
(в общем случае N—1) |
пуассоновских потоков, |
перево |
|||||||||
дящих систему из 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 (fe„ t) [8ц + |
Хц {t) At] ,= |
|
= . р « ( ' о . 0 + Д ' І ! |
h){t)Pü{ta,t). |
||
|
г=і |
|
|
Отсюда обычным образом получим |
систему уравнений |
||
d-M^n = |
j^h{t)Pil(t0lt), |
|
(2.37) |
решение которой показывает изменение вероятностей пе рехода во времени. Начальными условиями для системы (2.37) являются соотношения (2.34).