MATLAB--练习5（矩阵解方程组）

function [x ,y]=leqsolve(A,b)

    if nargin <2
        return
    end
    [m,n]=size(A);

    if norm(b)>0
        if rank(A)==rank([A,b])%系数矩阵的秩等于增广矩阵
            if rank(A)==m      %系数矩阵的秩等于行唯一解
                x=A\b;
            else%无穷解
                x=A\b;%特解
                y=null(x,'r');%基础解
            end
        else
            x=A\b;
        end

    else
        if rank(A)>=n
            x=zeros(m,1);
        else
            x=null(A,'r');
        end
    end
end 

 

。。。