MATLAB实现的Reed-Muller（RM码，里德-马勒编码）编码解码纠错以及BER分析

最新推荐文章于 2022-12-03 19:35:26 发布

H4de5

最新推荐文章于 2022-12-03 19:35:26 发布

阅读量2.5k

点赞数 2

分类专栏：编码文章标签： matlab 开发语言编码学

本文链接：https://blog.csdn.net/weixin_44480014/article/details/123205380

版权

MATLAB实现的Reed-Muller编码解码纠错以及BER分析

背景
代码
- 计算BER
计算不使用RM编码情况时的BER（模拟环境与理论情况）
与汉明编码做对比
总结

背景

本科时信息论与编码的作业，RM（2，4）编码

课程为

Information Theory & Coding
Vaibhav Kumar, PhD
School of Electrical & Electronic Engineering
University College Dublin – The Republic of Ireland

下述作业也是老师布置的，如果涉及到版权问题我会删掉该博客。流程逻辑如下：

数据->RM编码->BPSK调制->加入AWGN（模拟传输时的噪声）->BPSK解调->解码（Majority-Logic解码）->输出数据

给定生成矩阵为

v0                 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
v4                 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1
v3                 0 0 0 0 1 1 1 1 0 0 0 0 1 1 1 1
v2                 0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1
v1                 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1
v3v4      =        0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1
v2v4               0 0 0 0 0 0 0 0 0 0 1 1 0 0 1 1
v1v4               0 0 0 0 0 0 0 0 0 1 0 1 0 1 0 1
v2v3               0 0 0 0 0 0 1 1 0 0 0 0 0 0 1 1
v1v3               0 0 0 0 0 1 0 1 0 0 0 0 0 1 0 1
v1v2               0 0 0 1 0 0 0 1 0 0 0 1 0 0 0 1

接收到的数据可以写成r=v+e的格式

问题描述上
问题描述下

代码

clc,clear
Eb_N0_log=0:0.5:10; %dB
Eb_N0=10.^(Eb_N0_log/10);
R_coded=11/16; %R
BER_coded=zeros(1,21); %store BER in different Eb/N0
G=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1; %v0
   0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1; %v4
   0 0 0 0 1 1 1 1 0 0 0 0 1 1 1 1; %v3
   0 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1; %v2
   0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1; %v1
   0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1; %v3v4
   0 0 0 0 0 0 0 0 0 0 1 1 0 0 1 1; %v2v4
   0 0 0 0 0 0 0 0 0 1 0 1 0 1 0 1; %v1v4
   0 0 0 0 0 0 1 1 0 0 0 0 0 0 1 1; %v2v3
   0 0 0 0 0 1 0 1 0 0 0 0 0 1 0 1; %v1v3
   0 0 0 1 0 0 0 1 0 0 0 1 0 0 0 1; %v1v2
   ]; %generator matrix
for i=1:21 %calculate BER in different Eb/N0
    SNR_coded=2*R_coded*Eb_N0(i); %calculate SNR
    message=round(rand(1,1980000)); %generate original data u = (a0 a4 a3 a2 a1 a34 a24 a14 a23 a13 a12),
    px_w=1; %signal's power
    pn_w=px_w/SNR_coded; %noise's power
    error_number=0; %error