适合使用Gnuplot的时间序列

Question

我是Gnuplot的忠实粉丝，现在我想将fit-function用于时间序列。

我的数据集如下：

1.000000 1.000000 0.999795 0.000000 0.000000 0.421927 0.654222 -25.127700 1.000000 1994-08-12
1.000000 2.000000 0.046723 -0.227587 -0.689491 0.328387 1.000000 0.000000 1.000000 1994-08-12 
2.000000 1.000000 0.945762 0.000000 0.000000 0.400038 0.582360 -8.624480 1.000000 1995-04-19 
2.000000 2.000000 0.060228 -0.056367 -0.680224 0.551019 1.000000 0.000000 1.000000 1995-04-19
3.000000 1.000000 1.016430 0.000000 0.000000 0.574478 0.489638 -3.286880 1.000000 1995-07-15

我的拟合脚本：

set timefmt "%Y-%m-%d"
set xdata time
set format x "%Y-%m-%d"
f(x)=a+b*x
fit f(x) "model_fit.dat" u 10:($2==2?$4:1/0) via a,b

所以我对时间数据进行了条件拟合。 我的问题是，Gnuplot拟合函数不能处理时间数据。 我在这里发现了一个类似的问题：Linear regression for time series with Gnuplot但我不想使用其他软件。而且我也不知道如何将时间值更改为数字，然后再返回....

任何人都可以帮我解决这个问题吗？

非常感谢！

Answer 1

确实，gnuplot的拟合机制适用于时间数据。你必须只关注一些细节。

通常，可以精确地求解通过两个数据点的线性拟合。但是gnuplot通常是非线性拟合，所以初始值很重要。

用于该行的有效x值以秒为单位给出。使用初始值b = 1e-8可以正常工作：

set timefmt "%Y-%m-%d"
set xdata time
set format x "%Y-%m-%d"
f(x)=a+b*x
a = 1
b = 1e-8
fit f(x) "model_fit.dat" u 10:($2==2?$4:1/0) via a,b

plot "model_fit.dat" u 10:($2==2?$4:1/0) lt 1 pt 7 title 'data',\
     f(x) w l lt 1 title 'fit'