计算方法B_高斯消去

Posted 2021-01-08 wupenghao95

%计算方法No.1
%20180916 by wupenghao
%高斯消去
%！！！循环中的步长一定要设置准确，+1和-1等，一定要注意！！！
A=rand(10,10);
b=rand(10,1);
x=A;
%消元
cof=zeros(10,1);
root=zeros(10,1);
for k=1:1:9
    for i=k+1:length(A)
        cof(i)=A(i,k)/A(k,k);
        b(i)=b(k)*(-cof(i))+b(i);
        for j=k:length(A)
        A(i,j)=A(k,j)*(-cof(i))+A(i,j);
        end
    end
end
%回代
n=length(A);
root(n)=b(n)/A(n,n);
for k=(n-1):-1:1
    S=b(k);
    for j=(k+1):1:n
        S=S-A(k,j).*root(j);
    end
    root(k)=S/A(k,k);
end