Identifique processo Ornstein–Uhlenbeck regularmente amostrado como um processo autorregressivo 

O suporte aprimorado do Mathematica 10 para computação com fatias de processo permite que você use diretamente método de momentos para fatias de processo multivariado para estabelecer a equivalência de lei entre os dois processos.

Tanto o processo autorregressivo (AR) quanto o processo Ornstein–Uhlenbeck (OU) são processos Gaussianos, determinados por suas funções de média e covariância.

In[1]:=

X ar\[ScriptCapitalP] = ARProcess[c, {a1}, si2]; ou\[ScriptCapitalP] = OrnsteinUhlenbeckProcess[mu, si, th];

Use o método de momentos para distribuições de fatias para combinar os parâmetros do processo AR estacionário às restrições do processo OU estacionário a uma grade de tempo regular.

In[2]:=

X arSlice\[ScriptCapitalD] = ar\[ScriptCapitalP][{0, 1, 2, 3, 4}]; ouSlice\[ScriptCapitalD] = ou\[ScriptCapitalP][{dt, 2 dt, 3 dt, 4 dt, 5 dt}];

In[3]:=

X {sol1} = (Solve[{Mean[arSlice\[ScriptCapitalD]] == Mean[ouSlice\[ScriptCapitalD]], Covariance[arSlice\[ScriptCapitalD]] == Covariance[ouSlice\[ScriptCapitalD]]}, {c, a1, si2}] // FullSimplify)

Out[3]=

Confirme que as funções de covariância concordam com atraso simbólico.

In[4]:=

X Simplify[CovarianceFunction[ar\[ScriptCapitalP] /. sol1, h], dt > 0 && th > 0 && h \[Element] Integers]

Out[4]=

In[5]:=

X CovarianceFunction[ou\[ScriptCapitalP], h dt]

Out[5]=

Repita o exercício para processos AR e OU com condições iniciais.

In[6]:=

X ar\[ScriptCapitalP]2 = ARProcess[c, {a1}, si2, {x0ar}]; ou\[ScriptCapitalP]2 = OrnsteinUhlenbeckProcess[mu, si, th, x0];

In[7]:=

X arSlice\[ScriptCapitalD]2 = ar\[ScriptCapitalP]2[{0, 1, 2, 3, 4}]; ouSlice\[ScriptCapitalD]2 = ou\[ScriptCapitalP]2[{dt, 2 dt, 3 dt, 4 dt, 5 dt}];

In[8]:=

X {sol2} = (Solve[{Mean[arSlice\[ScriptCapitalD]2] == Mean[ouSlice\[ScriptCapitalD]2], Covariance[arSlice\[ScriptCapitalD]2] == Covariance[ouSlice\[ScriptCapitalD]2]}, {c, a1, si2, x0ar}] // FullSimplify)

Out[8]=

Use representação autorregressiva para amostrar o processo Ornstein–Uhlenbeck em uma grade de tempo regular.

In[9]:=

X RandomFunction[ ar\[ScriptCapitalP]2 /. sol2 /. {mu -> 1, si -> 2, th -> 2/3, x0 -> 10, dt -> 0.01}, {0, 200}, 5] // ListPlot[#, ImageSize -> 400] &

Out[9]=