rsc_table.m

用Matlab仿真turbocode的算法

function [state_table,output_table,ST]=rsc_table(g)
%to build output%state tables of outer encoder
%g=[1 1 1;1 0 1];


plylen=size(g,2);
plynum=size(g,1);
m = plylen-1;
state = zeros(1,m);
new_state=state;
input=[0 1];
state_table=zeros(2^length(state),2^(plynum-1));
output_table=state_table;
ST=zeros(4,2^(plynum-1)*2^length(state));

for state_skalar=1:2^length(state)					
   
    state=bin_state( state_skalar-1, length(state),'unipolar' );			
        
        
        for x=0:2^(plynum-1)-1
             d_k=input(x+1);
             a_k=mod(g(1,:)*[d_k state]',2);    %feedback
             v_k=mod(g(2,:)*[a_k state]',2);
             output_table(state_skalar,x+1)=int_state([d_k v_k],'unipolar');
             
             new_state=[a_k state(1:m-1)];
             state_table(state_skalar,x+1)=int_state(new_state,'unipolar')+1;
            
        end
           
end

ST(1,:)=mod([0:size(ST,2)-1],2^(plynum-1))+1;
ST_output=output_table';
ST(2,:)=ST_output(1:end)+1;
ST(3,:)=(floor([0:size(ST,2)-1]/(2^(plynum-1)))+1);
states=state_table';
ST(4,:)=states(1:end);