![](/user_photo/_userpic.png)
книги / Флуктуации в динамических системах под действием малых случайных возмущений
..pdfi 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 s£ — 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я |
|
|
|
|
|
|
|
|
|
|
Фг (г) |
= |
- - |
2зТаГJ |
Ф(r cos |
s> |
ro)s*n 5) s^n 5 ^ |
|
|
|
|
2л |
|
|
Фг (Г) ^ |
|
~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). Так как функция Ь(х, у) удовлетворяет условию
Липшица, то Ь(х) тоже удовлетворяет условию Липшица, так что решение уравнения существует и единственно.