Файл: Ганин, М. П. Прикладные методы теории марковских случайных процессов учебное пособие.pdf

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

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

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

Добавлен: 17.10.2024

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

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

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

jj

т-2

(t, x\ t +

At)

(29.33)

A t-*0

A t

 

 

 

lim m s ( t,

x ; t

+ A t )

(s ^ 3 , 4,...).

(29-34)

A t - 0

A t

 

 

 

Аналогично процессу броуновского движения со сносом может быть построен любой непрерывный марковский процесс, являю­ щийся суммой большого числа малых приращений. Можно также любой непрерывный марковский процесс рассматривать как реше­ ние так называемого стохастического дифференциального уравне­ ния, которое содержит процесс броуновского движения в качестве «вынуждающей силы». Вследствие этого непрерывный марковский процесс называется также диффузионным. Название «диффузион­ ный» происходит от того, что с помощью непрерывного марковского процесса достаточно точно описывается диффузионное движение частицы, т. е. медленное проникновение частиц одного вещества, например жидкости, в другое при их непосредственном соприкос­ новении. Перемещение такой частицы из точки х в момент t за ма­

лое время At можно представить в виде

суммы a(t, x)A£-f-6p-f-

+ 0 (Д/), где a(t, х) — некоторая средняя

скорость, а 6ц — случай­

ное перемещение, вызванное столкновениями частицы с находящи­ мися в тепловом движении молекулами жидкости. При этом

М,[(6ц)2] = ^(/, x)At + 0(At),

где b (t, х) — некоторая функция, пропорциональная средней кине­ тической энергии молекул жидкости в окрестности точки х в мо­ мент t. Функции a(t, х) и b(t, х), характеризующие непрерывней марковский процесс, называется коэффициентами сноса и диффу­ зии соответственно. Данные функции определяются формулами:

a(t, х) — lim--7 W?! (f, х;

t +

At)-,

(29.35)

At-*-0 A t

 

 

 

b(t, х) ~ lim —itn-zit, х\

£ +

Д£),

(29.36)

At-» 0

 

 

 

где rri\(t, х\ т) и ш2(i, х\ т) — условные моменты первого н второго порядка, для которых справедливы выражения (29.20) и (29.21).

Условный начальный момент первого порядка тх(/, х; т) яв­ ляется средним перемещением частицы за время т — t из фиксиро­ ванного положения х в момент t в результате случайных воздей­ ствий. Из (29.35) и (29.20) следует, что функция a(t, х) является средней скоростью изменения процесса в точке х в момент t или, что то же самое, средней скоростью изменения состояния х в мо­ мент t. Условный момент второго порядка tn2(t, х; т) пропорциона­ лен математическому ожиданию квадрата изменения случайного процесса и является мерой разброса значений Х(^) относительно

222


исходного значения х. Согласно (29.36) и (29.20) функция b(t, х) характеризует скорость изменения условной дисперсии марковского процесса в точке х в момент t. Если к описанию марковского про­

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

случайных толчках

и л и случайной силе, действующей на систему,

то функция b(t, х)

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

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

случае

a(t,

х )= а , b(t,

x) =

b. Кроме

того, для эдого процесса вы­

полняются условия

(29.34),

согласно

которым начальные моменты

ms (t, х;

t +

At) (s =

3,

4, ...) при At

-> 0 стремятся к нулю быст­

рее, чем М. С помощью условной плотности распределения равен­ ства (29.34) записываются в виде

со

 

 

lim 7 : Г(v — x)sf(t, х;

t-\~ М, у) dy 0

(29.37)

— со

 

 

(s = 3, 4,

...).

 

Чтобы определить любой другой непрерывный марковский про­ цесс X(t) , необходимо знать его основные характеристики, кото­ рыми являются непрерывные функции a{t, х) и b(t, х), т. е. коэф­ фициенты сноса и диффузии. Данные функции могут быть вычис­ лены по формулам (29.35) и (29.36), т. е. они выражаются через условную плотность распределения f(t, х; т, у). Достаточным усло­ вием непрерывности любого марковского процесса X(t), кроме су­ ществования непрерывных функций a(t, х) и b(t, у ), является вы­ полнение равенства

 

со

 

 

 

 

 

lim -^7 Г|у — x\2+hf(t, х;

t - { - l t , y ) d y

= 0

(29.38)

■Ц-0

 

 

 

 

 

при некотором положительном б, 0 < б ^ 1 .

 

 

Обозначим

через

РА вероятность

того, что |X(t -f- At) X{t) |

при условии

X(t) =

x

будет больше

заданного

положительного

числа е. Для этой вероятности имеем

 

 

 

 

Р * =

 

Г f ( t , x ; t + i t , y ) d y * C

 

 

 

 

I

У - Х | > Е

 

 

 

223


|у — JC|2+8 / (t, x ; t + At, у ) d y K

Iу-* l>s

< ~^+r 11 У — ■* l2+! / (*> ■*; H - A*, У) rfy.

(29.39)

oo

Сучетом (29.38) из (29.39) следует, что

lim ^ = 0.

(29.40)

4t-0 At

 

Следовательно, при выполнении условия (29.38) вероятность боль­ ших изменений |X(f-f-At) — X(t) |, когда X(t) = x, при At —» О стремится к нулю быстрее, чем At. Нужно отметить, что предельное

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

(29.38) не только при S > 0,

а уже при 8 > —2. Условие (29.38)

более жесткое, чем (29.40),

поэтому его называют условием усиленной непрерывности. Согласно

этому условию вероятность больших изменений

|^(^ + А/) — Я (/) |

при X{t) = x и At 0 настолько

быстро стремится к нулю, что

М [|X(t -f At) — X(t) \2+i/X{t) = x]

при At -*■0

убывает быстрее,

чем At. Это требование применительно к системе, подвергающейся воздействию случайных толчков, позволяет рассматривать случай ■ ный процесс X(t) как непрерывную функцию, т. е. как среднее за время, значительно превосходящее промежуток между случайными толчками.

§ 30. ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ ДЛЯ УСЛОВНОЙ И ОДНОМЕРНОЙ ПЛОТНОСТЕЙ РАСПРЕДЕЛЕНИЯ

Коэффициент сноса a(t, х) и коэффициент диффузии b(t, х) для непрерывного марковского случайного процесса X(t) опреде­ ляются формулами:

a(t,

х) — lim ±7M {[ X {t + bt) — X(t)\IX(t) = x\\

(ЗОЛ)

b{t,

x) = lim ^ M { [ X ( t + A t ) - X(t)]2/X(t)=x}.

 

(30.2)

 

At-*0

 

 

Если известна условная плотность распределения f(t, х\

т,

у), то

для вычисления коэффициентов сноса п диффузии можно

исполь­

зовать равенства:

 

 

оо

 

 

 

 

a{t, х ) =

lim -^

j (у — x)f{t,

х;

t +

It, y)dy;

(30.3)

 

 

— oo

 

 

 

 

 

 

00

 

 

 

 

b (t, x) =

U m ^

Г(у - x f f { t ,

x\

t +

At, y) dy.

(30.4)

224


При практических приложениях коэффициенты сноса и диффу­ зии обычно известны или находятся с помощью формул (30.1) и (30.2). Условная плотность распределения, условная функция рас­ пределения, одномерные законы распределения и другие функции, характеризующие непрерывный марковский случайный процесс, находятся как решение соответствующих дифференциальных урав­ нений в частных производных, коэффициенты которых выражаются через a(t, х) и b(t, х ). Чтобы получить одно из таких дифферен­ циальных уравнений относительно условной функции распределе­

ния F(t, х; т,

у),

воспользуемся уравнением (29.19), которое после

замены t на i — At и t' на i

записывается в виде

 

F(t — It, х;

-г, у) = J F{t,

z; -, y)f(t — At, х; t, z) dz.

(30.5)

 

 

— oo

 

 

Предположим,

что условная функция распределения F (t,

z\ т, у)

имеет непрерывные частные производные первого и второго по­

рядка по z при любых t,

z, у и т > t.

Тогда эту функцию

можно

разложить в ряд в окрестности точки

z = x, записав первые три

члена разложения в явном виде, а остаточный

член — в

форме

Пеаио. При этом получаем

 

 

 

 

 

 

z;

х, у) = F(t, х; т, у ) +

-

х) — ^

 

+

+

у ( г — X? -

- % Х1 Т’

У)

+ о [(z -

х )2].

(30.6)

Последнее слагаемое из (30.6) удовлетворяет условию

 

 

Пт

о [(* -

ху\

= 0,

 

(30.7)

 

Z->X

(z -

х )2

 

 

 

так что можно записать

 

 

 

 

 

 

 

0 \{z — х )2] = |z

— х |2+56 (х, z),

 

(30.8)

где б — некоторое положительное число, а ф (х, z) — ограниченная функция, т. е. |ф (х, z) |< С < оо. Подставляя (30.6) в (30.5), с учетом равенства

оо

 

 

j

х; t, z ) d z = 1

и(30.8) получаем

^[/=■(*, х; х, у) - F{t It, х; х, у)] +

15

225


 

+

dF(t, x;

т,

у)

1

 

(z

х) f(t At, х ; t, z) dz-f-

 

 

дх

 

 

At

 

 

 

 

 

 

 

 

 

 

 

 

 

_L

d2F{t, x\ -t,

у)

1

 

oo

 

 

 

 

 

 

(z — x )2 /(£ — Дt, x; t, z) dz -)-

 

 

2

дх2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

oo

 

 

 

 

 

 

 

 

 

-f-

Гty (x, z) |г — x |2+Г|f{t At,

x; t, z) dz =

0.

(30.9)

Первое слагаемое пз левой

частиэтого равенства при At -> 0

стре-

мптся

к

dF(t, х;

т,

у)

.

m

v

, .

 

dt

 

1

ак

как функция г|:(х,

г) ограничена,

то в

 

с

(29.38)

последнее

слагаемое после перехода

соответствии

кпределу обращается в нуль.

Сучетом этих замечаний и выражений (30.3), (30.4) из (30.8)

после перехода к пределу при -At —> 0 приходим к следующему диф­ ференциальному уравнению в частных производных относительно условной функции распределения:

dF(t, х; т, у)

 

a{t,

х)

dF(t, х;

т, у)

 

dt

 

 

 

 

 

дх

 

 

+

iyb(t,

х)

d2F(t, х; т,

у) =

0.

(30.10)

 

 

 

 

 

дх2

 

 

 

Если уравнение (30.10)

продифференцировать один раз по у и вос­

пользоваться равенством

(29.9), то получим

 

 

 

д/(*. х;т, у)

,f .

df{t,

х ; т,_у)

 

dt

 

\

 

 

 

 

 

,

1 , ,,

 

d2f(t, х;

т, у)

 

(30.11)

+

2 b{t’ x ) --------- Ш

---------- ° -

 

Данное равенство,

называемое

первым

(обратным)

уравнением

Колмогорова, характеризует изменение условной плотности распре­

деления f(t, х; т, у)

в зависимости от начального состояния, т. е.

от t и х. Конечное

состояние системы, т. е. параметр т (т > t) и

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

Наиболее часто известно не конечное, а начальное состояние си­ стемы. Поэтому необходимо иметь уравнение, в котором условная плотность распределения f(t, х; т, у) рассматривается как функция конечного состояния, т. е. в зайгсимости от т и у, а начальное со­

226