книги / Флуктуации в динамических системах под действием малых случайных возмущений
..pdfГАУССОВСКИЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ И ПОЛЯ |
141 |
Если a{zx z) > a(z2 z) для всех остальных точек z из D{JdD, то
lim pfsup |сХ2 — ф (2) |< б| еХ е |
= 1 |
при любом 6 ]> 0.
В случае однородного случайного поля Хъ уже нельзя указать функцию ср, вдоль которой с подавляющей вероят ностью при е 0 достигают уровня 1 реализации еХ£ при условии, что они его вообще достигают. Однако в этом случае функции a(z0, z) при разных z0получаются друг из друга сдвигом, и можно доказать, что
шах |
гХ |
_~ z ~~~2°^ |
1 |
|
= DU0D |
2 |
“ {0) |
||
|
||||
где б > 0, а(г) = |
a(z0, |
г0 + г). |
|
Г Л А В А 4
ГАУССОВСКИЕ ВОЗМУЩЕНИЯ ДИНАМИЧЕСКИХ СИСТЕМ.
ОКРЕСТНОСТЬ ПОЛОЖЕНИЯ РАВНОВЕСИЯ
§1. Функционал действия
Вэтой главе мы будем рассматривать возмущения ди намической системы
xt = b(xt), х0 = х, Ъ{х) = (Ь'(х), . . ., Ьг(х)) (1.1)
процессом белого шума и вообще гауссовским процессом. Функции 6*, если не оговорено противного, будем считать ограниченными и удовлетворяющими условию Липшица: |Ь(х) — Ъ(у) | К\х — 1/|, \Ь(х)\ К < оо. Основное вни мание мы уделяем здесь случаю, когда возмущенный про цесс имеет вид
X ? - b (X ? ) + e ^ f X l = x t |
(1.2) |
где wt — r-мерный винеровский процесс. В гл. 2 мы гово
рили о том, что при е -> 0 процессы X] сходятся по веро ятности равномерно на каждом конечном отрезке [О,Г] к траекториям динамической системы (1.1). На этот ре зультат можно смотреть как на вариант закона больших чисел. В гл. 2 имеется также утверждение типа центральной
предельной теоремы для процессов Х\: нормированная раз
ность (Xf — xt) сходится к гауссовскому процессу. Этот результат характеризует отклонения порядка е от предельной динамической системы. В этой главе мы изу чим асимптотику вероятностей больших (порядка 1) ук
лонений для семейства процессов Xf и рассмотрим ряд
§ 1) |
функционал действии |
143 |
задач, касающихся поведения возмущенного процесса на больших отрезках времени. В последнем параграфе мы
рассмотрим большие уклонения для процессов Х \ кото рые определяются уравнениями
|
|
|
|
|
|
|
XI = х , |
(1.3) |
|
где |
£* — гауссовский процесс в RJ\ b(x, |
у), х е RT, у е |
|||||||
е |
Л1, - непрерывная |
функция, |
для |
которой Ъ{х, 0) = |
|||||
= |
Ъ(х). |
\f)f — непрерывная функция на [0, Т\ со зна |
|||||||
|
|
Пусть |
|||||||
чениями в Rr. Рассмотрим в пространстве C0T{Rr) |
опера |
||||||||
тор |
Бх: ф |
и, |
где |
v = |
vt — решение уравнения |
|
|||
|
|
|
|
|
t |
|
|
|
|
|
|
|
Vt = |
х + |
J Ь(v8) ds + ф/, |
* ^ |
10, Т]л |
(1.4) |
|
|
|
Нетрудно |
|
о |
что при сделанных предположе |
||||
|
|
доказать, |
ниях относительно Ъ{х) решение уравнения (1.4) существу
ет и единственно для любой непрерывной функции ф и |
|
любого х е |
/Г . Оператор Вх имеет обратный |
|
t |
|
(З Г М / = Ф* = vt — X— j Ъ(vs) ds. |
|
о |
Л е м м а |
1.1. Предположим, что функция Ъ(х) удое* |
летворяет |
условию Липшица |
\Ь(х) — Ъ(у)| < К\х — у\.
Тогда оператор Вх в C0T(Rr) удовлетворяет условиюЛип шица'.
II#*Ф — |
< eKTh |
— # |
ф,Ф е |
C0T(Rr). |
Д о к а з а т е л ь с т в о . |
По |
определению опера |
||
тора Вх для и = |
Bxq>, v = Bxty имеем |
|
||
|
t |
|
|
|
|
<K -\\ut — vt \ds + 1(f — ф||, |
t e= [0 , Т]. |
||
|
0 |
|
|
|
144 ОКРЕСТНОСТЬ ПОЛОЖЕНИЯ РАВНОВЕСИЯ [ГЛ. 4
Отсюда на основании леммы 1.1 гл. 2 вытекает утвержде ние леммы 1 .1 .
На |
пространстве |
Сот (Rr) |
рассмотрим |
функционал |
||||||
5(ср) = |
S0T (ф), |
который |
для |
абсолютно |
непрерывных |
|||||
функций ф е Сот(Rr) |
определяется |
равенством |
|
|||||||
|
|
|
|
т |
|
|
|
|
|
|
|
S OT (ф) = j |
J |Ф* — Ъ(ф5) |ds; |
|
|
||||||
|
|
|
|
О |
|
|
|
|
|
|
для остальных |
ф е CoT(Rr) |
считаем |
5 0Т(ф) = 4 - 00. |
|
||||||
Т е о р е м а |
1.1. Функционал е~25(ф) является функ |
|||||||||
ционалом действия для семейства |
процессов |
X?, опреде |
||||||||
ляемых |
уравнением (1.2), |
в |
пространстве |
C0T(Rr) |
щи |
|||||
е —>■0 равномерно относительно начальной точки х е |
Rr. |
|||||||||
Д о к а з а т е л ь с т в о . |
|
При |
каждом |
фиксирован |
ном х утверждение вытекает из теоремы 3.1 гл. 3 и леммы 1 .1, если принять во внимание вид функционала действия для семейства процессов ewt.
Действительно, так как оператор Вх непрерывен в
С0Т (Rr) и имеет обратный, |
то функционал действия для |
семейства процессов X 8 = |
Bx(&w) имеет вид е~25(ф), |
т |
т |
|
°dt==Y j 1ф«— b(q>/)|2d#, |
|
О |
|
t |
если функция |
= (pt — х — [ Ъ(фв) ds абсолютно не- |
|
b |
прерывна. Ясно, что эта функция абсолютно непрерывна тогда и только тогда, когда абсолютно непрерывна функ ция ф*.
Теперь нам нужно проверить равномерность по х. Выполнение формул (3.11), (3.12) гл. 3 при малых значе
ниях параметра сразу для всех х е Rr вытекает из равно мерной по х непрерывности оператора Вх (условие Липши ца с константой, не зависящей от х). Остается проверить требование (0 К) § 3 гл. 3: полунспрерывность снизу функ ционала S0T и компактность U Ф*($) Для любого ком-
3 C G X
пакта К . Это |
выводится из того, что (х, ф) -+Bxty — |
го- |
меоморфизм, |
если рассматривать только ф с ф = |
0 . |
§ 1] |
ФУНКЦИОНАЛ ДЕЙСТВИЯ |
ш |
Отсюда, в частности, вытекает, что нижняя грань S()T{ср) по любому ограниченному замкнутому подмноже
ству Сот (Rr) достигается; причем значения, близкие к наименьшему, функционал S0T принимает только вблизи функций, на которых достигается минимум.
Отметим, что если SoT(ф) = 0, то функция ф на [О, Т| является траекторией динамической системы (1.1), так как в этом случае ф абсолютно непрерывна и почти всюду
на |0, Т ] удовлетворяет уравнению |
(pt ~ b(cpt). |
|
||||||
Рассмотрим некоторые простейшие применения тео |
||||||||
ремы 1.1. Пусть |
D — область в |
i?r, |
dD — ее граница, |
|||||
с(х) — непрерывная ограниченная функция в /?г, |
g{x) |
|||||||
ограниченная непрерывная функция, заданная па |
dD. |
|||||||
Обозначим |
|
т8 = |
min {£: Xf ф D) |
|
где Xf — решение |
|||
уравнения |
(1.6); |
HD{t, |
х) = |
{ф е |
CQT (Rr): ф0 = х, |
|||
фt s D{JdD}, |
f f D(t, |
х) = |
{ф <= Сот (Rr): ф0 = х, |
<р8ф D |
||||
при каком-то s е |
[О, |
Л). |
|
|
|
|
|
|
Т е о р е м а |
1.2. Предположим, что граница области |
|||||||
D совпадает с границей ее замыкания. Тогда |
|
|||||||
П т е2 lnP^ {X8 е О] = — |
min |
5 0г(ф), |
(1.Г>) |
|||||
е“*° |
|
|
|
фен0(/,х) |
|
|||
lim е2 lnP x [х8 <! £} = — |
min SQT{ф). |
(1.6) |
||||||
е'*° |
|
|
|
(peJJD«,*) |
|
Если экстремаль ф8, доставляющая минимум функционалу х), единственна и только при одном зна
чении s [0, t] она принимает значение из dD, то
е->° |
р * ь е< а |
|
|
= 8 (фг) охр j 1с (ф5) dsj. |
(1.7) |
|
Д о к а з а т е л ь с т в о . Так как &~2S0T (ф) |
есть |
функционал действия для семейства процессов Xf, то соотношения (1.5) и (1.6) вытекают из теоремы 3.4 гл. 3.
146 ОКРЕСТНОСТЬ ПОЛОЖЕНИЯ РАВНОВЕСИЯ (ГЛ. 4
При этом нужно еще воспользоваться регулярностью мно жеств ITD(t, х) и IID(t, х)у которая доказана в примере 3.5 гл. 3.
Чтобы доказать |
(1.7), рассмотрим |
события А\ = [т8^ |
Ро,(Хе, ф) < |
6}, Л® = |те< ^ , |
Ро*(Хв* ф )> б } . На |
множестве А\ выполняется неравенство
(1.8)
где Х&-> 0 при 6 | 0. Эта оценка вытекает из того, что время, проведенное кривой ф8 в 8-окрестности точки ФГ G= 0D, стремится к нулю при б \ 0, а после момента 5 экстремаль ф3 на 0D не попадает. Из соотношений AJ [} А\= = {т8 t) и (1.8) следует неравенство
f * (ф7)ехр | с ($ 3) ds — Хб
< |
£(ф ?) ехр |
j c ( ( p s) |
dsj-bX e P x O ^ o + U |
I W * ^ ) ’ |
|
|
|
|
(1.9) |
где |
||с|| = sup |
|с(а:)|, |
||^|| = sup ]g(x)\. Так |
как ‘ф — |
единственная |
экстремаль функционала действия на мно |
жестве HD(t, х), то значение функционала Sot на функ циях, достигающих 0D и находящихся на расстоянии не
менее б от ф, больше, чем ^ (ф ) + у, где у — какое-то положительное число. Пользуясь теоремой 1.1 (верхней оценкой с у/2 вместо у), получаем, что при достаточно ма лых е выполняется неравенство
Р* ( 4 ) < ехр [- <Г2 (So, (ф) + у12)}.
8 И |
ФУНКЦИОНАЛ ДЕЙСТВИЯ |
|
147 |
||
Отсюда, |
принимая во |
внимание соотношение |
(1.6), |
||
заключаем, что |
|
|
|
|
|
|
lim Рх {At)/Px [т * < |
*] = 0 . |
|
(1.10) |
|
|
e->Q |
|
|
|
|
Если теперь |
неравенство |
(1.9) разделить |
на Рх{тг |
£}, |
|
то, учитывая |
(1.10) и то, |
что lim |
= 0, |
получим |
(1.7). |
|
|
6 о |
|
|
|
Таким образом, вычисление главного члена логарифми |
|||||
ческой асимптотики вероятностей |
событий, связанных с |
процессом Хи свелось к решению некоторых вариацион ных задач. Эти задачи носят стандартный характер. Для экстремалей обычным образом пишется уравнение Эйлера, а сам минимум удобно находить с помощью уравнения Гамильтона — Якоби (см., например, Г е л ь ф а н д , Ф о м и н [1 ]). Если обозначить
|
|
V *1 У) = |
min |
Sot (ф), |
|
|
фо=*.фt^ V |
|
|
то |
min |
£о*(ф) = min |
v |
Х 1 У)у min Sot (ср) = |
|
q<=:HD (t,x) |
y^D\JdD |
<p^HD {Ux) |
= min V (s, x, у). Уравнение Гамильтона — Якоби для
V$D
V(t, x, у) имеет вид
4 г = Т • v »y & Х1 ^ I2 + М * Vj/^ {h у)), (1.11)
где Vy — оператор |
взятия градиента по переменным у. |
К уравнению (1.11) |
нужно еще добавить условияУТО, х. х)= |
= 0, V(t, х, у) > |
0. |
В заключение этого параграфа заметим, что функции и*(*. *)=» Р *(X ?e f i ) , ve(tt x) = Px [x* ^ t } i юв(< i*)“
представляют
148 |
|
ОКРЕСТНОСТЬ ПОЛОЖЕНИЯ РАВНОВЕСИЯ |
[ГЛ. 4 |
|||||
собой |
решения следующих |
задач: |
|
|
||||
~ |
— — Аи + (Ь(х), |
у л«е), |
^ > 0; |
|
||||
|
|
|
ие(0, х) — 1 при х е D, |
|
|
|||
|
|
|
и" (0, х) —0 при х ф D\ |
|
|
|||
— |
= |
~ |
Ai>e + |
(Ь(.г), у Л |
х <= D, |
t > 0; |
|
|
|
|
|
Vе (0, |
а;) — 0, |
г/ (t, |
х) |v.e aD = |
1; |
|
|
= |
-у |
Ди>е + |
(& (*), |
v X ) |
+ с (-О w\ xz=D, |
t > 0; |
w&(U, x) = 0,
we(<, ж) UeeD = g H ;
так что на теорему 1.2 можно смотреть как на некоторые утверждения, касающиеся поведения решений дифферен циальных уравнений с малым параметром при старших производных при стремлении этого параметра к нулю.
§ 2. Задача о выходе из области
Пусть D — ограниченная область в Rr, dD — ее гра ница, которую будем для простоты считать гладкой. Если траектория xt(x) системы (1.1), начинающаяся в точке х е D, за конечное время выходит из D (J dD, то траектории
процесса X], исходящие из х, тоже с вероятностью, близ кой к единице при малых е, покинут за это время область D, причем первый выход из области с большой вероят ностью произойдет вблизи точки выхода из D траектории xt(x) (см. гл. 2).
В этом параграфе мы будем считать, что (Ь(х), п(х))< 0 при х <= dD, где п(х) — внешняя нормаль к границе об ласти D, так что кривые xt(x) при х е D не могут покинуть
область D. Траектории процесса X], исходящие из точки х е D, и в этом случае тоже при каждом е ^ О с вероят ностью 1 выходят из D, однако точка выхода и время, необходимое для достижения границы, уже не определя ются при малых е траекторией динамической системы xt(x), а зависят, вообще говоря, от вида поля Ь{х) во всей
ЗАДАЧА О ВЫХОДЕ ИЗ ОБЛАСТИ |
149 |
области D . Здесь мы изучим задачу о выходе из области при простейшем устройстве поля Ъ(х) внутри D, совме стимом с условием (Ь(х)у п(х)) < 0. Более общий случай будет рассмотрен в гл. 6.
Пусть точка О е Дг есть асимптотически устойчивое положение равновесия системы (1.1), т. е. для всякой ок
рестности |
точки О найдется меньшая окрестность |
г |
такая, что |
траектории системы (1.1), начинающиеся |
в |
(д2, стремятся при t -> оо к пулю, не выходя из <gV
Мы говорим, что область D притягивается к О, если траектории xt(x), х е D при t -> оо стремятся к положе нию равновесия О, не покидая при этом D.
Квазипотенциалом динамической системы (1.1) от носительно точки О назовем функцию F(0, х), определен
ную |
равенством |
|
|
Отметим, что концы отрезка |
17\, Т21 не фиксируются. |
||
Смысл названия |
квазииотепццал будет выяснен в следую |
||
щем параграфе. Легко проверить, что У(0, х)^0, V(0, О) ~ |
|||
— 0 и функция |
V(Oy х) непрерывна. |
||
|
Т е о р е м а |
2.1. Пусть О — устойчивое положение |
|
равновесия системы (1.1), область D притягивается к О, |
|||
и (Ь(х), п(х)) < |
0 при х е dD. Предположим, что суще |
||
ствует единственная точка у0е |
dD, для которой V (О, у0) = |
||
= min V (О, у). Тогда для любого 6 > 0 и х G D |
|||
у^дО |
|
|
|
где |
т8 = inf{£: |
X8 е с?/?}. |
|
|
Д о к а з а т е л ь с т в о |
этой теоремы проведем по |
следующему плану. Сначала покажем (при помощи лем мы 2. 1), что, исходя из любой точки ж е й ,с вероятностью,
стремящейся к 1 при е —►0, марковские траектории Xf, прежде чем выйти на dD, попадут в малую окрестность
положения равновесия (рис. 3). Так как процесс Xet — строго марковский, то отсюда следует, что достаточно изу чить, как выходят из области траектории, начинающиеся в малой окрестности точки О.
150 |
ОКРЕСТНОСТЬ ПОЛОЖЕНИЯ РАВНОВЕСИЯ |
[ГЛ. 4 |
Пусть Г и у—две маленькие сферы радиусов ц и р/2 со ответственно с центром в положении равновесия О (рис. 4). Введем возрастающую последовательность марковских моментов т0, сг0, тх, ои т2,. . . следующим образом: т0 = 0>
an = |
inf |
{* > т„:.Х?«=Г}, т„ = inf |
{f > |
an_ x: X? |
<= |
e y (J |
dD} |
(если на каком-нибудь |
шагу |
процесс |
X® |
больше не достигает множества Г, полагаем соответствую
Рис. 3. |
Рис, 4, |
щий марковский момент и все следующие равны м и 00J заметим, впрочем, что рассмотрения бесконечных хп и оп можно избежать, если изменить подходящим образом
поле Ь(х) вне области D). Последовательность Zn = Х*п об
разует цепь Маркова на множестве y[jdD (вообще говоря, неконсервативную цепь: Zn не определено, если тп =оо; но это может быть только после выхода на dD). При малых е эта марковская цепь за один шаг с подавляющей вероят ностью переходит из любой точки х е y\JdD в множество у. Но оказывается, что если уж цепь совершит переход из # е у на 0D, то с вероятностью, стремящейся к 1 при е ->■ 0, этот переход произойдет в точку, лежащую в окрест-