如何在Brian2中对Quantity类型实例执行四舍五入操作？

• 原生Pythonround函数未做Brian2Quantity类型适配，直接调用会触发类型错误。
• 直接调用np.round得到0*mV的核心原因是：Brian2的Quantity底层统一存储SI基准单位的数值，例如77*mV底层实际存储值为0.077（对应伏特单位），无参数调用np.round时会对底层的0.077做整数取整得到0，转换回mV单位后就输出0*mV，和预期不符。

优先推荐剥离单位运算的方案，不会出现单位换算逻辑错误，不需要记忆SI单位量级：

方案1：剥离单位做数值运算后再绑定单位 #

Quantity除以对应单位后会得到无单位的纯数值，此时可以正常使用所有Python、NumPy的数值运算函数，运算完成后乘回单位即可得到合法的Quantity实例：

from brian2 import mV
import numpy as np

equilibrium_potential = my_func()
# 提取mV单位下的纯数值
ep_pure = equilibrium_potential / mV
# 四舍五入取整，也可以给round传参指定保留小数位
rounded_ep = round(ep_pure) * mV
# 若需保留1位小数，替换为numpy版本即可
# rounded_ep = np.round(ep_pure, decimals=1) * mV

方案2：直接对Quantity调用np.round（不推荐） #

如果不想剥离单位，需要根据目标精度对应SI单位的量级传入decimals参数，例如要取整到1mV精度（1mV = 1e-3 V，即SI基准单位下的1e-3量级），需要指定保留3位小数：

# 取整到1mV精度
rounded_ep = np.round(equilibrium_potential, decimals=3)

该方案需要手动计算不同单位对应的小数位，遇到μV、nS这类小单位时非常容易写错精度，非必要不使用。

判断电位是否在77*mV左右不需要先做四舍五入，直接设置合理容差做差值判断即可，灵活度更高：

# 设定允许的误差范围，例如上下1mV
tolerance = 1 * mV
if abs(equilibrium_potential - 77*mV) <= tolerance:
    # 数值符合预期的处理逻辑
    pass

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