求满足1-b≠1的最小IEEE浮点数b的技术咨询

我们以单精度IEEE 754浮点数为例推导（双精度逻辑一致，仅参数不同）：

1. IEEE单精度浮点数结构 #

单精度浮点数共32位，分为三部分：

• 符号位（S）：1位，0表示正数
• 指数位（E）：8位，实际指数 e = E - BIAS，其中BIAS=127
• 尾数位（M）：23位，规格化数的数值公式为：(-1)^S * (1 + M/2^23) * 2^e

浮点数1.0的表示为：S=0，E=127（e=0），M=0，对应数值1.0*2^0=1.0。

2. 关键：浮点数舍入规则 #

默认舍入规则为向最近的偶数舍入：

• 若计算结果落在两个相邻浮点数之间，舍入到更近的那个；
• 若刚好在两个数的中点，舍入到尾数最后一位为偶数的浮点数。

3. 寻找1.0的前一个相邻浮点数 #

我们需要找到比1.0小的最大可表示浮点数（记为F_prev）：

• 由于1.0的指数e=0，尾数M=0无法再减小，因此指数需减1：e=-1，对应E=127-1=126；
• 尾数取最大值（全1），即M=2^23-1；
• 计算数值：(1 + (2^23-1)/2^23) * 2^-1 = (2^24-1)/2^24 = 1 - 2^-24。

4. 确定舍入分界点 #

1.0和F_prev的中点为：
(1.0 + F_prev)/2 = (1 + 1 - 2^-24)/2 = 1 - 2^-25

根据舍入规则：

• 当1 - b > 1 - 2^-25时，结果离1.0更近，舍入后等于1.0；
• 当1 - b = 1 - 2^-25时，刚好在中点，因1.0的尾数最后一位是0（偶数），舍入后仍为1.0；
• 当1 - b < 1 - 2^-25时，结果离F_prev更近，舍入后等于F_prev，即1-b≠1。

5. 最小满足条件的b #

由上述结论，当b > 2^-25时，1-b≠1。我们需要找大于2^-25的最小正浮点数：

• 2^-25的表示为S=0，E=127-25=102，M=0；
• 大于它的最小浮点数为S=0，E=102，M=1，数值为(1 + 2^-23)*2^-25 = 2^-25 + 2^-48。

关于你提到的格式0|BIAS-1|111...1 #

这个格式对应的是F_prev（即1-b的结果，当b=2^-24时），而非b本身。b=2^-24确实满足1-b≠1，但它不是最小的满足条件的数——更小的b（如2^-25 + 2^-48）也能满足条件。

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