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

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

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

15)

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

151

пости точки у0 минимума квазппотенциала: для

любого

6 >

О

 

 

lim P2{|Z1- y 0|<6|Z1e ^ } = l

(2.1)

 

е-*0

 

равномерно по z G у. Из этого уже легко вывести утверж­ дение теоремы.

Доказательство (2.1) будет, разумеется, использовать функционал действия и его свойства.

Сформулируем и докажем вспомогательные утвержде­

ния, нужные

для доказательства.

Л е м м а

2.1. Пусть F компакт в Rr, T u b

положительные числа. Тогда существую п такие положи­ тельные е0, р, что при любом х е F и г < е0

Р , {рог (X е, * (z ))> 6 )< e x p {— e~2p)t

где xt(x) траектория динамической системы, выходя­ щая из точки х.

Д о к а з а т е л ь с т в о . Положим

G{x) = {(р е CGT: ф0 = х, роГ(ф, х(х)) > б}.

В силу следствия из теоремы 1.1 нижняя грань d функцио­

нала S0T (ф) на замкнутом множестве U G(x) достигается x&F

на какой-то функции из этого

множества.

Функционал

S0T обращается в нуль только

на траекториях динамиче­

ской системы, поэтому

d > 0.

 

и [} Фл (df)

Для любого d' < d

множества |J G(x)

 

 

x^F

x^F

не пересекаются; обозначим расстояние между ними (по­ ложительное в силу замкнутости первого и компактности второго) через б'. Воспользуемся теоремой 1.1: для про­ извольного у > 0 при достаточно малых е для всех х е F

р* {рот (х\ х (*)) >

6} Рх=

(X е е G (*)}

<

 

 

< Р, {Рот (Х 0е,

Фх (d') > б' ) < ехр { - 8

2 (i' -

Y)| *

Итак, утверждение

леммы

выполнено

для

р = d'

— у

(в качестве р можно, таким образом, взять любое число, меньшее d).

Ш

 

ОКРЕСТНОСТЬ ПОЛОЖЕНИЯ РАВНОВЕСИЯ

|ГЛ. 4

 

В дальнейшем через $ 6(а) будем обозначать б-окрест-

иость точки

а е /?г.

 

 

 

 

 

Следующая лемма будет нужна нам также в § 4.

 

Л е м м а

2.2. Предположим, что точка О устой­

чивое положение равновесия системы (1.1),

область D при­

тягивается к О и (Ь(х), п(х)) <С 0 при х е

дВ. Тогда для

любого а >

0:

 

 

 

 

 

 

а) существуют такие положительные константы а и

Т0, что для любой функции ф*, принимающей при

[0, Т ]

значения

в множестве В \JdB\^a(0), имеет место

нера­

венство:

S0T(ф) :> а(Т — Г0);

 

 

 

кие,

б) существуют положительные константы с, Т0 та­

что

при всех

достаточно

малых г > 0 для

любого

х е

D\JdD\g>a(0) справедливо

неравенство

 

 

 

Р*{£« >

П

< ехр { -

г~*с{Т -

Ти)},

 

где

£а =

inf {<: X*

 

D\«ra (0)).

а) Пусть

&а'{0) — такая

 

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

окрестность

точки 0,

что траектории xt(x) динамической

системы, начинающиеся в '(О), никогда не покидают <^а(0). Обозначим Г(а, х) время, которое затрачивает тра­ ектория xt(x) на достижение <^v(0). Так как область D притягивается к О, то Т{а, х) < оо при х ^В \]дВ . Функ­ ция Т(а, х) полунепрерывна сверху по х (из-за непрерыв­ ной зависимости xt(x) от х). Значит, она достигает своего

наибольшего значения Т0 =

шах Т(а, х) < оо.

 

 

xeDUdD

значения в

Множество функций из Сото» принимающих

/) U дВ\&(Х(0), замкнуто

в Сот0Последствию из теоре­

мы 1.1 функционал SOTO

достигает на этом

множестве

нижней грани. Эта нижняя

грань отлична от

нуля, так

как в противном случае в указанное множество входила бы траектория динамической системы.

Итак,

для

всех

таких функций

^от^ф) ^

А > 0.

В силу аддитивности функционала S для функций ф, про­

водящих

B D [}d D \ ga(0) время Г, большее Т0, 5 0Т(ф) ^

для функций,

проводящих в D[JdD\g>a(0) время

Т ^ 2Т0, S0T(ф) ^ 2Л, и т. д. Вообще,

 

 

^ог(ф) >

А [Т,То] > А (т1то “ !)

= л(Г -

Г0).

б) Из того, что область D притягивается к О и что (£>(д‘), п(х)) < 0 на границе этой области, вытекает, что

§ 2]

ЗАДАЧА

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

153

теми же свойствами будет обладать и S-окрестпость D при

достаточно малом 8 >

0. Будем считать 8 меньшим,

чем

а/2. В силу утверждения а) можно выбрать такие постоян­ ные Г0 и Л, что5'0То(ф) >> А для функций, которые на от­ резке от 0 до 70 не выходят за пределы замкнутой 8-окрест­

ности D и не заходят в /2(0). При x ^ D

функции из

мно­

жества

ФХ{А) = {ср: ср0 = х, *!?ОТо(ф) ^

А) за время

от 0

до TQдостигают ^а/г{0) или выходят за пределы

8-окрест­

ности

области D\ траектории процесса

X*, для

которых

находятся на расстоянии не менее 6 от этого мно­ жества. Отсюда в силу теоремы 1.1 получаем, что при ма­ лых е для всех ^ G D

Р * { £ « > Г о } < е * р

{ —е~2(/1 -

Y)>.

Далее пользуемся марковским свойством:

Р *{£ «> (« + 1 ) ? ’»} = мл|£«> «Л»; р д.,г

а * > Г0}| <

< Р лК «> и ^ }* 8 и р Р Л 5 «> ?’о};

по индукции получаем

 

 

P H fe a > n < P 3e{se > [ | - ] r 0} ’

«ирРи{?а > TJ 1Г/Г.1 <

ye=D

 

< exp

е-2

— 1j (Л — Y)J;

Итак, в качестве с можно взять — у)/Г0, где у — про­ извольно малое положительное число.

Еще одна очень простая лемма.

Л е м м а 2.3. Имеется положительная константа L

такая, что для любых х

и у е

ПТ существует гладкая

функция ср*, ср0 — х, Фг =

у, Т = у |,

для которой

^от(ф) < L -\х — у |.

 

 

 

Действительно, можно положить Ф* = х +

~r-~'x t t.

 

 

 

I у — хI

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

Пусть 8 > 0. Обозначим

 

 

 

d = miu{F(0, у): у е дО, |у -

у01> 8} -

V(0, у0).

154

 

ОКРЕСТНОСТЬ ПОЛОЖЕНИЯ РАВНОВЕСИЯ

 

 

[ГЛ. 4

Раз

I/Q — единственная точка минимума

функции

 

F, то

d > 0 .

 

положительное

число ц <

d!bL

так, чтобы

 

Выберем

сфера Г этого радиуса с центром в О лежала

внутри

D

(L — константа из леммы 2.3).

 

малых

г

 

 

 

 

Л е м м а

2.4.

При достаточно

 

 

 

 

 

Рx{Zt <= £ £ >}> exp { —8“2(F(0,

y0) +

0,45d)}

 

 

для всех x e

у* (Напомним, что у сфера радиуса р/2

с

центром О.)

 

 

Выберем

точку

г/х

вне

 

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

D [jdD на расстоянии не более ц/2 от г/(). Существует Т >

0

такое, что для любой точки х е

у найдется функция ф?,

0 <

* <

Г,

ср5 = я, Фг = Уь

^0т(фх) <

 

 

Уо) +

0,4d.

 

 

Действительно,

прежде

всего выберем функцию ф*1*,

0 <

t <

Гх,

ф(оП=

О, фг? =

г/„

такую,

что

5 0Т1(ф(г>) <

 

V(0 ,y0)-|-0,ld. «Отрежем»

ее

начало до точки xt = ф^

последнего пересечения ф^

с Г, т. е. введем новую

функ­

цию ф(2) = ф(,+(, 0 < t < Г 2= 7 ,1—tt. Имеем ф!>2)=ж1,Фг1 =*/„,

‘5'ог2(ф(2)) = ‘5'(1г1(ф(1)) < ЩО, Уо) +

0,Id.

Далее,

по

лем­

ме 2.3 выбираем функции ф((3), 0

£ <1 Г3 =

р,

ф(03>

= О,

фт,' =

X lt

SOT,(<P{3))<0,2d; Ф(4), 0 < г < 7 ’4, Фо4>= !/о. Фг.)==

= г/j,

SoTt (ф(4>) ^ 0, Id. Наконец,

но той же лемме

для

любого х е у

выбираем зависящую от х функцию

фх5),

0 < t

<

Ть =

 

р/2, ф(05) =

х,

ф^ =

О,

S0Tl (Ф(5)) <

0,Id.

Функцию

ф?

составляем

из

кусков ф(5>, ф(3), ф(2> и

ф<4>:

ф? = Ф(,5) при

0 < t <

Т ъ;

= фi-T, при Г6<г < ^ 5

+ Г3;

=Ф(—т,—г, при

Ть+

 

Ть +

Г3+

Z,*;=V<—г,—г,—г,

при Г5 +

Г3+ Г 2< « < 7 ’8 +

Гз +

г,2 +

7’4

(рис. 5).

Выберем положительное б' меньше, чем ц/4 и чем рас­

стояние от точки ух до dD, и воспользуемся

теоремой 1.1.

Получим для е, меньших какого-то е0, и для всех х ^ у

Р*{Рог ( * 8, Ф*) < б'} > exp { - e - W i Уо) + 0,4d +

+ 0t05d)}.

Но если траектория X\ проходит на расстоянии менее б'

§ 2]

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

155

 

 

от кривой ср*, то она, не задавая у после достижения Г, достигает б'-окрестности точки ух и по пути пересекает границу dD. Значит, вероятность того, что Zx принадле­ жит dD, не меньше чем ехр { —е~2(У(0, г/0) + О^Ы)}.

Л е м м а 2.5. При достаточно малых е

PX{ZXе= dD\&6(г/о)} < ехр { —е -2(У(0 ,у0) +

0,55d)}

для всех х е у ,

 

Д о к а з а т е л ь с т в о . Вспомним, что Zx =

Х\х, где

Tx=inf {t> o Q: X] е у U dD}. Введем обозначение т(у [} dD) —

== inf {t > 0: Xеs у U dD). Случайная величина Zx — не

что иное, как величина X®(vuaD), вычисленная для отрез­ ка траектории после момента а0, сдвинутого на а0 по оси времени. Воспользуемся строго марковским свойством относительно марковского момента оу, получим

Р * iz i е d D \& 6(у0)} =

[X'wem ^ d D \& &(у0)).

Так как Х®0 е Г, то эта вероятность для всех х^ у не пре­

восходит supPa{X®(Yua£))<= dD\&&(yQ)}] оценим послед-

х б Г

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

т

ОКРЕСТНОСТЬ ПОЛОЖЕНИЯ РАВНОВЕСИЯ

ГГЛ. 4

 

В силу леммы 2.2 для любого с^> О найдется такое

7\

что PX{T(V U9£)) > Т) < ехр { —8~2с} для всех х е

Г

п

8, меньших некоторого е0. Возьмем в качестве с, скажем, У(Оу у0) + d; чтобы получить нужную оценку, остается

оценить Р х {т (у U dZ>)<7\ X jyvdD)zzdD\&6 (у0)}. Мы полу­ чим эту оценку с помощью теоремы 1.1.

Рассмотрим замыкание р/2-окрестности множества dD\&6 (у0); обозначим его К. Никакая функция ф*, Т, ф0 е Г, такая, что 5 0Т(ф) <; V(0, у0) + 0,G5d, не задевает множества К. Действительно, предположим,

что ф^ е К при каком-то tx <; Т; тогда Sotjy) — SoT(ф) ^

^

У(0,

у0) + 0,65d.

Берем по лемме 2.3 функции

ф*1^

О <

t <7-!. ф^0 = О, Фу> = ф 0, с 5оГ,(ф('))<0,2^ и ф<2),

о <

<

< <

Т* Ф(о2) = ф,„

фг2) е= d D \ gb(y0), с £ 0г. (ф<2>) <

0,Id

и составляем из кусков ф^, ф, ф<2>новую функцию: ф(7) =

=Ф((° при 0 <17 <1 7V, =

ф(-т,

при

7\ <

г <

7\ -f 7t;

= ф(-r .-i, при fj+ fi <

t <!

7\ +

гх +

Тг. При этом фо= 0,

^T i+ tj+ r^

dD\&t,(yo),

а

^от.+^+тДф) ^

0,2d

+ V(0,

Уи) + 0,Id +

0,65(7. Это

меньше

нижней

грани

V(0, у)

по у е= dDXg’biijn), чего не

может быть.

 

 

Это означает, что все функции из [}

Фх (F (О, уо)+0,65(7)

проходят на расстоянии не менее ц/2 от дО \$>ь(у0).

Поль­

зуемся теоремой 1.1;

получаем

для

достаточно

малых

е и всех х е

Г:

 

 

 

 

 

 

Р* {т (у U dZ))< Т,

Ar®(VU9D) s dD\g>6(г/0)] <

 

< Р* {рог (А 8, Фх (V (О, у0) +

0,65d)) > |х/2) <

 

 

<

ехр ( -

е“ 2 (V (О, Уа) + 0,65d - 0,05(7)),

Р* (A^vuaD)e= 9D\&o (г/о)! < Рх {т (у U dD) > Т } +

 

+ Р* {* (Y U <№) ^

Т’Дцтиэо) е

dD\gb («/„)]<

<

ехр I -

e-2

(F (О, у0) + (7)) +

 

+ехр { - е~2 {V (О, у0) -f. 0,6(7)) <

<е х р { - в~2(Г {0^л) + 0J55(7)}.

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

Из лемм 2.4 и 2.5 вытекает, что

PX{ZX<= d D \ $ b{ (/„)} < РX{ZX<= dD} exp { —e -2-0,ld}

для достаточно малых е и всех х е

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

наименьшее га,

при котором Zn е

dD. Пользуясь строго

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

Р* { |К* ~ У* |>

6} e Р* (2v s

d D \ g6 (у0)} =

оо

 

 

 

— 2

(v = ni

dD\<g6 (у0)} =»

П=1

 

 

 

оо

= 2 Мж\ZXе V,.. Zn_) <= V.

П=1

Ргп^ { г х^ д О \ ^ йЫ } ) <

 

 

 

 

оо

 

 

 

 

 

 

 

 

П—1

 

 

 

P zn -J Z x ^ 3Z)}]-exp {— г~2-0,Ы\ =

 

 

2 P*{v— п) ехР I-

е

2•0,1с?] = ехр {— е

2-0,lcf}

0

 

П=1

 

 

 

 

 

 

 

при

8

0.

Таким образом,

теорема доказана для х е

у.

Если

о: — произвольная

точка

области

то

 

Р, { IГ?-

».| > 6} < р , [X;(vuaD) е

dD] +

 

 

 

 

 

+

Р* {^?(vuaD)е V’ |

|^

 

Первая вероятность стремится к нулю по лемме 2.1; вторую преобразуем, пользуясь строго марковским свойством, к виду

М * {^ (v U a D ) ^ У ' ^ z t(VU0D) { I^'Се ~~ Уо\ ^

что стремится к нулю по уже доказанному.

На языке теории дифференциальных уравнений теоре­ му 2.1 можно сформулировать следующим эквивалентным образом.

158

ОКРЕСТНОСТЬ ПОЛОЖЕНИЯ РАВНОВЕСИЯ

[ГЛ. 4

Т е о р е м а 2.2. Пусть g(x) непрерывная функция% определенная на границе 3D области D. Рассмотрим в об­ ласти D задачу Дирихле

Т Аи*(*> + 2 ъ1(*>

И = °i * е

UB(x)=g(x),

X<=:dD.

Тогда, если выполнены условия теоремы 2.1, mo Пт ue(.r) =

е~>0

Д о к а з а т е л ь с т в о немедленно следует из фор­ мулы ие(х) = M3cg(X®e) (см. § 5 гл. 1), если принять во

внимание непрерывность и ограниченность функции g(x). Из теоремы 2.2 с помощью той же формулы, в свою очередь^ можно получить теорему 2.1.

При некоторых дополнительных предположениях можно получить более точную информацию о том, как

траектория X® при малых е выходит из области D. Здесь нам будет удобнее пользоваться обозначениями

X 8(Z), ср(t) вместо Xf, cpt и т.

п.

Мы определили

У(0, у)

как нижнюю грань SoT (ф)

по всем функциям

ф(£), 0 ^

t ^ Г, ведущим из точки О

ву. Эта нижняя грань обычно не достигается (см. примеры

вследующем параграфе), но она достигается для функций,; определенных на бесконечном влево полуинтервале: су­

ществует функция ф(/),— оо ^

7\ такая, чтоф(—оо)=*

= О, ф(Г) = г/,

о,т (ф) = У(0,

у).

Не будем

этого

доказывать, а введем в качестве условия в теорему,

кото­

рую мы собираемся

сформулировать.

(Утверждение со­

держится в статье В е н т ц е л я , Ф р е й д л и п а [4], лемма 3.3, снамеченнымдоказательством). Экстремаль ф(£)не един­ ственна: вместе с ней экстремалью будет любой ее сдвиг

ф(0 — ф(* + я), —оо ^ t ^ Т а.

Далее введем следующее определение. Пусть G — окрестность точки О с гладкой границей 3G. Кривая ф(Z), ведущая из точки О на границу 3D области D , непременно где-то пересечет 3G. Обозначим 0aG (ф) последний момент,

когда ф(£) находится на dG: QQG(ф)

= sup

ф* s 3G}.

Если при этом для некоторого а >

0 при t е

[0эс(ф) —

§ 21

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

159

вдв{ф)) функция ф(t)

принимает значения внутри G, бу­

дем говорить,

что ф(t) выходит из G неособым

образом.

Т е о р е м а 2.3. Пусть выполнены условия

теоремы

2.1, и пусть

существует единственная с точностью до

сдвига экстремаль ф(£), —оо <

t ^ Т, ведущая

из О на

dG {а именно, в точку у0, рис.

6). Продолжим ф{t) произ­

вольным непрерывным

образом

на

 

значения t ^

Т.

некоторая

 

Пусть G a D —

 

окрестность

положения равнове­

 

сия, с гладкой границей, причем эк­

 

стремаль ф(£) выходит из этой ок­ рестности неособым образом. Обо­

значим O^G последний момент пре­ бывания траектории Xe(t) на dG

до выхода на

dD: 0*JG = max {t <

 

< те: Xе (t) е

dG}• Тогда для любо­

 

го 8 >

О и х е D

 

П т Р Л

шах

|XB(t) —ф (t

 

е->0 )

.

е

 

 

0вс(ф))|< б|= 1.

Рис. 6.

Иными словами, последний участок марковской тра­ ектории до выхода на dD, но после выхода и з б е вероят­ ностью, стремящейся к 1 при е -> 0, расположен в малой окрестности сдвинутой соответствующим образом экстре­ мали ф(£).

Д о к а з а т е л ь с т в о этой теоремы проводится по тому же плану, что и доказательство теоремы 2.1. Устанавливаем, что для функций ф(£), 0 ^ t ^ Т, веду­ щих из О на dD и лежащих далее чем на расстоянии 6/2 от

всех сдвигов экстремали, Si)T(ф) ^

V(0, у0) +

d, где

d — положительная константа. Затем выбираем ^

d/5L Д

Дб/2, сферы у и Г. Функции ф(£),

Т, начинающиеся

на сфере Г, лежащие далее чем на 8/2 от сдвигов экстремали

ф(0 и такие, чго £ог(ф) <1 V(0,

у0) +0,65 d, не достигают

р/2 -окрестности множества dD.

Отсюда1 так же как при

160

ОКРЕСТНОСТЬ ПОЛОЖЕНИЯ РАВНОВЕСИЯ

 

[ГЛ. 4

доказательстве теоремы 2.1, выводим, что

 

 

Р, 1т U dD)

< Т, X е (т U

dD)) е 3D,

 

 

 

max

|X е (t) — ср {t

(ф)) 1^ 6) ^

 

 

e^G« < тЕ

 

 

 

 

 

 

 

<

ехр [ — е~2 (V (О, у0) +

0,6J)}

для

всех х е

Г,

 

 

 

 

Рх (X е (т U

3D) ее 3D,

max |Хв (£) —

 

 

-

Ф (f - 0*G + 00G (ф)) |>

6] < Рх (X8 (т (Y и

dD)) s

 

 

 

 

е dZ)} exp 1—

в

2 • 0,1с?}

для х е V, откуда получаем утверждение теоремы для этих

ж, а затем и для всех х е

D.

 

 

Если динамическая система (1.1) имеет в Z) более одного

устойчивого положения равновесия, то прежде, чем выйти

из области, траектория X? может совершать переходы от одного устойчивого положения равновесия к другим. В этом случае точка выхода из области уже, вообще говоря, зависит от начальной точки. Такое более сложное устрой­ ство динамической системы мы рассмотрим в гл. 6 . Здесь мы остановимся еще на том случае, когда в области D си­ стема (1.1) имеет единственный устойчивый предельный цикл.

Итак, пусть для наглядности D — область на плоско­ сти, гомеоморфная кольцу; граница 3D, как и раньше, предполагается гладкой. Предположим, что система (1.1) имеет в D единственный устойчивый предельный цикл 11, и область D притягивается к циклу Г1, т. е. траектория xt(x) при х е D с ростом t стремится к П, не покидая при

этом области D.

На границе 3D области!)

предполагает­

ся выполненным

прежнее условие

(Ь(х),

тг(д:)) < 0 .

Квазипотенциалом системы (1.1)

относительно цикла

II назовем функцию

 

 

V (П, х) = inf {.Sr,Г, (ф): Фг, е

II, <рг, = *}•

Следующая теорема доказывается так же, как тео­ рема 2.1.