Suavize uma trajetória simulada de partículas

Use uma combinação de BlockMap e TrimmedMean para suavizar dados de uma trejetória de partículas.

Simule um comportamento estocástico acrescentando ruído a uma trajetória exata.

In[1]:=

exact[u_] := {u Cos[u], u Sin[u]} noise[u_] := RandomReal[CauchyDistribution[0, 1], 2] data = Table[exact[u] + noise[u], {u, 0, 6 \[Pi], 1/100}];

In[2]:=

ListPlot[data, AspectRatio -> Automatic]

Out[2]=

Suavize a trajetória simulada com BlockMap usando um TrimmedMean móvel para eliminar possíveis valores extremos.

In[3]:=

smooth[r_] := BlockMap[TrimmedMean, data, r, 1]

Aumente o tamanho da janela para obter uma trajetória mais suave.

In[4]:=

Table[ListPlot[smooth[r], AspectRatio -> Automatic], {r, {25, 50, 100, 200}}]

Out[4]=