我试图在numpy数值计算数组的二阶梯度.
a = np.sin(np.arange(0, 10, .01))
da = np.gradient(a)
dda = np.gradient(da)
这就是我的想法.应该这样做吗？
我问这个,因为在numpy中没有选项说np.gradient(a,order = 2).我担心这种用法是否错误,这就是为什么numpy没有实现这一点.
PS1:我确实知道有np.diff(a,2).但这只是单边估计,所以我很好奇为什么np.gradient没有类似的关键字.
PS2:np.sin()是玩具数据 - 真实数据没有分析形式.
谢谢!