matlab : R2018a 64bit
OS : Windows 10 x64
typesetting : Markdown
blog : my.oschina.net/zhichengjiu
gitee : gitee.com/zhichengjiu
code
clear
clc
% 系数矩阵
a=[-2 1 1;1 -2 1;1 1 -2];
b=[0;3;-3]
c=[a b] % 增广矩阵
r_a=rank(a);
r_c=rank(c);
d=size(a);
if(r_a==r_c)
fprintf('方程相容\n');
if(r_a<d(2))
e=d(2)-r_a;
fprintf('该齐次线性方程组具有无穷多解。\n任一解的通解式中含有%i个任意常数\n',e);
end
end
f=rref