книги / Флуктуации в динамических системах под действием малых случайных возмущений
..pdf§ 2) |
разложение по степеням малого параметра |
£1 |
Если функции Ь(х) и о(х) имеют ограниченные производ ные до порядка к + 1 включительно, то это стохастическое дифференциальное уравнение имеет единственное решение и определяет (к + 1) мерный марковский процесс. Как и в случае уравнения (1.1), нулевое приближение — детерми нированное движение вдоль траекторий иевозмущенной
динамической системы (1.2). Процесс Х*4) определяется из стохастического дифференциального уравнения, у ко
торого вектор переноса линейно зависит от X*1*, а коэф фициенты диффузии зависят только от t. Легко проверить (это следует, например, из того, что решение стохасти ческого дифференциального уравнения можно построить
методом последовательных приближений), что Х*4)— гауссовский процесс. Таким образом, решение уравнения (1.3) с точностью до величин порядка е2 есть гауссовский
марковский (неоднородный по времени) процесс Х*0) 4*
+ яХ\1\
Рассмотрим в качестве примера одномерное стохасти
ческое |
дифференциальное уравнение |
|
|||
|
X* = b ( X Et) + |
ewt1 |
Х% = х. |
|
|
Нулевое приближение |
Х*0) |
есть |
решение |
уравнения |
|
Х$0) = |
b (Х*0)), X(Q} == х. |
Для |
Х*(1) |
получаем |
уравнение |
Х\{) = V (Х[0)) Х(*0 + wt.
Если считать Х*0) известной функцией, то решение этого уравнения можно записать в виде
X(t° = J exp Jj* V (X (u0)) duj dws.
Таким образом, |
( t |
\ |
t |
X) = X\0) + e | exp IJ У (X£) du\ dws + о (e).
Сформулируем теорему относительно разложения по степеням малого параметра е решения стохастического дифференциального уравнения.
as |
ВОЗМУЩЕНИЯ НА КОНЕЧНОМ ОТРЕЗКЕ ВРЕМЕНИ |
(ГЛ. 2 |
||||||||
|
Т е о р е м а |
2.2. |
Предположим, |
что |
коэффициенты |
|||||
Ь1(х)иа)(х) |
имеют |
ограниченные |
частные |
производные |
||||||
до |
порядка k + |
i включительно. |
|
(1.3) |
|
|
|
|||
|
Тогда для |
решения X] |
уравнения |
имеет место |
||||||
разложение |
|
|
|
|
|
|
|
|
|
|
|
ХЧ = |
Х<0) + eXjl) + |
. . . + е*Х(» + |
|
(<), (2.7) |
|||||
где |
Х((0), Х)1), . . |
X\h) |
— определяются |
|
|
из |
уравне |
|||
ний (2.6). |
Случайный процесс Х *+1 = |
(Х*0), |
Х(/°) |
|||||||
определяется первыми k + |
1 уравнениями системы (2 .6). |
|||||||||
С точностью до величин порядка еа процесс Х\ |
прибли |
|||||||||
жается гауссовским процессом Х*0) + еХ*1*. |
Остаточный |
|||||||||
член в (2.7) удовлетворяет неравенству |
|
|
|
|||||||
|
sup |
(м | д ;+ !(*)|2) |/8 < с в к+‘, |
с < |
оо. |
|
|||||
|
Д о к а з а т е л ь с т в о |
этой теоремы |
лишь |
некото |
рыми техническими моментами отличается от доказа тельства теоремы 2.1. Эти технические моменты, впрочем, требуют довольно громоздких выкладок, связанных с доказательством дифференцируемости решения стохасти ческого дифференциального уравнения по параметру. Такие вопросы лежат в стороне от нашей основной темы, поэтому мы не будем приводить здесь доказательства тео
ремы |
2.5. Это |
доказательство |
можно извлечь из |
ста |
|
тей Б л а г о в е щ е н с к о г о и |
Ф р е й д л и н а |
11]1 |
|||
Б л а г о в е щ е н с к о г о |
[1]. |
|
|
||
Из |
разложения |
|
|
|
|
|
ХЧ = |
Х'(0) + еХ<° + |
. . . + |
thX[k) + о (е*) |
(2.8) |
для реализаций X* легко получить разложение по степе ням е для гладких функционалов от реализаций. Пусть функционал F дифференцируем по Фреше в точке Х<°>; производная F в этой точке — линейный функционал F'(X(°); h). В этом случае имеем
F(Xe) = F(X<°>) + еГ(Х<°>; Х<г>) + о(е) |
(2.9) |
при 8 -> Ох где о(г) понимается в том же смыслеАчто и в
§ 21 |
РАЗЛОЖЕНИЕ ПО СТЕПЕНЯМ МАЛОГО |
ПАРАМЕТРА |
83 |
(2.8) |
(равномерно по f G Ю, Т ] при почти всех со |
или в |
|
смысле сходимости по вероятности равномерно по t). |
|||
Из |
разложения (2.9) в случае, когда |
— гауссов |
ский случайный процесс, получаем, что значение функцио нала F(Xе) асимптотически нормально с дисперсией, пропорциональной е2. Коэффициент пропорциональности в этой асимптотической дисперсии выражается через про изводную F'(X<°>; h) и корреляционную функцию слу
чайного |
процесса |
X*1* |
(асимптотическое |
среднее |
равно |
|||||
F(X<0)), если |
только МХ*1} = 0 ). |
|
|
|
|
|||||
|
Рассмотрим теперь случай, когда функционал F диф |
|||||||||
ференцируем |
два |
раза, |
вторая производная — билиней |
|||||||
ный функционал |
F "(X (0); hXJ h2). Получаем |
разложение |
||||||||
с точностью до о(е2) при в |
0 : |
|
|
|
|
|||||
F (X е) = |
F (X (0)) + EF' (Х (0); Х(1)) + |
|
|
|
|
|||||
+ |
е2 [ 1 F" (Х (0); X(U, Х(1)) + F' (Х (0'; Х(2))] + |
о (е2); |
(2.10) |
|||||||
п т. Д. |
|
|
|
|
|
|
|
|
|
|
|
Например, для функционала F вида |
|
|
|
|
|||||
|
|
|
|
F (ф) = |
{ g (ф() dt |
|
|
|
(2.11) |
|
формулы |
(2.9), (2.10) принимают вид (для |
простоты обо |
||||||||
значений |
рассматриваем |
одномерный случай): |
|
|
||||||
1 |
1 |
|
|
|
|
1 |
|
|
|
|
j |
g (X?) dt = j g ( X H dt + |
e J *' (Х Г ) Xl1^ |
|
+ 0 |
(6); |
(2.12) |
||||
0 |
|
0 |
|
|
0 |
|
|
|
|
|
Jg (x?)dt = |
Jg (X)0)) dt + |
e fg' (X (t0)) X\l)dt + |
|
|
||||||
|
+ e2 i |
i r |
m |
(x i'o 2 dt+ j g’ ( x n x\»dt |
+ |
о(г2). |
||||
|
L |
о |
|
|
|
|
|
|
|
(2.13) |
|
|
|
|
|
|
|
|
|
|
В задачах, связанных со случайными процессами, час то приходится рассматривать функционалы^ определяемые
84 |
ВОЗМУЩЕНИЯ НА КОНЕЧНОМ ОТРЕЗКЕ ВРЕМЕНИ |
ТТЛ. 2 |
|
через первый момент выхода из области D. В случае об |
|||
ласти |
с |
гладкой границей функционал т(ф) = min {t: |
|
cpt ^ D) |
не будет дифференцируем по Фреше во всех |
||
точках |
ср |
пространства непрерывных функций; он даже |
|
не будет |
непрерывен. Однако он будет дифференцируем |
||
в тех |
точках ф, для которых ф4 имеет производную |
при |
|
t = т(ф), |
направление которой не касательно к границе. |
Не будем доказывать этого на языке производных в функ циональном пространстве, а сформулируем результат сра зу на языке разложений по степеням е.
Т е о р е м а 2.3. Пусть имеет место разложение (2.8) с к = 1. Пусть t0 — первый момент выхода траектории
Х\0) из области Z), те — первый момент выхода из этой области случайного процесса Х\. Пусть граница dD об ласти один раз дифференцируема в точке Х(^\ причем п —
вектор внешней нормали в этой точке. Пусть (X IV) > о . Тогда
|
|
( V |
- '0 |
+ |
о(е), |
(2-14) |
|
* о - е , ^ |
-4 |
||||
|
|
№ ») |
|
|
|
|
л |
___ v ( 0 ) |
е |
Х'П),-ш |
— |
V(у{ 0> № " ) |
+ о(г) (2.15) |
8 — Л.и + |
Atо |
Аи |
||||
|
|
|
|
|
(4 V J |
|
при 8 - » 0 |
(где о(г) |
понимается в смысле |
сходимости с |
вероятностью 1 или сходимости по вероятности в зависи мости от того, как понимается о(г) в разложении (2 .8)).
Д о к а з а т е л ь с т в о . |
Пользуемся разложением |
|
(2.8) |
на отрезке Ю, Г], где |
Т > t0. Прежде всего полу |
чаем, |
что тЕ-> t0 при е -> 0. |
Отсюда |
+sX^g* "j" О (б) =
=+ (т8 - О + о (Xе - g + eA'i*' + О(е). (2.16) Умножив (2.16) скалярно на /г, получим
( Х г% - Х%\ п) = (т8 - t0) (Х?,\ п) + о(т8 - О +
+ e(X};\i*) + 0 (e). (2.17)
§ 2] |
РАЗЛОЖЕНИЕ ПО СТЕПЕНЯМ МАЛОГО ПАРАМЕТРА |
85 |
С другой стороны, из-за гладкости 3D в точке Xi°\ скалярное произведение в левой части (2.17) будет беско
нечно |
малым по |
сравнению |
с Хг8е — |
Отсюда |
и из |
(2.16) |
получаем |
|
|
|
|
|
(Х'е - |
, п) = |
о (те - О |
+ о (е). |
(2.18) |
Из (2.17), (2.18) получаем разложение (2.14) для т8. Под ставляя его еще раз в (2.16), получаем (2.15).
Коэффициент при е в разложении (2.15) получается
проектированием вектора Х*^ параллельно Х*в0) на ка сательную в точке Х*о0).
Если имеет место разложение (2.8) с к = 2, причем функция Х\0) дифференцируема дважды, а случайная функция Х[{) — один раз, мы можем получить разложе ние те и Х8е с точностью до о(е2) (хотя соответствующий
функционал дважды дифференцируем только вдоль не которого подпространства). Если же случайная функция
х\» не дифференцируема (как это будет в случае диффу зионных процессов с малой диффузией, рассмотренном
в теореме 2 .2), |
разложение |
для т8 с точностью до о(е2) |
не получается. |
Поясним, с |
чем это связано. |
Дело в том, что в доказательстве теоремы 2.3 мы ни как не использовали, что т8 — это именно первый момент достижения границы, а только тог что это — какой-то
момент, когда Xf находится на границе, причем сходя щийся к t0. Если мы рассмотрим процесс X? простейшего вида: X® = х0 + t + ewt1 то его первый момент т8 дости жения точки Xi > х0 и последний момент о8 пребывания в точке хх отличаются на величину порядка е2. Действи тельно, с помощью строго марковского свойства относи
тельно марковского момента т8 получаем, |
что распределе |
||||||
ние |
а8 — х8 — такое |
же, как для |
случайной |
величины |
|||
С8 = |
max {t : t + ewt = 0}. Далее |
пользуемся |
тем, что |
||||
e~2(te2 + еш*8«) =* t + |
E-Xwtz' = t + |
wu |
где |
wt — снова |
|||
винеровский |
процесс, |
выходящий |
из нуля; |
и |
£8 = е2£, |
||
где |
£ = max |
{t: t + |
wt =* 0 }. |
|
|
|
|
83 |
ВОЗМУЩЕНИЯ НА КОНЕЧНОМ ОТРЕЗКЕ ВРЕМЕНЯ 1ГД. 2 |
§ 3. Эллиптические и параболические дифференциальные уравнения
смалым параметром при старших производных
Втеории дифференциальных уравнений эллиптического
ппараболического типов большое внимание уделяется
исследованию поведения при е -*■ 0 решений краевых
dv?’
задач для уравнений вида Ьвиъ+с{х)ив= {{х), — =
= LBvB+ c{x)vB+ g(x), где LB— эллиптический диф ференциальный оператор с малым параметром при старших производных:
L ~~~2 JU |
(.г) |
д* |
+ 2 ^ |
д |
д х 'д х ' |
д х i |
|||
|
|
|
i=1 |
|
Как было сказано в гл. 1, с каждым таким оператором LB (коэффициенты которого мы предполагаем достаточно
регулярными) связывается диффузионный процесс Xf’*; его можно задать с помощью стохастического уравнения
А’?'* = Ъ(X?**) + |
еа (Х?'х) wt, XS'* = х, |
(3.1) |
где о{х)о*(х) = (aVl{x)), |
Ь{х) = (6*(х), . . .,Ьг(х)). |
Для |
этого процесса мы будем иногда пользоваться обозначе
нием |
Х%’х, иногда Xf (х) (в |
рамках понятия |
марков |
ского семейства); а иногда — |
Xf, и тогда мы будем писать |
||
индекс х у вероятности: марковский процесс |
(Xf, Рх). |
||
В двух предыдущих параграфах этой главымы получили |
|||
ряд |
результатов относительно поведения при |
е -> 0 ре |
шений Xf * (со) уравнений (3.1). Так как решения краевых задач для оператора Ьв можно записать в виде средних значений некоторых функционалов от траекторий семей
ства (х?-*, Р), то, зная поведение Xf,x(co) при е 0, можно сделать те или иные выводы относительно поведения при е -^ 0 решений краевых задач. Этим вопросам и посвящен настоящий параграф.
Рассмотрим задачу Коши |
|
|
9V а ’-~ = |
х) + с (х) ve (t, х) + g (ж); |
* > 0, |
i>e(0,*) = |
/(*) |
(3.2) |
§ 3] |
|
УРАВНЕНИЙ С МАЛЫМ ПАРАМЕТРОМ |
87 |
||
при е > |
0 и вместе с ней задачу для оператора |
первого |
|||
порядка, |
который получается при |
е = |
0 : |
|
|
dv° (*< х) _ |
£ о^о |
t > |
02 х е |
Ягд |
|
|
= |
L°v° + с (х) и0+ g (х); |
i>° (0 , х) = f (я).
Мы |
считаем, |
что выполняются |
следующие условия: |
|
1) функции с(х), g(x) равномерно непрерывны и огра |
||||
ничены |
при х е |
Rr\ |
|
|
2 ) коэффициенты оператора L1 удовлетворяют условию |
||||
Липшица; |
Г |
|
|
|
3) к 2 2 |
|
при любых ве |
||
2 |
|
|||
щественных Хц 1 2, , , ,, |
и х е |
Rr, где к2 — положи |
||
тельная |
константа. |
|
|
|
При этих условиях решения задач (3.2) и (3.3) сущест |
||||
вуют и |
единственны. |
|
|
|
Все результаты этого параграфа сохраняются и в слу |
||||
чае, когда форма |
|
только неотрицательно опреде |
лена, но в случае вырождений нужно уточнить постановку краевых задач и ввести понятие обобщенного решения. Уточнения, необходимые при наличии вырождений, мы сделаем после разбора невырожденного случая.
Т е о р е м а 3.1. Если условия 1)—3) выполнены, то при
любой ограниченной "непрерывной |
начальной функции j(x), |
х е Rr, существует lim ye(f, х) = |
v°(t,x). Функция i>°(t, х) |
является решением задачи (3.3). |
заметим прежде всего, |
Для д о к а з а т е л ь с т в а |
что если выполнено условие 3), то существует матрица о(х) с удовлетворяющими условию Липшица элементами, для которой о(х)о*(х) = (а^(х)) (см. § 5 гл. 1)
Решение уравнения (3.2) может быть представлено следующим образом:
ve (t} *) = М / (XVх) е
88 |
ВОЗМУЩЕНИЯ На КОНЕЧНОМ ОТРЕЗКЕ ВРЕМЕНИ |
[ГЛ. 2 |
где XVх — марковское семейство, построенное с помощью уравнения (3.1). Из теоремы 1.2 следует, что процессы
Х*'х (со) сходятся по вероятности при е -> 0 равномерно на отрезке [0, t] к Х°’х— решению уравнения (1.2) с на
чальным условием X Q,X = х, Учитывая, что под знаком математического ожидания в (3.4) стоит непрерывный огра
ниченный функционал от Х®,х (со),по теореме Лебега о пре дельном переходе заключаем, что
Jc( * “ ’* ) * |
I |
(c (* 2 '* )d ,. |
О
Легко проверить подстановкой, что функция, стоящая в правой части последнего равенства есть решение задачи (3.3). Теорема 3.1 доказана.
Если предположить, что коэффициенты оператора L® имеют ограниченные производные до порядка к-\-1 вклю чительно, то матрицу о(х) можно выбрать так, чтобы ее элементы тоже имели к-\-1 ограниченную производную.
В этом случае в силу теоремы 2.2 для Xf,x можно написать разложение по степеням е до поряцка к. Если функции
/(#), с(х), |
g(x) |
имеют к-\-1 ограниченную производную, |
||||||
то, |
как |
это |
следует |
из (2.7), имеет место разложение |
||||
по |
степеням 8 до порядка к с |
остаточным членом поряд |
||||||
ка еЬ+Ч |
|
|
|
с(х) = 0 и г = 1, то реше |
||||
|
Так, например, если g(x) s = |
|||||||
ние задачи (3.2) можно записать в виде |
|
|
||||||
ve (ttx) = |
M j |
(X?) = |
|
|
|
|
|
|
|
= Мж/ (Х<{0) + еХ)” + . . . + |
е*Х,(М + |
Щ+1 (*)) = |
|
||||
|
|
|
|
|
к |
|
|
|
|
|
|
|
= |
S |
e‘M,G, + 0(8*+‘), |
(3.5) |
|
где |
Х(0), Xj^, |
Xfth) — коэффициенты |
разложения |
Xf |
||||
по |
степеням |
малого |
параметра, о которых шла |
речь |
§ 3] |
УРАВНЕНИЯ С |
МАЛЫМ ПАРАМЕТРОМ |
89 |
в тео{ еме 2 .2 ; |
|
|
|
б, = б, (Л’(А |
. . Х<п) = 1 |
^ / (Х Г + eXj0 + . . . |
|
|
|
•. . + |
8ftX<ft))|e=0. |
Из формулы (3.5) и уравнений, определяющих процессы
Х*г), можно вывести, что коэффициенты при нечетных сте пенях е обращаются в нуль. Коэффициенты при е2т яв ляются решениями некоторых уравнений в частных произ водных первого порядка; их можно, разумеется, найти также, решая системы обыкновенных дифференциальных уравнений.
Продемонстрируем вероятностный метод получения коэффициентов разложения ve(t, х) в простейшем случае— для размерности 1 — и разложения до членов порядка е2. Для коэффициентов разложения по степеям 8 решения стохастического уравнения
|
XI =» Ь(X?) + |
еа (X?) wu |
Х% = |
(3.6) |
|
выпишем первые три уравнения (2 .6): |
|
|
|||
Х(,0)= |
Ь (х (,0)), |
Х (о0) = х\ |
|
(3 .7) |
|
X\l) = |
V ( х И |
Х4° + |
а (х !0)) щ, |
Х(оП = 0, |
(3.8) |
х !2) = |
v (х Г ) |
х[2) + |
1 ъ" (х<(0)) ( х (4п ) 2 + |
|
|
|
+ а' (Х (,0)) |
Х(0 ) = 0. |
(3.9) |
Функция Х*0) — неслучайная, |
другое обозначение |
для |
||
нее — xt(x). |
|
|
|
|
Если / — дважды непрерывно дифференцируемая функ |
||||
ция, |
имеет место |
разложение |
|
|
/ (X?) |
= / (Х (0)) + |
гг (х Г ) |
+ |
|
+ |
е2|/' (X (t0)) X (t2) + 1 /" (Х<01) « ‘ >)2| + о (е2). |
(3.10) |
90 |
|
ВОЗМУЩЕНИЯ НА КОНЕЧНОМ ОТРЕЗКЕ ВРЕМЕНИ |
[ГЛ. 3 |
||||
Возьмем |
от |
обеих |
частей математическое ожидание |
||||
Vе ( t , |
X ) |
= |
М J |
( 4 ) = |
/ (Х ( (х )) + |
8 / ' (.Г, (X )) \Ax X'i + |
|
+ |
е2 |
Г (*» (X)) M.v.Al2) + | г |
(*, (*)) Мл ( а ;11)2] + |
О(е2). |
|||
|
|
|
|
|
|
|
( 3 .1 1 ) |
Процесс А*1*— гауссовский с нулевым средним, так что коэффициент при е обращается в нуль. Чтобы найти
Мх применяем формулу (3.8) и формулу Ито
^ -(А {0 ) 2 = 2V (А((0)) (А )0 ) 2 + 2о ( А(0)) 6' ( А |
+ |
+ о(А<,0))2. (3.12)
Беря математическое ожидание от обеих частей, получаем линейное неоднородное дифференциальное уравнение для
м* (А(1}) 2
|
д - М ( А * / ') 2 = 2 6 ' ( x t ( х ) ) М а ( А (, п ) 2 + а ( * , ( * ) ) |
( 3 .1 3 ) |
|
с начальным условием Mx.(X (0I)) 2 = |
0 . |
|
|
Решив это уравнение, находим и решение уравнения* |
|||
для |
Мх Х\г\ получаемое взятием математического ожида |
||
ния |
от (3.9): |
|
|
^ М хХ<(2) - И М * ) ) М яХ,< ,-(- |
|
|
|
|
+ y b "(* t (*))M ,(A l1))2, |
М,А(,2) = 0. |
(3.14) |
Итак, для нахождения коэффициентов в разложении vE(t, х) по е до порядка 2 достаточно решить нелинейное
уравнение xt(x) = b{xt(x)), и два линейных — (3.13) и (3.14).
Тот же результат проще получить стандартными ме тодами теории дифференциальных уравнений; однако ве роятностные методы могут быть применены также к менее стандартным асимптотическим задачам. Пусть, например,