np.gradient and matlab gradient

Question

我试图使用numpy.gradient函数在python中重现matlab梯度函数的结果，但我得到的结果与matlab提出的结果非常不同。

不幸的是，我没有再安装一个可以工作的Matlab，所以我无法详细检查它。

给定3个数字(沿三个方向递增(：

dx=1e-7
dy=1e-7
dz=1e-3,

我的matlab代码是：

[F1, F2, F3] = gradient(F, dx,dy,dz)

python中的等价表达式是什么？

我用过：

F1, F2, F3 = np.gradient(F, dx, dy, dz)

但我得到的图表远不能与之相比。

Answer 1

根据Numpy文档，numpy.gradient()使用内部点的中心差计算数值梯度，并允许您使用edge_order选项选择边界处的精度：

numpy.gradient(f, *varargs, axis=None, edge_order=1)

Matlab使用中心差分，其精度高达2阶，因此如果设置edge_order = 2，则应获得与Matlab的输出相同的结果。

[F1, F2, F3] = np.gradient(F, dx, dy, dz, edge_order=2)