• 约 6 分钟

完了😭,期末了来不及学矩阵分析了,算了就记住这些吧

根据B站知名UP主苑长,章节的重要性排序为:

1>3≈4>6>2≈51>3\approx4>6>2\approx5

第1章 线性空间与线性变换

前两章的概念和定理较多,着重理解记忆。

  • 定义
    • 封闭性
    • 交结零负,结单分分(😅)
  • 维数:可以理解为有多少影响因子
  • 基与基到基的过渡矩阵
    • 过渡矩阵
    • 一个定理
      • 如果 α\alpha 在基 α1,⋯ ,αn\alpha_1,\cdots,\alpha_n 下坐标为 (x1,⋯ ,xn)T(x_1,\cdots,x_n)^T ,在基 β1,⋯ ,βn\beta_1,\cdots,\beta_n 下坐标为 (y1,⋯ ,yn)T(y_1,\cdots,y_n)^T ,从 α1,⋯ ,αn\alpha_1,\cdots,\alpha_n 到 β1,⋯ ,βn\beta_1,\cdots,\beta_n 的过渡矩阵为 AA ,则(x1,⋯ ,xn)T=A(y1,⋯ ,yn)T (x_1,\cdots,x_n)^T=A(y_1,\cdots,y_n)^T
      • 且 AA 是可逆的
      • 务必注意顺序啊!
  • 子空间及其和、交、直和
    • 注意和的定义
    • 和的性质:交结
    • ⭐维数公式
      • dim(V1+V2)=dim(V1)+dim(V2)−dim(V1∩V2)dim(V_1+V_2)=dim(V_1)+dim(V_2)-dim(V_1\cap V_2)
  • 生成子空间
    • ⭐若 V1=Span(α1,⋯ ,αs),V2=Span(β1,⋯ ,βt)V_1=Span(\alpha_1,\cdots,\alpha_s),V_2=Span(\beta_1,\cdots,\beta_t) 则 V1+V2=Span(α1,⋯ ,αs,β1,⋯ ,βt)V_1+V_2=Span(\alpha_1,\cdots,\alpha_s,\beta_1,\cdots,\beta_t)
  • 线性变换及其在某一基下的矩阵
    • 在某一基下的矩阵相当于 α1,⋯ ,αn\alpha_1,\cdots,\alpha_n 到 T(α1),⋯ ,T(αn)T(\alpha_1),\cdots,T(\alpha_n) 的过渡矩阵
    • 线性变换的核与像
    • 一个定理
      • 如果 α\alpha 在基 α1,⋯ ,αn\alpha_1,\cdots,\alpha_n 下坐标为 (x1,⋯ ,xn)T(x_1,\cdots,x_n)^T , TT 为在这一基下的矩阵,则 T(α)T(\alpha) 在这一基下的坐标 (y1,⋯ ,yn)T=A(x1,⋯ ,xn)T(y_1,\cdots,y_n)^T=A(x_1,\cdots,x_n)^T
      • 注意与前面的那个定理区分
    • 线性变换是同一线性空间内的映射,因此其矩阵必须是方阵
    • ⭐线性变换在不同基下矩阵的关系
      • 线性变换 TT 在基 α1,⋯ ,αn\alpha_1,\cdots,\alpha_n 和基 β1,⋯ ,βn\beta_1,\cdots,\beta_n 下矩阵分别为 A,BA,B ,则 A∼BA\sim B ,若 α1,⋯ ,αn\alpha_1,\cdots,\alpha_n 到 β1,⋯ ,βn\beta_1,\cdots,\beta_n 过渡矩阵为 PP ,则 B=P−1APB=P^{-1}AP .

第2章 内积与欧氏空间

最小二乘解求法

ATAx=ATbA^TAx=A^Tb

解出x

求奇异值

  1. 求 AHAA^HA
  2. 求 AHAA^HA 的特征值
  3. 非零特征值开方就是奇异值了

第3章 矩阵的标准形

求Jordan标准形的方法

方法一

  1. λE−A\lambda E-A
  2. 行列式因子
  3. 不变因子
  4. 初级因子
  5. Jordan块
  6. J

方法二(据说更简单)

  1. λE−A\lambda E-A
  2. 化为等价得Smith标准形
  3. 不变因子->初级因子->…

Jordan标准形的1可以在左下也可以在右上

已知初级因子求不变因子

不同的根的最高次幂放在 dr(λ)d_r(\lambda)

⭐⭐⭐求 λ\lambda -矩阵Smith标准形

貌似只可意会不可言传,还是通过例题体会吧

求最小多项式方法

  1. 根据哈密顿-凯莱(Cayley-Hamilton)定理 f(λ)=∣λE−A∣=0f(\lambda)=|\lambda E-A|=0 是一个特殊的零化多项式
  2. 观察某些幂能不能低次

Hints

  1. Jordan标准形与块序无关
  2. 多项式中次数为非负整数
  3. 可以运用零化多项式简化矩阵运算
  4. λ\lambda -矩阵的秩:不为零子式的最大阶数
  5. 满秩 λ\lambda -矩阵不一定可逆,可逆的充要条件是 ∣A(λ)∣=C≠0|A(\lambda)|=C\not=0
  6. λ\lambda -矩阵的初等变换
    1. 互换
    2. 数乘非零k
    3. 乘多项式加到另一行

第4章 矩阵分解

LU分解步骤

  1. 写成分块矩阵形式 (A,E)(A,E)
  2. 只做第三类初等行变换,变为如下形式
[上三角∣主元为1的下三角]\begin{bmatrix}上三角&|&主元为1的下三角\end{bmatrix}
  1. 左侧记为 UU ,右侧记为 PP
  2. 求 PP 的逆 L=1∣P∣P∗L=\frac{1}{|P|}P^*
  3. A=LUA=LU

求Ax=b可先分解A=LU

LDU分解步骤

  1. LU分解
  2. 老U中的主对角元组成D,新U为分别除以主角元

⭐⭐QR分解步骤

  1. AA 写作 (α1,α2,⋯ ,αn)(\alpha_1,\alpha_2,\cdots,\alpha_n)
  2. 施密特 得 β1⋯\beta_1\cdots
  3. 单位化 得 γ1⋯\gamma_1\cdots
  4. 得
Q=(γ1,γ2,⋯ ,γn)Q=(\gamma_1,\gamma_2,\cdots,\gamma_n) R=[∣∣β1∣∣(α2,γ1)⋯(αn,γ1)∣∣β2∣∣⋯(αn,γ2)⋱⋮∣∣βn∣∣]R=\begin{bmatrix}||\beta_1||&(\alpha_2,\gamma_1)&\cdots&(\alpha_n,\gamma_1)\\\\&||\beta_2||&\cdots&(\alpha_n,\gamma_2)\\\\&&\ddots&\vdots\\\\&&&||\beta_n||\end{bmatrix}

满秩分解步骤

  1. 初等行变换为行最简形
  2. 每个首1在原矩阵中的列组成 PP
  3. 首1行组成 QQ
  4. A=PQA=PQ

奇异值分解步骤(这个简单,只要到时候不做就行了)

  1. 求出 AHAA^HA 的特征值 λ1⩾⋯⩾λr>0=λr+1⩾⋯⩾λn\lambda_1\geqslant\cdots\geqslant\lambda_r>0=\lambda_{r+1}\geqslant\cdots\geqslant\lambda_n ,确定非零奇异值 di=λid_i=\sqrt{\lambda_i}
  2. 求特征向量并单位正交化,得 γ1,⋯ ,γr,γr+1,⋯ ,γn\gamma_1,\cdots,\gamma_r,\gamma_{r+1},\cdots,\gamma_n 令
V=(V1,V2)V=(V_1,V_2) V1=(γ1,⋯ ,γr), V2=(γr+1,⋯ ,γn)V_1=(\gamma_1,\cdots,\gamma_r),\ V_2=(\gamma_{r+1},\cdots,\gamma_n)
  1. 计算 U1=AV1D−1U_1=AV_1D^{-1} ,其中 D=diag(d1,⋯ ,dr)D=diag(d_1,\cdots,d_r) ,取 U2∈Cm×(m−r)U_2\in C^{m\times(m-r)} 使得其 m−rm-r 个列向量和 U1U_1 的 rr 个列向量构成 CmC^m 的标准正交基
  2. 令 U=(U1,U2)U=(U_1,U_2) ,从而 A=UΣVHA=U\Sigma V^H ,其中
Σ=[DOOO]\Sigma=\begin{bmatrix}D&O\\\\O&O\end{bmatrix}

求加号逆的步骤

  1. 若行满秩则 A+=AH(AAH)−1A^+=A^H(AA^H)^{-1}
  2. 若列满秩则 A+=(AHA)−1AHA^+=(A^HA)^{-1}A^H
  3. 若非满秩则进行满秩分解 A=BCA=BC ,其中 BB 为列满秩, CC 为行满秩
  4. 求 B+B^+ 和 C+C^+
  5. A+=C+B+A^+=C^+B^+

第5章 范数理论及其应用

向量的范数计算方法

∣∣α∣∣1=∑i=1n∣xi∣||\mathbf\alpha||\substack 1=\sum_{i=1}^n|x_i| ∣∣α∣∣2=∑i=1n∣xi∣2||\mathbf\alpha||\substack 2=\sqrt{\sum_{i=1}^n|x_i|^2} ∣∣α∣∣∞=max⁡i∣xi∣||\mathbf\alpha||\substack \infty=\max_i|x_i|

都是要先算绝对值的

矩阵的范数计算方法

m-范数

∣∣A∣∣m1=∑i=1n∑j=1n∣aij∣||\mathbf{A}||\substack {m_1}=\sum_{i=1}^n\sum_{j=1}^n|a_{ij}| ∣∣A∣∣m2=∣∣A∣∣F=∑i=1n∑j=1n∣aij∣2=tr(AHA)||\mathbf A||\substack {m_2}=||\mathbf A||\substack F=\sqrt{\sum_{i=1}^n\sum_{j=1}^n|a_{ij}|^2}=\sqrt{tr(\mathbf A^H\mathbf A)} ∣∣A∣∣m∞=nmax⁡i,j∣aij∣||\mathbf A||\substack {m_\infty}=n\max_{i,j}|a_{ij}|

其他

∣∣A∣∣1=max⁡1⩽j⩽n∑i=1n∣aij∣||\mathbf A||\substack 1=\max_{1\leqslant j\leqslant n}\sum_{i=1}^n|a_{ij}|

即列模和最大者

∣∣A∣∣∞=max⁡1⩽i⩽n∑j=1n∣aij∣||\mathbf A||\substack \infty=\max_{1\leqslant i\leqslant n}\sum_{j=1}^n|a_{ij}|

即行模和最大者

∣∣A∣∣2=λn=max⁡iλi(AHA)||\mathbf A||\substack 2=\sqrt{\lambda_n}=\sqrt{\max_i\lambda_i(\mathbf A^H\mathbf A)}
  • λn\lambda_n 是矩阵 AHA\mathbf A^H\mathbf A 的最大特征值*

求极小范数最小二乘解

x0=A+b\mathbf x_0=\mathbf A^+\mathbf b

记住两个不等式

对于任意一种范数 ∣∣⋅∣∣||\mathbf\cdot|| ,有

ρ(A)⩽∣∣A∣∣\rho(\mathbf A)\leqslant||\mathbf A||

对 ∀ε>0\forall\varepsilon>0 ,存在 ∣∣⋅∣∣m||\cdot||\substack m 使得

∣∣A∣∣m⩽ρ(A)+ε||\mathbf A||\substack m\leqslant\rho(\mathbf A)+\varepsilon

其中 ρ(A)=max⁡j∣λj∣\rho(\mathbf A)=\max_j|\lambda_j| 为矩阵 A\mathbf A 的谱半径

第6章 矩阵分析及其应用

据说这章有30分

函数矩阵的导数&积分

就是每个元素求导或者求积分就完事儿了

矩阵的幂级数(收敛性判断)

  • 知识点
    • 矩阵的幂级数
      • 设 A=(aij)∈Cn×nA=(a_{ij})\in C^{n\times n} ,称形如 ∑k=0+∞akAk=a0E+a1A+⋯+akAk+⋯ \sum_{k=0}^{+\infty}a_kA^k=a_0E+a_1A+\cdots+a_kA^k+\cdots
      • 的矩阵级数为矩阵幂级数
    • 收敛性
      • 设幂级数 ∑k=0+∞akXk\sum_{k=0}^{+\infty}a_kX^k 的收敛半径为 RR , A∈Cn×nA\in C^{n\times n} ,则
        1. 当 ρ(A)<R\rho(A)<R 时, ∑k=0+∞akAk\sum_{k=0}^{+\infty}a_kA^k 绝对收敛
        2. 当 ρ(A)>R\rho(A)>R 时, ∑k=0+∞akAk\sum_{k=0}^{+\infty}a_kA^k 发散
    • 收敛半径 R=1/ρR=1/\rho
      • 其中 ρ=lim⁡k→∞∣ak+1ak∣\rho=\lim_{k\to\infty}|\frac{a_{k+1}}{a_k}|
  • 判断敛散性步骤
    1. 求收敛半径 RR
    2. 观察 ∣∣⋅∣∣1,∣∣⋅∣∣+∞||\cdot||\substack1,||\cdot||\substack{+\infty} 是否小于 RR ,若小于,则收敛
    3. 求特征值,然后 ρ(A)=max⁡jλj\rho(A)=\max_j\lambda_j
    4. 若 ρ(A)<R\rho(A)<R 则收敛
    5. 如果不好搞可以用定义判断: lim⁡k→+∞Ak=O\lim_{k\to+\infty}\mathbf A^k=\mathbf O
  • 求矩阵幂级数的和
    1. 记和为 SS
    2. 左乘 AA
    3. 错位相减

求矩阵函数步骤

  1. 求[[#求最小多项式方法|最小多项式]] φ(λ)\varphi(\lambda)
  2. 设多项式 r(λ)=a0+a1λ+⋯+arλrr(\lambda)=a_0+a_1\lambda+\cdots+a_r\lambda^r ,其次数比 φ(λ)\varphi(\lambda) 少1
  3. 设 f(λ)f(\lambda) 为原函数
  4. 将各特征值同时代入 r(λ)r(\lambda) 和 f(λ)f(\lambda) 中,如果约束不够,可以两个多项式分别对 λ\lambda 求导再代入
  5. 解得 a0,⋯ ,ara_0,\cdots,a_r
  6. f(A)=a0E+a1A+⋯+arArf(A)=a_0E+a_1A+\cdots+a_rA^r

微分方程组求解

齐次线性微分方程组

{dxdt=Axx(t0)=x0\begin{cases} \frac{dx}{dt}=Ax\\\\x(t_0)=x_0 \end{cases}

的解为

x(t)=eA(t−t0)x0x(t)=e^{A(t-t_0)}x_0

非齐次线性微分方程组

{dxdt=Ax+f(t)x(t0)=x0\begin{cases} \frac{dx}{dt}=Ax+f(t)\\\\ x(t_0)=x_0 \end{cases}

的解为

x(t)=eA(t−t0)x0+∫t0teA(t−s)f(s)ds=eA(t−t0)x0+eAt∫t0te−Asf(s)ds\begin{aligned} x(t)&=e^{A(t-t_0)}x_0+\int_{t_0}^te^{A(t-s)}f(s)ds \\\\ &=e^{A(t-t_0)}x_0+e^{At}\int_{t_0}^te^{-As}f(s)ds \end{aligned}

tt 是参数, t0t_0 是个数, x0,f(⋅)x_0,f(\cdot) 是向量

已知 eAte^{At} 求 AA

A=deAtdt∣t=0A=\frac{de^{At}}{dt}\bigg|\substack{t=0}
林威
林威 咖味十足的软件工程师