fai.m

LDPC 64QAM仿真程序

function E=Fai(belta)
index0=find(belta==0);
belta(index0)=1;
index1=find(belta>=30.0);
belta(index1)=1;
   %EObijindx(INDXN)=1e-100;
   %IIaijindx(INDXN)=1;
E=-log(tanh(belta/2));
E(index0)=50;
E(index1)=1e-10;