Моделирование сезонных данных 

Фиттинг модели 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]=