如何在gnuplot中对不同文件带行偏移的列做算术运算并绘图？

1. 前置准备 #

• 确认两个CSV文件格式：列分隔符为逗号（若为其他分隔符，后续脚本需调整），数据行无空行。
• 在VS Code中安装Gnuplot Integration扩展，方便直接运行脚本。

2. 编写gnuplot脚本 #

新建后缀为.gp的文件（比如plot_data.gp），写入以下内容，根据实际路径修改文件名：

# 设置CSV分隔符（若为制表符，改为"\t"）
set datafile separator comma

# 可选：设置绘图样式与标签
set style data linespoints  # 采用带点折线图
set title "File1第2列 / File2第2列（对应行）"
set xlabel "File1第1列数据"
set ylabel "计算结果"

# 步骤1：将file2第10-60行的第2列数据存入数组（共50行，对应file1的1-50行）
array file2_col2[50]
stats 'file2.csv' using (file2_col2[int($0)-9+1] = $2) every ::9::59 nooutput
# 说明：every ::9::59 取行索引9到59（对应实际第10到60行），gnuplot行索引从0开始，数组索引从1开始

# 步骤2：读取file1第1-50行，计算y轴值并绘图
plot 'file1.csv' every ::0::49 using 1:($2 / file2_col2[int($0)+1]) title "计算曲线"
# 说明：every ::0::49 取file1前50行（行索引0到49），$1为x轴（第1列），$2/file2_col2[...]为y轴计算值

关键参数说明 #

• every ::a::b：控制读取的行范围，a和b为行索引（从0开始），比如实际第1行对应索引0，第10行对应索引9。
• stats ... nooutput：仅提取数据到数组，不输出冗余统计信息。
• 若CSV含表头：需调整行索引范围。比如file1第1行是表头，要取第2-51行数据，将every ::0::49改为every ::1::50；file2同理，根据表头后的实际数据行调整every参数。

3. 在VS Code中运行脚本 #

• 打开.gp脚本文件。
• 右键点击编辑器区域，选择Run Gnuplot Script，或使用扩展快捷键（通常为Ctrl+Shift+G，以扩展说明为准）。
• gnuplot会自动弹出绘图窗口展示结果。

4. 常见问题排查 #

• 报错"array index out of range"：检查数组大小是否与读取行数一致（此处为50），以及行索引计算是否正确。
• 数据读取错误：确认set datafile separator设置的分隔符与CSV实际一致，比如分号分隔需改为set datafile separator ";"。
• 文件找不到：确保脚本中文件路径正确，Windows下可使用绝对路径，如'C:/Users/xxx/Documents/file1.csv'（正斜杠）或'C:\\Users\\xxx\\Documents\\file1.csv'（双反斜杠）。

内容的提问来源于stack exchange，提问作者ASD