作者Justinok (Sin)
看板comm_and_RF
标题[问题]J.Richardson的类下三角LDPC快速编码
时间Wed Mar 21 17:57:37 2012
小弟最近在验证由2001年 Thomas J.Richardson and Rudiger L. Urbanke所提出的
Efficient Encoding of Low-Density Parity-Check Codes这篇论文里头所提到的
类下三角LDPC快速编码,不晓得有没有先进验证过里头的公式可以分享一下经验,
小弟利用excel写出该公式流程,并输入 m = 1 0 0 0 0 0,经过运算之後的确会符合
Hu^T = 0的条件没错,但是如果我把 m 值改掉,例如:1 0 0 1 0 0 ,Hu^T 就不为 0了,
我验证过公式流程并没有错,H矩阵也按照论文上的经过高斯消去、行列转换後的图照
KEY,但出来就是不对,再请有经验的先进帮我解答。感谢
m是我传送讯息
u是我编码後的码字
补上公式:p1^T = (ET^-1A+C)m^T
p2^T = T^-1[Am^T + Bp1^T]
设:m = 1 0 0 1 0 0,
Am^T = 1 , Cm^T = 0 , T^-1 = 1 0 0 0 ET^-1 = 0 1 0 1
1 0 0 1 0 0 0 0 1 1
0 1 1 1 0
1 1 0 1 1
ET^-1(Am^T) = 0 , phi^-1*0 p1 = 1 1
1 1
Bp1^T = 1 , Tp2^T = 0 , p2 = 0 0 1 0
1 0
1 1
0 1
u = 1 0 0 1 0 0 1 1 0 0 1 0
Hu^T = 0 0 0 0 0 1 1
另外我有两个问题
1.PAPER里头有提到ordering 1,2,3,4,5,6,7,10,11,12,8,9放置到
H矩阵中,在PAPER范例中经过非奇异col"5",col"8"变换,之後可看为
1,2,3,4,10,6,7,5,11,12,8,9,我不晓得它这行数字代表甚麽? col"5"
以及"8"又是依据甚麽所以对换?
2.当D转为非奇异矩阵之後,他又可得一个等价矩阵,它是经过类似
互斥或转换後才得到全新的C D E吗? 感谢先进们的指导
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 120.96.59.197
1F:→ ato324:他的文章应该没错,我以前也验证过140.113.144.175 03/21 18:33
2F:→ Justinok:我推他的论文范例出来也对,可是码字改掉 120.96.59.197 03/21 18:37
3F:→ Justinok:就不行 120.96.59.197 03/21 18:37
4F:→ Justinok:留下我的msn,请有经验的提供您宝贵的 120.96.59.197 03/21 18:39
※ 编辑: Justinok 来自: 120.96.59.197 (03/21 20:01)
6F:→ lovewa:不好意思我没仔细看内容,但我想知道你所 140.113.144.28 03/22 11:21
7F:→ lovewa:谓的改掉u的意思是什麽?产生另外一组inf. 140.113.144.28 03/22 11:22
8F:→ lovewa:去encode吗?还是直接针对某组u直接改..? 140.113.144.28 03/22 11:22
9F:→ lovewa:如果是後者,那当然有可能发生Syndrom不为0 140.113.144.28 03/22 11:23
10F:→ lovewa:回在推文,可以麻烦你看一下吗? 140.113.144.28 03/22 11:23
※ 编辑: Justinok 来自: 120.96.59.197 (03/22 11:43)
11F:→ Justinok:对不起我上面没改到,我加上注解 120.96.59.197 03/22 11:44
12F:→ Justinok:原始传送讯息m是paper上的 1000000 120.96.59.197 03/22 11:44
13F:→ Justinok:我把原始讯息m改为 1 0 0 1 0 0 120.96.59.197 03/22 11:45
※ 编辑: Justinok 来自: 120.96.65.124 (03/22 14:26)
※ 编辑: Justinok 来自: 120.96.65.124 (03/22 15:45)
※ 编辑: Justinok 来自: 120.96.59.197 (03/28 20:35)
※ 编辑: Justinok 来自: 120.96.59.197 (03/30 10:58)
※ 编辑: Justinok 来自: 120.96.59.197 (03/30 11:02)