如何在Octave中通过fread读取17位小端无符号整数？

Octave的fread不像Matlab那样直接支持'bit17'这种自定义位宽的格式，所以我们需要手动解析文件的比特流来模拟Matlab的fread(fd, K, 'bit17', 'ieee-le')行为。下面是具体的实现步骤和代码：

核心思路 #

Matlab的'bit17'+'ieee-le'是按**LSB优先（小端比特顺序）**读取连续的比特流，每17位组成一个无符号整数。我们需要：

1. 读取足够的字节来覆盖K个17位整数的总比特数
2. 将字节转换为LSB优先的比特流
3. 把比特流分割为K个17位组，每个组转换为无符号整数

实现代码 #

高效向量化版本 #

function data = read_bit17_le(fd, K)
    total_bits = 17 * K;
    total_bytes = ceil(total_bits / 8);
    
    % 读取所有需要的字节
    bytes = fread(fd, total_bytes, 'uint8');
    
    % 将字节转换为LSB优先的比特流（向量化操作，避免循环提升效率）
    % 对每个字节，提取1到8位（对应LSB到MSB），再展平为一维数组
    bits = reshape(bitget(bytes, repmat(1:8, numel(bytes), 1))', 1, []);
    
    % 截取刚好需要的比特数，去掉末尾多余的填充比特
    bits = bits(1:total_bits);
    
    % 把比特流重塑为17行K列的矩阵，每列对应一个17位整数的比特（LSB到MSB）
    bits_matrix = reshape(bits, 17, K);
    
    % 计算每个整数：每个比特乘以对应的2的幂次后求和
    powers = 2.^(0:16); % 2^0 到 2^16，对应17位的权重
    data = sum(bits_matrix .* powers', 1);
end

使用示例 #

% 打开二进制文件（注意用'rb'模式以二进制读取）
fd = fopen('your_data_file.bin', 'rb');
if fd == -1
    error('无法打开文件');
end

K = 500; % 替换为你需要读取的17位整数个数
data = read_bit17_le(fd, K);

% 关闭文件
fclose(fd);

关键细节说明 #

• 比特顺序处理：用bitget(bytes, 1:8)提取每个字节的LSB到MSB，确保比特流和Matlab的'ieee-le'行为一致。
• 向量化优化：避免逐比特、逐整数的循环，用矩阵操作提升大K值下的运行效率。
• 边界处理：通过ceil(total_bits / 8)计算需要读取的总字节数，确保覆盖所有17位整数的比特，最后截取多余的比特避免干扰。

验证正确性 #

你可以用Matlab生成一个测试文件，再用Octave读取验证结果是否一致：

1. 在Matlab中生成测试数据并写入文件：

   K = 2;
   test_data = uint32([1, 2^16]); % 第一个整数是1，第二个是17位的最大值的一半（2^16）
   fd = fopen('test_bit17.bin', 'wb');
   fwrite(fd, test_data, 'bit17', 'ieee-le');
   fclose(fd);
   

2. 在Octave中读取并验证：

   fd = fopen('test_bit17.bin', 'rb');
   octave_data = read_bit17_le(fd, 2);
   fclose(fd);
   disp(octave_data); % 应输出 [1, 65536]
   

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