IEEE 754双精度浮点数区间插入N个均匀分布值的可行性判断

核心问题 #

要按公式 A + (B - A)/(N+1)*X（X∈[1,N]）生成N个严格递增、互不相同的双精度浮点数，关键是确保步长 step = (B - A)/(N+1) 足够大，避免舍入误差导致相邻计算值重合或逆序。

可靠判断逻辑 #

双精度浮点数的最小可分辨间隔（ULP，最后一位的单位）随数值变化：

• 正规数（绝对值≥2^(-1022)）：对于数x，ULP(x) = 2^(floor(log₂(|x|)) - 52)
• 亚正规数（绝对值<2^{(-1022)）：ULP固定为2}(-1074)

要保证生成的序列严格递增，需满足步长step > 区间内所有数的ULP的1/2。由于ULP在正规数区间随数值增大而递增，只需验证区间两端的条件即可：

1. 若B == A，直接判定无法划分（无中间数）
2. 计算step = (B - A)/(N+1)，若step == 0（下溢为0），判定无法划分
3. 对A和B分别验证：
   • 若为正规数：检查step > ULP(该数)/2
   • 若为亚正规数：检查step > 2^(-1075)（即亚正规ULP的1/2）
     两个端点条件都满足时，即可保证序列严格递增。

实用简化方案 #

如果不需要极致严谨的判断，可采用以下简化阈值：

• 对于常规数值区间（绝对值在1e-300到1e300之间的正规数）：确保 step > 2^(-50)（约8.88e-16），该值远大于常规数值的ULP（比如x=1时ULP为2^(-52)≈2.22e-16），留有足够余量避免舍入问题
• 对于亚正规数区间：确保 step > 2^(-1073)（约1.18e-323），即亚正规ULP的2倍，可保证相邻值可区分

注意事项 #

• 避免直接通过“(B-A)/N > 固定阈值”判断，因为不同数值范围的ULP差异极大（比如1e300的ULP约为1e247，而1e-300的ULP约为1e-352），固定阈值无法覆盖所有场景
• 若step接近区间内的ULP，建议改用递推方式生成序列并实时验证 x_{k+1} > x_k，确保严格递增

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