ssolve.m

Matlab中虽然有很多解方程的函数

%解一般线性方程组
function [X,phi]=ssolve(A,b)
[m,n]=size(A);
if nargin==1
    b=zeros(m,1);
end
AA=[A,b];
r=rank(AA);
if rank(A)~=r
    disp('无解！')
else
    [AA,o]=srref(AA);
    AA=AA(1:r,:);
    if r==n
        X=[];
        phi=permut(AA(:,n+1),o);
    else
        X=permut([-AA(:,r+1:n);eye(n-r)],o);
        phi=permut([AA(:,n+1);zeros(n-r,1)],o);
    end
end