Re: [问题] shell浮点数运算问题

楼主: rickieyang (Rickie Yang)   2018-12-04 23:26:17
忍不住回一下...
※ 引述《zchien (小建)》之铭言:
: 我用unix 系统
UNIX的商标权由国际开放标准组织所拥有
只有符合单一UNIX规范的UNIX系统才能使用UNIX这个名称
否则只能称为类UNIX(UNIX-like)。
既然你来到这个版, 我猜你是用 Linux, 就直接说用 Linux 吧
: 写sh时,遇到了一个bug ,就是浮点数无法运算
我猜你的意思是 写 shell script.
: 我使用ntpdate去校正时间, 并且记录校正了多少时间
: 例如:需要校正的时间是 0.000304 s
: 我写
: ntpdate -p 10.37.82.23 > time.txt
不知道你是用哪一个系统?
Linux 用 ntpdate -p 10.37.82.23 应该会报错吧?!
usage: ntpdate [-46bBdqsuv] [-a key#] [-e delay] [-k file] [-p samples] \
[-o version#] [-t timeo] server ...
或是你是指 ntpdate -q 10.37.82.23 ?
: NTP= grep delay time.txt | awk ' ' {print $6}
上面的语法应该也会报错吧?
NTP=$(grep delay time.txt | awk '{print $6}')
: # $NTP will show 类似 0.000304
: if [ $NTP -le 0.100001 ];then
: ntpdate -u 10.37.82.23
: fi
先不管 能不能用小数点的问题.
为什么 小于等于 0.100001 秒才更新?
如果 offset 大于 0.100001 秒就不更新?
: # 如果NTP server 连线异常, $NTP will show 0.000000
: if [ $NTP -eq 0.000000 ];then
: echo "NTP server not connect"
: fi
: 上面第一个if 会 不成立 , 可能类似浮点数运算, 他以为两个都是0
: 第二个if 也会显示错误, [0.000000 -eq 0.000000 ] 他竟然显示不相等
: 请问这是什么原因,我可以怎么修正
应该不会是不成立或是不相等吧?
-eq 应该会直接骂人, 类似:
./test.sh: line n: [: 0.000000: integer expression expected
首先, 既然你是用 > time.txt, 会一直被覆蓋, 所以假设 time.txt 并没有要保留.
那其实可以直接一路 | 到底
NTP=$(ntpdate -q 10.37.82.23 | grep delay | awk '{print $6}')
然后都已经用 awk 了, 可以不用特别再叫 grep 出来跑龙套
NTP=$(ntpdate -q 10.37.82.23 | awk '/delay/{print $6}')
另外, 不知道你有没有发现, 这样你得到的 NTP 后面其实会有一个 ,
$ ntpdate -q clock.stdtime.gov.tw
server 211.22.103.158, stratum 2, offset 0.000697, delay 0.03252
4 Dec 22:41:21 ntpdate[1947]: adjust time server 211.22.103.158 offset 0.000697 sec
这应该会影响你后面的判断式
把它去掉比较好一些
NTP=$(ntpdate -q 10.37.82.23 | awk '/delay/{print $6}'|sed -e 's/,$//')
或是可以在 awk print 时先乘 1000000, 然后往后都可以用 us(微秒) 来处理
这样 awk 会自己把 , 去掉, 也间接避开了小数点运算的限制
NTP=$(ntpdate -q 10.37.82.23 | awk '/delay/{print $6*1000000}')
if [ "$NTP" -le "100001" ];then
ntpdate -u 10.37.82.23
fi
if [ "$NTP" -eq "0" ];then
echo "NTP server not connect"
fi
先不管不能判断小数点, 逻辑判断上 0.000000 也是 -le 0.100001 呀...
把两个 if 合在一起吧.
NTP=$(ntpdate -q 10.37.82.23 | awk '/delay/{print $6*1000000}')
if [ "$NTP" -eq "0" ];then
echo "NTP server not connect"
elif [ "$NTP" -le "100001" ];then
ntpdate -u 10.37.82.23
fi
虽然机率不比民乐透高, 有考虑过如果刚刚好 offset 是 0.000000 秒吗?
man ntpdate:
DIAGNOSTICS
ntpdate’s exit status is zero if it finds a server and updates the
clock, and nonzero otherwise.
ntpdate 有错误时的 return code 不是 0, 建议用这个来判断比较好.
#!/bin/sh
ntp_response=$(ntpdate -q 10.37.82.23 2>&1)
rc=$?
if [ "${rc}" -ne "0" ];then
echo "ntpdate error, return code: ${rc}"
else
ntp_offset=$(echo ${ntp_response} | awk '/delay/{print $6*1000000}')
if [ "${ntp_offset}" -le "100001" ];then
echo "offset: ${ntp_offset}"
ntpdate -u 10.37.82.23
fi
fi
大概是酱~
最后有两个问题还是不解.
1. 为什么 offset -le 0.100001 才要 -u, 超过 0.100001 不处理吗?
2. 有考虑直接改用 ntpd 吗?
作者: kdjf (我抓得到什么呢?)   2018-12-05 00:57:00
看以前发的文很可能是solaris系列,所以参数之类不同不奇怪
作者: zchien (小建)   2018-12-05 23:35:00
感谢回复,solaris和你上面指令应该相同,很多指令都是硬记的,所以上面可能写的不完全,明天上班,我在试看看,感恩Server上面有跑程式,如果时间瞬间差太大,程式会当机,所以写crontab 让他每小时校时
作者: kdjf (我抓得到什么呢?)   2018-12-06 01:22:00
那就更适合ntpd,他本身就考虑到时间不能突然跳动的问题了只是你系统上不一定有现成的
作者: rexsony (雷克斯索尼)   2018-12-06 18:18:00
喔喔,高手过招果然很多学习的主体地方
作者: jamselee (jamselee)   2018-12-10 07:16:00
版上高手好多
作者: TPPCMAN (e?▄?▄?▄?▄?▄)   2018-12-12 11:31:00
这是 高手高手高高手了吧

Links booklink

Contact Us: admin [ a t ] ucptt.com