我正在尝试在 gnuplot 中拟合一些数据。数据是不同的曲线。我有三个.gnp文件：一个用于数学过程，一个用于调用文件进行分析，一个必须在 gnuplot 中启动。
我得到的是这个错误，但我不知道如何解决它。
gnuplot> load '0-one_fit.gnp'
gnuplot> id = $0
^
"0-one_fit.gnp" 第 5 行：从无效上下文调用的 column((
此文件应存储从 0 到变量 id 的序列号。序列号通过命令调用从外部 loop.gnp 传递给此脚本。
文件说：
id = $0
filename = sprintf("FeCoBSiNbCu%04d.chi",id)
fit [26:35.6] f110(x) filename u 1:2 via x110,A110,w110,a0,a1
fit [42.3:45.1] f200(x) filename u 1:2 via x200,A200,w200,b0,b1
fit [48:56.1] f23(x) filename u 1:2 via x2,A2,w2,x211,A211,w211,c0,c1
entry1 = sprintf("%dt%.3ft%.3ft%.1ft%.1ft%.3ft%.3ft",id,x110,x110_err,A110,A110_err,w110,w110_err)
entry2 = sprintf("%.3ft%.3ft%.1ft%.1ft%.3ft%.3ft",x200,x200_err,A200,A200_err,w200,w200_err)
entry3 = sprintf("%.3ft%.3ft%.1ft%.1ft%.3ft%.3ft",x2,x2_err,A2,A2_err,w2,w2_err)
entry4 = sprintf("%.3ft%.3ft%.1ft%.1ft%.3ft%.3f",x211,x211_err,A211,A211_err,w211,w211_err)
print entry1, entry2, entry3, entry4
更新： 我用id=ARG1更改了文件，但现在它给了我另一个错误
gnuplot> call '0-one_fit.gnp' 0
"0-one_fit.gnp" line 8: f_sprintf: attempt to print string value with numeric format
如果我现在用int(id)更改，它说：
load '0-one_fit.gnp'
"0-one_fit.gnp" line 8: Non-numeric string found where a numeric expression was expected
更新我尝试首先运行名为0-in-situ.gnp的文件，其中包含峰值拟合函数：
#### First peak: (110)
bg110(x) = a0 + a1*x
k110(x) = (x-x110)/w110
Gauss110(x) = A110*exp(-log(2)*k110(x)**2)
Cauchy110(x) = A110/(1 + k110(x)**2)
pV110(x) = alpha*Gauss110(x) + (1-alpha)*Cauchy110(x)
f110(x) = pV110(x) + bg110(x)
#### Second peak: (200)
bg200(x) = b0 + b1*x
k200(x) = (x-x200)/w200
Gauss200(x) = A200*exp(-log(2)*k200(x)**2)
Cauchy200(x) = A200/(1 + k200(x)**2)
pV200(x) = beta*Gauss200(x) + (1-beta)*Cauchy200(x)
f200(x) = pV200(x) + bg200(x)
#### Third peak & deconvolution (Q2 & 211):
bg(x) = c0 + c1*x
k2(x) = (x-x2)/w2
k211(x) = (x-x211)/w211
Gauss2(x) = A2*exp(-log(2)*k2(x)**2)
Cauchy2(x) = A2/(1 + k2(x)**2)
pV2(x) = gamma*Gauss2(x) + (1-gamma)*Cauchy2(x)
Gauss211(x) = A211*exp(-log(2)*k211(x)**2)
Cauchy211(x) = A211/(1 + k211(x)**2)
pV211(x) = delta*Gauss211(x) + (1-delta)*Cauchy211(x)
f23(x) = pV2(x) + pV211(x) + bg(x)
alpha = 0
beta = 0
gamma = 1
delta = 0
A110 = 15574
x110 = 31.1
w110 = 2.4
a0 = 1760
a1 = 6
A200 = 4153
x200 = 52.2
w200 = 3
b0 = 1500
b1 = 15
A2 = 1800
x2 = 51.6
w2 = 3
A211 = 1897
x211 = 76.6
w211 = 0.3
c0 = 1830
c1 = -2
#### Loop through all XRD patterns:
load "0-loop.gnp"
程序运行，但是当我打开应该包含数据的文件.dat时，显然它不会分析所有曲线，而只会分析第一条曲线。这是下面循环的问题吗？
call "0-one_fit.gnp" 447
call "0-one_fit.gnp" 448
call "0-one_fit.gnp" 449
call "0-one_fit.gnp" 450
call "0-one_fit.gnp" 451
call "0-one_fit.gnp" 452
您可以在 plot 命令或 plot 命令中使用的函数中使用$0或column(0)。在其他地方，我想这是"无效上下文"。 检查以下示例：
还要检查help using、help column、help pseudocolumns。
reset session
$Data <<EOD
11 12 13 14 15
21 22 23 24 25
31 32 33 34 35
41 42 43 44 45
EOD
# id = $0 # does not work, invalid context
# id = column(0) # does not work, invalid context
id(n) = column(n) # within a function, works
plot $Data u 0:1 w lp # works
plot $Data u ($0):1 w lp # works
plot $Data u (column(0)):1 w lp # works
plot $Data u (id(0)):1 w lp # works
编辑：
我想，我误解了你的问题。我以为您想为变量分配列值。
如果要将参数传递给 gnuplot 例程，则可以call带有参数的 gnuplot 文件(而不是load(。
call '0-one_fit.gnp' 123
在0-one_fit.gnp内，您可以通过id = ARG1访问此值。 在早期的 gnuplot 版本("旧式"(中，这是id=$0的。
检查help call和help call old-style(gnuplot 5.2(