BF比特翻轉譯ldpc碼

clc
clear all                          % y=[1,1,0,1,1,0,0]
load('H.mat') ; 
load('G.mat') ;                             %[1,0,1,1,0,0,0;0,1,0,0,1,0,0;1,1,0,0,0,1,0;0,0,1,0,0,0,1];
u=randi([0,1],1,84);
y=u*G;
Dim=size(H);
col=Dim(2);
iteration=20;   %迭代次數
for i=1:iteration
s=mod(H*y',2);
if s==0
    y=y;
    break
else 
    for n=1:col
    f(n)=s'*H(:,n);
    end
    for n=1:col
        if f(n)>=2
            y(n)=mod(y(n)+1,2);
        end
    end
end      
end