高级Bash脚本编程指南(31)：数学计算命令

将一个正数分解为多个素数.

[html]  view plain copy

1. root@ubuntu:~/resource/shell-study/0622-2013# factor 2374  
2. 2374: 2 1187  
3. root@ubuntu:~/resource/shell-study/0622-2013# factor 27417  
4. 27417: 3 13 19 37  
5. root@ubuntu:~/resource/shell-study/0622-2013#   

信不过他的话你可以自己去试试，O(∩_∩)O~

bc

Bash不能处理浮点运算, 并且缺乏特定的一些操作, 这些操作都是一些重要的计算功能. 幸运的是, bc可以解决这个问题.

bc不仅仅是个多功能灵活的精确计算工具, 且它还提供许多编程语言才具备的一些方便功能，bc比较类似于C语言的语法，因为它是一个完整的UNIX工具, 所以它可以用在pipe中, bc在脚本中也是很常用的.

这里有一个简单的使用bc命令的模版, 可以用来计算脚本中的变量. 这个模版经常用于命令替换中.

variable=$(echo "OPTIONS; OPERATIONS" | bc)

一个实例，不过不推荐你去看其中的算法，你只要只要怎么用bc去计算就够了

[html]  view plain copy

1. #!/bin/bash  
2. # 计算按月偿还贷款的数量.  
3.   
4. echo  
5. echo "Given the principal, interest rate, and term of a mortgage,"  
6. echo "calculate the monthly payment."  
7.   
8. bottom=1.0  
9.   
10. echo  
11. echo -n "Enter principal (no commas) "  
12. read principal  
13. echo -n "Enter interest rate (percent)"#如果是12%,那就键入"12",而不是".12".  
14. read interest_r  
15. echo -n "Enter term (months) "  
16. read term  
17.   
18. interest_r=$(echo "scale=9; $interest_r/100.0" | bc) # 转换成小数.  
19. # "scale"指定了有效数字的个数.  
20.   
21. interest_rate=$(echo "scale=9; $interest_r/12 + 1.0" | bc)  
22.   
23. top=$(echo "scale=9; $principal*$interest_rate^$term" | bc)  
24.   
25. echo; echo "Please be patient. This may take a while."  
26.   
27. let "months = $term - 1"  
28. # ====================================================================  
29. for ((x=$months; x > 0; x--))  
30. do  
31.    bot=$(echo "scale=9; $interest_rate^$x" | bc)  
32.    bottom=$(echo "scale=9; $bottom+$bot" | bc)  
33. #  bottom = $(($bottom + $bot"))  
34.  done  
35. # ====================================================================  
36.   
37. # --------------------------------------------------------------------  
38. #  Rick Boivie给出了一个对上边循环的修改方案,  
39. #+ 这个修改更加有效率, 将会节省大概2/3的时间.  
40.   
41. # for ((x=1; x <= $months; x++))  
42. # do  
43. #   bottom=$(echo "scale=9; $bottom * $interest_rate + 1" | bc)  
44. # done  
45.   
46.   
47. #  然后他又想出了一个更加有效率的版本,  
48. #+ 将会节省95%的时间!  
49.   
50. # bottom=`{  
51. #     echo "scale=9; bottom=$bottom; interest_rate=$interest_rate"  
52. #     for ((x=1; x <= $months; x++))  
53. #     do  
54. #          echo 'bottom = bottom * interest_rate + 1'  
55. #     done  
56. #     echo 'bottom'  
57. #     } | bc`       # 在命令替换中嵌入一个'for循环'.  
58. # -------------------------------------------------------------------------  
59. #  另一方面, Frank Wang建议:  
60. #  bottom=$(echo "scale=9; nterest_rate^$term-1)/($interest_rate-1)" | bc)  
61.   
62. #  因为 . . .  
63. #  在循环后边的算法  
64. #+ 事实上是一个等比数列的求和公式.  
65. #  求和公式是 e0(1-q^n)/(1-q),  
66. #+ e0是第一个元素, q=e(n+1)/e(n),  
67. #+ n是元素数量.  
68. # --------------------------------------------------------------------------  
69.   
70.   
71.  # let "payment = $top/$bottom"  
72.  payment=$(echo "scale=2; $top/$bottom" | bc)  
73.  # 使用2位有效数字来表示美元和美分.  
74.   
75.  echo  
76.  echo "monthly payment = \$$payment"  # 在总和的前边显示美元符号.  
77.  echo  
78.   
79.  exit 0  

结果：

[html]  view plain copy

1. root@ubuntu:~/resource/shell-study/0622-2013# ./test10.sh   
2.   
3. Given the principal, interest rate, and term of a mortgage,  
4. calculate the monthly payment.  
5.   
6. Enter principal (no commas) 50000  
7. Enter interest rate (percent)#如果是12%,那就键入12,而不是.12.5  
8. Enter term (months) 12  
9.   
10. Please be patient. This may take a while.  
11.   
12. monthly payment = $4280.37  
13.   
14. root@ubuntu:~/resource/shell-study/0622-2013#   

接着看一个数制转换的实例，挺不错的实例啊，可以研究一下

[html]  view plain copy

1. #!/bin/bash  
2. ##########################################################################  
3. # 脚本       :    base.sh - 用不同的数制来打印数字 (Bourne Shell)  
4. # 作者       :    Heiner Steven (heiner.steven@odn.de)  
5. # 日期       :    07-03-95  
6. # 类型       :    桌面  
7. # $Id: base.sh,v 1.2 2000/02/06 19:55:35 heiner Exp $  
8. # ==> 上边这行是RCS ID信息.  
9. ##########################################################################  
10. # 描述  
11. #  
12. # 修改纪录  
13. # 21-03-95 stv  fixed error occuring with 0xb as input (0.2)  
14. ##########################################################################  
15.   
16. # ==> 在本书中使用这个脚本通过了原作者的授权.  
17. # ==> 注释是本书作者添加的.  
18.   
19. NOARGS=65  
20. PN=`basename "$0"`                 # 程序名  
21. VER=`echo '$Revision: 1.2 $' | cut -d' ' -f2`  # ==> VER=1.2  
22.   
23. Usage () {  
24.     echo "$PN - print number to different bases, $VER (stv '95)  
25. usage: $PN [number ...]  
26.   
27. If no number is given, the numbers are read from standard input.  
28. A number may be  
29.     binary (base 2)     starting with 0b (i.e. 0b1100)  
30.     octal (base 8)      starting with 0  (i.e. 014)  
31.     hexadecimal (base 16)   starting with 0x (i.e. 0xc)  
32.     decimal         otherwise (i.e. 12)" >&2  
33.     exit $NOARGS  
34. }   # ==> 打印出用法信息的函数.  
35.   
36. Msg () {  
37.     for i   # ==> 省略[list].  
38.     do echo "$PN: $i" >&2  
39.     done  
40. }  
41.   
42. Fatal () { Msg "$@"; exit 66; }  
43.   
44. PrintBases () {  
45.     # 决定数字的数制  
46.     for i      # ==> 省略[list]...  
47.     do         # ==> 所以是对命令行参数进行操作.  
48.     case "$i" in  
49.         0b*)        ibase=2;;   # 2进制  
50.         0x*|[a-f]*|[A-F]*)  ibase=16;;  # 16进制  
51.         0*)         ibase=8;;   # 8进制  
52.         [1-9]*)     ibase=10;;  # 10进制  
53.         *)  
54.         Msg "illegal number $i - ignored"  
55.         continue;;  
56.     esac  
57.   
58.     # 去掉前缀, 将16进制数字转换为大写(bc命令需要这么做)  
59.     number=`echo "$i" | sed -e 's:^0[bBxX]::' | tr '[a-f]' '[A-F]'`  
60.     # ==> 使用":" 作为sed分隔符, 而不使用"/".  
61.   
62.     # 将数字转换为10进制  
63.     dec=`echo "ibase=$ibase; $number" | bc`  # ==> 'bc'是个计算工具.  
64.     case "$dec" in  
65.         [0-9]*) ;;           # 数字没问题  
66.         *)      continue;;       # 错误: 忽略  
67.     esac  
68.   
69.     # 在一行上打印所有转换后的数字.  
70.     # ==> 'here document'提供命令列表给'bc'.  
71.     echo `bc <<!  
72.         obase=16; "hex="; $dec  
73.         obase=10; "dec="; $dec  
74.         obase=8;  "oct="; $dec  
75.         obase=2;  "bin="; $dec  
76. !` | sed -e 's: :   :g'  
77.   
78.     done  
79. }  
80.   
81. while [ $# -gt 0 ]  
82. # ==>  这里必须使用一个"while循环",  
83. # ==>+ 因为所有的case都可能退出循环或者  
84. # ==>+ 结束脚本.  
85. # ==> (感谢, Paulo Marcel Coelho Aragao.)  
86. do  
87.     case "$1" in  
88.     --)     shift; break;;  
89.     -h)     Usage;;                 # ==> 帮助信息.  
90.     -*)     Usage;;  
91.     *)      break;;         # 第一个数字  
92.     esac   # ==> 对于非法输入进行更严格检查是非常有用的.  
93.     shift  
94. done  
95.   
96. if [ $# -gt 0 ]  
97. then  
98.     PrintBases "$@"  
99. else                    # 从stdin中读取  
100.     while read line  
101.     do  
102.     PrintBases $line  
103.     done  
104. fi  
105.   
106.   
107. exit 0  

结果：

[html]  view plain copy

1. root@ubuntu:~/resource/shell-study/0622-2013# ./test11.sh -h  
2. test11.sh - print number to different bases, 1.2 (stv '95)  
3. usage: test11.sh [number ...]  
4.   
5. If no number is given, the numbers are read from standard input.  
6. A number may be  
7.     binary (base 2)     starting with 0b (i.e. 0b1100)  
8.     octal (base 8)      starting with 0  (i.e. 014)  
9.     hexadecimal (base 16)   starting with 0x (i.e. 0xc)  
10.     decimal         otherwise (i.e. 12)  
11. root@ubuntu:~/resource/shell-study/0622-2013# ./test11.sh   
12. 6  
13. hex=6   dec=6   oct=6   bin=110  
14. 15  
15. hex=F   dec=15  oct=17  bin=1111  
16. D  
17. hex=D   dec=13  oct=15  bin=1101  
18. ^C  
19. root@ubuntu:~/resource/shell-study/0622-2013#   

当一个脚本需要将一个选项列表和多个命令传递到bc中时, 这种方法就显得非常合适了.

[html]  view plain copy

1. variable=`bc << LIMIT_STRING  
2. options  
3. statements  
4. operations  
5. LIMIT_STRING  
6. `  
7. ...or...  
8.   
9. variable=$(bc << LIMIT_STRING  
10. options  
11. statements  
12. operations  
13. LIMIT_STRING  
14. )  

其实上面的第一种方式在上面的实例中已经使用过了，很好，很强大

接着看一个实例：使用"here document"来调用bc

[html]  view plain copy

1. #!/bin/bash  
2. # 使用命令替换来调用'bc'  
3. # 并与'here document'相结合.  
4.   
5. var1=`bc << EOF  
6. 18.33 * 19.78  
7. EOF  
8. `  
9. echo $var1       # 362.56  
10.   
11. #  使用$( ... )这种标记法也可以.  
12. v1=23.53  
13. v2=17.881  
14. v3=83.501  
15. v4=171.63  
16.   
17. var2=$(bc << EOF  
18. scale = 4  
19. a = ( $v1 + $v2 )  
20. b = ( $v3 * $v4 )  
21. a * b + 15.35  
22. EOF  
23. )  
24. echo $var2       # 593487.8452  
25.   
26.   
27. var3=$(bc -l << EOF  
28. scale = 9  
29. s ( 1.7 )  
30. EOF  
31. )  
32. # 返回弧度为1.7的正弦.  
33. # "-l"选项将会调用'bc'算数库.  
34. echo $var3       # .991664810  
35.   
36.   
37. # 现在, 在函数中试一下...  
38. hyp=             # 声明全局变量.  
39. hypotenuse ()    # 计算直角三角形的斜边.  
40. {  
41. hyp=$(bc -l << EOF  
42. scale = 9  
43. sqrt ( $1 * $1 + $2 * $2 )  
44. EOF  
45. )  
46. # 不幸的是, 不能从bash函数中返回浮点值.  
47. }  
48.   
49. hypotenuse 3.68 7.31  
50. echo "hypotenuse = $hyp"    # 8.184039344  
51.   
52. exit 0  

结果：

[html]  view plain copy

1. root@ubuntu:~/resource/shell-study/0622-2013# chmod +x test12.sh   
2. root@ubuntu:~/resource/shell-study/0622-2013# ./test12.sh   
3. 362.56  
4. 593487.8452  
5. .991664810  
6. hypotenuse = 8.184039344  
7. root@ubuntu:~/resource/shell-study/0622-2013#   

一个模拟计算圆周率的实例，不过先说明，这个计算需要统计，执行的次数越多越好，我只让他跑了10次，差得十万八千里可以说，我只是测试，准不准我可不管，O(∩_∩)O~

[html]  view plain copy

1. #!/bin/bash  
2.   
3. # 这事实上是一个"Monte Carlo"蒙特卡洛模拟的非常简单的实例:  
4. #+ 蒙特卡洛模拟是一种由现实事件抽象出来的数学模型,  
5. #+ 由于要使用随机抽样统计来估算数学函数,所以使用伪随机数来模拟真正的随机数.  
6.   
7. #  想象有一个完美的正方形土地, 边长为10000个单位.  
8. #  在这块土地的中间有一个完美的圆形湖,  
9. #+ 这个湖的直径是10000个单位.  
10. #  这块土地的绝大多数面积都是水, 当然只有4个角上有一些土地.  
11. #  (可以把这个湖想象成为这个正方形的内接圆.)  
12. #  
13. #  我们将使用老式的大炮和铁炮弹  
14. #+ 向这块正方形的土地上开炮.  
15. #  所有的炮弹都会击中这块正方形土地的某个地方.  
16. #+ 或者是打到湖上, 或者是打到4个角的土地上.  
17. #  因为这个湖占据了这个区域大部分地方,  
18. #+ 所以大部分的炮弹都会"扑通"一声落到水里.  
19. #  而只有很少的炮弹会"砰"的一声落到4个  
20. #+ 角的土地上.  
21. #  
22. #  如果我们发出的炮弹足够随机的落到这块正方形区域中的话,  
23. #+ 那么落到水里的炮弹与打出炮弹总数的比率,  
24. #+ 大概非常接近于PI/4.  
25. #  
26. #  原因是所有的炮弹事实上都  
27. #+ 打在了这个土地的右上角,  
28. #+ 也就是, 笛卡尔坐标系的第一象限.  
29. #  (之前的解释只是一个简化.)  
30. #  
31. #  理论上来说, 如果打出的炮弹越多, 就越接近这个数字.  
32. #  然而, 对于shell 脚本来说一定会做些让步的,  
33. #+ 因为它肯定不能和那些内建就支持浮点运算的编译语言相比.  
34. #  当然就会降低精度.  
35.   
36.   
37. DIMENSION=10000  # 这块土地的边长.  
38.                  # 这也是所产生随机整数的上限.  
39.   
40. MAXSHOTS=10    # 开炮次数.  
41.                  # 10000或更多次的话, 效果应该更好, 但有点太浪费时间了.  
42. PMULTIPLIER=4.0  # 接近于PI的比例因子.  
43.   
44. get_random ()  
45. {  
46.     SEED=$(head -1 /dev/urandom | od -N 1 | awk '{ print $2 }')  
47.     RANDOM=$SEED                                  #  来自于"seeding-random.sh"  
48.                                                   #+ 的例子脚本.  
49.     let "rnum = $RANDOM % $DIMENSION"             #  范围小于10000.  
50.     echo $rnum  
51. }  
52.   
53. distance=        # 声明全局变量.  
54. hypotenuse ()    # 从"alt-bc.sh"例子来的,  
55. {                # 计算直角三角形的斜边的函数.  
56.     distance=$(bc -l << EOF  
57.     scale = 0  
58.     sqrt ( $1 * $1 + $2 * $2 )  
59. EOF  
60. )  
61. #  设置 "scale" 为 0 , 好让结果四舍五入为整数值,  
62. #+ 这也是这个脚本中必须折中的一个地方.  
63. #  不幸的是, 这将降低模拟的精度.  
64. }  
65.   
66. # 初始化变量.  
67. shots=0  
68. splashes=0  
69. thuds=0  
70. Pi=0  
71.   
72. while [ "$shots" -lt  "$MAXSHOTS" ]           # 主循环.  
73. do  
74.   
75.   xCoord=$(get_random)                        # 取得随机的 X 与 Y 坐标.  
76.   yCoord=$(get_random)  
77.   hypotenuse $xCoord $yCoord                  #  直角三角形斜边 =  
78.                                               #+ distance.  
79.   ((shots++))  
80.   
81.   printf "#%4d   " $shots  
82.   printf "Xc = %4d  " $xCoord  
83.   printf "Yc = %4d  " $yCoord  
84.   printf "Distance = %5d  " $distance         #  到湖中心的  
85.                                               #+ 距离 --  
86.                                               #  起始坐标点 --  
87.                                               #+  (0,0).  
88.   
89.   if [ "$distance" -le "$DIMENSION" ]  
90.   then  
91.     echo -n "SPLASH!  "  
92.     ((splashes++))  
93.   else  
94.     echo -n "THUD!    "  
95.     ((thuds++))  
96.   fi  
97.   
98.   Pi=$(echo "scale=9; $PMULTIPLIER*$splashes/$shots" | bc)  
99.   # 将比例乘以4.0.  
100.   echo -n "PI ~ $Pi"  
101.   echo  
102.   
103. done  
104.   
105. echo  
106. echo "After $shots shots, PI looks like approximately $Pi."  
107. # 如果不太准的话, 那么就提高一下运行的次数. . .  
108. # 可能是由于运行错误和随机数随机程度不高造成的.  
109.   
110. exit 0  

结果：

[html]  view plain copy

1. root@ubuntu:~/resource/shell-study/0622-2013# ./test13.sh   
2. #   1   Xc = 9684  Yc = 9797  Distance = 13775  THUD!    PI ~ 0  
3. #   2   Xc = 8460  Yc = 6413  Distance = 10615  THUD!    PI ~ 0  
4. #   3   Xc = 1111  Yc = 1189  Distance =  1627  SPLASH!  PI ~ 1.333333333  
5. #   4   Xc = 4230  Yc = 4986  Distance =  6538  SPLASH!  PI ~ 2.000000000  
6. #   5   Xc = 9403  Yc =   58  Distance =  9403  SPLASH!  PI ~ 2.400000000  
7. #   6   Xc = 9345  Yc = 5325  Distance = 10755  THUD!    PI ~ 2.000000000  
8. #   7   Xc = 6826  Yc = 4682  Distance =  8277  SPLASH!  PI ~ 2.285714285  
9. #   8   Xc = 3762  Yc = 4721  Distance =  6036  SPLASH!  PI ~ 2.500000000  
10. #   9   Xc = 7446  Yc = 5189  Distance =  9075  SPLASH!  PI ~ 2.666666666  
11. #  10   Xc = 1337  Yc = 9442  Distance =  9536  SPLASH!  PI ~ 2.800000000  
12.   
13. After 10 shots, PI looks like approximately 2.800000000.  
14. root@ubuntu:~/resource/shell-study/0622-2013#   

dc

dc(桌面计算器desk calculator)工具是面向栈的, 并且使用RPN(逆波兰表达式"Reverse Polish Notation"又叫"后缀表达式"). 与bc命令很相似, 但是这个工具具备好多只有编程语言才具备的能力.

正常表达式   逆波兰表达式
a+b                a,b,+
a+(b-c)          a,b,c,-,+
a+(b-c)*d      a,d,b,c,-,*,+

绝大多数人都避免使用这个工具, 因为它需要非直观的RPN输入, 但是, 它却有特定的用途.

一个实例：将10进制数字转换为16进制数字

[html]  view plain copy

1. #!/bin/bash  
2. # hexconvert.sh: 将10进制数字转换为16进制数字.  
3.   
4. E_NOARGS=65 # 缺少命令行参数错误.  
5. BASE=16     # 16进制.  
6.   
7. if [ -z "$1" ]  
8. then  
9.   echo "Usage: $0 number"  
10.   exit $E_NOARGS  
11.   # 需要一个命令行参数.  
12. fi  
13. # 练习: 添加命令行参数检查.  
14.   
15.   
16. hexcvt ()  
17. {  
18.     if [ -z "$1" ]  
19.     then  
20.       echo 0  
21.       return    # 如果没有参数传递到这个函数中的话就"return" 0.  
22.     fi  
23.   
24.     echo ""$1" "$BASE" o p" | dc  
25.     #                 "o" 设置输出的基数(数制).  
26.     #                   "p" 打印栈顶.  
27.     # 参考dc的man页来了解其他的选项.  
28.     return  
29. }  
30.   
31. hexcvt "$1"  
32.   
33. exit 0  

结果：

[html]  view plain copy

1. root@ubuntu:~/resource/shell-study/0622-2013# chmod +x test14.sh   
2. root@ubuntu:~/resource/shell-study/0622-2013# ./test14.sh 15   
3. F  
4. root@ubuntu:~/resource/shell-study/0622-2013# ./test14.sh 6  
5. 6  
6. root@ubuntu:~/resource/shell-study/0622-2013# ./test14.sh 20  
7. 14  
8. root@ubuntu:~/resource/shell-study/0622-2013#   

awk

在脚本中使用浮点运算的另一种方法是使用awk内建的数学运算函数

看一个用awk的实例：计算直角三角形的斜边

[html]  view plain copy

1. #!/bin/bash  
2. # hypotenuse.sh: 返回直角三角形的斜边.  
3. #               (直角边长的平方和,然后对和取平方根)  
4.   
5. ARGS=2                # 需要将2个直角边作为参数传递进来.  
6. E_BADARGS=65          # 错误的参数值.  
7.   
8. if [ $# -ne "$ARGS" ] # 测试传递到脚本中的参数值.  
9. then  
10.   echo "Usage: `basename $0` side_1 side_2"  
11.   exit $E_BADARGS  
12. fi  
13.   
14.   
15. AWKSCRIPT=' { printf( "%3.7f\n", sqrt($1*$1 + $2*$2) ) } '  
16. #             命令 / 传递给awk的参数  
17.   
18.   
19. # 现在, 将参数通过管道传递给awk.  
20. echo -n "Hypotenuse of $1 and $2 = "  
21. echo $1 $2 | awk "$AWKSCRIPT"  
22.   
23. exit 0  

结果：

[html]  view plain copy

1. root@ubuntu:~/resource/shell-study/0622-2013# ./test16.sh   
2. Usage: test16.sh side_1 side_2  
3. root@ubuntu:~/resource/shell-study/0622-2013# ./test16.sh 4 6  
4. Hypotenuse of 4 and