用shell计算小数

今天朋友给我处了一道难题

有一组数据，数据格式如下：每行一个数据，以列的形式存在，一共是140行

zhaoyj]# tail -7 num.txt
100.2457
100.2458
100.2459
100.2457
100.2456
100.2457
100.2455

用这些数据按照以下公式进行计算：

在讲解这个公式之前我先给大家介绍一下以下几个变量，帮助大家理解这个公式

Sum1： 所有数据之和

Count： 行数之和

Avg： 平均值（sum/count）

Minus:每个数据与平均值之差

Square：差值的平方

Sum2：平方之和

Qoot： 开方根

Sum1=$line1+$line2+$line3+.......

Avg=$Sum1/$Count

Minus1=$line1-$Avg ;Minus2=$line2-$Avg; Minus3=$line3-$Avg ; . . . . . .

Square1=$Minus1*$Minus1; Square2=$Minus2*$Minus2; Square3=$Minus3*$Minus3 ; . . . . . . .

Sum2=$Square1+$Square2+$Square3+......

Qoot=Sum2的开方根

实现以上功能的脚本如下：

zhaoyj]# cat haha.sh
#!/usr/local/bin/bash
sum=0
SUM=0
count=0
while read line
do
sum1=`echo $sum1+ $line | bc`
let count++
done < num.txt

while read line
do
avg=`echo "scale=4;$sum/$count" | bc`
minus=`echo $line-$avg|bc`
square=`echo "scale=9;($minus)^2"|bc`
sum2=`echo $sum2+$square|bc`
root=`echo "scale=9;sqrt($sum2)"|bc`
done < num.txt
printf "%2.8f/n" $root

运行结果如下：

[root@mby /home/sysadmin/zhaoyj]# ./haha.sh
0.00473286

shell本身不支持小数点的计算，所以得用到bc

最后需要将输出格式化以后才能显示正确结果，要不然不会显示整数部分，只显示小数点与小数点后面的部分

<!--EndFragment-->