numpy中的高阶有限差分

时间:2018-04-06 16:48:16

标签: python numpy

我在2D和3D numpy数组上有采样函数,我需要一种从这些数组中获取偏导数的方法。我可以手动编码有限差分方案,但我需要的不仅仅是二阶精度,可能是四阶甚至六阶。对于更高精度的订单编码,手动快速变得乏味,特别是因为我需要它用于不同尺寸的阵列。

是否有numpy或scipy中的函数或其他可以方便的包?

1 个答案:

答案 0 :(得分:1)

您可能需要查看findiff包。它让您方便地获取任何维度,任何衍生订单和任何所需精度订单的numpy数组的衍生物。项目网站说它的特点是:

  • 区分沿任意轴的任意数量维度的数组
  • 任何所需订单的部分衍生产品
  • 矢量微积分的标准算子,如梯度,发散和卷曲
  • 可以处理均匀和不均匀的网格
  • 可以处理具有常数和变量系数的导数的任意线性组合
  • 可以指定准确性顺序
  • 完全矢量化速度
  • 计算均匀和非均匀网格的任何顺序和精度的原始有限差分系数