ffsolvtion.m

《MATLAB 7.X程序设计》-王建卫-源代码及习题答案

function X=ffsolvtion(A,b,n)
%n方程组中未知变量的个数
B=[A,b];                              %增广矩阵
R_A=rank(A);                         %求系数矩阵的秩
R_B=rank(B);                         %求增广矩阵的秩
format rat
if R_A==R_B&R_A==n                  %判断是否有唯一解
   X=A\b;                             %求唯一解
elseif R_A==R_B&R_A<n              %判断是否有无穷解
   X=A\b;                             %求特解
   C=null(A,'r')                     %求AX=0的基础解系
else 
X='Equation has no solves';       %判断是否无解
end