Ldpc编译码方案MATLAB仿真

主要对LDPC码的编译码算法进行研究。首先,介绍LDPC的相关基本概念，其次,阐述了LDPC码的性能特点、发展应用以及LDPC码的基本理论知识。最后在译码算法上,简单对BP译码算法进行了介绍和仿真分析

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
% Bachelor of Science Graduation Dissertation
% ACE algorithm for Constructing LDPC Codes
% Copyright （C） 2009 Duanmu Fanyi BIT
% Code Designed by Student Duanmu Fanyi （Reg. No. 20050042）
% Directed by Fan Guangrong
% Class 01510500 Department of Electronic Engineering 
% School of Information Science Beijing Institute of Technology
% Reference:
% T Tao C Jones Villasenor R.D. and R Wesel. Construction of Irregular LDPC Codes with Low Error Floors [J]. IEEE
% Trans. Comm 2003 3: 3125-3129.
% For simplicity and efficiency the M and N is set small just for algorithm test

% Initialization
% “m” indicate the row number and “n” indicate the column number
% “dc” indicate the maximum number of “1”s in one column and “dr” indicate the maximum number of “1”s in one %row
% row_sum counts the number of “1”s in one row
% H1 and H2 are submatrix of final H matrix
% rank_num is used to check whether H2 is full-ranked or not
% cir is used to count the actual effective number of target column vectors
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

m=100;
n=200;
k=100;
R=0.5;
dc=3;
dr=6;

row_sum=zeros（m1）;       
Column_Vector=zeros（m1）; 
H1=zeros（mk）;
H2=zeros（mm）;
H=[H1H2];

row_num=m;
col_num=n;

rank_num=1;
cir=0;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
while（1）
Column=col_num-cir;
% redo:
% detect row degree
   row_sum=zeros（m1）;
   for c1=1:m
    for c2=1:n
     if H（c1c2）~=0
      %temp1=c1;
      %temp2=c2;
       row_sum（c11）=row_sum（c11）+1;
     end;
    end;
   end;

% Get Target Rows with smallest degree
    [v index]=sort（row_sum）;

% Create Vi
% check bits randomly generation
   if Column>m
      a=[1:m].‘;
      for i=1:m
       b=fix（rand（1）*m）+1;
       tmp=a（i）;
       a（i）=a（b）;
       a（b）=tmp;
      end
      Column_Vector=zeros（m1）;
      for k=1:3
       Column_Vector（a（k）1）=1;
      end;
% informationo bits uniformly generation
   else
      Column_Vector=zeros（m1）; 
       for k=1:3 
         Column_Vector（index（k）1）=1;
       end;
   end;
% obtain final H matrix
   for row_num=1:m
    H（row_numColumn）=Column_Vector（row_num1）;
    H2=H（:6:n）;
    H1=H（:1:m）;
   end;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% ACE processing
  if Column>m      %vi is a parity bit（check bit）    
% detect H2’ rank. If it is full continue to ACE detection; If it is not full go back to redo and generate another column % vector
    GE_num=rank（H2）;
    if GE_num~=rank_num 
      continue; %goto redo
    end;
  end;

% ACE detection for vi
% if ACE not applicable goto redo; if ACE applicable actual circulation increases and the H2’s rank increases
  cir=cir+1;
  rank_num=rank_num+1;
% If and only if circulation reaches the target column number stop; else continue to generate new column vectors
  if cir==n
    br

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-10-10 15:39  源码\
     目录           0  2016-03-06 11:15  源码\1\
     目录           0  2016-03-06 11:09  源码\1\H构造\
     文件        3076  2014-02-19 14:00  源码\1\H构造\ACE.m
     目录           0  2016-03-06 11:10  源码\1\LDPC\
     文件        1697  2007-05-04 15:06  源码\1\LDPC\decodeBitFlipping.m
     文件        2568  2007-05-04 15:09  源码\1\LDPC\decodeLogDomain.m
     文件        2169  2007-05-04 15:06  源码\1\LDPC\decodeLogDomainSimple.m
     文件        2468  2007-05-04 15:06  源码\1\LDPC\decodeProbDomain.m
     文件        2490  2012-06-12 15:51  源码\1\LDPC\ldpcBER.asv
     文件        2495  2013-05-27 09:30  源码\1\LDPC\ldpcBER.m
     文件        3498  2008-12-11 11:21  源码\1\LDPC\makeLdpc.asv
     文件        3468  2007-05-04 15:15  源码\1\LDPC\makeLdpc.m
     文件        3819  2007-05-04 14:25  源码\1\LDPC\makeParityChk.m
     文件         167  2016-03-06 11:15  源码\1\说明.txt
     目录           0  2018-10-10 15:39  演示\
     文件    12371456  2016-03-20 19:49  演示\9582录像2.avi
     文件          64  2018-10-10 16:34  说明.txt
     目录           0  2018-10-10 15:38  需求\
     文件        9962  2018-10-10 18:57  需求\MLFBO[T@OV]IZ2DZVAMOX2J.png