function?x_hat=bp_dec（rxHiterate）??????????%應該可以使用稀疏矩陣實現
%Simplified?log-BP?decoder?of?LDPC?
%rx?is?the?received?code?which?carries?the?soft?information
%H?is?the?LDPC?matrix
%iterate?is?the?iterating?times
%Author:LEE?Ming-yang100926?
%Mail:wo.lmy@163.com??QQ:905432530
lcj=-rx;?????????????????????????????????????%后驗似然信息
[M?N]=size（H）;
lrij=zeros（MN）;?????????????????????????????%初始化校驗信息
lqij=H.*repmat（lcjM1）;?????????????????????%關鍵，體現了校驗節點和變量節點在H陣中的關系
alphaij=zeros（MN）;??????????????????????????%sign（lqij）
betaij=zeros（MN）;???????????????????????????%lqij
lqj=zeros（1N）;
for?it=1:iterate
????alphaij=sign（lqij）;??????????????????????%i行j列，迭代之前賦值
????betaij?=abs（lqij）;
????%水平處理，得到校驗節點信息
????for?i=1:M????????????????????????????????%遍歷行
????????cl=find（H（i:））;?????????????????????%從某行中找非0列
????????for?j=1:length（cl）???????????????????%j都遍歷相關聯的變量節點
????????????min_beta=realmax;????????????????%每次使用之前初始化
????????????temp_beta=betaij（icl）;??????????%只要有用的，無關量是0，會干擾判斷
????????????temp_beta（j）=realmax;????????????%j個對應的那列不用于比較
????????????min_beta=min（temp_beta）;
????????????%betaij（ij）=min_beta;??????????%只是個標志，沒什么用
????????????temp_alpha=alphaij（icl）;
????????????temp_alpha（j）=1;
????????????prod_alpha=prod（temp_alpha）;
????????????lrij（icl（j））=prod_alpha*min_beta;
????????end
????end
????%垂直處理，得到變量節點信息
????for?j=1:N
????????rl=find（H（:j））;
????????for?i=1:length（rl）
????????????lqij（rl（i）j）=lcj（j）+sum（lrij（rlj））-lrij（rl（i）j）;
????????end
????????%累加處理
????????lqj（j）=lcj（j）+sum（lrij（rlj））;
????????if?lqj（j）>0
????????????x_hat（j）=0;
????????else
????????????x_hat（j）=1;
????????end
????end
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????6722??2015-10-19?21:52??H4_5.mat

?????文件???????1031??2010-09-24?13:02??inv_GF2.m

?????文件????????868??2015-10-19?23:34??Test_CCSDS_Matrix.m

?????文件???????1843??2010-10-12?10:30??bp_dec.m

?????文件????????253??2015-10-19?21:54??G_gen.m

?????文件??????14261??2013-03-11?14:53??H2_3.mat

-----------?---------??----------?-----??----

????????????????24978????????????????????6