Файл: Ганин, М. П. Прикладные методы теории марковских случайных процессов учебное пособие.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 зависит как

от разности 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