inv_gf2.m

来自「LDPC中的q进制编码算法」· M 代码 · 共 38 行

function [invA,stopcol]=inv_GF2(A);
%[invA,stopcol]=inv_GF2(A);

[rows cols]=size(A);

if rows~=cols
   error('Input matrix error: rows must equal cols!');
end

stopcol=0;

%AE为完全矩阵，如果AE为稀疏矩阵则当AE规模较大时速度会大大减慢
AE=[full(A) eye(rows)];

%对AE在GF(2)上进行初等行变换，前面的rows行 x rows列形成单位矩阵
for k=1:rows
   if AE(k,k)==0
      ind=find(AE(k+1:end,k)~=0)';
      if isempty(ind)
         invA=[];
         stopcol=k;
         return
      end
      AE(k,:)=rem(AE(k,:)+AE(k+ind(1),:),2);
   end

   ind=find(AE(:,k)~=0)';
   for x=ind
      if x~=k
         AE(x,:)=rem(AE(x,:)+AE(k,:),2);
      end
   end

end

%后面的E经同样的初等行变换得到的矩阵即为invA
invA=AE(:,rows+1:end);