matlab
优化模型算法
非齐次线性方程组的通解
A=[1 2 1; 2 -1 3; 3 4 -2];
b=[3 1 7]';
B=[A b];
n=4;
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='equition no solve' %判断无解
end
%结果展示
R_A =
3
R_B =
3
X =
1
1
-1/126884024284177488
C =
空的 3×0 double 矩阵