Марковские плотности вероятности — различия между версиями

Материал из synset
Перейти к: навигация, поиск
м (Защищена страница «Марковские плотности вероятности» ([edit=sysop] (бессрочно) [move=sysop] (бессрочно)))
Строка 7: Строка 7:
  
  
 +
<math>\textstyle \bullet</math> Вернёмся к винеровскому процессу с нулевым сносом <math>\textstyle \mu=0</math> и единичной волатильностью. Так как случайная функция <math>\textstyle x(t)</math> зависит от гауссовой переменной <math>\textstyle \varepsilon</math>:
 +
 +
:<center><math>x=x_0+\varepsilon\cdot \sqrt{t-t_0}\;\;\;\;\;\;\;или\;\;\;\;\;\;\;\varepsilon=(x-x_0)/\sqrt{t-t_0},</math></center>
 +
 +
то, воспользовавшись распределением Гаусса (см. стр. \pageref{prob_P_from_f}), можно записать условную плотность вероятности в виде:
 +
 +
{| width="100%"
 +
| width="90%" align="center"|<math> P(x_0,t_0 \Rightarrow x,t) = \frac{1}{\sqrt{2\pi\cdot(t-t_0)}} \cdot \exp \left\{ -\frac{1}{2}\,\frac{(x-x_0)^2}{t-t_0}\right\}. </math>
 +
| <div width="10%" align="right" style="color:#0000CC">'''(EQN)'''</div>
 +
|}
 +
 +
Чем меньше разница <math>\textstyle t-t_0</math>, тем более высоким и узким будет колокол гауссианы, стремясь в пределе <math>\textstyle t\to t_0</math> к ''дельта-функции Дирака'':
 +
 +
{| width="100%"
 +
| width="90%" align="center"|<math> P(x_0,t_0 \Rightarrow x,t) = \delta (x-x_0)\;\;\;\;\;при\;\;\; t\to t_0. </math>
 +
| <div width="10%" align="right" style="color:#0000CC">'''(EQN)'''</div>
 +
|}
 +
 +
Она равна бесконечности при <math>\textstyle x=x_0</math> и нулю в других точках, так, что интеграл по <math>\textstyle x</math> в окрестности <math>\textstyle x_0</math> равен единице (см. Приложение М, стр. \pageref{math_delta_dirac_int}). Функции Дирака равна любая условная плотность вероятности при <math>\textstyle t=t_0</math>. Действительно, в бесконечно близкий к <math>\textstyle t_0</math> момент времени отлична от нуля только вероятность в окрестности начального значения <math>\textstyle x\approx x_0</math>.
 +
 +
К дельта-функции Дирака при <math>\textstyle t\to t_0</math> стремится также условная плотность вероятности Коши:
 +
 +
{| width="100%"
 +
| width="90%" align="center"|<math> P(x_0,t_0 \Rightarrow x,t) = \frac{(t-t_0)/\pi}{(x-x_0)^2+(t-t_0)^2}. </math>
 +
| <div width="10%" align="right" style="color:#0000CC">'''(EQN)'''</div>
 +
|}
 +
 +
Интеграл от этой функции по <math>\textstyle x</math> равен единице, среднее значение &mdash; <math>\textstyle x_0</math>. Однако моменты второго и более высоких порядков равны бесконечности. Соответственно равна бесконечности волатильность. В результате становятся вероятными очень большие выбросы случайных чисел. Подобные процессы называются процессами со скачками. У колокола распределения существует типичная ширина, пропорциональная <math>\textstyle t-t_0</math>. По мере удаления от начального момента времени происходит "расплывание" распределения вероятностей и очень быстрый уход процесса от начального значения <math>\textstyle x_0</math>. Поэтому в теории диффузных процессов мы не рассматриваем распределение Коши, хотя, как мы увидим чуть ниже, оно является марковским.
 +
 +
<math>\textstyle \bullet</math> Плотность вероятности марковских процессов должна удовлетворять определённым уравнениям. Рассмотрим три последовательных момента времени <math>\textstyle t_1<t_2<t_3</math>, в которых <math>\textstyle x(t)</math> принимает значения <math>\textstyle x_1</math>, <math>\textstyle x_2</math> и <math>\textstyle x_3</math>. Совместная плотность вероятности для <math>\textstyle x_1</math> и <math>\textstyle x_3</math> равна:
 +
 +
{| width="100%"
 +
| width="90%" align="center"|<math> P(x_1,x_3) = \int P(x_1,x_2,x_3) dx_2, </math>
 +
| <div width="10%" align="right" style="color:#0000CC">'''(EQN)'''</div>
 +
|}
 +
 +
где для краткости опущены времена <math>\textstyle t_i</math>. В () мы суммируем все возможные реализации "промежуточного" значения <math>\textstyle x_2</math>. В результате из трехточечной совместной плотности вероятности получается двуточечная. Подставим в левую часть <math>\textstyle P(x_1,x_3)=P(x_1)\cdot P(x_1\Rightarrow x_3)</math> определение условной вероятности, а в правую, с учётом ''марковости'' процесса, трёхточечную плотность вероятности (см. (), на стр. \pageref{chain_markov}):
 +
 +
:<center><math>P(x_1,x_2,x_3)=P(x_1)\cdot P(x_1\Rightarrow x_2) P(x_2\Rightarrow x_3).</math></center>
 +
 +
Восстанавливая времена, получаем:
 +
 +
{| width="100%"
 +
| width="90%" align="center"|<math> P(x_1,t_1 \Rightarrow x_3, t_3) = \int P(x_1,t_1 \Rightarrow x_2, t_2)\cdot P(x_2,t_2 \Rightarrow x_3, t_3) \;dx_2. </math>
 +
| <div width="10%" align="right" style="color:#0000CC">'''(EQN)'''</div>
 +
|}
 +
 +
Это интегральное ''уравнение Чепмена-Колмогорова''. В качестве упражнения (<math>\textstyle \lessdot</math> H) имеет смысл проверить, что этому уравнению удовлетворяет гауссова плотность вероятности (). Второе упражнение (<math>\textstyle \lessdot</math> H) состоит в записи уравнения Чепмена-Колмогорова для характеристических функций, если <math>\textstyle P(x_0,t_0\Rightarrow x,t)=P(x-x_0, t-t_0)</math>, и проверке марковости распределения Коши ().
 +
 +
Уравнению Чепмена-Колмогорова должны удовлетворять любые вероятности марковских процессов. Правда в таком виде оно слишком общее, и нам нужны его более конкретные представления. В уравнении () времена <math>\textstyle t_1,</math> <math>\textstyle t_2</math> и <math>\textstyle t_3</math> могут быть удалены друг от друга как угодно далеко. Однако особый интерес представляет ситуация бесконечно близких времён. В результате глобальные свойства <math>\textstyle P(x_0, t_0\Rightarrow x,t)</math> определяются из решений локальных дифференциальных уравнений. Так как условная плотность вероятности имеет две пары аргументов, то возможны, по крайней мере, два уравнения относительно <math>\textstyle \{x_0,t_0\}</math> и <math>\textstyle \{x,t\}</math>. Из () в следующем разделе мы получим уравнение относительно <math>\textstyle \{x_0,t_0\}</math>, которое называется первым уравнением Колмогорова. Аналогично выводится уравнение Фоккера-Планка, или второе уравнение Колмогорова относительно <math>\textstyle \{x,t\}</math>. Мы найдём его при помощи стохастического дифференциального уравнения. Этот вывод покажет непосредственную связь двух математических аппаратов.
  
 
----
 
----

Версия 18:55, 21 февраля 2010

Квазидетерминированное приближение << Оглавление >> Уравнение для плотности вероятности


Вернёмся к винеровскому процессу с нулевым сносом и единичной волатильностью. Так как случайная функция зависит от гауссовой переменной :

Невозможно разобрать выражение (синтаксическая ошибка): {\displaystyle x=x_0+\varepsilon\cdot \sqrt{t-t_0}\;\;\;\;\;\;\;или\;\;\;\;\;\;\;\varepsilon=(x-x_0)/\sqrt{t-t_0},}

то, воспользовавшись распределением Гаусса (см. стр. \pageref{prob_P_from_f}), можно записать условную плотность вероятности в виде:

(EQN)

Чем меньше разница , тем более высоким и узким будет колокол гауссианы, стремясь в пределе к дельта-функции Дирака:

Невозможно разобрать выражение (синтаксическая ошибка): {\displaystyle P(x_0,t_0 \Rightarrow x,t) = \delta (x-x_0)\;\;\;\;\;при\;\;\; t\to t_0. }
(EQN)

Она равна бесконечности при и нулю в других точках, так, что интеграл по в окрестности равен единице (см. Приложение М, стр. \pageref{math_delta_dirac_int}). Функции Дирака равна любая условная плотность вероятности при . Действительно, в бесконечно близкий к момент времени отлична от нуля только вероятность в окрестности начального значения .

К дельта-функции Дирака при стремится также условная плотность вероятности Коши:

(EQN)

Интеграл от этой функции по равен единице, среднее значение — . Однако моменты второго и более высоких порядков равны бесконечности. Соответственно равна бесконечности волатильность. В результате становятся вероятными очень большие выбросы случайных чисел. Подобные процессы называются процессами со скачками. У колокола распределения существует типичная ширина, пропорциональная . По мере удаления от начального момента времени происходит "расплывание" распределения вероятностей и очень быстрый уход процесса от начального значения . Поэтому в теории диффузных процессов мы не рассматриваем распределение Коши, хотя, как мы увидим чуть ниже, оно является марковским.

Плотность вероятности марковских процессов должна удовлетворять определённым уравнениям. Рассмотрим три последовательных момента времени , в которых принимает значения , и . Совместная плотность вероятности для и равна:

(EQN)

где для краткости опущены времена . В () мы суммируем все возможные реализации "промежуточного" значения . В результате из трехточечной совместной плотности вероятности получается двуточечная. Подставим в левую часть определение условной вероятности, а в правую, с учётом марковости процесса, трёхточечную плотность вероятности (см. (), на стр. \pageref{chain_markov}):

Восстанавливая времена, получаем:

(EQN)

Это интегральное уравнение Чепмена-Колмогорова. В качестве упражнения ( H) имеет смысл проверить, что этому уравнению удовлетворяет гауссова плотность вероятности (). Второе упражнение ( H) состоит в записи уравнения Чепмена-Колмогорова для характеристических функций, если , и проверке марковости распределения Коши ().

Уравнению Чепмена-Колмогорова должны удовлетворять любые вероятности марковских процессов. Правда в таком виде оно слишком общее, и нам нужны его более конкретные представления. В уравнении () времена и могут быть удалены друг от друга как угодно далеко. Однако особый интерес представляет ситуация бесконечно близких времён. В результате глобальные свойства определяются из решений локальных дифференциальных уравнений. Так как условная плотность вероятности имеет две пары аргументов, то возможны, по крайней мере, два уравнения относительно и . Из () в следующем разделе мы получим уравнение относительно , которое называется первым уравнением Колмогорова. Аналогично выводится уравнение Фоккера-Планка, или второе уравнение Колмогорова относительно . Мы найдём его при помощи стохастического дифференциального уравнения. Этот вывод покажет непосредственную связь двух математических аппаратов.


Квазидетерминированное приближение << Оглавление >> Уравнение для плотности вероятности

Стохастический мир - простое введение в стохастические дифференциальные уравнения