如何使用MATLAB从奇偶校验矩阵48X96导出LDPC生成矩阵

Question

我已经生成了LDPC稀疏奇偶校验矩阵，其中n = 96，M = 48，k = 48，速率为1/2。 enter link description here（N = 96，K = 48，M = 48，R = 0.5）

以下链接将显示H矩阵声明enter link description here

生成矩阵的MATLAB源代码

function [b]=inv_GF2(A)
 %Ainv=inv_GF2(A)
 dim=size(A);
 rows=dim(1);
 cols=dim(2);

 for i=1:rows
   for j=1:rows
     unity(i,j)=0;
   end

   unity(i,i)=1;

 end

 for i=1:rows
   b(1:rows,i)=gflineq(A,unity(1:rows,i));
 end

 H=[h1 h2];
 [m,n]=size(H);

 A(1:m,1:m)=H(1:m,1:m);
 B(1:m,1:m)=H(1:m,m+1:n);

 d=mod(inv_GF2(B)*A,2); 
 I=eye(m);
 G=[ I d'];

上述MATLAB代码显示的结果不正确。