Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
лекция_сп1_проектор.doc
Скачиваний:
0
Добавлен:
01.04.2025
Размер:
162.82 Кб
Скачать

2. Винеровский процесс.

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

Дискретным аналогом этого процесса может служить следующая модель случайного блуждания. Частица меняет свое положение лишь в дискретные моменты времени, кратные t. Изменение положения происходит таким образом, что находясь в точке х, частица, независимо от предшествующего поведения, переходит с равными вероятностями в одну из соседних точек х+х или х-х (речь идет лишь об одной координате частицы.) В пределе t0, х0, получается непрерывное случайное движение, характерное для процесса броуновского движения.

Обозначим (t) – положение частицы в момент времени t. Пусть (0)=0. При дискретном блуждании частица за время t совершает n=t/t шагов. Пусть Sn - число шагов, совершаемых в положительном направлении.

Тогда (t)=Snх-(n- Sn)х=(2Sn-n)х.

Поскольку (0)=0, то (t)=[(t)-(s)]+ [(s)-(0)] для любого s в интервале [0,t].

В описанной модели случайного блуждания с.в. [(t)-(s)] и [(s)-(0)] являются независимыми. Причем распределение вероятностей приращения [(t)-(s)] такое же, как и приращения [(t-s)-(0)]. Поэтому

D(t)=2(t)=2(s+(t-s))=D[(t)-(s)]+D[(s)-(0)]=2(s)+2(t-s).

Следовательно, дисперсия есть линейная функция времени:

D(t)=2t,

где 2>0 - некоторая постоянная, называемая коэффициентом диффузии.

С другой стороны, дисперсия с.в. (t), описывающей смещение за n=t/t шагов, имеет вид

D(t)=D((2Sn-n)х)= D(2Snх -nх)=4(х)2DSn=(х)2n=(х)2t/t.

В итоге получаем следующее соотношение между х и t: 2=(х)2/t.

Положение частицы описывается формулой (t)=(2Sn-n)х, где Sn – число успехов в n испытаниях Бернулли.

Согласно теореме Бернулли с.в. имеет стандартное нормальное распределение и .

Т.О. и

Т.О. (t) – нормальная с.в. N(0,t2).

Т.О., для моделирования броуновского движения можно определить следующий процесс.

Определение. Винеровским процессом называется случайный процесс с независимыми приращениями, обладающий следующими свойствами:

  1. (0)=0.

  2.  t>s с.в. [(t)-(s)] имеет нормальное распределение N(0,2(t-s)).

Моментные функции Винеровского процесса.

1. mt=0.

2. Если t>s, то

R(t,s)= M(t)(s) = M((t)-(s))(s)+M((s))2 = M((t)-(s))M(s)+M((s))2= D(s)=2s.

Если t<s, то R(t,s)=2t.

Окончательно получим R(t,s)=2min(t,s).

Т.О. корреляционные функции Пуассоновского и Винеровского процесса совпадают, но это совершенно разные процессы.

Пример на Винеровский процесс.

Пусть x – время первого достижения частицей точки х, то есть (x)=х.

Найдем распределение для с.в. x.

Если (t)х, то xt поскольку траектории частицы непрерывны.

P((t)х|xt)= P((t)х)/P(xt).

Если в момент времени x частица находилась в точке х, то далее она могла смещаться как вправо, так и влево, и поэтому могла с равной вероятностью оказаться как справа, так и слева от точки х.

Поэтому P((t)х|xt)=1/2= P((t)х)/P(xt) и

P(xt)= P((t)х)/(1/2)=2P((t)х) .

Плотность вероятности с.в. x найдем, продифференцировав последнее выражение по t:

.

Замечание. Винеровский процесс является представителем еще одного класса процессов, так называемых Гауссовских с.п.

Гауссовским с.п. называется процесс, у которого все распределения являются гауссовскими многомерными распределениями.

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

Гауссовский и винеровские процессы имеют очень широкие области применения, включающие радиотехнику, теорию диффузии, модели шумов и.т.д.

17

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]