基于MATLAB的线性代数 判断齐次方程组有无穷多解

镇场诗：吾心为他人，此智不私藏。学术有小成，愿帮大家忙。诚心说愿望，愿见此文者得启发。

/////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////

%  x1+2*x2+2*x3+x4=0
%  2*x1+x2-2*x3-2*x4=0
%  x1-x2-4*x3-3*x4=0

a=[1 2 2 1;2 1 -2 -2;1 -1 -4 -3];
b=rref(a)

% b =
%    1 0 -2 -5/3
%    0 1  2  4/3
%    0 0  0   0

c=rank(a);
d=size(a);

if(d(2)>c)
e=d(2)-c;
fprintf('该齐次线性方程组具有无穷多解。\n任一解的通解式中含有%i个任意常数\n',e);
end

%   x1=2*x3+5/3*x4
%   x2=-2*x3-4/3*x4
%
%1 、让x3=0,x4=1  2、让x3=1,x4=0

//////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////

这个是一位同学在寒假中做的，比较系统。由于希望自己做的东西能够得到发扬，努力能获得认同。所以写到博客中，

不可否认的是，部分代码一定存在错误，如果不加研究，没有具体问题具体分析，胡乱复制粘贴，那么会造成严重的后果的。

如果有错误的话，请您指出来，一定改正。代码功能比较基础，比较好理解与分析。

在进行中，ilovematlab论坛及百度上热爱分享知识的各位给了很大的帮助。感恩。

希望大家看到可以得到启发。