Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Флуктуации в динамических системах под действием малых случайных возмущений

..pdf
Скачиваний:
4
Добавлен:
12.11.2023
Размер:
21.91 Mб
Скачать

§ 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 :

 

d(*< х) _

£ о^о

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).

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