消息传递算法mp的matlab代码

% % 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;