免费注册 查看新帖 |

ChinaUnix.net

  平台 论坛 博客 文库
最近访问板块 发新帖
查看: 1386 | 回复: 3

[文本处理] 两个文件,对应位置数据求和 [复制链接]

论坛徽章:
1
程序设计版块每日发帖之星
日期:2016-05-04 06:20:00
发表于 2018-06-01 10:17 |显示全部楼层
问题描述:
第一步,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



论坛徽章:
21
申猴
日期:2014-10-10 15:56:39射手座
日期:2014-10-10 15:57:18黑曼巴
日期:2018-05-14 11:05:122016科比退役纪念章
日期:2018-05-14 11:05:0715-16赛季CBA联赛之北控
日期:2018-05-14 11:05:0015-16赛季CBA联赛之江苏
日期:2017-02-27 18:11:0715-16赛季CBA联赛之上海
日期:2018-08-15 09:48:5415-16赛季CBA联赛之佛山
日期:2018-07-20 17:14:23
发表于 2018-06-01 11:48 |显示全部楼层
  1. [root@AnsibleC14 tmp]# cat aa.py
  2. #!/usr/bin/python

  3. import os
  4. import sys

  5. File="sed -n '1,12p' 1.txt"
  6. File1="sed -n '13,$p' 1.txt"
  7. File2="sed -n '13,$p' 1.txt"

  8. String=""
  9. List1=os.popen(File1).read().strip('\r\n').split()
  10. List2=os.popen(File2).read().strip('\r\n').split()

  11. print os.popen(File).read().strip('\r\n')

  12. if len(List1)==len(List2) :
  13.     i=0
  14.     while i<len(List1) :
  15.         x=float(List1[i])+float(List2[i])
  16.         String=String+str(x)+" "
  17.         i+=1
  18.         if i%10 == 0 :
  19.            String=String+"\n"
  20.     print String        

  21. [root@AnsibleC14 tmp]# ./aa.py
  22. VBM
  23.    1.00000000000000
  24.      3.844636    0.000000    0.000000
  25.      1.922318    3.329553    0.000000
  26.      1.922318    1.109851    3.139132
  27.    Si
  28.      2
  29. Direct
  30.   0.000000  0.000000  0.000000
  31.   0.250000  0.250000  0.250000

  32.    56   56   56
  33. -4.8474 -4.4954 -3.4792 -1.91258 0.02979 2.1388 4.204 6.0466 7.5444 8.638
  34. 9.3252 9.6442 9.656 9.4292 9.0306 8.5192 7.9338 7.3002 6.6422 5.9836
  35. 5.345 4.7448 4.1968 3.7122 3.2998 2.968 2.724 2.5744 2.524 2.5744
  36. 2.724 2.968 3.2998 3.7122 4.1968 4.7448 5.345 5.9836 6.6422 7.3002
  37. 7.9338 8.5192 9.0306 9.4292 9.656 9.6442 9.3252 8.638 7.5444 6.0466
  38. 4.204 2.1388 0.02979 -1.91258 -3.4792 -4.4954 -4.4954 -3.6928 -2.1894 -0.173258
  39. 2.1156 4.4258 6.5376 8.2936 9.6108 10.4754 10.9238 11.0204 10.8398 10.4542
  40. 9.9282 9.311 8.6306 7.9112 7.1772 6.4518 5.7566 5.109 4.524 4.0132
  41. 3.5876 3.2568 3.03 2.9146 2.9146 3.03 3.2568 3.5876 4.0132 4.524
  42. 5.109 5.7566 6.4518 7.1772 7.9112 8.6306 9.311 9.9282 10.4542 10.8398
  43. 11.0204 10.9238 10.4754 9.6108 8.2936 6.5376 4.4258 2.1156 -0.173258 -2.1894
复制代码

论坛徽章:
6
CU十四周年纪念徽章
日期:2017-08-28 17:16:1915-16赛季CBA联赛之北控
日期:2018-03-29 11:39:0615-16赛季CBA联赛之上海
日期:2018-05-21 10:46:10技术图书徽章
日期:2018-09-03 13:48:54技术图书徽章
日期:2018-09-03 13:54:18技术图书徽章
日期:2018-09-03 13:54:28
发表于 2018-06-01 14:33 |显示全部楼层
本帖最后由 龙牙地主天 于 2018-06-01 14:53 编辑

回复 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
复制代码

论坛徽章:
21
程序设计版块每日发帖之星
日期:2016-05-03 06:20:0015-16赛季CBA联赛之深圳
日期:2018-06-15 14:59:3715-16赛季CBA联赛之八一
日期:2018-07-03 16:56:4615-16赛季CBA联赛之八一
日期:2018-07-05 10:34:09黑曼巴
日期:2018-07-06 15:19:5015-16赛季CBA联赛之佛山
日期:2018-08-03 13:19:3315-16赛季CBA联赛之山西
日期:2018-08-07 19:46:2315-16赛季CBA联赛之广夏
日期:2018-08-08 19:31:5015-16赛季CBA联赛之青岛
日期:2018-11-26 15:21:5015-16赛季CBA联赛之青岛
日期:2018-06-08 13:45:2815-16赛季CBA联赛之同曦
日期:2018-06-04 19:42:2015-16赛季CBA联赛之山东
日期:2018-05-30 12:44:59
发表于 2018-06-05 16:44 |显示全部楼层
回复 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
复制代码
您需要登录后才可以回帖 登录 | 注册

本版积分规则 发表回复

数据风云,十年变迁
DTCC 第十届中国数据库技术大会已启航!

2019年5月8日~5月10日,由IT168旗下ITPUB企业社区平台主办的第十届中国数据库技术大会(DTCC2019),将在北京隆重召开。大会将邀请百余位行业专家,就热点技术话题进行分享,是广大数据领域从业人士的又一次年度盛会和交流平台。与SACC2018类似,本届大会将采用“3+2”模式:3天传统技术演讲+2天深度主题培训。大会不仅提供超100场的主题演讲,还会提供连续2天的深度课程培训,深化数据领域的项目落地实践方案。
DTCC2019,一场值得期待的数据技术盛会,殷切地希望您报名参与!

活动入口>>
  

北京盛拓优讯信息技术有限公司. 版权所有 16024965号-6 北京市公安局海淀分局网监中心备案编号:11010802020122
中国互联网协会会员  联系我们:huangweiwei@it168.com
感谢所有关心和支持过ChinaUnix的朋友们 转载本站内容请注明原作者名及出处

清除 Cookies - ChinaUnix - Archiver - WAP - TOP