JavaScript 浮點數(shù)陷阱及解法

JavaScript 浮點數(shù)陷阱及解法

眾所周知兜蠕,JavaScript 浮點數(shù)運算時經(jīng)常遇到會 0.0000000010.999999999 這樣奇怪的結(jié)果,如 0.1+0.2=0.30000000000000004锋拖、1-0.9=0.09999999999999998,很多人知道這是浮點數(shù)誤差問題,但具體就說不清楚了燕差。本文幫你理清這背后的原理以及解決方案潜叛,還會向你解釋JS中的大數(shù)危機和四則運算中會遇到的坑秽褒。

浮點數(shù)的存儲

首先要搞清楚 JavaScript 如何存儲小數(shù)。和其它語言如 Java 和 Python 不同威兜,JavaScript 中所有數(shù)字包括整數(shù)和小數(shù)都只有一種類型 — Number销斟。它的實現(xiàn)遵循 IEEE 754 標(biāo)準(zhǔn),使用 64 位固定長度來表示椒舵,也就是標(biāo)準(zhǔn)的 double 雙精度浮點數(shù)(相關(guān)的還有float 32位單精度)蚂踊。計算機組成原理中有過詳細介紹,如果你不記得也沒關(guān)系笔宿。

這樣的存儲結(jié)構(gòu)優(yōu)點是可以歸一化處理整數(shù)和小數(shù)悴势,節(jié)省存儲空間。

64位比特又可分為三個部分:

  • 符號位S:第 1 位是正負數(shù)符號位(sign)措伐,0代表正數(shù)特纤,1代表負數(shù)
  • 指數(shù)位E:中間的 11 位存儲指數(shù)(exponent),用來表示次方數(shù)
  • 尾數(shù)位M:最后的 52 位是尾數(shù)(mantissa)侥加,超出的部分自動進一舍零

64 bit allocation

實際數(shù)字就可以用以下公式來計算:

$ V = (-1)^{S}\times M \times 2^{E} $

注意以上的公式遵循科學(xué)計數(shù)法的規(guī)范捧存,在十進制是為0<M<10,到二進行就是0<M<2担败。也就是說整數(shù)部分只能是1昔穴,所以可以被舍去,只保留后面的小數(shù)部分提前。如 4.5 轉(zhuǎn)換成二進制就是 100.1吗货,科學(xué)計數(shù)法表示是 1.001*2^2,舍去1后 M = 001狈网。E是一個無符號整數(shù)宙搬,因為長度是11位,取值范圍是 0~2047拓哺。但是科學(xué)計數(shù)法中的指數(shù)是可以為負數(shù)的勇垛,所以再減去一個中間數(shù) 1023,[0,1022]表示為負士鸥,[1024,2047] 表示為正闲孤。如4.5 的指數(shù)E = 1025,尾數(shù)M為 001烤礁。

最終的公式變成:

$ V = (-1)^{S}\times (M+1) \times 2^{E-1023} $

所以 4.5 最終表示為(M=001讼积、E=1025):

4.5 allocation map

(圖片由此生成 http://www.binaryconvert.com/convert_double.html)

下面再以 0.1 例解釋浮點誤差的原因肥照, 0.1 轉(zhuǎn)成二進制表示為 0.0001100110011001100(1100循環(huán)),1.100110011001100x2^-4勤众,所以 E=-4+1023=1019建峭;M 舍去首位的1,得到 100110011...决摧。最終就是:

0.1 allocation map


轉(zhuǎn)化成十進制后為 0.100000000000000005551115123126亿蒸,因此就出現(xiàn)了浮點誤差。

為什么 0.1+0.2=0.30000000000000004掌桩?

計算步驟為:

// 0.1 和 0.2 都轉(zhuǎn)化成二進制后再進行運算
0.00011001100110011001100110011001100110011001100110011010 +
0.0011001100110011001100110011001100110011001100110011010 =
0.0100110011001100110011001100110011001100110011001100111

// 轉(zhuǎn)成十進制正好是 0.30000000000000004

為什么 x=0.1 能得到 0.1边锁?

恭喜你到了看山不是山的境界。因為 mantissa 固定長度是 52 位波岛,再加上省略的一位茅坛,最多可以表示的數(shù)是 2^53=9007199254740992,對應(yīng)科學(xué)計數(shù)尾數(shù)是 9.007199254740992则拷,這也是 JS 最多能表示的精度贡蓖。它的長度是 16,所以可以使用 toPrecision(16) 來做精度運算煌茬,超過的精度會自動做湊整處理斥铺。于是就有:

0.10000000000000000555.toPrecision(16)
// 返回 0.1000000000000000,去掉末尾的零后正好為 0.1

// 但你看到的 `0.1` 實際上并不是 `0.1`坛善。不信你可用更高的精度試試:
0.1.toPrecision(21) = 0.100000000000000005551

大數(shù)危機

可能你已經(jīng)隱約感覺到了晾蜘,如果整數(shù)大于 9007199254740992 會出現(xiàn)什么情況呢?
由于 M 最大值是 1023眠屎,所以最大可以表示的整數(shù)是 2^1024 - 1剔交。這就是能表示的最大整數(shù)。但你并不能這樣計算這個數(shù)字改衩,因為從 2^1024 開始就變成了 Infinity

> Math.pow(2, 1023)
8.98846567431158e+307

> Math.pow(2, 1024)
Infinity

那么對于 (2^53, 2^63) 之間的數(shù)會出現(xiàn)什么情況呢岖常?

  • (2^53, 2^54) 之間的數(shù)會兩個選一個,只能精確表示偶數(shù)
  • (2^54, 2^55) 之間的數(shù)會四個選一個葫督,只能精確表示4個倍數(shù)
  • ... 依次跳過更多2的倍數(shù)

下面這張圖能很好的表示 JavaScript 中浮點數(shù)和實數(shù)(Real Number)之間的對應(yīng)關(guān)系竭鞍。我們常用的 (-2^53, 2^53) 只是最中間非常小的一部分,越往兩邊越稀疏越不精確候衍。

floating-point numbers map to real numbers

在淘寶早期的訂單系統(tǒng)中把訂單號當(dāng)作數(shù)字處理笼蛛,后來隨意訂單號暴增洒放,已經(jīng)超過了
9007199254740992蛉鹿,最終的解法是把訂單號改成字符串處理。

要想解決大數(shù)的問題你可以引用第三方庫 bignumber.js往湿,原理是把所有數(shù)字當(dāng)作字符串妖异,重新實現(xiàn)了計算邏輯惋戏,缺點是性能比原生的差很多。所以原生支持大數(shù)就很有必要了他膳,現(xiàn)在 TC39 已經(jīng)有一個 Stage 3 的提案 proposal bigint响逢,大數(shù)問題有問徹底解決。

toPrecision vs toFixed

數(shù)據(jù)處理時棕孙,這兩個函數(shù)很容易混淆舔亭。它們的共同點是把數(shù)字轉(zhuǎn)成字符串供展示使用。注意在計算的中間過程不要使用蟀俊,只用于最終結(jié)果钦铺。

不同點就需要注意一下:

  • toPrecision 是處理精度,精度是從左至右第一個不為0的數(shù)開始數(shù)起肢预。
  • toFixed 是小數(shù)點后指定位數(shù)取整矛洞,從小數(shù)點開始數(shù)起。

兩者都能對多余數(shù)字做湊整處理烫映,也有些人用 toFixed 來做四舍五入沼本,但一定要知道它是有 Bug 的。

如:1.005.toFixed(2) 返回的是 1.00 而不是 1.01锭沟。

原因: 1.005 實際對應(yīng)的數(shù)字是 1.00499999999999989抽兆,在四舍五入時全部被舍去!

解法:使用專業(yè)的四舍五入函數(shù) Math.round() 來處理族淮。但 Math.round(1.005 * 100) / 100 還是不行郊丛,因為 1.005 * 100 = 100.49999999999999。還需要把乘法和除法精度誤差都解決后再使用 Math.round瞧筛±魇欤可以使用后面介紹的 number-precision#round 方法來解決。

解決方案

回到最關(guān)心的問題:如何解決浮點誤差较幌。首先揍瑟,理論上用有限的空間來存儲無限的小數(shù)是不可能保證精確的,但我們可以處理一下得到我們期望的結(jié)果乍炉。

數(shù)據(jù)展示類

當(dāng)你拿到 1.4000000000000001 這樣的數(shù)據(jù)要展示時绢片,建議使用 toPrecision 湊整并 parseFloat 轉(zhuǎn)成數(shù)字后再顯示,如下:

parseFloat(1.4000000000000001.toPrecision(12)) === 1.4  // True

封裝成方法就是:

function strip(num, precision = 12) {
  return +parseFloat(num.toPrecision(precision));
}

為什么選擇 12 做為默認精度岛琼?這是一個經(jīng)驗的選擇底循,一般選12就能解決掉大部分0001和0009問題,而且大部分情況下也夠用了槐瑞,如果你需要更精確可以調(diào)高熙涤。

數(shù)據(jù)運算類

對于運算類操作,如 +-*/,就不能使用 toPrecision 了祠挫。正確的做法是把小數(shù)轉(zhuǎn)成整數(shù)后再運算那槽。以加法為例:

/**
 * 精確加法
 */
function add(num1, num2) {
  const num1Digits = (num1.toString().split('.')[1] || '').length;
  const num2Digits = (num2.toString().split('.')[1] || '').length;
  const baseNum = Math.pow(10, Math.max(num1Digits, num2Digits));
  return (num1 * baseNum + num2 * baseNum) / baseNum;
}

以上方法能適用于大部分場景。遇到科學(xué)計數(shù)法如 2.3e+1(當(dāng)數(shù)字精度大于21時等舔,數(shù)字會強制轉(zhuǎn)為科學(xué)計數(shù)法形式顯示)時還需要特別處理一下骚灸。

能讀到這里,說明你非常有耐心慌植,那我就放個福利吧甚牲。
遇到浮點數(shù)誤差問題時可以直接使用

https://github.com/dt-fe/number-precision

完美支持浮點數(shù)的加減乘除、四舍五入等運算蝶柿。非常小只有1K鳖藕,遠小于絕大多數(shù)同類庫(如Math.js、BigDecimal.js)只锭,100%測試全覆蓋著恩,代碼可讀性強,不妨在你的應(yīng)用里用起來蜻展!

參考

數(shù)學(xué)公式無法正常顯示喉誊,請查看原文:camsong/blog 轉(zhuǎn)載請注明出處。
想不想加入阿里巴巴一起玩 ES7纵顾,React伍茄,F(xiàn)RP 等最新技術(shù),歡迎簡歷到 neosoyn@gmail.com

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末施逾,一起剝皮案震驚了整個濱河市敷矫,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌汉额,老刑警劉巖曹仗,帶你破解...
    沈念sama閱讀 212,454評論 6 493
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場離奇詭異蠕搜,居然都是意外死亡怎茫,警方通過查閱死者的電腦和手機,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 90,553評論 3 385
  • 文/潘曉璐 我一進店門妓灌,熙熙樓的掌柜王于貴愁眉苦臉地迎上來轨蛤,“玉大人,你說我怎么就攤上這事虫埂∠樯剑” “怎么了?”我有些...
    開封第一講書人閱讀 157,921評論 0 348
  • 文/不壞的土叔 我叫張陵掉伏,是天一觀的道長缝呕。 經(jīng)常有香客問我澳窑,道長,這世上最難降的妖魔是什么岳颇? 我笑而不...
    開封第一講書人閱讀 56,648評論 1 284
  • 正文 為了忘掉前任照捡,我火速辦了婚禮颅湘,結(jié)果婚禮上话侧,老公的妹妹穿的比我還像新娘。我一直安慰自己闯参,他們只是感情好瞻鹏,可當(dāng)我...
    茶點故事閱讀 65,770評論 6 386
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著鹿寨,像睡著了一般新博。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上脚草,一...
    開封第一講書人閱讀 49,950評論 1 291
  • 那天赫悄,我揣著相機與錄音,去河邊找鬼馏慨。 笑死埂淮,一個胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的写隶。 我是一名探鬼主播倔撞,決...
    沈念sama閱讀 39,090評論 3 410
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼慕趴!你這毒婦竟也來了痪蝇?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 37,817評論 0 268
  • 序言:老撾萬榮一對情侶失蹤冕房,失蹤者是張志新(化名)和其女友劉穎躏啰,沒想到半個月后,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體耙册,經(jīng)...
    沈念sama閱讀 44,275評論 1 303
  • 正文 獨居荒郊野嶺守林人離奇死亡丙唧,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 36,592評論 2 327
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發(fā)現(xiàn)自己被綠了觅玻。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片想际。...
    茶點故事閱讀 38,724評論 1 341
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖溪厘,靈堂內(nèi)的尸體忽然破棺而出胡本,到底是詐尸還是另有隱情,我是刑警寧澤畸悬,帶...
    沈念sama閱讀 34,409評論 4 333
  • 正文 年R本政府宣布侧甫,位于F島的核電站,受9級特大地震影響,放射性物質(zhì)發(fā)生泄漏披粟。R本人自食惡果不足惜咒锻,卻給世界環(huán)境...
    茶點故事閱讀 40,052評論 3 316
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望守屉。 院中可真熱鬧惑艇,春花似錦、人聲如沸拇泛。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,815評論 0 21
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽俺叭。三九已至恭取,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間熄守,已是汗流浹背蜈垮。 一陣腳步聲響...
    開封第一講書人閱讀 32,043評論 1 266
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留裕照,地道東北人攒发。 一個月前我還...
    沈念sama閱讀 46,503評論 2 361
  • 正文 我出身青樓,卻偏偏與公主長得像牍氛,于是被迫代替她去往敵國和親晨继。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點故事閱讀 43,627評論 2 350

推薦閱讀更多精彩內(nèi)容