在 gnuplot 脚本中使用 eval 和 awk 表达式

Question

我正在尝试在 gnuplot 脚本中绘制参数表达式，其系数存储在文本文件的最后一行。为此，首先我尝试了这个：

plot "<awk 'END{print $1"*cos(t)*cos("$2")-"$3"*sin(t)*sin("$4"), "$1"*cos(t)*sin("$2")+"$3"*sin(t)*cos("$4")"}' manip_file.csv"

但 gnuplot 说 undefined variable: t。所以接下来我尝试了以下方法：

plotCMD = 'awk 'END{print "plot " $1"*cos(t)*cos("$2")-"$3"*sin(t)*sin("$4"), "$1"*cos(t)*sin("$2")+"$3"*sin(t)*cos("$4")"}' manip_file.csv'
eval(plotCMD)

但是这次 gnuplot 说的是';' expected。如果我在命令行中运行awk 命令，它会给我一个正确的方程，gnuplot 绘制它没有问题。因此，丢失一些单/双引号不是问题。试图避开美元符号 (\$1) 也没有解决问题。有什么想法吗？

Answer 1

您将 gnuplot 和 awk 语法完全混合在一起。这是使用 gnuplot 的解决方案：

# 1. Get the last row of your text file
p = system('tail -1 manip_file.csv')

# 2. Get the four parameters, assuming they are white-space separated
p1 = word(p, 1)*1.0
p2 = word(p, 2)*1.0
p3 = word(p, 3)*1.0
p4 = word(p, 4)*1.0

# 3. Define the functions
x(t) = p1*cos(t)*cos(p2) - p3 * sin(t)*sin(p4)
y(t) = p1*cos(t)*sin(p2)+p3*sin(t)*cos(p4)

# 4. Plot them
set parametric
plot x(t), y(t)