导数求解

[provance]

各位好，我是一个Matlab初学者。现在需要将一列数据 V 对 z 座标求二阶偏导数Vzz，我是根据导数定义进行求解的，即Vz ＝ diff(V) ./ diff(z), Vzz = diff(Vz) ./ diff(z)。
结果如附件中的Vzz所示，非常的不平滑；原因可能是一阶导数 Vz 已经不是光滑的函数，如附件中Vz所示。或许在求解此类数值微分时不宜采用定义的方法，但目前不知道有其它的方法可以求解，不知各位能不能提供一定建议？感激不尽。