如何用Python基于大气噪声（MP3/WAV）实现真随机数生成器并转比特作密钥？

一、音频转比特的核心逻辑 #

不管是WAV还是MP3，本质是要先拿到PCM原始采样数据——这是音频的数字表示，每个采样点是固定位宽的整数（常见16位、24位）。转比特的完整流程：

1. 解码音频：WAV是未压缩格式，直接读取采样值；MP3需要解码为PCM（有损压缩格式需还原成原始采样数据）。
2. 采样值转二进制：将每个整数采样值转换为对应位宽的二进制字符串（补前导零保证位数一致），拼接后得到原始比特流。
3. 熵优化：大气噪声可能存在硬件采样偏移等微小偏差，需通过去偏算法（如von Neumann）生成无偏比特流，才能用作密码学安全的密钥。

二、Python实现步骤与代码 #

1. 依赖准备 #

安装处理音频的基础库：

pip install wave pydub numpy

注意：pydub依赖FFmpeg，需提前安装（Ubuntu用apt install ffmpeg，Windows下载后配置环境变量）。

2. 读取并解码音频到PCM采样 #

读取WAV文件

import wave
import numpy as np

def read_wav(file_path):
    with wave.open(file_path, 'rb') as wf:
        n_channels, sampwidth, _, n_frames = wf.getparams()[:4]
        frames = wf.readframes(n_frames)
        # 根据采样位宽匹配数据类型
        dtype_map = {2: np.int16, 3: np.int32, 4: np.int32}
        if sampwidth not in dtype_map:
            raise ValueError("不支持的采样位宽")
        samples = np.frombuffer(frames, dtype=dtype_map[sampwidth])
        # 立体声取单声道（左声道）
        if n_channels > 1:
            samples = samples[::n_channels]
        return samples

读取MP3文件（解码为PCM）

from pydub import AudioSegment

def read_mp3(file_path):
    audio = AudioSegment.from_mp3(file_path)
    samples = np.array(audio.get_array_of_samples())
    # 立体声取单声道
    if audio.channels > 1:
        samples = samples[::audio.channels]
    return samples

3. 采样数据转比特流并去偏 #

def samples_to_bits(samples, bit_width=16):
    bits = []
    mask = (1 << bit_width) - 1
    for sample in samples:
        # 处理负数补码，转成固定位宽的二进制字符串
        bin_str = format(sample & mask, f'0{bit_width}b')
        bits.extend(list(bin_str))
    return bits

def von_neumann_debias(bits):
    # 相邻比特不同则保留第一个，相同则丢弃，实现去偏
    debiased = []
    for i in range(0, len(bits)-1, 2):
        b1, b2 = bits[i], bits[i+1]
        if b1 != b2:
            debiased.append(b1)
    return debiased

4. 生成真随机数或流密码密钥 #

def generate_key(bits, key_length_bytes=32):
    required_bits = key_length_bytes * 8
    if len(bits) < required_bits:
        raise ValueError("比特流长度不足，无法生成指定长度密钥")
    key_bits = bits[:required_bits]
    # 按8位分组转成字节串
    key_bytes = bytes(int(''.join(key_bits[i:i+8]), 2) for i in range(0, required_bits, 8))
    return key_bytes

# 示例流程
if __name__ == "__main__":
    # 读取录制的大气噪声WAV文件
    samples = read_wav("atmospheric_noise.wav")
    # 转成原始比特流
    raw_bits = samples_to_bits(samples, bit_width=16)
    # 去偏得到无偏比特流
    debiased_bits = von_neumann_debias(raw_bits)
    # 生成256位（32字节）流密码密钥
    stream_key = generate_key(debiased_bits, key_length_bytes=32)
    print("流密码密钥（十六进制）:", stream_key.hex())
    
    # 生成32位随机整数
    random_int = int(''.join(debiased_bits[:32]), 2)
    print("随机整数:", random_int)

三、密码学安全关键注意事项 #

• 熵验证：生成的比特流必须通过NIST SP 800-22等标准熵测试，确保随机性满足密码学要求。可使用ent工具或cryptography库辅助验证。
• 密钥唯一性：流密码中同一个密钥绝对不能加密多个明文，否则会导致信息泄露。
• 采样选择：优先提取采样值的低位比特（如最低4-8位），这部分的随机性通常更高，能减少硬件采样的系统偏差。
• 音频纯净度：必须使用未经降噪、压缩处理的原始大气噪声录音，任何后期处理都会破坏随机性。

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