关于numerical_approx函数行为不符合预期的技术咨询

核心原因 #

你碰到的问题源于SymPy numerical_approx（即evalf）的双重逻辑：

• 当指定digits=2时，它会为底层数值设置对应两位十进制有效数字的二进制存储精度（约7位）。但512是2的整数次幂，在7位二进制精度下可以被精确存储，因此变量t的实际值仍是512，并未发生舍入。
• 打印输出时，SymPy会严格按照你指定的两位有效数字格式显示，将512格式化输出为510，但这只是显示层面的处理，不会改变内部存储的真实数值。这就导致print(t)显示510，而t-510计算结果为2.0。

获取可计算误差的近似值 #

如果你需要真正经过舍入的近似值来计算误差，可以用以下方法：

• 使用Python内置round函数：根据有效数字位数确定舍入位置。对于512（三位数字），两位有效数字需要舍入到十位，即：

  approx = round(512, -1)
  print(approx)  # 输出510
  print(512 - approx)  # 输出2.0
  

• SymPy内结合转换处理：将numerical_approx的结果转为Python浮点数后，再按有效数字舍入：

  from sympy import *
  t = 512.numerical_approx(digits=2)
  t_float = float(t)
  # 计算有效数字对应的舍入位数：总位数 - 有效数字位数
  round_pos = len(str(int(t_float))) - 2
  approx = round(t_float, -round_pos)
  print(approx)  # 输出510
  print(t_float - approx)  # 输出2.0
  

内容的提问来源于stack exchange，提问作者Ilan