基于MATLAB的线性代数 非齐次方程组的基础解系和通解

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

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

%系数矩阵
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(c)

% f =
%
%      1     0    -1    -1
%      0     1    -1    -2
%      0     0     0     0
%有了这个后续工作就好求了

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

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

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

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

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

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