兩個文件，對應(yīng)位置數(shù)據(jù)求和

vcuteym

問題描述：
第一步，file1和file2的前12行一樣，將file1和file2合并成file3，前12行保持不變。
第二步，file1和file2的第13行至最后一行，對應(yīng)數(shù)據(jù)求和，合成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

復(fù)制代碼

龍牙地主天

回復(fù) 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

復(fù)制代碼

wh7211

回復(fù) 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

復(fù)制代碼