![](/user_photo/_userpic.png)
книги / Флуктуации в динамических системах под действием малых случайных возмущений
..pdfi п |
НЕ ОЧЕНЬ БОЛЬШИЕ УКЛОНЕНИЯ |
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 ®и дру