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

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

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

§ И

ЛОКАЛЬНЫЕ ТЕОРЕМЫ, ТОЧНАЯ АСИМПТОТИКА

391

где

(ciij(x)) = (а^(.г)) 1. Устанавливается, что

 

 

lim h In ph(t, х, у) = — V (t, х, у).

(1.4)

 

hi О

 

Эта грубая локальная теорема о больших уклонениях по­

лучается применением грубых теорем из

работы

В е н т-

ц е л

я, Ф р е й д л и н а

[4] (интегральных) и оценок из

работы А р о н с о н а

[1 ]:

 

 

 

 

А„

 

'

cn\y— x(t,

а Л I 2 )

 

 

Рь (*. х. У) < (hty'Z cxpj

 

Ы

 

j

 

Ph(*,

А,

exp |f

ci\lf — x(t,x) |

 

 

У )> (ht)Tli

 

l

ht

 

c2\(y — x(t,

 

 

 

 

( / ■ о

" " -

M p

Ы ' ]

для достаточно малых t, где A h ct и а — положительные

константы, a x(t, х)

— решение системы х = Ъ(х) с началь­

ным условием #(0, х) — х.

 

Далее,

если

(Xf, P j)— диффузионный

процесс,

соответствующий

оператору Lh в области

D с глад­

кой границей 3D, исчезающий по достижении этой гра­

ницы, то его плотность вероятностей

перехода qh(t, х,

у) —■функция Грина для уравнения

^ = Lhu

с гра­

ничным условием и = 0 на 3D. Доказывается, что

lim h In qk (t, x, y) = — VD (t, x, y)t

(1.5)

hi 0

 

 

где VD(J, х1 у) определяется как нижняя грань значений функционала (ф) на кривых ср, соединяющих точки х и у за время г, не выходя из области D. Отсюда, в частности, вытекает, что

Ип1 l ib .. *<.Vl = о,

(1.

h i O p n (t, х , г / )

 

если все экстремали, на которых

достигается

мини­

мум (1.2), выходят из D U 3D.

Доказывается

также

противоположный результат: если все экстремали прохо-

.92 УТОЧНЕНИЯ И ОБОБЩЕНИЯ 1ГЛ. 9

дят

внутри D,

то

 

 

 

 

 

 

 

 

 

 

= 1.

(1.7)

 

 

h i О р

(t,

X,

у)

 

 

(По

выражению

М. Каца,

процесс X

при малых

h

«не

ощущает» границы.)

схему

естественным образом

В рассматриваемую

почти в точности) включаются также вопросы, касающие­ ся асимптотики переходной плотности диффузионного процесса (фундаментального решения параболического уравнения), не зависящего от параметра, при малых зна­ чениях времени. Действительно, пусть p(t, х, у) — плот­

ность

вероятностей перехода

диффузионного

процесса

(X t,

Рх), соответствующего эллиптическому

оператору

 

Lu =

1 V

д~и

,

 

, i / ч д и

 

 

 

(х) ---:---- +

ZJЬ (X) ---

 

 

 

 

д х 'д х *

 

д х

1

 

Рассмотрим семейство

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

процессов

Хъ, =

= Xht, /г > 0 . Переходная плотность

p{h, х, у)

процесса

Х г за время /г

равна

переходной

плотности

рл(1, х, у)

процесса Xht за время 1. Процесс (X*, Р£) управляется диф­

ференциальным

оператором

 

 

 

 

L lu hLu =

aij(x) д2 а

 

 

0-8)

 

 

д х ' д х Т + ' * 2 '> '« £ •

Коэффициенты

переноса

к-Ь{(х)

стремятся

к нулю при

/г | 0, так что семейство операторов

(1.8) — это почти то

же, что семейство (1.1) с заменой Ь1(х)

на нуль. Во всяком

случае, функционал действия

выписывается без

труда:

он имеет вид h~lSQ (ф),

 

 

 

 

 

 

 

t

 

 

 

 

^оНф)

о

 

 

 

(1.9)

 

 

 

 

 

 

Задача о нахождении

минимума

такого

функциона­

ла решается в терминах римановой метрики р(х, у), свя­ занной с матрицей (atj) :

t

р{х, у) = min | [2 Х Д ф8) фЖ ] ,/2* -

Фо— I—У о

§ И

ЛОКАЛЬНЫЕ ТЕОРЕМЫ, ТОЧНАЯ АСИМПТОТИКА

393

Легко проверить, что минимум функционала (1.9) по всем параметризациям <рв, 0 ^ s <1 t, данной кривой равен квадрату римановой длины кривой, умноженному на (2£)~*; а значит, минимум (1.2) равен (2г)—1 р(.г, г/)2.

Применение формулы (1.4) дает

lim h ln p (h, х, у) =

-----^-p(x, i/)3.

(1.10)

hi 0

^

 

Соответственно для переходной плотности диффузионного процесса с исчезновением на границе области D полу­ чается

limbing (/г, х, у) =

— 4 - р D (ж, У?,

(1-11)

hi 0

*

 

где pD (х, у) — нижняя грань римановых длин кривых, соединяющих точки х и у, не выходя из D . Получаются также результаты, касающиеся отношения плотностей q и р при малых значениях временного аргумента, соот­ ветствующие (1.6), (1.7); в частности, «принцип неощути­ тельности границы» принимает такой вид:

lim £lL£iJd

= 1,

( 1. 12)

Н0 Р(*, У)

если все кратчайшие геодезические, соединяющие точки х

п у, лежат целиком в D .

были получены

в статьях

Результаты

(1.10) — (1.12)

В а р а д а н а

[2], [3] ((1.12)

для

частного случая

вине-

ровского процесса — у Ц и е с е л

ь с к о г о

[1]).

[1 ]

В статьях

К и ф е р а [1 ],

[3],

М о л ч а н о в а

были получены точные варианты этих результатов. При их получении существенную роль играют соответствую­ щие грубые результаты,—несущественно, в локальной или интегральной: форме: они дают возможность исключить из рассмотрения всё, за исключением окрестности экстремали (экстремалей) функционала действия; после этого задача становится локальной.

Точная асимптотика плотности вероятностей перехода из точки х в точку у оказывается различной в зависимости от того, являются точки х и у сопряженными вдоль какойлибо из связывающих их экстремалей или нет (см. Г е л ь ­ фанд и Ф о м и н [1]). В случае, когда точки х и у не

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

894

УТОЧНЕНИЯ И ОБОБЩЕНИЯ

1ГЛ. 9

оказываются сопряженными, а коэффициенты оператора гладки, получается не только асимптотика для плотно­ стей с точностью до эквивалентности, но и асимптоти­ ческие разложения по степеням малого параметра: для семейства процессов, соответствующих операторам (1.1),

Ph(t,

х. У) = (2яht)~r/* ехр

{ —

 

х, y)}\KQ(t, х, у) +

 

 

+

hK1{t, х,у)

+ . . .

+h™Km{t, х, y)+o(h™)] (1.13)

при

h |

0 ( К и ф е р

[3]); для процесса с

производящим

оператором

(1.8)

 

 

 

 

 

 

 

P(t,

х,

у)

= (2л£)"г/2 ехр

{ —р(я, y)2/2t}

1К0(х,

у) +

 

 

 

+ tK1(z,

у) + . . .

+ tmKm{x, у) + o{tm)] (1.14)

при

t

|

0 ( М о л ч а н о в

[1]).

Вероятностные методы

в этих работах комбинируются с

аналитическими.

Наме­

тим

доказательство

разложения

(1.13) в самом простом

случае,

притом для т = 0;

т. е. речь идет о доказательст­

ве

существования

конечного

предела

 

 

 

 

lim (2nht)rl2ph (t, х, у) ехр [hT'V (t, ххг/)).

(1.15)

 

 

hi 0

 

 

 

 

 

 

 

Пусть

оператор Lh имеет вид

 

 

 

 

 

 

 

 

Ь Ы -.-^ Ь и + '^ Ь ^ х )^ .

 

(1.16)

Процесс,

соответствующий этому оператору, можно за­

дать прп

помощи

стохастического уравнения

 

 

 

 

 

X j = b ( X j ) +

/*1/2i^,

X h0 = xt

(1.17)

где ws ■— r-мерный винеровский процесс. Рассмотрим на­ ряду с X* неоднородный по времени диффузионный про­ цесс Yg, задаваемый стохастическим уравнением

Yh,=<Va + hU2wn Yh0 = x t

(1.18)

где фв, 0 ^ s ^ экстремаль функционала действия, ведущая из точки х в у за время t (предполагаем., что она

§ И

ЛОКАЛЬНЫЕ ТЕОРЕМЫ, ТОЧНАЯ АСИМПТОТИКА.

395

единственна). Плотность распределения Yht легко

вы­

писывается: это — плотность нормального распределения со средним срt и матрицей ковариаций ЫЕ; в точке у она

равна (2л/^)“ г/2-

 

 

 

Отношение плотностей

вероятностей

Xht и У?

равно

пределу отношения

 

 

 

P[X?£=i/ +

D)/Ptr?<=i/ +

Z)}

(1.19)

при стремлении к нулю диаметра окрестности D начала координат. Воспользуемся тем, что меры в пространстве

Cot(Rr), соответствующие случайным процессам

Xj, Yhs,

абсолютно непрерывны друг* относительно друга;

плот­

ность имеет

вид

 

 

d

(Ул) =

ехр А 1/2 j (iЬ(У*) — <ps, dwa)

 

 

rfPyft

 

 

 

 

 

 

 

 

(2h)~l 11 b (Y a) — cps |ds 1.

(1.20)

 

 

0

J

 

Абсолютная непрерывность позволяет выразить вероят­ ность любого события, связанного с процессом Х л, в виде интеграла от функционала от Yh\ в частности,

Р1Х?ег/ + £} =м{у?е=г/ + Я; ^ (y * )J .

Выражение (1.19) принимает вид

 

 

 

М ( exp |а_ 1/2 j (У?) - ф„ dwa) -

 

 

 

- (2А)-1 J |b (У?) -

% |ds

У? s у +

D

=

= M jexp |A 1/2 J (b (<ps + hU2ws) — <pS) dws)

— (2A)“ 1j |b (cp, +

hl,2ws) -

cps |c?sJ

wt s

h~l,iDJ

1 4 *

396 УТОЧНЕНИЯ И ОБОБЩЕНИЯ 1ГЛ. Я

(пользуемся тем, что w0 = 0, и Y* = фв +

W 2w9 при

О < 5 <

t).

ней Л” 1/2/))

При

стягивании окрестности D (а вместе с

к 0 получаем условное математическое ожидание при ус­ ловии, что wt = 0. Итак,

(;2п% ~ г/*

= М { ехр Г

1/21 ^

+ hUZw^

~ ф*’ dw$) ~

— (2h)~l J |Ъ(cps +

hl,2ws) — <р. |ds

wt = 0 . (1.21)

(Формула

(1.21) выписывалась фактически еще в работе

Х а н т а

[1 ].)

существование предела (1.15) при

Чтобы

установить

h | 0, прежде всего урезаем математическое ожидание (1.21) умножением экспоненты на индикатор события

{шах |Л1/2и>в|< б};

то,

что отброшенной

частью можно

О<8<*

0,

устанавливается с

помощью гру­

пренебречь при h j

бых оценок,

связанных с функционалом действия. Затем

преобразуем

выражение под знаком экспоненты с учетом

гладкости функции Ъ. В первом интеграле берем разложе­

ние Ь(ф3 + Л1/2^ ) по степеням h1/2w3 до членов первого

порядка, а во втором — до членов второго порядка. Глав­

ный член, возникающий из второго

интеграла,

равен

t

 

 

 

 

 

 

—(2A)_1J l % s) — %\2ds,

 

 

 

О

 

 

 

 

 

 

и при подстановке в формулу

(1.15) он сокращается

с

/i"1F(^, я, у). Если мы проверим, что сокращаются друг

с

другом также члены порядка й -1/2, то

чтобы

завершить

доказательство, останется

только установить, что сходи­

мости не помешают члены

порядка 1 и бесконечно малые

при h | 0 (это удается установить в случае,

когда

точки

х и у не сопряжены).

 

fe-1/2,

возникающими

из

Займемся членами порядка

первого интеграла в (1.21);

произведем интегрирование по

§ И

ЛОКАЛЬНЫЕ ТЕОРЕМЫ, ТОЧНАЯ АСИМПТОТИКА

397

частям:

 

 

 

 

 

t

 

 

 

 

h ~ in J (b (ф8) -

<ps, dws) = h~m (b ((ft) -

<p„ wt) -

 

 

0

 

t

 

 

 

 

 

 

 

 

 

— h~ i/2 j

i

Ф (ф.) -

ч>.) )

Виеинтегральный член исчезает в силу условия

wt = 0;

а

выражение

^ Ф(ц>8) —*Ф3) преобразуется с

учетом

уравнения Эйлера для экстремали, и, как легко проверить, интеграл сокращается с членами порядка ft-1/2, возникаю­ щими из второго интеграла в (1.21).

Из намеченного доказательства видно, что коэффициент K0(t, х, у) в разложении (1.13) имеет смысл условного ма­ тематического ожидания экспоненты от квадратичного функционала от винеровского процесса при условии wt = 0.

Положение не сложнее и в случае, когда точки х и у соединены конечным числом экстремалей, причем они не сопряжены вдоль никакой из них. В случае сопряженных точек можно выразить переходную плотность за время t при помощи уравнения Чэпмена — Колмогорова через переходную плотность за меньшее время; в полученном интеграле основную роль будут играть плотности в не­ сопряженных точках, и асимптотика получается приме­ нением метода Лапласа (конечномерного). Для плотности

p(t,

х, у) при t | 0 это проведено в работе М о л ч а н о ­

ва

[1 ]; в частности, при

разной структуре множества ми­

нимальных геодезических, соединяющих х и у,

возни­

кают асимптотические

выражения вида р(£, х,

у) ~

~ C t —a e~~p(‘x'y)*/2t с различными а.

Что касается точной асимптотики в задачах о больших уклонениях, не сводящихся к одномерным или конечно­ мерным распределениям, то здесь пока сделано мало. Полученные результаты относятся не к вероятностям Ph{Xh e i ) и не к плотностям, а к математическим ожи­ даниям вида

М* e x p {ft-Ч Д Х * )} ,

(1.22)

398

Ут о ч н е н и я и о б о б щ е н и я

[ГЛ. 9

где F — гладкий функционал (нормирующий коэффици­ ент предполагается равным Ь,—1). Рассмотрение задач та­ кого рода в качестве первого шага естественно, потому что даже в случае больших уклонений для сумм независимых двумерных случайных векторов точная асимптотика ин­ тегралов, аналогичных (1.22), находится существенно легче, чем точная асимптотика вероятности попадания в область.

Выражение (1.22) логарифмически эквивалентно explfe^max [F — 5]}, где S — нормированный функ­ ционал действия. Если экстремаль ф, доставляющая этот максимум, единственна, математическое ожидание (1.22) отличается от

ММР(*Л.Ф) < б; ехР { й - Ч т } }

(1.23)

на величину, экспоненциально малую по сравнению с (1.22) или (1.23). Это позволяет локализовать задачу.

План дальнейшего исследования в большой степени аналогичен тому, что делалось в статье К р а м е р а [1 ]: производится обобщенное преобразование Крамера, за­

меняющее меру РЛ на новую вероятностную меру РЛ такую, что при малых h «наиболее вероятной» относитель­

но

РЛтраекторией процесса X ht оказывается экстремаль

Ф*.

Случайный процесс fe-1/2 [Х^1 — ф*]

относительно

новой вероятностной меры оказывается асимптотически гауссовским с характеристиками, которые легко найти. Е:ли функционалы F и S дважды дифференцируемы в точке ф (требование гладкости S сводится к требованиям гладкости локальных характеристик семейства процессов

X?) и квадратичный функционал, соответствующий вто­ рой производной (F" — 5")(ф), строго отрицательно определен, для средних (1.23), (1.22) получаем точную асимптотику: при h | О

М* expffe-^X *)} - К0 ехр{^-1[7г(Ф) - 5(Ф)]}. (1.24)

Константа К0 выражается в виде математического ожида­ ния определенного функционала от гауссовского случай­ ного процесса.

5 21

БОЛЬШИЕ УКЛОНЕНИЯ ДЛЯ СЛУЧАЙНЫХ МЕР

399

Если F и S дифференцируемы v + 2 раза, для (1.22), (1.23) получается асимптотическое разложение вида

exp {h-4F(Ф) - S(ср)]}

(К0 + KJi + . . .

 

................................................+

K W 2 ] h W n + 0(fev/2)).

(1.25)

Аналогичные результаты получаются для математи­ ческих ожиданий функционалов вида G(Xh)exр xF(Xh)}.

Эта программа реализована в работе Ш и л ьд е р а

[1 ]

для

семейства случайных

процессов X? = Ъг^Ри)^

где

wt — винеровский процесс,

и в работах Д у б р о в ­

с к о г о

[1], [2], [3] для семейств локально безгранично

делимых марковских процессов, рассмотренных нами в

гл. 5.

Заметим, что в случае, рассмотренном у Ши л ь -

д е р а

[1 ], из намеченной нами схемы выпадает (ввиду ее

тривиальности) часть, касающаяся асимптотической гаус-

совости /г-1/2 [Xht — cpj относительно Ph. Точно так же в этой простой ситуации нет надобности в осознании свя­

зи применяемого метода с

методом

Г. Крамера.

§ 2. Большие уклонения для случайных мер

Рассмотрим винеровский

процесс

на отрезке [0, 1 ]

с отражением в концах отрезка. Для каждого борелевско-

го

множества Г cz [0,

1 ]

можно

рассмотреть случайную

величину лг(Г) — долю

времени,

проведенную процессом

l t

в множестве Г при

i е

[О, Т\:

 

 

 

 

т

 

 

 

Я г ( Г ) - И х г ( Ы * .

(2-1)

 

 

 

О

 

 

Случайная функция лт(Г) = лт(Г, со) является

вероят­

ностной мерой по Г. Пусть в пространстве мер на отрезке [О, 1 ] задана некоторая метрика р; например, определен­ ная равенством

р (р, V) =

sup I ц [0, х] — v [О, *] |.

(2.2)

 

 

0<х<1

 

 

Известно, что

для

любой

начальной точки

£0 = х е

^ [0,11 меры

лг при Т-+оо

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

1 сходятся

400

УТОЧНЕНИЯ И ОБОБЩЕНИЯ

1ГЛ. 9

к

лебеговой мере I на [0, 11. Если А — множество в

пространстве мер, находящееся на положительном рас­

стоянии

от лебеговой меры

Z,

то

Рх{яг е А }

0 при

Т

с»;

т. е. событие {со : яг

G

4 } относится к области

больших

уклонений.

определяет случайную

меру

 

Формула (2.1), конечно,

я г (Г) для любого измеримого случайного процесса

£*(со),

t ^

0. Если этот процесс, например, стационарный и эр-

годический, то по теореме

Биркгофа яг (Г) стремится

при Т -> оо к некоторой неслучайной мере иг(Г) (к одно­ мерному распределению процесса), и не стремящиеся к нулю отклонения меры яг от тпотносятся к области боль­ ших уклонений.

Семейства случайных мер, сходящиеся по вероятности к неслучайной мере, возникают и в других ситуациях, в частности, в задачах о пересечении уровня. Пусть,

например,

£* — стационарный

процесс

с достаточно

ре­

гулярными

траекториями, r\J = %т-г

>

0).

С

про­

цессом h можно связать семейство случайных мер ят

на

отрезке [0, 1 ], где ят(Г) — нормированное делением на

Т

число пересечений процессом r|f некоторого

заданного

уровня, происшедших при t е

Г с : [0,

1 ]. При

некото­

рых требованиях регулярности

на процесс

меры

 

лт

будут сходиться при Т —>* оо к мере Лебега,

умноженной

на среднее число пересечений на единицу

времени.

Для случайных мер можно рассматривать предельные теоремы различного рода. Здесь мы приведем некоторые результаты относительно поведения вероятностей боль­ ших уклонений. Грубая асимптотика вероятностей боль­ ших уклонений для мер (2.1), связанных с марковскими

процессами, изучалась в серии работ Д о н с к е р а

и

В а р а д а н а [1 ],

[2]. Независимо от этих авторов

общие результаты

были получены Г е р т и е р о м

[2],

[3]; Гертиер также применил эти результаты к марков­ ским и к диффузионным процессам и рассмотрел ряд при­ меров. Наше изложение близко к работам Гертнера.

Пусть , 3!) — измеримое пространство; обозначим В

пространство

ограниченных измеримых

функций на Е\

V — пространство всех конечных счетно-аддитивных функ­

ций множества (зарядов) на а-алгебре

Введем обозна­

чение < ц,

/ > = J f(x) р, (dx).

 

Е