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