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

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

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

§ 5]

ЗАДАЧА О ВЫХОДЕ ИЗ ОБЛАСТИ

261

Траектории динамической системы, выходящие из то­ чек х из левой части области D (левее сепаратрис, входя­ щих в точку if4), за исключением неустойчивого положе­ ния равновесия Kv притягиваются к циклу К2; точки пра­ вой части притягиваются к К3, а точки на линии разде­ ла — к К±. Значит, из точек левой половины D при малых

е процесс Х\ будет выходить на dD в малой окрестности У2, из точек правой половины — вблизи У3, а из точек разделительной линии — вблизи У2 или Из­

боли мы увеличим VD {Kt, dD) так, чтобы VD{K2, dD)r=

= 16, VD(K3, dD) = 10, VD (tf4, 0D) = 9, то { д и ­ графов, минимизирующих сумму (5.3), тоже два: Кх->■ К2,

К2->- К31 К3-+ dD,

 

К2 или К±-+-К3, для всех i имеем

M(i) = {3}. Значит, из всех точек области при

малых

е

выход на границу будет осуществляться вблизи У3.

 

Обратимся к вопросу о времени, проводимом процес­

сом X* в D до выхода на границу.

 

 

 

. .

Кz,

Рассмотрим графы

на множестве L = {/£1,

х, dD). Положим

 

 

 

 

 

 

 

 

 

 

M D(*)=-

 

min

2

^ D («,

P).

 

(5.8)

7

 

 

 

g<=G(x+ldD})

< a - > F ) e g

 

 

 

4

Обозначение G(a

 

W) введено в § 3 после формулировки

леммы 3.4.

5.3. Минимум (5.8)

можно записать также

Л е м м а

CAedywupiM

образом:

 

 

 

 

 

 

 

 

MD (х) =

min

2

Уо(а, Р);

 

 

(5.9)

 

g^G(xMdD)) ( a - P ) G = £

 

 

 

 

 

 

Л/ D И

= WD Д min [FD («А Kt) +

MD(#*)],

 

(5.10)

 

 

 

 

i

 

 

 

 

 

 

 

ede MD(Kt)

onpedeAnemcH как

м иним ум

 

 

 

 

MD (Kt) =

 

min

 

2

FD(a, P),

 

(5.11)

 

 

 

 

 

(a - > p ) € = £

 

 

 

 

в кот ором

участ вуют

графы

на м нож ест ве

{Kv . . .

. . ., Кх, d D ) (a

WD onpedeAnemcH

ф орм улой

(5.3)).

(5.10)

При нaxoжdeнuu минимумов

(5.9)

или

(5.11),

можно выбросить из рассмотрения все неустойчивые ком­ пакты Kt.

262

ВОЗМУЩЕНИЯ НА БОЛЬШИХ ОТРЕЗКАХ ВРЕМЕНИ 1ГЛ. б

Д о к а з а т е л ь с т в о опять аналогично доказатель­ ству лемм 4.1, 4.3.

Т е о р е м а 5.3. Равномерно по х в пределах каждого компактного подмножества F области D

lim е2 In М*хе -

W D - MD (x).

(5.12)

e->0

 

 

 

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

Выбираем у', p0,

plt

p2,

как при доказательстве предыдущей теоремы,

но

еще

с условием, чтобы среднее время выхода из р0-окрестно-

стей компактов K t не превосходило ехр {е~2у '}

(см. лем­

му

1.7);

рассматриваем

марковские моменты

т0(=0),

т1?

т2, . .

и цеиъХп

Обозначим через v номер

шага, на котором эта цепь впервые выходит на dD, т. е. первое п, для которого те = тп. Пользуясь строго мар­ ковским свойством, мы можем написать:

М^те = S MUZn<£dD;

М|втЛ.

Леммы 1.7, 1.8, 1.9 дают нам, что

М !^ в этой сумме

не превосходит (при малых е) 2 ехр {е -2*у'}, но больше

ехр {—е~2-у'}. Итак,

с точностью до множителя, заклю­

ченного в пределах

[2 ехр {е -2^ '}!* 1, М* те совпадает

оо

M£v. Но это математическое ожида-

с 2 Рl { Z n ^ dD} =

п—О ние оценивается с помощью леммы 3.4.

Сначала для ^ e U Gt получаем1

 

пользуясь оценками

(2.3) — (2.5),

г

е:

 

 

при малых

 

 

ехр U~2{WD -

MD {К,) - (4г + 1) у')\ <

< Млетх < 2 ехр (е“ 2 (W D -

MD {К,) +

(4' + l) у % (5.13)

Затем для начальной точки х е

i

 

 

 

мЬе= мгл + 2 м*1zie=dgti м|,т£}.

4=1

§

5]

ЗАДАЧА О ВЫХОДЕ ИЗ ОБЛАСТИ

263

Учитывая

то, что М/гх ^

2ехр (в

2у'1, оценки

(2.8), (2.10)

и

(5.13),

получаем, что

М/с8

заключено в

пределах

exp {&~2{WD—MD(X) =т=у)}. Это справедливо и для ^ G U Gt.

Так

как у > 0

 

i

произвольно мало, получаем утвержде­

ние

теоремы.

 

 

Вернемся к

рассмотренному

примеру (с VD(Kt, Kj),

заданными матрицей (2.1), и

VD(K2, dD) == 8, F^ (К3,

dD) = 2, VD(Kv dD) = 1). Вычислим асимптотику мате­ матического ожидания времени тЕ выхода на границу для траекторий, начинающихся в положении устойчивого рав­

новесия К3. Находим: WD =

10 (минимум (5.3) достига­

ется на двух графах

Кх->• К2, К2

dD, К3-> dD, КА

-> К2 или /£4->- К3);

Мп(К3) = 6 (минимум (5.11) дости­

гается на

графах

Ко,

/v4,

К2; К ^ К2,

КВ- * К 2,

К2-+ К 2\ К г + К 2,

К3- > К 2, К ^ К 3). Итак,

математическое ожидание времени выхода логарифмиче­ ски эквивалентно exp {e -2(tFD — Л^ГХ>(АГ3))} = ехр{4е-2}.

Рассматривая цепь Z n, мы можем понять, за счет чего получается такое среднее время выхода. Начинаясь в /i3, эта цепь при малых е с вероятностью, близкой к 1, прове­ дет на dg3 число шагов порядка exp {&~2VD (К3, dD)} = = exp {2е~2}, затратив время того же порядка; после

этого

она с вероятностью, близкой к 1,

выйдет

на dD,

а с

вероятностью

порядка

ехр{—&~2(VD(K3,

К2)

VD(KA4 dD))} =

ехр{—4е~2}

перейдет

к устойчивому

циклу К2 (может быть, задержавшись на относительно небольшое число шагов вблизи неустойчивого положения

равновесия К4). Если

это произошло, то

Znс большой ве­

роятностью

в

течение времени порядка

ехр{е~2*VD(K2,

dD)} =

exp

{8е~2}

будет

совершать

переходы

в пре­

делах

dg2

и

(примерно

в ехр{е~2

VD(K2,

Кг)} =

= ехр {е~2} раз реже) dgv а после этого выйдет на грани­ цу. Итак, математическое ожидание порядка ехр{4е~2} образуется за счет маловероятных — с вероятностью по­ рядка ехр {—4е~2} — значений порядка ехр {8е~2}.

Вспомним, что в случае области, притягивающейся к одному устойчивому положению равновесия, среднее

время выхода М/с8 имеет тот же порядок, что границы диапазона наиболее вероятных значений т8 (теорема 4.2 гл. 4); в частностщ любая квантиль распределения т

264 ВОЗМУЩЕНИЯ НА БОЛЬШИХ ОТРЕЗКАХ ВРЕМЕНИ [ГЛ. 6

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

Это сильно снижает ценность доказанной теоремы как результата, характеризующего предельное поведение рас­

пределения

те.

 

 

результат, сформулиро­

Приведем

соответствующий

ванный на языке дифференциальных уравнений.

Т е о р е м а

5.4.

Пусть g(x) положительная не­

прерывная функция в D (J dD;

vz{x) — решение уравнения

{х) = g{x)

в D

с нулевыми граничными условиями

на dD. Тогда равномерно по х в пределах любого компакт­

ного подмножества при г

О

V8 (X) X exp {е-2

(WD — MD(я))}.

§ 6. Разбиение на циклы. Субпредельные распределения

В задачах, которым посвящена эта глава, имеется два больших параметра: е~2 и t — время, в течение которого рассматривается возмущенная динамическая система. Естественно исследовать, что будет, если стремление этих параметров к бесконечности согласовано тем или иным образом. Будем интересоваться предельным поведением

мер Р £ {Х ^ еГ }; ограничимся случаем компактного многообразия (как в § 4).

Наиболее простой случай — когда сначала устремляет­ ся к бесконечности с”2, а затем t. Здесь все определяется поведением невозмущенпой динамической системы. Ясно,

что lim lim Р* {X] G T) = 1, если открытое множество

f - ю о Е—>0

Г целиком содержит (о-пределыюе множество траектории xt(x), начинающейся в точке х0(х) = х, и этот предел ра­ вен нулю, если Г находится на положительном расстоянии от этого множества.

В § 4 рассмотрен случай, когда сначала к бесконечно­ сти устремляется t, а затем е-2. Теорема 4.1 дает воз-

§ 61 ЦИКЛЫ. СУБПРЕДЕЛЬНЫЕ РАСПРЕДЕЛЕНИЯ 265

можпость

установить,

что lim lim P* {Xf е Г) = 1 для

открытых

множеств Г,

£—*0 /-»оо

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

на которых достигается минимум W(Kt). В случае общего положения этот минимум достигается на одном компак­ те; если к тому же на этом компакте сосредоточена единст­ венная нормированная инвариантная мера ц0 динами­

ческой системы, то lim lim Р* (Xf е Г} = р0 (Г) для всех £~>С t~*ОО

Г с нулевой (10-мерой границы.

Займемся исследованием поведения процесса Х\ на отрезках времени длины Z(e~2), где £(е“ 2) — функция, монотонно возрастающая при возрастании е-2. Ясно, что если £(е“ 2) растет достаточно медленно, то за это вре­

мя траектория процесса Х\ не успеет отойти далеко от устойчивого компакта, в области притяжения которого находится начальная точка. За большие отрезки времени осуществится переход из окрестности этого компакта в окрестности других, сначала к «ближайшему» (в смысле функционала действия) компакту, затем — к все более и более «далеким». Установим прежде всего, в каком по­

рядке процесс

Xet будет

достигать

окрестностей

различ­

ных компактов

Ki.

 

2, . . ., Z};

Q ка­

Т е о р е м а

6.1. Пусть L = {1,

кое-то подмножество L. Рассмотрим момент тQ первого

достижения процессом

Х\ границ

dgj р-окрестностей

gj компактов Kj с номерами из L \ Q. Пусть процесс начинается в giljdgi, i е Q. Тогда при достаточно малых

р с вероятностью, стремящейся к 1 при 8 -^ 0 , X^ бу­

дет принадлежать одному из множеств dgj таких, что в одном из (L \ 0)-графов g, на которых достигается

A (Q) = min

2

V(Km,K n),

(6.1)

geG(L\Q) (m->n)€Eg

 

цепочка стрелок ведет из i в ] е

L \ Q.

помо­

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

легко проводится при

щи лемм 2.1 и 3.3.

 

 

 

В формулировке теоремы можно также V(Km, К п) заменить на V(Km, К п), и в этом случае мы можем также выбросить все неустойчивые компакты и рассматривать

266

ВОЗМУЩЕНИЯ НА БОЛЬШИХ ОТРЕЗКАХ ВРЕМЕНИ [ГЛ. 6

только переходы от одного устойчивого компакта к дру­ гому.

Рассмотрим пример. Пусть K h i = 1, 2, 3, 4, 5,— устойчивые компакты, содержащие со-предельные мно­ жества; значения V(Ki4 Kj) задаются матрицей

О

4

9 13

12\

7

0

5 10

11 \

6

8

0 17

15

1

(3

6

8 0

2

 

5

7

10 3

 

0/

Пусть процесс начинается вблизи Кх. Находим {2, 3, 4, 5}-граф, минимизирующий сумму значений V. Этот граф состоит из одной стрелки 1 -> 2; значит, первым из

других компактов,

к которому подойдет процесс, бу­

дет К2.

куда мы пойдем из окрестности К2.

Посмотрим далее,

Полагаем Q = {2}. Находя минимизирующий сумму (6.1) {1, 3, 4, 5}-граф 2 3, получаем, что следующим компак­ том, к которому приблизится процесс, будет с подавляю­ щей вероятностью К3. Далее граф 3 -> 1 показывает, что процесс вернется к компакту Кх. Переходы от Кх к К2ч от К2к К3 и от К3 обратно к Кх будут происходить много раз, но в конце концов процесс подойдет к одному из ком­ пактов А 4, К ъ. Чтобы понять, к какому именно, пользуем­

ся теоремой 6.1 для Q

{1 ,2 , 3}.

Находим

{4, 5}-граф,

на

котором достигается

минимум

(6.1);

это

граф 3 -> 1 ,

1

2, 2 -> 4. Итак, с вероятностью, близкой при малых

е к 1, процесс, начинающийся вблизи

 

К2 или

К3,

достигнет окрестности А 4 раньше, чем

окрестности

Къ.

(Мы не утверждаем, что до этого процесс будет совершать переходы между компактами Кх, К2, К3только в наиболее

вероятном

порядке Кх

К2-> К3-> Кг\

более того,

в данном

случае можно

доказать, что до

достижения

окрестности А 4 будут с вероятностью,- близкой к 1, про­ исходить и переходы в противоположном направлении.)

Далее процесс будет ходить от Ki к Къ и обратно (это показывают графы 4 -> 5 и 5 -> 4); затем он вернется в окрестности компактов Кх, К2, К 3, причем вероятнее все­

го сперва к компакту Кх (как показывает граф 5 4,

§ 61

ЦИКЛЫ. СУБПРЕДЕЛЬНЫЕ

РАСПРЕДЕЛЕНИЯ

267

переходы между «циклами» Кг

К2-*■ К3

Кг и К х-у-

-*■

4, причем чаще всего

они будут

происходить

так, как описано выше. Таким образом, получается «цикл из циклов»— цикл второго ранга.

Переходы между окрестностями компактов К г и в об­ щем случае описываются с помощью иерархии циклов (Ф р е й д л и н [9], [10]). Пусть Ки . . ., К1о—- устой­ чивые компакты; Q — какое-то подмножество множества L — ( 1 , 2 , . . . , 10}. Предположим, что существует единст­ венный (£\<2)-граф g*, на котором достигается минимум (6.1). Определим /?Q(0> * ^ Q» как тот элемент из L\Q, в котором заканчивается цепочка стрелок в графе ведущая из i в множество L\Q.

Опишем теперь разбиение множества L тта иерархию циклов. Начнем с циклов первого ранга. Для каждого

г0е

L рассмотрим

последовательность

i0,

ii,

i2, . . .

. .

In, ••м

В которой in — 11{in_i) i’n-i)-

Пусть п

наименьший номер, при котором начнутся

повторения:

in =

о <

т <

гг, причем для /с,

меньших,

чем гг,

все гл различны. Тогда циклами первого ранга, порожден­

ными

элвхментом

i0 е L,

мы

будем называть

группы

{*о}»

{ * l } * * * •»

{ * m —l} »

{*m

^ * m + i

l **

-> zm}; причем последняя группа рассматривается вместе

суказанным циклическим порядком. Циклы, порожден­ ные различными начальными точками г0 е L, либо не пересекаются, либо совпадают, и в этом случае цикли­ ческий порядок на них — один и тот же. Итак, циклы первого ранга выделены (некоторые из них — только из одной точки).

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

циклы — 1)-го

ранга

(короче — (к — 1)-циклы)

я?” 1» я*-1 , •••,

 

уже определены;

это — множест­

ва, состоящие

из

— 2)-циклов,

причем снабжен­

ные определенным циклическим порядком. В конечном счете каждый цикл состоит из отдельных точек — элемен­ тов L; будем обозначать множество точек, из которых составляется цикл, так же, как и сам цикл. Будем гово­

рить,

что

цикл лj" '1 следует

за

и

писать

-*■

л^~1, ‘ если R ла

при

w е л*” 1. Проверя­

ется,

что

функция RЛ;ь_1 (т)

принимает

одно

и то же

268

ВОЗМУЩЕНИЯ НЛ

БОЛЬШИХ

ОТРЕЗКАХ

ВРЕМЕНИ

|ГЛ 6

значение при всех т е

п\ \ так

что

приведенное

опре­

деление корректно.

 

 

 

 

 

 

 

 

Рассмотрим теперь какой-то цикл я*^1

и

начинаю­

щуюся

из него

последовательность

циклов

я?0—1—>

я?” 1-> n tr1

 

я*” 1 "-*•••• В этой последователь­

ности

с

некоторого

номера

начнутся

повторения;

пусть

п — наименьший

такой

номер:

я^~1= я £ ~ 1>

О ^ т <

п. Тогда будем

говорить, что

цикл

я^"”1

по­

рождает

циклы

к-то ранга {я ?"1}, {я*-1 }, .

.

{ я ? ^ ]

циклов к-то ранга , состоящих каждый из

одного цикла

предыдущего ранга) и (я -" 1-*- я*“_^-* ...

 

я я - " 1}.

Беря

всевозможные начальные

— 1)-циклы я*” 1, раз­

биваем на /с-циклы все циклы — 1)-го

ранга.

 

Циклами нулевого ранга назовем отдельные точки множества L; при некотором к все — 1)-циклы войдут в один fe-цикл, который охватит все L.

Разбиение на циклы полностью определяет наиболее

вероятный при малых е порядок обхода траекторией А'? окрестностей устойчивых компактов (все это, разумеется, относится к случаю «общего положения», когда каждый минимум (6.1) достигается только на одном графе). Теперь обратим внимание на время, которое процесс будет про­ водить в пределах того или иного цикла.

Т е о р е м а 6.2. Пусть я — какой-то цикл. Положим

С (л) = А (л) — min

min

2

v (Km, K n), (6.2)

ien

g ^ G n { i )

(jn-»n)eg

где Л(я) определяется формулой (6.1),

a Gn{i} — совокуп­

ность всех {1}-графов над множеством я. Тогда при доста­ точно малых р >> О

lim е2 In М^Тд =

С (я)

 

(6.3)

с—►О

 

 

 

равномерно по х из р-окрестности U Кь \для любого у > О

 

ien

 

 

lim Р* {es-2(c;*)-v> < т„ <

е6

= 1

(6.4)

е-»0

 

 

 

равномерно по всем указанным х.

§ б]

ЦИКЛЫ. СУБПРЕДЕЛЬНЫЕ РАСПРЕДЕЛЕНИЯ

269

 

Д о к а з а т е л ь с т в о . Соотношение (6 3)

выводит­

ся из теоремы 5.3; напомним, что для произвольного мно­

жества Q С

L (не цикла) lim в2 In M*TQ зависит, вооб-

ще говоря,

е-*0

(см. § 5). Доказа-

от выбора точки х е (J

 

ieQ

 

тельство утверждения (6.4), аналогичного утверждению теоремы 4.2 гл. 4, опускаем.

Разбиение на циклы и теорема 6.2 позволяют ответить

на вопрос о поведении процесса X\ на отрезках времени длины £(е~2), где £(в”2) имеет показательную асимптотику

при в”2 -> оо: £ (е~~2) жеСг 2. Пусть л, л',

. . .,

л<*> —

циклы

предпоследнего

ранга,

объединяющиеся

непо­

средственно в последний цикл,

охватывающий все

точки

из L.

Если константа

С больше, чем С(л),

С(л'), . . .

.. ., С(л(5)), то за время порядка ехр {Се”2} процесс

успеет

уже много раз пройти все эти циклы (а внутри них — все циклы меньшего ранга), и для него уже успеет установить­ ся предельное распределение, найденное нами в § 4. В тер­ минах иерархии циклов это предельное распределение будет сосредоточено на том из циклов л, л', . . ., л(в), для которого соответствующая константа С(л), С(л'), . . ., или С(л(г5)) — наибольшая; в пределах этого цикла — на том из его подциклов, для которого соответствующая константа, определяемая по формуле (6.2), больше всего, и т. д. вплоть до отдельных точек i и соответствующих им устойчивых компактов К ь. То же будет, если не £(в”2) ж Ж ехр {Се-2 }, а только lim в2 1п£ (в-2 ) > max {С (л),

С ( л С (л<*>)}.

За времена £(в 2), имеющие меньший порядок, процесс не успеет пройти всех циклов, и общее для всех начала ных точек х предельное распределение еще не успеет уста­ новиться; при этом на циклах не максимального ранга, которые процесс уже успеет пройти, будут устанавливать­ ся некоторые субпредельные распределения, зависящие от начальной точки.

Вернемся к примеру, рассмотренному после теоремы 6.1. Изобразим иерархию циклов на чертеже (рис. 17):

циклы нулевого ранга

(отдельные

точки) объединяются

в циклы первого ранга

{1 2 -> 3

1} и {4 -> 5 -> 4),

а они — в единственный цикл второго ранга. На стрел­

270 ВОЗМУЩЕНИЯ НА БОЛЬШИХ ОТРЕЗКАХ ВРЕМЕНИ {ГЛ. б

ках, выходящих из каждого цикла я, указаны значения

констант С(я). Если lim е2 In t (е~2) между 0 и 2, за вре-

е->0

мя £(е~2) процесс еще не успеет отойти от того компакта К ь вблизи которого он начинался, и субпредельное распределение, соответствующее начальной точке х, бу­

 

 

 

дет сосредоточено на том К ь

 

 

 

к области притяжения

кото­

 

 

 

рого

принадлежит

х.

За

 

 

 

время

t{е”2),

для которого

 

 

 

lim е2 In t (е-2) между

2 и 3,

 

 

 

е-*0

 

 

схождение с

 

 

 

произойдет

 

 

 

0-цикла

{4)\ т. е. единствен­

 

 

 

ное,

что

произойдет,— это,

 

 

 

если

процесс

был

вблизи

 

 

 

компакта К 4, то он перейдет к

 

 

 

компакту Къ. При этом суб­

будет сосредоточено

на

пределыюе

распределение

КR и для начальных точек

в об-

ласти

притяжения

К 4,

в области притяжения

К*.

Если

lim е2 In t (е~2)

между 3 и 4,

за время t(e

2) уста-

 

е-»0

 

 

 

 

 

5

4},

новится предельное распределение на цикле {4

но не произойдет ничего больше; так как это распределе­ ние сосредоточено на точке 5, соответствующей компакту

Къ, результат будет

тот

же, что и для

lim е2 In t (е~2)

между 2 и 3.

 

 

 

е-*0

 

 

между 4 и 5 произойдет схож­

При

lim е2 In t (в-2 )

дение

с

его

{4

5

4}, процесс попадет в окрест­

цикла

ность

компакта

Кг и перейдет оттуда к

компакту К

но не дальше; субпредельное распределение будет сосре­ доточено на К2 для начальных точек, притягивающихся к любому компакту, кроме К3 (а для точек, притягиваю­

щихся к К3,— на К3). Наконец, при lim в2 In t (е~2) > 5 е-»0

устанавливается уже предельное распределение (хотя схождение с цикла { 1 - + 2 - + 3 - * ! } и даже с цикла {5}, может быть, еще не произошло).

Теперь мы можем указать субпредельные распределе­ ния для начальных точек из областей притяжения всех устойчивых компактов Kt. Например, для х%притягиваю-