encode1513.m

turbo 码matlab编译码程序

function [Y1,Y2,X]=encode1513(X,I)
% turbo encode (15,13).force zero
L=length(X);
m=3;
T1=0;
T2=0;
T3=0;

for  i=1:L
    Recur_gener=mod(X(i)+T1+T3,2);
    Y1(i)=mod(Recur_gener+T2+T3,2);
    T3=T2;
    T2=T1;
    T1=Recur_gener;
end
for i=L+1:L+m
    if mod(T1+T3,2)==0
        X(i)=0;
    else
        X(i)=1;
    end
    Recur_gener=mod(X(i)+T1+T3,2);
    
    Y1(i)=mod(Recur_gener+T2+T3,2);
    
    T3=T2;
    T2=T1;
    T1=Recur_gener;
end
X_inter=X(I);
for  i=1:length(X_inter)
    Recur_gener=mod(X_inter(i)+T1+T3,2);
    Y2(i)=mod(Recur_gener+T2+T3,2);
    T3=T2;
    T2=T1;
    T1=Recur_gener;
end