"ref_1.gif"

"ref_2.gif"

Случайный марковский процесс:  блуждание частиы по целочисленной решётке.

1. Так как мы желаем только проиллюстрировать случайный процесс, матрицу переходных вероятностей  полагаем случайной. Генерация матрицы переходных вероятностей P.

In[35]:=

"accidental_error_1.gif"

Out[35]=

"accidental_error_2.gif"

In[36]:=

"accidental_error_3.gif"

Out[36]=

"accidental_error_4.gif"

In[37]:=

"accidental_error_5.gif"

Out[37]//MatrixForm=

"accidental_error_6.gif"

Проверяем, равна ли 1 ссумма элементов каждой строки.

In[38]:=

"accidental_error_7.gif"

Out[38]=

"accidental_error_8.gif"

2. Определяем сучайную величину Direct, принимающую значения up, right, down, left, которая определяет направление движения блуждающей частицы.

In[5]:=

"accidental_error_9.gif"

In[6]:=

"accidental_error_10.gif"

In[7]:=

"accidental_error_11.gif"

In[8]:=

"accidental_error_12.gif"

In[9]:=

"accidental_error_13.gif"

3, Генерируем случайную последовательность направлений.

In[51]:=

"accidental_error_14.gif"

4. Вычисление последовательности частичных сумм элементов списка lst.

In[52]:=

"accidental_error_15.gif"

In[53]:=

"accidental_error_16.gif"

In[61]:=

"accidental_error_17.gif"

5  Визуализация случайного блуждания (50 шагов).

In[63]:=

"accidental_error_18.gif"

Out[63]=

"accidental_error_19.gif"

6. Вычисление стационарной матрицы переходов "accidental_error_20.gif"

In[65]:=

"accidental_error_21.gif"

In[66]:=

"accidental_error_22.gif"

In[67]:=

"accidental_error_23.gif"

Out[67]//MatrixForm=

"accidental_error_24.gif"

Блуждание частицы в стационарном режиме (матрица переходных вероятностей "accidental_error_25.gif".

1, Генерируем случайную последовательность направлений.

In[68]:=

"accidental_error_26.gif"

2. Вычисление последовательности частичных сумм элементов списка laccd.

In[69]:=

"accidental_error_27.gif"

In[70]:=

"accidental_error_28.gif"

In[73]:=

"accidental_error_29.gif"

5  Визуализация случайного блуждания (50 шагов).

In[74]:=

"accidental_error_30.gif"

Out[74]=

"accidental_error_31.gif"

"ref_1.gif"

"ref_2.gif"

Spikey Created with Wolfram Mathematica 7.0
Hosted by uCoz