gauss.m

来自「用mathlab实现Gauss算法,gauss算法是科学计算中一个非常著名的常用」· M 代码 · 共 34 行

function [result]=Gauss(A,b,x0)
%计算矩阵大小
n=length(b);
%初始化x1
x1=zeros(n,1);
while true
    %计算x1(1)
    x1(1)=b(1)/A(1,1);
    for j=2:n
        x1(1)=x1(1)-A(1,j)*x0(j)/A(1,1);
    end
    %计算x1(2)到x1(n-1)
    for i=2:n-1
        x1(i)=b(i)/A(i,i);
        for j=1:i-1
            x1(i)=x1(i)-A(i,j)*x1(j)/A(i,i);
        end
        for j=i+1:n
            x1(i)=x1(i)-A(i,j)*x0(j)/A(i,i);
        end
    end
    %计算x1(n)
    x1(n)=b(n)/A(n,n);
    for j=1:n-1
        x1(n)=x1(n)-A(n,j)*x1(j)/A(n,n);
    end
    
    if norm(x1-x0)<10^-6
        result=x1;
       break;
    else
        x0=x1;
    end
end