%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?
%?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