Файл: Ганин, М. П. Прикладные методы теории марковских случайных процессов учебное пособие.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 17.10.2024
Просмотров: 232
Скачиваний: 0
ГЛАВА 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 зависит как |
||
от разности t2— t\. Из (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 х(ш) = |
2к |
е - 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