作者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)