【代码】LDPC校验矩阵转换到alist格式

title: converting ldpc parity matrix to alist format (for gnuradio ldpc encoder use)

LDPC(低密度奇偶校验码)是通信中常用的FEC(前向纠错)编码之一。

GnuRadio的LDPC编解码器接受的输入格式是alist格式的校验矩阵(通常记作H矩阵),本文提供了将H矩阵转换为alist格式的工具。

alist是一种通用的稀疏矩阵存储格式,具体定义见 www.inference.org.uk/mackay/codes/alist.html

可用以下matlab代码将矩阵x转换为alist,并保存到名字为filename的文件。

function write_alist(x, filename)
    n = size(x,1);
    m = size(x,2);
    num_nlist = sum(x,1)';
    num_mlist = sum(x,2);
    biggest_num_m = max(num_mlist);
    biggest_num_n = max(num_nlist);
    for i = 1:n
        nlist{i} = find(x(i, :) > 0);
    end
    for i = 1:m
        mlist{i} = find(x(:, i) > 0)';
    end

    fd = fopen(filename, 'w');
    
    fprintf(fd, '%d %d\n', n, m);
    fprintf(fd, '%d %d\n', biggest_num_n, biggest_num_m);
    
    fprintf(fd, '%d ', num_nlist);
    fprintf(fd, '\n');
    fprintf(fd, '%d ', num_mlist);
    fprintf(fd, '\n');
    
    for i=1:length(nlist)
        fprintf(fd, '%d ', nlist{i});
        fprintf(fd, '\n');
    end
    for j=1:length(mlist)
        fprintf(fd, '%d ', mlist{j});
        fprintf(fd, '\n');
    end
    
    fclose(fd);
end

调用示例:

% convert sparse matrix to alist format
% http://www.inference.org.uk/mackay/codes/alist.html

clear all;
x= [0 0 1 0 0 0 0 1 0 1 0 0 1 0 0 0;
0 0 0 1 0 0 1 0 1 0 0 0 1 0 0 0;
0 1 0 0 1 0 1 0 0 1 0 0 0 0 0 0;
0 0 0 1 0 1 0 0 0 0 1 0 0 1 0 0;
0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 1;
1 0 0 0 0 1 0 0 1 1 0 0 0 0 0 0;
0 0 0 1 0 0 0 1 0 0 0 1 0 0 1 0;
0 1 0 0 0 1 0 0 0 0 0 1 0 0 0 1;
1 0 0 0 0 0 1 0 0 0 0 0 0 1 0 1;
0 0 1 0 1 0 0 0 0 0 0 1 0 1 0 0;
0 1 0 0 0 0 0 0 0 0 1 0 1 0 1 0;
1 0 0 0 1 0 0 1 0 0 1 0 0 0 0 0;];

write_alist(x, 'x.alist');

输出结果与原文链接中的例子完全一致。


======TODO:补充CCSDS C2 LDPC 7/8的例子=====

来自:电子信息 / 电子技术
4

想参与大家的讨论?现在就 登录 或者 注册

warmonkey
专家 学者 机友 笔友
文章
317
回复
6339
学术分
17
2008/10/11注册,2 天前活动

科创航天主席

插入资源
全部
图片
视频
音频
附件
全部
未使用
已使用
正在上传
空空如也~
上传中..{{f.progress}}%
处理中..
上传失败,点击重试
{{f.name}}
空空如也~
(视频){{r.oname}}
{{selectedResourcesId.indexOf(r.rid) + 1}}
插入表情
我的表情
共享表情
Emoji
上传
注意事项
最大尺寸100px,超过会被压缩。为保证效果,建议上传前自行处理。
建议上传自己DIY的表情,严禁上传侵权内容。
点击重试等待上传{{s.progress}}%处理中...已上传
空空如也~
草稿箱
加载中...
此处只插入正文,如果要使用草稿中的其余内容,请点击继续创作。
{{fromNow(d.toc)}}
{{getDraftInfo(d)}}
标题:{{d.t}}
内容:{{d.c}}
继续创作
删除插入插入
{{forum.displayName}}
{{forum.countThreads}}
篇文章,
{{forum.countPosts}}
条回复
{{forum.description || "暂无简介"}}
ID: {{user.uid}}
{{submitted?"":"投诉"}}
请选择违规类型:
{{reason.description}}
支持的图片格式:jpg, jpeg, png