为什么十进制计算也需要精确

Souloss

公告

欢迎来到我的博客！这是一条示例公告

Learn More

标签

Souloss

公告

欢迎来到我的博客！这是一条示例公告

Learn More

标签

Souloss

公告

欢迎来到我的博客！这是一条示例公告

Learn More

标签

397 字

1 分钟

为什么十进制计算也需要精确

2023-08-07

技术科普

/

底层原理

在为什么浮点数计算不精确中，我们了解到 0.1 + 0.2 = 0.300000004。但如果是在金融场景，即使是 0.00000001 的误差也是不可接受的。

一、金融计算的特殊性#

1.1 浮点数的误差在金融场景的后果#

1
# 银行利息计算错误示例
2
balance = 1000.00
3
interest_rate = 0.03  # 3% 年利率
4

5
# 错误的浮点数计算
6
calculated_interest = balance * interest_rate
7
print(f"利息: {calculated_interest}")  # 可能输出 30.000000000004
8

9
# 错误累计
10
total = 0.0
11
for _ in range(1000000):
12
    total += 0.0000001
13

14
print(f"总计: {total}")  # 不是 0.1！

1.2 金融计算的要求#

场景	精度要求
银行利息	精确到分（0.01）
外汇交易	精确到小数点后 8 位
加密货币	精确到聪（satoshi）
股票交易	精确到 0.001 元

二、Decimal 的工作原理#

2.1 Decimal vs Float#

1
from decimal import Decimal, getcontext
2

3
# Decimal 使用字符串初始化
4
a = Decimal('0.1')
5
b = Decimal('0.2')
6
print(a + b)  # 0.3
7

8
# 浮点数初始化 Decimal（不推荐）
9
c = Decimal(0.1)  # 会带入浮点误差
10
print(c)  # 0.1000000000000000055511151231257827021181583404541015625

2.2 Decimal 的表示#

1
Decimal = 符号 × 系数 × 10^指数
2

3
例如：Decimal('3.14') = (+1) × 314 × 10^(-2)

2.3 可配置的精度#

1
from decimal import Decimal, getcontext
2

3
# 设置精度（默认 28 位）
4
getcontext().prec = 50
5

6
a = Decimal('1') / Decimal('3')
7
print(a)  # 0.33333333333333333333333333333333333333333333333333

三、金融计算的最佳实践#

3.1 用字符串创建 Decimal#

1
from decimal import Decimal, ROUND_HALF_UP
2

3
# 始终用字符串创建 Decimal
4
price = Decimal('19.99')
5
tax_rate = Decimal('0.08')
6
quantity = Decimal('3')
7

8
# 计算
9
subtotal = price * quantity
10
tax = subtotal * tax_rate
11
total = subtotal + tax
12

13
# 四舍五入到分
14
total = total.quantize(Decimal('0.01'), rounding=ROUND_HALF_UP)
15
print(f"总计: {total}")  # 64.77

3.2 避免浮点数混合#

1
# 错误：混合使用
2
result = Decimal('0.1') + 0.2  # 0.2 是浮点数！
3

4
# 正确：全部使用 Decimal
5
result = Decimal('0.1') + Decimal('0.2')
6

7
# 如果必须从浮点数转换
8
from decimal import Decimal
9
result = Decimal(str(0.2)) + Decimal('0.1')

3.3 货币计算示例#

1
from decimal import Decimal, ROUND_HALF_UP
2

3
class Money:
4
    def __init__(self, amount, currency='CNY'):
5
        self.amount = Decimal(str(amount))
6
        self.currency = currency
7

8
    def __add__(self, other):
9
        if self.currency != other.currency:
10
            raise ValueError("货币类型不匹配")
11
        return Money(self.amount + other.amount, self.currency)
12

13
    def __mul__(self, multiplier):
14
        return Money(self.amount * Decimal(str(multiplier)), self.currency)
15

16
    def round(self, places=2):
17
        return Money(
18
            self.amount.quantize(
19
                Decimal(10) ** -places,
20
                rounding=ROUND_HALF_UP
21
            ),
22
            self.currency
23
        )
24

25
# 使用
26
price = Money('19.99')
27
tax = price * Decimal('0.08')
28
total = (price + tax).round()
29
print(f"总价: {total.amount} {total.currency}")  # 21.59 CNY

四、Python 中的 Decimal 性能#

4.1 性能对比#

1
# Decimal 比 float 慢约 10-100 倍
2
# 因为是软件模拟，而非硬件支持
3

4
# benchmark
5
import timeit
6

7
float_time = timeit.timeit('0.1 + 0.2', number=100000)
8
decimal_time = timeit.timeit(
9
    "Decimal('0.1') + Decimal('0.2')",
10
    setup="from decimal import Decimal",
11
    number=100000
12
)
13

14
print(f"Float: {float_time:.4f}s")
15
print(f"Decimal: {decimal_time:.4f}s")
16
# Float: 0.0023s
17
# Decimal: 0.1847s