高斯消元法

ID:334306 · 发表于 2018-5-20 17:03

function x=gauss_elim(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

帐号		自动登录	找回密码
密码			立即注册