计算场的时间导数

Question

我在三个不同时间有三个维度为 (4x4) 的字段：

varT1 = np.array([[1,6,2,5],[9,4,2,1], [8,2,9,5], [9,3,2,1]])
varT2 = np.array([[9,4,1,6],[9,8,6,9], [8,4,1,1], [7,2,4,4]])
varT3 = np.array([[1,3,4,9],[8,7,7,4], [7,3,1,3], [3,7,0,6]])

是否可以使用

np.gradient

计算时间导数？

Answer 1

我将沿着一个新轴堆叠三个字段，我们可以将其视为时间，然后使用

np.gradient

沿该轴进行计算

var_stack = np.stack([varT1, varT2, varT3], axis=0)
temporal_derivative = np.gradient(var_stack, axis=0)