Gnuplot直方图与框不同的结果

时间:2014-06-14 16:05:33

标签: plot gnuplot histogram binning

我在自己的例程中进行了一个非常简单的测试来计算直方图(用IDL编写)和用gnuplot计算的那个。

clear
reset

width=0.02
set terminal pngcairo size 800,500 enhanced font 'Verdana,14'
set output "GNUvsMY.png"
set boxwidth width
set style fill transparent solid 0.5 border #fillstyle
set grid
set ytics 0.,5,40
set xrange [-0.09:0.09]
set yrange [0:40]
set xlabel "x"
set ylabel "Frequency [N]"
##########
hist(x,width)=width*floor(x/width)+width/2.0
#########
plot "randomIDL.dat" u 2:1 w boxes lc rgb "blue" title 'test3 my hist',\
     "random.dat" u (hist($1,width)):(1.0) smooth freq w boxes lc rgb "red" title 'test3 GNU shift'

为了使测试非常简单,我在每个bin中生成特定数量的对象(bin宽度等于0.02)。所以我现在完全知道每个箱子里有多少物体。

现在我在IDL程序中计算bin,以便计算最小值和最大值内的数字,实际上它会根据需要计算它们。但是当我尝试使用gnuplot的hist函数时,我可以重现我的IDL直方图的(正确)结果。我附加randomIDL.dat,其中您可以在每个bin中找到对象的数量,并在random.dat文件中找到所有对象(这是要传递给GNUplot hist的文件)

这里有不同的输出:

myhistogram

此外,如果我移动hist(x,width)功能,我可以重现正确的IDL直方图。

我点了帖子:Histogram using gnuplot?但我无法理解如何告诉GUNplot如何正确移动内部bin的最小值和最大值。

1 个答案:

答案 0 :(得分:3)

计算计数的间隔不同。您的中心位于0,0.02,0.04,依此类推。 Gnuplot的中心位于0.01,0.03,0.05,依此类推。将bin函数更改为:

hist(x,width)=width*floor(x/width+0.5)