基于matlab的线性分组码的编码译码程序
一、实验目的
1.熟悉线性分组码的编码译码原理
2.掌握编写线性分组码的编码译码程序的要点
3.掌握使用matlab的仿真要点
二、实验内容
1.完成对任意信息序列的编码
2.根据生成矩阵，形成监督矩阵
3.根据得到的监督矩阵，得到伴随式，并根据它进行译码
4.验证工作的正确性
三、实验代码
G=[1 0 0 1 0 1; 0 1 0 1 1 0; 0 0 1 0 1 1] H=gen2par(G) disp('监督矩阵为:H='); disp(H); I=[0 0 0;0 0 1;0 1 0;0 1 1;1 0 0;1 0 1;1 1 0;1 1 1]; C=rem(I*G,2); disp('所得许用码组结果为:C='); disp(C) %译码并判别 clear all; close all; G=[1 0 0 1 0 1; 0 1 0 1 1 0; 0 0 1 0 1 1]; R=input('请输入接收到的码组R:'); [a,b]=size(R) E=[0 0 0 0 0 0;1 0 0 0 0 0; 0 1 0 0 0 0;0 0 1 0 0 0; 0 0 0 1 0 0;0 0 0 0 1 0;0 0 0 0 0 1]; H=gen2par(G); S=rem(R*H',2); disp('所得伴随式为:S='); disp(S); i=1; for i=1:1:a M(i,1)=S(i,1).*4+S(i,2).*2+S(i,3); end for i=1:1:a switch(M(i,1)) case 0 A(i,:)=R(i,:)+E(1,:); case 5 A(i,:)=R(i,:)+E(2,:); case 6 A(i,:)=R(i,:)+E(3,:); case 3 A(i,:)=R(i,:)+E(4,:); case 4 A(i,:)=R(i,:)+E(5,:); case 2 A(i,:)=R(i,:)+E(6,:); case 1 A(i,:)=R(i,:)+E(7,:); end end for i=1:1:a switch(M(i,1)) case 0 disp('没有出现错误!'); case 1 disp('注意:第1位出现一个错误!请纠正!'); case 2 disp('注意:第2位出现一个错误!请纠正!'); case 4 disp('注意:第3位出现一个错误!请纠正!'); case 3 disp('注意:第4位出现一个错误!请纠正!'); case 6 disp('注意:第5位出现一个错误!请纠正!'); case 5 disp('注意:第6位出现一个错误!请纠正!'); end end A=rem(A,2); disp('检纠错后的码组A='); disp(A); j=1; while j<=3 I(:,j)=A(:,j); j=j+1; end disp('译出的信息序列I='); disp(I);
