LDPC仿真程序

%BCJR算法G=[1 （1+D+D2）/（1+D2）]
%x：2*n接收码字
%La：先验概率
%isinterweave：是否收尾
function Le=BCJR（xLaLcisinterweave）
    n=length（La）;
    next_net0=[1 3 2 4];
    next_net1=[2 4 1 3];
    last_net0=[1 3 2 4];
    last_net1=[3 1 4 2];
    net_out0=[-1 1 -1 1];
    net_out1=[1 -1 1 -1];
    %计算分支度量
    for i=1:n
        for j=1:4
            gama0（ij）=-La（i）/2+Lc/2*（x（2i）*net_out0（j）-x（1i））;%从j状态出发码字为0的路径
            gama1（ij）=La（i）/2+Lc/2*（x（2i）*net_out1（j）+x（1i））;%从j状态出发码字为1的路径
        end
    end
    %计算前向度量
    a（11）=0;
    a（21）=gama0（11）;
    a（22）=gama1（11）;
    a（31）=a（21）+gama0（21）;
    a（32）=a（21）+gama1（21）;
    a（33）=a（22）+gama0（22）;
    a（34）=a（22）+gama1（22）;
    for i=4:n+1
        for j=1:4
            a（ij）=max（a（i-1last_net0（j））+gama0（i-1last_net0（j））a（i-1last_net1（j））+gama1（i-1last_net1（j）））; 
            %到达j状态的前向度量
        end
    end
    if（isinterweave==0）
        a（n2）=0;
        a（n3）=0;
        a（n4）=0;
        a（n-12）=0;
        a（n-14）=0;
        %计算后向度量
        b（n+11）=0;
        b（n1）=gama0（n1）;
        b（n3）=gama1（n3）;
        b（n-11）=b（n1）+gama0（n1）;
        b（n-12）=b（n3）+gama0（n3）;
        b（n-13）=b（n1）+gama1（n1）;
        b（n-14）=b（n3）+gama1（n3）;
    else
        b（n+1:）=a（n+1:）;
        for j=1:4
            b（nj）=max（b（n+1next_net0（j））+gama0（nj）b（n+1next_net1（j））+gama1（nj））;
        end
    end
    for i=n-1:-1:1
        for j=1:4
            b（ij）=max（b（i+1next_net0（j））+gama0（ij）b（i+1next_net1（j））+gama1（ij））;
            %由j状态出发的后向度量
        end
    end
    b（12）=0;
    b（13）=0;
    b（14）=0;
    b（23）=0;
    b（24）=0;
    for i=1:n
        for j=1:4
            m（1j）=a（ij）+gama0（ij）+b（i+1next_net0（j））;
            m（2j）=a（ij）+gama1（ij）+b（i+1next_net1（j））;
        end
        Le（i）=max（m（2:））-max（m（1:））;
    end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2016-12-08 20:58  LDPC_Test
     文件        1993  2016-11-28 16:44  LDPC_TestBCJR.m
     文件         700  2016-11-30 16:15  LDPC_TestBuildHG_80211n.m
     文件         537  2016-12-01 08:47  LDPC_TestBuildHG_CCSDS.m
     文件         378  2016-11-30 16:17  LDPC_TestBuildHG_random.m
     文件         822  2016-11-28 16:44  LDPC_TestH2G.m
     文件        1450  2016-12-02 16:55  LDPC_TestLDPC.m
     文件        1240  2016-12-06 15:08  LDPC_TestTurbo.m
     文件         364  2016-12-05 10:23  LDPC_Testconvcoding.m
     文件         214  2016-11-28 16:44  LDPC_TestexchHG.m
     文件        4998  2016-12-01 08:47  LDPC_TestgenBM_CCSDS.m
     文件        4111  2016-11-28 16:44  LDPC_TestgetH.m
     文件        3596  2016-12-02 15:27  LDPC_Testldpc_decode.m
     文件        1162  2016-12-02 16:47  LDPC_Testldpc_decode_LLR.m
     文件          51  2016-11-28 16:44  LDPC_Testldpc_encode.m
     目录           0  2016-12-01 08:47  LDPC_Testmatrix
     文件     2927403  2016-11-30 14:28  LDPC_TestmatrixG10240_25.mat
     文件       17096  2016-11-29 09:50  LDPC_TestmatrixG1296_12.mat
     文件       34211  2016-11-29 09:51  LDPC_TestmatrixG1944_12.mat
     文件       35967  2016-11-29 18:59  LDPC_TestmatrixG1944_23.mat
     文件       35467  2016-11-29 19:00  LDPC_TestmatrixG1944_34.mat
     文件       31573  2016-11-29 19:01  LDPC_TestmatrixG1944_56.mat
     文件       68111  2016-11-30 09:08  LDPC_TestmatrixG2560_25.mat
     文件      157209  2016-11-29 09:31  LDPC_TestmatrixG4512_12.mat
     文件      739721  2016-11-30 09:24  LDPC_TestmatrixG5120_25.mat
     文件        4877  2016-12-01 08:47  LDPC_TestmatrixG640_25.mat
     文件        6558  2016-11-29 09:46  LDPC_TestmatrixG648_12.mat
     文件       17373  2016-11-28 18:14  LDPC_TestmatrixG648_12r.mat
     文件      352052  2016-11-30 14:28  LDPC_TestmatrixH10240_25.mat
     文件        9133  2016-11-29 09:50  LDPC_TestmatrixH1296_12.mat
     文件       17559  2016-11-29 09:51  LDPC_TestmatrixH1944_12.mat
............此处省略55个文件信息