5

最近移民到 Python 和使用 Python 进行科学计算。这是一个避免任何可能已经存在的代码重复的问题。

我有一个字段在常规网格中作为 x 和 y 的函数进行采样。我想对数据进行插值,不仅可以获得网格上任意点的字段值,还可以获得一阶和二阶导数。如果我使用 interp2d 对它进行双三次插值,我可以获得该字段的值。

有没有人建议如何使用现有的 numpy 或 scipy 函数获得该字段的一阶和二阶导数?

谢谢!

4

1 回答 1

3

scipy.interpolate.interp2d.__call__方法具有选项dxdy可以在点评估更高的导数(至少从版本 0.14.0 开始)。

于 2015-04-02T14:49:22.223 回答