使用已知积分Python进行曲线拟合

Question

我有一些数据是箱中未知曲线的积分。 为了您的利益，数据为海浪能量，垃圾箱为方向，例如0-15度。 如果可能的话，我想在数据上拟合一条曲线，以节省分档内的积分。 我试着用铅笔在记事本上绘制草图，似乎有可能。 有谁知道Python中有任何曲线拟合工具可以做到这一点，例如在scipy插值子包中？

提前致谢

编辑：

谢谢您的帮助。 如果这样做，看来我将尝试本文第4节中推荐的方法： http : //journals.ametsoc.org/doi/abs/10.1175/1520-0485%281996%29026%3C0136%3ATIOFFI ％3E2.0.CO％3B2 。 从理论上讲，它基本上是使用矩阵从每个波段之间的已知积分中得出一些“伪”数据。 绘制后，该数据将生成一个插值线图，该图保留了积分。

Answer 1

距离我的bailiwick有点远 ，但是我建议您看看SciKits ，看看是否有任何有用的东西。 要浏览的其他软件包是pandas和StatsModels 。 祝好运！

Answer 2

如果您有一条曲线f(x)近似于另一条曲线g(x)的积分，即f=int(g,x)则两者通过微积分的基本定理进行关联，即您的原始函数是第一条曲线的导数g = df/dx 。 这样，您可以使用numpy.diff或任何高阶方法来近似df/dx以获得原始曲线的估计值。

Answer 3

一种可能性是：计算仓体积的累积总和（np.cumsum），对它进行插值样条拟合，然后采用导数获得曲线。

scipy样条曲线具有计算导数的方法。

如果与您的情况相关，则是唯一的限制，累积总和的样条曲线可能不是单调的，并且在某些时间间隔上导数可能为负。

我猜想有关平滑直方图的文献对积分/箱的体积有类似的约束，但是我还没有任何参考资料。

Answer 4

1 / fit2直方图

您的问题是关于拟合直方图。 我刚看过一些文档样本，其中包含用于多变量模式分析的Python软件包，PyMVPA，并提出了一些用于直方图拟合的功能。 一个例子在这里： PyMVPA 。

但是，我想可用的发行版集仅限于著名的发行版。

2 /积分计算

如前所述，下一个解决方案是近似积分值，并使模型适合所得数据集。 您要么知道导数的显式表达式，要么使用计算导数：有限差分，解析方法。