GNUPLOT:从平滑累积中保存数据

时间:2014-06-08 12:02:24

标签: plot statistics gnuplot cumulative-sum

我制作了这个简单的累积和直方图,显示了实数的均匀随机分布(n = 1000):

http://www.filedropper.com/random1_1:random1.dat

宏是:

unset key
clear
reset


n=120 #number of intervals
max=4. #max value
min=1. #min value

width=(max-min)/n #interval width
#function used to map a value to the intervals
bin(x,width)=width*floor(x/width)+width/2.0 # cosi viene centrato in mezzo
set xtics min,(max-min)/10,max
set boxwidth width
set style fill solid 0.5 border


set ylabel 'Frequency'
set y2label 'Cumulative frequency'
set y2tics 0,100,1000
set ytics nomirror

set xrange [0.9:4.1]
set yrange [0:25]

set terminal pngcairo size 800,500 enhanced font 'Verdana,14'
set output "testCum.png"

plot 'random1.dat' using (bin($1,width)):(1.0) smooth frequency with boxes title 'histogram',\
'' using (bin($1,width)):(1.0) smooth cumulative axis x1y2 w l lt 2 lw 2 lc rgb 'green' title 'cumul'

现在output.png是:

enter image description here

我如何告诉Gnuplot我不仅要记录累积的情节,还要记录保存在特定file.dat中的数字?

1 个答案:

答案 0 :(得分:4)

使用smooth应用set table ...后,您可以保存数据。在最简单的情况下,如果您只需要累积数据,只需使用:

set table 'random1-smoothed.dat'
plot 'random1.dat' using (bin($1,width)):(1.0) smooth cumulative
unset table

为了更好地包含在您的脚本中,您还可以将整个现有plot命令包装在set table中:

...
set table 'random1-smoothed.dat'
plot 'random1.dat' using (bin($1,width)):(1.0) smooth frequency with boxes title 'histogram',\
'' using (bin($1,width)):(1.0) smooth cumulative axis x1y2 w l lt 2 lw 2 lc rgb 'green' title 'cumul'
unset table

set terminal pngcairo size 800,500 enhanced font 'Verdana,14'
set output "testCum.png"
replot