2012-07-05 200 views
1

我有一些数据是分箱内未知曲线的积分。为了您的兴趣,数据是海浪能量,并且箱子用于方向,例如, 0-15度。如果可能的话,我想对曲线进行拟合,以保存箱内积分的数据。我试着用铅笔在笔记本上画草图,看起来好像可能。有没有人知道Python中的任何曲线拟合工具来做到这一点,例如在scipy插值子包中?已知积分的曲线拟合Python

在此先感谢

编辑:

感谢您的帮助。如果我这样做,看起来我会尝试本文第4部分推荐的方法:http://journals.ametsoc.org/doi/abs/10.1175/1520-0485%281996%29026%3C0136%3ATIOFFI%3E2.0.CO%3B2。从理论上讲,它基本上使用矩阵来根据每个频段之间的已知积分制作一些“假”数据。绘制时,这些数据会生成一个内插线图,保留积分。

回答

0

这是我的住所以外的一点,但我可以建议看看SciKits看看是否有什么可能有用。其他要浏览的软件包将是pandasStatsModels。祝你好运!

0

如果你有一个曲线f(x)这是一个近似的积分另一条曲线g(x),即f=int(g,x)那么这两个由Fundamental theorem of calculus是相关的,也就是说,你的原始功能是第一曲线g = df/dx的衍生物。因此,您可以使用numpy.diffhigher order methods中的任何一个来近似df/dx以获得原始曲线的估计值。

+0

谢谢,但问题是我没有曲线功能 - 这是一组数据。我应该更清楚一点。 – LaurieW 2012-07-26 15:54:56

+0

@LaurieW我明白 - 你不知道'g',但你知道'g'的近似值,这是你的分类数据。你对积分的估计当然是一个近似值,但它是你能做的最好的! – Hooked 2012-07-26 17:17:23

0

一种可能性:计算仓体积的累积和(np.cumsum),对其进行插值样条拟合,然后取导数得到曲线。

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

唯一的限制,如果它与您的情况相关,通过累计和的样条可能不是单调的,并且导数可能在某些时间间隔内是负的。

我想平滑直方图的文献着眼于积分/容积的相似约束,但我没有任何参考准备就绪。

+0

我认为我有类似的想法,但是我不知道在重新绘制曲线时是否会完全保留仓体积。 – LaurieW 2012-07-27 17:15:00

0

1/fit2histogram

你的问题是关于装修的直方图。我刚刚通过一些用于多变量模式分析的Python软件包PyMVPA的文档,并提出了一些用于直方图拟合的函数。一个例子是:PyMVPA

但是,我猜那套可用分布仅限于着名的分布。

2 /积分运算

如已经mentionned,下一个解决方案是在接近的整数值,以适应模型,以所得到的一组数据。要么你知道导数的明确表达式,要么使用计算推导:有限差分分析方法。