我试图用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)但我得到的图表远不是可比较的。
发布于 2020-10-15 21:45:34
对于每个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)
https://stackoverflow.com/questions/64370598
复制相似问题