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

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

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

* 9]

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

351

щем отрезку [&Д, +

1)

Д]

 

 

М IYt -

У<| =

 

 

 

 

 

 

 

t

 

 

 

 

 

 

ад

 

 

 

 

 

t

 

 

 

 

+ 7 = f

[С {XI

п ) - с (Л1д, Щ dw,

 

< ^ - 1 ( 1

М |X® — Л’аД|<fs +

J Ml v! y f |ds ) +

 

'АД

 

 

АД

/

 

/ t

 

 

t

 

4

4— “

( j* M |A- — АЕД

+ J

M f Y$ — YE\~ds

 

 

'АД

 

t

kA

 

 

 

^

(p ' ^ "e) J M I

— А^д |d$ +

 

 

 

 

АД

 

 

 

 

 

 

t

 

 

 

 

+ C4 (|L +

i ) J MI t f - Yes 12ds.

(9.5)

Здесь и в дальнейшем мы обозначаем через С* постоянные, зависящие только от константы Липшица коэффициентов Ъ{{х, у), В'(х, у), а) (х, I/), Cj (.г, у), максимума модулей этих коэффициентов и от размерности пространства.

Из ограниченности коэффициентов стохастического

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

А® вытекает при Д <

1 оценка

 

 

 

М |X? -

А ЕД| < Сь|« -

ЛД |

(9.6)

для s s

[АгЛ, (А:-}-1)Д].

Из этого неравенства и (9.5)

при

t е [/сД,

+

1)Д] получаем

 

 

352

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

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

ГГЛ. 7

откуда

приходим к соотношению

 

М\Y\ — Y\|*

+ |f) Д ехр [С6 ( § + 4)'- .

Отсюда заключаем, что (9.4) выполняется, если положить

д =

д(е) = t V In в-1 -

 

 

 

 

 

 

 

Покажем теперь,

что

для

любого б >

О

 

 

р ( sup

| Х ? - Я ? | > в \ — О

(9.7)

при е - ^ 0 и Д = Д(е) =

е / . а

, -

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

/?г,

у е

R1. Действительно,

из определения

процессов

X*

и Х\ следует, что

 

 

 

 

 

 

 

P f

sup |Х<- X f | >

б ) <

 

 

 

 

 

\о«<г

J

 

 

 

 

 

 

 

< Р |j

|ь {XI, У*) -

Ь( Х ^ д ,

У ,) I ds >

6

Оценивая вероятность, стоящую в правой части,с помощью неравенства Чебышёва и принимая во внимание (9.4) и (9.6), получаем (9.7). Из оценок (9.4) и (9.6) легко полу­ чить также, что

sup М| Х ? -£ ? | 2-* 0

(9.8)

0<*<Т

 

при е -*• 0.

 

Теперь убедимся в том, что sup |X® —

| стремится

о<s<T

 

к нулю по вероятности вместе с 8. Отсюда и из (9.7) будет следовать, очевидно, утверждение теоремы.

Прежде всего заметим, что из определения процесса

У® вытекает,

что процесс Z8=

У*д4-« ПРИ

$ ^ [0, Д ]

совпадает по

распределению о

процессом

Yf£*Yh*1

определяемым уравнением (9.2). Следует только в (9.2) выбрать винеровский процесс wt не зависящим от X\^t Y\

§ 91

 

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

353

Отсюда,

учитываяг что е-1Д(е) -> оог на основании

(9.3)

получаем

 

 

(ft+l)AГА Л*

 

 

 

М

J

Ъ(Х®д, Y?) dt Д6 (Х®д)

 

 

I

 

 

 

 

 

АД

 

 

 

 

 

 

 

 

Д /8

 

 

=

ДМ i

j [ H x ;Atz , ) - b ( x j A)]ds <Д -х(Д /е).

Используя

эту

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

 

М j sup

\b (Xf8/A]A,

1<К«т

i

1

 

Ml

max

f t ) ds — [ b (X 8e) ds

I

 

0

J

|

(M-i)A

 

2

J ■[ b C x L t f ! ) -

 

fcA

 

 

 

 

-

b{XtA)] ds

+

с 7д <

 

 

<

[Т/Д1

м

(HГ1/<-l)A-»

 

 

 

 

S

J

[Ь (хгд, У8е) -

Ц х ^ )1 ds

+ С7Д <

 

k=0

 

ЙД

<

С7Д + 7и (Д/е)

0

(9.9)

 

 

 

 

при

е -> 0 ,

так как Д(е) -> 0^ х(Д(е)/е)->0

при е

0.

О цент

тг(t) = М|Х*— Xt |* Из определения

про­

цессов Xt и Х\

вытекает

 

 

 

 

 

 

t

 

b (xOl ds +

 

 

Xt -

X, =

j

[ь (xf,M1Al ft ) -

 

 

 

 

0

 

 

 

 

 

 

 

 

t

 

 

 

t

 

 

 

+

J

[b{xt) - b (X ,)] ds +

J [a (X 8e) -

a (X .)] dip..

 

 

6

 

 

 

0

 

 

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

т

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

[Ел. г

придем

к соотношению

 

 

 

i

г

 

те(t) ^

Cst J те (s) ds +

С9J те (s) ds +

 

 

+ ЗМ

 

 

 

 

\b (X l)d s

 

Из этого

соотношения

при

i e

[О,

Т] получаем

 

 

t

 

 

t

 

eCJ0(T+T*)#

 

тп (f) < ЗМ J* Ь (Xfs/A]A, Y s) d s - fF (X 3e)] ds

 

 

0

 

 

0

 

 

 

Отсюда в силу (9.9) вытекает,

что

m?(t) -*■ 0 при е

О

равномерно на отрезке [О, Г].

 

неравенство

Далее1

для любого

б >

0

справедливо

<^Р{

sup

’ [ь(А1,/д]л, У 5 ) - Ь М ] ds

> 5/ej +

|о<^г

 

 

 

( Т ____

 

\

+

 

 

+ PH \b(Xs) - b (X * )\ d s > №

J

 

 

 

 

16

 

 

 

 

 

+

р

|j

б|(Х 8е) -

I (Х|) |ds > б,б| +

 

+

Р {

sup

( [ст (X f) — а (Х 8)] dws

>

8/61

+

 

I <к«т

 

 

 

I

 

 

 

 

+ Р

sup

J [o (X 8e) - o ( X

8)] dwa > 8/6 .

Первое слагаемое в правой части этого неравенства стре­ мится к нулю в силу (9.9). Чтобы показать стремление к нулю второго и третьего слагаемых, нужно воспользовать­ ся неравенством Чебышева, соотношением (9.8) и тем, что me(t) -> 0 при е | 0. Четвертое и пятое слагаемые оцени­ ваются с помощью неравенства Колмогорова и тоже

4 91

УСРЕДНЕЙЙЕ В (УГОХАСТЙЧЕСКЙХ УЁАВНЁНЙЯХ

355

стремятся к нулю.

Таким образом, мы получаем,

что

s u p J x f - X j - ^ O

при е\ 0 по вероятности. Отсюда и

пз (9.7) вытекает утверждение теоремы 9.1.

 

 

Скажем теперь несколько слов относительно случая,

когда элементы матрицы о зависят от а: и от у. Предполо­ жим, что выполнено условие (9.3) и, кроме того, существу­

ет

матрица

а(х) =

{а^(х))

такая2 что при любых t ^ О*

х е

Rr> у е

R1

 

 

 

 

 

 

 

t+T

 

 

 

 

maxM

|ji

(

2

о* (*i

al (х, Y*v)ds — а^{х) |<х(7т)1

i.i

 

t

h

 

 

 

 

 

 

 

 

 

 

 

 

(9.10)

где

х(Г) ->■ 0

при

Т —►оо, Ygy

решение

уравнения

(9.2).

Пусть

X t — решение стохастического

дифферен­

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

 

 

 

 

 

 

%t = b(X t) +

o(Xt)w t,

Х0 = xi

 

где а(х) = (a(x))lf2, wt —■г-мернып винеровский процесс. Оказывается, что и в этом случае будет сходимость X® к Xt.

Однако эту сходимость нужно понимать в несколько более слабом смысле, а именно, меры в пространстве

траекторий, соответствующие процессам Xf, будут при г | 0 слабо сходиться к мере, соответствующей процессу

X*. Мы не будем здесь приводить доказательство этого

утверждения,

отсылая

читателя к

работе (X а с ь м и н-

с к и й [6]),

где оно

содержится.

Отметим, что в этой

работе принцип усреднения доказан в предположениях, допускающих некоторый рост коэффициентов. Условия (9.3) и (9.10) также заменены менее жесткими.

Рассмотрим пример. Пусть (ге, фр) — двумерный мар­ ковский процесс, управляемый дифференциальным опера­ тором

L * = 1 • £ + b (Г1 Ч>) Т г + т [ В <Г1 Ч>) - Щ + Т С2

Ч>)

356

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

ИГЛ. 7

Этот процесс можно описать также с помощью стохасти­ ческих уравнений

г8 = Ь(г*3 <ре) +

ф8 = ъ-1В{г*1 фе) -J-

фe)w2.

Предположим, что функции Ь(г, ф), В(г, ф), С(г, ф) перио­ дические по ф с периодом 2л; С(г, ф) ^ с0 > 0. В этом

случае процесс ф(*Го)2

который получается из процесса

t

г

ф((го> = ф(оГо) + j в (г0, ф(; о)) ds + J с(г0, ф(; о)) dw,

о

о

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

 

 

 

 

 

К о / ($tr,)) — j / (ф) т (г0, ф) dф < С е м

sup |/ (ф) |

 

 

 

 

 

0<Ф^2л

(см., например, Ф р е й д л и н

[3 ]). Отсюда вытекает соот-

 

 

 

 

 

ношение

(9.3). Действительно^

положив b{r) = J Ь(гг ф)х

X т(г,

cp)dcp, получим

 

о

 

 

 

 

л

 

 

 

М,

Y

J

ь(г, Ф(3Г)) ds — b (г)

 

 

'фо

 

 

 

 

 

 

< М ,фо Y 'J &(г1 Ф(/0 ds — b (г)

 

 

т

т

 

 

 

 

j

j Мф. (Ь (г1ф!г)) - Ь(г)) (б (г, ф!г)) -

6 (г)) ds

§ 9]

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

357

ТТ

«= jj j* dsJ мф. (b(ri ф*г)) - b(г))мф3 (г-ф»-.) —

Об

т т

при Т оо.

Таким образом, в силу теоремы 9.1, г* равномерно на отрезке 0 ^ i ^ Т сходится по вероятности при е \ О

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

 

Если имеется некоторый непрерывный в пространстве

Сот{йг)

функционал F [rj,

0

Т]

от

процесса г®,

то

из сказанного следует,

что

F [г®, 0 ^

t ^

71]

0 ^

£

Т] по вероятности при е \ 0.

Эта сходимость

сохранится, если функционал F имеет разрывы, но мно­ жество функций, на которых F разрывен, имеет меру нуль относительно меры в пространстве функций индуциро­

ванной предельным процессом rt.

Результаты, приведенные в этом параграфе, можно применить для исследования поведения при е | 0 реше­ ний некоторых эллиптических и параболических уравне­

ний с малым

параметром.

задачу Дирихле

Рассмотрим,

например,

 

Ь*иЕ(гг <р)

=

0,

г <= (rv

г ,);

u8(ri,

ср) =

Cv

и*(г2, ср) =

(9.11)

С2

в области D = {(г,

ср):

 

0 <

гх < г <

г2}. Если (г, ср)

интерпретировать как полярные координаты на плоскости, то эта область представляет собой кольцо, ограниченное двумя концентрическими окружностями радиусов гг и г2 соответственно с центром в начале координат. Как известно2 решение этой задачи можно записать в виде

(Гг ф) “ С]Рг,чр ( г \е ■*= Fj.)

C JP г,<f { г те = Гг ),

358 ДИНАМИЧЕСКИЕ СИСТЕМЫ С УСРЕДНЕНИЕМ [ГЛ. 7

где те =

inf {^:

rf ф (гь г2)}.

Обозначим

т =

inf {£ >

0:

rt ф (гА,

г2)}. Легко проверить,

что max

Рг{ т > Г } - > 0

при

 

оо

и

что

граничные

 

г,<г<г,

 

[rlf

г2]

 

точки

отрезка

регулярны

для

процесса

rt

в

[г*,

г2],

т.

е.

что

Рг|{т = 0} =

1,

i =

1, 2

(см.

В е н т ц

е л ь [1 ]).

Отсюда

и из равномерной

сходимости

по

вероятности

т\ к rt

на

любом

конечном

отрезке

[0,

 

Т]

вытекает,

что

П т

и¥(г}

ф) =

u(r)

= С1Р г(гТ =

 

 

 

 

 

Функ-

е 10

_

 

 

 

 

 

 

 

 

 

 

 

 

 

ция и(г)

может быть пайдена как решение задачи

 

 

 

-j- и" (г) -J- Ь (г) и1(г) — 0,

г е ( г „ гг);

 

 

 

 

 

 

и (^i) = ^1» и (^2) =

Сг*

 

 

 

 

Решая эту задачу* окончательно получаем

Пт иг (гх ф) = и (г) =* е;о

Некоторые примеры более общего характера можно найтв

в работе Х а с ь м и н с к о г о

[6].

Рассмотрим теперь большие уклонения в системах типа

(9.1). Мы ограничимся случаем,

когда в медленном движе­

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

менных.

пусть Л/, Е —- два

римановых

многообразия

Итак,

класса С00, многообразие Е

компактно,

dim

М = г,

dim Е =

I. Обозначим ТМх и ТЕУкасательные простран­

ства к этим многообразиям в точках х е М , у е

Е соот­

ветственно. Рассмотрим семейство векторных полей Ь(х, у) на Л/, зависящих от у е Е как от параметра, и семейство В(х, у) полей па зависящих от х е М. Кроме того, на многообразии Е рассмотрим эллиптический дифферен­ циальный оператор L второго порядка* переводящий

§ 9 ) Усреднение б стохастических уравнениях 359

константы в нуль. Функции Ъ(х, у), В ( х , у ), так же как коэффициенты оператора L, предполагаются бесконечно

дифференцируемыми

по совокупности

переменных.

На прямом произведении М X Е рассмотрим семейство

марковских процессов

Z\ = (X Y])t

управляемых опе­

раторами

 

 

У) = {Ь(х, у), ? J(x , у)) +

+ B-1lLgf(x, у) 4-(В(х, у), V yf(x, y))\t

где V x, V у — операторы градиента на М и Е соответст­

венно. В координатах траектории процесса Ъ\ = (X ft Yf) можно задать системой стохастических уравнений

л f =

г>(xf, yf),

 

(9.12)

У? = ® * [fi (Xf, у О + 8 (у?)] + г~ т С (yf)

ГДе

= (gl(l/)i ...» #4*/)) — коэффициенты при первых

производных в операторе L, матрица С(у) связана с коэф­ фициентами aij(y) при старших производных в операторе L соотношением С(у)С*(у) = {а^(у))\ wt — /-мерный винеровский процесс.

Если многообразие М не компактно, то нужно еще сделать предположение о том, что поле Ъ{х, у) не слишкОхМ

быстро растет: для любого Т >

0 и х е

М можно указать

компакт

F czM ,

x ^ F , такой,

что

Рxy\X t^F при

/ е [О,

Т]} = 1

для любых

у е

Е,

г > 0.

Пусть а — элемент пространства Т*МХ, сопряженного к ТМХ. Введем в рассхмотрение дифференциальный опера­ тор R = R{x, z2а), действующий на функции j(y)x у ^ Ег по формуле

R{x, z, а)/(у) =

= £/(*/) + (^(2 , y),Vyf(y)) + (а, Ь(х, y))f(y)\,х

х, z е М, а е — параметры. При произвольных значениях параметров Д(я, z, а) — эллиптический диф­ ференциальный оператор в пространстве функций, опре­ деленных на кохмпакте Е. Этот оператор является сужени-

360

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

[ГЛ. 7

ем на гладкие функции производящего оператора некото­ рой положительной полугруппы. Аналогично тому, как это делалось в § 4, отсюда можно вывести, что опера­ тор R(xy z, а) имеет однократное собственное значение ц(.г, z, а) с наибольшей вещественной частью. Это собст­ венное значение вещественно и в силу однократности диф­ ференцируемо по параметрам х, z, а.

Введем в

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

диффузионный

процесс Y\y

z е М, на

многообразии

Е% управляемый

оператором

 

N * = L +

(B(z, у), Vy).

 

 

Л е м м а

9.1. Пусть а е Т*МхУ и

F некоторый

компакт в М. Равномерно по х, z e F , у е

Е

существует

предел

 

 

 

 

 

т

 

 

 

lim -=-1пМу exp (а, Ъ(х, У")) dsj =

\х (х, z, а).

Г-.оо 1

 

 

 

 

Функция ц(.г, z, а) выпукла вниз по переменным а.

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

Обозначим

V(xy уу а) =

=* (а, Ь(ху у)). Тогда семейство операторов

действую­

щих по формуле

 

 

 

T jf

(у) = Му/ (Я ) exp j j У (х, Yzst а) d*j

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

Обозначим L(xy zy Р) (ху z <= Му $ ^

ТМХ) преобразо­

вание Лежандра функции

z1

а)

по последнему аргу­

менту:

 

 

 

 

L (xi zi Р) = SUP Ц а 3

Р)

(« 1

2 , а ) ] .

а

 

 

 

 

Мы будем иногда рассматривать функцию L(xy z, Р) при совпадающих первых двух аргументах. Будем обозначать Ь(ху х, Р) = L(xy Р); эта функция, очевидно, является преобразованием Лежандра от р,(;г, ху а). Заметим,, что функция Ь(х± z, Р) полунепрерывна снизу.