Gnuplot:如何为不同单位添加y2轴刻度

时间:2012-12-20 10:00:00

标签: gnuplot

我正在绘制文件中的数据。数据点以公制单位表示。我想在右边显示第二个刻度(y2),它是标准单位。

该文件代表火箭发动机的推力。数据以牛顿为单位。我想在左边显示牛顿(这自然会自然发生),右边显示磅力。转换是一个简单的因素(将N乘以0.2248以获得lbf)。

我可以设置y2tics,如果我手动设置y2range,它们会出现在右侧。我不知道怎么做是自动将y2range设置为y1range *因子。

我最终的解决方案是绘制两次,在y1上以牛顿为单位,在y2上以磅为单位绘制一次,并使y2图几乎不可见:

plot '-' using 1:($2*0.2248) with dots axes x1y2 lc rgb 'white' notitle, \
      '' using 1:2 with lines lc rgb '<color>' title '<title>'

上面的解决方案通常会产生稍微不同的y刻度:使用autoragne,gnuplot会将范围向上舍入,因此每个轴上的顶部刻度是一个圆数,当然对于不同的单位,舍入是不同的。

最终我最终得到了每个图中找到最高推力值的Python代码,然后我明确地将yrange设置为该数字,并将y2range设置为该数字* 0.2248:

f.write("set yrange [0:%s]; set y2range[0:%s]\n" % (peak_thrust, peak_thrust*NEWTON_LBF));

以下是最终结果:http://www.lib.aero/hosted/motors/cesaroni_12-15-12.html(下面的示例图表) enter image description here

3 个答案:

答案 0 :(得分:16)

在我看来,最简单的方法是简单地缩放数据:

set y2tics
plot sin(x) w lines, 5*sin(x) w lines axes x1y2

当然,您正在绘制文件中的数据,因此它看起来更像是:

set y2tics
FACTOR=0.2248  #conversion factor from newtons to lbf
plot 'datafile' u 1:2 w lines, '' u 1:(FACTOR*$2) w lines

如果您明确设置了yrange(您可能需要这样做):

set yrange [ymin:ymax]
set y2range [ymin*FACTOR:ymax*FACTOR]

最后,如果你真的想依靠自动缩放,你需要做一些“体操”。

首先,设置一个虚拟终端,这样我们就可以plot而无需制作情节:

set term unknown
plot 'datafile' u 1:2  #collect information on our data

现在我们收集了有关数据的信息,我们可以设置真实的y2range

FACTOR=0.2248
set y2range [FACTOR*GPVAL_Y_MIN : FACTOR*GPVAL_Y_MAX]
set y2tics nomirror
set ytics nomirror

现在设置终端并绘制数据:

set term ...
set output ...
plot 'datafile' u 1:2 w lines

答案 1 :(得分:11)

5.0版增加了对yy2(或xx2)轴之间此类关系的支持:

set xrange[0:370]
set ytics nomirror
set y2tics
set link y2 via 0.2248*y inverse y/0.2248
plot x

example for <code>set link</code>

答案 2 :(得分:2)

我知道这是一个老问题,答案已被接受,但我认为值得分享我的方法。

我只是为x2axis使用修改后的标签。在你的情况下,这将是

set y2tics ("10" 10/0.2248, "20" 20/0.2248 etc etc...

可以这样循环

do for [i=0:1000:10] { set y2tics add (sprintf("%i",i) i/0.2248) }

应根据您的数据调整for范围(您可以使用stats和变量GPVAL_DATA_Y_MAX完全放心)。 别忘了

set ytics nomirror

这将准确地给出你在寻找什么,(几乎)一个班轮: enter image description here

如果你想使用一个网格,并在x2axis上有转换因子,那么例如标签y = 50 N将对应y2 = 11.2(如果使用网格,它会保持整洁)你可以做到

do for [i=0:1000:50] { set y2tics add (sprintf("%5.1f",i*0.2248) }

这是结果: enter image description here