补充一点,关于生成矩阵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);