☰

AWK는 인라인 나누기를 시뮬레이션하고 절대값을 취하는 데 사용됩니다.

tag-icon • tag-icon shell-script awk syntax
AWK는 인라인 나누기를 시뮬레이션하고 절대값을 취하는 데 사용됩니다.

다음 형식의 데이터 파일이 있습니다. (윗줄은 참고용입니다.)

    E          T               N         D                 S             s                 R          k
   0.45847  300.0000      9.99979156  0.44254427E+02  0.49221658E-04  0.98763353E+19  0.73239256E-09  0.24892539E+15
   0.45947  300.0000      9.97603847  0.51763106E+02 -0.24189016E-04  0.97909193E+19 -0.24467623E-08  0.23855274E+15

이제 다음을 포함하는 동일한 열 형식의 파일을 생성해야 합니다.

N         (k/(s/300.0000))         Abs(S)/E         exp(Abs(S)/E)

즉, 먼저 두 열을 나누고 Column(5)의 절대값을 취한 다음 다시 거듭제곱을 수행해야 합니다.

나는 코드를 사용하고 있습니다 :

while read -a line; do echo -e " ${line[2]}\t" `awk "BEGIN {print ${line[7]}/${line[5]}/300}"` ; done < filename

하지만 이렇게 하면 두 번째 열만 볼 수 있습니다. CI는 다음 두 열과 같이 절대값을 취해야 합니까?

나는 이미 도착했습니다:

 while read -a line; do echo -e ${line[4]} | awk ' { if($line[4]>=0 ) {print $line[4] } else {print $line[4]*(-1)} }'; done <

하지만 마지막 열 단위 나누기 코드와 결합할 수는 없습니다. 일부 구문 BEGIN이나 END사용법이 누락되었나요 awk? 결과를 10진수 형식으로 표시하는 것도 혼란스럽습니다. ㅏ

PS: awk이 작업을 수행하는 데 적합한 도구인가요? 나는 그것을 시도하고 싶습니까 python? 이와 같은 파일이 20개가 넘으므로 스프레드시트를 사용할 수 없다는 것을 알고 있습니다.

답변1

예, awk이것이 올바른 도구입니다. 다음과 같이 할 수 있습니다.

$ awk 'function abs(x) {return x<0 ? -x : x}
     NR>1{printf ("%e %e %e %e\n", $3, $NF/($6/300.0000), abs($5)/$1, exp(abs($5)/$1) );
}' infile

9.999792e+00 7.561268e-03 1.073607e-04 1.000107e+00
9.976038e+00 7.309408e-03 5.264547e-05 1.000053e+00

위와 같이 abs() 함수를 시뮬레이션하거나 함수를 정의하지 않고 삼항 연산자를 사용하여 다음과 같이 생성할 수 있습니다.

$ awk 'NR>1{ $5*=($5<0?-1:1);
    printf("%e %e %e %e\n", $3, $NF/($6/300.0000), $5/$1, exp($5/$1) )}' infile
9.999792e+00 7.561268e-03 1.073607e-04 1.000107e+00
9.976038e+00 7.309408e-03 5.264547e-05 1.000053e+00

답변2

Python 스크립트를 사용하여 이 작업을 수행하려면 다음을 수행하십시오.

#!/usr/bin/env python
##
# use: ingest.py infile1 infile2 ... infilen

import math,re,argparse

parser = argparse.ArgumentParser()
parser.add_argument('infiles', type=str, help='infile names', nargs='+')
args = parser.parse_args()

for infile in args.infiles:
    with open(infile) as in_file, open(infile+".out", "w+") as out_file:
        column_names = in_file.readline().split()
        out_file.write("N (k/(s/300.0000)) Abs(S)/E exp(Abs(S)/E)\n")
        for line in in_file.readlines():
            column_values = [float(x) for x in line.split()]
            value_dict = dict(zip(column_names,column_values))
            div = value_dict["k"]/(value_dict["s"]/300.0000)
            ABS = abs(value_dict["S"])/value_dict["E"]
            EXP = math.exp(abs(value_dict["S"])/value_dict["E"])
            out_file.write("{N} {DIV} {ABS} {EXP}\n".format(N=value_dict["N"],DIV=div,ABS=ABS, EXP=EXP))

관련 정보