比特翻转译码

tic
clc
clear all
count=1e4;  %测试次数，可自己选择
%初始化
load（‘H.mat‘）;
load（‘G.mat‘）;
for SNR=1:6

    total_err=0;  %错误比特数初始化为0
    for loop=1:count
        x = randi（[01]1252）;
        x_send=mod（x*G2）;     %编码
        x_bpsk=2*x_send-1;    %BPSK调制
        re_y=awgn（x_bpskSNR‘measured‘）;    %进入高斯信道
        max_ite=20;      %最大迭代次数，可自己选择
        [rowcol]=size（H）;
        iteration=0;  
        hard=0.5*（sign（re_y）+1）;
        y_re = hard;
        syn = mod（y_re*H‘2）;          % 计算伴随式（校正子）syn
        while （sum（syn） ~= 0） & （iteration < max_ite）   %循环条件syn=0 或达到最大迭代次数
            iteration = iteration + 1;
            
            %Step 2.  计算每一比特参与校正子为1的方程的个数fn
            fn=zeros（11）;
            for i = 1:col