GNUPLOT适合循环

Question

我有一个包含N列的文本文件，我想使用GNUPLOT使用custum f（x）的命令拟合来拟合每列的数据，例如f（x）= exp（x / t），通过t拟合。文本文件的第一列包含x的值。 我无法在列上进行循环。 任何人都可以帮助我吗？

非常感谢， 弗朗西斯

Answer 1

也许这不是最好的代码，但它有效。

    #Given a txt file, where the first colum is x, this script calculates
#the best fit for each column starting from column 2
#N=number of column
#(R0:R1)=range of x values over which calculate fit
#The script saves parameters fit in file parameters.txt.
#Remember to erase this file before relaunching because appends results.

N=28
R0=0
R1=7

do for [i=2:N:1] {
   set terminal postscript enhanced color solid eps
   set output graph(i)
   set xrange[0:20]

   f(x) = exp(-x/tau)
   set fit errorvariables
   fit [R0:R1] f(x) 'myfilename.txt' u 1:i via tau

   set print "parameters.txt" append
   print tau,tau_err
   plot"myfilename.txt" u 1:i, f(x)
   reset  
}