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

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

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

i п

НЕ ОЧЕНЬ БОЛЬШИЕ УКЛОНЕНИЯ

331

няется

соотношение

 

 

П т

In М ехр е*~‘

(as,/ (Is/e)) tfej =•

 

по

0

т

 

 

 

 

 

= ^ ( C a s,a s)ds,

(7.1)

 

 

О

 

где С — некоторая симметричная постоянная матрица порядка г. Легко видеть, что из (7.1) при as == а вытекает равенство

lim Г2к

ЧпМехр

(а, /(5,)) ds

- у

(Са, а).

(7.2)

00

 

 

 

 

 

При

определенных

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

которые

мы

уточнять не будем, матрица С является матрицей вторых производных по переменным а от функции Н(х, а), введен­ ной в § 4, при х = 0, а = 0.

В дальнейшем мы будем предполагать для простоты, что det С Ф 0. Если матрица С вырождается, но какое-то собственное значение отлично от нуля, то по существу все результаты сохранятся, только их формулировки ста­ нут более громоздкими.

На функциях из CoT(Rr) определим функционал

г

S (ф) = S OT (ф) = - у j (С"‘(ф, - Я Ф ^ Ф. - £ ф 4)d st

О

если cps абсолютно непрерывна; на остальных <p е CoT(Rr) полагаем S(ср) = +оо. С таким функционалом мы уже встречались в гл. 4; S(ср) — нормированный функционал

действия для семейства гауссовских процессов ^

опреде­

ляемых дифференциальным уравнением

 

+

(7.3)

где Wt — процесс белого шума. Ниже мы поясним, какое отношение к нашей задаче имеет уравнение (7.3).

332

 

д п н а 'м и ч е с к и е си стем ы с

у сре д н е н и е м

 

1ГЛ. 7

Т е о р е м а

7.1.

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

что выполнено

усло­

вие Fx,

и det С Ф 0.

Пусть для

какого-то

у >

1 — 2х

при

любом

б >

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

 

 

л

 

_

 

 

 

lim evln Р

sup

Г*

( Ы

-

В) ,f e(S_U)B/ (|ц/е) du ds >

е-*0

 

О < t< T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

> б | = — 00,

 

(7.4)

и пусть, кроме того, существуют

t0 е (О,

Т ]

и функ­

ция 0(«), 0 (0 >

0 , Пт 0 (0 =

0 ,

такие, что

 

 

 

 

 

 

 

по

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8х-

1

ft-ft

 

 

< 00.

lim

sup

е1

2х In М exp

 

 

/( U )

 

е|0

е ^t<l„

 

 

 

в(<)•j

 

 

 

 

0<h < T —t

 

 

 

 

 

 

 

 

 

 

Тогда

&2*~lSoT(q>) будет функционалом действия при

е 1 0

в

пространстве CQT{R)

для

семейства процессов

Ъ\ =

e"~HXf,

к е (0,

1/2),

где

Xf — решение

уравнения

(1.2)

с

начальным условием

X Q =

0

(напоминаем,

что

6(0)

=

0).

 

 

 

этой

теоремы представляет

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

собой комбинацию некоторых рассуждений, примененных при доказательстве теорем 3.1 и 4.1, поэтому мы его толь­ ко наметим. Во-первых, доказывается, что е2х_1*5(ср) будет функционалом действия для семейства процессов

если e2x_1 S(ср) — функционал действия для про­

цессов ~Z\= e“”xXf, где Xf — решение линеаризованной системы

х ! = п ъ/ г) + в ( 11/г) х 1 X Q = о.

Далее, используя (7.4), доказывается, что оценки для Z®, которые входят в определение функционала действия, выполняются, если они справедливы для процесса

Z\ = гГ’>лХ]1 где X® — решение уравнения

Х1 = /(Ы + Ъх1

( 7 .5)

§ 11

НЕ ОЧЕНЬ БОЛЬШИЕ УКЛОНЕНИЯ

333

Чтобы получить оценки для Xf, сначала вычисляется

функционал

действия для

семейства

процессов

тр =

= е~к j* / (£*/е) ds. Отсюда с

помощью

теоремы 3.1

гл. 3

о

 

 

 

 

находится функционал действия для процесса Xf.

Нако­

нец, вычисление функционала действия для процессов

г), проводится так же, как это делалось в лемме 4.3 (см. также (Г е р т н е р 11])). При этом используется послед­ нее условие теоремы.

Таким образом, функционал S((р) и нормирующий коэффициент /(г) = е2х-1 характеризуют отклонения

порядка 8х, х е (0, 1/2). Процесс K_1/2Xf в соответствии с теоремой 3.1 при е->- О сходится к гауссовскому процес­

су

В нашем случае этот гауссовский процесс

как

легко

проверить, удовлетворяет уравнению

(7.3)

при

X = 1.

Отношение

X f/ex можно

записать

в

виде

е1/2—х-e~~1/2Xf и

интерпретировать

результат

теоре­

мы 3.1 следующим образом: вероятности уклонений по­ рядка 8х, х е (0, 1/2), имеют ту же асимптотику, что и уклонения порядка 1 под действием гауссовского шума

XCw, X = 81/2~х. Этот факт, конечно, находится в пол­ ном соответствии с тем, что большие, по не очень большие, уклонения для сумм независимых слагаемых имеют ту же асимптотику (в главных членах), что и соответствующие уклонения нормального приближения (см. Л и н н и к,

И б р а г и м о в [1 ]).

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

Т е о р е м а 7.2. Пусть Е£ — стохастически непрерыв­ ный марковский процесс с конечным числом состояний;

Pij(t) его переходные вероятности; qi}

dt

о

 

Предположим, что все qtj Ф 0.

Обозначим Х(а) собствен­

ное значение матрицы Qa =

(g?j), =

а , Ь(0, i))9]

имеющее

наибольшую вещественную часть.

Тогда условия

теоремы

7.1 при х е (0,

1/2)

выполнены,

и С =

где

=

д-1 (а)

 

 

 

да. да . ,а=о

 

 

 

334

ДИНАМИЧЕСКИЙ СЙСТЕМЫ С УСРЕДНЕНИЕМ

[ГД. ?

 

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

этой

теоремы использует

конструкцию, примененную

при

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

тео­

ремы 4.2. Прежде всего, так же, как это сделано в лемме

4.3, доказывается, что условие F* в случае процесса ^ выте­ кает из соотношения (7.2). Чтобы проверить (7.2), заме­ тим, что

М,- ехр \ г* j (а, / (5,)) ds\ = (т\ *al) (i)t

(7.6)

где T f— полугруппа операторов, действующая в прост­ ранстве функции g(i) на фазовом пространстве процесса It по формуле

№ ) ( 0 =M ,g (&,)«“

Используя обозначения, введенные при доказательстве теоремы 4.2, получим

г2х_1 in с + г2*-1 in (гГх<х1) (0<

< t2*~l In (Т\ *°e)(l) = t2*\ (Г « а) + t2*~l In е, <

 

<

i 2K“ 1 l n (т‘,

Заметим, что 1 — наибольшее по модулю собственное

значение оператора Т? при а =

0,

поэтому Х(0) = 0.

Принимая во внимание, что 6(0) =

0, нетрудно проверить,

что все первые производные функции Х(а) при а = 0

обращаются в нуль, поэтому X (а) =

--- ^

д а 'доГ^ aiaJJr

+ o(|a|2) при |a|~>0. С учетом

ij

*

J

последнего

равенства

щз (7.7) вытекает, что

 

 

 

Нгп^-Чп ( Г р а1)(0 =

 

 

 

L "V д2К

(0) а;0С/ =

-2 " (Со* а).

2

§ 7] НЕ ОЧЕНЬ БОЛЬШИЕ УКЛОНЕНИЙ 335

Отсюда и из (7.6) следует (7.2), а стало быть, и условие Fx. Проверку остальных условий теоремы 7.1 мы предо­ ставляем читателю.

Аналогичный результат верен, конечно, и для некото­ рых других марковских процессов, например, если — невырожденный диффузионный процесс на компакте. Близкие по существу рассуждения позволяют проверить условия F и Fx и для некоторых немарковских процессов с хорошими свойствами перемешивания (ср. С и н а й [1 ]).

Выпишем уравнение Гамильтона — Якоби для функ­ ции u(t, х) = inf {SоДср): фо = 0, ф* = х}\ как неодно­ кратно говорилось, через эту функцию выражаются глав­

ные члены

асимптотики

многих

интересных

величин.

Вычисляя

преобразование

Лежандра

от

 

функции

L(x,

$ ) = ± . { ( Г '{ $ - В х \ $ - В х ) , придем к следую­

щему уравнению для u(t, х):

 

 

 

 

 

 

 

\

(CVxu>Vxu) +

(Вх, У*и).

 

 

Можно выписать и уравнения Эйлера для

экстремален;

в рассматриваемом

случае

эти

уравнения

линейны.

Из теоремы 7.7, в частности, следует, что при к е

(0, 1/2)

lime1-2*]nP { sup

| X f| > ex*d\ =

 

 

 

 

e-tO

 

 

 

 

/

 

 

 

 

 

 

 

= — min (£, x) : |x |= d1 i e

[0, Г]}.

Уклонения порядка ex определяют среднее время, ко­

торое

понадобится

траектории

процесса

X8, X Q= 0,.

для того, чтобы впервые

выйти из

области Z)e,

содержа­

щей точку 0, если De получается сжатием в ех раз из не­

которой

фиксированной области D: Dh= 8х •/).

Пусть

 

 

т8 = inf

U iX f e O 8}, V(x) =

inf u{t,x).

 

 

 

 

 

t>о

 

Т е о р е м а

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

7.2,

матрица С не вырождается,

причем Ds — ех-£>,

х е

(0,

1/2),

где D ограниченная область с

гладкой

ш

ДИНАМИЧЕСКИЕ

СИСТЕМЫ

С УСРЕДНЕНИЕМ

[ГЛ. 7

границей. Положим F0 =m inF (;r).

Тогда

 

 

 

 

x&OD

 

 

 

 

 

lim g1

2х In М0т8 =

F0,

 

 

 

£i О

 

 

 

 

 

 

lim Р0 lee2x~ 1(r<l_v) <

т£ <

ее2х_1(У«Н)] = 1

 

ею

 

 

 

 

 

 

при любом у > 0.

 

этой

 

теоремы

аналогично

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

 

доказательству соответствующих

результатов

из

гл. 4,

и мы его

опустим.

слов

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

вычисления

Скажем несколько

функции V(x). Аналогично тому, как это делалось в § 3 гл. 4, можно доказать, что функция V(x) есть решение задачи R0 для уравнения

~ (C W (X), VV (х)) + (Ъх, VV (X)) = 0.

Непосредственной проверкой нетрудно убедиться в том, что если матрица С~1В симметрична, то F(:r) =

=*— (С^Вх, х).

Взаключение этого параграфа заметим, что аналогич­ ные оценки можно получить не только для отклонений порядка ех, х е (0, 1/2), от положения равновесия, но п для отклонений порядка 8х от любой усредненной траек­ тории.

§8. Примеры

1.Рассмотрим сначала случай, когда правые части уравнения (1.2) не зависят от х и процесс £* стационарный.

Тогда X] можно представить в виде

 

 

 

t

ае

 

 

X] = х +

j* Ъ(S,/e) els = х +

t (fye)-1 j

b (Ss) ds.

 

 

 

 

о

0

 

 

Обозначим m =

Mb(£s), В1>(т) =

М(й*(£5+Т) — т')(Ь1(£3)

mi)

II

 

г

в н (т) -*■ 0

при т

оо.

предположим, что

Отсюда

с

 

i=l

 

 

любом

помощью неравенства Чебышева прн

§ 8]

 

 

 

ПРИМЕРЫ

 

 

837

б > 0

получаем

 

 

 

 

t+T

 

 

 

 

 

■jr

j*

ь (|s) ds — m > 8

 

 

 

 

 

1+ T l+T

I M(*‘(5.) ~ " W

 

 

 

Иг*

j

j

(5u) -

«О* <*» -

 

 

t

«

i=1

 

 

 

 

 

 

 

И-T t+ T

T

 

 

 

 

 

 

= 5 ^ ] j ^ В п ( и - s)duds~* 0

 

 

 

 

< *

i=1

 

 

с ростом

T равномерно no t^ O . Таким образом,

в этом

случае выполняются условия теоремы 2.1,

и sup

! Х\

 

 

 

 

 

 

о<«г

 

х — иг£| -> 0 при е -> 0 по вероятности. Если процесс Е* обладает достаточно быстро убывающим коэффициен­ том перемешивания а(т), и функции Ь(у) растут не слиш­ ком быстро, например, если они ограничены, то примени­ ма теорема З.Т. В рассматриваемом случае эта теорема сводится к утверждению, что семейство процессов

 

 

 

 

 

 

V е

 

= ] / ё | (6 (| .) - т е ]*

 

 

 

о

 

слабо сходится при е -> 0 к гауссовскому процессу

имеющему

нулевое

среднее

и

матрицу

ковариаций

 

 

 

 

 

оо

(MlkD = (t Л s)-K = (t A s) (Klj),

где Kij =

J Bi}(x) dx.

Очевидно,

что

процесс

£*

имеет

—oo

независимые

приращения. Утверждение о сходимости распределения

(при фиксированном I) к гауссовскому составляет содержание центральной предельной теоремы для слу­ чайных процессов (см., например, Р о з а н о в [1], И б р а г и м о в и Л и н н и к [1]).

Если для процесса X] выполняется условие F из § 4, то применима теорема 4.1, которая дает возможность оце­

нить большие (порядка 1) уклонения Xet от линейной функции mt + х, t е [О, Т]. Так как правая часть урав­ нения (1.2) не зависит от х, то функция Н(х7 ос), опреде-

1 2 А А . Д , В е н т ц е л ь , м . И . Ф р е й д л и а

838 ДИНАМИЧЕСКИЕ СИСТЕМЫ С ОСРЕДНЕНИЕМ [ГЛ. 1

ляемля условием F, тоже пе зависит от х, и нормирован­

ный функционал действия £ оТ(ф) имеет вид: *5от(ф) ^

г

= J* L (cps) ds, где L(P) — преобразование Лежандра функ-

и

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

U(t, х) = ini{Sot{y): ф0 = х0, ср* = х}.

В случае, рассматриваемом здесь, эта нижняя грань и экстремаль, на которой она достигается, просто вычис­ ляются. Действительно, будем считать, что функция L(P) строго выпукла. Тогда уравнения Эйлера для функцио­ нала 5(ф)

- i r t W - o

показывают, что экстремалями являются только прямые

Ф = с, с е Rr. С учетом условий на концах отрезка [О, t ] получаем, что нижняя грань достигается на функции

Ф* = *0 +

и и (*> х) =

tL

 

 

Допустим, мы хотим найти асимптотику при е \0 для

 

t

 

 

 

 

InP {Xf е D], где

Xf = j b(^/£) ds, D — некоторое огра-

 

o

границей 0D.

Из

теоремы

ниченное множество в й г с

4.1 следует, что

 

 

 

 

 

limelnP { X f e D }

=

— £• inf

\ 1

 

ею

 

 

зеениан

J

если только эта нижняя грань совпадает с нижней гранью, взятой по совокупности внутренних точек множества D. Если, кроме того, нижняя грань достигается только в

одной точке х е D[JdD, то нетрудно доказать, что для любого б > О

H m P l i X ? - £ | < 6 | X ? s f l U dD) = 1 .

о

§ 8]

 

ПРИМЕРЫ

333

Пусть, например, процесс l t

определяется равенством

h =

Лi ПРИt G

U\ * +

1), * — целое,

(8.1)

где г|0, гц, . .

г\п, . . .

— последовательность

независи­

мых величин с общей функцией распределения F(x). Тогда условие F выполняется и

со

# ! (а) = In e{a'J * v))d F ( y ) f

—оо если только интеграл, стоящий под логарифмом, сходит­

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

но течению процесса Х\ на целом отрезке t е [О, Г], а не только для событий, относящихся к одному фикси­

рованному моменту времени.

имеет вид

 

2.

Пусть теперь уравнение (1.2)

 

 

Х? = Ь ( Х ? ) + о ( Х ? ) Ь /е,

Х8о =

ze= i?r,

(8.2)

где

Ъ{х) = (^(х), . .

Ьг(х)),

а(х)

= (а) (а:)),

г-мерный случайный

процесс, М£* =

0.

Функции

Ь1(х),

о) (х)

предполагаются

ограниченными

и достаточно

глад­

кимиЕсли диагональные элементы корреляционной мат­

рицы Z?(s, t)

процесса £* стремятся к нулю при 11 s\->

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

процесс X] при в \0 сходится по вероятности равномерно

на отрезке

0 ^ t <1 Т

к решению

дифференциального

уравнения

 

 

 

 

xt =

b(xt), х0 = х.

(8.3)

Если процесс обладает хорошими свойствами перемеши­ вания, то с помощью теоремы 3.1 можно оценить нормаль­ ные отклонения от xt: вычислить характеристики гаус­ совского процесса предельного для процессов = в ~1/2х

х ( х ? - 5 , ) .

12*

840

ДИНАМИЧЕСКИЕ СИСТЕМЫ С УСРЕДНЕНИЕМ

[ГЛ. 7

Предположим теперь, что для процесса l t выполняется

условие

F: существует функция Я^(а): Rr

R1 такая,

что для

любой ступенчатой функции as :

[0, Т ]

Rr

lim e In М exp ti о

и пусть функция Н|(сс) дифференцируема по а. Тогда, как легко проверить, условие F выполняется для уравне­ ния (8.2), и

Н(х, а) = (Цх), а) + Нъ(о*(х)а).

(8.4)

Преобразование Лежандра Ь(х, р) функции Н(х, а) прос­ то выражается через Ь^ф) — преобразование Лежандра для Н^(а):

Цх, р)

= /^(ст-*(д:)(Р — Ь{х))),

 

 

если только матрица сг(#) не вырождается.

 

 

Пусть, например,

— марковский процесс, принимаю­

щий два значения ех, е2 е

Rr; (pu{t)) — матрица

вероят­

ностей перехода,

ен . .

Как

доказано

в §

4,

 

0).

для процесса \t условие F выполняется, и функция

а)

равна наибольшему

собственному

значению матрицы

 

/?и +

(a. ei)

Я12

\

 

 

\Я21 Ягг + («, e2)J

Рассмотрим случай, когда дп = д22 = —д, ех = —е2 =

= е е Rr. Решая характеристическое уравнение, нахо­ дим, что

(а) = - q + у q%+ (а, ef,

и с помощью соотношения (8.4) получаем функцию Н(х, а) для семейства процессов X].

Предположим, что 0 — асимптотически устойчивое по­ ложение равновесия для системы (8.3). Для нахождения асимптотики среднего времени выхода из области, содер­ жащей точку 0, точки, через которую этот выход происхо­

дит, асимптотики инвариантной меры процесса X ®и дру­