平滑仿真粒子轨道
用 BlockMap 和 TrimmedMean 的组合平滑来自粒子轨道的数据.
通过对精确轨道添加噪声仿真一个随机路径.
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]=
平滑带有 BlockMap 的仿真轨道，使用移动 TrimmedMean 来消除可能的异常值.
In[3]:=
smooth[r_] := BlockMap[TrimmedMean, data, r, 1]增加窗口大小获取更平滑轨道.
In[4]:=
Table[ListPlot[smooth[r],
AspectRatio -> Automatic], {r, {25, 50, 100, 200}}]Out[4]=
