Go語言浮點數(shù)的存儲方式

本文將為大家詳細介紹Go語言浮點數(shù)的存儲方式，內(nèi)容詳細步驟清晰，細節(jié)處理妥當，希望大家通過這篇文章有所收獲，我們先來看看浮點數(shù)如何在程序中被使用的：

創(chuàng)新互聯(lián)-專業(yè)網(wǎng)站定制、快速模板網(wǎng)站建設、高性價比奉節(jié)網(wǎng)站開發(fā)、企業(yè)建站全套包干低至880元,成熟完善的模板庫,直接使用。一站式奉節(jié)網(wǎng)站制作公司更省心,省錢,快速模板網(wǎng)站建設找我們，業(yè)務覆蓋奉節(jié)地區(qū)。費用合理售后完善，10余年實體公司更值得信賴。

下面的一段簡單程序 0.3 + 0.6 結(jié)果是什么？

1 var f1 float64 = 0.32 var f2 float64 = 0.63 fmt.Println(f1 + f2)

有人會天真的認為是0.9，但實際輸出卻是0.8999999999999999（go 1.13.5）

問題在于大多數(shù)小數(shù)表示成二進制之后是近似且無限的。
以0.1為例。它可能是你能想到的最簡單的十進制之一，但是二進制看起來卻非常復雜：0.0001100110011001100…
其是一串連續(xù)循環(huán)無限的數(shù)字（涉及到10進制轉(zhuǎn)換為2進制，暫不介紹）。
結(jié)果的荒誕性告訴我們，必須深入理解浮點數(shù)在計算機中的存儲方式及其性質(zhì)，才能正確處理數(shù)字的計算。
golang 與其他很多語言（C、C++、Python…）一樣，使用了IEEE-754標準存儲浮點數(shù)。

IEEE-754 如何存儲浮點數(shù)

IEEE-754規(guī)范使用特殊的以2為基數(shù)的科學表示法表示浮點數(shù)。

32位的單精度浮點數(shù) 與 64位的雙精度浮點數(shù)的差異

符號位：1 為 負數(shù)， 0 為正數(shù)。
指數(shù)位：存儲 指數(shù)加上偏移量，偏移量是為了表達負數(shù)而設計的。
小數(shù)位：存儲系數(shù)的小數(shù)位的準確或者最接近的值。

以 數(shù)字 0.085 為例。

小位數(shù)的表達方式

以0.36 為例:
010 1110 0001 0100 0111 1011 = 0.36 (第一位數(shù)字代表1/2,第二位數(shù)字是1/4…，0.36 是所有位相加)
分解后的計算步驟為:

go語言顯示浮點數(shù) - 驗證之前的理論

接下來用一個案例有助于我們理解并驗證IEEE-754 浮點數(shù)的表示方式。

math.Float32bits 可以為我們打印出32位數(shù)據(jù)的二進制表示。(注：math.Float64bits可以打印64位數(shù)據(jù)的二進制)

下面的go代碼將輸出0.085的浮點數(shù)二進制表達，并且為了驗證之前理論的正確性，根據(jù)二進制表示反向推導出其所表示的原始十進制0.085

 

輸出：表明我們對于浮點數(shù)的理解正確。

1 Starting Number: 0.0850002 Bit Pattern: 0 | 0111 1011 | 010 1110 0001 0100 0111 10113 Sign: 0 Exponent: 123 (-4) Mantissa: 0.360000 Value: 0.085000

經(jīng)典問題：如何判斷一個浮點數(shù)其實存儲的是整數(shù)

下面是一個有趣的問題，如何判斷一個浮點數(shù)其實存儲的是整數(shù)？

思考10秒鐘…

下面是一段判斷浮點數(shù)是否為整數(shù)的go代碼實現(xiàn)，我們接下來逐行分析函數(shù)。
它可以加深對于浮點數(shù)的理解

 

1、要保證是整數(shù)，一個重要的條件是必須要指數(shù)位大于127，如果指數(shù)位為127，代表指數(shù)為0. 指數(shù)位大于127，代表指數(shù)大于0， 反之小于0.

下面我們以數(shù)字234523為例子：

第一步,計算指數(shù)。由于 多減去了23，所以在第一個判斷中 判斷條件為 exponent < -23
exponent := int(bits >> 23) - bias - 23

第二步，
(bits & ((1 << 23) - 1)) 計算小數(shù)位。

| (1 << 23) 代表 將1加在前方。
1 + 小數(shù) = 系數(shù)。

 第三步，計算intTest 只有當指數(shù)的倍數(shù)可以彌補最小的小數(shù)位的時候，才是一個整數(shù)。
如下，指數(shù)是17位，其不能夠彌補最后6位的小數(shù)。即不能彌補1/2^18 的小數(shù)。
由于2^18位之后為0.所以是整數(shù)。

golang decimal 包詳解

要理解decimal包，首先需要知道兩個重要的概念，Normal number、denormal (or subnormal) number 以及精度。

1.概念：Normal number and denormal (or subnormal) number

wiki的解釋是：

 

什么意思呢？在IEEE-754中指數(shù)位有一個偏移量，偏移量是為了表達負數(shù)而設計的。比如單精度中的0.085，實際的指數(shù)是 -3， 存儲到指數(shù)位是123。
所以表達的負數(shù)就是有上限的。這個上限就是2^-126。如果比這個負數(shù)還要小，例如2^-127,這個時候應該表達為0.1 * 2 ^ -126. 這時系數(shù)變?yōu)榱瞬皇?為前導的數(shù)，這個數(shù)就叫做denormal (or subnormal) number。
正常的系數(shù)是以1為前導的數(shù)就叫做Normal number。

2.概念：精度

精度是一個非常復雜的概念，在這里筆者討論的是2進制浮點數(shù)的10進制精度。
精度為d表示的是在一個范圍內(nèi)，如果我們將d位10進制（按照科學計數(shù)法表達）轉(zhuǎn)換為二進制。再將二進制轉(zhuǎn)換為d位10進制。數(shù)據(jù)不損失意味著在此范圍內(nèi)是有d精度的。
精度的原因在于，數(shù)據(jù)在進制之間相互轉(zhuǎn)換時，是不能夠精準匹配的，而是匹配到一個最近的數(shù)。如圖所示：

                                          精度轉(zhuǎn)換

 

在這里暫時不深入探討，而是給出結(jié)論：（注：精度是動態(tài)變化的，不同的范圍可能有不同的精度。這是由于 2的冪 與 10的冪之間的交錯是不同的。）

float32的精度為6-8位，

float64的精度為15-17位

目前使用比較多的精準操作浮點數(shù)的decimal包是shopspring/decimal。鏈接:https://github.com/shopspring/decimal

decimal包使用math/big包存儲大整數(shù)并進行大整數(shù)的計算。

比如對于字符串 “123.45” 我們可以將其轉(zhuǎn)換為12345這個大整數(shù)，以及-2代表指數(shù)。參考decimal結(jié)構體：

在本文不會探討math/big是如何進行大整數(shù)運算的，而是探討decimal包一個非常重要的函數(shù)：

NewFromFloat(value float64) Decimal

其主要調(diào)用了下面的函數(shù)：

 

此函數(shù)會將浮點數(shù)轉(zhuǎn)換為Decimal結(jié)構。
讀者想象一下這個問題：如果存儲到浮點數(shù)中的值（例如0.1）本身就是一個近似值，為什么decimal包能夠解決計算的準確性？
原因在于，deciimal包可以精準的將一個浮點數(shù)轉(zhuǎn)換為10進制。這就是NewFromFloat為我們做的事情。
下面我將對此函數(shù)做逐行分析。

 

第5行：剝離出IEEE浮點數(shù)的指數(shù)位
exp := int(bits>>flt.mantbits) & (1<

第6行：剝離出浮點數(shù)的系數(shù)的小數(shù)位
mant := bits & (uint64(1)<

第7行：如果是指數(shù)位為0，代表浮點數(shù)是denormal (or subnormal) number；
默認情況下會在mant之前加上1，因為mant只是系數(shù)的小數(shù)，在前面加上1后，代表真正的小數(shù)位。
現(xiàn)在 mant = IEEE浮點數(shù)系數(shù) * 2^53

第13行：加上偏移量，exp現(xiàn)在代表真正的指數(shù)。
第14行：引入了一個中間結(jié)構decimal

第15行：調(diào)用d.Assign(mant) , 將mant作為10進制數(shù)，存起來。
10進制數(shù)的每一位都作為一個字符存儲到 decimal的byte數(shù)組中

第16行：調(diào)用shift函數(shù)，這個函數(shù)非常難理解。

 

此函數(shù)的功能是為了獲取此浮點數(shù)代表的10進制數(shù)據(jù)的整數(shù)位個數(shù)以及小數(shù)位個數(shù)，此函數(shù)的完整證明附后。（注1)
exp是真實的指數(shù)，其也是能夠覆蓋小數(shù)部分2進制位的個數(shù)。（參考前面如何判斷浮點數(shù)是整數(shù)）
exp - int(flt.mantbits)代表不能被exp覆蓋的2進制位的個數(shù)
如果exp - int(flt.mantbits) > 0 代表exp能夠完全覆蓋小數(shù)位 因此 浮點數(shù)是一個非常大的整數(shù)，這時會調(diào)用leftShift(a, uint(k))。否則將調(diào)用rightShift(a, uint(-k)), 常規(guī)rightShift會調(diào)用得更多。因此我們來看看rightShift函數(shù)的實現(xiàn)。

第5行：此for循環(huán)將計算浮點數(shù)10進制表示的小數(shù)部分的有效位為 r-1 。
n >> k 是一個重要的衡量指標，代表了小數(shù)部分與整數(shù)部分的分割。此函數(shù)的完整證明附后。（注1)

第21行：此時整數(shù)部分所占的有效位數(shù)為a.dp -=（r-1）
第24行：這兩個循環(huán)做了2件事情：
1、計算10進制表示的有效位數(shù)
2、將10進制表示存入bytes數(shù)組中。例如對于浮點數(shù)64.125，現(xiàn)在byte數(shù)組存儲的前5位就是64125

 

繼續(xù)回到newFromFloat函數(shù)，第18行，調(diào)用了roundShortest函數(shù)，
此函數(shù)非常關鍵。其會將浮點數(shù)轉(zhuǎn)換為離其最近的十進制數(shù)。
這是為什么decimal.NewFromFloat(0.1)能夠精準表達0.1的原因。

參考上面的精度，此函數(shù)主要考察了2的冪與10的冪之間的交錯關系。四舍五入到最接近的10進制值。
此函數(shù)實質(zhì)實現(xiàn)的是Grisu3 算法,有想深入了解的可以去看看論文。筆者在這里提示幾點：
1、2^exp <= d < 10^dp。
2、10進制數(shù)之間至少相聚10^(dp-nd)
3、2的冪之間的最小間距至少為2^(exp-mantbits)
4、什么時候d就是最接近2進制的10進制數(shù)？
如果10^(dp-nd) > 2^(exp-mantbits)，表明 當十進制下降一個最小位數(shù)時，匹配到的是更小的數(shù)字value - 2^(exp-mantbits)，所以d就是最接近浮點數(shù)的10進制數(shù)。

 

 

繼續(xù)回到newFromFloat函數(shù)，第19行 如果精度小于19，是位于int64范圍內(nèi)的，可以使用快速路徑，否則使用math/big包進行賦值操作，效率稍微要慢一些。
第36行，正常情況幾乎不會發(fā)生。如果setstring在異常的情況下會調(diào)用NewFromFloatWithExponent 指定精度進行四舍五入截斷。

注一：快速的獲取一個浮點數(shù)代表的十進制

以典型的數(shù)字64.125 為例 ， 它可以被浮點數(shù)二進制精準表達為：
Bit Patterns: 0 | 10000000101 | 0000000010000000000000000000000000000000000000000000
Sign: 0 | Exponent: 1029 (6) | Mantissa: 0.001953

即 64.125 = 1.001953125 * 2^6
注意觀察浮點數(shù)的小數(shù)位在第九位有1, 代表2^-9 即 0.001953125.

我們在浮點數(shù)的小數(shù)位前 附上數(shù)字1，10000000010000000000000000000000000000000000000000000 代表其為1 / 2^0 .

此時我們可以認為這個數(shù)代表的是1.001953125. 那么這樣長的二進制數(shù)變?yōu)?0進制又是多少呢:4512395720392704。

即 1.001953125 = 4512395720392704 * 2^(-52)

所以64.125 = 4512395720392704 * 2^(-52) * 2^6 = 4512395720392704 * 2^(-46)
在這里，有一種重要的等式。即 (2 ^ -46) 等價于向左移動了46位。并且移動后剩下的部分即為64,而舍棄的部分其實是小數(shù)部分0.125。
這個等式看似復雜其實很好證明，即第46位其實代表的是2^45。其除以2^46后是一個小數(shù)。依次類推…

因此對于數(shù)字 4512395720392704 ， 我們可以用4，45，451，4512 … 依次除以 2 ^ 46. 一直到找到數(shù)451239572039270 其除以2^46不為0。這個不為0的數(shù)一定為6。
接著我們保留后46位，其實是保留了小數(shù)位。

假設 4512395720392704 / 2^46 = (6 + num)
64.125 =(6 + num) * 10 + C = 60 + 10* num + C

當我們將通過位運算保留后46位，設為A, 則 A / 2^46 = num
所以 (A * 10 + C) / 2 ^46 =(num * 10 +C) = 4.125
此我們又可以把4提取出來。

關于Go語言浮點數(shù)的存儲方式就分享到這里了，希望以上內(nèi)容可以對大家有一定的幫助，可以學到更多知識。如果覺得文章不錯，可以把它分享出去讓更多的人看到。

名稱欄目：Go語言浮點數(shù)的存儲方式
文章路徑：http://m.jiaotiyi.com/article/jhppjo.html