两个文件，对应位置数据求和

vcuteym

问题描述：
第一步，file1和file2的前12行一样，将file1和file2合并成file3，前12行保持不变。
第二步，file1和file2的第13行至最后一行，对应数据求和，合成file3。例如file1的第13行是  1.4  2.3  2.1  ....，file2第13行是 3.3 5.2  5....。合并成file3后，第13行是 4.7  7.5  7.1  .....

file1和file2样例如下，file1和file2有很多行，没有列全
=========file1================
VBM
   1.00000000000000
     3.844636    0.000000    0.000000
     1.922318    3.329553    0.000000
     1.922318    1.109851    3.139132
   Si
     2
Direct
  0.000000  0.000000  0.000000
  0.250000  0.250000  0.250000

   56   56   56
-2.4237     -2.2477     -1.7396     -.95629     0.14895E-01  1.0694      2.1020      3.0233      3.7722      4.3190
  4.6626      4.8221      4.8280      4.7146      4.5153      4.2596      3.9669      3.6501      3.3211      2.9918
  2.6725      2.3724      2.0984      1.8561      1.6499      1.4840      1.3620      1.2872      1.2620      1.2872
  1.3620      1.4840      1.6499      1.8561      2.0984      2.3724      2.6725      2.9918      3.3211      3.6501
  3.9669      4.2596      4.5153      4.7146      4.8280      4.8221      4.6626      4.3190      3.7722      3.0233
  2.1020      1.0694     0.14895E-01 -.95629     -1.7396     -2.2477     -2.2477     -1.8464     -1.0947     -.86629E-01
  1.0578      2.2129      3.2688      4.1468      4.8054      5.2377      5.4619      5.5102      5.4199      5.2271
  4.9641      4.6555      4.3153      3.9556      3.5886      3.2259      2.8783      2.5545      2.2620      2.0066
  1.7938      1.6284      1.5150      1.4573      1.4573      1.5150      1.6284      1.7938      2.0066      2.2620
  2.5545      2.8783      3.2259      3.5886      3.9556      4.3153      4.6555      4.9641      5.2271      5.4199
  5.5102      5.4619      5.2377      4.8054      4.1468      3.2688      2.2129      1.0578     -.86629E-01 -1.0947

==========================================================================

========file2===================
VBM
   1.00000000000000
     3.844636    0.000000    0.000000
     1.922318    3.329553    0.000000
     1.922318    1.109851    3.139132
   Si
     2
Direct
  0.000000  0.000000  0.000000
  0.250000  0.250000  0.250000

   56   56   56
-2.4237     -2.2477     -1.7396     -.95629     0.14895E-01  1.0694      2.1020      3.0233      3.7722      4.3190
  4.6626      4.8221      4.8280      4.7146      4.5153      4.2596      3.9669      3.6501      3.3211      2.9918
  2.6725      2.3724      2.0984      1.8561      1.6499      1.4840      1.3620      1.2872      1.2620      1.2872
  1.3620      1.4840      1.6499      1.8561      2.0984      2.3724      2.6725      2.9918      3.3211      3.6501
  3.9669      4.2596      4.5153      4.7146      4.8280      4.8221      4.6626      4.3190      3.7722      3.0233
  2.1020      1.0694     0.14895E-01 -.95629     -1.7396     -2.2477     -2.2477     -1.8464     -1.0947     -.86629E-01
  1.0578      2.2129      3.2688      4.1468      4.8054      5.2377      5.4619      5.5102      5.4199      5.2271
  4.9641      4.6555      4.3153      3.9556      3.5886      3.2259      2.8783      2.5545      2.2620      2.0066
  1.7938      1.6284      1.5150      1.4573      1.4573      1.5150      1.6284      1.7938      2.0066      2.2620
  2.5545      2.8783      3.2259      3.5886      3.9556      4.3153      4.6555      4.9641      5.2271      5.4199
  5.5102      5.4619      5.2377      4.8054      4.1468      3.2688      2.2129      1.0578     -.86629E-01 -1.0947

chengchow

1. [root@AnsibleC14 tmp]# cat aa.py
   
2. #!/usr/bin/python
   
3. 
   
4. import os
   
5. import sys
   
6. 
   
7. File="sed -n '1,12p' 1.txt"
   
8. File1="sed -n '13,$p' 1.txt"
   
9. File2="sed -n '13,$p' 1.txt"
   
10. 
    
11. String=""
    
12. List1=os.popen(File1).read().strip('\r\n').split()
    
13. List2=os.popen(File2).read().strip('\r\n').split()
    
14. 
    
15. print os.popen(File).read().strip('\r\n')
    
16. 
    
17. if len(List1)==len(List2) :
    
18.     i=0
    
19.     while i<len(List1) :
    
20.         x=float(List1[i])+float(List2[i])
    
21.         String=String+str(x)+" "
    
22.         i+=1
    
23.         if i%10 == 0 :
    
24.            String=String+"\n"
    
25.     print String        
    
26. 
    
27. [root@AnsibleC14 tmp]# ./aa.py
    
28. VBM
    
29.    1.00000000000000
    
30.      3.844636    0.000000    0.000000
    
31.      1.922318    3.329553    0.000000
    
32.      1.922318    1.109851    3.139132
    
33.    Si
    
34.      2
    
35. Direct
    
36.   0.000000  0.000000  0.000000
    
37.   0.250000  0.250000  0.250000
    
38. 
    
39.    56   56   56
    
40. -4.8474 -4.4954 -3.4792 -1.91258 0.02979 2.1388 4.204 6.0466 7.5444 8.638
    
41. 9.3252 9.6442 9.656 9.4292 9.0306 8.5192 7.9338 7.3002 6.6422 5.9836
    
42. 5.345 4.7448 4.1968 3.7122 3.2998 2.968 2.724 2.5744 2.524 2.5744
    
43. 2.724 2.968 3.2998 3.7122 4.1968 4.7448 5.345 5.9836 6.6422 7.3002
    
44. 7.9338 8.5192 9.0306 9.4292 9.656 9.6442 9.3252 8.638 7.5444 6.0466
    
45. 4.204 2.1388 0.02979 -1.91258 -3.4792 -4.4954 -4.4954 -3.6928 -2.1894 -0.173258
    
46. 2.1156 4.4258 6.5376 8.2936 9.6108 10.4754 10.9238 11.0204 10.8398 10.4542
    
47. 9.9282 9.311 8.6306 7.9112 7.1772 6.4518 5.7566 5.109 4.524 4.0132
    
48. 3.5876 3.2568 3.03 2.9146 2.9146 3.03 3.2568 3.5876 4.0132 4.524
    
49. 5.109 5.7566 6.4518 7.1772 7.9112 8.6306 9.311 9.9282 10.4542 10.8398
    
50. 11.0204 10.9238 10.4754 9.6108 8.2936 6.5376 4.4258 2.1156 -0.173258 -2.1894

复制代码

龙牙地主天

回复 1# vcuteym

1. awk 'NR==FNR && NR<=12{print $0}NR==FNR && NR>12{for(i=1;i<=NF;i++) a[NR,i]=$i ; next}FNR>12{for(j=1;j<=NF;j++) printf strtonum($j+a[FNR,j])"\t" ; printf "\n"}' file1 file2

复制代码

wh7211

回复 1# vcuteym

1. awk 'FILENAME==ARGV[1]{if(NR<=12){print}else{for(i=1;i<=NF;i++){a[NR" "i]=$i}};next}FNR>=13{for(i=1;i<=NF;i++){$i=a[FNR" "i]+$i};print}' file1 file2

复制代码