线性代数 MATLAB求解线性方程组的解

clear all;clc
%% 2.36矩阵求逆
A=[3 2 1;3 1 5;3 2 3];
if det(A)==0
    disp('A矩阵不可逆');
else
    disp('A矩阵可逆');A1=inv(A)
end

B=[0 2 -1;1 1 2;-1 -1 -1];
if det(B)==0
    disp('B矩阵不可逆');
else
    disp('B矩阵可逆');B1=inv(B)
end

C=[2 2 3;1 -1 0;-1 2 3];
if det(C)==0
    disp('C矩阵不可逆');
else
    disp('C矩阵可逆');C1=inv(C)
end
%% 2.37线性方程组求解
D=[1 1 -1;0 2 2;1 -1 0];
E=[1 -1;1 1;2 1];
if det(D)==0
    disp('D矩阵不可逆');
else
    disp('D矩阵可逆');
    X1=inv(D)*E;
    disp(X1);
end

F=[4 1 -2;-1 0 1;1 -1 0];
E=[1 -3;2 2;3 -1];
if det(F)==0
    disp('F矩阵不可逆');
else
    disp('F矩阵可逆');
    X2=inv(F)*E;
    disp(X2);
end

%% 2.38线性方程组求解
G=[1 1 -1;0 2 2;1 -1 0];
E=[1 -1 1;1 1 0];
if det(G)==0
    disp('G矩阵不可逆');
else
    disp('G矩阵可逆');
    X3=E*inv(G);
    disp(X3);
end

H=[0 2 1;2 -1 3;-3 3 -4];
E=[1 2 3;2 -3 1];
if det(H)==0
    disp('H矩阵不可逆');
else
    disp('H矩阵可逆');
    X4=E*inv(H);
    disp(X4);
end