有没有办法在gnuplot中绘制平均值?
假设我有一个包含两列数据的文件有没有办法在gnuplot中绘制平均值?,gnuplot,numerical-analysis,Gnuplot,Numerical Analysis,假设我有一个包含两列数据的文件file.dat。 我通常会用 plot "file.dat" u 1:2 我想平均10(例如)个以上的点和10个以下的点,然后在同一个图上绘制。 我可以使用一些外部脚本轻松地完成这项工作,其中我制作了另一个专栏: for(i=-10;i<=10;++i) $3[j] += $2[j-i] for(i=-10;i这可能令人惊讶地不是内置于gnuplot中的。由于gnuplot是如何将数据作为流处理的,因此没有好的方法来操作gnuplot中的各个数据点,
file.dat
。
我通常会用
plot "file.dat" u 1:2
我想平均10(例如)个以上的点和10个以下的点,然后在同一个图上绘制。
我可以使用一些外部脚本轻松地完成这项工作,其中我制作了另一个专栏:
for(i=-10;i<=10;++i)
$3[j] += $2[j-i]
for(i=-10;i这可能令人惊讶地不是内置于gnuplot中的。由于gnuplot是如何将数据作为流处理的,因此没有好的方法来操作gnuplot中的各个数据点,也没有数据点的范围
gnuplot的一个优点是调用外部脚本和工具非常容易。如果您想使用外部脚本来处理来自gnuplot的数据,可以这样做:
plot "<script.py data.dat" u 1:2
看起来Bézier曲线是内置的,可以用作plot”file.dat“u 1:2 s B
。有某种平均值是很好的,但是自定义的会更好,所以问题仍然悬而未决。
#!/usr/bin/python2.7
import sys
if (len(sys.argv) > 6):
print ""
print "This script takes one mandatory argument, the name of a file containing"
print "data to be plotted. It takes up to four optional arguments as follows:"
print " 1) the number of points before a data point to add into average."
print " 2) the number of points after a data point to add into average."
print " 3) the column number of y data (first column is column 1)"
print " 4) the column number of x data (first column is column 1)"
print ""
exit()
# set variable defaults
box_back = 10 # number of points before current point to add into average
box_front = 10 # number of points after current point to add into average
y_col = 2 # column number of y data (first column is column 1)
x_col = 1 # column number of x data (first column is column 1)
# assign variables from command line arguments
inputFileName = str(sys.argv[1])
if (len(sys.argv) > 2):
box_back = int(sys.argv[2])
if (len(sys.argv) > 3):
box_front = int(sys.argv[3])
if (len(sys.argv) > 4):
y_col = int(sys.argv[4])
if (len(sys.argv) > 5):
x_col = int(sys.argv[5])
# open input file
f = open(inputFileName)
# make list from lines in file
lines = f.readlines()
# make sure boxcar average will work
if ((box_back + box_front + 1) > len(lines)):
print ""
print "ERROR: too many points for boxcar averaging."
print ""
exit()
# this is the number of points encompassed in the boxcar average
num_points = box_back + box_front + 1
# this variable is the running sum.
sum_vals = 0
# add up values for first boxcar average
for i_ in range(0,num_points):
sum_vals += float(lines[i_].split()[y_col-1])
print float(lines[box_back].split()[x_col-1]),sum_vals/num_points
# each subsequent average differs only in the first and last points from the
# previous average.
for i_ in range(box_back+1,len(lines)-box_front):
sum_vals += float(lines[i_+box_front].split()[y_col-1])
sum_vals -= float(lines[i_-box_back-1].split()[y_col-1])
print float(lines[i_].split()[x_col-1]),sum_vals/num_points