Новое в системе Wolfram
Mathematica
9
◄
предыдущая
|
следующая
►
Новое в системе Wolfram
Mathematica
9
›
Временные ряды и стохастические дифференциальные уравнения
Моделирование сезонных данных
Фиттинг модели
SARMAProcess
по почасовым замерам температуры в августе.
In[1]:=
X
data = Select[WeatherData["Champaign", "Temperature", {2012, 8}], FreeQ[#, Missing] &];
In[2]:=
X
DateListPlot[data, Filling -> Bottom, Joined -> True]
Out[2]=
Нахождение параметров процесса.
In[3]:=
X
eproc = EstimatedProcess[data[[All, 2]], SARMAProcess[{}, {}, {24, {c}, {d}}, v], ProcessEstimator -> "MethodOfMoments"]
Out[3]=
Проверка, что процесс является стационарным в широком смысле.
In[4]:=
X
WeakStationarity[eproc]
Out[4]=
Сравнение корреляционных функций.
In[5]:=
X
ListPlot[CorrelationFunction[#, {100}] & /@ {eproc, data[[All, 2]]}, PlotStyle -> {PointSize[Large], PointSize[Medium]}, PlotRange -> All, PlotLegends -> Placed[{"Model", "Data"}, {Right, Top}], ImageSize -> 300]
Out[5]=