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

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

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

§ 4]

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

171

= 1 +

TYKmin Р2 {те < Тх + Т2})~{ <

 

 

гео

 

 

< 2 ( 7 \ + r 2)e x p ie -2(F0 +

d)l ,

если только е достаточно мало. Отсюда вытекает утверж­ дение а).

Докажем теперь утверждение б). Введем в рассмотре­ ние марковские моменты т&, ok и марковскую цепь Zn, определенные при доказательстве теоремы 2.1. Фазовым

пространством

цепи

Ъп

служит множество у [ j

3D,

где

у =

{х ^

Rr: \х — 0\ =|л/2}. Для вероятностей

перехода

этой

цепи

за один

шаг

имеется оценка

 

 

Р , 3D) ^

шах Ру {тх =

т8} =

шах [Ру{т8 = тг <

Т) +

 

 

 

уеГ

 

 

 

уеГ

 

 

 

 

 

 

 

 

+ Pj,{xe = t i > 7’}].

(4.4)

Как следует из леммы 2,2,

Т можно выбрать столь

большим, что вторая вероятность допускает оценку

 

Ру {*8 =

>

Т) < 4 - exp I - е "2 (V0 - h)}‘

(4-5)

Чтобы оценить первую вероятность, заметим, что траекто­

рии

Xf,

0 ^ t ^ Г, для которых т8 =

< Г, при лю­

бом

h >

0 находятся на

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

расстоянии от

множества функций {ф е

С^Т(ЯГ): ф0 = у,

50Г(ф) < F0—

— /г/2}, если только ц достаточно мало. Отсюда на основа­ нии свойства функционала действия получаем, что

Ру{тЕ = *1 < Т) < ехр {—e-2(F0 — Щ

при у е Г и достаточно

малых е, > 0.

и (4.5) сле­

Из последнего неравенства, оценок (4.4)

дует:

 

 

Р(я, 3D) ^

ехр { —e -2(F0 — fe)},

(4.6)

если е, ц достаточно малы.

Обозначим, как и в теореме 2.1, v наименьшее п, при

котором Zn =

Х\пе

3D. Из (4.6) следует, что

P*(v >

П} >

и - ехр { —е -2(У0 - fe)}

172 ОКРЕСТНОСТЬ ПОЛОЖЕНИЯ РАВНОВЕСИЯ [ГЛ. 4

при х е у.

Очевидно, т8 = (х1 — т0) +

(т2 — тх) + . . .

 

00

 

 

. . .+(TV—тv-j), поэтому Мд-Т8 = 2 Мд. ( v >

п; т„ — хп-

{).

Пользуясь

п— 1

 

Xf

строго марковским свойством процесса

относительно момента on- v получим, что

 

 

Мх {V > щ хп— тп_ А} > Мх {v > п; хп— ап_ !} >

^ Рх { v ^ n}* inf Млтх. асег

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

Собирая все полученные оценки вместе, при достаточно малых 8 и [х для х ^ у получим

Мхте >

* i 2 ( l -

n *ev

n

- exp { - e " 2 (F0 -

h)})n~ l = t, exp {e~2 (F0 - A)] .

Отсюда вытекает утверждение б) при ж е у . Для лю­ бого учитывая, что Ра{т8 > хх) -> 1 при 8 -> О, будем иметь

МдТ8 = Мд. (т8 < Tj; т8} + М* {те > тх; т8} >

> М Я{х8> Xi, M8t т8} >

exp {(F0 — h) e- 2 ]x

X

Р Д т8>

t ! } > -j-exp {e_ 2 (F0 — A)}.

Тем самым доказано утверждение б) для

любых х е D.

З а м е ч а н и е .

Проанализировав

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

этой теоремы, можно заметить, что предположения о глад­

кости многообразия

3D и

о том, что (Ь(х), ?г(.г)) <

0 при

х е 3D можно ослабить.

Вместо них достаточно

пред­

положить, что границы D и замыкания D совпадают и что

при любом х е 3D

траектория динамической

системы

xt(x) для всех t > О расположена в D.

распреде­

Приведем еще один результат, касающийся

ления случайной величины т8 — момента первого выхода из D.

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

173

Т е о р е м а 4.2. Предположим, что выполнены усло­ вия теоремы 4.1. Тогда для всякого а 0 при х е U

lim Рж(ее~ 2(У‘ - а) < т8 < ее~2(^ +<х)} = 1.

 

 

е->0

 

 

 

 

 

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

Если UmPxltE> 6e

2(у«+а)}>

> 0

для

какого-то

а > 0, то lim

8->0

 

,

е2 In МЛте ^

V0 ccf

что

противоречит

теореме

4.1.

Поэтому

при люоом

а >

0 и

х е В

 

 

 

 

 

 

 

lim Рх {т8 < exp (е~2 (F0 + а )]} =

1.

$.7)

 

 

8—»О

 

 

 

 

 

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

Рх |тЕ< eB_2(V“-a )} <

М* {ti < Т8,

S

РХ\ W -

»,

 

 

 

 

 

v

 

П=1

 

 

 

 

те < exp [е~2 (F0 — a ))l} +

Р*{1е =

Ti}-

(4-8)

 

Последняя вероятность в правой части (4.8) стремится

к нулю.

Оценим

остальные

слагаемые. Пусть

те =

=

exp {e -2(F0 — а )}]; постоянную С мы выберем поз­

же.

Тогда

при

 

 

 

 

 

 

 

2

Ржlv = «, тЕ< exp {е-2 (F0 — a)} <

 

 

n=l

 

 

ОО

 

 

 

 

 

 

< P * {v <

тг) +

p x {v =га, т„ < ex p {e _2(F0 — a )ll<

2

 

 

n=m e

 

 

 

 

 

 

< P x{v < m e} +

Px {Tme< e x p {e -2(F0 - a ) } / .

(4.9)

Используя неравенство Px{v =

1} <; exp { —e -2(F0 —h)},

которое выполняется при x e

у,

h >

0 и достаточно ма­

лых е, получим

 

 

 

 

 

 

 

P«{v < mt} ^ 1 — (l — exp {— e

2 (F0 — h)})me-+ 0

(4.10)

174

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

1ГЛ. 4

вместе с е при любых С, а > 0 , если только h достаточно мало.

Оценим второе слагаемое в правой части (4.9). Можно выбрать такое 0 > 0, что Рл.{тг > 0} 1/2 при всех х е= у, 8 > 0. Для числа успехов Smв т испытаниях Бер­ нулли с вероятностью успеха 1/2 выполняется неравен­ ство

 

P{Sm >771/3} > 1 - 6

 

при т > /п0.

Тогда,

так

как тт =

(тх — т0) + (т2 —

— тх) + . . . +

(тт — Тщ-j),

используя

свойство строгой

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

получим

 

 

р* К е < ee~2(V“” а)} = Р*

< 4-] < б>

<4Л1)

0

 

1

и тг достаточно велико.

 

если -g- >

 

 

Собирая вместе оценки (4.8) — (4.11), приходим к со­

отношению

 

 

 

 

 

lim P jx 6 < ее_2(Уо_а)) =

0, x<=D.

(4.12)

 

 

Е—>О

 

 

 

Из (4.7) и (4.12) вытекает утверждение теоремы.

е 0

Перейдем

теперь к изучению

поведения при

инвариантной меры процесса Xf, определяемого уравнением (1.2). Для существования конечной инвариантной меры процесса Xf необходимо сделать некоторые предположения о поведении поля Ь(х) в окрестности бесконечности. Если

никаких предположений не делать, то траектории Xf могут, например, с вероятностью 1 уходить на бесконеч­ ность, в этом случае конечной инвариантной меры не су­ ществует. Мы будем предполагать, что вне достаточно большого шара с центром в начале координат проекция поля Ь(х) на радиус-вектор г(я), ведущий в точку х, отри­ цательна и отделена от нуля, т. е. существует столь боль­ шое число Nr что (&(#), г(я)) < — 1IN при \х\ > N.

Это условие, которое в дальнейшем называется услови­ ем Л, обеспечивает достаточно быстрое возвращение про­

цесса Xf в окрестность начала координат и существова­ ние инвариантной меры. Доказательство можно найти

§ 4)

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

175

у Х а с ь м и н с к о г о [1]; там же содержатся и более общие условия, обеспечивающие существование конечной инвариантной меры.

Если инвариантная мера рЕ(*) процесса Х\ суще­ ствует, то она абсолютно непрерывна по лебеговой мере,

и плотность те (х) = удовлетворяет стационарному

прямому уравнению Колмогорова. В нашем случае это уравнение имеет вид

Т

X АтЕ (*) - 2 ^1

(* )mS(*)) = °-

(4ЛЗ)

Вместе с дополнительными

условиями J rne(x)

dx —1,

 

R r

 

гпг(#)> О, это уравнение единственным образом определяет функцию т? (х).

Рассмотрим сначала случай потенциального поля Ь(х): b(x) = —VU(x). Условия существования конечной ин­ вариантной меры в этом случае означают, что потенциал U(x) достаточно быстро растет вместе с \х\, например, быстрее некоторой линейной функции а\х\ + р. Оказы­ вается, что если поле Ь{х) потенциально, то плотность ин­ вариантной меры может быть найдена в явном виде. Не­ посредственной подстановкой в уравнение (4.13) проверя­ ется, что

тг (х) = сг exp { —2s-2 £/(#)},

(4.14)

где с6 — нормирующий

множитель, определяемый из

условия нормировки се=

^ J^exp {— 2е 2U (#)} =

dxj- 1.

Сходимость интеграла, входящего в определение се, является необходимым и достаточным условием существования конечной инвариантной меры в потенциальном случае.

Пусть D — некоторая

область в Rr.

Тогда р,8 (D)

= с8 J ехр { —2e-2U(x)}dx.

Используя

это представле-

D

 

 

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

0.

Пусть U(х)

0, и в некоторой точке О потенциал

обра­

176 ОКРЕСТНОСТЬ ПОЛОЖЕНИЯ РАВНОВЕСИЯ (ГЛ. 4

щается в нуль: U(0) =

0. Тогда нетрудно проверить, что

е2 In р8 (D) = е2 In с£+

е2 In \ехр (— 2е“ 2£/ (я)) dx-+

 

Ъ

 

inf 2U(x) (4.15)

 

x& D

при 8 —> 0. С помощью метода Лапласа можно найти и бо­

лее

точную асимптотику р8 («) при

е -> 0

(см. Б е р н ­

 

ш т е й н

[1],

Н е в е л ь -

 

с о н

[1 ]).

Если

поле

 

Ь(х) не

потенциально,

то,

 

вообще говоря, нельзя вы­

 

писать

 

явное

выражение

 

для плотности инвариант­

 

ной меры. Тем не менее

 

оказывается,

что

соотно­

 

шение

(4.15) сохраняется,

 

если под 2 U(x) понимать

 

квазипотенциал поля Ъ(х).

 

Т е о р е м а

4.3. Пусть

 

точка

 

О единственное

 

устойчивое положение рав­

 

новесия

системы (1.1), и

 

все пространство

Rr при­

 

тягивается к

точке

О.

 

Предположим,

что

выпол­

нено

условие А. Тогда процесс Х\ при каждом 8 >

0

име­

ет единственную инвариантную меру

р8,

и

для

любой

области D CL Rr с компактной границей dD,

общей для

D и замыкания Z),

 

 

 

 

 

lim 82 in р8 (D) = -

inf V (О, х),

(4.16)

 

е-»0

 

 

D

 

 

где V(0,

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

поля

Ъ(х) относительно

точки О:

V(0,

х) =

inf {£ от(ср): Ф е

C0T{Rr),

ф0 = О*

Фг х, Т > 0}

(рис.

8).

 

 

 

Наметим доказательство этой теоремы. Как уже отме­ чалось, из условия А следует существование и единствен­ ности конечной инвариантной меры. Чтобы доказать (4.16), достаточно проверить, что для любого h > 0 найдется е0 = eQ(h) такое, что при е < е0 выполняются

« 41

ВРЕМЙ

ВЫХОДА

Й ИЙВАРНАЙТЙЛЙ МЕРА

1?7

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

 

 

 

 

а)

>

ехр{—е~2(У0 + Л)},

 

 

б)

|.i* (D) <

ехр{—e -2(F0 — /г)},

 

где F0 = inf V(0, z).

Если У0 = 0, то неравенства а) и б) очевидны. Мы оста­ новимся на случае VQ> 0. Ясно, что VQ> 0, только если р(0, D) = р0 > 0. Для доказательства неравенств а) и б) воспользуемся следующим представлением инвариант­ ной меры. Пусть у и Г, как и раньше,— сферы радиуса р/2 и р соответственно с центром в положении равновесия О; ji < р0. Как и при доказательстве теоремы 2.1, рассмот­ рим возрастающую последовательность марковских мо­ ментов т0, а0, т1? а1? т2, . . . Из условия А вытекает, что все эти моменты конечны с вероятностью 1.

Последовательность X®t, Х®8, . . . , Х®п, ... образует

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

нормированную меру Iе (dx). Как

вытекает,

например,

из

работы

Х а с ь м и н с к о г о

[1],

нормированная

инвариантная мера ре (•) процесса

X]

следующим обра­

зом выражается через инвариантную меру цепи

(Х®п)

на

у:

 

 

 

 

 

 

 

 

 

 

^ (D) = Се J М* ( %D{XI) ds I« (dx)t

 

(4.17)

 

 

 

V

о

 

 

 

 

 

где

%D{x) — индикатор

множества

D,

а

множитель с8

определяется из условия нормировки р8 (Rr) =

1.

Поло­

жим

TD =

min {£: Xf

е D у dD}.

Из

(4.17)

имеем

це (D) =

Сг j

Мх j xD (х*) ds Iе(dx) <

 

 

 

 

 

 

 

v

о

 

 

 

 

 

 

 

 

 

^ с ешах PX{^D <

ti}* шах

i.

(4.18)

y€.dD

178 окрестность положения Равновесия [гл. 4

Из условия А и компактности 8D следует, что равномерно пэ е < е0 < I

шах Мутх < аг < оо.

(4.19)

V&D

 

Далее, из анализа доказательства теоремы 2.1 вытекает, что при достаточно малых ц и е и ж е ^

РХ{Ъ < T*i} <

exp { - e - 2(F0 - hi2)}.

(4.20)

Принимая во внимание

соотношение

 

(dx))-\

v

нетрудно проверить, что при достаточно малых р, и е

0 < 1 п с е < А

(4.21)

Из оценок (4.18)—(4.21) заключаем, что при достаточно малых 8 выполняется соотношение

In pe (D) <

^

 

V o-h

 

 

Е2

*

 

 

 

 

 

 

 

Чтобы доказать б), введем в рассмотрение множество

Z)_p = <= D:

р(я,

dD) >

Р}.

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

малых

Р это множество не пусто и

inf

F(0, х) <

F0+ft/4. Если

я е у и е , р >

 

 

зсео—р

 

 

0 достаточно малы, то

 

 

Рж(min {t: Х\ е

D _p( < til >

exp (— е-2 (F0 +

hl2)\,

 

 

 

Ti

 

 

 

 

inf

M* f XD (X s8) ds> a2 > 0.

 

 

oceD__p

0

 

 

 

 

Собирая эти оценки вместе, из (4.17) при достаточно малых 8 получим

ре (Z)) ^ сг•min Рх {min [t:X*<=D$} < TJ х

х^у

Tl

X inf М А XD ( X S) ds > се exp {— e“ 2(F0 + h/2)} •a2.

о

8 5]

ГАУССОВСКИЕ ВОЗМУЩЕНИЯ ОЁЩЕГО ВИДА

179

Отсюда,

если

учесть

(4.21)г вытекает утвержде­

ние б).

 

 

 

 

Отметим, что если поле Ъ(х) имеет устойчивые положе­

ния равновесия,

отличные от нуля, то утверждение теоре­

мы 4.3, вообще

говоря, неверно. Поведение инвариантной

меры в случае поля Ъ(х) более общего вида рассматривает­ ся в гл. 6.

§ 5. Гауссовские возмущения общего вида

 

Пусть функция Ъ(х, у), х е

# г, у е

R1, со значения­

ми в Rr такова, что |Ь{хц уг)— Ь(х21 у2)\ < К-{\хг х2\+

+

1 УгI)-

 

X* — X* (х), t ^ О,

 

Рассмотрим случайный процесс

который удовлетворяет дифференциальному уравнению

 

X? = b {X le Q ,

Xl=*x,

(5.1)

где

Lt — гауссовский случайный процесс в R1.

среднее

 

Будем считать, что процесс

£* имеет нулевое

и непрерывные с вероятностью 1 траектории. Для непре­ рывности, как известно, достаточно, чтобы матрица кор­

реляций a(s, t) = (aij(s, t)), ali(s, t)

=

имела дваж­

ды дифференцируемые элементы.

Обозначим b(x) = Ъ(х%

0). Процесс Х\ можно рассматривать как результат ма­ лых случайных возмущений динамической системы xt =*

= xt(x):

*

xt = b(xt), x0 = x.

В гл. 2 было показано,* что при е -* 0 равномерно на каждом конечном отрезке времени X* -* xt. Там же было изучено разложение X? по степени малого параметра е.

В этом параграфе мы займемся большими уклонениями X® от траекторий динамической системы xt.

Для непрерывных функций cps, s > 0, со значениями в R1 определим оператор S^cp) = и по формуле

t

u = ut = [b {u s, фa)ds + x.

о

180

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

[ГЛ. 4

Иными словами,

ut = В^Дф) — это

решение

уравнения

(5.1), где

вместо

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

и берется

начальное

условие

и0 =

х.

Через оператор

Вх можно

записать:

х? =

обозначать | ||с и | \\L*

нормы в

C0T{Rr)

и

Будем

LQT(R1) соответственно.

 

 

i/)

Л е м м а

5.1.

Предположим, что функции Ь(х,

удовлетворяют условию Липшица с константой К\ и =

= В^ф), v = В^ф), гдг ф, ф е C0T{Rr).

Тогда

 

 

 

||И- н||с < к

У Г гкт ||ф -

41к>.

 

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

По

определению

операто

ра Вх

имеем

 

 

 

 

 

VA «

(

 

 

 

 

 

J [&(“ .. Ф5) ~

b(vt, Tt>,)]ds

 

 

 

 

 

 

 

t

 

 

 

 

v, |ds +

*J| ф.

4 . 1ds-

 

 

 

 

 

о

 

Отсюда с помощью леммы 1.1 гл. 2 заключаем, что

 

 

т

 

 

 

 

I U -

v||с <

ектк J I ф, -

t , I ds <

К ут ект\Ф -

ф Ц..

 

 

о

 

 

 

 

Обозначим А корреляционный

оператор процесса £3.

Он действует в пространстве L Q T ( R 1)' по формуле

т

Aq>t = ja (s1t) ysds.

о

Как доказано в предыдущей главе, функционал дейст­

вия для

семейства процессов е£* в пространстве LOT (R1)

при с

0 имеет вид

 

<5'2>

Если А 1/2 фне определено, то считается, что 5ог(ф) — + 00•