gallager_f.m

来自「density evolution for LDPC codes design 」· M 代码 · 共 8 行

% gallager_f.m
function y = gallager_f(x)
if x < 0.00001
  y = 12.21; % 12.21=gallager_f(0.00001)
else
  y = log((exp(x)+1.0)/(exp(x)-1.0));
end;