How to generating LDPC parity-check matrices with arbitrary sizes?

Question

Hamid Reza Barzegar 2017-1-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/318834-how-to-generating-ldpc-parity-check-matrices-with-arbitrary-sizes

评论： han xiangqun 2024-5-13

I want to use comm.LDPCEncoder/comm.LDPCDecoder and i know to use this encoder decoder in Matlab the H matrices should be Systematic (support Systematic codes only) otherwise i will receive this error.

"the last (N-K) columns of the parity-check matrix must be invertible in GF(2)."

I know two methods from MATLAB that will generate parity-check matrices:

H = dvbs2ldpc(r)
h = hammgen(m)

However, these methods are restricted to a certain ratio between rows and columns. I tried some other third-party scripts to generate the matrix, but some of them don't seem to return a valid matrix, since "comm.LDPCDecoder" will fail with the same error.

I used some predefined H matrices of Mackay, but no way.

lent of my code-word is not fix it's about 1000 to 3000 it depend to the situation of my simulation(for example code-word length is equal to 1536 or another case 2048), that's why i should create a spars matrix then use it.

Now how to generate the LDPC parity-check matrices with arbitrary sizes[depend to length of my code-word]?

How to generating LDPC parity-check matrices with arbitrary sizes?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

David Ding 2017-1-5

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/318834-how-to-generating-ldpc-parity-check-matrices-with-arbitrary-sizes#answer_249510

在 MATLAB Online 中打开

Hi Hamid,

I believe you might find the following link useful:

http://mathforum.org/kb/message.jspa?messageID=9818971

Essentially, you just need to do a bit of pre-processing on the parity check matrix in permuting its columns, such that the last (n-k) columns of the new H is always invertible in F_2. To do this, as the author in the link suggested, you may use:

% Permute H to obtain a systematic code
% This is always possible if H has full rank and if infobits_loc is
% correctly defined
paritybits_loc = setdiff(1:size(H,2),infobits_loc);
H2 = H(:,[infobits_loc paritybits_loc]);

Hope this helps.

David