algdecode.m

TPC的Chase迭代译码时测试图样的程序

%------------------------------硬判决译码---------------------
%------------------从伴随式s中找出错误图样e--------------------
function e=algdecode(H,r)
% r=[0     0     1     0     1     1     0     0     1     1     1     1     1     0 1];
% H=[1     0     0     1     1     0     1     0     1     1     1     1     0     0 0];
[M,N]=size(H);
n_k=2^M;
s(n_k,:)=de2bi(n_k-1);
len=length(s(n_k,:));
for i=0:n_k-2
    s(i+1,:)=de2bi(i,len);
end
for i=1:n_k
    ee(i,:)=s(i,:)*(pinv(H))';
end
s_r=r*H';
len_e=length(ee(1,:));
e=zeros(1,len_e);
for i=1:n_k
    if s_r==s(i,:)
        e=ee(i,:);
    end
end