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

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

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

ГАУССОВСКИЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ И ПОЛЯ

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 .

SQT (ф) на В Д

§ 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;

 

 

 

 

(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, этот переход произойдет в точку, лежащую в окрест-