不同结果的格点直方图

Question

我在我自己的计算直方图的例程(用IDL编写)和用gnu图计算的直方图之间做了一个非常简单的测试。

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的文件)。

在这里，不同的输出：

​

​

此外，如果我移动hist(x,width)函数，我可以复制正确的IDL直方图。

我对文章进行了修改：使用gnuplot的直方图？，但我无法理解如何告诉GUNplot如何正确地移动内部计算元素的最小和最大值。

Answer 1

计算计数的间隔是不同的。你的中心是0，0.02，0.04，等等。侏儒图形的中心是0.01，0.03，0.05，等等。将bin函数更改为：

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