【代码】LDPC校验矩阵转换到alist格式
warmonkey2018/11/27电子技术 IP:广东

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

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

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

alist是一种通用的稀疏矩阵存储格式,具体定义见 XXXXXXXXXXXXXXXXXXXX/mackay/codes/XXXXXXXXml

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

<code>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</code>

调用示例:

<code>% convert sparse matrix to alist format
% <a href=/l?t=aHR0cDovL3d3dy5pbmZlcmVuY2Uub3JnLnVrL21hY2theS9jb2Rlcy9hbGlzdC5odG1s target="_blank">http://www.inference.org.uk/mackay/codes/alist.html</a>

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');</code>

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


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

来自:电子信息 / 电子技术
1
4
已屏蔽 原因:{{ notice.reason }}已屏蔽
{{notice.noticeContent}}
~~空空如也

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

所属专业
所属分类
上级专业
同级专业
warmonkey
学者 机友
文章
357
回复
7605
学术分
14
2008/10/11注册,16分26秒前活动

Cubesat

主体类型:个人
所属领域:无
认证方式:手机号
IP归属地:未同步
文件下载
加载中...
{{errorInfo}}
{{downloadWarning}}
你在 {{downloadTime}} 下载过当前文件。
文件名称:{{resource.defaultFile.name}}
下载次数:{{resource.hits}}
上传用户:{{uploader.username}}
所需积分:{{costScores}},{{holdScores}}下载当前附件免费{{description}}
积分不足,去充值
文件已丢失

当前账号的附件下载数量限制如下:
时段 个数
{{f.startingTime}}点 - {{f.endTime}}点 {{f.fileCount}}
视频暂不能访问,请登录试试
仅供内部学术交流或培训使用,请先保存到本地。本内容不代表科创观点,未经原作者同意,请勿转载。
音频暂不能访问,请登录试试
支持的图片格式:jpg, jpeg, png
插入公式
评论控制
加载中...
文号:{{pid}}
投诉或举报
加载中...
{{tip}}
请选择违规类型:
{{reason.type}}

空空如也

加载中...
详情
详情
推送到专栏从专栏移除
设为匿名取消匿名
查看作者
回复
只看作者
加入收藏取消收藏
收藏
取消收藏
折叠回复
置顶取消置顶
评学术分
鼓励
设为精选取消精选
管理提醒
编辑
通过审核
评论控制
退修或删除
历史版本
违规记录
投诉或举报
加入黑名单移除黑名单
查看IP
{{format('YYYY/MM/DD HH:mm:ss', toc)}}