计算每列数值平均

shihyu

1. #!/usr/bin/python
   
2. # -*- coding: utf-8 -*-
   
3. import re
   
4. import sys
   
5. 
   
6. def get_nodes(inputline):
   
7.     p = re.compile("g\=\(([-+]?\d+\.\d+)\,\s([-+]?\d+\.\d+),\s([-+]?\d+\.\d+).*a\=\(([-+]?\d+\.\d+),\s([-+]?\d+\.\d+),\s([-+]?\d+\.\d+)")
   
8.     lines = inputline.read()
   
9.     lines = lines.strip().split('\n')
   
10.     lens = len(lines)
    
11. 
    
12.     print(lens)
    
13. 
    
14.     s0 = 0
    
15.     s1 = 0
    
16.     s2 = 0
    
17.     s3 = 0
    
18.     s4 = 0
    
19.     s5 = 0
    
20. 
    
21.     for i in range(0, lens, 1):
    
22.         n1 = p.findall(lines[i])[0]
    
23.         s0 = s0 + float(n1[0])
    
24.         s1 = s1 + float(n1[1])
    
25.         s2 = s2 + float(n1[2])
    
26.         s3 = s3 + float(n1[3])
    
27.         s4 = s4 + float(n1[4])
    
28.         s5 = s5 + float(n1[5])
    
29. 
    
30.     print(s0 / lens)
    
31.     print(s1 / lens)
    
32.     print(s2 / lens)
    
33.     print(s3 / lens)
    
34.     print(s4 / lens)
    
35.     print(s5 / lens)
    
36. 
    
37. if __name__ == "__main__":
    
38.     fin = open(sys.argv[1])
    
39.     get_nodes(fin)
    

复制代码

g=(0.00000, -0.00195, 0.00000), a=(9.13281, -0.03906, 0.56836)
g=(0.00000, -0.00195, 0.00000), a=(9.12305, -0.05859, 0.57813)
g=(0.00000, -0.00195, 0.00000), a=(9.06445, -0.04883, 0.57813)
g=(0.00000, -0.00195, 0.00000), a=(9.04492, -0.02930, 0.57813)
g=(0.00000, -0.00195, 0.00000), a=(9.07422, -0.05859, 0.59766)
g=(0.00000, -0.00195, 0.00000), a=(9.02539, -0.04883, 0.58789)
g=(0.00000, -0.00195, 0.00000), a=(9.05469, -0.06836, 0.52930)
g=(0.00000, -0.00195, 0.00000), a=(9.11328, -0.03906, 0.63672)
g=(0.00000, -0.00195, 0.00000), a=(9.08398, -0.06836, 0.55859)
g=(0.00000, -0.00195, 0.00000), a=(9.06445, -0.06836, 0.57813)
g=(0.00000, -0.00195, 0.00000), a=(9.12305, -0.06836, 0.62695)
g=(0.00000, -0.00195, 0.00000), a=(9.10352, -0.05859, 0.57813)
g=(0.00000, -0.00195, 0.00000), a=(9.08398, -0.07813, 0.57813)



我目前是用 python regular expression 计算出每列数值总和的平均～ 不知道用 sed or awk 一行代码有办法作到同样事情嘛？

谢谢

jason680

回复 1# shihyu

How about this ...

$ awk -F, '{gsub(".=\\(","");for(n=1;n<=NF;++n)a[n]+=$n}END{for(n=1;n<=NF;++n)s=s a[n]/NR " ";print s}' FILE
0 -0.00195 0 9.08398 -0.05634 0.582635

$ cat FILE
g=(0.00000, -0.00195, 0.00000), a=(9.13281, -0.03906, 0.56836)
g=(0.00000, -0.00195, 0.00000), a=(9.12305, -0.05859, 0.57813)
g=(0.00000, -0.00195, 0.00000), a=(9.06445, -0.04883, 0.57813)
g=(0.00000, -0.00195, 0.00000), a=(9.04492, -0.02930, 0.57813)
g=(0.00000, -0.00195, 0.00000), a=(9.07422, -0.05859, 0.59766)
g=(0.00000, -0.00195, 0.00000), a=(9.02539, -0.04883, 0.58789)
g=(0.00000, -0.00195, 0.00000), a=(9.05469, -0.06836, 0.52930)
g=(0.00000, -0.00195, 0.00000), a=(9.11328, -0.03906, 0.63672)
g=(0.00000, -0.00195, 0.00000), a=(9.08398, -0.06836, 0.55859)
g=(0.00000, -0.00195, 0.00000), a=(9.06445, -0.06836, 0.57813)
g=(0.00000, -0.00195, 0.00000), a=(9.12305, -0.06836, 0.62695)
g=(0.00000, -0.00195, 0.00000), a=(9.10352, -0.05859, 0.57813)
g=(0.00000, -0.00195, 0.00000), a=(9.08398, -0.07813, 0.57813)