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

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

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

§ 91

УСРЕДНЕНИЕ В СТОХАСТИЧЕСКИХ УРАВНЕНИЯХ

361

 

Т е о р е м а 9.2 (Ф р ей д л и н [И]). Пусть

X? —

первая компонента марковского процесса Z\ на М X Ег

управляемого оператором j ? 8. Положим

т

SOT (ф) = j L((p„ q>s) ds

 

 

 

О

 

 

 

 

 

 

для абсолютно

непрерывных

функций

ср е

Cor(Af);

для

остальных

ср е

CQT{M) считаем

SnlT(<p) =

-foo.

 

Функционал е-1 50:г(ф) является функционалом действия

для семейства

процессов Xf,

^ е

[О,

Т]х

в

пространстве

СцтЩ) при е | 0.

 

 

 

 

 

 

 

Д о к а з а т е л ь с т в о .

Вместе

с

процессом

Z] =

= (Xfj y f)

рассмотрим процесс

Z\ = (Xf,

Y\) :

 

 

х? =

ь(Х?, у?),

 

 

 

 

 

 

 

Tf =

e_1^(Tf) +

e“ 1/2C(yf) wt.

 

Обозначим *(#, г/) = С-1(у) В(х, у). Процесс Zet отличается

от процесса Z\ изменением вектора переноса по перемен­ ным, по которым имеется невырожденная диффузия, так что меры, соответствующие этим процессам в пространстве траекторий, абсолютно непрерывны друг относительно друга. Учитывая это обстоятельство, для любой функции

ср : [0, Т ] ->• М и б >

0 получим

 

 

 

Px,v lpur (X 8, ф) < 8} = Мж>!/ (рот (.Vе, ф) < б,

 

 

ехр { в - ,/2j

(е {X l

Yt),

div,) - (2г)~1j |в (X*,

Vf) pds}.

 

 

 

 

 

 

(9.13)

Пусть

[0,

Т ]->

М — ступенчатая функция такая,

что роТ(ср, ф(п)) <1 1In. Для любых у, С > 0

при достаточ­

но малых б

и 1In выполняется неравенство

 

 

 

 

 

т

 

 

У_.

 

У.е), dwt) — [ (е (\|>*П)У?),

dws)

<

 

 

 

О

 

 

4е *

 

 

 

 

 

 

 

 

 

Рог (Х е3 ф) < 6j <

ехр [— Сг *} .

13 А. Д. Вентцель, М. И. Фрейдлип

362

ДИНАМИЧЕСКИЕ СИСТЕМЫ С УСРЕДНЕНИЕМ

[ГЛ. 7

Эта оценка проверяется с помощью экспоненциального неравенства Чебышева. Отсюда и из (9.13) получаем, что

для любой функции ф еС 0т(Л-/) и любого у > 0 при доста­ точно малых б и 1!п

м*,„(Рот (X е, ф) < б; exp |e- 1/2 j

?!), dws) -

- ( 2 в Г ‘j |е

Yl) 14 s - - jjJ

<

рх>!),1рог ( X е, Ф) < б )<

< М*,у |рот (X е, ф) < б; exp je~,/2 j (е (ф*п), У|), dws)

-(28)“ * j |в ( й п), Yt) |?ds + - ^ j . (9.14)

Введем в рассмотрение еще процесс Z? = (Xf, У/), который в координатной форме определяется стохасти­ ческими уравнениями

 

А 1 = Ь ( Х ? ,

У?),

 

 

 

 

 

 

 

 

 

 

Ь

=

г-'g (У?) + е~1Б № " \

у?) + е— 1/2С (У?) щ .

Из неравенства (9.14) с учетом абсолютной непрерыв­

ности мер, соответствующих процессам Z et

и Z ft вытекает,

что при достаточно малых 6 и 1!п

 

 

 

 

 

 

 

—Y/28 ^

^х,у (Рог

*

ф) < б]

 

у!2г

(9.15)

 

 

 

 

 

Px.HPorUC, Ч ')< б)< е

 

 

 

 

 

 

 

 

 

 

Пусть

tx<

t2<

...<^m- i — точки,

в

которых функ-

ция

 

\j)(fn)

имеет

скачки,

t0 =

0,

tm =

Г;

ф(,п) =

ф(/с) при

t е

\th, tk+1), к =

0,1,...,

т — 1. На

каждом полуинтер­

вале

[£ft, £ft+1) процесс X]

удовлетворяет условиям теоре­

мы

4.1.

Роль процесса

£t/e

играет процесс

Уf, который

можно

представить

в

виде

Y] =

У}/е,

У ? = ? (У }) +

+ 5

(^(/с), У() + С (У?) wt.

Выполнение условия

F выте­

§ 9] УСРЕДНЕНИЕ В СТОХАСТИЧЕСКИХ УРАВНЕНИЯХ 363

кает из лемм 4.2 и 9.1. При этом соответствующий функ­

ционал действия на абсолютно

непрерывных функциях

имеет вид

lk+i

 

 

 

 

L (ф (А*), ф8, cps) ds.

Отсюда

и

из (9.15)

вытекает,

что при любом у >

0 для достаточно

малых б

и 1In найдется такое е0 > 0,

что при е <

е0 справедливы

оценки

 

 

 

 

 

^ Рэе{Рот (Xе, ф) < б } <1

Заключительная часть доказательства этой теоремы про­ водится так же, как окончание доказательства теоремы 4.1, и мы его предоставим читателю.

З а м е ч а н и е 1. Утверждение теоремы 9.2 остается в силе, если многообразие Е заменить компактным мно­

гообразием с краем Е, а в качестве процессаZ f= (X 8, Y 8) выбрать процесс в М X Е, который во внутренних точках управляется оператором i ? 8, а при у, принадле­ жащих дЕ — краю многообразия Е — удовлетворяет не­ которым граничным условиям, например условию отра­ жения вдоль поля п(у), где п(у) — векторное поле на дЕ,

касательное к Е, но не касательное к дЕ. Изучение такого процесса можно свести к рассмотрению процесса на мно­ гообразии М X Е\ где Е' — многообразие без края,

которое получается путем склейки двух экземпляров Ё

вдоль края дЕ (см.

Ф р е й д л и н [3]).

З а м е ч а н и е

2. Нетрудно выписать функционал

действия, характеризующий уклонения порядка ех, к е=

е (0, 1/2), процесса Xf, определенного уравнениями (9.12), от усредненной системы. Пусть, например, начало координат 0 — положение равновесия усредненной систе­

мы. Тогда функционал действия для процесса Z? =*

13*

364 ДИНАМИЧЕСКИЕ СИСТЕМЫ С УСРЕДНЕНИЕМ [ГЛ. 7

= е

*Xf

имеет такой же вид, как и в теореме 7.1; в ка-

честве

 

^

(д2\х(0, 0, а) \

матрицы С нужно

взять матрицу^— да~да.— J

при

а

=

0.

 

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

мер. Пусть y f = |//е, где It — винеровский

процесс па

отрезке [—1,1]

с отражением на концах.

Процесс X®

в R T определим

дифференциальным уравнением: X? =

— b{Xzt lY^). Чтобы выписать функционал действия для

семейства процессов Xf, в соответствии с теоремой 9.2 и за­ мечанием 1 (или согласно теореме 4.1; здесь быстрое и медленное движение разделяются), нужно рассмотреть задачу на собственные значения

Nu (у) = 4~ Jji + («> ъ (*> У)) и = (У),

= 0.

дУ У=± 1

Пусть X = Х(х, а) — собственное значение оператора N с максимальной вещественной частью и Ь(х, (5) — преобра­ зование Лежандра функции Цх, а) по второму аргументу. Тогда нормированный функционал действия для процес­

сов Х* при в \ 0 на абсолютно непрерывных функциях

т

имеет вид Sor (ф) — ^ Ь (ф8, фS)ds. Чтобы найти асимптоти-

о

ку Рх,у [X* е О}, D а Дг, при г | 0 и асимптотику вероятностей целого ряда других событий, связанных с

процессом

X®, нужно, как объяснялось ранее,

вычислить

ux(t, z) =

inf{Sot(cp): ср е= Ht(x, z)},

где Ht{x,

z) — мно­

жество функций ф таких, что ф0 =

х, Ф< = z.

Начальная

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

§ 91

УСРЕДНЕНИЕ В СТОХАСТИЧЕСКИХ УРАВНЕНИЯХ

363

что в нашем случае усредненная система имеет вид

 

 

 

 

 

 

1

 

 

=

Ь(х),

Ъ(х) =

- j-

j Ь(х, у) dy,

(9.17)

 

 

 

 

 

—1

 

так

что, если

точка z

лежит

на

траектории x t системы

(9.17), для которой х0 ~ х и хи = z, то u(t0, z) = 0. Для нахождения функции u(t, z) можно воспользоваться урав­ нением Гамильтона — Ячоби. В рассматриваемом случае оно имеет вид

•i— Чг, V*")-

(9-18)

Так как N — самосопряженный полуограииченный оператор, то для его наибольшего собственного значения Jc(#, а) справедливо представление

X {х, а) = nun |-Е j [/' (i/)]2 dy J (a, b (x, y)) f2 (y) cfyj,

где S = {/ : J |/ (y) | d y= 1, f (1) = /' ( - 1) = oj. Отсюда

вытекает, что решения уравнения (9.18) удовлетворяют следующему нелинейному дифференциальному уравнению

IF (*i z) = mj n {х ] [f'

dy ~

 

 

 

~

J (VzU,

b(z, y ))/2(t/)dj/J.

Таким образом, уравнени

для

н(£, z)

можно написать,

не находя первого собственного значения оператора N. Допустим, нам нужно найти асимптотику при е 1 0

функции Vе (t, X,

у) = МХ'У}(Х * ), где f{x) :

RT i ? —

неотрицательная

гладкая функция, отличная

от нуля на

множестве G (Z Rr. Легко видеть, что vz(tx х, у) -> f(xt(x))f

3G6 ДИНАМИЧЕСКИЕ СИСТЕМЫ С УСРЕДНЕНИЕМ [ГЛ. 7

где xt(x) — решение усредненной системы (9.17) с началь­ ным условием XQ{X) = х. В частности, если xt(x) ф G,

то vE(t1х, у) -> 0 при в \ 0. Скорость сходимости ve(t, х, у)

к нулю можно оценить с

помощью теоремы о больших

уклонениях для семейства

процессов Xf :

 

lim eln ve (t} x, у)

=

— inf

ux (txz).

(9.19)

 

г 10

 

 

 

 

 

zeGUdG

 

 

 

Как известно, функция v*(t, x, у) есть решение следую­

щей

задачи:

 

 

 

 

 

 

 

 

 

ди*

1 d*ve

V

гЛ/

ч Я1*

 

 

 

 

 

dt

 

 

i—1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x<=Rr, i/<=(— 1, 1),

t> 0 ,

 

(9.20)

 

ve (0, x, y)= }(x),

di’

(tlX ly)

 

= 0.

 

 

J= ± I

 

 

 

 

 

 

 

 

 

 

 

 

Функция

v(t, x)=

f{xt(x)),

очевидно, удовлетворяет урав-

нению

 

 

 

 

 

 

 

 

 

 

 

dv

(*,*) =

2

(*) гт»

г > 0 >

 

 

 

dt

 

, Q o n

 

 

 

 

i=l

 

fa

 

 

 

 

(9.21)

 

v (0, x) =

/ (x)t

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

где

b' (x) = -E

I

b1(x, y) dy.

Поэтому из

сходимости Xf

к траекториям

—l

 

(9.17)

вытекает сходимость при

системы

6 | 0 решения задачи (9.20) к решению задачи (9.21). Соотношение (9.18) позволяет оценить скорость этой схо­ димости.

Как мы видели в гл. 4 и 6, в стационарных задачах с малым параметром большие уклонения могут определять главный член асимптотики, а не только члены, стремя­ щиеся к нулю вместе с е. Рассмотрим стационарную зада­ чу! соответствующую уравнению (9.20). Для краткости

§ 9]

УСРЕДНЕНИЕ В СТОХАСТИЧЕСКИХ УРАВНЕНИЯХ

367

будем считать г = 1:

 

 

 

 

ж е

(— 1,1),

} е ( -

1,1),

(9.22)

 

 

 

у=±i

- О,

и;®(1, */)U r+ =

Ф(1, »),

 

 

 

 

 

, У2/) U|уегr _ = ^ ( - 1 ? »),

 

 

где

Г+ =

{у е=

[—1,

1]: 6(1,

у)

>

0},

Г_

= (I/ G

[—1,

1]:

Ь(— 1,

у) <

0},

\р(1,

у) и

ф(—1,

у) — непрерывные

функции,

определенные при у е

[—1,

1 ]. ‘Решение зада­

чи (9.22), по крайней мере обобщенное, всегда существует, но без дополнительных предположений оно не единственно

(см.

Ф р е й д л и н [1 ], [3]).

Потребуем,

чтобы для каж­

дого

х

[—1, 1]

нашлось

у0 = у0(х0)

такое, что либо

Ь(х, у0) >> 0 при х ^

с0, либо 6(ж, у0) < 0

при х ^ х0, Это

условие обеспечивает единственность решения задачи

(9.22). Если Zf = (Xf, yf) — процесс в полосе |^ 1 с отражением по нормали к границе, управляемый при

\у\ <

1, х е

(—оо,

оо)

оператором

L?y

тЕ=

min {t:

I Х\ |=

1},

то единственное решение задачи

(9.22)

можно

I

 

виде

wF(x,

у) =

. g

Уте).

 

записать в

(Х те,

 

Рассмотрим теперь усредненную динамическую систе­ му (9.17) на прямой. Предположим, что траектории усред­ ненного уравнения, исходя из любой точки х е [—1, 1 ],

выходят из этого отрезка. Тогда, очевидно, функция Ь(х)

не меняет знака. Пусть для определенности

b(x) J> 0 при

х е

[— 1,

1]. Если

дополнительно предположить,

чго

6(1,

у) ^

0 при у е

[—1, 1], то нетрудно

доказать,

что

Если не предполагать, что в (1.9) 6(1, у)

0 при у е

е[—1, 1], то положение существенно усложняется. См.

по этому

поводу

С а р а ф я н , С.а фа

р ян, Ф р е й -

д л и н [1 ]. В этой

же работе рассмотрен

случай, когда

траектории

усредненного движения, не

выходят из ин­

тервала (— 1, 1).

 

 

Г Л А В А 8

УСТОЙЧИВОСТЬ ПРИ СЛУЧАЙНЫХ ВОЗМУЩЕНИЯХ

§1. Постановка задачи

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

x t = b(xt)

(1.1)

в Rr имеет в точке О положение равновесия: Ь(0) — 0. Положение равновесия О называется устойчивым (устойчивым по Ляпунову), если для всякой окрестности

Ux э

О найдется окрестность U2такая, что решения урав­

нения

(1.1) с начальными

условиями х0 = х е

U2 ни­

когда не покидают

множество

Uv при положительных t.

Если,

кроме того,

lim xt

=

О для траекторий,

исходн­

щих из точек я0 =

ое»

 

 

 

х, достаточно близких к О, то положе­

ние равновесия О называется асимптотически устойчивым. С устойчивостью относительно возмущения начальных условий тесно связана задача об устойчивости при постоян­ но действующих возмущениях. Чтобы пояснить, в чем

состоит эта задача^ рассмотрим наряду с

уравнением

(1.1) уравнение

 

= Kxt) + £,,,

(1.2)

§ 1]

ПОСТАНОВКА ЗАДАЧИ

3G9

где Z>t — некоторая ограниченная непрерывная функция

на полупрямой [0, оо) со значениями в Rr. Задача об устойчивости при постоянно действующих возмущениях может быть сформулирована так: при каких условиях на поле Ь(х) решение задачи (1.2) с начальным условием

х0 = х будет при О| +

sup

|^|

0 сходиться рав-

номерно на полупрямой

Q<.t<oo

к

постоянному ре­

[0, оо)

шению хг = О? Можно доказать, что если положение рав­ новесия устойчиво в достаточно сильном смысле относи­ тельно малых возмущений начальных условий, то оно устойчиво и при постоянно действующих возмущениях. Например, если положение равновесия равномерно асимп­ тотически устойчиво, т. е. Пш xt = О равномерно по

I—> оо

х0 — х из некоторой окрестности точки О, то это положе­ ние равновесия устойчиво и относительно малых постоян­ но действующих возмущений (см. М а л к и н [1 ]).

Положение существенно меняется, если отказаться от ограниченности функции £* на полупрямой [0, оо). В этом случае решения уравнения (1.2) могут, вообще говоря, выходить из любой окрестности положения равновесия, даже если это положение равновесия устойчиво относи­ тельно возмущений начальных условий в самом сильном смысле. Да и само понятие малости постоянно действую­ щих возмущений в этом случае требует уточнения.

Пусть теперь £* в (1.2) — случайный процесс: £< —

со). Если возмущения £*(со) равномерно малы по ве­

роятности,

т. е. sup |£t(со) |-> 0 по вероятности,

то поло-

жеиие

не

 

0 < г < о о

от

детерминированного

случая:

отличается

если точка О устойчива для системы (1.1)

в достаточно

сильном

смысле,

то

при

стремлении

по

вероятнос­

ти sup |£*(со)|

к

нулю,

а

начального

условия х0 = х

0«<ОО

 

равновесия

О имеем: Р {sup

\xt О|^

к положению

^ 6 } - > 0 д л я

любого

б >

0.

О

о о

 

 

Однако

в большинстве

задач

предположение

о том,

что £* стремится к нулю

равномерно на

всей

полупрямой

[0, оо),

слишком огра­

ничительно. Более естественно предположение о стрем­

лении к нулю sup М |^I2 или каких-либо других харак- t

теристик процесса £*, делгюцее маловероятным большое

370

УСТОЙЧИВОСТЬ ПРЛ СЛУЧАЙНЫХ ВОЗМУЩЕНИЯХ (ГЛ 8

значение |£*| в каждый фиксированный момент t, но позволяющее функции |£*(о>)| принимать большие зна­ чения в какие-то моменты времени, зависящие от со. При предположениях такого рода траектории процесса

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

ни было мало М|£*|2= а Ф 0, траектории xt свероятностыо 1 как угодно далеко отклоняются от положения рав­ новесия.

Чтобы пояснить характер тех задач, касающихся устой­ чивости при случайных возмущениях, о которых будет идти речь в этой главе, рассмотрим некоторый объект, состояние которого описывается точкой х е Rr. Предпо­ ложим, что эволюция во времени этого объекта описыва­ ется в отсутствие случайных возмущений уравнением (1.1), а случайные возмущения С*(о>) приводят к тому, что изменение состояния описывается уравнением

Xt = b(Xt, £,)•

( 1 . 3 )

Мы допускаем в качестве £* и некоторые обобщенные слу­ чайные процессы, лишь бы уравнение (1.3) имело смысл, и для любой начальной точки х0 = х е Rr существовало единственное решение уравнения (1.3). Так, мы будем рассматривать случай, когда Ь(х, у) = Ь(х) -]-о(х)у, а £* — процесс белого шума или производная от пуассоновского процесса.

Предположим, что можно указать некоторую область D d Rr такую, что пока фазовая точка, характеризую­ щая состояние нашего объекта, принадлежит Z), объект не претерпевает существенных изменений, а когда фазо­ вая точка покидает!), объект разрушается. Такую область

D мы будем называть критической областью.

Пусть X? — решение уравнения (1.3) с начальным ус­ ловием XQ = х. Введем случайную величину т* = min{£:

X? ф. D] — время до разрушения объекта. В качестве меры устойчивости системы относительно случайных воз­ мущений при критической области D естественно выби­ рать различные вероятностные характеристики случайной