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

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

1>34>6>251>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(V1V2)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 ,则 ABA\sim B ,若 α1,,αn\alpha_1,\cdots,\alpha_n β1,,βn\beta_1,\cdots,\beta_n 过渡矩阵PP ,则 B=P1APB=P^{-1}AP .

第2章 内积与欧氏空间

最小二乘解求法

ATAx=ATbA^TAx=A^Tb

解出x

求奇异值

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

第3章 矩阵的标准形

求Jordan标准形的方法

方法一

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

方法二(据说更简单)

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

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

已知初级因子求不变因子

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

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

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

求最小多项式方法

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

Hints

  1. Jordan标准形与块序无关
  2. 多项式中次数为非负整数
  3. 可以运用零化多项式简化矩阵运算
  4. λ\lambda -矩阵的:不为零子式的最大阶数
  5. 满秩 λ\lambda -矩阵不一定可逆,可逆的充要条件是 A(λ)=C0|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=1PPL=\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
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=AV1D1U_1=AV_1D^{-1} ,其中 D=diag(d1,,dr)D=diag(d_1,\cdots,d_r) ,取 U2Cm×(mr)U_2\in C^{m\times(m-r)} 使得其 mrm-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=1nxi||\mathbf\alpha||\substack 1=\sum_{i=1}^n|x_i| α2=i=1nxi2||\mathbf\alpha||\substack 2=\sqrt{\sum_{i=1}^n|x_i|^2} α=maxixi||\mathbf\alpha||\substack \infty=\max_i|x_i|

都是要先算绝对值的

矩阵的范数计算方法

m-范数

Am1=i=1nj=1naij||\mathbf{A}||\substack {m_1}=\sum_{i=1}^n\sum_{j=1}^n|a_{ij}| Am2=AF=i=1nj=1naij2=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)} Am=nmaxi,jaij||\mathbf A||\substack {m_\infty}=n\max_{i,j}|a_{ij}|

其他

A1=max1jni=1naij||\mathbf A||\substack 1=\max_{1\leqslant j\leqslant n}\sum_{i=1}^n|a_{ij}|

即列模和最大者

A=max1inj=1naij||\mathbf A||\substack \infty=\max_{1\leqslant i\leqslant n}\sum_{j=1}^n|a_{ij}|

即行模和最大者

A2=λn=maxiλ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 使得

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

其中 ρ(A)=maxjλ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 ACn×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
      • 其中 ρ=limkak+1ak\rho=\lim_{k\to\infty}|\frac{a_{k+1}}{a_k}|
  • 判断敛散性步骤
    1. 求收敛半径 RR
    2. 观察 1,+||\cdot||\substack1,||\cdot||\substack{+\infty} 是否小于 RR ,若小于,则收敛
    3. 求特征值,然后 ρ(A)=maxjλj\rho(A)=\max_j\lambda_j
    4. ρ(A)<R\rho(A)<R 则收敛
    5. 如果不好搞可以用定义判断: limk+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(tt0)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(tt0)x0+t0teA(ts)f(s)ds=eA(tt0)x0+eAtt0teAsf(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=deAtdtt=0A=\frac{de^{At}}{dt}\bigg|\substack{t=0}