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

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

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

§ 7)

ЗАДАЧИ

О СОБСТВЕННЫХ ЗНАЧЕНИЯХ

271

щихся

к i f 4,

HmP£{A'®/ _ 2\е

Г} = 1

Для

открытых

Г ZD К 4, если

е-*0

д

'

____

 

 

0 <

lim е2 In t (г~~2) ^ lim e 2 In t (г~2) < 2;

lim Р* {Х^е_ 2)е Г } =

е-*0

для

8">0

Т ZD К& при

1

открытых

2 < lim е2 In t (е~2) <1 lime2 In t (е—2) < 4; если эти нижний

Т-*о

Е_>°

 

 

и верхний пределы больше 4 п меньше 5, то lim Р* [Хг,

2)^

 

е->0

v

е Г } = 1 для открытых множеств, содержащих К2. На­ конец, если lim е2 In £ (е—2) > 5, то распределение Х^е__2)

е-*о

К3 независимо от начальной точ­

будет притягиваться к

ки X.

 

схему доказа­

Общую формулировку результата и

тельства можно найти

в статье Ф р е й

д л п н а [10].

Субпредельные распределения за экспоненциальное вре­

мя устанавливаются и в случае процесса X? в

области,

останавливающегося

по

достижении границы

(ситуа­

ция § 5).

что вопрос

о предельном распределении

Заметим,

Агде_ 2) тесно

связан

с задачей о стабилизации

решений

параболических дифференциальных уравнений с малым параметром (см. ту же статью).

§ 7. Задачи о собственных значениях

Пусть L — эллиптический дифференциальный опера­ тор в ограниченной области D с гладкой границей dD. Наименьшее по модулю собственное значение Ах операто­ ра — L с нулевыми граничными условиями, как извест­ но,— действительное, положительное и однократное. Оно допускает следующую вероятностную характериза­ цию. Пусть (X*, Рх) — диффузионный процесс с произво­

дящим

оператором L, т — момент

первого выхода X t

из области

D. Тогда

значение

Ах — пограничное между

теми А,

для которых

<

оо,

и теми, для которых

М л^ х =

оо

(см. Х а с ь м и н с к и й [2]).

Результаты, касающиеся функционала действия для

семейства

процессов

(Х*г Р£)г

соответствующих опера-

272

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

[ГЛ. б

торам

 

 

 

L* ~ \ 2 atl (*) dxidx) + 2 b*(*)дх1

(7.1)

могут быть применены к нахождению асимптотики собст­ венных значений — Ьг при е 0. Здесь возникают два качественно различных случая в зависимости от того, вы­

ходят все траектории динамической системы xt — b(xt) из D\JdD, или же в D имеются устойчивые со-предельные

множества

этой

системы.

случай рассмотрен

в

статье

Первый,

более

простой

В е н т ц е л я

[6].

Здесь,

как было

установлено

в

§ 1

(лемма 1.9), вероятность того, что траектория Xf

проведет

в D больше времени Г, оценивается сверху выражением

ехр { —г~2-с(Т — Г0)}, с >

0.

Отсюда вытекает,

что

МЕеЯтЕ < оо при Х<&~2-с и,

таким

образом,

 

 

2с.

Более точно скорость стремления ЯЕк бесконечности

ука­

зывает

 

 

 

 

 

 

 

 

 

Т е о р е м а

7.1. При е ->■ 0 имеем: Xе= (сх +

о (1)) х

X 8 2, где

 

 

 

 

 

 

 

 

 

сг — lim Т~хm in{S0T (ф): ф* е

D U dD,

 

(7.2)

T->OQ

 

 

 

 

 

 

 

 

 

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

семейст­

ва диффузионных процессов (Xf, Р*).

всего установим

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

 

Прежде

существование предела (7.2). Обозначим минимум,

участ­

вующий в (7.2), через а(Т). Прежде всего, ясно, что а(Т/п) ^ а(Т)/п при любом натуральном п (потому что значение функционала £(ф) хотя бы на одном из отрезоч-

ков [кТ/п,

+

i)T!n\ меньше или равно

S^T(q)ln).

Отсюда для любых 71, Т >

0

 

 

 

а(Т) >

[Т/Т1а(Т).

(7.3)

Чтобы получить оценку противоположного знака, вос­

пользуемся

тем,

что существуют положительные Т0 и А

такие, что любые две точки х и у из D{JdD можно соеди­ нить кривой y t(x, у), 0 < * < Т(х, у) ^ Т0, со значением

§ 71

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

273

функционала SoT(XiV) (q)(x, у ) ) ^ А . Пусть

при каком-

то большом Т минимум а(Т) достигается на функции ф.

Полагаем х =

фт,

у =

ф0 и составляем

функцию ф* из

кусков:

ф* =

ф, при

0

£ <1 Т \

ф* =

ф*_т (я, У)

при

Т ^ t ^

Т +

Т(х,

у);

 

а далее продолжаем ф периоди­

чески с периодом Т +

 

Т(х, у). Получаем для произволь­

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

Т

 

 

 

S 0T (Ф) <

( [ Г /( ? +

Т (*,?))] + l ) ( S f f ($ )

+

 

 

 

 

 

 

+

S0T(X,v) (<р(х, у))),

(7.4)

a(T)^([T/T] +

1 )(а (Г ) - М ) .

 

 

(7.5)

Деля (7.3), (7.5) на Т и переходя к пределу при Т -> оо, получаем

? “ ‘ e(7 ,) < llm Г ~ 'а (Г )< П Ж T~i а (Т) 7’-1 (а (?)+ Л ). г—

(7.6)

Переходя затем к пределу при Т -> оо, получаем сущест­

вование предела (притом конечного).

выхода

процесса

Пусть теперь т8 — первый момент

X] из D. Докажем, что М* ехр {е~26те} = сх> при

доста­

точно малых е, если b >

 

Легко видеть, что

 

 

lim T~l

sup

min {£ оТ (ф): Ф0 = х, Фг =

у, ф( <= £> (J

<?£>,

Т-+оо

K,y^D

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

также равен

сг.

Выбираем

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

меньшее

— сх)/3, п такое

Г, что

 

 

 

 

 

T~l sup

min {SoT(ф): ф0 =

х, Фт = У, Ф< ^ D U dD,

 

x,yeD

 

 

 

О

T } < c L+ K. (7.7)

 

 

 

 

 

Теперь выберем такое 8 >

О, чтобы

 

 

 

min (5 оГ (ф): ф0 =

х1 фт = у,

Ф ,е В - { U dD-б, 0 < i <

 

 

 

 

 

< T } < T ( Cl +

2x)

(7.8)

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

при всех х, у ^ 0 6, где 0 _ 5 — множество точек в 0 на расстоянии, большем чем 6, от границы. Это можно сде­

лать в силу леммы (1.4) (Т = Т, потому что функции ф* определены на одном и том же отрезке [О, Т ]; см. доказа­

тельство

этой

леммы). Для

произвольных

х е

0 _ й,

у е

0 26

выбираем соединяющую их кривую

ф*,

0 ^

^

t <1 Т,

целиком проходящую в D-b\J dD-b, для

кото­

рой SoT(ф) ^

Т(сг + 2х). По

теореме 3.2 гл. 5

получаем

при достаточно малых е для всех х е 0 _ 6*.

 

 

Р* (те > Т, Хгт<= Z)_e} > Ре {рйГ (X е, ф) < 6} >

> е х р { — е - 2 [50Г (ф) + K r il^ e x p f — г~2Т (сх+3х)}.

 

 

 

 

(7.9)

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

(п — 1)-кратное

применение мар­

ковского свойства

дает

 

 

Р* {те > пТ} ^ ехр { — пг~~2Т (сх+ Зх)}.

Отсюда при малых е для всех п

 

 

Мл-ехр{е“ 2Ьте}^ е х р {е “ 2б7гГ}Р®{те >

п Т }^

 

 

^ exp {m ~2T (b с1— Зх)},

что стремится к бесконечности при

п->- оо. Поэтому

М*ехр (е—26те} =

оо.

 

Теперь докажем, что если Ъ<

то

М*ехр {е—26те} <

< оо при малых е. Опять берем 0 <

х < (сг Ъ)!3; бе­

рем такие Т и б >

0,

что

 

 

min {£ от (ф): ф, е

0+6 U 50 ■6, 0 < t < Г } > Т (с1— 2х)

 

 

 

 

(М О )

(еще раз применяем лемму 1.4). Расстояние между мно­ жеством функций ф*, целиком лежащих в D при 0 ^ t ^

^

Г, и любым

из множеств

2 х)) = {ф.*фо =

=

х, SoT(ф) ^

T(ci — 2х)} не меньше б; поэтому в силу

верхней оценки теоремы 3.2 гл. 5 при достаточно малых е для любого х е 0

Р* {те > Г ) < Ре {рог (X е, Фх (Т (с, - 2х))) > 6} <

< е х р { - е - * 2 , (е1 - 3 х ) } . (7.11)

§ 7] ЗАДАЧИ О СОБСТВЕННЫХ ЗНАЧЕНИЯХ 275

Марковское свойство дает нам Р* {т8> /г Г }^ ехр { —пг~~2Х

X Т {сх— Зх)}, п

М£ехр {е—2Ьте} ^

с»

<

2 ехр {г~ 2Ь{п + 1) Т} Р1{пТ < тЕ<; (га -f

1) Г } <

 

п=0

 

 

 

 

 

<

2

ехр {г~2Ь(п + 1) Т) Рх {тЕ> пТ}^.

 

 

п—О

 

 

 

 

 

 

 

оо

 

 

^

ехр {е,~2ЪТ} 2 ехр {пе~2Т (Ь— сх +

Зх)} < с».

 

 

 

 

п=0

 

Теорема

доказана.

 

Получить широкий класс примеров, в которых предел

(7.2)

вычисляется,

помогает

 

Т е о р е м а

7.2.

Пусть поле Ь{х) потенциально отно­

сительно римановой метрики, связанной с матрицей диф­ фузии, т. е. существует гладкая в D[JdD функция U(x)

такая,

что

 

 

 

 

 

 

Ь1(*) = - 2

aii (*) ^ г -

 

(7.12)

 

 

i

 

 

 

Тогда

равио

 

 

 

 

 

min

 

i

(7.13)

где матрица (аи(х)) = (а^(х)Уг.

 

 

 

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

Обозначим

минимум

(7.13)

через с*. Пусть минимум а(Т) достигается

на функции ф:

 

т

(«Р») ( ф*- V(ф() )(ф/ — Ь; (ф,)) Л.

 

асП = SУ И

(7.14)

 

О

 

 

 

 

Преобразуем этот интеграл:

 

 

 

 

т

 

 

 

 

а (Г) =

1 т 2

(ф*) Ф Wtdt —

 

 

 

о

 

 

 

 

Г

 

(ф4) фIdt +

т

 

(Ф,) dt.

J Z аи (Ф«)

f1 2 в» (Ф,) 6* (Ф()

о

276

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

ВРЕМЕНИ

[ГЛ. 6

Первый

интеграл

здесь

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

второй

в

силу

(7.12)

равен

£7(ф0) — £/(фт),

а третий не меньше чем с\у

умноженное

на Т. Отсюда

вытекает,

что

Т

[а(Т) ^

^с\ Т~*6\

где

С — максимум U(x) U(y)

по

всем

I , I/ G

D(JdD. Переходя к пределу при Т

 

оо, получаем

функции

Неравенство

 

вытекает

из

того,

что

для

cpt == х0,

где х0 — точка, в которой

достигается

минимум

(7.13),

T ^S QT (ф) = с*.

 

случай — ког.

Теперь рассмотрим

противоположный

да в D есть со-предельные множества системы xt = b(xt)- Предположим, что выполнено требование А § 2 (существо­ вание конечного числа компактов Ки . . ., К г, состоящих из эквивалентных друг другу точек, содержащих все со-предельные множества). В этом случае, как мы уже гово­ рили, наш процесс с малой диффузией хорошо прибли­ жается цепью Маркова с (I + 1) состояниями, отвечаю­ щими компактам K t и границе области dD, с вероятностя­ ми перехода порядка

ехр { -г -* У в (Кь К;)},

(7.15)

exp { - ^ V D{Kb dD)}

(вероятности перехода из 0D в Kt полагаются рапными О, диагональные элементы матрицы вероятностей перехода такие, что суммы но строкам равны 1).

Представляется правдоподобным, что условие конеч­ ности М х^ 8 близко к требованию конечности математи­

ческого ожидания е^8, где Vе — число шагов в цепи с указанными переходными вероятностями до первого достижения состояния dD. Пограничным между теми Я, для которых это математическое ожидание конечно и бес­ конечно, оказывается взятый с обратным знаком логарифм наибольшего после единицы собственного значения мат­ рицы вероятностей перехода.

Асимптотика собственных значений матрицы с эле­ ментами порядка (7.15) находится в заметке В е н т ц е л я 13]; обоснование перехода от цепи к диффузионному про­ цессу Xt — у В е н т ц е л я 12 ].

Т е о р е м а

7.3. Пусть Я8, Я8, ... Д? — занумеро­

ванные в порядке

убывания модуля собственные значения

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

277

(iисключая собственное значение 1) стохастической матри­ цы с элементами, логарифмически эквивалентными вы­

ражениям (7.15) при е -> 0. Определим константы

FW,

к — 1,

. . .,

Z, l

1,

формулой

 

 

 

 

 

 

 

VW =

min

2

VD(a,p),

 

(7.16)

 

 

 

gE-GW

 

 

 

 

 

 

где G^

множество

всех

W-графов

над множеством

L — {К!,

К t, dD) таких, что в множестве W содер­

жится к элементов. (Заметим,

что

при к = I +

1

мно­

жество

W =

L;

при

этом

имеется

ровно один

IF-граф,

а именно пустой, и сумма в (7.16) равна нулю как и F^+1).)

Тогда при е

 

0

 

 

 

 

 

 

 

 

Re (l — Xl) ж; exp ( — (F (W— V(ft+1)) б—2j.

 

(7.17)

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

аналогично

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

леммы 3.1 и других лемм § 3: коэффициенты характеристи­ ческого многочлена матрицы вероятностей перехода, из которой вычтена единичная матрица, выражаются в ви­ де сумм произведений n(g) по PF-графам, откуда находит­ ся асимптотика этих коэффициентов; из нее выводится

асимптотика

корней.

 

Заметим,

что

из утверждения теоремы вытекает, что

7 <а- 1) — учю

 

V(b) — F ^*1) < V(k+v — F<A+2F

В слу­

чае «общего положения», когда эти неравенства

стро­

гие, собственное

значение к\ при малых в оказывается

вещественным и

однократным (старшее собственное зна­

чение Х\ всегда вещественно и однократно).

Что касается обоснования перехода от асимптотической задачи для цепи Маркова к задаче для диффузионного про­ цесса, оно разбивается на две части: конструкцию, связы­ вающую диффузионные процессы с определенными диск­ ретными цепями, не зависящую от наличия малого пара­ метра е; и оценки, связанные с цепями Маркова Z n, вве­ денными в § 2.

Для первого оказывается существенной следующая лемма.

Пусть (Хи Рх) — диффузионный процесс с производя­ щим оператором L в ограниченной области D с гладкой

278

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

границей; dg—-какая-то поверхность внутри D, Г— повер­ хность какой-то окрестности dg, отделяющая поверхность

dg

от

dD.

Введем случайные

моменты

времени

а0=

min

{ t : X t е

Г}, Тх =

min {t ^

а0: X t е

dg\JdD}.

Определим для а >

О операторы </а, действующие па огра­

ниченные функции

на dg:

 

 

 

 

 

 

qaf ( x) = Mx [XXt(=dg-,eax'f(XXi)}l

x<=dg. (7.18)

Л е м м а

7.1.

Наименьшее собственное значение

Xi

оператора L с

нулевыми

граничными

условиями

на

dD пограничное между теми а, при которых наиболь­ шее собственное значение оператора qa меньше 1, и теми, при которых оно больше 1.

Д о к а з а т е л ь с т в о этого утверждения и его обобщение на другие собственные значения можно найти в статье В е н т ц е л я [5].

Чтобы применить это к диффузионным процессам с ма­

лой диффузией, мы выбираем dg = |Jdg*ii Г =

(JI\ так,

г

г

как в § 2, и доказываем для соответствующих операторов

ql следующие оценки.

 

 

 

 

 

Л е м м а

7.2. Для любого у > 0 можно выбрать ра­

диусы

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

компактов

K t настолько

малыми,

что при всех а ^ е~уг~2 и достаточно малых е

 

1 + ае~уг~2

<

1 +

aev°~2,

x s % ;

(7.19)

exp { -

( Ь (Kh Kj) +

у) e“ 2) <

qhogj (x) =

 

= К

(X?, e

dgj; eaxt) < exp 1 -

D (Kt,

K3) - y) e~2},

 

 

 

 

 

xe£dgt;

(7.20)

exp { -

(VD(Klt dD) +

V) e~ 2) <M * {X?, s

dD; eaXi] <

< e x p { — (VD(Ki, dD) — v)e~ 2),

x ^ d g t.

(7.21)

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

аналогично

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

лемм 2.1 и 1.7; верхние оценки получаются с несколько большим -трудом.

§ 7]

Задачи о собственных значениях

279

Окончательно получается такой результат:

 

Т е о р е м а 7.4. Наименьшее собственное

значение Я®

оператора

Ьг при е

0 логарифмически эквивалентно

 

ехр { —e"2(V(D - F<2>)},

(7.22)

где F<x>, F<2) определяются формулой (7.16).

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

7ч ^ ехр ( — е~2 min V (О, г/)},

(7.23)

y^dU

 

где F — квазипотенциал.

Теоремы 7.1, 7.4 не дают ничего содержательного в случае, когда внутри!) есть только неустойчивые предель­ ные множества или когда устойчивое предельное множест­ во находится на границе области. Некоторые результаты, относящиеся к этим случаям, содержатся в статье Д е- в и н а ц а, Э л л и с а , Ф р и д м а н а [1 ].

Результаты, содержащиеся в теоремах 7.3, 7.4, родст­ венны результатам § 6. В частности, константы, при пере­

ходе через которые

предела

lim е2 In £(е“2)

происходит

смена

одного

субпределытого

е-+0

другим,—

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

не что

иное,

как F

— F^H).

 

Поясним механизм этой связи. Пусть (Xf, Р£) — диф­ фузионный процесс с малой диффузией на компактном многообразии М или процесс в замкнутой области D\JdD, останавливающийся по достижении границы; f(x) — не­ прерывная функция, обращающаяся в нуль на дБ. Допус­

тим, что функция ие х) = Ml-/ (А7*6) — решение соот-

д н г

ветствующей задачи для уравнения — = Ьгиг — раз­

лагается в ряд по собственным функциям е\{х) опера­ тора —

игх) = 2 с\е Xlt£ (х)

(7-24)

i

 

с коэффициентами cf, получаемыми интегрированием Д умноженного на соответствующие собственные функции

е\ сопряженного оператора. (Так как оператор Ь г не

280

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

самосопряжен, предположение о возможности такого разложения не выполняется автоматически.) Допустим, что первые 10 собственных значений вещественны, неотри­ цательны и экспоненциально быстро убывают:

ж ехр {— С # " 2}, + °о > Сх> С2> ... > Cto.

Пусть функция £(е“2) такова, что Cj+i<lim e2ln t

2) ^

 

е-»и

 

< Time2 In t(s~2) < C }; тогда

Hm<>- ^ '(e 2) = 1

При

е-^0

е-*0

 

и этот предел равен нулю при i > /. Если бесконечный хвост суммы (7.24) не помешает этому, то получим

Пт иг (t (е~2), х) =

2 Пт f / {у) е* (у) dy-lim е\(ж). (7.25)

е-*0

 

г—1 е-*0 J

е<-*0

 

Иначе

говоря,

характер

предельного

поведения

цр(£(е“2), х) меняется при переходе lim е2 In t (е“ 2) через С]. е->0

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