用matlab求解非齐次线性方程组

function X = f(A,b)

% A 方程组的系数矩阵

% b 常数项对应的列向量

n = size(A,2);

B =[A,b];

RA = rank(A);

RB = rank(B);

if RA==RB & RA ==n %判断是否有唯一解

X =A\b;

elseif RA == RB & RA<n %判断是否有无穷解

[L,U] = lu(A);

X = U\(L\b);

C = null(A,'rational')

else

X = 'error:no solution'; %无解

end

end