extended_hamming_coding.m

matlab编写的turbo码示例程序

function code = extended_hamming_coding(msg)

%%%% (31,26)
%% x^5+x^2+1
poly=[1 0 0 1 0 1];
L=length(msg);
divisor = [msg,zeros(1,5)];
for k=1:L
    if divisor(k)==1
        divisor(k:k+5)=mod(divisor(k:k+5)+poly,2);
    end
end

code = [msg,divisor(end-4:end)];

one_num=sum(code);

if mod(one_num,2)==1
    code=[code,1];
else
    code=[code,0];
end