正确因子表回代程序.m

本程序用于电力系统计算的因子表形成、前代、回代

%本程序的功能是利用因子表对常数项进行回代
A=input('请输入矩阵：A=');
B=input('请输入常数项矩阵：B=');
[n,m]=size(A);
for i=1:n
    A(i,i)=1./A(i,i);
    for j=i+1:n
        A(i,j)=A(i,j)*A(i,i);
    end
    for k=i+1:n
        for j=i+1:n
            A(k,j)=A(k,j)-A(k,i)*A(i,j);
        end
    end
end
disp('矩阵A的因子表为：A=');
disp(A)
for i=1:n
    B(i)=B(i)*A(i,i);
    for j=i+1:n
        B(j)=B(j)-A(j,i)*B(i);
    end
end
for i=n-1:-1:1
    for j=i+1:-1:2
        B(j-1)=B(j-1)-A(j-1,i+1)*B(i+1);
    end
end
disp('在因子表的基础上求解线性方程组的解为：X=');
disp(B)