0

我有如下输入

*KEYWORD
$TIME_VALUE = 9.9999993e-004
$STATE_NO = 2
$Output for State 2 at time = 0.001
*END
$NODAL_RESULTS
$RESULT OF Resultant Displacement
721810 1.7188E-2
721812 6.1973E-2
721825 1.1481E+0
721827 1.0962E+0
721852 5.1831E-1
721854 1.3085E-2
721867 1.1077E+0
 .        .
 .        .
 .        .

我需要找到第 2 列中的最大值及其平均值。然后我还需要输出代表最大值的第一列中的数字。

我使用以下代码计算最大值和平均值，但是除以零。

awk: cmd. line:5: fatal: division by zero attempted

代码如下

# 1.k is the input file name.
sed  -n '/^[0-9]\{1\}/p' 1.k  > 2.k   # delete all lines not starting with number
mv 2.k 1.k
sed  -i -e '/^$/d' 1.k                # delete all lines that are empty
#sed -i -e  's/^[ \t]*//;s/[ \t]*$//' 1.k  

awk 'BEGIN{min=999}
{a[NR]=$0;if($2<min){min=$2;m[1]=NR;}if($2>max){max=$2;m[2]=NR;}m[2]+=$2;}
END{print "Min:"a[m[1]];
print "Max:"a[m[2]];
print "Number Of Nodes:" NR;
print "Avg:"m[3]/NR}' 1.k

有人可以帮我解决这个问题吗？

问候，

4

5 回答 5

0

计算.awk：

{ 
  sum += $2  
  if (NR == 1)  {   
    min = max = $2
    minv= maxv= $1
  }
  if (min > $2) { min = $2; minv = $1 }
  if (max < $2) { max = $2; maxv = $1 }
}
END { 
  print "Min: " minv ", " min 
  print "Avg: " sum / NR  
  print "Max: " maxv ", " max 
  print "# Nodes: " NR
}
于 2012-08-24T12:41:27.600 回答
0

经过尝试和尝试，我找到了一个可行的方法，但我认为这不是问题的最佳解决方案。

sed -i -e  's/^[ \t]*//;s/[ \t]*$//' 1.k
sed  -n '/^[0-9]\{1\}/p' 1.k  > 2.k
mv 2.k 1.k
sed  -i -e '/^$/d' 1.k  



awk 'BEGIN{min=999}

{a[NR]=$0;if($2<min){min=$2;m[1]=NR;}if($2>max){max=$2;m[2]=NR;}m[3]+=$2;}

END{     print "Max:"a[m[2]];     
     print "Min:"a[m[1]];     
    print "Number Of Calls:" NR;  
    print "Avg:"m[3]/NR}' 1.k   > result

感谢您宝贵的建议。

于 2012-08-24T14:01:16.647 回答
0

尝试简单的-r解决方案：

perl -walne 'print $F[1] if /^\d/' | r summary -

Simple-r 是一个 R 包装器，用于在命令行中进行快速统计分析。可以在以下位置找到：

https://code.google.com/p/simple-r/

于 2013-10-01T16:59:28.910 回答
0

如果您首先过滤掉非数字信息，那么此awk脚本应该执行以下操作：

awk 'BEGIN{max=-999}\
{\
  col1[NR]=$1;\
  col2[NR]=$2;\
  if($2>max){max=$2;imax=NR};\
  sum+=$2\
}\
END{print col1[imax]" "col2[imax]" average: "sum/NR}' yourinputfile
于 2012-08-24T12:57:08.987 回答
0

Perl 解决方案：

<1.k perl -ne 'next unless ($key, $val) = /^([0-9]+)\s+([-+E.0-9]+)/;# Only process the important lines.
               if ($val > $max) {                                    # New maximum.
                   $max = $val;
                   @maxk = $key;
               } elsif ($max == $val) {                              # The maximum appears more than once.
                   push @maxk, $key;
               }
               $sum += $val;
               $count++;
           } {
               print "MAX: $max at @maxk, AVG: ", $sum / $count ,"\n"; '
于 2012-08-24T12:38:51.167 回答