Оценка нерегулярных выборок случайных процессов
Сформируем нерегулярную выборку с помощью функции OrnsteinUhlenbeckProcess.
In[1]:=
![Click for copyable input](assets.ru/estimation-of-irregularly-sampled-random-processes/In_81.png)
sample = TimeSeriesResample[
RandomFunction[
OrnsteinUhlenbeckProcess[0, .1, .3], {0, 100, .1}], {Sort[
RandomReal[100, 1000]]}]
Out[1]=
![](assets.ru/estimation-of-irregularly-sampled-random-processes/O_60.png)
код на языке Wolfram Language целиком
Out[2]=
![](assets.ru/estimation-of-irregularly-sampled-random-processes/O_61.png)
Вычислим параметры процесса на основе нерегулярных выборок данных.
In[3]:=
![Click for copyable input](assets.ru/estimation-of-irregularly-sampled-random-processes/In_83.png)
EstimatedProcess[sample,
OrnsteinUhlenbeckProcess[\[Mu], \[Sigma], \[Theta]]]
Out[3]=
![](assets.ru/estimation-of-irregularly-sampled-random-processes/O_62.png)
Возьмем цены на акции компании General Electric с 1 января 2013 года, и конвертируем их в временные ряды с помощью функции TemporalData .
In[4]:=
![Click for copyable input](assets.ru/estimation-of-irregularly-sampled-random-processes/In_84.png)
price = TemporalData[FinancialData["GE", "Jan. 1, 2013"]]
Out[4]=
![](assets.ru/estimation-of-irregularly-sampled-random-processes/O_63.png)
код на языке Wolfram Language целиком
Out[5]=
![](assets.ru/estimation-of-irregularly-sampled-random-processes/O_64.png)
Штамп времени данных о ценах акций неоднороден.
In[6]:=
![Click for copyable input](assets.ru/estimation-of-irregularly-sampled-random-processes/In_86.png)
MinMax[Differences[price["Times"]]]
Out[6]=
![](assets.ru/estimation-of-irregularly-sampled-random-processes/O_65.png)
Предположим, что логарифм цены соответствует характеристикам дробного броуновского движения (см. FractionalBrownianMotionProcess) и вновь вычислим параметры процесса.
In[7]:=
![Click for copyable input](assets.ru/estimation-of-irregularly-sampled-random-processes/In_87.png)
EstimatedProcess[Log[price],
FractionalBrownianMotionProcess[\[Mu], \[Sigma], h]]
Out[7]=
![](assets.ru/estimation-of-irregularly-sampled-random-processes/O_66.png)