rsc_encoding.m

用Matlab仿真turbocode的算法

function [output1_bin,state1_int]=rsc_encoding(g1,st_t1,out_t1,info)
%RSC coding for each channel
%to get output and the last state of rsc encoder
n=length(info);
mem1=size(g1,2)-1;
state1_int=1;  
output1_int=zeros(1,n+mem1);  
for input=1:n   
    output1_int(input)=out_t1(state1_int,info(input)+1);
    state1_int=st_t1(state1_int,info(input)+1); %next state
end   
info_ter=zeros(1,length(mem1));
for input=1:mem1  %mem1 termination bits
    info_ter(input)=mod(g1(1,2:end)*bin_state(state1_int-1,mem1,'unipolar')',2);   
    output1_int(n+input)=out_t1(state1_int,info_ter(input)+1);
    state1_int=st_t1(state1_int,info_ter(input)+1);
end
buffer=bin_state(output1_int,2,'unipolar')';
output1_bin=buffer(1:end); 
clear buffer;
info=[info,info_ter];
output1_bin=-2*output1_bin+1; %to map unipolar(0/1) to bipolar(+1/-1)