![](/user_photo/_userpic.png)
книги / Флуктуации в динамических системах под действием малых случайных возмущений
..pdf§ 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 естественно выби рать различные вероятностные характеристики случайной