如何用Shell处理两个文档中多行数据的加法？

Osirix

文档1格式是这个样子的：

1. 
   
2. [...]
   
3. #下面是需要处理的数据，前三列是三维坐标，最后一列是“标号”
   
4.   0.9986  0.9231  0.8750    1
   
5.   0.4986  0.9231  0.8750    2
   
6.   0.9986  0.4231  0.8750    3
   
7.   0.4986  0.4231  0.8750    4
   
8.   0.9986  0.9231  0.3750    5
   
9.   0.4986  0.9231  0.3750    6
   
10.   0.9986  0.4231  0.3750    7
    
11.   0.4986  0.4231  0.3750    8
    
12.   0.5013  0.5768  0.6250    9
    
13.   0.0013  0.5768  0.6250   10
    
14. #end of file1
    

复制代码

文档2的格式是这样的:

1. 
   
2. #number of files will be created
   
3.     7
   
4. #下面每一行中前三个数据（即三维坐标）需要分别加到文档1中同“标号”所在行的三维坐标中后将文档1生成一个新的文档
   
5. #上一行中的7，标示下面总共有7行，即需要产生7个新的文档
   
6. 0.1000   0.0000  0.0000  2
   
7. 0.0000   0.0000 -0.1100  2
   
8. 0.2000   0.0000  0.0000  5
   
9. 0.0000  -0.1200  0.0000  5
   
10. 0.5000   0.0000  0.0000  8
    
11. 0.0000   0.5000  0.5000  8
    
12. 0.9000   0.0000  0.0000  10
    
13. #end of file2
    

复制代码

这个有啥好法子可以实现呢？

Osirix

譬如根据file2中的下面这一行:

1. 
   
2. 0.1000   0.0000  0.0000  2
   

复制代码

表示应该对file1中的标号为2的这一行

1. 
   
2.   0.4986  0.9231  0.8750    2
   

复制代码

进行加法,
从而产生下面一个新的文档

1. 
   
2. [...]
   
3. #下面是需要处理的数据，前三列是三维坐标，最后一列是“标号”
   
4.   0.9986  0.9231  0.8750    1
   
5.   0.5986  0.9231  0.8750    2
   
6.   0.9986  0.4231  0.8750    3
   
7.   0.4986  0.4231  0.8750    4
   
8.   0.9986  0.9231  0.3750    5
   
9.   0.4986  0.9231  0.3750    6
   
10.   0.9986  0.4231  0.3750    7
    
11.   0.4986  0.4231  0.3750    8
    
12.   0.5013  0.5768  0.6250    9
    
13.   0.0013  0.5768  0.6250   10
    
14. #end of file1
    

复制代码

lastart

1. 
   
2. #!/usr/bin/mawk -f
   
3. BEGIN{
   
4.     total=1;
   
5.     while(getline temp<"file2"){
   
6.         if (match(temp,"^ *#")) continue;
   
7.         tNF=split(temp,t);
   
8.                 if(tNF!=4)continue;
   
9.         for(i=1;i<total;i++) if (data[i,4]==t[4]) {
   
10.             for(j=1;j<4;j++) data[i,j]+=t[j];
    
11.             break;}
    
12.         if(i==total){for(j=1;j<=4;j++) data[i,j]=t[j];
    
13.             total++;}
    
14.     }
    
15. }
    
16. {for(i=1;i<total;i++) if (data[i,4] && $4 == data[i,4]){$1+=data[i,1];$2+=data[i,2];$3+=data[i,3];delete data[i,4]}
    
17. print
    
18. }
    

复制代码

运行awk -f file1. 输出的格式会乱, 但应该影响不在