np.gradient与matlab梯度

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 documentation，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)