前言

在上一篇文章《8位16位64位等任意数量用户CDMA直接序列扩频通信系统的Matlab仿真》中，介绍了一种多用户CDMA传输模型，但该模型存在一个缺陷，那就是无论信噪比多高，误码率始终无法为0，为了引入纠错码从而解决这个缺陷，在空闲时间研究了汉明码与交织编码技术，做出了Matlab模型，整理成文章给大家参考。

原理

汉明码

代码中编码与解码均采用奇校验。

编码过程

冗余位数量计算

∀ x ⊂ Z + 且满足 2 x ≥ x + m + 1 ， i = m i n ( x ) \forall x \subset Z^+且满足2^x \geq x+m+1，i = min(x) ∀x⊂Z+且满足2x≥x+m+1，i=min(x)
其中，m为每一帧数据的二进制位数，i为冗余位数量。

校验位位置计算

校验位位置都出现在2的幂次方位，所以根据冗余位数量计算即可。如果冗余位数量为4，则校验位位置就是[1,2,4,8]，注意，以上位置对应数组索引从一开始，即1表示第一位。

计算校验相关位

汉明码就是奇偶校验的升级版，它在编码结果中加入冗余位数量个奇偶校验位。每一个校验位都有不同的校验相关位，第一位校验位，它就跟索引在二进制下第一位是1的所有数据位有关。
假设传输的数据7位一帧，那么它的冗余位数量就是4，第一个冗余位出现在编码后的第一位，这一位要与[3 5 7 9 11]这些位做奇偶校验。
第2位校验位则跟索引在二进制第2位是1的所有数据进行奇偶校验；
。。。
第k位校验位则跟索引在二进制第k位是1的所有数据进行奇偶校验。
代码中用了按位与的方法，把每一个校验位对应的多个相关位找出来。

开始编码

对于每一帧数据，首先把用户的数据填入到编码结果的非校验位，然后对每个校验位及其校验相关位做奇校验，把校验结果填入校验位，即完成一帧的编码。
例如传输的数据7位一帧，那么它的冗余位数量就是4，冗余位序号，索引以及校验相关位如下图所示：

解码过程

解码同样要计算校验位位置与校验相关位，然后对每一个校验位及其校验相关位进行奇校验，校验正确输出0，校验失败输出1。仍是上面的例子，假设接收端接收时出错，上面的编码结果：10111000101的第6位变成了1，即接收到了：10111100101，接收端的解码过程如下：
把接受结果依次与四种冗余相关位进行逻辑与，然后分别统计结果中”1“的个数是否为奇数个，是则输出0，否则输出1。

实验结果

实验参数：

1. 汉明编码一帧对应8个用户码元
2. 交织编码一帧对应16个用户码元
3. 用户码元数设置为1e5 * 8个
4. 信道误码率范围：10^(-4) 到10^(-0.3)
   
   实验结果显示，在信道误码率小于0.128时，使用纠错编码传输对于通信质量有显著改善，甚至在信道误码率小于4e-4时，可以0误码传输，在信道误码率大于0.128时，使用纠错码无实际意义。

仿真代码

可以修改的参数

1. 汉明编码一帧的码元数
2. 交织编码一帧的码元数
3. 测试的码元数量
4. 误码率开始数
5. 误码率结束数
6. 误码率尝试数

下载链接

GitHub: https://github.com/highskyno1/hamming_en-decode
度盘：https://pan.baidu.com/s/13YSblZQoAWjCibgU7Mtmmw
提取码：jrgi

主函数

main.m

%{本代码实现了对用户二进制码元的纠错编码与解码。核心功能采用的汉明码的方法。2020/8/4 23:10加入了交织编码技术,一定程度上避免了突发连续错误带来的误码2020/8/5 11:19
%}
%汉明编码一帧的码元数,对数据冗余度有影响
hamming_frame_size = 8;
%交织编码一帧的码元数,越大对于抵抗突发连续错误越好,但会增加处理时延
mixed_frame_size = 16;
%测试的码元数量,必须是frame_size以及mixed_frame_size的整数倍
code_num = 1e5 * 8;%生成用户码元
user_code = genBipolar(code_num);
%交织编码
mix_code = mixed_encode(user_code,mixed_frame_size);
%发送出去的编码后的码元
send_code = hamming_encode(mix_code,hamming_frame_size);
fprintf('冗余度：%f\n',length(send_code)/length(user_code));
%误码率开始数
error_start = -4; % -4指10^-4
%误码率结束数
error_end = -0.3; % -0.3指10^(-0.3)
%误码率尝试数
error_size = 100;
%记录信道误码率的数组
channel_error_rec = zeros(1,error_size);
%记录实际误码率的数组
error_rec = zeros(1,error_size);
parfor error_index = 1:error_size%当前信道误码率error_rate = 10^(error_start + (error_index-1) / error_size * (error_end - error_start));channel_error_rec(error_index) = error_rate;%随机误码rec_code = randErrorCode(send_code,error_rate);%解码后的码元de_code = hamming_decode(rec_code,hamming_frame_size);%交织解码res_code = mixed_decode(de_code,mixed_frame_size);%计算实际误码率[~,accuracy] = compare(user_code,res_code);error_rec(error_index) = 1 - accuracy;fprintf('信道误码率：%f,实际误码率：%f\n',error_rate,1-accuracy);
end
%绘制实际误码率与信道误码率关系图
figure();
loglog(channel_error_rec,error_rec);
hold on;
loglog(channel_error_rec,channel_error_rec);
xlabel('信道误码率');
ylabel('实际误码率');
legend('使用纠错码传输','直接传输');

汉明码编解码测试模块

hamming_test.m

user_code = [1 0 0 1 1 0 1];
en_code = encode(user_code,7);
%误码位
en_code(11) = ~en_code(11);
de_code = decode(en_code,7);
%计算正确率
[~,accuracy] = compare(user_code,de_code);
fprintf('误码率：%f\n',1-accuracy);

汉明码编码器

hamming_encode.m

function res = hamming_encode(input,frame_size)
%encode 对输入的二进制数组进行汉明码编码,奇校验
%input:需要编码的二进制数组
%frame_size:每一帧数据的二进制位数
len_input = length(input);
if mod(len_input,frame_size) ~= 0error('输入的数组数必须是分组数的整数倍')
end
%计算冗余位数量
red_bit_num = red_bit_cal(frame_size);
%总的位数
bit_num = red_bit_num + frame_size;
%帧数量
frame_num = len_input / frame_size;
%计算校验位的位置
pos_red_bit = 2.^(0:red_bit_num-1);
%计算校验位的相关位
pos_relat_bit = uint8(zeros(red_bit_num,bit_num));
for i = 1:red_bit_numfor j = 1:bit_numif bitand(2^(i-1),j) ~= 0pos_relat_bit(i,j) = 1;endend
end
%初始化一个结果单元
res = uint8(zeros(frame_num,bit_num));
%填满每一帧的数据位
index = 1;
for i = 1:frame_numfor j = 1:bit_numif ismember(j,pos_red_bit)continueendres(i,j) = input(index);index = index + 1;end
end
%对每一帧加入校验位
for i = 1:frame_numfor j = 1:red_bit_numfoo = bitand(res(i,:),pos_relat_bit(j,:));foo = mod(sum(foo)+1,2);res(i,pos_red_bit(j)) = foo;end
end
res = res';
res = res(:);
end

汉明码解码器

hamming_decode.m

function res = hamming_decode(input,frame_size)
%decode 对输入的二进制数组进行汉明码解码,奇校验
%input:需要编码的二进制数组
%frame_size:每一帧数据的二进制位数
len_input = length(input);
%计算冗余位数量
red_bit_num = red_bit_cal(frame_size);
%总的位数
bit_num = red_bit_num + frame_size;
if mod(len_input,bit_num) ~= 0error('输入的数组数必须是总的位数的整数倍');
end
%帧数量
frame_num = len_input / bit_num;
%对输入进行重构
input = uint8(input);
input = reshape(input,bit_num,frame_num);
input = input';
%计算校验位的位置
pos_red_bit = 2.^(0:red_bit_num-1);
%计算校验位的相关位
pos_relat_bit = uint8(zeros(red_bit_num,bit_num));
for i = 1:red_bit_numfor j = 1:bit_numif bitand(2^(i-1),j) ~= 0pos_relat_bit(i,j) = 1;endend
end
%开始校验,注意是奇校验
for i = 1:frame_num%该帧的校验结果res_verify = uint8(zeros(1,red_bit_num));for j = 1:red_bit_numfoo = sum(bitand(input(i,:),pos_relat_bit(j,:)));if mod(foo,2) == 0%奇校验失败res_verify(j) = 1;endend%计算出错位%转换为小端序res_verify = fliplr(res_verify);error_bit_pos = 0;for j = 1:red_bit_numerror_bit_pos = error_bit_pos*2 + res_verify(j);end%出错位取反if error_bit_pos <= bit_num && error_bit_pos ~=0%如果为0，则没有发生传输错误input(i,error_bit_pos) = ~input(i,error_bit_pos);end
end
%提取结果
res = uint8(zeros(1,frame_num*frame_size));
res_index = 1;
for i = 1:frame_numfor j = 1:bit_numif ismember(j,pos_red_bit)continueendres(res_index) = input(i,j);res_index = res_index + 1;end
end
end

冗余位计算模块

red_bit_cal.m

function res = red_bit_cal(group_num)
%red_bit_cal 计算冗余位数量
%group_num:每一帧数据的二进制位数
res = 0;
while 2^res - res < group_num + 1res = res + 1;
end
end

交织编码器

mixed_encode.m

function res = mixed_encode(input,frame_size)
%mined_encode 本函数用于对输入信号进行交织编码
%input：用户数据
%frame_size：一帧数据所包含的二进制码元数
len_input = length(input);
if mod(len_input,frame_size) ~= 0error('用户数据二进制数必须为帧数的整数倍');
end
frame_num = len_input / frame_size;
res = reshape(input,frame_num,frame_size);
res = res';
res = res(:);
end

交织解码器

mixed_decode.m

function res = mixed_decode(input,frame_size)
%mined_encode 本函数用于对输入信号进行交织解码
%input：用户数据
%frame_size：一帧数据所包含的二进制码元数
len_input = length(input);
if mod(len_input,frame_size) ~= 0error('用户数据二进制数必须为帧数的整数倍');
end
frame_num = len_input / frame_size;
res = reshape(input,frame_size,frame_num);
res = res';
res = res(:);
end

随机误码模块

randErrorCode.m

function res = randErrorCode(input,error_rate)
%randErrorCode 对输入的码元数组添加随机误码
%input:数组
%error_rate:误码率(%)
res = input;
for i = 1:length(input)if rand() < error_rateres(i) = ~res(i);end
end
end

比较模块

compare.m

%比较两个数组
%input1:数组1
%input2:数组2
%res:正确码元数量
%accuracy:正确率
function [res,accuracy] = compare(input1,input2)len1 = length(input1);len2 = length(input2);len = min(len1,len2);right_code = 0;for i = 1:lenif input1(i) == input2(i)right_code = right_code + 1;endendres = right_code;accuracy = right_code / len;
end

单极性码生成模块

genBipolar.m

%产生单极性码
%num:单极性码的规模
%res:满载单极性码的数组
function res = genBipolar(num)res = rand(1,num);res = value2Bipolar(res);
end

随机码转单极性码模块

value2Bipolar.m

% 本函数用户把随机生成的0~1之间的double数转换成双极性码
function res = value2Bipolar(input)res = zeros(1,length(input));for index = 1:length(input)if(input(index) > 0.5)res(index) = 1;endendres = uint8(res);
end

后语

终于把纠错码模型做出来了，接下来可以去修改任意用户CDMA传输模型了。前几天趁着有空，去了张家界旅游，自由行，一共5天。那边空气非常清新，但是景区的设计非常坑爹，经常都是至少坐一次缆车或者天梯，才能在景区下班前离开景区，但风景是真的不错。昨天收到通知，初步确定9月24号才能回到学校，这个暑假又延长了！

【通信系统仿真系列】基于Matlab的汉明码(Hamming Code)纠错传输以及交织编码(Interleaved coding)仿真相关推荐

1. 【OFDM通信】基于matlab OFDM通信系统仿真【含Matlab源码 2521期】

   ⛄一.获取代码方式 (附课程作业报告) 获取代码方式1: 完整代码已上传我的资源:[OFDM通信]基于matlab OFDM通信系统仿真[含Matlab源码 2521期] 点击上面蓝色字体,直接付费下 ...

2. 【数字信号】基于matlab GUI多音双频（DTMF）拨号音频解码仿真系统【含Matlab源码 1084期】

   ⛄一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[数字信号]基于matlab GUI多音双频(DTMF)拨号音频解码仿真系统[含Matlab源码 1084期] 点击上面蓝色字体,直接付费 ...

3. 基于Matlab模拟、检测和跟踪飞机着陆进场中异常的仿真（附源码）

   目录 一.介绍 二.生成和标记轨迹 三.定义方案 四.运行方案并检测异常轨道 五.将跟踪异常报告与事实进行比较 六.总结 七.程序 该示例显示了如何自动检测最终接近机场跑道的飞机的偏差和异常.在此示例 ...

4. 【通信系统仿真系列】基于延迟相乘鉴相的2QPSK调相通信系统仿真

   基于延迟相乘鉴相的2QPSK调相通信系统仿真 1 前言 2 原理 2.1 码元生成 2.2 差分处理 2.2.1 不使用差分码 2.2.2 使用差分码 2.3 载波 2.4 调制 2.5 经过高斯信道 ...

5. 基于matlab的am系统仿真论文,基于Matlab的AM通信系统仿真

   基于Matlab的AM通信系统仿真 分析当前通信类课程实验教学存在的问题,提出采用Matlab仿真来弥补实验室实验设备等的不足.利用Mat-lab的Simulink工具箱建立AM系统的两种仿真模型,包 ...

6. matlab控制系统仿真论文,基于MATLAB的过程控制系统仿真毕业设计论文.doc

   基于MATLAB的过程控制系统仿真毕业设计论文 毕业设计论文 基于MATLAB的过程控制系统仿真研究 摘 要 在水箱液位控制系统中,通过建立数学模型以及实验中对实验数据的分析,分别确定了单容.双容.三 ...

7. 电力系统仿真：基于MATLAB的电力系统稳态和暂态仿真，分析电力系统的稳定性、损耗和负载特性

   电力系统仿真报告: 一.问题定义 本次仿真案例拟对一个电力系统进行暂态和稳态分析,以验证电力系统的稳定性和其它特性.具体目标包括: 建立电力系统的暂态和稳态模型 分析电力系统在不同负荷和拓扑下的稳定性 ...

8. 基于matlab的风力发电系统仿真研究,基于MATLAB的风力发电系统仿真研究

   基于MATLAB的风力发电系统仿真研究 向恺X刘永前 (华北电力大学动力工程系,北京102206) 摘要:本文建立了风力发电系统风轮的数学模型.传动系统模型.三相异步发电机的数学模型,并用MATLAB ...

9. 基于matlab 系统仿真学报,基于MATLAB/Simulink的混合动力汽车正向仿真软件的开发...

   系 统 仿 真 学 报 Vol. 17 增刊 JOURNAL OF SYSTEM SIMULATION Dec. 2005 16 基于 MATLAB/Simulink 的混合动力汽车正向仿真软件的开发 ...

最新文章

1. Reactor模式与Proactor模式
2. 二分图带权最大匹配费用流_一文掌握阻抗匹配
3. Python 随机森林特征重要度
4. 七夕了！！聊聊《最受欢迎的男友职业排行榜Top10》
5. Flex与javascript通信
6. Tomcat 内存调大
7. setseed_Java Random setSeed（）方法与示例
8. Oracle11g的安装和使用
9. oracle数据库直方图,数据库 – Oracle伪造直方图
10. vue项目运用繁体字
11. 屏幕录制专家——录制视频没声音的解决办法
12. 1688淘口令链接API接口-item_password-获得淘口令真实url 接口，淘口令API接口
13. k8s ceph rbd动态挂载
14. 用vmware安装gho文件
15. 第114课:SparkStreaming+Kafka+Spark SQL+TopN+Mysql+KafkaOffsetMonitor电商广告点击综合案例实战(详细内幕版本)
16. 妈妈计算机英语怎么说,妈妈的英文翻译,妈妈英语怎么说
17. 境外服务器网站能挂支付宝么,支付宝如何境外改为境内？在国外可以使用吗？...
18. 我国最高山峰是珠穆朗玛峰，8848米。现在我有一张足够大的纸，它的厚度是0.01米。请问，我折叠多少次，可以折成珠穆朗玛峰的高度?
19. 知识汇总：Python办公自动化应该学习哪些内容
20. C语言练习题——函数

热门文章

1. 信号完整性之浅谈理解（一）
2. Android View截图
3. 模拟角频率和数字角频率
4. 【UML】——活动图
5. 介绍一款web端的截图工具
6. [论文笔记] Google Spanner Distributed Database
7. 彻底解决：fake_useragent.errors.FakeUserAgentError: Maximum amount of retries reached
8. 强化学习：Actor-Critic、SPG、DDPG、MADDPG
9. 如何将CSDN博客下载为PDF文件
10. DOSBOX使用经验