如何在Python中使用有限差分法处理非均匀网格？

Question

我有一个函数

f = dp/de

（压力变化作为能量变化的函数），我想用有限差分计算这个导数。我有一个

p

值列表和一个

e

值列表。然而，

p

与

e

插值形成的网格并不均匀。

我早些时候尝试按如下方式计算 f：

f.append((p[i+1] - p[i])/(e[i+1]- e[i]))

但是，至少在第一点上我需要更好的改进。

如何在Python中对非均匀网格使用有限差分法？基本上网格有下面这种格式（照片）

网格

Answer 1

我有一个名为 uts（不均匀时间序列）的 python 库，它允许计算一阶和二阶的中心差。 重要的代码是here。