%?%?MP?algorithm
%?Qij（x）=1./（1+exp（-2*x*y/sigma^2））?for?each?ij?such?that?hij=1.
%?in?code?0->-11->+1such?as?q0?means?q-1;p0?means?P（X=-1|Y=y）;
%?BPSK??0----->+11----->-1;
function?[outsuccessiter]=ldpc_mp（ysigmaH）
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））;
%?%?initilization
for?i=1:s（2）
????j=find（H（:i）==1）;
????q0（ji）=p0（i）;
????q1（ji）=p1（i）;
end

iter=0;
success=0;
max_iter=10;
while?（（success==0）&（iter????iter=iter+1;
for?jj=1:s（1）
????init_q=1;
????ii=find（H（jj:）==1）;
????for?k=1:length（ii）
????init_q=init_q*（1-2*q0（jjii（k）））;
????end
????r1（jjii）=0.5+0.5*init_q./（1-2*q0（jjii））;
????r0（jjii）=1-r1（jjii）;
end
for?iii=1:s（2）
????init_r1=1;
???