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

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

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

§ 3]

СВОЙСТВА КВАЗИПОТЕНЦИАЛА. ПРИМЕРЫ.

161

Т е о р е м а 2.4. Предположим, что область D при­ тягивается к устойчивому предельному циклу П, и (Ъ{х)г п(х))< 0 при х е 9D. Предположим, что существует един­ ственная точка уЛе 3D, Зля которой

V{U, у0)= тщ У (П , у).

V&dD

Тогда для любого 6>0 при х е D

е “ Р “ ^ ^

- J/o I < б 1 = 1-

Утверждение теоремы 2.4 остается верным и в том слу­ чае, когда не вся область D притягивается к циклу П,) а только часть, внешняя по отношению к этому циклу, при условии, что внутренняя часть не имеет выхода на границу dD. При этом цикл может быть устойчив только снаружи.

Если экстремаль функционала S(cp), ведущая с П на dD, единственна с точностью до сдвига, то можно доказать аналог теоремы 2.3.

§3, Свойства квазииотенциала. Примеры

Впредыдущем параграфе было показано, как, зная функцию V(0, х) — квазипотенциал динамической си­ стемы, найти точку у0 на границе области, вблизи которой

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

тотики среднего времени М^т8, затрачиваемого траекто­ рией X8, чтобы достигнуть границы области. Поведение

при е

О инвариантной меры процесса X 8 также опи­

сывается

функцией V(0, х).

В этом параграфе мы изучим вопрос о вычислении ква­ зипотенциала V(0, х), установим некоторые свойства экстремалей функционала 5(ср) и рассмотрим примеры. Мы будем заниматься только квазипотенциалом динами­ ческой системы относительно устойчивого положения равновесия; случай квазипотенциала относительно устой­ чивого предельного цикла рассматривается аналогично.

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

162

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

[ГЛ. 4

Т е о р е м а 3.1. Предположим, что векторное

поле

Ъ(х)

допускает разложение

 

 

Ъ(х) = -S/U(x) + 1{х),

(3.1)

где функция U(x) непрерывно дифференцируема eD\jdD,

U(О) = О, U(x) >

0 и VU(x) Ф 0

при х ф О ,

причем

(l(x), VU{x)) = 0.

Тогда квазипотенциал V(0, я) динамиче­

ской системы (1.1) относительно

точки О совпадает

с 2U(x) во всех точках х

3ля которых

U(x) ^

^ U0 = min С/(i/). Z?c./m функция U(x) дважды непрерывно y&dD

дифференцируема, то единственная экстремаль функцио­

нала 5(ф) на множестве функций

фв,

—оо ^

s

Г,

ведущих из О в х, задается уравнением

 

 

 

 

 

ф, =

VU{(fs) +

/(ф8), s <= (—ОО,

Т),

 

фг =

X.

(3.2)

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

Если

функция

(ps

при

$ ее [7^,

Т2]

не выходит ii3D\jdD,

то

из

соотношения

 

 

г,

 

 

 

 

 

 

 

 

U (срг8) — U (фгх) = [

(V t% s),

Фв)^5

следует

неравенство

 

 

Тг

 

 

 

 

 

 

 

 

 

 

Т!

 

 

 

 

 

 

 

 

STiT, (Ф) = 4

f IФ. “

W (ф.) - I (ф.) | ds +

 

 

 

 

 

Тг

Тх

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

2 f (ф„ VU (ф,)) ds >

2 [U (Фг8) -

U (фГ,)].

 

(3.3)

 

Ti

 

 

 

 

 

 

 

 

 

Отсюда выводится, что STtTi(ф) ^ 2Щх)

для любой кри­

вой ф5, фГ1 =

О, Фг2—

гдео: таково, что U(x)^UQ.Действи­

тельно, если эта кривая за время от TLдо

Т2 не выходит

из D\JdD, то утверждение сразу следует из (3.3) (напом­ ним, что U(0) = 0). Если же ф5 выходит из D\JdD, то эта

функция в какой-то момент Т ^ (7\, Т2) проходит через поверхность уровня U(x) = U0. В силу неотрицательно­ сти и аддитивности функционала S(ф) получим

(ф) > ST%J(ф) > 2и (ф* ) = 2С/0 > 2и (X).

§ 3]

СВОЙСТВА КВАЗИПОТЕНЦИАЛА. ПРИМЕРЫ

163

 

С другой стороны,

если q)s — решение задачи (3.2),

то

Ф_оо = О.

Это

вытекает из

того, что dU{y^!ds

=

=

|VС/(Фз) | >

0 при

Ф О, и точка О — единственный

нуль функции

Щх). Отсюда

получаем:

 

 

 

 

г

 

 

 

5 _ 00,г ( ф) =

2

J (ф „ V £/($s))^=2[C /(x)-C /(0)]= 2C /(x).

 

 

—оо

 

 

 

Таким образом, для любой кривой ф5, соединяющей точки

О и

х,

имеем:

£(ф) ^

2U(x) = £(ф);

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

7(0,

х) =

inf £(ф) = 2U(x), и функция ф8

— экстремаль.

Если функция

U(x)

дважды непрерывно дифференци­

руема, то решение уравнения (3.2) единственно, а значит, и экстремаль, ведущая из О в #, нормированная условием Фт = х, только одна.

Из доказанной теоремы следует, что если поле Ь(х) по­

тенциально, т. е. b(x) = —VU(x), то функция 7(0, х) лишь множителем отличается от потенциала U(x). Поэтому мы и называем функцию 7 квазипотенциалом.

 

Если векторное поле Ь(х) допускает разложение (3.1),

то

из условия ортогональности VU(x) и 1(х) = Ь(х) +

+

VU(x) для квазипотенциала получаем уравнение

 

-Y (VF (О, х)х W (V, х)) + (Ъ(х), VV (V, х)) = 0, (3-4)

которое, по существу, есть уравнение Якоби для вариаци­ онной задачи, определяющей квазипотенциал. Теорема 3.1 говорит, что решение этого уравнения, удовлетворя­

ющее условиям 7(0, О) = 0, 7(0, х) > 0 и V7(0, х) Ф О при х ф 0, есть квазипотенциал. Можно доказать, что на­ оборот, если квазипотенциал 7(0, х) непрерывно диффе­ ренцируем, то для него выполняется уравнение (3.4), т. е. имеет место разложение (3.1). Однако легко привести пример, показывающий, что функция 7(0, х) может не быть дифференцируемой.

Совсем плохой функция 7(0, х) все же не может быть: она удовлетворяет условию Липшица; это легко выводит­ ся из леммы 2.3.

Теперь обратимся к исследованию экстремалей функ­ ционала S(ф), ведущих из О в х. Для них обычным образом

6*

164

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

[ГЛ. 4

выписываются уравнения Эйлера:

Если векторное поле Ъ(х) имеет разложение (3.1), то для экстремалей выписываются более простые урав­ нения (3.2). Из них можно сделать ряд качественных выводов о поведении экстремалей. В каждой точке области скорость движения по экстремали в силу (3.2) равна

VU(x) + 1(х), тогда как b(x) = —SjU(x) + 1(х). Отсюда заключаем, что скорость движения вдоль экстремали в точ­ ке равна по величине скорости движения траекторий сис­

темы: Ivf7(s) + 1(х)\ = I—VU(x) + 1(х)| = \Ь(х)\. Пос­ леднее утверждение остается справедливым, и если не предполагать существования разложения поля Ь(х). Это вытекает из следующей леммы.

Л е м м а 3.1. Пусть £(ср) < оо. Обозначим ср функ­ цию, получающуюся из ф путем изменения параметриза­

ции таким образом, что |ф5|= |b(cps)| при почти всех s.

Тогда £(ф)

5(ф),

причем равенство достигается только

тогда,

когда |ф,| =

|Ь(ф8)| при почти всех 8.

 

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

При

замене

$ =* $(£),

ф* =

*= фз(^)

получим

 

 

 

 

 

<5(ф) =

f |Ь(фа)—фз I2^

 

t(T2)

 

 

 

j 1&(ф< ) -

 

 

т,

 

ЦТ,)

 

«г,)

 

 

 

 

 

 

 

 

 

- <pt« (о-1 12«(оd t = 4 ~

Iъ(фОI2« (o+fol2®(о-1) ^ -

 

t(T,)

_

^

«т,)

^

^

 

 

- J

(Ь(ф<). 4>t)d t >

J

|б(ф|)Иф»|Л-

 

 

 

 

 

«Г.)

 

 

 

 

 

 

 

ЦТ,)

^

 

 

 

 

 

~

S ф^ Ф * )# *

(3-5)

ЦТ, )

§ 3]

СВОЙСТВА ЦВАЗИПОТЕНЦЙАЛА. ПРИМЕРЫ

165

где t(s) — функция, обратная к s(t). }В (3.5) мы воспользо­ вались неравенством ах2 + a “ V 2 ^ 2ху, справедливым при любом положительном а. Определим функцию s(t) из равенства

 

8 (0

 

* =

J |ф«1|ь(Фи)Г^«-

 

 

О

 

Эта функция монотонно возрастает, а ср* =

абсолют­

но непрерывна по t.

Для этой функции |ср*| =

|5(ф*)| при

почти всех t, и для нее неравенство (3.5) обращается в ра­ венство

ЦТ,) ^ ^ ЦТ,)

s ( ф) =

J |ь (ф<)М ф* 1 ^ -

J

 

t(Ti)

 

ЦТ,)

 

Это доказывает лемму.

некоторые

примеры.

 

Рассмотрим

теперь

систему

П р и м е р

3.1. Рассмотрим

динамическую

xt = b(xt) на

прямой,

Ь(0) = 0,

Ь(х) > 0 при х < 0 и

Ь(х) < 0 при х > 0; D — интервал (ах, а2) С

содер-

 

 

 

X

 

жащий точку 0. Если

обозначить

U (х) = — J b (у) dy, то

 

 

 

о

 

b(x)~—dU/dx, так что в одномерном случае любое поле по­

тенциально. Из теорем 2.1 и 3.1 следует, что первый вы- t

ход процесса Xf = х + J b(Xs)ds +ew t из интервала (alf a2)

о

при малых е с вероятностью, близкой к единице, происхо­ дит через тот конец а*, в котором функция U(x) принимает

меньшее значение.

Пусть

для

определенности U(aх) <

<

U(a2).

Уравнение для

экстремали

имеет

вид

ф5

=

=

—Ь(ф8);

следует

брать

то решение

ф8, —оо

<

s ^

Г,

этого уравнения, которое при s

—с» стремится к нулю,

ифт = <хх.

Водномерном случае функцию vz(x)=Px {Хте = аг} можно, решая соответствующую краевую задачу, вычис­

166

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

IfЛ. 4

лить

явно:

 

 

V* (*) = j e"*~h w dy •( ] e - ^ ^ d y

 

х\а х

Из

этой

формулы

нетрудно,

койечно,

усмотреть,

что

lim Рх [Xl* = ах} =

1,

если

U(a1) <

U(а2).

 

е-*0

 

 

 

 

 

 

 

 

 

 

 

 

 

Интересно, используя явный вид функции 1>е (я), ис­

следовать, как выходит процесс Хг

из

интервала

(а1э

а2), если U(a1) =

С/(а2). Теорема 2.1 в этом случае непри­

менима.

 

 

 

 

 

 

 

 

асимптотического

оце­

Используя метод Лапласа для

нивания интегралов (см., например,

Е в г р а ф о в

[1]),

получим,

что

для

х е

(ах,

а2) при е

О

 

 

 

 

а2

 

 

 

ea

exp{2e

2U(аг)}

 

 

 

 

]'

ezs~ 2u(v)dy

 

 

 

 

2

 

| V (a2) |

*

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

e z

2m)dy

e2

exp{2e“ V (a ,)}

 

 

 

 

2

 

|U' («!)|

 

 

 

 

 

CCi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда

находим,

что

если

Ща^ — U(a2) и £/'(а() Ф О,

i =

i,

2,

то

при

в-> О

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

—1

 

 

 

 

P *U xee =

а 4}-*-

 

(ai)l

 

 

 

 

 

 

 

 

|У'(а»)1

 

 

 

 

 

 

 

 

и'(о1) Г ‘ +

 

Таким образом, в этом случае траектории X f в пределе

при

е -*• 0 с

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

вероятностью выходят из

отрезка через оба конца, причем вероятность выхода че­ рез конец ос* обратно пропорциональна |С/'(осj) |= |Ь(а*)|.

П р и м е р 3.2. Рассмотрим линейную однородную

систему xt = Axt с постоянными коэффициентами в Rr. Предположим, что матрица А нормальна, т. е. А А* = = А*А, и что симметричная матрица А + А* отрица­ тельно определена. В этом случае начало координат О является асимптотически устойчивым положением равно­ весия. Действительно, решение нашей системы можно

§ 31 СВОЙСТВА КВАЗИПОТЕНЦИАЛА. ПРИМЕРЫ 167

представить в виде xt =* eAix }, где х0 — начальное усло­ вие. Из нормальности матрицы А вытекает нормальность матрицы ем . Используя это заАчечание и соотношение

eMeA 4 = e U+A*)l'

получЩ1

 

 

 

| = (И Ч

eAtx0) =

(е(А+А*»х,О’

+>)<!

е,А+А*»х0|.

Так как матрица

А +

Л*

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

определена, то

16(A+A*)fT01

Q, и значит,

|.г*|2-> 0 при

любом

началь­

ном условии

х0.

 

дифференцированием легко

прове­

Непосредственным

рить, что векторное поле Ах допускает разложение

А

 

ril

А А- А*

\ . А — А*

/0

Ах = — V I--------j— х, х \-|------^— х.

(3.6)

При этом векторные поля V ^-----~(Л + Л*) х, х^ =

-----~-х

X (4 + 4 * ) х H-ij- (4 —4*)х

ортогональны: I-----— (4 + 4 * ) х,

-i-(4 — 4 * )х ) =

— -^ [(4 х , 4 х)—(4*х,

4*х)] =

— J - x

X [(А*Ах, х) (АА*х, .г)] = 0.

Пусть правые части нашей системы возмущаются про­ цессом белого шума:

А* = АХ\

и мы интересуемся тем, как траектории процесса Х\ выходят из ограниченной области D, содержащей положе­ ние равновесия О. Из теоремы 3.1 и формулы (3.6) заклю­ чаем, что квазипотенциал V(0, х) нашей динамической системы относительно положения равновесия О равен —V2((^ + А*)х, х). Чтобы указать на границе dD об­ ласти D точку, вблизи которой с вероятностью, стремя­

щейся к 1 при б — 0, траектории X? впервые выходят из D, нужно найти минимум V(0, х) на 6D. Уравнение для экстремалей имеет вид

Ф, = - ^ - ( 4 + 4*) + ± (4 - 4*) Ф, = — 4*Ф{.

168

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

[ГЛ. 4

Если у0 ^

6D — единственная точка, где У(0, х) достига­

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

по времени записать в виде ф4

* АНУо

 

Пусть,

например, на плоскос­

ти R2 задана

динамическая си­

стема

 

Л

 

 

x t =

Xt,

 

Xf

 

X? = х\

Xt.

Матрица этой системы нормаль­ на, и начало координат асимпто­ тически устойчиво. Траектории системы представляют собой лога­ рифмические спирали, наворачи­ вающиеся на начало координат но часовой стрелке (рис. 7). Квази­ потенциал равен [(я1)2 + (я2)2],

так что

выход

процесса Xf из

области D происходит с подавляю­

щей при

малых

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

вблизи точки y0^.dD, ближайшей к началу координат. Легко проверить, что экстремали то­ же представляют собой логарифмические спирали, но сво­ рачивающиеся с точки О по часовой стрелке.

§4. Асимптотика среднего времени выхода

иинвариантной меры для окрестности положения равновесия

Пусть, как и в § 2, D — ограниченная область

в Rr,

D э О — устойчивое положение равновесия системы (1.1),

т8 — первый момент выхода процесса X? из

обла­

сти D. В этом параграфе, в предположении, что область D притягивается к О, вычисляется главный член In М^т8 при 8 -> 0, а также главный член In тг (D), где т* (•) —

инвариантная мера процесса X?, D = Rr\D. Так как существование инвариантной меры и ее свойства зависят

§ 41

ВРЕМЯ ВЫХОДА И ИНВАРИАНТНАЯ МЕРА

1б§

не только от поведения векторного поля Ъ(х) при х е

D,

то при изучении те (D) нам придется сделать некоторые

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

относительно

поля во

всем

простран­

стве

Rr.

4.1. Пусть О асимптотически устой­

Т е о р е м а

чивое

положение

равновесия

системы

(1.1),

и область

D CZ Rr притягивается к О. Предположим, что граница dD области D есть гладкое многообразие, и пусть при х е= dD скалярное произведение(Ь(х), п(х)) < 0, где п(х) внешняя нормаль к границе области D. Тогда при х е D

lim е2 lnM^x8 =

F0 = min F (О, у).

(4.1)

е-»0

у е а о

 

Здесь функция F(0, у)

— квазипотенциал

динамиче­

ской системы (1.1) относительно точки О.

Д о к а з а т е л ь с т в о этой теоремы использует рас­ суждения, конструкции и обозначения, которые приме­ нялись при доказательстве теоремы 2.1.

Чтобы доказать (4.1), достаточно убедиться в том, что

для любого

d > 0 найдется

е0

такое, что при е <; е0 вы­

полняются

неравенства:

d,

б)

е2 In М^т8 > F0 — d.

а) е2 In М^т8 < F0 +

Докажем сначала неравенство а). Выберем положи­

тельные числа ц, h, Тъ

Т2 так, чтобы выполнялись сле­

дующие условия: во-первых, все траектории невозмущен­ ной системы, начинающиеся в точках х ^ D \}dD, до мо­ мента Тг попадают в ц/2-окрестность положения рав­ новесия, и никогда после момента Т1 этой окрестности не

покидают;

во-вторых,

для каждой точки х,

лежащей

в шаре G =

е Rr:

0\ ^ ц} существует функция

Ф* такая, что ф* = х,

Ф? достигает внешности

A-окрест­

ности области D в момент Т(х) ^ Г2; при этом G? не по­ падает в р,/2-окрестность О после выхода из G и

S QT(X)(4>X) < V Q + - y *

Первое из этих условий может быть удовлетворено, так как О — асимптотически устойчивое положение рав­ новесия, область D притягивается к О, а для х е dD вы­

полняется

неравенство (Ь(я),

п(х))

0.

Функции

ф*,

входящие

во

второе условие, были

построены

в

§

2.

 

 

170

 

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

[ГЛ. 4

Для

i / e G

из определения функционала

действия

получаем

 

 

 

Ру\

sup

\Xt — Ф ?| < М >

 

to<t<T(y)

J

 

>

exp j — e-2

^SoT(y) (q>»0 + -§-J} > exp ( - e-2 (У0 + d)} t

если только 8 достаточно мало.

Так как точка у\у) не принадлежит /г-окрестности области D, то из последнего неравенства заключаем, что

при достаточно малых e n i / e G

 

Ру 1те < Тг\>

Ру W8 < Т (1/)1 > ехр { - е-2 (У0 +

d)}. (4.2)

Обозначим

а

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

шара G:

а0 = min [ t : X] е

С). Используя свойство строгой мар­

ковости процесса

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

при доста­

точно малых е будем иметь

 

 

 

Р* {*Е<

+ Т г) > М х { а < Г1;

(т* <

Г,}} >

 

 

> Р х{о < Fi}-exp {— е~2 (Уо +

<?)}>

 

 

> -| -ех р

[— е_ 2 (У0 -|-d)].

(4.3)

Здесь мы воспользовались неравенством (4.2) для оценки Р е (тЕ< Т2} и тем, что при е -> 0 траектории Х\ равно-

мерно на отрезке [0, Тг] сходятся по вероятности к xt.

Далее, используя марковское свойство процесса Хьи из (4.3) получаем

М*те< S (п + 1)(7\ + Тй) Рх {п (Т! + Т2) < ' п-0

< т е < ( п + 1)(7’1+ 7 ’2)} =

=

+

Тг) £

Рх

> п (Тг + Т2)} <

 

 

п=0

 

 

<

(Тг +

Тг) 2

1 -

min Рг (те < Тх+ Тг)

n=o L