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?????????????iteration?=?iteration?+?1;
????????????
????????????%Step?2.??計算每一比特參與校正子為1的方程的個數fn
????????????fn=zeros（11）;
????????????for?i?=?1:col
????????