numpy中的高阶有限差分

Question

我已经采样了2D和3D numpy阵列上的函数，我需要一种从这些数组中获取偏导数的方法。我可以手动编码有限差分格式，但我需要的不仅仅是2阶精度，可能是4阶甚至6阶。随着更高的精度顺序编码，它很快就变得乏味，特别是因为我需要它的数组的不同尺寸。

在numpy或scipy或其他包中是否有功能可以方便地做到这一点？

Answer 1

您可能想看看终距包。让我们方便地获取任意维数的numpy数组的导数，任何导数顺序，以及任何期望的精度顺序。该项目的网站说，它的特点是：

• 沿任意轴区分任意维数的数组
• 任意阶偏导数
• 向量演算中的标准算子，如梯度、散度和卷曲
• 能处理均匀和非均匀网格。
• 可以处理常系数和变系数导数的任意线性组合。
• 可以指定精度顺序。
• 完全矢量化的速度
• 计算均匀和非均匀网格的任意阶和精度的原始有限差分系数。