大家好，又见面了，我是你们的朋友全栈君。
function [ out ] = inv_bin( in )
%INV_BIN Summary of this function goes here
% Detailed explanation goes here
% 计算二进制稀疏矩阵的逆
[m,n] = size(in);
if(m~=n)
fprintf('m~=n\n');
return ;
end
E = eye(m);
%%
%做行变换，变成下三角阵
for i = 1:m
noneZerosIndex = find(in(:,i)); %i到end行，第i列非零元
noneZerosIndex = noneZerosIndex(find(noneZerosIndex>=i));
if(length(noneZerosIndex)==0) %该行为0
randIndex = randi([i+1,m],1);
%列交换
temp = in(:,i);
in(:,i) = in(:,randIndex);
in(:,randIndex) = temp;
%E交换
temp = E(:,i);
E(:,i) = E(:,randIndex);
E(:,randIndex) = temp;
end
id1 = noneZerosIndex(1);
%in交换
temp = in(i,:);
in(i,:) = in(id1,:);
in(id1,:) = temp;
%E交换
temp = E(i,:);
E(i,:) = E(id1,:);
E(id1,:) = temp;
noneZerosIndex = find(in(:,i)); %第i列非零元,用其中的第一行消其他行
for cc = 1:length(noneZerosIndex)
if(noneZerosIndex(cc)~=i) %跳过第i行
temp = mod(in(noneZerosIndex(cc),:)+in(i,:) , 2);
in(noneZerosIndex(cc),:) = temp;
temp = mod(E(noneZerosIndex(cc),:)+E(i,:) , 2);
E(noneZerosIndex(cc),:) = temp;
end
end
end
out = E;
end
版权声明：本文内容由互联网用户自发贡献，该文观点仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请联系我们举报，一经查实，本站将立刻删除。
发布者：全栈程序员-站长，转载请注明出处：https://javaforall.net/129146.html原文链接：https://javaforall.net
