項目二：LDPC碼高斯信道譯碼

clc;clear;
degree_sum = 100;
for sum_loop = 1:degree_sum
    load PEGH_3_6_1008.mat;
    [m,n] = size(H);
    C = zeros(1,n);
    for i = 1:n
        if C(i)==1
            C1(i) = -1;
        else
            C1(i) = 1;
        end
    end
    C1 ;
    SNR = 2;
    b2 = 10^((-0.1)*(SNR));                           %方差
    Y = C1+(sqrt(b2))*randn(1,n);     %通过高斯信道被污染的信道
    %******************软判决begin*****************************
    Lch = 2*Y/b2;
    degree = 30;                       %迭代次数
    %*****************************迭代开始*************************************
    for Loop = 1:degree
        %**************************初始化开始***********************************
        if Loop == 1
            for j = 1:n
                for i = 1:m
                    if H(i,j)~=0
                        Lch_1(j,i) = Lch(j);
                    end
                end
            end
            Lch &