LDPC源代码

代码详细，步骤清晰，想学习LDPC的可以参考。

function [H] = alist2sparse（fname）
% reads binary parity check matrix in “alist“ format from file FNAME and
% converts it to sparse matrix used in MATLAB routines.
% This is an interface to matrices at http://wol.ra.phy.cam.ac.uk/mackay/codes/
%
% Example
%        [H] = alist2sparse（‘A‘）;   % A is the ascii file in alist format


%   Copyright （c） 1999 by Igor Kozintsev igor@ifp.uiuc.edu
%   $Revision: 1.1 $  $Date: 2000/03/23 $ Bug fixed by Hatim Behairy

fid = fopen（fname）;
n = fscanf（fid‘%d‘1）;
m = fscanf（fid‘%d‘1）;
maxinrow = fscanf（fid‘%d‘1）; 
junk = fscanf（fid‘%d‘1）; % no need
num = fscanf（fid‘%d‘[1 n]）; % number of elements in rows

num2（1:n）=maxinrow;    
junk = fscanf（fid‘%d‘[1 m]）; % no need

position = zeros（nmaxinrow）;
for i=1:n
   for j=1:num2（i）    
      position（ij） = fscanf（fid‘%d‘1）;
   end
end

ii = zeros（1sum（num））;
jj = ii;
k = 1;
for i=1:n
      for j=1:num（i）
      jj（k） = i;
      ii（k） = position（ij）;
      ss = 1;
      k = k+1 ; 
   end
end
H = sparse（iijjssmn）;
fclose（fid）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        446  2000-05-09 18:30  ldpc_h2g.m

     文件       1079  2000-05-09 18:30  alist2sparse.m

     文件      12536  2000-06-05 13:11  ldpc_decode.m

     文件       1258  2000-05-09 18:30  ldpc_encode.m

     文件       8984  2000-05-09 18:30  ldpc_generate.c

     文件        712  2000-05-09 18:30  ldpc_generate.m

     文件      12505  2000-05-09 18:30  ldpc_h2g.c

     文件       9721  2008-12-16 19:13  使用帮助：新手必看.htm

     文件        183  2008-09-03 11:09  Matlab中文论坛--助努力的人完成毕业设计.url

----------- ---------  ---------- -----  ----

                47424                    9