7-4线性分组码的matlab程序

补充一点，关于生成矩阵G和校验矩阵H的关系：生成矩阵G是4*7矩阵，分为两块，前4列组成的4*4矩阵为单位矩阵，后3列组成的4*3矩阵我称它监督矩阵，校验矩阵是3*7矩阵，前4列组成的3*4矩阵是监督矩阵的转置得到的，后面三列组成的3*3矩阵是单位矩阵。

%设本组（7,4）线性分组码信息位和检验位的约束关系为：

%c5=c1+c2+c3;c6=c2+c3+c4;c7=c1+c2+c4；

%由约束关系可以写出生成矩阵G.

clear all;

G1=[1,0,0,0;0,1,0,0;0,0,1,0;0,0,0,1];

G2=[1,0,1;1,1,1;1,1,0;0,1,1];

G=[G1,G2];

C=[1,0,1,1];  %C=[c1,c2,c3,c4]要进行编码信息码字 

R1=C*G;

R=mod(R1,2);

fprintf('输出的编码为：R=')

disp(R);

 

 

 

%已知生成矩阵G，可以求出校验矩阵H。

clear all;

E=[1,1,1,1,1,1,1];

H1=[1,1,1,0;0,1,1,1;1,1,0,1];

H2=[1,0,0;0,1,0;0,0,1];

H=[H1,H2];%校验矩阵。

R=[1,0,1,1,0,0,1];%接收到的码字。

S1=R*(H'); %S为校阵子；

S=mod(S1,2);

 

 

for i=1:7; %用改for循环取出H中每一列，然后与S相加。

 T=H(:,[i]);

 %disp(T);

 B1=S+T';

 B=mod(B1,2);

 

    if (all(B(:)==0)); %如果S与H的第i列之和B为0矩阵，则表示R中第i个码字有误。

      fprintf('R中序列中错误码位是第：');

      disp(i) 

    else

     E(1,i)=0;

     

    end;

   

 

end;

C=mod((R+E),2);

fprintf('纠错后的码字应该为：C=');

disp(C);