hammingdecode.m

AWGN信道下高斯噪声干扰下GPS编码译码系统性能仿真

function t=HammingDecode(T)
%函数用于纠错
%T为输入一个500行的15列的（15，11）编码后的矩阵
%t为输出一个1行7500列的（15，11）编码纠错后的矩阵
%函数只能纠正一位错码,当有一位以上错码时视为超出纠错范围
T=reshape(T,1,32);
H=[1 0 1 0 0 1;
   0 1 0 1 1 0;
   1 0 1 0 1 0;
   1 1 0 1 0 0;
   1 1 1 0 1 1;
   0 1 1 1 0 0;
   1 0 1 1 1 1;
   1 1 0 1 1 1;
   0 1 1 0 1 0;
   0 0 1 1 0 1;
   0 0 0 1 1 1;
   1 0 0 0 1 1;
   1 1 0 0 0 1;
   1 1 1 0 0 0;
   1 1 1 1 0 1;
   1 1 1 1 1 0;
   0 1 1 1 1 1;
   0 0 1 1 1 0;
   1 0 0 1 1 0;
   1 1 0 0 1 0;
   0 1 1 0 0 1;
   1 0 1 1 0 0;
   0 1 0 1 1 0;
   0 0 1 0 1 1;
   1 0 0 1 0 1;
   0 1 0 0 1 1;
   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];
t=[];
R=T;
S=T*H;
S=mod(S,ones(1,6)*2);
if S==[0 0 0 0 0 1]
  R=xor(R,[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1]);
elseif S==[0 0 0 0 0 0]
      R=xor(R,[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]);
elseif S==[0 0 0 0 1 0]
      R=xor(R,[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0]);
elseif S==[0 0 0 1 0 0]
      R=xor(R,[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0]);
elseif S==[0 0 1 0 0 0]
      R=xor(R,[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0]);
elseif S==[0 1 0 0 0 0]
      R=xor(R,[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0]);
elseif S==[1 0 0 0 0 0]
      R=xor(R,[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0]);
elseif S==[0 1 0 0 1 1]
      R=xor(R,[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0]);
elseif S==[1 0 0 1 0 1]
      R=xor(R,[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0]);
elseif S==[0 0 1 0 1 1]
      R=xor(R,[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0]);
elseif S==[0 1 0 1 1 0]
      R=xor(R,[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0]);
elseif S==[1 0 1 1 0 0]
      R=xor(R,[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0]);
elseif S==[0 1 1 0 0 1]
      R=xor(R,[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0]);
elseif S==[1 1 0 0 1 0]
      R=xor(R,[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0]);
elseif S==[1 0 0 1 1 0]
      R=xor(R,[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0]);
elseif S==[0 0 1 1 1 0]
      R=xor(R,[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0]);
elseif S==[0 1 1 1 1 1]
      R=xor(R,[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]);
elseif S==[1 1 1 1 1 0]
      R=xor(R,[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]);
elseif S==[1 1 1 1 0 1]
      R=xor(R,[0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]);
elseif S==[1 1 1 0 0 0]
      R=xor(R,[0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]);
elseif S==[1 1 0 0 0 1]
      R=xor(R,[0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]);
elseif S==[1 0 0 0 1 1]
      R=xor(R,[0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]);
elseif S==[0 0 0 1 1 1]
      R=xor(R,[0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]);
elseif S==[0 0 1 1 0 1]
      R=xor(R,[0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]);
elseif S==[0 1 1 0 1 0]
      R=xor(R,[0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]);
elseif S==[1 1 0 1 1 1]
      R=xor(R,[0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]);
elseif S==[1 0 1 1 1 1]
      R=xor(R,[0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]);
elseif S==[0 1 1 1 0 0]
      R=xor(R,[0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]);
elseif S==[1 1 1 0 1 1]
      R=xor(R,[0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]);
elseif S==[1 1 0 1 0 0]
      R=xor(R,[0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]);
elseif S==[1 0 1 0 1 0]
      R=xor(R,[0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]);
elseif S==[0 1 0 1 0 1]
      R=xor(R,[0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]);
elseif S==[1 0 1 0 0 1]
      R=xor(R,[1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]);
end; 
if R(2)==0
        R=R;
end;    
if R(2)==1
        R=xor(R,[0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 1 0 1 1 0]);
end;  
t=R(3:26);