clear?x;
x=demodz;

clear?prob0;
clear?prob1;
prob0=zeros（1lieH）;
prob1=zeros（1lieH）;

for?a=1:lieH
????prob1（1a）=1/（1+exp（x（1a）*2/sigma^2））;
????prob0（1a）=1-prob1（1a）;
end

q0=zeros（hangHlieH）;
q1=zeros（hangHlieH）;
r0=zeros（hangHlieH）;
r1=zeros（hangHlieH）;

for?i=1:lieH
????for?jo=1:liezhong
????????q0（ncol（joi）i）=prob0（1i）;
????????q1（ncol（joi）i）=prob1（1i）;
????end
end




for?lp=1:20
????for?i=1:hangH
????????for?jo=1:hangzhong
????????????delta=1;
????????????for?l=1:hangzhong
????????????????if?l~=jo
????????????????????delta=delta*（1-2*q1（inrow（il）））;
????????????????end
????????????end
????????????r0（inrow（ijo））=0.5+0.5*delta;
????????????r1（inrow（ijo））=1-r0（inrow（ijo））;
????????end
????end

????for?i=1:lieH
????????for?jo=1:liezhong
????????????alpha1=1;
????????????alpha2=1;
????????????alpha3=1;
????????????alpha4=1;
????????????for?l=1:liezhong
????????????????if?l~=jo
????????????????????alpha1=alpha1