f.m

ldpc编码

function [y] = F(x)
if (x == 0.0) 
	y = 1.0e+30;               
elseif (x > 30.0)                          
	y = 0.0;                              
else                                         
	y = log ( (exp(x)+1.0)/(exp(x)-1.0) );
end