Марковские плотности вероятности — различия между версиями
WikiSysop (обсуждение | вклад) |
WikiSysop (обсуждение | вклад) |
||
Строка 15: | Строка 15: | ||
{| width="100%" | {| 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> | | 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">'''( | + | | <div width="10%" align="right" style="color:#0000CC">'''(4.1)'''</div> |
|} | |} | ||
Строка 22: | Строка 22: | ||
{| width="100%" | {| width="100%" | ||
| width="90%" align="center"|<math> P(x_0,t_0 \Rightarrow x,t) = \delta (x-x_0)\;\;\;\;\;при\;\;\; t\to t_0. </math> | | 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">'''( | + | | <div width="10%" align="right" style="color:#0000CC">'''(4.2)'''</div> |
|} | |} | ||
Строка 31: | Строка 31: | ||
{| width="100%" | {| 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> | | 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">'''( | + | | <div width="10%" align="right" style="color:#0000CC">'''(4.3)'''</div> |
|} | |} | ||
Строка 40: | Строка 40: | ||
{| width="100%" | {| width="100%" | ||
| width="90%" align="center"|<math> P(x_1,x_3) = \int P(x_1,x_2,x_3) dx_2, </math> | | 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">'''( | + | | <div width="10%" align="right" style="color:#0000CC">'''(4.4)'''</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> определение условной вероятности, а в правую, с учётом ''марковости'' процесса, трёхточечную плотность вероятности (см. ( | + | где для краткости опущены времена <math>\textstyle t_i</math>. В (4.4) мы суммируем все возможные реализации "промежуточного" значения <math>\textstyle x_2</math>. В результате из трехточечной совместной плотности вероятности получается двуточечная. Подставим в левую часть <math>\textstyle P(x_1,x_3)=P(x_1)\cdot P(x_1\Rightarrow x_3)</math> определение условной вероятности, а в правую, с учётом ''марковости'' процесса, трёхточечную плотность вероятности ([[Модель аддитивного блуждания|см. (1.42)]]: |
:<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> | :<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> | ||
Строка 51: | Строка 51: | ||
{| width="100%" | {| 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> | | 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">'''( | + | | <div width="10%" align="right" style="color:#0000CC">'''(4.5)'''</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 \lessdot</math> H) имеет смысл проверить, что этому уравнению удовлетворяет гауссова плотность вероятности (4.1). Второе упражнение (<math>\textstyle \lessdot</math> H) состоит в записи уравнения Чепмена-Колмогорова для характеристических функций, если <math>\textstyle P(x_0,t_0\Rightarrow x,t)=P(x-x_0, t-t_0)</math>, и проверке марковости распределения Коши (4.3). |
− | Уравнению Чепмена-Колмогорова должны удовлетворять любые вероятности марковских процессов. Правда в таком виде оно слишком общее, и нам нужны его более конкретные представления. В уравнении () времена <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>. Мы найдём его при помощи стохастического дифференциального уравнения. Этот вывод покажет непосредственную связь двух математических аппаратов. | + | Уравнению Чепмена-Колмогорова должны удовлетворять любые вероятности марковских процессов. Правда в таком виде оно слишком общее, и нам нужны его более конкретные представления. В уравнении (4.5) времена <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>. Из (4.5) в следующем разделе мы получим уравнение относительно <math>\textstyle \{x_0,t_0\}</math>, которое называется первым уравнением Колмогорова. Аналогично выводится уравнение Фоккера-Планка, или второе уравнение Колмогорова относительно <math>\textstyle \{x,t\}</math>. Мы найдём его при помощи стохастического дифференциального уравнения. Этот вывод покажет непосредственную связь двух математических аппаратов. |
---- | ---- |
Текущая версия на 18:14, 15 марта 2010
Квазидетерминированное приближение << | Оглавление | >> Уравнение для плотности вероятности |
---|
Вернёмся к винеровскому процессу с нулевым сносом и единичной волатильностью. Так как случайная функция зависит от гауссовой переменной :
Невозможно разобрать выражение (синтаксическая ошибка): {\displaystyle x=x_0+\varepsilon\cdot \sqrt{t-t_0}\;\;\;\;\;\;\;или\;\;\;\;\;\;\;\varepsilon=(x-x_0)/\sqrt{t-t_0},}
то, воспользовавшись распределением Гаусса (см. стр. \pageref{prob_P_from_f}), можно записать условную плотность вероятности в виде:
(4.1)
|
Чем меньше разница , тем более высоким и узким будет колокол гауссианы, стремясь в пределе к дельта-функции Дирака:
Невозможно разобрать выражение (синтаксическая ошибка): {\displaystyle P(x_0,t_0 \Rightarrow x,t) = \delta (x-x_0)\;\;\;\;\;при\;\;\; t\to t_0. } | (4.2)
|
Она равна бесконечности при и нулю в других точках, так, что интеграл по в окрестности равен единице (см. Приложение М, стр. \pageref{math_delta_dirac_int}). Функции Дирака равна любая условная плотность вероятности при . Действительно, в бесконечно близкий к момент времени отлична от нуля только вероятность в окрестности начального значения .
К дельта-функции Дирака при стремится также условная плотность вероятности Коши:
(4.3)
|
Интеграл от этой функции по равен единице, среднее значение — . Однако моменты второго и более высоких порядков равны бесконечности. Соответственно равна бесконечности волатильность. В результате становятся вероятными очень большие выбросы случайных чисел. Подобные процессы называются процессами со скачками. У колокола распределения существует типичная ширина, пропорциональная . По мере удаления от начального момента времени происходит "расплывание" распределения вероятностей и очень быстрый уход процесса от начального значения . Поэтому в теории диффузных процессов мы не рассматриваем распределение Коши, хотя, как мы увидим чуть ниже, оно является марковским.
Плотность вероятности марковских процессов должна удовлетворять определённым уравнениям. Рассмотрим три последовательных момента времени , в которых принимает значения , и . Совместная плотность вероятности для и равна:
(4.4)
|
где для краткости опущены времена . В (4.4) мы суммируем все возможные реализации "промежуточного" значения . В результате из трехточечной совместной плотности вероятности получается двуточечная. Подставим в левую часть определение условной вероятности, а в правую, с учётом марковости процесса, трёхточечную плотность вероятности (см. (1.42):
Восстанавливая времена, получаем:
(4.5)
|
Это интегральное уравнение Чепмена-Колмогорова. В качестве упражнения ( H) имеет смысл проверить, что этому уравнению удовлетворяет гауссова плотность вероятности (4.1). Второе упражнение ( H) состоит в записи уравнения Чепмена-Колмогорова для характеристических функций, если , и проверке марковости распределения Коши (4.3).
Уравнению Чепмена-Колмогорова должны удовлетворять любые вероятности марковских процессов. Правда в таком виде оно слишком общее, и нам нужны его более конкретные представления. В уравнении (4.5) времена и могут быть удалены друг от друга как угодно далеко. Однако особый интерес представляет ситуация бесконечно близких времён. В результате глобальные свойства определяются из решений локальных дифференциальных уравнений. Так как условная плотность вероятности имеет две пары аргументов, то возможны, по крайней мере, два уравнения относительно и . Из (4.5) в следующем разделе мы получим уравнение относительно , которое называется первым уравнением Колмогорова. Аналогично выводится уравнение Фоккера-Планка, или второе уравнение Колмогорова относительно . Мы найдём его при помощи стохастического дифференциального уравнения. Этот вывод покажет непосредственную связь двух математических аппаратов.
Квазидетерминированное приближение << | Оглавление | >> Уравнение для плотности вероятности |
---|
Стохастический мир - простое введение в стохастические дифференциальные уравнения