抱歉,您的浏览器无法访问本站
本页面需要浏览器支持(启用)JavaScript
了解详情 >

CodingStudio

努力进步

引言

  • 本文为硕士课程中矩阵理论的课堂笔记
  • 能力有限,不喜勿喷

矩阵理论

  • 重根的特征向量若不正交,则使用许米正交化方法正交
  • 秩1阵的分解应用
    秩1阵的根
    例题

预备知识

  • AH=AA^H=A, A为Hermite阵
    • AH=AA^H=-A, A为斜Hermite阵
  • Hermite阵的对角线元素为实数
    • 实对称阵是Hermite矩阵
    • 若B为斜Hermite阵,则iB,B/i为Hermite阵特征根都为虚数或者0(利用B/i为Hermite阵证明)
    • 若A为Hermite阵,则iA为斜Hermite阵
  • Herimte分解定理:若A为Hermite阵,则存在优阵Q使得A=QDQ1=Q(λ100λn)Q1,λ1,...λn为实数A =\mathbf{Q} \mathbf{D} \mathbf{Q}^{-1}=\mathbf{Q} \left(\begin{array}{ccc}\lambda_{1} & & 0 \\ & \ddots & \\ 0 & & \lambda_{n}\end{array}\right) \mathbf{Q}^{-1},且\lambda_{1} ,...\lambda_{n}为实数
    • Hermite阵A必有n个正交特向
    • Hermite阵特征根都是实数
    • AHA=0A=0A^HA=0 \leftrightarrow A=0(利用秩进行证明,或利用迹公式证明)

迹公式

  • 矩阵模平方(F范数):tr(AHA)=tr(AAH)=ai,j2\operatorname{tr}\left(A^{H} A\right)=\operatorname{tr}\left(A A^{H}\right)=\sum\left|a_{i, j}\right|^{2}
  • tr(ABT)=tr(BTA)=ai,jbi,j\operatorname{tr}\left(A B^{T}\right)=\operatorname{tr}\left(B^{T} A\right)=\sum a_{i, j} b_{i, j}
  • tr(ABH)=tr(BHA)=aijbij\operatorname{tr}\left(A B^{H}\right)=\operatorname{tr}\left(B^{H} A\right)=\sum a_{i j} \overline{b_{i j}}

内积公式

  • (X,Y)=YHX=tr(XYH)=x1y1++xnyn,(X, Y)=Y^{H} X=\operatorname{tr}\left(X Y^{H}\right)=x_{1} \overline{y_{1}}+\cdots+x_{n} \overline{y_{n}}, \quad for X,YCnX, Y \in \mathrm{C}^{\mathrm{n}}
  • 模长(X,X)=tr(XXH)=XHX=x1x1++xnxn=x12++xn2=X2(X, X)=\operatorname{tr}\left(X X^{H}\right)=X^{H} X=x_{1} \overline{x_{1}}+\cdots+x_{n} \overline{x_{n}}=\left|x_{1}\right|^{2}+\cdots+\left|x_{n}\right|^{2}=|X|^{2}
  • 复空间矩阵内积:
  • (A,B)=YHX=tr(ABH)=a1b1++anbn,(A, B)=Y^{H} X=\operatorname{tr}\left(A B^{H}\right)=a_{1} \overline{b_{1}}+\cdots+a_{n} \overline{b_{n}}, \quad for A,BCmnA, B \in \mathrm{C}^{\mathrm{m*n}}
  • 矩阵模长公式:A2=tr(AAH)=tr(AHA)=aij2\|A\|^{2}=\operatorname{tr}\left(A A^{H}\right)=\operatorname{tr}\left(A^{H} A\right)=\sum\left|a_{i j}\right|^{2}
  • 向量垂直:XY(Y,X)=0(X,Y)=0X \perp Y \Leftrightarrow(Y, X)=0 \Leftrightarrow(X, Y)=0
    • 性质:1.满足勾股定理;2.可推广至n个向量互垂直

优阵_预优阵

  • 预优阵:方阵各列互垂直
    • 预半优阵:矩阵各列互垂直
    • 引理:A为预优阵(预半优阵)\Leftrightarrow当且仅当AHAA^{H} A为对角形
  • 优阵:方阵各列互垂直,且各列的模长为1
    • 半优阵:矩阵各列互垂直,且各列的模长为1
    • 引理:A为优阵(半优阵)\Leftrightarrow当且仅当AHA=IA^{H} A = I(AAH=IA A^{H} = I),则AH=A1A^{H} = A^{-1}
    • **性质:A为优阵(半优阵)
        1. AX2=X2|AX|^{2} = |X|^{2};保模长
        1. xyAxAyx \perp y \Rightarrow A x \perp A y;保正交
        1. (Ax,Ay)=(x,y)(A x, A y)=(x, y)

许尔公式

  • 补充公式:设n阶可逆阵P可以按列写成P=(α1,α2,,αn)\mathrm{P}=\left(\alpha_{1}, \alpha_{2}, \cdots, \alpha_{n}\right),则P1α1=e1,P1α2=e2,,P1αn=en\mathrm{P}^{-1} \alpha_{1}=e_{1}, \mathrm{P}^{-1} \alpha_{2}=e_{2}, \cdots, \mathrm{P}^{-1} \alpha_{n}=e_{\mathrm{n}}
  • 二阶许尔公式的引理:
    • 二阶许尔公式
  • 许尔公式:每个方阵A,存在可逆阵P使得P1AP=D=(λ10λn)\mathbf{P}^{-1} \mathbf{A P}=\mathbf{D}=\left(\begin{array}{ccc}\lambda_{1} & \ldots & * \\ & \ddots & \vdots \\ 0 & & \lambda_{n}\end{array}\right)为上三角
  • 优相似三角化(许尔公式2):每个方阵A,存在优阵Q,使Q1AQ=D=(λ10λn)\mathbf{Q}^{-1} \mathbf{A Q}=\mathbf{D}=\left(\begin{array}{ccc}\lambda_{1} & \ldots & * \\ & \ddots & \vdots \\ 0 & & \lambda_{n}\end{array}\right)为上三角
    • 推论:每个方阵都优相似与上三角阵

换位公式

  • A=An×p\mathbf{A}=\mathbf{A}_{\mathrm{n} \times p}B=Bp×n\mathbf{B}=\mathbf{B}_{\mathrm{p} \times n},且np\mathrm{n} \geq \mathrm{p},则λInAB=λnpλIpBA\left|\lambda \mathbf{I}_{\mathrm{n}}-\mathbf{A B}\right|=\lambda^{n-p}\left|\lambda \mathbf{I}_{p}-\mathbf{B} \mathbf{A}\right|
    • AB为n阶方阵,BA为p阶方阵,AB与BA的非0根相同,只相差n-p个0根
      • tr(AB)=tr(BA)\operatorname{tr}\left(A B\right)=\operatorname{tr}\left(B A\right)
      • 行列式降阶公式:In+AB=Ip+BA\left|\mathbf{I}_{\mathrm{n}}+\mathbf{A B}\right|=\left|\mathbf{I}_{\mathrm{p}}+\mathbf{B} \mathbf{A}\right|
    • 若n=p,则AB与BA的特征根相同

秩1阵

  • 方阵的秩 rank(A)=1 ,则A叫做秩1阵(比例阵)
    • 特征根为λ(A)={tr(A),0,,0}\lambda(\mathbf{A})=\{\operatorname{tr}(\mathbf{A}), 0, \cdots, 0\}
    • λ1=tr(A)\lambda_1 = \operatorname{tr}(\mathbf{A}),则Aα=λ1αA\alpha = \lambda_1 \alpha,且A2=λ1AA^2 = \lambda_1 A
    • 若方阵满足A2=λ1AA^2 = \lambda_1 A,则A中各列都是特征向量
    • 秩1阵必有分解:A=αβ=αn×1β1×n\mathbf{A}=\boldsymbol{\alpha} \boldsymbol{\beta}=\boldsymbol{\alpha}_{\mathrm{n} \times 1} \boldsymbol{\beta}_{1 \times \mathrm{n}},且Aα=λ1αA\alpha = \lambda_1 \alphaλ1=tr(A)\lambda_1 = \operatorname{tr}(\mathbf{A})
      秩1阵证明

遗传公式

  • 根遗传公式
    • 满足平移公式、倍公式、幂公式、逆公式
  • 向量遗传定理:A的特征向量为X1,X2,...,XnX_1,X_2,...,X_n也是f(A)的特征向量
    • 使用平移法求解矩阵的特征值和特征向量

补充定理:若A=A(np)A = A_(n*p),B=B(pn)B = B_(p*n),且n>p,则AB=0|A B | = 0
证明

秩公式

  • 定理1: 任意复矩阵A,AHAA^HA,AAHAA^H为Hermite阵,且半正定
    • 证明:显然AHAA^H A为Hermite阵,利用二次型XH(AHA)X=(AX)HAX=AX20X^H (A^H A) X = (AX)^H AX = | AX |^2 \geq 0,则AHAA^H A为半正定的
  • 定理2: 任意矩阵A,则(AHA)X=0(A^H A)X = 0AX=0AX=0具有相同解
    • 证明:若(AHA)X=0(A^HA)X = 0成立,则AX2=(AX)HAX=XH(AHA)X=0|AX|^2=(AX)^H AX=X^H (A^H A)X=0,AX=0AX=0\rightarrow |AX|=0 \rightarrow AX=0
  • 定理3:任意m*n矩阵,r(AHA)=r(A)r(A^HA)=r(A)
    • 证明:由定理2可知,(AHA)X=0(A^H A)X = 0AX=0AX=0同解,则r(AHA)=r(A)r(A^HA)=r(A)
  • 秩的第二定义:m*n阶矩阵A,r(A)=r的等价条件使A中必有一个r阶子式非0,所有r+1阶子式都为0
    • r(AHA)=r(A),r(AHA)=r(AH),r(A)=r(AH)=r(Aˉ)r(A^HA) = r(A),r(A^HA)=r(A^H),r(A)=r(A^H)=\operatorname{r}(\bar{A})

正规阵

  • 若方阵A满足AHA=AAHA^HA=AA^H,则A叫正规阵
    • 对角阵必正规
    • Hermite阵与斜Hermite阵必正规
      • 若A正规,则A与AHA^H必有相同特征向量,AX=cXAHX=cˉXAX=cX\leftrightarrow A^HX = \bar{c}X
      • 若A的特根为λ(A)={λ1,...λn}\lambda(A)=\{\lambda_1,...\lambda_n\},则λ(AH)={λ1ˉ,...,λnˉ}\lambda(A^H)=\{\bar{\lambda_1},...,\bar{\lambda_n}\}
    • 实对称阵与实反对称阵都正规
    • 优阵必正规
    • 已知正规阵,可知正规阵的方法
      • 正规阵的倍数依然正规
      • 正规阵平移后依然正规
  • 优相似定理:正规阵的优相似必正规
  • 多项正规定理:正规阵A的多项式f(A)也正规
  • 三角正规定理:三角正规阵一定是对角形(同时适用于分块矩阵)
    • 严格三角阵不是正规阵

正规分解定理

  • 若方阵A正规,则存在优阵Q(QH=Q1Q^H=Q^{-1}),使Q1AQ=D=(λ100λn)(对角形),其中λ(A)={λ1,...λn}\mathbf{Q}^{-1} \mathbf{A Q}=\mathbf{D}=\left(\begin{array}{ccc}\lambda_{1} & & 0 \\ & \ddots & \\ 0 & & \lambda_{n}\end{array}\right)(对角形),其中\lambda(A) = \{\lambda_1,...\lambda_n\},Q1AQ=D\mathbf{Q}^{-1} \mathbf{A Q}=\mathbf{D}可知,Q中的列都是特向量
    证明

    • 推论:正规阵A中恰有n个正交特向,不同特征根对应的特征向量正交(利用内积证明)
    • 正规阵一定优相似与对角形
    • 若A优相似与对角形,则A必正规
  • 分解方法:

    1. 求特征根
    2. 求正交特征向量X1X2...XnX_1 \perp X_2 \perp ...X_n
    3. Q=(X1X1,,XnXn)Q=\left(\frac{X_{1}}{\left|X_{1}\right|}, \cdots, \frac{X_{n}}{\left|X_{n}\right|}\right),D=(λ100λn)\mathbf{D}=\left(\begin{array}{ccc}\lambda_{1} & & 0 \\ & \ddots & \\ 0 & & \lambda_{n}\end{array}\right)

结论:A=(0110),λ(A)={i,i},X1=(1,i)T,X2=(i,1)TA=\left(\begin{array}{cc}0 & -1 \\ 1 & 0\end{array}\right),\lambda(A) = \{-i,i\},X_1=(1,i)^T,X_2=(i,1)^T为正规阵

QR分解

  • QR公式:若A为列无关(高阵,列满秩,r(A)=列数pr(A)=列数p),则A=QR,Q=Qnp为列U(QHQ=Ip)R=(r10rn)为上三角,对角元为正数A=QR,Q=Q_{n*p}为列U阵(Q^HQ=I_p),\mathbf{R}=\left(\begin{array}{ccc}r_{1} & \ldots & * \\ & \ddots & \vdots \\ 0 & & r_{n}\end{array}\right)为上三角,对角元为正数,且有公式R=QHAR=Q^HA
    • 若A为可逆方阵,则有A=QR,Q=Qnn为优阵A=QR,Q=Q_{n*n}为优阵
  • 许米正交公式:
    • 设3向量α1,α2,α3\alpha_1,\alpha_2,\alpha_3线性无关,则令β1=α1,β2=α2(α2,α1)α12α1,β3=α3(α3,α1)α12α1(α3,β2)β22β2\beta_{1}=\alpha_{1},\beta_{2}=\alpha_{2}-\frac{\left(\alpha_{2}, \alpha_{1}\right)}{\left|\alpha_{1}\right|^{2}} \alpha_{1},\beta_{3}=\alpha_{3}-\frac{\left(\alpha_{3}, \alpha_{1}\right)}{\left|\alpha_{1}\right|^{2}} \alpha_{1}-\frac{\left(\alpha_{3}, \beta_{2}\right)}{\left|\beta_{2}\right|^{2}} \beta_{2},则β1,β2,β3\beta_1,\beta_2,\beta_3互正交
  • 分解方法
    1. 使用许米特公式求解优阵(或半优阵)Q
    2. 使用R=QHAR=Q^HA求三角阵R
    3. 写出分解

满秩分解(高低分解)

  • A=AmnA=A_{m*n}的秩r(A)=r>0r(A)=r>0,则有分解A=BCA=BC,其中B为列满秩(高阵),C为行满秩(低阵)
  • 分解方法:
    • 将A进行行变换化简为简化阶形式,A 行变换 (Ir00),r=r(A)A \stackrel{\text { 行变换 }}{\longrightarrow}\left(\begin{array}{cc}I_{r} & * \\ \cdots & \cdots \\ 0 & 0\end{array}\right), r=r(A),从A中去除前r列即为B,C=(Ir,)C=(I_r,*)
    • 秩1分解法:若r(A)=1(各列成比例)r(A)=1(各列成比例),则A=(a1an)(b1bn)=αβA=\left(\begin{array}{c}a_{1} \\ \vdots \\ a_{n}\end{array}\right)\left(\begin{array}{lll}b_{1} & \cdots & b_{n}\end{array}\right)=\alpha \beta

乔利斯分解与平方根公式

  • 乔利斯定理:若A>0(正定),则A=RHRA=R^HR,其中R=(b10bn)\mathbf{R}=\left(\begin{array}{ccc}b_{1} & \ldots & * \\ & \ddots & \vdots \\ 0 & & b_{n}\end{array}\right)为上三角,且b1,...bn>0b_1,...b_n>0
    • A=RHRA=R^HR为乔利斯分解
  • 平方根分解:
    • 若A>0(正定),则有A=B2,B>0(正定)A=B^2,B>0(正定),且矩阵B唯一矩阵B唯一,可写B=A,A=(A)2B=\sqrt{A},A=(\sqrt{A})^2
    • A0A\geq0半正定,则存在B0B\geq0半正定,使A=B2A=B^2,且BH=B,λ(B)={λ1,...,λn}B^H=B,\lambda(B)=\{\sqrt{\lambda_1},...,\sqrt{\lambda_n}\}全为非负根

补充知识

  • 正定阵:设A是Hermite阵,若f(x)xHAx>0f(x)-x^HAx>0对一切非0向量成立,则A为正定阵(A>0A>0)
    • 等价条件
      • A为正定
      • A有分解A=PHPA=PPHA=P^HP或A=PP^H,P可逆
      • 特根λ(A)={λ1>0,...,λn>0}\lambda(A)=\{\lambda_1 >0,...,\lambda_n >0\}(可用特商公式λ1=xHAxx2\lambda_{1}=\frac{x^{H} A x}{|x|^{2}}证明)
        证明
  • 半正定阵:设A是Hermite阵,若f(x)xHAx0f(x)-x^HAx \geq0对一切非0向量成立,则A为正定阵(A0A\geq 0)
    • 等价条件
      • A为半正定
      • A有分解A=PHPA=PPHA=P^HP或A=PP^H
      • 特根λ(A)={λ10,...,λn0}\lambda(A)=\{\lambda_1 \geq 0,...,\lambda_n \geq 0\}都非负
  • 定理
    • A=AmnA=A_{m*n},则AHA,AAHA^HA,AA^H都半正定,特征根非负(利用二次型证明)
    • A=AmnA=A_{m*n}为高阵,或方阵A可逆,则AHAA^HA正定
    • Hermite阵的特征根都是实数

奇异值

  • 正奇异值定义:给定 AmnA_{m*n},则 AHAAAHA^HA与AA^H 有相同正根λ1λ2...λr>0,r=r(A)\lambda_1 \geq \lambda_2 \geq ... \lambda_r >0,r=r(A),λ1,...λr\sqrt{\lambda_1},...\sqrt{\lambda_r}叫做A的正奇异值,全体正奇异值记作s+(A)={λ1,...λr}s_+(A)=\{\sqrt{\lambda_1},...\sqrt{\lambda_r}\},其中λ1\sqrt{\lambda_1}为A的最大奇异值
  • 奇异值定义:给定 AmnA_{m*n},则 AHAAAHA^HA与AA^H 有相同非负根λ1λ2...λr0\lambda_1 \geq \lambda_2 \geq ... \lambda_r \geq 0,λ1,...λr\sqrt{\lambda_1},...\sqrt{\lambda_r}叫做A的奇异值,全体奇异值记作s(A)={λ1,...λr}s(A)=\{\sqrt{\lambda_1},...\sqrt{\lambda_r}\}

奇异值分解

  • 正SVD分解(短SVD分解):设A=Amn,r=r(A)>0A=A_{m*n},r=r(A)>0,正奇异值为λ1,λ2,...,λr\lambda_1,\lambda_2,...,\lambda_r,则有分解A=PΔQHA=P\Delta Q^H,其中Δ=(λ100λr)\Delta=\left(\begin{array}{ccc}\sqrt{\lambda_{1}} & \cdots & 0 \\ \vdots & \ddots & \vdots \\ 0 & \cdots & \sqrt{\lambda_{r}}\end{array}\right),P,Q为半优阵,PHP=I,QHQ=IP,Q为半优阵,P^HP=I,Q^HQ=I
  • 将P,Q扩大为U阵W=(P1,P2),V=(Q1,Q2)W=(P_1,P_2),V=(Q_1,Q_2),可验证WDVH=W(Δ000)VH=PΔQH=AW \mathrm{D} V^{H}=W\left(\begin{array}{ll}\Delta & 0 \\ 0 & 0\end{array}\right) V^{H}=P \Delta Q^{H}=A
  • 奇异分解(SVD):设A有正奇异值λ1,λ2,...,λr>0\lambda_1,\lambda_2,...,\lambda_r >0,则有A=WDVHA=W \mathrm{D} V^{H},其中D=(Δ000),W=Wmm,V=Vnn为两个优阵D = \left(\begin{array}{ll}\Delta & 0 \\ 0 & 0\end{array}\right),W=W_{m*m},V=V_{n*n}为两个优阵
  • 分解方法:
    1. AHAA^HA的特征根,正奇值
    2. 求特征根对应的正交特向
    3. 令列优阵Q=(X1X1,,XrXr);Q=\left(\frac{X_{1}}{\left|X_{1}\right|}, \cdots, \frac{X_{r}}{\left|X_{r}\right|}\right) ; \quadP=(AX1AX1,,AXrAXr)\quad P=\left(\frac{A X_{1}}{\left|A X_{1}\right|}, \cdots, \frac{A X_{r}}{\left|A X_{r}\right|}\right),正SVD分解A=PΔQHA=P \Delta Q^{H}
    4. 可用观察扩充法求两个U阵,W=(P,Y),V=(Q,X)W=(P,Y),V=(Q,X),可得SVD公式A=WDVHA=W D V^{H},D=(Δ000)D = \left(\begin{array}{ll}\Delta & 0 \\ 0 & 0\end{array}\right)
    • 备注:
      • 对于V=(Q,X)V=(Q,X),可解AX=0AX=0XX
      • 对于W=(P,Y)W=(P,Y),可解AHY=0A^HY=0YY
  • 转置法:
    • BH=A=WDVHB=VDHWHB^H=A=WDV^H \rightarrow B=VD^HW^H

正规阵谱公式

正规分解复习

  • 正规阵谱公式:若A=AnnA=A_{n*n},全体互异根为λ1,...λk\lambda_{1} ,...\lambda_{k},则有A=λ1G1+...+λkGkA=\lambda_1G_1+...+\lambda_kG_k,叫做A的谱分解
    1. G1+G2+...GK=IG_1+G_2+...G_K=I
    2. G1G2=0,...,GiGj=0G_1G_2=0,...,G_iG_j=0
    3. G12=G1,...,Gk2=GkG_1^2=G_1,...,G_k^2=G_k(幂等),G1H=G1,...,GkH=GkG_1^H = G_1,...,G_k^H=G_k
    4. AP=λ1PG1+...+λkPGkA^P=\lambda_1^PG_1+...+\lambda_k^PG_k
    5. f(A)=f(λ1)G1+...+f(λk)Gkf(A)=f(\lambda_1)G_1+...+f(\lambda_k)G_k
    6. G1=(Aλ1I)(AλkI)(λ1λ1)(λ1λk)G_{1}=\frac{\left(A-\lambda_{1} I\right) \cdots \cdots\left(A-\lambda_{k} I\right)}{\left(\lambda_{1}-\lambda_{1}\right) \cdots \cdots\left(\lambda_{1}-\lambda_{k}\right)},G2=(Aλ1I)(Aλ2I)(AλkI)(λ2λ1)(λ2λ2)(λ2λk)G_{2}=\frac{\left(A-\lambda_{1} I\right)\left(A-\lambda_{2} I\right) \cdots \cdots\left(A-\lambda_{k} I\right)}{\left(\lambda_{2}-\lambda_{1}\right)\left(\lambda_{2}-\lambda_{2}\right) \cdots \cdots\left(\lambda_{2}-\lambda_{k}\right)},\cdots \cdots,Gk=(Aλ1)(Aλk)(λkλ1)(λkλk)G_{k}=\frac{\left(A-\lambda_{1}\right) \cdots \cdots\left(A-\lambda_{k}\right)}{\left(\lambda_{k}-\lambda_{1}\right) \cdots \cdots\left(\lambda_{k}-\lambda_{k}\right)}
      • 若A正规,只有两个不同根,公式G1=(Aλ2)λ1λ2G_{1}=\frac{\left(A-\lambda_{2}\right)}{\lambda_1-\lambda_2},G2=(Aλ1)λ2λ1G_{2}=\frac{\left(A-\lambda_{1}\right)}{\lambda_2-\lambda_1}G1+G2=IG_1+G_2=I
      • 若A正规,有三个不同根,公式G1=(Aλ2)(Aλ3)(λ1λ2)(λ1λ3)G_{1}=\frac{(A-\lambda_{2})(A-\lambda_{3})}{(\lambda_1-\lambda_2)(\lambda_1-\lambda_3)},G1=(Aλ1)(Aλ3)(λ2λ1)(λ2λ3)G_{1}=\frac{(A-\lambda_{1})(A-\lambda_{3})}{(\lambda_2-\lambda_1)(\lambda_2-\lambda_3)},G3=(Aλ1)(Aλ2)(λ3λ1)(λ3λ2)G_{3}=\frac{(A-\lambda_{1})(A-\lambda_{2})}{(\lambda_3-\lambda_1)(\lambda_3-\lambda_2)}
    7. AG1=λ1G1,...AGk=λkGkAG_1=\lambda_1G_1,...AG_k=\lambda_kG_k,G1,G2,...,GkG_1,G_2,...,G_k中各列都是A的特征向量
  • 分块法求Ak=(A100A2)k=(A1k00A2k)A^k=(\begin{array}{ll}A_1 & 0 \\ 0 & A_2\end{array})^k=(\begin{array}{ll}A_1^k & 0 \\ 0 & A_2^k\end{array})

单阵

  • 单阵(单纯阵,可对角阵)定义: 若A=AnnA=A_{n*n}相似于对角形,即Q1DQ=D=(λ100λn)Q^{-1}DQ=D=\left(\begin{array}{ccc}\lambda_{1} & & 0 \\ & \ddots & \\ 0 & & \lambda_{n}\end{array}\right),称A为单纯阵
    • 正规阵一定是单纯阵
    • 充分条件:设n阶方阵A恰有n个不同根,则A为单阵
    • 单阵的谱公式:若A为单阵,全体不同根为λ1,...,λk\lambda_1,...,\lambda_k,则有A=λ1G1+...+λkGkA=\lambda_1G_1+...+\lambda_kG_k
      • 单阵的谱阵G1,...GkG_1,...G_k是幂等阵,但是不一定为Hermite阵(正规阵谱分解)
  • 单阵判定法:设λ1,...,λk\lambda_1,...,\lambda_k为A的全体不同根
    • (Aλ1I)...(Aλk)I=0(A-\lambda_1I)...(A-\lambda_k)I=0则A为单阵(相似于对角形)
    • (Aλ1I)...(Aλk)I0(A-\lambda_1I)...(A-\lambda_k)I \not ={0}则A为单阵(相似于对角形)

*极小式(0化式)

  • λ1,...,λk\lambda_1,...,\lambda_k为不同根,且(Aλ1I)...(Aλk)I=0(A-\lambda_1I)...(A-\lambda_k)I=0,则A必是单阵,此时**m(x)=(xλ1)(xλ2)...(xλk)m(x)=(x-\lambda_1)(x-\lambda_2)...(x-\lambda_k)为A的极小式**
  • 0化式:若多项式f(x)符合f(A)=0,称f(x)为A的0化式
  • 极小式:若多项式m(x)符合m(A)=0且m(x)具有最小次数,称m(x)为A的极小式
    • 特征多项式一定是0化式
    • 极小式鄙视特征式的因式
    • 极小式必为每个0化式f(x)的因子
  • 极小式m(x)的求法:
    • xIA=(xa)2(xb)|xI-A|=(x-a)^2(x-b)(AaI)(AbI)=0是否成立(A-aI)(A-bI)=0是否成立,若成立则m(x)=(xa)(xb)m(x)=(x-a)(x-b),否则为m(x)=(xa)2(xb)2m(x)=(x-a)^2(x-b)^2
    • |xI-A|=(x-a)(x-b)(x-c),则极小式为m(x)=(xa)(xb)(xc)m(x)=(x-a)(x-b)(x-c)
      补充证明题

补充验单法

  • 已知n阶方阵A
    • 若每个k>1重根λ1\lambda_1使得rank(Aλ1)=nkrank(A-\lambda_1) = n-k,则A为单阵
    • 若有k>1重根λ1\lambda_1使得rank(Aλ1)nkrank(A-\lambda_1) \not ={n-k},则A非单阵

补充引理

  • (Aλ1I)P=0(A-\lambda_1I)P=0,则P中非0列都是λ1\lambda_1的特向

镜面阵

镜面阵
引理1
例题1

3个矩阵函数

  • 常见解析函数f(x)与矩阵函数A,其中A为方阵
    • ex=1+x+12!x2++1k!xk+e^{x}=1+x+\frac{1}{2 !} x^{2}+\cdots+\frac{1}{k !} x^{k}+\cdots \cdots,则eA=1+A+12!A2++1k!Ak+e^{A}=1+A+\frac{1}{2 !} A^{2}+\cdots+\frac{1}{k !} A^{k}+\cdots \cdots
      • eA=k=01k!Ake^{\mathbf{A}}=\sum_{k=0}^{\infty} \frac{1}{k !} \mathbf{A}^{k}
    • sinx=xx33!+x55!x77!+\sin x=x-\frac{x^{3}}{3 !}+\frac{x^{5}}{5 !}-\frac{x^{7}}{7 !}+\cdots \cdots,则sinA=AA33!+A55!A77!+\sin A=A-\frac{A^{3}}{3 !}+\frac{A^{5}}{5 !}-\frac{A^{7}}{7 !}+\cdots \cdots
      • sinA=k=0(1)kA2k+1(2k+1)!\sin \mathbf{A}=\sum_{\mathrm{k}=0}^{\infty}(-1)^{k} \frac{\mathbf{A}^{2 k+1}}{(2 k+1) !}
    • cosx=xx22!+x44!x66!+\cos x=x-\frac{x^{2}}{2 !}+\frac{x^{4}}{4 !}-\frac{x^{6}}{6 !}+\cdots \cdots,则cosA=AA22!+A44!A66!+\cos A=A-\frac{A^{2}}{2 !}+\frac{A^{4}}{4 !}-\frac{A^{6}}{6 !}+\cdots \cdots
      • cosA=k=0(1)kA2k(2k)!\cos \mathbf{A}=\sum_{\mathrm{k}=0}^{\infty}(-1)^{k} \frac{\mathbf{A}^{2 k}}{(2 k) !}
  • 为了方便应用引入参数t
    • etA=1+tA+12!tA2++1k!tAk+e^{tA}=1+tA+\frac{1}{2 !} tA^{2}+\cdots+\frac{1}{k !} tA^{k}+\cdots \cdots
    • sintA=tAtA33!+tA55!tA77!+\sin tA=tA-\frac{tA^{3}}{3 !}+\frac{tA^{5}}{5 !}-\frac{tA^{7}}{7 !}+\cdots \cdots
    • costA=tAtA22!+tA44!tA66!+\cos tA=tA-\frac{tA^{2}}{2 !}+\frac{tA^{4}}{4 !}-\frac{tA^{6}}{6 !}+\cdots \cdots
  • 性质:
    • e0=I,sin0=0,cos0=0e^0=I,sin0=0,cos0=0
    • cos(A)=A,sin(A)=sin(A)cos(-A)=A,sin(-A)=-sin(A)
    • eiA=cosA+isinAe^{iA}=cosA+isinA,cosA=12(eiA+eiA)\cos A=\frac{1}{2}\left(e^{i A}+e^{-i A}\right),sinA=12i(eiAeiA)\sin A=\frac{1}{2 i}\left(e^{i A}-e^{-i A}\right),eitA=costA+isintAe^{itA}=costA+isintA
    • 交换公式,若AB=BAAB=BA,则eAeB=eA+B=eBeAe^Ae^B=e^{A+B}=e^Be^A
    • 特别公式:eAeA=eAeA=Ie^Ae^{-A}=e^{-A}e^A=I
    • 可逆公式:(eA)1=eA(e^A)^{-1}=e^{-A}
    • 推论:令n方阵A=(ai,j)A=(a_{i,j}),则f(A)=eAf(A)=e^A的行列式为det(eA)=eA=etr(A)det(e^A)=|e^A|=e^{tr(A)}
  • 幂等阵,若A2=AA^2=A,则etA=I+(et1)Ae^{tA}=I+(e^t-1)A

幂零阵

  • Ak=0A^k=0,A称为幂零阵
  • 幂零阵公式1:若Ak=0A^k=0,f(x)为任一解析式,则f(A)=f(0)I+f(0)A+f(0)2!A2++fk1(0)(k1)!Ak1f(A)=f(0) I+f^{\prime}(0) A+\frac{f^{\prime \prime}(0)}{2 !} A^{2}+\cdots+\frac{f^{k-1}(0)}{(k-1) !} A^{k-1}
  • 幂零阵公式0:若(AaI)k=0(A-aI)^k=0,f(x)为任一解析式,则f(A)=f(a)I+f(a)(Aa)+f(a)2!(Aa)2++fk1(a)(k1)!(Aa)k1f(A)=f(a) I+f^{\prime}(a) (A-a)+\frac{f^{\prime \prime}(a)}{2 !} (A-a)^{2}+\cdots+\frac{f^{k-1}(a)}{(k-1) !} (A-a)^{k-1}
  • 备注1:若Ak=0A^k=0为幂零阵,则A的全体根λ(A)=0,0,...,0\lambda(A)={0,0,...,0}
  • 备注2:若(AaI)k=0(A-aI)^k=0为幂零阵,则A的全体根λ(A)=a,a,...,a\lambda(A)={a,a,...,a}

Jordan块

  • n阶上三角A=(a10a11a)n,nA=\left(\begin{array}{lllll}a & 1 & & & 0 \\ & a & 1 & & \\ & & \ddots & \ddots & \\ & & & \ddots & 1 \\ & & & & a\end{array}\right)_{n, n}叫n阶Jordan块
  • 可知λ=a\lambda=a为n重根
  • 若a=0,可得n阶0根Jordan块,A=(0100110)n,n=(0,e1,e2,...,en1A=\left(\begin{array}{lllll}0& 1 & & & 0 \\ & 0 & 1 & & \\ & & \ddots & \ddots & \\ & & & \ddots & 1 \\ & & & & 0\end{array}\right)_{n, n}=(0,e_1,e_2,...,e_{n-1},eje_{j}为单位阵I的各列
    • Dn+1=0D^{n+1}=0
  • 定理:n阶Jordan块A=(a10a11a)n,nA=\left(\begin{array}{lllll}a & 1 & & & 0 \\ & a & 1 & & \\ & & \ddots & \ddots & \\ & & & \ddots & 1 \\ & & & & a\end{array}\right)_{n, n}对任一解析式f(x)有公式:f(A)=(f(a)f(a)12!f(a)1(n1)!f(n1)(a)f(a)f(a)12!f(a)f(a)f(a)f(a))f(A)=\left(\begin{array}{ccccc}f(a) & f^{\prime}(a) & \frac{1}{2 !} f^{\prime \prime}(a) & \cdots & \frac{1}{(n-1) !} f^{(n-1)}(a) \\ & f(a) & f^{\prime}(a) & \cdots & \vdots \\ & & \ddots & \ddots & \frac{1}{2 !} f^{\prime \prime}(a) \\ & & & f(a) & f^{\prime}(a) \\ & & & & f(a)\end{array}\right)

广谱公式:非单阵的f(A)公式

  • 若A满足(AaI)2(AbI)=0,ab(A-aI)^2(A-bI)=0,a \geq b,即A有0化式(xa)2(xb)(x-a)^2(x-b),则有广谱公式f(A)=f(a)G1+f(b)G2+f(a)D1,G1+G2=If(A)=f(a)G_1+f(b)G_2+f^{'}(a)D_1,G_1+G_2=I,其中G1,G2,D1为固定矩阵(广谱阵)G_1,G_2,D_1为固定矩阵(广谱阵)
    例题
    例题

张量积

  • 定义:设A=(ai,j)mn,B=(bi,j)pqA=(a_{i,j})_{m*n},B=(b_{i,j})_{p*q},称下面的分块矩阵为A与B的直积(张量积),(a11Ba12Ba1nBa21Ba22Ba2nBam1Bam2BamnB)mp×nq\left(\begin{array}{llll}a_{11} B & a_{12} B \cdots \cdots & a_{1 n} B \\ a_{21} B & a_{22} B \cdots \cdots & a_{2 n} B \\ \cdots \cdots \cdots \cdots & & \\ a_{m 1} B & a_{m 2} B \cdots \cdots & a_{m n} B\end{array}\right)_{m p \times n q},记作AB=(aijB)mp×nqA \otimes B=\left(a_{i j} B\right)_{m p \times n q}
  • 定理1:
    • 两个上三角阵的直积意识上三角阵
    • 两个对角阵的直积仍是对角阵
    • InIm=ImIn=ImnI_n \otimes I_m = I_m \otimes I_n = I_{m*n}
  • 推论1:设A=(α1,...,αt)nt,b=(β1,...,βt)pq(按列分块)A=(\alpha_1,...,\alpha_t)_{n*t},b=(\beta_1,...,\beta_t)_{p*q}(按列分块),则AB=(α1β1,...,αtβp)nptqA \otimes B=(\alpha_1\otimes\beta_1,...,\alpha_t\otimes\beta_p)_{np*tq},记为AB=(全列αiβj)按字典次序A\otimes B = (全列\alpha_i\otimes\beta_j)--按字典次序
  • 性质:
    • 直积满足分配律,结合律,倍数关系
    • 吸收律:(AB)(CD)=(AC)(BD)(A\otimes B)(C \otimes D)=(AC)\otimes(BD)
    • 转置公式:(AB)T=ATBT,(AB)H=AHBH(A\otimes B)^T = A^T \otimes B^T,(A \otimes B)^H=A^H \otimes B^H
    • (AB)1=A1B1(A \otimes B)^{-1}=A^{-1} \otimes B^{-1}
    • 秩公式:rank(AB)=(rankA)(rankB)rank(A \otimes B) = (rankA)(rankB)
    • 设A为m阶方阵,B为n阶方阵,则tr(AB)=trAtrB,det(AB)=(detA)n(detB)mtr(A \otimes B) =trA*trB,det(A \otimes B)=(detA)^n(detB)^m
  • 推论2:若A为m阶方阵,B为n阶方阵,则
    • (AB)k=AkBk(A\otimes B)^k=A^k \otimes B^k
    • (AIn)(ImB)=(ImB)(AIn)=AB(A \otimes I_n)(I_m \otimes B)=(I_m \otimes B)(A \otimes I_n) = A \otimes B
  • 推论3:若A,B为优阵,则ABA \otimes B为优阵
  • 推论4:设X1,...,XpX_1,...,X_p式p个线性无关的列向量,Y1,...,YqY_1,...,Y_q式q个线性无关的列向量,则全体pq个向量XiYjX_i \otimes Y_j线性无关

张量积与特征值

  • 定理1:设m阶方阵A的特征根为λ1,...λm\lambda_1,...\lambda_m,n阶方阵B的特征根为t1,...,tnt_1,...,t_n,则ABA \otimes B的全体特征根恰为mn个λktj\lambda_k t_j
  • 定理2:设m阶方阵A的特征根为λ1,...λm\lambda_1,...\lambda_m,n阶方阵B的特征根为t1,...,tnt_1,...,t_n
    • AIn+ImBA \otimes I_n + I_m \otimes B的全体特征根恰为mn个λk+tj\lambda_k +t_j
    • AInImBA \otimes I_n - I_m \otimes B的全体特征根恰为mn个λktj\lambda_k -t_j

**矩阵的拉直

  • 按行拉直为一个列向量,A=(a11,a12,,a1n,a21,a22,,a2n,,am1,am2,,amn)T\vec{A}=\left(a_{11}, a_{12}, \cdots, a_{1 n}, a_{21}, a_{22}, \cdots, a_{2 n}, \cdots, a_{m 1}, a_{m 2}, \cdots, a_{m n}\right)^{T}
  • 性质:
    • 线性公式A+b=A+B\vec{A+b}=\vec{A}+\vec{B},kA=kA\vec{kA}=k\vec{A}
    • 拉直公式:ABC=(ACT)B\vec{ABC} = (A \otimes C^T)\vec{B}
  • 推论:
    • AX=(AIn)X\vec{AX} = (A \otimes I_n)\vec{X}
    • XB=(ImBT)X\vec{XB} = (I_m \otimes B^T)\vec{X}
    • AX+XB=(AIn+ImBT)X\vec{AX+XB} = (A \otimes I_n + I_m \otimes B^T)\vec{X}

方程求解

  • AX+XB=CAX+XB=C:,利用拉直公式AX+XB=(AIn+ImBT)X\vec{AX+XB} = (A \otimes I_n + I_m \otimes B^T)\vec{X}将方程拉直为(AIn+ImBT)X=C(A \otimes I_n + I_m \otimes B^T)\vec{X} = \vec{C}

广义逆

  • A+A^+定义:若m*n矩阵A与n*m矩阵适合四个条件,则X叫A的加号逆(伪逆,广义逆),即为X=A+X=A^+

    1. AXA=A
    2. XAX=X
    3. (AX)H=AX(AX)^H=AX
    4. (XA)H=XA(XA)^H=XA
  • A+A^+唯一定理:矩阵A存在A+A^+且唯一

  • 性质

    • AA+A=AAA^+A=A
    • A+AA+=A+A^+AA^+=A^+
    • (AA+)H=AA+(AA^+)^H=AA^+
    • (A+A)H=A+A(A^+A)^H=A^+A
    • 幂等公式:(AA+)2=AA+,(A+A)2=A+A(AA^+)^2=AA^+,(A^+A)^2=A^+A
  • 例题

    • 若A为优阵,则A+=AH=A1A^+=A^H=A^{-1}(假设A+=AHA^+=A^H,验证四个条件)
    • 若A为半优阵,则A+=AHA^+=A^H(假设A+=AHA^+=A^H,验证四个条件)
      例题
      例题
  • 引理1(半优公式):

    • 若A为列优阵AHA=IA^HA=I,则A+=AHA^+=A^H
    • AHA^H为列优阵AAH=IAA^H=I(A为行优阵),则A+=AHA^+=A^H
  • A+A^+第一公式:若A=AmnA=A_{m*n}有正SVD,A=PΔQH,(P,Q为半优阵)A=P\Delta Q^H,(P,Q为半优阵),则有公式A+=QΔ1PH,Δ+=Δ1A^+=Q\Delta^{-1}P^H,\Delta^+=\Delta^{-1}

  • A+A^+秩1公式:若r(A)=1,各列成比例r(A)=1,各列成比例,则A+=1aij2AH=1tr(AHA)AHA^{+}=\frac{1}{\sum\left|a_{i j}\right|^{2}} A^{H}=\frac{1}{tr(A^HA)} A^H

  • 分块公式:A=(A100A2)m×nA=\left(\begin{array}{cc}A_{1} & 0 \\ 0 & A_{2}\end{array}\right)_{\operatorname{m} \times n},则A+=(A1+00A2+)m×nA^+=\left(\begin{array}{cc}A_{1}^+ & 0 \\ 0 & A_{2}^+\end{array}\right)_{\operatorname{m} \times n}

  • 特注:

    • 一般情况下A+AI,AA+I,A+AAA+A^+A \not ={I},AA^+ \not ={I},A^+A \not ={AA^+}
    • 一般情况下(AB)+A+B+(AB)^+ \not ={A^+B^+}
    • 一般情况下(AB)+B+A+(AB)^+ \not ={B^+A^+},只有满秩分解时成立
  • 引理2(优分解公式):设P,Q为优阵,则(PDQ)+=QHD+PH=Q+D+P+(PDQ)^+=Q^HD^+P^H=Q^+D^+P^+

  • 引理3:设D为对角形,则D+D=DD+D^+D=DD^+

  • 定理:若A正规,则A+A=AA+A^+A=AA^+,且(A+)k=(Ak)+(A^+)^k=(A^k)^+

  • QR公式:设A为高阵,且有QR分解,A=QR,Q为列优阵,则A+=R1QHA^+=R^{-1}Q^H

  • A+A^+的谱分解公式:A+=λ1+G1++λk+GkA^{+}=\lambda_{1}^{+} G_{1}+\cdots+\lambda_{k}^{+} G_{k},A+A^+非正规时不一定成立,A为单阵且可逆时一定成立

    • 推论:若A正规,则A+A=AA+,(A+)p=(Ap)+A^+A=AA^+,(A^+)^p=(A^p)^+
    • 推论:若A正规,则A+正规,(A+)p=(Ap)+A^+正规,且(A^+)^p=(A^p)^+
  • 补充A+A^+公式:

    • A+=(AHA)+AHA^+=(A^HA)^+A^H
    • A+=AH(AAH)+A^+=A^H(AA^H)^+

应用

  • 高低阵引理:

    • 若B为高阵(列满秩),则存在左逆阵BLB_L使得BLB=IB_LB=I,其中BL=(BHB)1BHB_L=(B^HB)^{-1}B^H
    • 若C为低阵(行满秩),则存在右逆阵CRC_R使得CCR=ICC_R=I,其中CR=CH(CCH)1C_R=C^H(CC^H)^{-1}
  • 高低广逆公式:

    • 若A为高阵(列满秩),加号逆A+=(AHA)1AHA^+=(A^HA)^{-1}A^H
    • 若A为低阵(行满秩),加号逆A+=AH(AAH)1A^+=A^H(AA^H)^{-1}
  • A+A^+第二公式:若有高低分解A=BC(满秩分解),则有公式A+=C+B+A^+=C^+B^+,其中C+=CH(CCH)1C^+=C^H(CC^H)^{-1},B+=(BHB)1BHB^+=(B^HB)^{-1}B^H,B+B=I,CC+=IB^+B=I,CC^+=I

  • A+A^+第三公式:A+=(AHA)+AHA^+=(A^HA)^+A^H

  • 张量的广逆公式:(AB)+=A+B+(A \otimes B)^+=A^+ \otimes B^+

  • 常见分块公式:

    • (A100A2)+=(A1+00A2+)\left(\begin{array}{cc}A_{1} & 0 \\ 0 & A_{2}\end{array}\right)^{+}=\left(\begin{array}{cc}A_{1}^{+} & 0 \\ 0 & A_{2}^{+}\end{array}\right)
    • (0A1A20)+=(0A2+A1+0)\left(\begin{array}{cc}0 & A_{1} \\ A_{2} & 0\end{array}\right)^{+}=\left(\begin{array}{cc}0 & A_{2}^{+} \\ A_{1}^{+} & 0\end{array}\right)

正规方程

  • AHAx=AHbA^HAx=A^HbAx=bAx=b的正规方程
  • 正规方程AHAx=AHbA^HAx=A^Hb必有特解x0=A+bx_0=A^+b,使AHAx0=AHbA^HAx_0=A^Hb
  • 推论:正规方程AHAx=AHbA^HAx=A^Hb必有解,方程Ax=bAx=b可能无解
  • 特解定理:若Ax=bAx=b有解(相容),则必有特解x0=A+bx_0=A^+b
  • 无解定理:若x0=A+b使Ax0bx_0=A^+b使Ax_0 \not ={b},则Ax=bAx=b无解(不相容)

子空间

  • 子空间定义:设WCn(非空集合)W \subset \mathrm{C}^{\mathrm{n}}(非空集合),若对α,βW\forall \boldsymbol{\alpha}, \boldsymbol{\beta} \in W, 有 α+βW\boldsymbol{\alpha}+\boldsymbol{\beta} \in W(对加法封闭),对αW,kC\forall \boldsymbol{\alpha} \in W, \boldsymbol{k} \in C, 有 kαW\boldsymbol{k\alpha} \in W(对数乘封闭),称W为CnC^n中的子空间(简称空间)

  • 推论1:子空间W一定含有0W\vec{0} \in W,且子空间对线性组合封闭

  • 推论2:子空间W一定含有0W\vec{0} \in W;若W不含0向量,则W不是子空间

  • 生成子空间:设α1,...αmCn,W=x=k1α1+...+kmαmkjC,j=1,...,m\alpha_1,...\alpha_m \in C^n,W={x=k_1\alpha_1+...+k_m\alpha_m|k_j \in C,j=1,...,m},称W为由α1,...αm\alpha_1,...\alpha_m生成的子空间,记为W=L(α1,...,αm)W=L(\alpha_1,...,\alpha_m)

  • 基定义:若子空间W中有r个向量α1,...,αr\alpha_1,...,\alpha_r线性无关,W中任一α\alpha可由α1,...,αr\alpha_1,...,\alpha_r表示,称[α1,...,αr][\alpha_1,...,\alpha_r]为W的一个基,r叫做W的维数,记作dimW=rdim W =r

    • 基性质1:设W是r维子空间(dim W=r),则W中任r+1个向量必线性相关
    • 基性质2:W的基[α1,...,αr][\alpha_1,...,\alpha_r]必是向量组W的一个极大无关组
    • 基性质3:若dim W=r,则W中任取r个线性无关向量都是W的基
  • 备注:矩阵A=AmnCmnA=A_{m*n} \in C^{m*n}产生两个子空间:值域R(A)CmR(A) \subset C^m,核空间N(A)Cn\mathrm{N}(A) \subset \mathrm{C}^{\mathrm{n}}

    • A=Am,nA=A_{m,n}的值域为R(A)=y=AXXCnCmR(A)={y=AX|X \in C^n} \subset C^m也叫列空间,令A=(α1,...αn)A=(\alpha_1,...\alpha_n),可写R(A)=L(α1,...,αn)R(A)=L(\alpha_1,...,\alpha_n)
      • 维数=秩数:dimR(A)=rank(A)dim R(A) = rank(A)
      • 值域可以用CmC^m中过原点的一个平面表示
    • A=Am,nA=A_{m,n}的核空间为N(A)=XCmAX=0CnN(A)={X \in C^m | AX=0} \subset C^n
      • 维数公式:dimN(A)=nrank(A)dim N(A) = n- rank(A)
      • 核空间可以用CmC^m中过原点的一个平面表示
  • 维数公式:dim N(A) + dim R(A) = n

  • 正交引理1:

    • CnC^n中正交引理:任取向量bCnb \in C^n,令x0=A+bCnx_0=A^+b \in C^n,则有x0N(A)x_{0} \perp \mathrm{N}(A),即x0x,xN(A)=xAx=0x_{0} \perp x ,\forall x \in N(A)={x|Ax=0}
      • 推论1:若x0=A+b,x0+x2=x02+x2X02(x02最小)x_0=A^+b,则|x_0+x|^2=|x_0|^2+|x|^2 \geq |X_0|^2(|x_0|^2最小)
      • 推论2:若Ax=bAx=b有特解x0x_0,则通解公式为x=x0+w,wN(A)x=x_0+w,\forall w \in N(A)
      • 最小范数解定理:若Ax=bAx=b有解(相容),则x0=A+bx_0=A^+b为最小范数解
  • 正交引理2:

    • CmC^m中正交引理:任取向量bCmb \in C^m,令x0=A+bCnx_0=A^+b \in C^n,则有(bAx0)R(A)(b-Ax_{0}) \perp \mathrm{R}(A)
    • 推论1:若x0=A+b,Axb2Ax0b2(为最小)x_0=A^+b,则|Ax-b|^2 \geq |Ax_0-b|^2(为最小)
    • 结论:仅当A(xx0)=0A(x-x_0)=0Axb2|Ax-b|^2取到最小值Axb2=Ax0b2|Ax-b|^2 = |Ax_0-b|^2

最小二乘解

  • 定义:若Ax0b2\left|A x_{0}-b\right|^{2}取到最小:Ax0b2=min{Axb2xCn}\left|A x_{0}-b\right|^{2}=\min \left\{|A x-b|^{2} \mid x \in C^{n}\right\}最小值,则x0x_0Ax=bAx=b的一个最小二乘解(极小二乘解), 简称“小二解”
  • 定理1:若Ax=b无解(不相容)Ax=b无解(不相容),则x0=A+bx_0=A^+b恰是一个最小二乘解,且全体最小二乘解公式为x=x0+w,wN(A)x=x_0+w,w \in N(A),即x=x0+w=x0+(t1Y1+...+tkYk),w=t1Y1+...+tkYk,k=nr(A),其中Y1,...,YkAx=0的基本解x=x_0+w=x_0+(t_1Y_1+...+t_kY_k),w=t_1Y_1+...+t_kY_k,k=n-r(A),其中Y_1,...,Y_k为Ax=0的基本解
  • 定理2:若Ax=b无解(不相容)Ax=b无解(不相容),则x0=A+bx_0=A^+b是最小范数的最小二乘解(最佳小二解)
  • 备注:
    • A=AmnA=A_{m*n}列满秩(高阵),则有N(A)=0N(A)={\vec{0}},即AX=0AX=0只有零解X=0X=\vec{0}
    • Ax=bAx=b不相容,且A为高阵,则x0=A+bx_0=A^+b时唯一最小二乘解(也是最佳二乘解)

**其他通解公式

  • 定理1:齐次方程Ax=0Ax=0必有通解公式x=(InA+A)y,yCnx=(I_n-A^+A)y, \forall y \in C^n
    • 备注:可写核空间公式N(A)=w=(InA+A)yyCnN(A)={w=(I_n-A^+A)y|y \in C^n}
  • 定理2:非齐次方程Ax=bAx=b必有解,则有通解公式x=(A+b)+(InA+A)y,yCnx=(A^+b)+(I_n-A^+A)y, \forall y \in C^n
  • 定理3:Ax=bAx=b的全体小二解公式为x=x0+w,wN(A),x0=A+b,x=(A+b)+(InA+A)yx=x_0+w,w \in N(A),x_0=A^+b,x=(A^+b)+(I_n-A^+A)y
  • 引理:A=Amn,bCmA=A_{m*n},b \in C^m,则x0=A+bx_0=A^+b适合AHAx0AHb=0A^HAx_0-A^Hb=0
  • 推论:A=Amn,binCm,AHAx=AHb\forall A=A_{m*n}, \forall b in C^m,\rightarrow A^HAx=A^Hb必有解

广逆补充**

  • 定理1:若矩阵方程AXB=DAXB=D有解(相容),则有特解x0=A+DB+x_0=A^+DB^+
    • 推论(无解定理):若x0=A+DB+x_0=A^+DB^+使得AX0BDAX_0B \not ={D},则无解(不相容)
    • 备注1:齐次方程AXB=0AXB=0必有通解公式X=YA+AYBB+X=Y-A^+AYBB^+,Y是适当阶数的任一矩阵
  • 定理2:若AXB=DAXB=D有解,令X0=A+DB+X_0=A^+DB^+,必有通解公式X=X0+(YA+AYBB+)X=X_0+(Y-A^+AYBB^+)
  • 特别结论矩阵方程AXA=AAXA=A必有解,且特解x0=A+AA+=A+x_0=A^+AA^+=A^+,且有通解X=X0+(YA+AYAA+)X=X_0+(Y-A^+AYAA^+)

其他广逆简介

  • 1号逆或减号逆定义:若AXA=A,则称X为A的一个减号逆,记作AA(1)A^-或A^(1)

  • 减号逆满足AAA=AAA^-A=A

  • 求法:任意矩阵的全体减号逆通解公式为A=X=A++(YA+AYAA+)A^-=X=A^++(Y-A^+AYAA^+)

  • 引理1:若A=(Ir000)m,nA=\left(\begin{array}{cc}I_{r} & 0 \\ 0 & 0\end{array}\right)_{m, n}为标准行,则全体A=(IrBCD)n,mA^{-}=\left(\begin{array}{ll}I_{r} & B \\ C & D\end{array}\right)_{n, m}

  • 补充公式1:若PAQ=(Ir000)m,nP A Q=\left(\begin{array}{ll}I_{r} & 0 \\ 0 & 0\end{array}\right)_{m, n},则全体A=Q(IrBCD)n,mPA^{-}=Q\left(\begin{array}{ll}I_{r} & B \\ C & D\end{array}\right)_{n, m} P
    求法
    补充定理1

  • 补充定理2:N(A)AY=0N(A)或AY=0的通解为Y=(IAA)y,yCnY=(I-A^-A)y,\forall y \in C^n

  • 补充定理3:若Ax=bAx=b有解,令x0=Abx_0=A^-b,则有通解公式x=x0+(InAA)y,yCnx=x_0+(I_n-A^-A)y,\forall y \in C^n

  • 性质:

    • (AA)2=AA(A^-A)^2=A^-A
    • (AA)2=AA(AA^-)^2=AA^-
    • rank(AA)=rank(A)rank(AA^-)=rank(A)
    • rank(AA)=rank(A)rank(A^-A)=rank(A)
  • 定义1:A(1,3)=xAXA=A,(AX)H=AXA^{(1,3)}={x|AXA=A,(AX)^H=AX}

  • 定义2:A(1,4)=xAXA=A,(XA)H=XAA^{(1,4)}={x|AXA=A,(XA)^H=XA}

  • 引理1:A(1,3)=xAHAX=AHA^{(1,3)}={x|A^HAX=A^H}

  • 引理2:A(1,4)=xXAAH=AHA^{(1,4)}={x|XAA^H=A^H}

  • 引理3:A(1,3)=xAX=AA+A^{(1,3)}={x|AX=AA^+},A(1,4)=xXA=A+AA^{(1,4)}={x|XA=A^+A}

矩阵函数求导公式

  • A(x)=(ai,j(x))mnA(x)=(a_{i,j}(x))_{m*n}的元素都是x的可导函数,定义A(x)A(x)关于x的求导为A(x)=ddxA(x)=(ddxaij(x))m×nA^{\prime}(x)=\frac{d}{d x} A(x)=\left(\frac{d}{d x} a_{i j}(x)\right)_{m \times n},同理可得积分式
  • 常见三个矩阵函数
    • f(At)=etA,cos(tA),sin(tA)f(At)=e^{tA},cos(tA),sin(tA)
  • 定理:设方阵ACnnA \in C^{n*n},则
    • ddxeAt=AeAt=eAtA\frac{d}{d x} e^{At} = Ae^{At}= e^{At}A
    • ddxsinAt=AcosAt=cosAtA\frac{d}{d x} sin{At} = Acos{At}= cos{At}A
    • ddxcosAt=AsinAt=(sinAt)A\frac{d}{d x} cos{At} = -Asin{At}= -(sinAt)A
  • 备注
    • 公式ddxeAt=AeAt=eAtA\frac{d}{d x} e^{At} = Ae^{At}= e^{At}A用-A替换A可得ddxeAt=AeAt\frac{d}{d x} e^{-At} = -Ae^{-At}
    • ddt(ettX(t))=AetAX(t)+etAX(t)=etA(X(t)AX(t))\frac{d}{d \mathrm{t}}\left(e^{-t t} \mathrm{X}(\mathrm{t})\right)=-A e^{-t A} \mathrm{X}(\mathrm{t})+e^{-t A} \mathrm{X}^{\prime}(\mathrm{t})=e^{-t A}\left(\mathrm{X}^{\prime}(\mathrm{t})-A \mathrm{X}(\mathrm{t})\right)

求解微分方程

  • 定理1:齐次方程dxdt=Ax,x(0)=c\frac{dx}{d t}= Ax,x(0)= \vec{c},其中x=(x1(t),...,xn(t))T,A=Annx=(x_1(t),...,x_n(t))^T,A=A_{n*n}为常数矩阵,有唯一解公式x=eAtc,x=eAtx(0)x=e^{At}\vec{c},即x=e^{At}x(0)

  • 定理2:齐次方程dYdt=AY+YB,Y(0)=D,Y=(yij(t))n×p,A=An,n,B=Bp,p\frac{d Y}{d t}=A Y+Y B, Y(0)=D, Y=\left(y_{i j}(t)\right)_{n \times p}, A=A_{n, n}, B=B_{p, p}有唯一解公式Y=e(At)De(Bt)Y=e^(At)De^(Bt),即Y=eAtY(0)eBtY=e^{At}Y(0)e^{Bt}
    例题解2阶常微分方程

  • 备注:非齐次方程{dx(t)dt=Ax(t)+f(t),x(t0)=c( 初始条件 )\left\{\begin{array}{l}\frac{d x(t)}{d t}=\mathrm{A} x(t)+f(t), \\ x\left(t_{0}\right)=\vec{c}(\text { 初始条件 })\end{array}\right.有唯一解x(t)=eA(tt0)c+eAtt0teAuf(u)dux(t)=e^{A\left(t-t_{0}\right)} \vec{c}+e^{A t} \int_{t_{0}}^{t} \mathrm{e}^{-A u} f(u) d u

例题

  • 备注公式:
    • A=(0110)A=\left(\begin{array}{cc}0 & 1 \\ -1 & 0\end{array}\right) 则有 etA=(costsintsintcost)e^{tA}=\left(\begin{array}{cc}\cos t & \sin t \\ -\sin t & \cos t\end{array}\right)
    • e(0atat0)=(cosatsinatsinatcosat)e^{\left(\begin{array}{cc}0 & a t \\ -a t & 0\end{array}\right)}=\left(\begin{array}{cc}\cos a t & \sin a t \\ -\sin a t & \cos a t\end{array}\right)

范数理论

  • 谱半径定义:称ρ(A)=max{λ1,λ2,,λn}\rho(A)=\max \left\{\left|\lambda_{1}\right|,\left|\lambda_{2}\right|, \cdots,\left|\lambda_{n}\right|\right\}为方阵的谱半径,其中方阵A的特征根λ(A)={λ1,...,λn}\lambda(A)=\{\lambda_1,...,\lambda_n\}
  • 备注:任一方阵A=AnnA=A_{n*n},必有ρ(A)0(非负性)\rho(A) \geq 0 (非负性)
  • 谱半径性质:(齐次公式)ρ(kA)=kρ(A)\rho (kA) =|k| \rho (A)
  • 幂公式:ρ(Ak)=kρ(A)\rho(A^k)=|k| \rho(A)
  • 谱范不等式:ρ(A)A\rho(A) \leq ||A||,对一起矩阵范数||A||成立

向量范数

  • 向量空间CnC^n中的模长公式:x=(x,x)=xHx=tr(xHx)=x12+...+xn2|x|=\sqrt{(x,x)}=\sqrt{x^Hx}=\sqrt{tr(x^Hx)}=\sqrt{|x_1|^2+...+|x_n|^2}
  • 性质:正性,齐性,三角性
  • 线性空间(向量空间)中范数定义
    • 定义1:设V是数域F(实数或复数域)上线性空间,若对于任一xVx \in V,对应一个非负数,记为x||x||,满足一下3个条件,则称x||x||为空间V上一个向量范数:1.正性,2.齐次性,3.三角不等式
      • 备注:空间V上一个向量范数就是V上一个非负函数ϕ(x)=x,xV\phi(x)=||x||,x \in V满足三个条件:1.正性,2.齐次性,3.三角不等式
    • 定义2:若线性空间V上有一个函数ϕ(x),xV\phi (x),x \in V适合1.正性,2.齐次性,3.三角不等式,则称ϕ(x)\phi (x)为V上的一个范数,记作ϕ(x)=x\phi (x) = ||x||
  • 常用范数,令向量x=(x1,...,xn)TCnx=(x_1,...,x_n)^T \in C^n
    • 1 范数: x1=xj=x1++xn\|\mathrm{x}\|_{1}=\sum\left|x_{j}\right|=\left|x_{1}\right|+\cdots+\left|x_{n}\right|
    • 2 范数: x2=x=(x,x)=x12++xn2\|\mathrm{x}\|_{2}=|\mathrm{x}|=\sqrt{(\mathrm{x}, \mathrm{x})}=\sqrt{\left|x_{1}\right|^{2}+\cdots+\left|x_{n}\right|^{2}} \infty 范数: x=max{x1,,xn}\|\mathrm{x}\|_{\infty}=\max \left\{\left|x_{1}\right|, \cdots,\left|x_{n}\right|\right\}
    • pp-范数: xp=(i=1nxip)1p,p1;\|\mathrm{x}\|_{p}=\left(\sum_{i=1}^{n}\left|x_{i}\right|^{p}\right)^{\frac{1}{p}}, \quad p \geq 1 ;
  • 向量范数的性质:
    • 单位化公式
    • x=x||-x||=-||x||
    • xyxy||x-y|| \geq ||x||-||y||
  • CnC^n上范数等价性:
    • 定理CnC^n上任两个范数xa,xb||x||_a,||x||_b存在任意正数k1>0,k2>0k_1>0,k_2>0使k_{1}|\mathrm{x}|{b} \leq|\mathrm{x}|{a} \leq k_{2}|\mathrm{x}|_{b}对一切x成立,简称这两个范数等价
  • 收敛定义:设CnC^n中向量序列x(k)=(x1(k),...,xn(k)),α=(α1,...,αn)Tx^{(k)}=(x_1^{(k)},...,x_n^{(k)}),\alpha=(\alpha_1,...,\alpha_n)^T,若 x1(k)a1,x2(k)a2,,xn(k)an,(k)x_{1}^{(k)} \rightarrow a_{1}, x_{2}^{(k)} \rightarrow a_{2}, \cdots, x_{n}^{(k)} \rightarrow a_{n},(k \rightarrow \infty),称 x(k)α(k),x^{(k)} \rightarrow \alpha(k \rightarrow \infty), \quadlimx(k)=α\lim x^{(k)}=\alpha
  • 收敛引理:x(k)αx(k)a0x^{(k)} \rightarrow \alpha \leftrightarrow ||x^{(k)}-a|| \rightarrow 0

矩阵范数

  • 矩阵范数定义1:对于任意方阵A,矩阵范数A||A||表示按某个法则与A对应的非负函数ϕ(A)\phi (A)满足四个条件:正性,齐性,三角性,相容性(次乘性),满足四个条件的矩阵范数A||A||叫做相容范数
  • 常用矩阵范数:
    • 列范数 (最大列和): A1=maxji=1naijj=1,2,,n;\|\mathbf{A}\|_{1}=\max_{j} \sum_{i=1}^{n}\left|a_{i j}\right| ( j=1,2, \cdots, n ;
    • 行范数 (最大行和): A=maxij=1naij(i=1,2,,,n)\|\mathbf{A}\|_{\infty}=\max_{i} \sum_{j=1}^{n}\left|a_{i j}\right|(i=1,2,, \cdots, n);
    • 谱范数: A2=(λ1(AHA))1/2,λ1(AHA)\|\mathbf{A}\|_{2}=\left(\lambda_{1}\left(\mathbf{A}^{H} \mathbf{A}\right)\right)^{1 / 2}, \lambda_{1}\left(\mathbf{A}^{H} \mathbf{A}\right) 表示 AHA\mathbf{A}^{H} \mathbf{A} 的最大特征值,即 A2\|\mathbf{A}\|_{2}A\mathbf{A} 的最大奇异值;
    • 总和范数 (求总和): AM=aij\|\mathbf{A}\|_{M}=\sum\left|a_{i j}\right|;
    • F\mathrm{F}- 范数: AF=(aij2)12=tr(AHA)\|\mathbf{A}\|_{\mathrm{F}}=\left(\sum\left|a_{i j}\right|^{2}\right)^{\frac{1}{2}}=\sqrt{\operatorname{tr}\left(\mathbf{A}^{H} \mathbf{A}\right)};
    • GG-范数: AG=nmaxi,j{\|\mathbf{A}\|_{G}=n \cdot \max_{i, j}\left\{\right. 最大的 aij}\left.\left|a_{i j}\right|\right\}
  • 备注
    • 以上6个范数满足矩阵范数的4个条件
    • 任两个范数都是等价的
  • 可用如下记号:
    • \infty 范数: A=max{A11,A21,,An1}\|A\|_{\infty}=\max \left\{\left\|A_{1}\right\|_{1},\left\|A_{2}\right\|_{1}, \cdots,\left\|A_{n}\right\|_{1}\right\} \quad (行范数)
    • 1 范数: A1=max{α11,α21,,αn1}\|A\|_{1}=\max \left\{\left\|\alpha_{1}\right\|_{1},\left\|\alpha_{2}\right\|_{1}, \cdots,\left\|\alpha_{n}\right\|_{1}\right\} \quad (列范数)
    • 2 范数(谱范数): A2= 最大的 λ1(AHA)=λ1\|A\|_{2}=\sqrt{\text { 最大的 } \lambda_{1}\left(A^{H} A\right)}=\sqrt{\lambda_{1}} (最大奇异值)
    • 都有
      • (1)正性;
      • (2)齐性;
      • (3)三角性;
      • (4)相容性: ABAB\|A B\| \leq\|A\| \cdot\|B\|
  • 备注*(公式):
    • A=AH1\|A\|_{\infty}=\left\|A^{H}\right\|_{1}
    • A1=AH\|A\|_{1}=\left\|A^{H}\right\|_{\infty}
    • A2=AH2\|A\|_{2}=\left\|A^{H}\right\|_{2}
    • AF=AHF=(tr(AHA))1/2|\mathbf{A}\|_{\mathrm{F}}=\left\|\mathbf{A}^{\mathrm{H}}\right\|_{\mathrm{F}}=\left(\mathrm{tr}\left(\mathbf{A}^{H} \mathbf{A}\right)\right)^{1 / 2}
  • 补充定理
    • U,VU, V 为西阵, 则 UAF=AVF=UAVF=AF\|\mathrm{UA}\|_{\mathrm{F}}=\|\mathbf{A} \mathrm{V}\|_{\mathrm{F}}=\|\mathrm{UAV}\|_{\mathrm{F}}=\|\mathbf{A}\|_{\mathrm{F}}
    • AvCn×n,xCn,\mathbf{A} \in \mathbb{v} \in \mathbb{C}^{\mathrm{n} \times n}, \mathrm{x} \in \mathrm{C}^{\mathrm{n}}, \quadAx2AFx2\|\mathbf{A x}\|_{2} \leq\|\mathbf{A}\|_{\mathrm{F}}\|\mathrm{x}\|_{2}
  • 新范数公式:
    • 已知CnnC^{n*n}上一个矩阵范数.||.||,P是可逆阵,令φ(A)=P1AP\varphi(\mathbf{A})=\left\|P^{-1} \mathbf{A} P\right\|,或φ(A)=PAP1,ACn×n\varphi(\mathbf{A})=\left\|P \mathbf{A} P^{-1}\right\|, \mathbf{A} \in \mathrm{C}^{\mathrm{n} \times n},可记这个新矩阵范数为φ(A)=Ap(P有关)\varphi(A)=||A||_p(与P有关)
  • 备注
    1. A±BA+B\|A \pm B\| \leq\|A\|+\|B\|, 推广有 A1±A2±±AkA1+A2++Ak\left\|A_{1} \pm A_{2} \pm \cdots \pm A_{k}\right\| \leq\left\|A_{1}\right\|+\left\|A_{2}\right\|+\cdots+\left\|A_{k}\right\|
    2. ABAB\|A B\| \leq\|A\| \cdot\|B\| 推广有 A1A2AkA1A2Ak\left\|A_{1} A_{2} \cdots \cdots \cdot A_{k}\right\| \leq\left\|A_{1}\right\| \cdot\left\|A_{2}\right\| \cdots \cdots \cdot\left\|A_{k}\right\|
    3. 幂公式: AkAk(k=1,2,3)\left\|A^{k}\right\| \leq\|A\|^{k} \quad(k=1,2,3 \cdots \cdots)
    4. 幕公式: ρ(Ak)=[ρ(A)]k,k=1,2,3\rho\left(A^{k}\right)=[\rho(A)]^{k}, k=1,2,3 \cdots \cdots

矩阵范数产生向量范数

  • 定理1:CnnC^{n*n}上任一矩阵范数.||.||都产生一个向量范数Φ(X)=XV\varPhi(X) = ||X||_V,使φ(AX)Aφ(X), 即 AXVAXV,ACn,n,XCn\varphi(A X) \leq\|A\| \varphi(X),\text { 即 }\|A X\|_{\mathrm{V}} \leq\|A\|\|X\|_{\mathrm{V}}, \quad \forall A \in \mathrm{C}^{\mathrm{n}, \mathrm{n}}, \forall X \in \mathrm{C}^{\mathrm{n}}
  • 备注(定义):若矩阵范数Am||A||_m与向量范数xv||x||_v适合AXVAmXV\|A X\|_{\mathrm{V}} \leq\|A\|_m\|X\|_{\mathrm{V}},则说Amxv相容||A||_m与||x||_v相容
  • 向量范数生成公式:Cn×n\mathrm{C}^{n \times n} 上任一矩阵范数 A\|A\| 产生一个向量范数公式:φ(X)=Xv(x100xn00),X=(x1xn)Cn\varphi(X)=\|X\|_{\mathrm{v}} \triangleq\left\|\left(\begin{array}{cccc}x_{1} & 0 & \cdots & 0 \\\vdots & \vdots & \ddots & \vdots \\x_{n} & 0 & \cdots & 0\end{array}\right)\right\|, X=\left(\begin{array}{c}x_{1} \\\vdots \\x_{n}\end{array}\right) \in \mathrm{C}^{\mathrm{n}},且Xv\|X\|_{\mathrm{v}}满足相容性:AXvAXv,ACn,n,XCn\quad\|A X\|_{\mathrm{v}} \leq\|A\|\|X\|_{\mathrm{v}}, \quad A \in \mathrm{C}^{\mathrm{n}, \mathrm{n}}, X \in \mathrm{C}^{\mathrm{n}}
  • 小范数定理:设 ACn×n\mathbf{A} \in \mathrm{C}^{\mathrm{n} \times n} 固定, 任取很小正数 ε>0\varepsilon>0, 则有矩阵范数 ()ε\|(\bullet)\|_{\varepsilon} 使得 St:Aερ(A)+ε\mathrm{St}:\|\mathbf{A}\|_{\varepsilon} \leq \rho(\mathbf{A})+\varepsilon
  • 特别推论:若ρ(A)<1\rho (A) < 1,则有某个范数Aε<1||A||_{\varepsilon} < 1

收敛阵

  • 定义:若方阵A满足Ak0A^k \rightarrow 0,即lim(k)Ak=0lim_{(k \rightarrow \infin)} A^k=0,称A为收敛阵
  • 推论1:lim(k)Ak=0Ak0lim_{(k \rightarrow \infin)} A^k=0 \leftrightarrow ||A^k|| \rightarrow 0
  • 定理1:
    • ρ(A)<1Ak0(k)Ak0(k)\rho(A)<1 \Leftrightarrow\left\|A^{k}\right\| \rightarrow 0(k \rightarrow \infty) \Leftrightarrow A^{k} \rightarrow 0(k \rightarrow \infty)
    • 某一范数 A<1Ak0(k)\|A\|<1 \Rightarrow\left\|A^{k}\right\| \rightarrow 0(k \rightarrow \infty)
  • 牛曼公式(收敛公式):
    • ρ(A)<1\rho(A)<1, 则 I+A+A2++Ak+=(IA)1I+A+A^{2}+\cdots+A^{k}+\cdots=(I-A)^{-1}
    • 若某范数 A<1\|A\|<1, 则 I+A+A2++Ak+=(IA)1I+A+A^{2}+\cdots+A^{k}+\cdots=(I-A)^{-1}
    • ρ(A)1\rho(A) \geq 1, 则 I+A+A2++Ak+发散I+A+A^{2}+\cdots+A^{k}+\cdots发散
      发散例题

谱半径估计

  • 定义:称ACnnA \in C^{n*n},称A的n个特征值的模的最大者为A的谱半径,即为p(A)p(A)

  • 定理:设ACnnA \in C^{n*n},则ρ(A)\rho (A)不大于A的任一矩阵范数,即为ρ(A)A\rho(A) \leq ||A||

  • 特别:ρ(A)A\rho(A) \leq ||A||_{\infty},且ρ(A)A1\rho(A) \leq ||A||_{1},即ρ(A)AT\rho (A) \leq ||A^T||_\infty

  • 正矩阵定义1:一个实矩阵A=(aij)RmnA=(a_{ij}) \in R^{m*n}

    • 若每个元素aij0a_{ij} \geq 0,称A为非负矩阵,记作A0A \geq 0
    • 若每个元素aij>0a_{ij} > 0,称A为正矩阵,记作A>0A > 0
  • 正矩阵定理:设非负阵A=(aij)nn0A=(a_{ij})_{n*n} \geq 0,令h=(A的最小行和),l=(A的最小列和)h=(A的最小行和),l=(A的最小列和),则

    • hρ(A)A(A的最大行和)h \leq \rho(A) \leq ||A||_\infty(A的最大行和)
    • lρ(A)A1(A的最大列和)l \leq \rho(A) \leq ||A||_1(A的最大列和)
    • 特别,若A=Ann>0A=A_{n*n} > 0为正矩阵,且h<Al<A1h<||A||_\infty或l<||A||_1
      • h<ρ(A)<Ah < \rho(A) < ||A||_\infty
      • l<ρ(A)<A1l < \rho(A) < ||A||_1

盖尔圆估计

盖尔圆定义

  • 圆盘定理1:方阵A的全体特征根都在A的n个Ger圆的并集中
    • 定义:若A的k个Ger圆相连在一起,且与其他n-k个圆分离,称此k个圆盘为一个连通分支,简称分支
  • 圆盘定理2:设A的k个Ger圆构成一个连通分支D,则在D中恰有k个特征根(含重复),特别一个孤立圆中恰有一个根
  • 由于A与转置ATA^T有相同特征根,可用ATA^T的Ger半径代替A的半径,可得A的列(圆盘定理)
  • 注:实矩阵A的n个Ger圆中心都在x轴上
  • 复习定理
    • 实系数方程的虚根一定在共轭出现
    • 实矩阵A的虚特征根必成对出现,实矩阵A的孤立圆中恰有一个实根
  • 推论1:对方阵A,若原点Oi=1nGiO \notin \bigcup_{i=1}^{n} G_{i},则A为可逆阵
  • 推论2:若A行对角占优:ajj>Rj|a_{jj}|>R_j,则A为可逆阵
  • 推论3:若A的n个Ger圆互相分离(都是孤立圆),则A是单阵(可对角化).特别,若实矩阵A的n个Ger圆互相分离,则特征根全为实数
  • 其他推论:
    • 推论1:对ACnnA \in C^{n*n},n个盖尔圆盘G1,...GnG_1,...G_n,若原点Oi=1nGiO \notin \bigcup_{i=1}^{n} G_{i},则A为非奇异阵
    • 推论2: A=(aij)Cn×nA=\left(a_{i j}\right) \in \mathbb{C}^{n \times n}. 若A对角占优, 即 aij>j=1jinaij(i=1,2,,n)\left|a_{i j}\right|>\sum_{j=1 \atop j \neq i}^{n}\left|a_{i j}\right|(i=1,2, \cdots, n) (行对角占优)或aij>j=1jinaji(i=1,2,,n)\left|a_{i j}\right|>\sum_{j=1 \atop j \neq i}^{n}\left|a_{j i}\right|(i=1,2, \cdots, n) (列对角占优), 则A为非奇异阵

评论