bp1.m

Mackey1随机构造LDPC码的校验矩阵进行编译码的仿真过程（其中译码算法采用置信传播算法）

%归一化BP算法
function [v]=BP1(y,H,sigma_2,maxiter)
% H=[1 1 1 0 0 0 0 0;0 0 0 1 1 1 0 0;1 0 0 1 0 0 1 0;0 1 0 0 1 0 0 1];
% y=[0.2 0.2 -0.9 0.6 0.5 -1.1 -0.4 -1.2];
% sigma_2=0.5;
% maxiter=50;
[m,n]=size(H);
success=0;k=0;                           
Lc=2*y/sigma_2;                                        %初始化信息节点的信息Lc.
Lg=repmat(Lc,m,1);                                     %初始化矩阵Lg(i,j)
Lg=Lg.*H;
[hj,hi]=find(H==1);
while((success==0)&(k<maxiter))                        %c*h'~=0或k(迭代次数)未达到最大迭代次数maxiter,继续进行迭代译码.
    for j=1:m                                          %计算校验节点向信息节点传递的消息Lh（j,i)
    for i=1:n
        if H(j,i)==1
            A=1; t=0;
            for ii=1:n
                if (ii~=i)&(H(j,ii)==1)
                    A=A*sign(Lg(j,ii)); b(j,ii)=abs(Lg(j,ii));t=t+1;
                    if t==1
                        B0=b(j,ii);
                    else B=b(j,ii);
                        if B0>B
                            B0=B;
                        end
                    end
                 end
            end
            x=1.25;Lh(j,i)=A*B0/1.25;
        end
    end
end
for i=1:n                                            %计算信息节点向校验节点传递的信息 Lg(i,j)
      rowind=find(hi==i);
      temp=Lh(hj(rowind),i);
      Q=sum(temp);
      Q1=Q-temp;
      Lg(hj(rowind),i)=Lc(i)+Q1;
      LQ=Lc(i)+Q;                                   %软判决
      if LQ<0                                       %硬判决
       v(i)=1;
      else v(i)=0;
      end
end
  k= k+1;
  if mod(v*H',2)==0                                  %判决译码是否成功,若成功(success=1),则退出循环运算,否则继续.
     success=1;
  else  success=0;
  end
end