6-3coding.rar_分组码

function y=linblkcoderx(x) % here x is recieved vector 6 bits long % channel coding illustration % (6,3) linear block code H=[ ]; % PARITY CHECK MATRIX P=[1 1 0; 0 1 1; 1 0 1]; % P is parity matrix H=[eye(3) P']; Ht=H'; %transpose of H S=zeros(1,3); %syndrome of recieved vector x for i=1:3 S(i)=x(1) & Ht(1,i); for j=2:6 S(i)=bitxor(S(i), (x(j) & Ht(j,i))); end end 'Syndrome of vector is' S %%****SYNDROME LOOK UP TABLE************ %%************************************** if S==[0 0 0] e=[0 0 0 0 0 0]; y=bitxor(x,e); end if S==[1 0 1] e=[0 0 0 0 0 1]; y=bitxor(x,e); end if S==[0 1 1] e=[0 0 0 0 1 0]; y=bitxor(x,e); end if S==[1 1 0] e=[0 0 0 1 0 0]; y=bitxor(x,e); end if S==[0 0 1] e=[0 0 1 0 0 0]; y=bitxor(x,e); end if S==[0 1 0] e=[0 1 0 0 0 0]; y=bitxor(x,e); end if S==[1 0 0] e=[1 0 0 0 0 0]; y=bitxor(x,e); end if S==[1 1 1] e=[0 1 0 0 0 1]; y=bitxor(x,e); end 'error pattern is' e 'correct bits are'