Hi, I am trying to make a parity check matrix from non-systematic to systematic. Hence I am attaching my code below. Somewhat it is correct, but there are some problems. It would be really great if someone could help me in this. Thanks

10 次查看（过去 30 天）
sb
sb 2016-6-9
评论： SP22 2017-11-20
Information theory and coding. I am working on LDPC coding and decoding. Please check the code below
  2 个评论
sb
sb 2016-6-10
编辑：Walter Roberson 2017-5-25
For e.g.
This is my H matrix
H =
1 0 1 1 0
0 0 1 0 1
1 0 0 1 0
1 0 1 1 1
I want to have an identity matrix inside the matrix. the dimension on H matrix here is (mxn) which is (4x5).
I should have matrix as this in the result:
Hsys =
0 1 0 0 0
0 0 1 0 0
1 0 0 1 0
0 0 0 0 1
I should have an identity matrix of dimension 'm'

请先登录，再进行评论。

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

回答（1 个）

freebil
freebil 2016-6-10
You have to do gauss jordan elimination to convert a parity check matrix to upper triangular form. For example,
H=[1 1 0 1 1 0 0 1 0 0;
0 1 1 0 1 1 1 0 0 0;
0 0 0 1 0 0 0 1 1 1;
1 1 0 0 0 1 1 0 1 0;
0 0 1 0 0 1 0 1 0 1];
There is rref() in matlab and you have to do it in GF(2). So,
HH = mod(rref(H),2)
gives
HH = 1 0 0 0 0 0 1 1 1 0
0 1 0 0 0 1 0 1 0 0
0 0 1 0 0 1 0 1 0 1
0 0 0 1 0 0 0 1 1 1
0 0 0 0 1 1 1 0 0 1
  2 个评论
SP22
SP22 2017-11-20
This might help you to get parity check matrix in form H=[A In-k].
temp=HH(:,1:5) %Access the identity matrix through column
temp2=HH(:,6:10) %Acces the parity through column
Hsyst=horzcat(temp2,temp) %Actual systematic matrix in the form H((n-k)*n)=[P' :I(n-k)]

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 AI for Wireless 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by