qinv_gf2.m

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

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);

stopcol=AE2EA(AE);

if stopcol==0
   invA=AE(:,rows+1:end);
else
   invA=[];
end