gnuplot：for循环中数据文件的可变路径

如何解决gnuplot：for循环中数据文件的可变路径

我想使用for循环在同一张图上绘制多条曲线。每个数据文件（名为stat_coupe）位于不同的文件夹（fwal055wal055 / rep16 /和fwal055wal055_c2 / rep20 /）中。 fwal055wal055和fwal055wal055_c2对应于仿真的名称。首先，我需要在其他文件（名为file_fwal055wal055和file_fwal055wal055_c2）中得到一个先前的结果，即一个数字（Utau）。借助命令awk，成功完成了此操作。结果取决于文件：Utaufwal055wal055 = 10.5和Utaufwal055wal055_c2 = 12.2。

然后，我需要将对应于路径fwal055wal055 / rep16 /的文件stat_coupe的第一列除以Utaufwal055wal055的值，并对对应于路径fwal055wal055_c2 / rep20 /的文件stat_coupe进行相同的处理，并将其值的Utaufwal055wal055_c2 。此外，每个图应具有特定的格式，具体取决于模拟运行的类型（fwal055wal055或fwal055wal055_c2）。

提出的问题简化为2个模拟fwal055wal055和fwal055wal055_c2和1个图，但是我要绘制约20个模拟和15个各种图形，这就是为什么我要使用for循环的原因。

总结一下，我每次迭代：

• 特定格式
• 特定路径
• Utau的特定值

我想在for循环的每次迭代中指出Utau的wright格式，路径和值。我在下面提出的解决方案成功地允许为每个模拟获取Utau的值，但是代码@ path_.i和@ format_.i不起作用。

#!/bin/bash
for elem in fwal055wal055 fwal055wal055_c2;
do
  Utau[${elem}]=$(awk 'FNR==5{print $1}' file_$elem)
done

gnuplot -persist  <<-EOFMarker
format_fwal055wal055='pt 1 ps 1.0 lc 0 title "WALE"'
format_fwal055wal055_c2='pt 2 ps 1.0 lc 0 title "WALE c2"'

path_fwal055wal055='"fwal055wal055/rep16/stat_coupe"'
path_fwal055wal055_c2='"fwal055wal055_c2/rep20/stat_coupe"'

list="fwal055wal055 fwal055wal055_c2"
plot for [i in list] @path_.i u 1:(\$2/${Utau[${i}]}) @format_.i
EOFMarker

我想获得等同于以下内容的东西：

plot @path_fwal055wal055 u 1:(\$2/${Utau[${i}]}) @format_fwal055wal055,\
@path_fwal055wal055_c2 u 1:(\$2/${Utau[${i}]}) @format_fwal055wal055_c2

有人可以帮我解决这个问题吗？

非常感谢你， 马丁

解决方法

选中help sprintf，help words和help word。 我将创建两个具有相同数量项的字符串，然后将它们与sprintf()组合在一起。从gnuplot 5.2开始，您还可以使用数组。

# Version 1
PATHS = '"fwal055wal055/rep16/stat_coupe" "fwal055wal055_c2/rep20/stat_coupe"'
FILES = "fwal055wal055 fwal055wal055_c2"

plot for [i=1:words(FILES)] sprintf("%s_%s",word(PATHS,i),word(FILES,i)) u 1:2 

或者您可以为文件名定义一个函数，以使plot命令简短易读。

# Version 2
PATHS = '"rep16/stat_coupe" "rep20/stat_coupe"'
FILES = "fwal055wal055 fwal055wal055_c2"

myFilename(i) = sprintf("%s/%s_%s",i))

plot for [i=1:words(FILES)] myFilename(i) u 1:2 

添加（经过一些澄清...）

如果我现在能正确理解您的问题，则以下代码可以完成这项工作。 为了提取UTAUS，您需要在绘制之前进行单独的循环，然后将提取的值存储在字符串中。在绘制期间，您可以通过word(UTAUS,i)返回这些值。由于您执行了数学运算column(2)/word(UTAUS,i)，因此gnuplot会将其解释为数字。选中help words，help word，help sprintf，help every。

代码：

### extract and normalize in a loop with individual files and directories
reset session

FILES  = 'fwal055wal055 fwal055wal055_c2'
DIRS   = 'rep16 rep20'
TITLES = '"WALE" "WALE c2"'   # if you have spaces you need to put it into double quotes
UTAUS   = ''

# define functions for better readability
myExtractionFile(i) = sprintf("file_%s",i))
myDataFile(i) = sprintf("%s/%s/stat_coupe",word(DIRS,i))
myTitle(i) = word(TITLES,i)

# define point or line appearance. Add more if you have more files
set style line 1 pt 1 ps 1.0 lc 0
set style line 2 pt 2 ps 1.0 lc 1

# extract the UTAUs
do for [i=1:words(FILES)] {
    set table $Dummy
        plot myExtractionFile(i) u (utau=$1) every ::4::4 w table   # extract value row 5,column 1 (not counting header lines)
    unset table
    UTAUS = UTAUS.sprintf(" %g",utau)   # append the extracted value as string
}

plot for [i=1:words(FILES)] myDataFile(i) u 1:(column(2)/word(UTAUS,i)) ls i title myTitle(i)
### end of code