0

我正在尝试使用数学中的移动平均线平滑 3D 直方图。我知道有一个叫做 smoothhistogram3D 的函数，它接近我想要的，但是，它似乎只能选择使用分布函数来平滑曲线。我能够通过修改这个 stackoverflow 答案来创建一个平滑 2D 直方图的函数，以包含一个 interpolationOrder 和移动平均功能。

MovAvgHistoPlot[MovAvg_, dx_] := Module[{histList, transposedHistList, histListAvg},
  histList = HistogramList[data, {dx}];

  transposedHistList = Transpose[{histList[[1]],ArrayPad[histList[[2]], {0, 1}, 
    "Fixed"]}];
  histListAvg = MovingAverage[transposedHistList, MovAvg];
  histPlot = ListPlot[histListAvg, InterpolationOrder -> 3, Joined -> True, 
    AxesOrigin -> {histListAvg[[1, 1]], 0}, PlotRange -> All, InterpolationOrder -> 3,
    PlotStyle -> Black, FrameLabel -> {"kTh", "Ion Intensity"},Frame -> {{True, False}, 
    {True, False}}, Axes -> False, ImageSize -> Large]
]

Manipulate[MovAvgHistoPlot[MovAvg, dx], {{MovAvg, 1, "Moving Average"}, 1, 500, 1}, 
  ContinuousAction -> False]

我尝试使用下面的代码将其扩展到第三维，但没有成功。

MovAvgHistPlot3D[MovAvg_] := Module[{HistList3D, XAndZGroupedValues, XValues, ZValues, 
  XAndZValues, YValues, ListPlot3DPoints},

  HistList3D = HistogramList[data];
  XAndZValues = Flatten[Partition[Table[Riffle[HistList3D[[1, 2]], 
    HistList3D[[1, 1, i]], {1, -1, 2}], {i,Length[HistList3D[[1, 1]]]}], {1, 2}]];
  YValues = Flatten[ArrayPad[Map[ArrayPad[##, {0, 1}, "Fixed"] &, HistList3D[[2]]], 
    {0, 1}, "Fixed"]]; 
  ListPlot3DPoints = Partition[Riffle[XAndZValues, YValues, {3, -1, 3}], 3];
  ListPlot3DPointsAvg = MovingAverage[ListPlot3DPoints, MovAvg];
  ListPlot3D[ListPlot3DPointsAvg, InterpolationOrder -> 3, Joined -> True]
  ]

Manipulate[MovAvgHistPlot3D[MovAvg], {{MovAvg, 1, "Moving Average"}, 1, 1000, 1},
  ContinuousAction -> False]

但是，3D 函数使用我的数据集输出此图像：http: //imgur.com/MJeBbwW

我首先尝试使用与此类似的方法，但可以选择使用移动平均线对其进行平滑处理：

ListPlot3D[HistogramList[filteredData1]]

但是，它输出了这样的图像：

http://imgur.com/Bkj0R9W
(Sorry, can't post more than two links due to lack of reputation points)

我想要一个与 smoothhistogram3D 的输出非常相似的数据集，但可以选择使用移动平均线进行平滑。

http://imgur.com/NRj6V2R

有什么建议么？有没有更简单的方法我没有意识到？

对不起，我意识到代码，尤其是第二部分，几乎不可读。我是mathematica的新手，只是想让它工作。

这也是我第一次在堆栈溢出上发帖，所以请原谅任何格式或指南错误。

4

0 回答 0