0

我有三个文件！

coord.xvg
veloc.xvg
force.xvg

这些文件中的每一个都有多个数字的行让我们说 10000 我想构建一个脚本

打开三个文件读取列并对每一行进行算术运算。

例如

if every file has 4 words

coord.xvg >>  Time  x  y  z
veloc.xvg >>  Time vx vy vz
force.xvg >>  Time fx fy fz

c,v,f 代表 coord.xvg, veloc.xvg,force.xvg

if I write the operation 2*v*v+c*f*c the output should be
column1   Column2            Column3           Column4
Time      2*vx*vx+cx*fx*cx   2*vy*vy+cy*fy*cy  2*vz*vz+cz*fz*cz

我在互联网上找到了以下内容

awk '{
       { getline < "coord.xvg" ; if (FNR==90307) for(i=1;i<=2;i+=1) c=$i} 
       { getline < "veloc.xvg" ; if (FNR==90307) for(i=1;i<=2;i+=1) v=$i}
       { getline < "force.xvg" ; if (FNR==90307) for(i=1;i<=2;i+=1) f=$i}
     }
END {print c+v+f}' coord.xvg

它代表我想在 90307 行之后开始阅读的文件。

但它对我没有多大帮助，因为它只返回每个变量的最后一个值

任何想法？？

4

1 回答 1

1

如果我理解正确，可以帮助您入门

$ cat *.xvg
Time 1 2 3
Time 4 5 6
Time 7 8 9
Time 10 11 12
Time 13 14 15
Time 16 17 18
Time 19 20 21
Time 22 23 24
Time 25 26 27

以下 awk 脚本

{ if (FNR>=1) {
    { getline < "coord.xvg" ; c1=$2;c2=$3;c3=$4}
    { getline < "veloc.xvg" ; v1=$2;v2=$3;v3=$4}
    { getline < "force.xvg" ; f1=$2;f2=$3;f3=$4}
    print c1,c2,c3,v1,v2,v3,f1,f2,f3
    print $1, c1+v1+f1, c2+v2+f2, c3+v3+f3
}}

从每个文件中读取一行并将数据放入变量中，如此处所示

$ awk -f s.awk coord.xvg 
1 2 3 19 20 21 10 11 12
Time 30 33 36
4 5 6 22 23 24 13 14 15
Time 39 42 45
7 8 9 25 26 27 16 17 18
Time 48 51 54

该if (FNR>=1)部分控制显示哪些行。计数从 1 开始，根据需要进行更改。我留给你的实际计算:-)

于 2013-09-26T09:44:14.420 回答