B站台大MATLAB教学视频自学整理---线性方程与线性系统【11】

B站视频链接：台大郭彦甫教授的教学视频

1.线性方程

1.1 高斯消除法

• 增广矩阵：就是把线性方程的常数项加入矩阵的最右侧。
• 高斯消除法： 就是求矩阵的秩。
• rref([A b]) : 高斯消除法消除行

1.2 LU消除法

• LU分解是指矩阵 可以分解为 $L^{-1}U$乘积的形式，其中$L$是单位下三角矩阵,$U$是单位上三角矩阵。
• lu(A) : 得到A的L矩阵和U矩阵
• 在matlab里用左除算矩阵：Ax=b --> x=A\b;
• mldivide(A ， B)：左除

练习1

%% 练习1
syms R1 R2 R3 R4 R5 V1 V2
A = [R1 0   0   R4  0  ;
     0  R2  0  -R4  R5 ;
     0  0  -R3  0   R5 ;
     1 -1   0  -1   0  ;
     0  1  -1   0  -1 ];

b = [V1 ; 0 ; V2 ; 0 ; 0];

x = A\b

矩阵分解函数

函数	释意
qr(A)	返回 QR 分解 A = Q*R 的上三角 R 因子
ldl()
ilu()
lu()
chol()
gsvd
svd

1.3矩阵求逆

• inv(A) : 求矩阵A的逆（A必须满秩）
• det(A) : 求矩阵A的行列式值。
  练习2
  

%% 练习2
[x,y] = meshgrid(0:1:40 , -40:1:100);%meshgrid 二维或三维网格
z1 = -x-y;
z2 = y-x;
z3 = (-1/3)*x;
mesh(x,y,z1)
hold on
mesh(x,y,z2)
hold on
mesh(x,y,z3)

1.4矩阵条件数

• cond : 计算矩阵条件数
• rank : 矩阵的秩

2.线性系统

$$y=Ab$$
$y$未知，为系统的输出。

• eig(A) : 算出A的特征值和特征向量
• expm(A) : 求以e为底的矩阵指数函数.