本文为一些矩阵相关知识的自习笔记。
矩阵特征值
若存在标量 δ 和非零向量 x 使矩阵满足方程 Ax=δx(x=0),则 δ 和 x 为矩阵的特征值和特征向量。
特征值 δ 通过方程 det(A−δI)=0 求解。
性质
- 行列式:矩阵的行列式等于所有特征值的乘积(det(A)=∏δi);
- 迹:矩阵的 迹(对角线元素之和) 等于所有特征值的和;
- 相似矩阵:若 A 和 B 为相似矩阵(即存在 P 使 P−1AP=B),则 A 和 B 具有相同的特征值;
- 特殊矩阵:三角矩阵的特征值就是其主对角线元素。
共轭转置矩阵
记作 AH 或 A∗。共轭转置即先将矩阵进行转置,再对每个元素取复共轭的操作。对于实数矩阵,共轭转置等同于普通转置。
矩阵的奇异值
A 的奇异值是 A∗A 的特征值 δi 的非负平方根 σi=δi。
矩阵范数
矩阵范数用于度量矩阵的大小。常见有 1-范数,∞-范数,Frobenius范数,2-范数。
MATLAB计算代码:norm(A, 1),norm(A, inf),norm(A, 'fro'),norm(A, 2)。
1-范数
∣∣A∣∣1:矩阵列向量的绝对值之和的最大值(列和最大)。
∞-范数
∣∣A∣∣∞:矩阵行向量的绝对值之和的最大值(行和最大)。
Frobenius范数
∣∣A∣∣F:矩阵元素绝对值平方和的平方根,即 ∑i,j∣aij∣2。
2-范数
∣∣A∣∣2:矩阵的最大奇异值,即 A∗A 最大特征值的平方根。
舒尔补
舒尔补通常用于将大的复杂的矩阵分解成小的易于解决的矩阵,更具体地,通常用于把含有未知变量相乘(非线性)的复杂不等式,拉伸成一个维度更高、但完全线性的矩阵不等式,从而能够直接求解。
舒尔补的定义
假设一个 (p+q)×(p+q) 的矩阵 M 被分为 A,B,C,D 四个部分,分别是 p×p、p×q、q×p 以及 q×q 的矩阵,也就是说:
M=[A CBD]
若 D 可逆,则 D 在矩阵内的舒尔补为 A−BD−1C;同理,若 A 可逆,则 A 在矩阵内的舒尔补为 D−CA−1B。
舒尔补的性质
以下以 A 可逆,S=D−CA−1B 为 A 的舒尔补为例。
行列式关系如下:
det(M)=det(A)×det(D−CA−1B)
正定性关系如下:
矩阵 M 正定,当且仅当 A 和其舒尔补 S=D−CA−1B 正定。
舒尔补在解方程中的应用
考虑一个方程 Mx=b,将 M 进行划分:
M=[A CBD],x=[x1 x2],b=[b1 b2]
通过第一行得到
x1=A−1(b1−Bx2)
代入到第二行得到
(D−CA−1B)x2=b2−CA−1b1
通过 A 和它的舒尔补解出方程。
舒尔补在求矩阵的逆时的应用
如果 M 可逆且 A 可逆,则:
M−1=[A CBD]−1=[(A−BD−1C)−1 −D−1C(A−BD−1C)−1−(A−BD−1C)−1BD−1D−1+D−1C(A−BD−1C)−1BD−1]
如果 M 可逆且 A,D 可逆,则:
M−1=[A CBD]−1=[(A−BD−1C)−1 −D−1C(A−BD−1C)−1−A−1B(D−CA−1B)−1(D−CA−1B)−1]
更特殊的,如果 A,B,C,D 都是系数(数字)时,有:
M−1=AD−BC1[D −C−BA]