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

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

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

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

Добавлен: 17.10.2024

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

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

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

С учетом (32.13) и (32.15) из (32.11) и (32.12) получаем следую­ щие выражения для коэффициентов сноса и диффузии стохастиче­ ского дифференциального уравнения (32.1):

a(t, *) = <р(х, 0 ; b(t, *) = [ф(*, t)f.

(32.16)

Данные выражения называются коэффициентами сноса и диффу­ зии по Ито.

Равенства (32.16) получены при замене функции ф[Х(^), t'] из

(32.6), когда t

t -f- н Д^-9-0, на ip[^(^), f\. Подобная замена

производится при

определении общего стохастического интеграла

по Ито. Указанный стохастический интеграл содержит недифференцируемую случайную функцию (процесс броуновского движе­ ния), а потому правила преобразования и вычисления этого ин­ теграла отличаются от соответствующих правил для дифференци­ руемых функций. Например, если w(t) — дифференцируемая функ­ ция, то

 

t

 

J

w(x)dw(i) = — \\w(t)Y— [®>^o)]2}-

(32.17)

to

 

 

Если же W(t)

— процесс броуновского движения, то при определе­

нии стохастического интеграла по Ито получается

 

t

'

 

j V СО dW(x) = ± { { w (О]2 - \W(*0)]2 - ( t - t0)) .

•(32.18)

to

 

 

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

жения.

Такое определение

дано Р.

JI. Стратоновичем. При этом

в (32.6)

подынтегральная

функция

ф Р ф '), t'] заменяется

не

на

ф [*(0 , А. а на разложение в ряд с учетом первых двух членов, т.

е.

принимается ty[X(t'),

П =

^ + ^х[^(0 > (ДО,

 

 

где

 

 

 

 

 

 

 

+

 

l(n d t"- + Q(t' - t) .

(32.19)

Тогда вместо (32.10)

будет

t

 

 

 

 

t+ At

 

 

 

bZ, (t) =

Ф[X (t),

 

 

 

t\ j \ (t') dt’ -f-

 

 

 

 

 

 

t

 

 

245


 

 

 

 

 

t+At

 

 

 

 

 

+ < ?;!*(*),

*]| ? [* (* ),

t\ j ( t ' - t m n d t ' ±

 

 

 

t + At

 

 

 

 

 

 

-h Ц* [ЛГ(0 ,

t]

j

и

 

Л '1 +

0(Л*).

(32.20)

 

 

 

t

 

 

 

 

 

Подставляя (32.9)

и (32.20)

в (32.7), получаем

 

 

 

 

 

 

 

 

 

t + A t

t

 

а (t, х)

= <р(х, t) +

{х,

t) <Ь{х, t) liin

1

\

\ Ь{ f -

t") dt" d f .

 

 

 

 

 

At-^O

 

 

 

 

 

 

 

 

 

t

t

 

(32.21)

 

 

 

 

 

 

 

 

 

Полагая f

tf/ т, с

учетом (32.14) находим

 

 

 

 

 

i

i 8(T>di

dt' = 0,5 Дt.

 

 

(32.22)

 

 

t

о

 

 

 

 

 

 

Следовательно, для коэффициента сноса справедливо следующее выражение:

a(t, х) = <?{х, * )+ - £ - <Н*, !f)'I/(x, t).

(32.23)

Нетрудно убедиться, что дополнительные члены в правой части

(32.20) не влияют

на коэффициент диффузии, а потому справед­

ливо равенство

 

 

 

b(t, х) = [ф(х, t)f.

(32.24)

Выражения (32.23)

и (32.24) называются

коэффициентами сноса

и диффузии для стохастического дифференциального уравнения

(32.1) по

Стратоповичу.

Для функций <р (*,-/)

и ф (х,

t)

при этом

справедливы следующие выражения:

 

 

 

9 (х, t) =

a (t, х) - -1

дЬ{^ / ~ : Ф(*. *) =

/ Н Г * )

(32.25)

Таким образом, если марковский случайный процесс X(t) опре­ делен стохастическим дифференциальным уравнением (32.1), то ко­ эффициенты сноса и диффузии по Стратоповичу находятся с по­ мощью соотношений (32.23) и (32.24). Если a(t, х) и b(t, х) — коэффициенты сноса и диффузии марковского случайного процесса X(t), то коэффициенты <р(х, t) и ф(х, t) стохастического дифферен­ циального уравнения (32.1) определяются формулами (32.25).

Из сравнения (32.16) с (32.23) и (32.24) следует, что в обоих рассмотренных выше случаях определения стохастического инте­ грала коэффициенты диффузии одни и те же, а коэффициенты

246


сноса

совпадают, только если функция г|з(-*, t)

не зависит от х, т. е.

когда

ф' (х, t) = 0. При коэффициентах сноса

и диффузии

(32.16)

по Ито вместо (32.25) справедливы равенства:

 

 

 

<р(х, t)=a(t, х); ф(х, t) = \fb(t, х) .

(32.26)

Рассмотрим систему автоматического регулирования с запазды­ ванием, функционирование которой описывается следующим диф­ ференциальным уравнением:

X(t +

r]) = <p[X(t),t] + q[X(t),t}Ut),

(32.27)

где использованы те же обозначения, что и в (32.1),

а г) — малое

время запаздывания

(ц > 0 и ц -* 0). При этом вместо

(32.4) будет

X (t + п +

ДО — X (t -f n) = Д^1 (t) + AZ2(t),

(32.28)

гдй AZi(t) и AZ3(t) определяются равенствами (32.5) и (32.6). Положим:

Если для AZi(t) и AZ2(t)

использовать выражения

(32.9) и (32.10),

то по аналогии с (32.16)

получим

 

 

 

 

 

an (t, x) =cp(x,

t)\

bn(t, x ) = ['[>

(x,

7)]3.

 

(32.33)

Коэффициенты сноса

a(t, х)

и

диффузии

b(t, х)

применительно

к уравнению

(32.27)

определяются как предельные значения функ­

ций а-,](/, х)

и brt(t, х)

при т)

0.

Из (32.33) следует, что при этом

справедливы равенства (32.16). Убедимся,

что

данные выражения

справедливы и в том случае,

если случайную

функцию ДZ2(t) оп­

ределить другим рассмотренным

выше способом.

Вместо

(32.19)

в соответствии с (32.28) имеем

 

 

 

 

 

о

X ( n - X ( t ) = q i X ( t - n ) , t - r \ ] ( t ' - t ) +

 

 

 

 

t ' - T )

 

 

 

 

 

+ ф [Х(г — -n), i

t]]

U t")dt"+0(t'-t).

(32 .34)

247


Так как время г| мало, то можно считать

 

 

 

X{t') - X(t) = q>[X(t),

t](t' — t) +

 

 

 

t'

 

 

 

 

 

(32.35)

 

+ 4 W ) , п | б ( ^ - л ) Л "

+

0 ( ^ - 0 + 0(л).

 

 

t

 

 

 

 

 

 

С учетом этого выражения по аналогии с (32.21) находим

 

 

a A t, х) = <?(х,

£)-fO(r()-f

 

 

 

 

t + At

 

V

 

 

+ ^ (* .

* )1, т . - Л

 

8

{t' -

t" г г,) dt" dt'.

(32-36)

 

At—*0

-*£•

t

 

t

 

 

 

 

 

 

 

 

 

 

Полагая во внутреннем интеграле

t' t" = т, получаем

 

 

t'

 

 

 

t’ - t

8 ( i +■»))<&.

 

 

o ( t ' - t " + ^)dt" =

 

j

 

 

t

 

 

 

oJ

 

 

Так как

t' — t^ 0 , то

т +

г )> 0 ,

а потому 8 (т + г]) = 0. Следова­

тельно,

a-4t, х) — у(х,

0 +

0 (Tl))

а

 

потому коэффициент

сноса

a{t, х) =

<р(х, t).

 

 

 

 

 

 

 

Таким образом, если в системе, функционирование которой опи­ сывается стохастическим дифференциальным уравнением (32.27), время запаздывания ц -> 0, то коэффициенты сноса и диффузии оп­ ределяются формулами (32.16); при этом справедливы обратные равенства (32.26)

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

248


§ 33. С Т А Ц И О Н А Р Н А Я П Л О Т Н О С Т Ь Р А С П Р Е Д Е Л Е Н И Я

Одномерная, плотность распределения f\(y, т) марковского слу­ чайного процесса является решением дифференциального уравне­ ния Колмогорова (30.24), или, что то же самое, уравнения (30.26), т. е.

<?М у; Я>

, d S ,( y ; х) _ п

(33.1)

дг

+ '

ду

 

где

 

 

 

5 , Су ; х) = а ( х , у) А (у; т ) -------У ) М у ; Т)1 •

(33.2)

В общем случае функция f\(y; т) зависит от обоих аргументов, т. е. от ординаты у и от времени т. Если с увеличением т одномерная плотность распределения f\(y, т) стремится к определенному пре­ делу f(y), то говорят, что существует стационарное решенпе урав­ нения (33.1). Функция f(y) называется стационарной плотностью распределения, причем

f(y) = Umfi(y; т).

(33.3)

Ясно, что стационарное решение уравнения (33.1)

существует не

всегда. Необходимым условием существования такого решения яв­

ляется

независимость коэффициентов

сноса а(т, у) и диффузии

Ь(т, у)

от времени т.

плотность, распределения

При

указанном условии условная

f(t, х ; т, у), являющаяся решением дифференциального уравнения (30.23), по прошествии большого промежутка времени т — t обычно не зависит не только от t и т, но и от начального состояния х. Пре­ дельное значение этой условной плотности распределения назы­ вается стационарным решением уравнения (30.23). Так как коэф­ фициенты уравнений (30.23) и (30.24) совпадают, то стационарным решением уравнения (30.23) является та же стационарная плот­ ность распределения f(y). Следовательно, если стационарное реше­ ние существует, то

f ( y ) = Hmf(/, х; т, у).

(33.4)

Т—t-~

 

От начального условия стационарная плотность распределения не зависит, а от вида граничных условий зависит.

Дифференциальное уравнение для [(у) получается из (33.1) и (33.2) при замене функций fi (у; т), а (г, у) и Ь (т, у) соответственно

на f(y), а(у) и Ь(у). Так как ^ ^ = 0, то из (33.1) следует, что

dSj

f~= 0, а потому ноток вероятности Si (у; т) постоянный, т. е. один dy

249