matlab入门 线性方程组求解

Ax=B类型方程组

一、方程组有唯一解

A = [2 1 -5 1;1 -3 0 -6;0 2 -1 2;1 4 -7 6]
B = [8;9;-5;0]

%矩阵的秩
rankA = rank(A)
rankB = rank([A B])
%rankA = rankB = 4 = 4(未知数) 所以唯一解

%使用克拉默法则求解
D = det(A)                          %求系数行列式的值，不等于0时有唯一解
D1 = det([B A(:,2:4)])
D2 = det([A(:,1) B A(:,3:4)])
D3 = det([A(:,1:2) B A(:,4)])
D4 = det([A(:,1:3) B])
x1 = [D1/D;D2/D;D3/D;D4/D]

%逆矩阵求解
invA = inv(A);
x2 = invA*B

%使用左除
x3 = A\B;

%使用linsolve
x4 = linsolve(A,B)

二、方程组有无穷多解

A = [1 2 3 4;2 2 1 1;2 4 6 8;4 4 2 2]
B = [1;3;2;6]

%矩阵的秩
rankA = rank(A)
rankB = rank([A B])
%rankA = rankB = 2 < 4(未知数) 所以有无穷多解

%求齐次方程组的基础解系
null(A,'r')     %齐次解
pinv(A)*B       %特解

%求齐次方程组的基础解系(使用符号变量)
syms a1 a2
z = null(sym(A))
x0 = sym(pinv(A))*B
x = a1* z(:,1) + a2*z(:,2) + x0