trellis_out.m

LDPC码的编解码Matlab仿真 希望对大家有用

function [state,state_out]=trellis_out(coef_expand,N_sample)

fltr='filt_0.125b_7t.txt'; 
fltr_coef=textread(fltr);

for i=1:2^(coef_expand)
    ind=dec2bin((i-1),coef_expand);
    for j=1:length(ind)
        if ind(j)=='0'
            state(i,j)=0;
        else
            state(i,j)=1;
        end
    end
    
    wav(i,:)=(rcos_pulse_mod(state(i,:),N_sample))';
    filt_out(i,:)=filter(fltr_coef,1,wav(i,:));
    state_out(i)=filt_out(i,length(filt_out(i,:)));
    
end