t_tpc2d.m

B3g_phase2_C语言_Matlab程序及说明

%test for 2D tpc code
clear;
%define the code structure
M=[6,6];%3 -> (8,4,4); 4-> (16, 11, 4); 5->(32,26,4);6->(64,57,4)
N=2.^M;
K=N-M-1;
msg_L=K(1)*K(2);
code_L=N(1)*N(2);

iter_num=4;
delta=4;

Rate=msg_L/code_L

snr1=3.5; snr2=3.5;
for snr=snr1:0.2:snr2
    snr=snr+10*log10(Rate)+3;
    errors=0;
    for kk=1:1000
        kk
        msg=randint(1,msg_L);
        code=tpc_encode2d(msg, M);
        s=2*code-1;
        r=awgn(s,snr);
        [decod llr_all]=tpc_decode2d(r,M,iter_num, delta);
        err=sum(abs(decod-msg));
        errors=errors+err;
        ber=errors/kk/msg_L;
        errors
        ber
    end
end