td=TemporalData[];

DateListPlot[td,Filling->Bottom,Joined->True]

eproc=EstimatedProcess[td,SARMAProcess[2,1,{24,0,2}],ProcessEstimator->"MaximumConditionalLikelihood"]

WeakStationarity[eproc]

ListPlot[CorrelationFunction[#,{100}]&/@{td,eproc},PlotLegends->Placed[{"data","model"},{Right,Top}],PlotStyle->{PointSize[Medium],PointSize[Large]},ImageSize->300]