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