ber512.m

turbo 均衡的经典文献matlab实现

%TE 主程序
%corpyright:flyingxiang

clear all;
close all;
eb=2;
ebno=10.^([0:5]/10);
No=eb./ebno;
var=No./2;


limiar_err=1000;%传输中最小错误总数






for i=1:length(ebno)
    err0=0;err1=0;err2=0;err3=0;err4=0;
    k=0;j=0;l=0;g=0;p=0;
    while  err4<limiar_err
        %______________________________________________________________________
        %初始状态,迭代为零
        
        [bits_trans,simb_cod,Ldc0,Ldb0, Lec_deint0,y,numero0,taxa0]=it0_512(var(i));
        err0=err0+numero0;
        j=j+1;num_bit0=j;             
        %______________________________________________________________________
        %迭代一次
        [Ldc1,Ldb1,Lec_deint_apriori_dec1,numero1,taxa1]=itN_512(var(i),Ldc0,Lec_deint0,y,bits_trans);
        err1=err1+numero1; 
        k=k+1;num_bit1=k;
        
        %______________________________________________________________________
        %迭代两次
        [Ldc2,Ldb2,Lec_deint_apriori_dec2,numero2,taxa2]=itN_512(var(i),Ldc1,Lec_deint_apriori_dec1,y,bits_trans);
        err2=err2+numero2 ;
        l=l+1;num_bit2=l;
        
        %______________________________________________________________________
        %迭代三次
        [Ldc3,Ldb3,Lec_deint_apriori_dec3,numero3,taxa3]=itN_512(var(i),Ldc2,Lec_deint_apriori_dec2,y,bits_trans);
        err3=err3+numero3; 
        g=g+1;num_bit3=g;
        %______________________________________________________________________
        %迭代四次
        [Ldc4,Ldb4,Lec_deint_apriori_dec4,numero4,taxa4]=itN_512(var(i),Ldc3,Lec_deint_apriori_dec3,y,bits_trans);
        err4=err4+numero4; 
        p=p+1;num_bit4=p;
        end
        eblog(i)=10*log10(ebno(i));
        %计算BER
        ber0(i)=err0/(num_bit0*512); 
        ber1(i)=err1/(num_bit1*512);  
        ber2(i)=err2/(num_bit2*512);
        ber3(i)=err3/(num_bit3*512);
        ber4(i)=err4/(num_bit4*512);
        
        
        %生成BER--SNR图像
        semilogy(eblog,ber0,'rx-',eblog,ber1,'bx-',eblog,ber2,'kx-',eblog,ber3,'go-',eblog,ber4,'mo-');grid;
        %semilogy(eblog,ber0,'rx-');grid;
        legend('iterative0','iterative1','iterative2','iterative3','iterative4');
        xlabel('Eb/No','FontSize',16);
        ylabel('BER','FontSize',16);
        pause(.0001);
    end