如何将十进制数转换为IEEE754格式？求ibm2float32反向实现方法

IEEE 754单精度浮点数是32位二进制格式，结构分为三部分：

• 符号位（第31位）：0代表正数，1代表负数
• 指数位（第30-23位，共8位）：采用偏移127的表示方式，存储的指数值 = 实际指数 + 127
• 尾数位（第22-0位，共23位）：存储规格化后有效数字的小数部分，隐含一个最高位的1（规格化后有效数字形式为1.xxxxx...）

转换步骤：

1. 将目标十进制数转换为二进制数，区分正负
2. 对二进制数进行规格化，得到1.xxxxx × 2^e的形式（e为实际指数）
3. 计算存储指数：存储指数 = e + 127，转为8位二进制
4. 截取有效数字的小数部分，保留23位（不足补0，超出则截断或按规则舍入）
5. 按「符号位+指数位+尾数位」的顺序拼接成32位二进制，再转为十六进制或无符号整数

IBM单精度格式为32位结构：1位符号位 + 7位指数位（偏移量64，基数为16） + 24位尾数位（形式为0.xxxxx，无隐含最高位）。以下是反向转换的实现代码：

import numpy as np

def float32_to_ibm(f):
    # 处理0值
    if f == 0.0:
        return np.uint32(0)
    
    # 确定符号位
    sign_bit = 0 if f > 0 else 1
    abs_f = abs(f)
    
    # 规格化为 0.m × 16^e，其中 0.0625 ≤ m < 1
    exponent = 0
    while abs_f >= 1:
        abs_f /= 16
        exponent += 1
    while abs_f < 0.0625:
        abs_f *= 16
        exponent -= 1
    
    # 计算偏移后的指数（7位，范围0-127）
    shifted_exponent = exponent + 64
    if shifted_exponent < 0 or shifted_exponent > 127:
        raise ValueError("数值超出IBM单精度浮点数的表示范围")
    
    # 计算24位尾数：将0.m乘以2^24后取整
    mantissa = int(abs_f * (1 << 24))
    # 确保尾数仅保留24位
    mantissa &= 0xFFFFFF
    
    # 拼接成32位无符号整数
    ibm_value = (sign_bit << 31) | (shifted_exponent << 24) | mantissa
    return np.uint32(ibm_value)

# 测试示例：反向验证原转换
test_num = ibm2float32(np.uint32(0xc1180000))
print(f"原转换结果：{test_num}")
# 转回IBM格式
converted_back = float32_to_ibm(test_num)
print(f"转回后的十六进制：{hex(converted_back)}")  # 输出应为0xc1180000

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