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