log_map_demod.asv

来自「短波信道抗多音干扰的性能分析及其仿真」· ASV 代码 · 共 34 行

function pb=log_map_demod(rho_in_dB,BPH,number_of_states)

N=1000;
fanout=2^BPH;
rho=10^(rho_in_dB/10);
L=floor(log(number_of_states)/log(fanout));
source=[randint(1,N,fanout),zeros(1,L)];
dsource=zeros(1,N*BPH);
if(BPH~=1)
    for i=1:N
        dsource((i-1)*BPH+1:i*BPH)=deci2change(source(i),BPH,2);
    end
else
    dsource=source(1:N);
end

nextstate=zeros(number_of_states,fanout); 
formerstate=number_of_states.*ones(number_of_states,fanout,fanout);
for i=0:number_of_states-1
    for j=0:fanout-1
        next_state=G_func(i,j,L,fanout);
        nextstate(i+1,j+1)=next_state;
        k=1;
        if(formerstate(next_state+1,1,j+1)==number_of_states)
           formerstate(next_state+1,1,j+1)=i; 
       elseif(formerstate(next_state+1,2,j+1)==number_of_states)
           formerstate(next_state+1,2,j+1)=i; 
       elseif(formerstate(next_state+1,3,j+1)==number_of_states)
           formerstate(next_state+1,3,j+1)=i; 
       else
           formerstate(next_state+1,4,j+1)=i; 
    end
end