Loops Gnuplot参数拟合

Loops Gnuplot参数拟合,loops,unix,gnuplot,physics,curve-fitting,Loops,Unix,Gnuplot,Physics,Curve Fitting,我有50个文件,其中包含不同温度下的能量与体积。我有一个脚本,它可以工作,并将使用一些标准多项式对这些文件进行循环完全拟合 F(x)=a+bx+cx2+d*x3 我的能量-体积曲线的导数是-压强,单位为eV/Angstrom,dF(T)/dV=-p(T)。我想找到压力为360GPA的50个温度下的相应体积。以eV/Angstrom为单位的360 GPA为2.246943285 所以我想找到每个文件的x(体积),其中能量体积曲线的梯度 dF(x)=-2.246943285 在哪里 dF(x)=b+

我有50个文件,其中包含不同温度下的能量与体积。我有一个脚本,它可以工作,并将使用一些标准多项式对这些文件进行循环完全拟合

F(x)=a+bx+cx2+d*x3

我的能量-体积曲线的导数是-压强,单位为eV/Angstrom,dF(T)/dV=-p(T)。我想找到压力为360GPA的50个温度下的相应体积。以eV/Angstrom为单位的360 GPA为2.246943285

所以我想找到每个文件的x(体积),其中能量体积曲线的梯度

dF(x)=-2.246943285

在哪里

dF(x)=b+2cx+3dx**2

然后我想将这些值存储到新文件中,例如,100度的Pressure result.100

F(x) = a + b*x + c*x**2 + d*x**3
dF(x)= b + 2*c*x + 3*d*x**2


do for [n=100:5000:100]{
   input_file = sprintf('Gibbs.%.f',n)
   fit F(x) input_file via a,b,c,d

   # output file for print command
   output_file = sprintf('Pressure.%.f.result', n)
   set print output_file

   # Here I want to show the file and it's temperature and the corresponding volume at 360GPA
   print(input_file)
   print(dF(x)=-2.246943285)

}