example.m

码率自适应LDPC的MATLAB源程序

% Author:  David Varodayan ([e-mail])
% Date:    May 9, 2006

n = 396;                            % blocklength
ladderFile = '396_regDeg3.lad';     % regular degree 3 codes
pCrossover = 0.11;                  % BSC crossover probability

source = double(rand(1, n)>0.5);

accumSyndrome = encodeBits(source, ladderFile);

sideinfo_BSC = mod(source + double( rand(1, n)<pCrossover ), 2);
pCond = (1-pCrossover).*(1-sideinfo_BSC) + pCrossover.*sideinfo_BSC; % P(source=0|sideinfo)
LLR_intrinsic = log( pCond./(1-pCond) ); % log( P(source=0|sideinfo)/P(source=1|sideinfo) )

[decoded, rate, numErrors] = decodeBits( LLR_intrinsic, accumSyndrome, source, ladderFile );