Skip to Content
Matrix Theory

矩阵论

2019-01-24legacy assets may be incomplete

概念补充

  • 同构算子:A:VV\mathcal{A}:V\rightarrow V'
    • 全映射
    • 可逆映射
  • 同构:VVVV'是同构的线性空间
  • 数域PP上两个有限线性空间同构    \iff两空间维数相同
  • 相似ABA\sim B:基 1 下矩阵为AA,基 2 下矩阵为BB,从基 1 到基 2 的过渡矩阵为CC,则B=C1ACB=C^{-1}AC
  • 相抵ABA\simeq B: A,Bm*n阶矩阵,阶矩阵,m阶方阵阶方阵Dn阶方阵阶方阵C满足满足B=DAC$
    • 不同维线性空间VnV^nVmV^m下,同一个线性变换在不同基偶下所对应的矩阵AABB之间的关系
  • 非奇异    \iff非退化    \iff满秩    \iff可逆
  • AA有 n 个线性无关特征向量时,AA有完备的特征向量系;否则AA为亏损矩阵

内积空间上的等积变换

  • 度量矩阵(Gram 矩阵): 基的度量矩阵AA,Aij=(ei,ej)A_{ij}=(\vec e_i,\vec e_j)

    • 对称正定矩阵
    • 不同基的度量矩阵相合
  • (x,y)=XTAY(\vec x,\vec y)=X^TAY

  • xyx\perp y: (x,y)=0(x,y)=0, x 与 y 正交

  • 施密特正交规范法

    • ym=xm+km1ym1++km,m1y1\vec y_m'=\vec x_m + k_{m1}\vec y_{m-1}'+\cdots+k_{m,m-1}\vec y_1'
    • kmi=(xm,ymi)(ymi,ymi)k_{mi}=-\frac{(\vec x_m,\vec y_{m-i}')}{(\vec y_{m-i}',y_{m-i}')}
    • yi=1yiyi\vec y_i=\frac{1}{|\vec y_i'|}\vec y_i'
  • 标准正交基    \iff度量矩阵为单位矩阵

  • 复数域内积

    • (x,y)=(y,x)(x,y)=\overline{(y,x)}
    • (x+y,z)=(x,z)+(y,z)(x+y,z)=(x,z)+(y,z)
    • (kx,y)=k(x,y)(kx,y)=k(x,y)
    • (x,x)0(x,x)\geq0, (x,x)=0    x=0(x,x)=0\iff x=0
  • 酉空间:定义了内积的复数域上的线性空间

  • 转置共轭:yH=yTy^H=\overline{y}^T

    • (x,y)=yHx(\vec x,\vec y)=\vec y^H\vec x
      • 夹角:cos2x,y=(x,y)(y,x)(x,x)(y,y)cos^2\langle\vec x,\vec y\rangle=\frac{(x,y)(y,x)}{(x,x)(y,y)}
  • 柯西施瓦兹不等式:(x,y)(y,x)(x,x)(y,y)(x,y)(y,x)\leq(x,x)(y,y)

  • 单位向量: xTx=1x^Tx=1

  • 正交变换: (A(x),A(y))=(x,y)(\mathcal{A}(\vec x),\mathcal{A}(\vec y))=(\vec x,\vec y)

  • 正交变换的充要条件:

    • xV, (A(x),A(x))=(x,x)\forall \vec x\in V,\ (\mathcal{A}(\vec x),\mathcal{A}(\vec x))=(\vec x,\vec x)
    • 任一组标注正交基变换后仍是一组标准正交基
    • 任一组标准正交基下为正交矩阵 AT=A1A^T=A^{-1}
  • 正交矩阵    \iff行(列)向量为标准正交向量

  • 初等旋转变换(Givens 变换)

    • 初等旋转矩阵RijR_{ij}:

      Li=[10000100cs001000s0c0001], c2+s2=1L^i=\left[\begin{matrix} 1 & 0 &\cdots & 0 & \cdots & 0\newline 0 & 1 &\cdots & 0 & \cdots & 0\newline \vdots & \vdots & c & \vdots & s & \vdots \newline 0 & 0 & \cdots & 1 & \cdots & 0 \newline 0 & 0 & -s & 0 & c & \vdots\newline \vdots & \vdots & \ddots & \vdots & \ddots & \vdots\newline 0 & 0 & \cdots & 0 & \cdots & 1 \end{matrix}\right],\ c^2+s^2=1
      • detRij=1\det R_{ij}=1
      • RijR_{ij}正交
      • Y=RijXY=R_{ij}X, 则仅两个分量变化
        • yi=cxi+sxjy_i=cx_i+sx_j
        • yj=sxi+cxjy_j=-sx_i+cx_j
      • 化第jj个分量为 0
        • s=xjxi2+xj2s=\frac{x_j}{\sqrt{x_i^2+x_j^2}}
        • c=xixi2+xj2c=\frac{x_i}{\sqrt{x_i^2+x_j^2}}
        • yi=xi2+xj2y_i=\sqrt{x_i^2+x_j^2}
  • 镜像变换(Householder 变换)

    • 镜像变换:线性变换将向量ξ\xi映射为与单位向量ω\omega正交的 n-1 维子空间的对称像η\eta,且η=(I2ωωT)ξ=Hξ\eta=(I-2\omega\omega^T)\xi = H\xi

    • 初等反射矩阵:H=I2ωωTH=I-2\omega\omega^T

    • HH对称正交

    • QQ为初等旋转矩阵之积

      H=Q1[100011001]QH=Q^{-1} \left[\begin{matrix} -1 & 0 & \cdots & 0\newline 0 & 1 & \cdots & 1 \newline \vdots & \vdots & \ddots & \vdots\newline 0 & 0 & \cdots & 1 \end{matrix}\right] Q

埃尔米特矩阵

kronecker product

(AB)p(r1)+v,q(s1)+w=arsbvw(AB)p(r1)+v,q(s1)+w=arsbvw{\displaystyle (A\otimes B)_{p(r-1)+v,q(s-1)+w}=a_{rs}b_{vw}} {\displaystyle (A\otimes B)_{p(r-1)+v,q(s-1)+w}=a_{rs}b_{vw}}

求导

  • xTax=aTxx=a\frac{\partial x^Ta}{\partial x}=\frac{\partial a^Tx}{\partial x}=a
  • ABx=AxB+ABA\frac{\partial AB}{\partial x}=\frac{\partial A}{\partial x}B+A\frac{\partial B}{\partial A}
  • A1x=A1AxA1\frac{\partial A^{-1}}{\partial x}=-A^{-1}\frac{\partial A}{\partial x}A^{-1}
  • bTXTXcX=X(bcT+cbT)\frac{\partial b^TX^TXc}{\partial X}=X(bc^T+cb^T)
  • (Bx+b)TC(Dx+d)x=BTC(Dx+d)+DTC(Bx+b)\frac{(Bx+b)^TC(Dx+d)}{\partial x}=B^TC(Dx+d)+D^TC(Bx+b)

奇异值分解

  • A=UΣVTA=U\Sigma V^T
  • UTU=IU^TU=I, VTV=IV^TV=I

瑞利商

  • R(A,x)=xHAxxHxR(A,x)=\frac{x^HAx}{x^Hx}, AA 为 Hermitan 矩阵
    • 最大值:AA 最大的特征值
    • 最小值:AA 最小的特征值
  • 广义瑞利商:R(A,B,x)=xHAxxHBx=R(B12AB12,B12x)R(A,B,x)=\frac{x^HAx}{x^HBx}=R(B^{-\frac{1}{2}}AB^{-\frac{1}{2}},B^{-\frac{1}{2}}x)

LU 分解

  • 消元矩阵: LiL^{i} 为高斯消元第 i 次对应的矩阵,由倍加阵相乘得到
Li=[10000100001000a(i+1)iiaii00aniiaiii1]L^i=\left[\begin{matrix} 1 & 0 &\cdots & 0 & \cdots & 0\newline 0 & 1 &\cdots & 0 & \cdots & 0\newline \vdots & \vdots &\ddots & \vdots & \ddots & \vdots \newline 0 & 0 & \cdots & 1 & \cdots & 0 \newline 0 & 0 & \cdots & -\frac{a^i_{(i+1)i}}{a^{ii}} & \cdots & \vdots\newline \vdots & \vdots & \ddots & \vdots & \ddots & \vdots\newline 0 & 0 & \cdots & -\frac{a^i_{ni}}{a^i_{ii}} & \cdots & 1 \end{matrix}\right]
  • 高斯消元过程: 上三角矩阵U=An=Ln1L2L1A1U=A^n=L^{n-1}\cdots L^{2}L^{1}A^{1}
  • 高斯消元过程可以进行到底    \iff AA的前n1n-1个顺序主子式都不为 0
  • L=(L1)1(Ln1)1L=(L^1)^{-1}\cdots(L^{n-1})^{-1} 为单位下三角矩阵,Lij=aijjajjj,i<jL_{ij}=-\frac{a^j_{ij}}{a^j_{jj}},i<j
  • 三角分解(LU 分解): 方阵A=LUA=LU
  • 杜立特分解(Doolittle Decomposition): LL为单位下三角矩阵
  • 克劳特分解(Crout Decomposition): LL为单位上三角矩阵
  • LDULDU基本定理:n 阶方阵AA可唯一分解为A=LDU    A=LDU\iff AA的前n1n-1个顺序主子式Ak=A_k=\not0, L,UL,U为单位下、上三角矩阵,DD为对角矩阵diag(A1A0,A2A1,...,AnAn1),A0=1diag(\frac{A_1}{A_0},\frac{A_2}{A_1},...,\frac{A_n}{A_{n-1}}),A_0=1
    • 杜立特和克劳特分解存在则唯一
  • LDU 分解
    • 计算L1,,Ln,A1,,AnL^1,\cdots,L^n,A^1,\cdots,A^n
    • L=(L1)1(Ln1)1L=(L^1)^{-1}\cdots(L^{n-1})^{-1}
    • An=DUA^n=DU
    • A=LAn=LDUA=LA^n=LDU
  • 克劳德分解
    • lij=aijk=1j1likukjl_{ij} = a_{ij}-\sum_{k=1}^{j-1}l_{ik}u_{kj}
    • uij=(aijk=1i1likukj)/liiu_{ij} = (a_{ij}-\sum_{k=1}^{i-1}l_{ik}u_{kj})/l_{ii}
    • 从上至下,从左至右
  • 杜立特分解
    • uij=aijk=1i1likukju_{ij}=a_{ij}-\sum_{k=1}^{i-1}l_{ik}u_{kj}
    • lij=(aijk=1j1likukj)/ujjl_{ij}=(a_{ij}-\sum_{k=1}^{j-1}l_{ik}u_{kj})/u_{jj}
  • n 阶对称正定矩阵AA,唯一存在一个实的非奇异的下三角矩阵LLA=LLTA=LL^T,称为乔累斯基分解(平方根分解)
  • 乔累斯基分解
    • lij=(aijk=1j1)/ljjl_{ij}=(a_{ij}-\sum_{k=1}^{j-1})/l_{jj}
  • n 阶对称正定矩阵AA可唯一分解为A=LDLTA=LDL^T,其中,对角矩阵DD元素为LL对角线元素倒数
  • 不带平方根的乔累斯基分解
    • lij=aijk=1j1likljklkkl_{ij}=a_{ij}-\sum_{k=1}^{j-1}\frac{l_{ik}l_{jk}}{l_{kk}}
  • 置换矩阵P=(eπ(1),,eπ(n))TP=(e_{\pi(1)},\cdots,e_{\pi(n)})^T
  • 任意矩阵AA可进行PLUPLU分解,A=PLUA=PLU,其中PP为置换矩阵

QR 分解

  • QRQR分解(正交三角分解):实非奇异矩阵AA能化成正交矩阵QQ和实非奇异上三角矩阵RR的乘积,A=QRA=QR
  • 任意非奇异 n 阶矩阵,存在三角分解,A=QRA=QR且相差一个对角线元素之差绝对值全等于 1 的对角矩阵因子
    • 利用施密特正规化正交向量矩阵AA和对应的上三角矩阵BB
  • Givens 方法=-9
  • Householder 方法

最大秩分解

奇异值分解

谱分解