Skip to Content

§1. Броуновское движение

В §6 гл. 14 была введена модель симметричного случайного блуждания. Допустим, что для заданного значения nn выполнено следующее преобразование графика траектории блуждания (рис. 1 у Лагутина): масштаб по оси абсцисс сжат в nn раз, а по оси ординат — в n\sqrt n раз. При этом кусок графика на [0,n][0, n] преобразуется в (случайную) ломаную, заданную на [0,1][0, 1]. Обозначим её через Wn(t)W_n(t).

Известно, что в пределе при n→∞n \to \infty на [0,1][0, 1] возникает случайный процесс W(t)W(t), называемый стандартным броуновским движением или винеровским процессом (обозначение W(t)W(t) введено в честь Н. Винера, внёсшего значительный вклад в построение теории броуновского движения). Сходимость Wn(t)W_n(t) к W(t)W(t) можно понимать как сходимость конечномерных распределений процессов: для любых 0≤t1<t2<⋯<tk≤10 \le t_1 < t_2 < \dots < t_k \le 1 последовательность (Wn(t1),…,Wn(tk))→d(W(t1),…,W(tk))\bigl(W_n(t_1), \dots, W_n(t_k)\bigr) \xrightarrow{d} \bigl(W(t_1), \dots, W(t_k)\bigr) (см. П5).

По кирпичикам: почему по вертикали n\sqrt n, а не nn? Сумма SnS_n из nn шагов ±1\pm1 имеет дисперсию nn, т. е. типичное отклонение n\sqrt n. Сожми в nn раз — получишь ноль (закон больших чисел), сожми в n\sqrt n — получишь нормальную величину с дисперсией 1 в момент t=1t = 1 (центральная предельная теорема). А в момент tt ломаная прошла tntn шагов, и Wn(t)≈Stn/n∼N(0,t)W_n(t) \approx S_{tn}/\sqrt n \sim \mathcal N(0, t).

Смысл по-пацански: пацан идёт по улице и на каждом шаге бросает монетку: орёл — шаг вправо, решка — шаг влево. Посмотри на его маршрут издалека — когда шагов миллион, а ты сжал картинку правильно, — увидишь не ломаную, а дрожащую непрерывную кривую. Это и есть броуновское движение.

Процесс W(t)W(t) характеризуется следующими свойствами:

  1. W(0)=0W(0) = 0;
  2. приращения W(t2)−W(t1),…,W(tk)−W(tk−1)W(t_2) - W(t_1), \dots, W(t_k) - W(t_{k-1}) независимы;
  3. W(t)−W(s)∼N(0,t−s)W(t) - W(s) \sim \mathcal N(0, t - s) при 0≤s<t0 \le s < t;
  4. траектории W(t)W(t) непрерывны с вероятностью 1.

Броуновское движение как предел блужданий

3,0−3,001t

Три симметричных блуждания по ±1 из n шагов (§6 гл. 14), сжатые по оси t в n раз, а по оси y — в √n раз. При n = 4 это грубые ломаные, при n = 2048 — уже «пушистые» кривые: в пределе возникает стандартное броуновское движение W(t). Почему по y именно √n: к моменту n блуждание уходит от нуля на величину порядка √n (дисперсия суммы n шагов ±1 равна n).

Нигде не дифференцируемые траектории

Однако траектории процесса W(t)W(t) нигде не дифференцируемы с вероятностью 1. Этот факт соответствует физической природе процесса, а математически объясняется разной скоростью сжатия в горизонтальном и вертикальном направлениях ломаной случайного блуждания, вследствие чего каждый локальный максимум (минимум) траектории становится с ростом nn всё более «острым».

По кирпичикам: по свойству 3 приращение за время Δt\Delta t имеет порядок Δt\sqrt{\Delta t}, а точнее M∣W(t+Δt)−W(t)∣=2Δt/π\mathsf M|W(t + \Delta t) - W(t)| = \sqrt{2\Delta t/\pi}. Разностное отношение ∣ΔW∣/Δt|\Delta W|/\Delta t — порядка 1/Δt→∞1/\sqrt{\Delta t} \to \infty при Δt→0\Delta t \to 0. У гладкой функции оно стремилось бы к ∣f′(t)∣|f'(t)|, а здесь — улетает в бесконечность.

Проверка

Какова средняя «крутизна» M∣ΔW∣/Δt\mathsf M|\Delta W|/\Delta t на сетке с шагом Δt=1/1024\Delta t = 1/1024? (Два знака после запятой.)

У броуновской траектории нет производной

mсредняя |ΔW|/Δt√(2m/π)
163,23,2
645,86,4
25612,712,8
102425,525,5
409650,951,1

Средняя крутизна |ΔW|/Δt на сетке с шагом Δt = 1/m. Для гладкой функции она стабилизировалась бы около |f′|, а у W(t) растёт как √m (теория: √(2m/π)): приращение за время Δt порядка √Δt, а делим на Δt. Отсюда — траектории нигде не дифференцируемы, и каждый локальный максимум с ростом n становится всё «острее».

Смысл по-пацански: под микроскопом гладкая кривая становится прямой. А броуновская — нет: увеличь в 100 раз, и она будет такой же рваной, как была. Касательную к ней не проведёшь ни в одной точке.

Ковариационная функция

Из свойств 1 и 3 вытекает, что MW(t)=0\mathsf MW(t) = 0 при всех tt. Вычислим ковариационную функцию R(s,t)=MW(s)W(t)R(s, t) = \mathsf MW(s)W(t). При s<ts < t, используя свойства 2 и 3, находим, что

R(s,t)=M[W(s)−W(0)][W(t)−W(s)]+MW2(s)=M[W(s)−0]⋅M[W(t)−W(s)]+s=s.R(s, t) = \mathsf M\bigl[W(s) - W(0)\bigr]\bigl[W(t) - W(s)\bigr] + \mathsf MW^2(s) = \mathsf M[W(s) - 0]\cdot\mathsf M[W(t) - W(s)] + s = s.

В общем случае, очевидно, R(s,t)=min⁡{s,t}R(s, t) = \min\{s, t\}.

По кирпичикам: разбили W(t)=W(s)+[W(t)−W(s)]W(t) = W(s) + [W(t) - W(s)]. Первое слагаемое — «то, что уже прошли», второе — «то, что будет», и оно от прошлого не зависит (свойство 2) и в среднем равно нулю. Поэтому ковариация W(s)W(s) с W(t)W(t) — это просто дисперсия общего куска пути, DW(s)=s\mathsf DW(s) = s.

Смысл по-пацански: где был наш пацан с монеткой в 10 вечера и где он в полночь — связано ровно настолько, насколько полночь «помнит» 10 вечера: весь путь до 10 у них общий, а дальше — новые шатания, к прошлому отношения не имеющие.

Проверь моделированием:

Загрузка редактора…

Проверь себя

Тест

Почему при переходе от блуждания к броуновскому движению вертикальный масштаб сжимают в √n, а не в n раз?

Тест

Чему равна ковариация W(0,2) и W(0,9)?
Обновлено