delgauss.m

来自「matlab 多个计算方程租的算法汇集 主要讲的是一般方法 还有正定矩阵的解法」· M 代码 · 共 29 行

function x=DelGauss(a,b)
% Gauss消去法
[n,m]=size(a);
nb=length(b);
det=1;%存储行列式值
x=zeros(n,1);
for k=1:n-1
    for i=k+1:n
        if a(k,k)==0
            return
        end
     m=a(i,k)/a(k,k);
     for j=k+1:n
         a(i,j)=a(i,j)-m*a(k,j);
     end
     b(i)=b(i)-m*b(k);
 end
 det=det*a(k,k);
end
det=det*a(n,n);
    

for k=n:-1:1  %回代
    for j=k+1:n
        b(k)=b(k)-a(k,j)*x(j);
    end
    x(k)=b(k)/a(k,k);
end