isexist.m

MATLAB科学计算与工程分析源代码源程序4

function y=isexist(A,b)
% 该函数用来判断线性方程组Ax=b的解的存在性
% 若方程组无解则返回0,若有唯一解则返回1,若有无穷多解则返回Inf
% 编写日期：2007-5-14

[m,n]=size(A);
[mb,nb]=size(b);
if m~=mb
    error('输入有误！');
    return;
end
r=rank(A);
s=rank([A,b]);
if r==s&r==n
    y=1;
elseif r==s&r<n
    y=Inf;
else
    y=0;
end