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

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

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

i 8)

ПРИМЕРЫ

341

гих интересных

характеристик нужпо, как это следует

 

 

т

из § 6, вычислить функцию V (я) = inf { JL (фв, cps) ds: Фо35

о

= 0 , фг = я, Т > 0}. Эту функцию можно найти как ре­ шение задачи /?0 (см. § 4 гл. 5) для уравнения

Ф (я), VF (я)) — q + V V + (а* (з) VF (з), е)2 = 0.

Если система (8.3) имеет асимптотически устойчивый пре­ дельный цикл Г, то отклонения от этого цикла описывают­ ся квазипотенциалом Рг(я), который можно найти как решение задачи Яг (см. § 4 гл. 5) для того же уравнения.

Рассмотрим в этом примере отклонения порядка ех, х е (0, 1/2), от положения равновесия 0. Если для про­

цесса

выполнено

условие Fx,

т. е. если

существует

предел

 

 

 

 

Ит е 1

2х In М ехр

 

 

as) dst

е ; о

 

 

 

 

где Сi — некоторая

симметричная

матрица,

а8 — произ­

вольная ступенчатая функция на [0, Г], то, как легко проверить, условие F* выполняется также и для процес­

са

Xf,

определяемого

уравнением

(8.2),

и

матрица С =*

=

а(0)С^а*(0).

Пусть

— марковский

процесс с

двумя

состояниями, о котором шла речь выше. Тогда, как сле­

дует из теоремы 7.2, условие F*

выполняется, и

=*

 

д^Нг (а)

= -7 - ( eV ), где е1,

е\

 

ет— коор­

 

д а . д а .

 

 

i

J а=0

 

Таким образом,

для

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

динаты

 

вектора е.

цессов

Х\ в этом случае получаем С =

а (0) (еV ) а* (0).

 

Легко проверить, что условия F и Fx для уравнения

(8.2) выполняются и в случае процессов

£*, заданных ра­

венством

(8.1).

 

 

 

 

 

 

 

3.

 

Рассмотрим уравнение Ван дер Поля со случайными

возмущениями:

 

 

 

 

 

 

 

 

 

х +

со2# =

е [/(я, з, vt) +

ф(з, я)£*).

( 8 . 5 )

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

342

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

№Л. 1

Здесь /(я, х, vt)

— достаточно гладкая периодическая по

t с

частотой v

функция,

— стационарный процесс в

R1 с нулевым математическим ожиданием и корреляцион­

ной

функцией

/5Г(т), ф(я,

х) — гладкая ограниченная

функция.

Вводя, как и в детерминированном случае, переменные Ван дер Поля (г, 0), которые определяются из соотно­ шений

х = г cos (соt + 9), х = —гсо sin (со£ + 0),

и «медленное» время s = е£, уравнение (8.5) можно пере­ писать в виде системы двух уравнений

= Fi («s/е + 0s, vs/г, гг, |?/Е) г.

г

(8.6)

- ^ 7 = ^2 (“ s/e + 6е, vs/г, ге, &/е),

где

(т, vt, г, I) = — -i- [/ (г cos т, — rco sin т, vt) +

+ ср (г cos т, — rco sin т) £] sin т*

F2(т, vt1гД) — — I/ cos т, — rco sin т, vt) +

-J- Ф (г cos т, — rco sin т) £] cos т.

Если корреляционная функция К(т) убывает с ростом т, то правые части уравнений (8.6) удовлетворяют условиям теоремы 7.1. При этом, как легко проверить,

 

т

П т

4- I ф (г cos (cos+0),

г-~

о

т— rco sin (cos+0)) £ssin (сОЗ+ 0)<& —

=lim ^ ip(r cos (cos + 0),

— rco sin (cos + 0)) I, cos (cos + 0) ds = 0,

примерь!

 

343

так что слагаемые, содержащие множители

в выраже­

ниях для функций Fly F2 <,

после

усреднения

исчезнут,

н усредненные уравнения

будут

иметь такой же вид,

как при отсутствии стохастических членов. Если отноше­ ние co/v иррационально (этот случай называется нерезо­ нансным), то усредненная система имеет вид

где функции Fx(r), F2(r) задаются формулами (1.9). Предположим теперь дополнительно, что процесс £*

имеет конечные моменты до седьмого порядка включи­ тельно и удовлетворяет условию сильного перемешива­ ния с достаточно быстро убывающим коэффициентом а(т). Тогда в силу^еоремы 3.1 семейство случайных про­ цессов г~~х/2(г? — ru 0f — t е [О, Г], слабо сходит­ ся к марковскому гауссовскому процессу (р*, Ф*), который соответствует дифференциальному оператору

Lu(Pl ^ ) = 4 ( >

где функции А11 (г) определяются равенствами

 

2зх

oo

X cp (г cos s,

— rco sin s) sin s •sin tt

 

 

 

 

 

A*2(r) =

j dt I

dsK ( 'Ч

г ') ф (r cos ~ ra sin

x

 

0

—oo

'

'

 

 

 

 

 

X cp (r cos — rco sin s) cos t

cos s,

 

2 л

00

 

 

 

 

 

 

 

Xq)(rcoss1

— rco sin s) cos t sin $.

12*

ш

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

[ГЛ. 7

 

Как указывалось в § 1, если г0 — единственный корень

уравнения Fx(r) = 0 и функция Fx(r) при переходе через г0 меняет знак с плюса на минус, то, каковы бы ни были начальные условия, при достаточно малых е в системе, описываемой уравнением (8.5) без случайных возмущений (ср = 0), установятся периодические колебания с амплиту­ дой, близкой к г0, и частотой, близкой к со. Если случай­ ные возмущения присутствуют, то, как вытекает из теоре­

мы 7.1, с вероятностью, стремящейся к 1 при е

0, фа­

зовые траектории (X®, Xf) за время порядка е -1 прибли­

зятся к предельному

циклу вида

ГГо =

{(х,

х): х =

= r0cos (со£ “Ь 0), х =

—r0co sin(a)/ +

9)}, и

будут

совер­

шать движение, близкое к колебаниям невозмущенной

системы вдоль цикла ГГо,

отклоняясь время

от времени

от этих колебаний и снова

к ним возвращаясь.

Предполо­

жим, что в начальный момент система уже находилась на цикле ГГо. Тогда на отрезках времени порядка е -1 откло­ нения от периодических колебаний с большой вероят­

ностью будут в

силу теоремы 3.1

иметь порядок ]/е .

С помощью этой

теоремы можно

вычислить различные

вероятностные характеристики таких отклонений. Напри­ мер, чтобы найти вероятность того, что за время [0, Т/г|

амплитуда г® хотя бы раз отклонится от г0 больше чем на h ]/е Анужно решить краевую задачу

- Т АП (го)

 

(Го) г

— h < r < h -,

и (0* г) = 02

и ($2

— h) ~ и (s2 К) =

1.

Тогда искомая вероятность

будет при

е < 1

близка к

и(Т, 0).

 

 

 

 

Отклонения на расстояние порядка 1 происходят за времена, большие по порядку чем е -1. Чтобы изучить эти отклонения, нужно воспользоваться результатами § 4. Пусть, например, l t = при [а - яАг/со, а + я/с/оо), где Цо, тц, . . ., г)А, . . одинаково распределенные не­ зависимые случайные величины; а — случайная величина, не зависящая от {ц*} и распределенная равномерно на

§ В]

ПРИМЕРЫ

345

[О, 2л/со). Тогда легко проверить, что условие F для системы (8.6) выполнено, и

Я(г, 0, <»!, а2) = Я(г, ах, а2) =

где

-= Fx(r)<*i +

+

# ri(ai<Pi(r) + «^Ps(r))^

 

 

 

 

 

 

 

 

 

 

 

Фг (г)

=

- -

2зТаГJ

Ф(r cos

s>

ro)s*n 5) s^n 5 ^

 

 

 

 

 

 

Фг (Г) ^

 

~2 л г ф

J Ф (r cos 5>

гсо sin 5) cos 5 d s 9

 

 

 

 

О

 

 

Ял(т) =

1пМех\

 

 

 

Функции

Fj(r)

и F2(r) определяются равенствами (1.9).

С помощью функции Я(г, a t, a2) можно, как показано в §§ 4 и 6, вычислить асимптотику различных вероятностных характеристик больших уклонений от невозмущенного

движения.

Пусть (а, 6 ) э г 0 и

т8 =

min

{*: rf ф (а, Ь)\.

Вычислим

Пте1пМ г,ет0 при

г е

(а, Ь).

Из результат

е i О

тов § 6 вытекает, что этот предел равен min (и(а), ЦЬ)), где функция и(г) на полупрямой г ^ 0 есть квазипотен­ циал случайных возмущений. Ее можно найти как реше­ ние задачи R ro для уравнения

в (г .

о ) = ? . М 7 ? + Я ,

и £ ) - 0 . (8 .7 )

Решение этой задачи сводится, очевидно, к отысканию ненулевого корня уравнения Fx(r)z + Я Т1(ф1(г)г) =* 0 и последующему интегрированию. Пусть, например, ве­ личины r\k имеют гауссово распределение, Мщ = 0, T>\]h =* о2. Результаты § 4 в этом случае применимы (см.

Г р и н ь [3]). Тогда Я „(т ) = - ~ - ,

z ( r ) = : - .. 2.р (г) И

*

a2<pf (г)

Г

 

Го

346

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

1ГЛ. Т

Предположим теперь, что уравнение Fx(r) = 0 имеет несколько корней г0 < гх < . . г2п, причем в корнях

с четными индексами функция F^r) меняет знак с плюса на минус, в точках r2h + г — с минуса на плюс. Тогда,

если го е (r2/t_1, г2*+1), то при малых е в системе без случайных возмущений установятся колебания с амплиту­ дой, близкой к г2к. Случайные возмущения приведут, вообще говоря, к переходам между устойчивыми предель­

ными циклами.

Пусть

u2h(r)

— решение

задачи i?r2ft

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

(8.7)

на

отрезке

[r2h _ i,

r2ft+1],

Щk(r) < оо и u2k{r2k + г) < и2к(г2к_г). Тогда

с

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

близкой к 1 при малых

е, с цикла ГГ2А =

{(г, 0):

г = г2к}

произойдет переход на

цикл

ГГ2(Л+1),

причем

среднее

время, необходимое для этого перехода, логарифмически

эквивалентно

exp

на

(0,

оо)

равенствамп

Определим функцию F(r)

У 00 = “ го ОО

ПРИ Г е (°> rll;

Г (г) »

V (raft_x)+ Wr2h(r) —

 

— ur2k(r2k- 1)

при

re=

[r^ .j, r2ft+1].

Функция V(r) имеет локальные минимумы в точках г0, г2, . . ., г2п. Предположим, что абсолютного минимума функции У(г) достигает в единственной точке г2^*. Тогда,

как вытекает из результатов гл. 6, предельный цикл ГГ2к* является «самым устойчивым»: подавляющее при

малых е время траектории (г®, 0*) будут проводить в окрестности цикла ГГ2Л*.

4. Рассмотрим линейную систему

Х? = Л(Ь/е) ^ + 6(Ь/.).

(8.8)

Элементы матрицы А (у) = (Л) (у)) и координаты вектора Ъ(у) = (Ь*(у), . . br(i/)) предполагаются ограниченными. Относительно процесса предположим, ■что он имеет достаточно быстро убывающий коэффициент перемешива­

ния; МЛ (£t) = Л = (Л)), Мб(^) = 6. На основании тео­

ремы 2.1 заключаем, что Х\ сходится по вероятности при е \0 равномерно на интервале 0 ^ t ^ Т к решению

§ 8]

ПРИМЕРЫ

347

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

 

 

xt = Axt + 6, х0=

xf

которое можно записать в виде

 

 

t

_

xt =

\exp (Л (t — .5)} 6 ds-f eAiz-

 

о

 

Чтобы оценить нормальные отклонения X\ от хи нужно воспользоваться теоремой 3.1. Предположим для просто­ ты, что процесс h — стационарный и Ь(у) = 0. Обозна­

чим К)п%(т) взаимную корреляционную функцию процес­ сов Л] &) и Anm(tt)y

 

 

 

ео

 

 

 

 

 

 

 

К% =

f irjm (т) dr,

G'n (х) = 2

 

V ".

 

 

 

 

-оо

 

 

i.™

 

 

Тогда

по

теореме

3.1

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

разность

£? =

= е”"1/2 (X® — xt) слабо

сходится при 8 - > 0 к гауссовско­

му процессу

| ^(<~ s)Adrjs,

где г|5 —

гауссовский

про-

цесс с

 

 

0

 

приращениями

и

нулевым Сред-

независимыми

 

 

 

 

 

 

't

 

 

 

ним, для

которого

Мг]*1!? =

|Gin (х$) ds

(см. X

а с ь-

м и н с к и й

[4]).

 

 

о

 

 

 

 

 

 

 

 

 

Нетрудно привести примеры, показывающие, что сис­ тема, получающаяся из (8.8) после усреднения, может иметь асимптотически устойчивые точки покоя там, где векторные поля А(у)х + Ь(у) ни при каких значениях параметра у не имеют положений равновесия или имеют неустойчивые положения равновесия. В окрестности та­

ких точек процесс X®при в 1 проводит большое время или даже с большой вероятностью притягивается к ним. Пусть для определенности, — марковский процесс с конечным числом состояний, рассматриваемый в теоре­ ме 4.2; Ъ= 0, и все собственные значения матрицы А имеют отрицательные вещественные части. Тогда начало

$48

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

СИСТЕМЫ

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

1ГЛ. 7

координат

0 — асимптотически

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

равновесия

усредненной

системы.

Пусть D — ограничен­

ная область с границей dD, содержащая начало коорди­

нат,

и

V(x, dD) — функция,

введенная в

§ 6.

Если

F(0,

dD) <

со, то в соответствии с теоремой 6.1

траекто­

рии

X?, XI = х ^ D,

выходят

из D за время

т8,

стре­

мящееся

по

вероятности к

бесконечности

при

е -> 0;

In Млт8 ~

е -^ О , dD).

Случай

F(0, dD) =

+ оо

рас­

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

— процесс с двумя состоя­

ниями, и матрица А (у) в этих состояниях обращается соот­ ветственно в

Рассмотрим систему Xf = А (£*/е) X?. Если бы процесс It не переходил из одного состояния в другое, то при любом начальном состоянии £0 начало координат было бы неустойчивой точкой покоя — седлом. Пусть матрица Q для процесса \t имеет вид

Стационарное распределение такого процесса — равно­ мерное — (1/2, 1/2). Усредненная система

х1= у а) х11 х2= —•а) х%

при а > а имеет в точке 0 асимптотически устойчивое положение равновесия. Нетрудно доказать, что в этом

случае траектории процесса Xjr при любом е > 0 с вероят­ ностью 1 при оо стремятся к 0, т. е. за счет случайных переходов процесса £* система приобретает устойчивость. С помощью результатов настоящей главы можно вычис­

лить логарифмическую асимптотику

Р ^ т8

< оо}

при

е-^-0,

где T8« m i n { £ :

X f ^ D ] ,

(D — некоторая

ок­

рестность положения

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

а также

асимптоти­

ку этой

вероятности

при

х ->• 0, е s

const.

 

 

$ 0>

УСРЕДНЕНИЕ В СТОХАСТИЧЕСКИХ УРАВНЕНИЯХ

849

§ 9. Принцип усреднения для стохастических дифференциальных уравнений

Рассмотрим систему дифференциальных уравнений

Xе =

Ь(Xе, Уе) + а (X е, Уе) W,

XI =

Xt

Уе =

е

( X е, Уе) + г~шС (X е, Ке) ^ У^ =

у,

где г е / Г ,

У е= Rl,

Ь(х,-

у) *=

(Ъ\х,

у),

Ьг(*, у)),

В(х, у) =

(/?*(#, у), ...,

/?*(#,

г/)),

гг;* —

/г-мерный

винеров-

ский процесс; а (х, у) — (а- (ж, у)), С {х, у) = (cj (х, у))

матрицы, переводящие Вп в Лг и R1соответственно. Функции Ь{(х,у),В{{х, у), а] (о:, у), Cj(xt у) предполагаются огра­ ниченными и удовлетворяющими условию Липшица. На этом примере мы рассмотрим уравнения типа (1.5), когда скорость быстрого движения зависит от медленных пере­ менных. Заметим также, что, в отличие от предыдущих параграфов, медленные переменные в (9.1) даже при фик­

сированном

Yfl

t е [0А

t], образуют

случайный

процесс.

 

 

 

 

 

 

Введем в

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

случайный

процесс

У*у,

х е i f , у е

R1,

который

определяется

стохастическим

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

 

 

 

Y?v = В (х, YГ) + С (х, Y f) wu

JT =

У-

(9.2)

Решения этого уравнения образуют марковский процесс в

R\ зависящий от х е i f , как от параметра.

Сначала мы сформулируем и докажем принцип усред­ нения в случае, когда элементы матрицы а(х, у) не зависят от г/, а затем укажем на те изменения, которые необходи­

мо сделать, чтобы рассмотреть общий случай.

 

Предположим, что существует

функция Ъ{х) =

(Ьг(х),

Щх), .... 1r(x)),

are Rr, такая*

что при любых

t ^ О*

X <=

R\ у <= R1

 

 

 

 

 

1+Т

 

 

 

 

М

Y

J

b (х>Y*v) ^ ~

Ь(х)

(9.3)

где

х(Г)->-0

при

Г -*• оо.

 

 

350

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

[ГЛ. 7

=

Т е о р е м а 9.1. Пусть элементы матрицы а(.г,

у)

а(.г) не зависят от у и выполнено условие (9.3). Обозна­

чим X t случайный процесс в Rrt определяемый дифферен­ циальным уравнением *)

Xt=b (Xt)+ а (Xt) wt, Х0 = X.

Тогда для произвольных

Т > 0,

б

0, х е i?r, у е R1

l i m p( sup

\X*t - X

t\>

б) = 0 .

е->0

 

 

>

Д о к а з а т е л ь с т в о . Рассмотрим разбиение отрез­ ка [0, Т] на одинаковые интервалы длины Д. Построим

вспомогательные процессы Yt и Xt с помощью следую­ щих соотношений:

г

И = П д

+ - ^ В{Х1&

Yt) ds + - L

f С (Хеш

Yt) dw.

 

 

feA

 

 

Уе L

 

 

 

 

 

 

 

 

< е [ Ц

(Л+1)Д],

 

 

 

 

t

 

 

 

 

%

= * +

I b (XfsMU) f t ) ds + S O (Xt) dwt.

 

 

 

0

 

0

 

 

 

 

что

Покажем, что интервалы Д =

Д(е) можно выбрать так,

е_1Д(е) -> оо2

Д(е)

0 при

е -^ 0

и

 

 

 

 

М |К ?-У ?|2-+ 0

 

(9.4)

равномерно по х е

К ,

у <= Яг,

i s

(0, Т]. Из определе­

ния процессов Yt и Yf

следует, что прп t, принадлежа­

 

*) В

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

wf — тот же

самый

винеровский про­

цесс, что и в (9.1). Так как функция Ь(х, у) удовлетворяет условию

Липшица, то Ь(х) тоже удовлетворяет условию Липшица, так что решение уравнения существует и единственно.