置信度传播算法

% % MP algorithm% Qi,j(x)=1./(1+exp(-2*x*y/sigma^2)) for each i,j such that hi,j=1.% in code 0->-1,1->+1,such as q0 means q-1;p0 means P(X=-1|Y=y);% BPSK 0----->+1,1----->-1;function [out,success,iter]=ldpc_mp(y,sigma,H)s=size(H);q0=zeros(s);q1=zeros(s);r0=zeros(s);r1=zeros(s);p0=1./(1+exp(2*y/sigma^2));p1=1./(1+exp(-2*y/sigma^2));% % initilizationfor i=1:s(2) j=find(H(:,i)==1); q0(j,i)=p0(i); q1(j,i)=p1(i);enditer=0;success=0;max_iter=10;while ((success==0)&(iter0.5 out(i)=1; else out(i)=0; endendif rem(H*out.',2)==0 success=1; endend