gnuplot中使用for循环绘制多函数失败的原因及解决办法

问题原因 #

你用sprintf("f%d(x)",i)生成的是字符串（比如"f1(x)"），gnuplot在plot命令中遇到字符串时，默认会将其视为数据文件名尝试读取，而非解析为函数表达式执行，这就是出现“无法找到文件”警告的核心原因。

解决方法 #

方法1：使用eval解析字符串为函数表达式 #

通过eval命令将拼接出的字符串转换为可执行的函数调用，修改后的代码如下：

set multiplot layout 2,2 
set xrange [-pi:pi] 
set yrange [-1:1] 
f1(x)=sin(x)
f2(x)=cos(x)
f3(x)=sinh(x)
f4(x)=cosh(x)

plot for [i=1:4] eval(sprintf("f%d(x)",i))
unset multiplot

eval会把字符串形式的函数调用转换成gnuplot能识别的表达式，从而正常绘制函数曲线。

方法2：使用函数数组（Gnuplot 5.0及以上版本支持） #

如果你的gnuplot版本在5.0及以上，可以直接定义函数数组，循环调用更简洁直观：

set multiplot layout 2,2 
set xrange [-pi:pi] 
set yrange [-1:1] 
f(1,x)=sin(x)
f(2,x)=cos(x)
f(3,x)=sinh(x)
f(4,x)=cosh(x)

plot for [i=1:4] f(i,x)
unset multiplot

这种方式无需字符串拼接和eval，直接通过数组索引调用对应函数，代码可读性更高。

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