Question

我有一个长3列的数值数据表。例如，它的一小部分看起来像这样：

-5.3986816409999996e+00 8.3394692357616169e-23 2.2891221151890116e-23
-5.3984375003749996e+00 8.1834931317429596e-23 2.7962314960022732e-23
-5.3981933597499996e+00 7.9967086078053667e-23 3.2928354334362533e-23
-5.3979492191249996e+00 7.7798165392355974e-23 3.7770639137136620e-23
-5.3977050784999996e+00 7.5336312057849817e-23 4.2470933118491389e-23
-5.3974609378749996e+00 7.2590772420102232e-23 4.7011532627697938e-23
-5.3972167972499996e+00 6.9571861716728761e-23 5.1375333330996674e-23
-5.3969726566249996e+00 6.6290925391538841e-23 5.5545894684658167e-23
-5.3967285159999996e+00 6.2760296523720059e-23 5.9507501919987393e-23
-5.3964843753749996e+00 5.8993249531280359e-23 6.3245225306694521e-23
-5.3962402347499996e+00 5.5003950322868833e-23 6.6744976470308050e-23

要清楚，第一列显示时间（t），另外两列表示一些其他变量作为时间的函数，比如说f（t）和g（t）。我想计算每个函数的时间导数并将其保存在新的数据文件中。我不确定我该如何处理这个数字。任何的想法？

Answer 1

一个简单的解决方案。使用numpy包，numpy.diff为您提供衍生产品。

Answer 2

如果你只是在寻找一阶估计，你可以使用δx/δt≈Δx/Δt。

# import file
t = []; f = []; g = []
for line in open('in_file.dat'):
    line = line.split()
    t.append(float(line[0]))
    f.append(float(line[1]))
    g.append(float(line[2]))

# calculate approximate derivatives
dt = [] # times for derivative values
df = [] # approximate derivative of f
dg = [] # approximate derivative of g
for i in range(len(t) - 1):
    delta_t = t[i + 1] - t[i]  # time between points
    dt.append(t[i] + delta_t * 0.5)  # midway between points
    df.append((f[i - 1] - f[i]) / delta_t)
    dg.append((g[i - 1] - g[i]) / delta_t)

# save to new file
f = '{:23.16e} {:23.16e} {:23.16e}\n'
out_lines = [f.format(dt[i], df[i], dg[i]) for i in range(len(dt))]
open('out_file.dat', 'w').writelines(out_lines)

请注意，衍生工具在点之间的中间有时间指数。因此，新的时间和衍生物列表将比输入短1个元素。这是正常的，虽然一些更高级的方法将具有1对1次。另请注意，输出文件格式由您决定。我刚刚使用了与上面大致相同的格式。

Answer 3

for i in range(len(t) - 1):
    delta_t = t[i + 1] - t[i]  # time between points
    dt.append(t[i] + delta_t * 0.5)  # midway between points
    df.append((f[i + 1] - f[i]) / delta_t)
    dg.append((g[i + 1] - g[i]) / delta_t)

从 - 到+ ==＆gt;的校正f [i + 1] - f [i]）/ delta_t

.dat文件中的数字时间导数

3 个答案: