![](/user_photo/_userpic.png)
книги / Флуктуации в динамических системах под действием малых случайных возмущений
..pdf* 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, что для любой огра ниченной измеримой функции на окружности /(ф)
|
|
|
2Л |
|
|
К о / ($tr,)) — j / (ф) т (г0, ф) dф < С е м |
sup |/ (ф) | |
||||
|
|
|
|
|
0<Ф^2л |
(см., например, Ф р е й д л и н |
[3 ]). Отсюда вытекает соот- |
||||
|
|
|
|
|
2л |
ношение |
(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, Р) полунепрерывна снизу.