bertest.m.bak

采用 DVB反向链路（DVB-RCS）标准中的双二元Turbo码matlab仿真程序和仿真结果。支持1/2和1/3码率。支持MAP

%双二元Turbo码的顶层误码测试程序
clear
clc
snr_db=2:0.25:2.5;
% snr_db=0.5;
% snr_db=[2.25 2.5];
snr=10.^(snr_db/10);
% N=212;   %帧长,对应于ATM包
 N=440;
% N=228;
%N=64;
IntTimes=8; %迭代次数

for k=1:length(snr)
    InVar=1.0/snr(k);  %方差,码率1/2时使用
%     InVar=1.5/snr(k);  %方差,码率1/3时使用
    ErrorNum=0; %误码个数
    FrameNum=0; %帧数
    FrameErrorNum=0;    %误帧
%     if snr_db(k)<=.5
%         MaxFrameErrorNum=400;
%     elseif snr_db(k)<=1.0
%         MaxFrameErrorNum=100;
%     else
%         MaxFrameErrorNum=30;
%     end
    if snr_db(k)<=.5
        MaxFrameErrorNum=400;
    elseif snr_db(k)<=1.0
        MaxFrameErrorNum=200;
    elseif snr_db(k)<=1.5
        MaxFrameErrorNum=75;
    elseif snr_db(k)<=1.75
        MaxFrameErrorNum=30;
    else
        MaxFrameErrorNum=15;
    end
    while FrameErrorNum<MaxFrameErrorNum
        sourse=randint(1,2*N);
        %编码
        [VecOutA,VecOutB,Sc1,Sc2]=DouBinTurboEncode(sourse,N,1/2);
        %通过信道
        Ar=VecOutA(1:N)+sqrt(InVar)*randn(1,N);
        Yr1=VecOutA(N+1:2*N)+sqrt(InVar)*randn(1,N);
        Wr1=zeros(1,N); %码率1/2，全部W被凿孔
%         Wr1=VecOutA(2*N+1:3*N)+sqrt(InVar)*randn(1,N);
        Br=VecOutB(1:N)+sqrt(InVar)*randn(1,N);
        Yr2=VecOutB(N+1:2*N)+sqrt(InVar)*randn(1,N);
        Wr2=zeros(1,N); %码率1/2，全部W被凿孔
%         Wr2=VecOutB(2*N+1:3*N)+sqrt(InVar)*randn(1,N);
        %解码
        decoded=DouBinTurboDecode(Ar,Br,Wr1,Yr1,Wr2,Yr2,N,IntTimes,InVar,Sc1,Sc2);
        %统计误码
        [Num ber0]=biterr(sourse,decoded);
        if Num~=0
            FrameErrorNum=FrameErrorNum+1;
        end
        ErrorNum=ErrorNum+Num;
        FrameNum=FrameNum+1;
        if mod(FrameNum,10)==0
            sprintf('ML_MAP,snr=%fdB,%d Frames transmitted, %d bit errors,%d frame errors\n',snr_db(k),FrameNum,ErrorNum,FrameErrorNum)
        end
    end
    ber(k)=ErrorNum/(FrameNum*N*2);
    Fer(k)=FrameErrorNum/FrameNum;
    sprintf('Finally,ML_MAP,snr=%fdB,%d Frames transmitted, %d errors,bit error ratio is %f\n',snr_db(k),FrameNum,ErrorNum,ber(k))
end
% figure
% semilogy(snr_db,ber)