逆矩阵.txt

来自「逆矩阵 行列式的值 ｍａｔｌａｂ 数值计算方法 源程序」· 文本 代码 · 共 22 行

function  B=ni(A)
[n,m]=size(A);
I=eye(n);
if n==m
    for i=1:n
        p=A(i,i);
        A(i,:)=A(i,:)/p;
        I(i,:)=I(i,:)/p;
        for k=1:n
            if A(i,i)~=0
            L=A(k,i)/A(i,i);
            if k~=i
        A(k,:)=A(k,:)-A(i,:)*L;
        I(k,:)=I(k,:)-I(i,:)*L;
            end
        end
        end
    end
B=I;
else
    B=zeros(n);
end