列依次求比值

angjzy

菜鸟求助，

源文件，7行8列

X1	X2	X3	X4	X5	X6
1	717.2	10.6	61.2	9.383018	1.30037
1.3	299.8	36.7	59.9	30.50118	0.229287
0.5	418	10	38.3	44.2	0.061
0.8	580	11	42.4	46	0.481
0.5	395	31.84	41.5	41.9	0.131
0.8	441	10.6	36.5	27.7	0.164
1.8	357.8	44	48.2	53.01714	0.354479

处理生成新文件得到 7行15列 数据
X1/X2 X1/X3 X1/X4 X1/X5 X1/X6    X2/X3 X2/X4 X2/X5 X2/X6  。。。直到X5/X6

谢谢高手

jason680

$ awk -vOFS="\t" 'NR>1{_=1;_OO=NF;for(_0=_;_0<_OO;++_0){_O0=_;for(_O=_0+_;_O<=_OO;++_O)_O0*=$_0/$_O;$(_0+_OO)=_O0}print}' FILE
1        717.2        10.6        61.2        9.383018        1.30037        1.76154e-07        3.34268e+07        1.59499        306.968        7.21565
1.3        299.8        36.7        59.9        30.50118        0.229287        8.05557e-07        525457        117.998        513.047        133.026
0.5        418        10        38.3        44.2        0.061        7.23974e-08        2.95634e+07        9.68387        544.058        724.59
0.8        580        11        42.4        46        0.481        5.4747e-08        1.09661e+07        1.41876        81.251        95.6341
0.5        395        31.84        41.5        41.9        0.131        1.0908e-08        3.35646e+06        141.705        313.77        319.847
0.8        441        10.6        36.5        27.7        0.164        4.22755e-07        2.15195e+07        7.18292        293.266        168.902
1.8        357.8        44        48.2        53.01714        0.354479        1.325e-06        411202        94.0382        123.62        149.564

jason680

Resend ....

angjzy

回复 4# jason680

谢谢高手，但是数值好像不太对，第七列第一个数应该是1/717.2 （X1/X2）, 第八列第一个数是1/10.6(X1/X3)

jason680

回复 5# angjzy

$ awk -vOFS="\t" 'NR>1{_=1;_OO=NF;for(_0=_;_0<_OO;++_0){for(_O=_0+_;_O<=_OO;++_O)$(++__+_OO)=$_0/$_O;}print}' FILE
1    717.2    10.6    61.2    9.383018    1.30037    0.00139431    0.0943396    0.0163399    0.106576    0.769012    67.6604    11.719    76.436    551.535    0.173203    1.1297    8.15153    6.52242    47.0635    7.21565
1.3    299.8    36.7    59.9    30.50118    0.229287            0.00433622    0.0354223    0.0217028    0.0426213    5.66975    8.168945.00501    9.82913    1307.53    0.612688    1.20323    160.061    1.96386    261.245    133.026
...

chengchow

1. [root@zbxagent01 temp]# cat a.py
   
2. #!/usr/bin/env python
   
3. 
   
4. import os
   
5. import re
   
6. import sys
   
7. 
   
8. File="/data/temp/a"
   
9. 
   
10. Column=int(os.popen("cat " + File + " | wc  -l").read().strip('\r\n'))
    
11. Column+=1
    
12. Str=os.popen("awk 'NR!=1' " + File).read().strip('\r\n')
    
13. 
    
14. for i in range(2,Column) :
    
15.         Str=os.popen("awk -v num=" + str(i) + " 'NR!=1&&NR==num' " + File).read().strip('\r\n')
    
16.         Lst=Str.split()
    
17.         LstLength=len(Lst)
    
18.         for m in Lst :
    
19.                 BeginElementIndex=Lst.index(m)
    
20.                 for n in range(BeginElementIndex+1,LstLength) :
    
21.                         print m+"/"+Lst[n],
    
22.                 print "\t",
    
23.         print "\n",
    
24. [root@zbxagent01 temp]# ./a.py
    
25. 1/717.2 1/10.6 1/61.2 1/9.383018 1/1.30037         717.2/10.6 717.2/61.2 717.2/9.383018 717.2/1.30037         10.6/61.2 10.6/9.383018 10.6/1.30037         61.2/9.383018 61.2/1.30037         9.383018/1.30037                 
    
26. 1.3/299.8 1.3/36.7 1.3/59.9 1.3/30.50118 1.3/0.229287         299.8/36.7 299.8/59.9 299.8/30.50118 299.8/0.229287         36.7/59.9 36.7/30.50118 36.7/0.229287         59.9/30.50118 59.9/0.229287         30.50118/0.229287                 
    
27. 0.5/418 0.5/10 0.5/38.3 0.5/44.2 0.5/0.061         418/10 418/38.3 418/44.2 418/0.061         10/38.3 10/44.2 10/0.061         38.3/44.2 38.3/0.061         44.2/0.061                 
    
28. 0.8/580 0.8/11 0.8/42.4 0.8/46 0.8/0.481         580/11 580/42.4 580/46 580/0.481         11/42.4 11/46 11/0.481         42.4/46 42.4/0.481         46/0.481                 
    
29. 0.5/395 0.5/31.84 0.5/41.5 0.5/41.9 0.5/0.131         395/31.84 395/41.5 395/41.9 395/0.131         31.84/41.5 31.84/41.9 31.84/0.131         41.5/41.9 41.5/0.131         41.9/0.131                 
    
30. 0.8/441 0.8/10.6 0.8/36.5 0.8/27.7 0.8/0.164         441/10.6 441/36.5 441/27.7 441/0.164         10.6/36.5 10.6/27.7 10.6/0.164         36.5/27.7 36.5/0.164         27.7/0.164                 
    
31. 1.8/357.8 1.8/44 1.8/48.2 1.8/53.01714 1.8/0.354479         357.8/44 357.8/48.2 357.8/53.01714 357.8/0.354479         44/48.2 44/53.01714 44/0.354479         48.2/53.01714 48.2/0.354479         53.01714/0.354479

复制代码

1cpuer

回复 6# jason680

# [ /home/soio/1bs/awks ] {2018-06-11 09:26:57}
: 1528680417:0;➜  sed ';s/_0/i/g; s/\<_O\>/j/g' jason.6
$ awk -vOFS="\t" 'NR>1{_=1;_OO=NF;for(i=_;i<_OO;++i){for(j=i+_;j<=_OO;++j)$(++__+_OO)=$i/$j;}print}' FILE