LA@分块矩阵@初等变换@初等矩阵#逆矩阵计算@初等变换法

news/2024/11/27 19:27:45/

文章目录

  • LA@分块矩阵@初等变换@初等矩阵#逆矩阵计算@初等变换法
    • 分块矩阵
      • 乘法
      • 准对角矩阵
      • 转置
    • 矩阵的初等变换🎈
      • 等价矩阵
      • 行`阶梯形`矩阵🎈
        • 变换步骤
      • 行简化阶梯形矩阵@最简矩阵
      • 等价标准形矩阵
        • 变换步骤
      • 小结
    • 初等矩阵
      • 初等矩阵的结论🎈
      • 小结
    • 逆矩阵计算@初等变换法
      • 步骤
      • 利用初等行变换直接求解简单`矩阵方程`

LA@分块矩阵@初等变换@初等矩阵#逆矩阵计算@初等变换法

分块矩阵

乘法

  • 注意对相乘的两个高阶矩阵的划分方法,子矩阵之间要能够相乘

    • Am×nBn×p=Cm×pA_{m\times{n}}B_{n\times{p}}=C_{m\times{p}}Am×nBn×p=Cm×p

      • m=∑i=1tmim=\sum_{i=1}^{t}m_{i}m=i=1tmi
      • n=∑i=1lnin=\sum_{i=1}^{l}n_{i}n=i=1lni
      • p=∑i=1rpip=\sum_{i=1}^{r}p_{i}p=i=1rpi
      • A=(Aij)t×lA=(A_{ij})_{t\times{l}}A=(Aij)t×l
      • B=(Bij)l×rB=(B_{ij})_{l\times{r}}B=(Bij)l×r
      • C=((AB)ij)t×rC=((AB)_{ij})_{t\times{r}}C=((AB)ij)t×r
      • (Aij)mi×nj,(Bij)ni×pj(A_{ij})_{m_i\times{n_j}},(B_{ij})_{n_{i}\times{p_j}}(Aij)mi×nj,(Bij)ni×pj
    • Cij=∑k=1lAikBkj,(i=1,2,⋯,t;j=1,2,⋯,r)C_{ij}=\sum_{k=1}^{l}A_{ik}B_{kj},(i=1,2,\cdots,t;j=1,2,\cdots,{r})Cij=k=1lAikBkj,(i=1,2,,t;j=1,2,,r)

  • 借助分块乘法,可以将高阶矩阵求逆转为低阶矩阵求逆

  • n阶矩阵:D=(A0CB)D=\begin{pmatrix}A&\bold{0}\\C&B\end{pmatrix}D=(AC0B),A,B均可逆

    • 待定系数法:D−1=(X11X12X21X22)D^{-1}=\begin{pmatrix}X_{11}&{X_{12}}\\X_{21}&X_{22}\end{pmatrix}D1=(X11X21X12X22)

    • DD−1=(A0CB)(X11X12X21X22)=(Ek00Er)=En=k+rDD^{-1}=\begin{pmatrix}A&\bold{0}\\C&B\end{pmatrix} \begin{pmatrix}X_{11}&{X_{12}}\\X_{21}&X_{22}\end{pmatrix} =\begin{pmatrix}E_k&{\bold{0}}\\\bold{0}&E_r\end{pmatrix}=E_{n=k+r} DD1=(AC0B)(X11X21X12X22)=(Ek00Er)=En=k+r

    • AX11=EkAX12=0CX11+BX21=0CX12+BX22=Er对前两个式子两边左乘A−1,得到X11=A−1;X12=0分别带入第3,4式中,得到:X21=−B−1CA−1BX22=Er,左乘B−1,X22=B−1得:D−1=(A−10−B−1CA−1B−1)AX_{11}=E_k\\ AX_{12}=\bold{0}\\ CX_{11}+BX_{21}=0\\ CX_{12}+BX_{22}=E_r \\ 对前两个式子两边左乘A^{-1},得到X_{11}=A^{-1};X_{12}=\bold{0} \\分别带入第3,4式中,得到: \\X_{21}=-B^{-1}CA^{-1} \\BX_{22}=E_r,左乘B^{-1},X_{22}=B^{-1} \\得: \\D^{-1}=\begin{pmatrix}A^{-1}&\bold{0}\\ -B^{-1}CA^{-1}&B^{-1}\end{pmatrix} AX11=EkAX12=0CX11+BX21=0CX12+BX22=Er对前两个式子两边左乘A1,得到X11=A1;X12=0分别带入第3,4式中,得到:X21=B1CA1BX22=Er,左乘B1,X22=B1:D1=(A1B1CA10B1)

    • 特别的,当C=0C=\bold{0}C=0

      • D=(A00B)D−1=(A−100B−1)D=\begin{pmatrix}A&\bold{0}\\\bold{0}&B\end{pmatrix} \\ D^{-1}=\begin{pmatrix}A^{-1}&\bold{0}\\ \bold{0}&B^{-1}\end{pmatrix} D=(A00B)D1=(A100B1)
    • 类似的:

      • (AC0B)−1=(A−1−A−1CB−10B−1)(0AB0)−1=(0B−1A−10)\\ \begin{pmatrix}A&C\\\bold{0}&B\end{pmatrix}^{-1} =\begin{pmatrix}A^{-1}&-A^{-1}CB^{-1}\\ \bold0&B^{-1}\end{pmatrix} \\ \begin{pmatrix}\bold{0}&A\\B&\bold{0}\end{pmatrix}^{-1} =\begin{pmatrix}\bold{0}&B^{-1}\\A^{-1}&\bold{0}\end{pmatrix} (A0CB)1=(A10A1CB1B1)(0BA0)1=(0A1B10)

      • 注意区分行列式中的拉普拉斯展开

准对角矩阵

  • A=(A110⋯00A22⋯0⋮⋮⋱⋮00⋯Ann)A−1=(A11−10⋯00A22−1⋯0⋮⋮⋱⋮00⋯Ann−1)将两个矩阵相乘,容易验证A=\begin{pmatrix} A_{11}& 0& \cdots & 0 \\ 0& A_{22}& \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0& 0& \cdots & A_{nn} \end{pmatrix} \\ A^{-1}=\begin{pmatrix} A_{11}^{-1}& 0& \cdots & 0 \\ 0& A_{22}^{-1}& \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0& 0& \cdots & A_{nn}^{-1} \end{pmatrix} \\将两个矩阵相乘,容易验证 A=A11000A22000AnnA1=A111000A221000Ann1将两个矩阵相乘,容易验证

转置

  • Ap×q=(A11A12⋯A1qA21A22⋯A2q⋮⋮⋱⋮Ap1Ap2⋯Apq)Aq×pT=(A11TA21T⋯Ap1TA12TA22T⋯Ap2T⋮⋮⋱⋮A1qTA2qT⋯ApqT)A_{p\times{q}}=\begin{pmatrix} A_{11}& A_{12}& \cdots & A_{1q} \\ A_{21}& A_{22}& \cdots & A_{2q} \\ \vdots & \vdots & \ddots & \vdots \\ A_{p1}& A_{p2}& \cdots & A_{pq} \end{pmatrix} \\ A^T_{q\times{p}}=\begin{pmatrix} A_{11}^T& A_{21}^T& \cdots & A_{p1}^T \\ A_{12}^T& A_{22}^T& \cdots & A_{p2}^T \\ \vdots & \vdots & \ddots & \vdots \\ A_{1q}^T& A_{2q}^T& \cdots & A_{pq}^T \end{pmatrix} Ap×q=A11A21Ap1A12A22Ap2A1qA2qApqAq×pT=A11TA12TA1qTA21TA22TA2qTAp1TAp2TApqT

矩阵的初等变换🎈

  • 初等变换分为行变换和列变换,仅以行变换为例:

    • 交换矩阵的某两行位置ri↔rjr_i\leftrightarrow{r_j}rirj

    • 用一个非零常数k乘以矩阵的某一行kri,k≠0kr_i,k\ne{0}kri,k=0

    • 将矩阵的某一行乘以非零常数k后加到另一行:kri+rj,k≠0kr_i+r_j,k\ne{0}kri+rj,k=0

等价矩阵

  • 如果矩阵B可以由矩阵A经过有限次初等变换得到,则称A和B为等价矩阵,记为A≅BA\cong{B}AB
    • 其中≅\cong式全等(congruent)号,在此处是等价号
  • 等价矩阵的性质:
    • A≅AA\cong{A}AA
    • 若A≅B则B≅A若A\cong{B}则B\cong{A}ABBA
    • 若A≅B,B≅C则A≅C若A\cong{B},B\cong{C}则A\cong{C}AB,BCAC

阶梯形矩阵🎈

  • Row Echelon Form,记为REF-Matrix
  1. 若有零行(矩阵中元素全为0的行),则零行都位于矩阵下方(特点1)
  2. 从第一行起,每行的主元(第一个非零元素)前面零的个数逐行增加(严格增加). (特点2)
    • 或者说,每行主元的列指标严格增加

变换步骤

  • 为了得到行简化矩阵,可以执行一系列的初等变换(这些初等变换要么全是行变换,要么全是列变换)

  • 将矩阵行中最接近全零行的行调整到第一行

  • 利用第一行将第一列到最后一列尽可能的零化(使得结果满足特点1,2)

  • A=(31561−13−221351111)→r1↔r4(11111−13−221353156)→r2−r1r3−2r1r4−3r1(11110−22−30−1130−223)→r2↔r3(11110−1130−22−30−223)→r3−2r2r4−2r2(11110−113000−9000−3)→−19r3(11110−1130001000−3)→r4+3r3(11110−11300010000)A=\left( \begin{matrix} 3& 1& 5& 6\\ 1& -1& 3& -2\\ 2& 1& 3& 5\\ 1& 1& 1& 1\\ \end{matrix} \right) \xrightarrow{r_1\leftrightarrow r_4}\left( \begin{matrix} 1& 1& 1& 1\\ 1& -1& 3& -2\\ 2& 1& 3& 5\\ 3& 1& 5& 6\\ \end{matrix} \right) \xrightarrow{\begin{array}{c} r_2-r_1\\ r_3-2r_1\\ r_4-3r_1\\ \end{array}} \\ \left( \begin{matrix} 1& 1& 1& 1\\ 0& -2& 2& -3\\ 0& -1& 1& 3\\ 0& -2& 2& 3\\ \end{matrix} \right) \xrightarrow{r_2\leftrightarrow r_3}\left( \begin{matrix} 1& 1& 1& 1\\ 0& -1& 1& 3\\ 0& -2& 2& -3\\ 0& -2& 2& 3\\ \end{matrix} \right) \xrightarrow{\begin{array}{c} r_3-2r_2\\ r_4-2r_2\\ \end{array}} \\ \left( \begin{matrix} 1& 1& 1& 1\\ 0& -1& 1& 3\\ 0& 0& 0& -9\\ 0& 0& 0& -3\\ \end{matrix} \right) \xrightarrow{-\frac{1}{9}r_3}\left( \begin{matrix} 1& 1& 1& 1\\ 0& -1& 1& 3\\ 0& 0& 0& 1\\ 0& 0& 0& -3\\ \end{matrix} \right) \xrightarrow{r_4+3r_3}\left( \begin{matrix} 1& 1& 1& 1\\ 0& -1& 1& 3\\ 0& 0& 0& 1\\ 0& 0& 0& 0\\ \end{matrix} \right) A=3121111153316251r1r41123111113351256r2r1r32r1r43r11000121212121333r2r31000112211221333r32r2r42r2100011001100139391r31000110011001313r4+3r31000110011001310

  • 行阶梯形矩阵的条件比较宽松,同一个矩阵化为行阶梯型矩阵的结果可以是多种多样的

    • 从一般行阶梯形矩阵化为行简化矩阵(最简矩阵)的过程(中间结果矩阵)都是行阶梯形矩阵
  • 但是最简形式是一样的(行简化阶梯形矩阵)

行简化阶梯形矩阵@最简矩阵

  • 简化列阶梯形矩阵简约行梯形式矩阵(reduced row echelon form),也称作行规范形矩阵(row canonical form),记为RREF-Matrix如果满足额外的条件:

    • 阶梯形矩阵中非零行主元(第一个非零元素)为1

    • 且主元所在的的其他元素都为0

  • (11110−11300010000)→−r2(111101−1−300010000)→r1−r2(102401−1−300010000)→r2+3r3;r1−4r3(102001−1000010000)\left( \begin{matrix} 1& 1& 1& 1\\ 0& -1& 1& 3\\ 0& 0& 0& 1\\ 0& 0& 0& 0\\ \end{matrix} \right) \xrightarrow{-r_2} \left( \begin{matrix} 1& 1& 1& 1\\ 0& 1& -1& -3\\ 0& 0& 0& 1\\ 0& 0& 0& 0\\ \end{matrix} \right) \xrightarrow{r_1-r_2} \left( \begin{matrix} 1& 0& 2& 4\\ 0& 1& -1& -3\\ 0& 0& 0& 1\\ 0& 0& 0& 0\\ \end{matrix} \right) \xrightarrow{r_2+3r_3;r_1-4r_3} \\ \left( \begin{matrix} 1& 0& 2& 0\\ 0& 1& -1& 0\\ 0& 0& 0& 1\\ 0& 0& 0& 0\\ \end{matrix} \right) 1000110011001310r21000110011001310r1r21000010021004310r2+3r3;r14r31000010021000010

等价标准形矩阵

  • Matrix equivalence - Wikipedia

  • Canonical form

  • The rank property yields an intuitive canonical form for matrices of the equivalence class of rank k as

    • (100⋯0010⋯000⋱0⋮1⋮0⋱0⋯0){\displaystyle {\begin{pmatrix}1&0&0&&\cdots &&0\\0&1&0&&\cdots &&0\\0&0&\ddots &&&&0\\\vdots &&&1&&&\vdots \\&&&&0&&\\&&&&&\ddots &\\0&&&\cdots &&&0\end{pmatrix} } } 100001000100000
  • where the number of 1s on the diagonal is equal to k. This is a special case of the Smith normal form, which generalizes this concept on vector spaces to free modules over principal ideal domains.

  • 设A为n阶矩阵

    • 其等价矩阵(行简化阶梯形矩阵)中,位于左上角的子矩阵是一个r阶单位阵(r⩽nr\leqslant{}nrn)

    • 其余子块都是零矩阵(如果有的话)

    • note:当A是可逆方阵时,r=nr=nr=n

    • C=(Er0r×s0s×r0s)r+s=n(r⩽n),当A可逆时取等号C= \begin{pmatrix} E_r&0_{r\times{s}}\\ 0_{s\times{r}}&0_{s} \end{pmatrix} \\ r+s=n(r\leqslant{n}),当A可逆时取等号 C=(Er0s×r0r×s0s)r+s=n(rn),A可逆时取等号

    • 对A经过一系列的初等变换得到标准形矩阵C,可以用初等矩阵乘法表示为PAQ=CPAQ=CPAQ=C

      • 对两边取行列式:∣P∣∣A∣∣Q∣=∣C∣|P||A||Q|=|C|P∣∣A∣∣Q=C

        • 由拉普拉斯展开

          • ∣AmROBn∣=∣AmORBn∣=∣Am∣⋅∣Bn∣\begin{vmatrix} A_m&R \\ O&B_n \end{vmatrix} = \begin{vmatrix} A_m& O\\ R&B_n \end{vmatrix} =|A_m|\cdot|B_n| AmORBn=AmROBn=AmBn
        • ∣P∣∣A∣∣Q∣=∣Er∣∣0s∣|P||A||Q|=|E_r||0_s|P∣∣A∣∣Q=Er∣∣0s

      • 若A可逆,则∣A∣≠0|A|\neq{0}A=0

        • 初等矩阵P,QP,QP,Q均可逆(∣P∣∣Q∣|P||Q|P∣∣Q不为0),等式∣P∣∣A∣∣Q∣=∣Er∣∣0s∣|P||A||Q|=|E_r||0_s|P∣∣A∣∣Q=Er∣∣0s左边不为0,所以等式右边也不为0

        • 可见,只有当r=n时,等式右边不为0,此时C是n阶单位矩阵

变换步骤

  • 假设矩阵A经过初等行变换得行简化阶梯形矩阵B=RREF(A)

  • 对B再所若干次初等列变换,得到A的等价标准型矩阵
    (102001−1000010000)→r3−2r1(100001−1000010000)→r3+r2(1000010000010000)→r3↔r4(1000010000100000)\left( \begin{matrix} 1& 0& 2& 0\\ 0& 1& -1& 0\\ 0& 0& 0& 1\\ 0& 0& 0& 0\\ \end{matrix} \right) \xrightarrow{r_3-2r_1} \left( \begin{matrix} 1& 0& 0& 0\\ 0& 1& -1& 0\\ 0& 0& 0& 1\\ 0& 0& 0& 0\\ \end{matrix} \right) \xrightarrow{r_3+r_2} \left( \begin{matrix} 1& 0& 0& 0\\ 0& 1& 0& 0\\ 0& 0& 0& 1\\ 0& 0& 0& 0\\ \end{matrix} \right) \\ \xrightarrow{r_3\leftrightarrow{r_4}} \left( \begin{matrix} 1& 0& 0& 0\\ 0& 1& 0& 0\\ 0& 0& 1& 0\\ 0& 0& 0& 0\\ \end{matrix} \right) 1000010021000010r32r11000010001000010r3+r21000010000000010r3r41000010000100000

  • 将过上述变换,将一个一般的方阵A变换为一个与之等价的标准形矩阵

小结

  • 任意矩阵都可以通过若干次初等变换化为"行阶梯形矩阵"和"行简化阶梯形矩阵"
  • 任意矩阵都可以通过初等变换(包括行变换和列变换)化为等价的标准形矩阵

初等矩阵

  • 单位矩阵E做一次初等行变换得到的矩阵称为初等矩阵

    • 易知单位矩阵本身也是初等矩阵

    • 每个初等变换都有一个与之对应的初等矩阵

    • 交换第i,j两行(列)的位置,记为P(i,j)P(i,j)P(i,j)

      • 交换第i,j行和交换第i,j列效果一样

      • (1⋱1(0)[1]1⋱1(1)[0]1⋱1)\begin{pmatrix} 1& & & & & & & & & & \\ & \ddots& & & & & & & & & \\ & & 1& & & & & & & & \\ & & & (0)& & & & [1]& & & \\ & & & & 1& & & & & & \\ & & & & & \ddots& & & & & \\ & & & & & & 1& & & & \\ & & & (1)& & & & [0]& & & \\ & & & & & & & & 1& & \\ & & & & & & & & & \ddots& \\ & & & & & & & & & & 1\\ \end{pmatrix} 11(0)(1)11[1][0]11

    • k倍乘第i行(列),记为P(i(k))P(i(k))P(i(k))

      • 第i行乘以k和第i列乘以k效果一样

      • (1⋱1k1⋱1)\begin{pmatrix} 1& & & & & & \\ & \ddots& & & & & \\ & & 1& & & & \\ & & & k& & & \\ & & & & 1& & \\ & & & & & \ddots& \\ & & & & & & 1 \\ \end{pmatrix} 11k11

    • k倍乘第j行,加到第i行,记为P(i,j(k))P(i,j(k))P(i,j(k))

      • 第j行的k倍加到第i行(ri+krj)(r_i+kr_j)(ri+krj)效果等价于第i列乘以k倍加到第j列(cj+kci)(c_j+kc_i)(cj+kci)

      • i,ji,{j}i,j相等,此时变为倍乘效果

      • (1⋱1⋯k⋱⋮1⋱1)\begin{pmatrix} 1 & & & & & & \\ & \ddots & & & & & \\ & & 1 & \cdots & k & & \\ & & & \ddots & \vdots & & \\ & & & & 1 & & \\ & & & & & \ddots & \\ & & & & & & 1 \\ \end{pmatrix} 11k11

初等矩阵的结论🎈

  • 初等矩阵有三类:P(i,j),P(i(k)),P(i,j(k))P(i,j),P(i(k)),P(i,j(k))P(i,j),P(i(k)),P(i,j(k))

  • 初等矩阵均可逆,其逆矩阵仍然为初等矩阵(用逆矩阵定义可以求证下列结论)

    • 对初等矩阵取行列式,都不为0
    • P(i,j)−1=P(i,j)P(i,j)^{-1}=P(i,j)P(i,j)1=P(i,j)
    • P(i(k))−1=P(i(k−1))P(i(k))^{-1}=P(i(k^{-1}))P(i(k))1=P(i(k1))
    • P(i,j(k))−1=P(i,j(−k))P(i,j(k))^{-1}=P(i,j(-k))P(i,j(k))1=P(i,j(k))
  • 对矩阵Am×nA_{m\times{n}}Am×n

    • 做一次初等行变换,相当于对A左乘相应的mmm阶初等矩阵
    • 做一次初等列变换,相当于对A右乘相应的nnn阶初等矩阵.
    • 初等矩阵是由单位矩阵经过初等变换得到的.
      • 假设初等矩阵M是由单位阵E经过初等变换P(f)P(f)P(f)得到的E→P(f)ME\xrightarrow{P(f)}MEP(f)M
        • P(f)P(f)P(f)变换可细分为
          • 初等行变换Pr(f)P_r(f)Pr(f)
          • 初等列变换,Pc(f)P_c(f)Pc(f)
          • 能够满足E→ME\to{M}EM的变换即可
        • 可以初等矩阵M详细记为M=MP(f)M=M_{P(f)}M=MP(f)
      • 那么
        • C=MAC=MAC=MA相当于:对矩阵A做一次Pr(f)P_r(f)Pr(f)变换得到矩阵C
        • D=AMD=AMD=AM相当于:对矩阵A做一次Pc(f)P_c(f)Pc(f)变换得到矩阵D
  • 证明:(以初等行变换)为例

    • 做如下矩阵分块(预计算),相比直接计算原始矩阵更加简洁

    • 一般的,设m阶矩阵Bm×m=(b11b12⋯b14b21b22⋯b24⋮⋮⋮b41b42⋯b44)对矩阵Am×n进行分块(按行分块),Ai为第i行(i=1,2,⋯m)由矩阵的分块乘法得到:Cm×n=BA=(∑k=1mb1kAk∑k=1mb2kAk⋮∑k=1mbmkAk)其中∑k=1mbikAk和Ak都是1×n的矩阵(行向量)一般的,设m阶矩阵 B_{m\times{m}}=\begin{pmatrix} b_{11} &b_{12} &\cdots &b_{14} \\ b_{21} &b_{22} &\cdots &b_{24} \\ \vdots &\vdots & &\vdots \\ b_{41} &b_{42} &\cdots &b_{44} \\ \end{pmatrix} \\ 对矩阵A_{m\times{n}}进行分块(按行分块),A_i为第i行(i=1,2,\cdots{m}) \\由矩阵的分块乘法得到: \\C_{m\times{n}}=BA=\begin{pmatrix} \sum\limits_{k=1}^{m}b_{1k}A_{k} \\ \sum\limits_{k=1}^{m}b_{2k}A_{k} \\ \vdots \\ \sum\limits_{k=1}^{m}b_{mk}A_{k} \\ \end{pmatrix} \\其中\sum\limits_{k=1}^{m}b_{ik}A_{k}和A_k都是1\times{n}的矩阵(行向量) 一般的,m阶矩阵Bm×m=b11b21b41b12b22b42b14b24b44对矩阵Am×n进行分块(按行分块),Ai为第i(i=1,2,m)由矩阵的分块乘法得到:Cm×n=BA=k=1mb1kAkk=1mb2kAkk=1mbmkAk其中k=1mbikAkAk都是1×n的矩阵(行向量)

    • B=P(i,j)B=P(i,j)B=P(i,j)
      C=BA=P(i,j)A=(A1⋮Aj⋮Ai⋮Am)C=BA=P(i,j)A=\begin{pmatrix} A_{1} \\ \vdots \\ \boxed{A_{j}} \\ \vdots \\ \boxed{A_{i}} \\ \vdots \\ A_{m} \\ \end{pmatrix} C=BA=P(i,j)A=A1AjAiAm

      • 这相当于把矩阵A的第i行和第j行互换
    • 类似的,令B=P(i(k))B=P(i(k))B=P(i(k))有:

      • C=BA=P(i(k))A=(A1⋮kAi⋮Am)C=BA=P(i(k))A=\begin{pmatrix} A_{1} \\ \vdots \\ \boxed{kA_{i}} \\ \vdots \\ A_{m} \\ \end{pmatrix} C=BA=P(i(k))A=A1kAiAm

      • 这相当于把矩阵A的第i行乘以k倍

    • B=P(i,j(k))B=P(i,j(k))B=P(i,j(k))

      • C=BA=P(i,j(k))A=(A1⋮Ai+kAj⋮Aj⋮Am)C=BA=P(i,j(k))A=\begin{pmatrix} A_{1} \\ \vdots \\ \boxed{A_{i}}+kA_j \\ \vdots \\ \boxed{A_{j}} \\ \vdots \\ A_{m} \\ \end{pmatrix} C=BA=P(i,j(k))A=A1Ai+kAjAjAm

      • 相当于将矩阵A第j行乘以k倍加到第i行

    • 右乘初等矩阵和列变换关系的证明

      • 设Am×n,Bn×n对A,B进行分块:A=(A1A2⋮Am)=(α1α2⋯αn)其中Ai(i=1,⋯,m)是1×n的行向量αj(j=1,⋯,n)是m×1的列向量B=(B1B2⋯Bn)其中Bi(i=1,⋯,n)是n×1的列向量C=AB=(A1A2⋮Am)(B1B2⋯Bn)=(A1B1A1B2⋯A1B4A2B1A2B2⋯A2B4⋮⋮⋮A4B1A4B2⋯A4B4)Cij=AiBj是1×1的矩阵设A_{m\times{n}},B_{n\times{n}} \\对A,B进行分块: \\A=\begin{pmatrix} A_{1} \\ A_{2} \\ \vdots \\ A_{m} \\ \end{pmatrix} =\begin{pmatrix} \alpha_{1}&\alpha_{2}&\cdots&\alpha_{n} \end{pmatrix} \\其中A_i(i=1,\cdots,m)是1\times{n}的行向量 \\ \alpha_{j}(j=1,\cdots,n)是m\times{1}的列向量 \\\\ B=\begin{pmatrix} B_{1}&B_{2}&\cdots&B_{n} \end{pmatrix} \\其中B_i(i=1,\cdots,n)是n\times{1}的列向量 \\ C=AB=\begin{pmatrix} A_{1} \\ A_{2} \\ \vdots \\ A_{m} \\ \end{pmatrix} \begin{pmatrix} B_{1}&B_{2}&\cdots&B_{n} \end{pmatrix} =\begin{pmatrix} A_{1}B_{1}&A_{1}B_{2}&\cdots &A_{1}B_{4} \\ A_{2}B_{1}&A_{2}B_{2}&\cdots &A_{2}B_{4} \\ \vdots &\vdots & &\vdots \\ A_{4}B_{1}&A_{4}B_{2}&\cdots &A_{4}B_{4} \\ \end{pmatrix} \\C_{ij}=A_{i}B_{j}是1\times{1}的矩阵 Am×n,Bn×nA,B进行分块:A=A1A2Am=(α1α2αn)其中Ai(i=1,,m)1×n的行向量αj(j=1,,n)m×1的列向量B=(B1B2Bn)其中Bi(i=1,,n)n×1的列向量C=AB=A1A2Am(B1B2Bn)=A1B1A2B1A4B1A1B2A2B2A4B2A1B4A2B4A4B4Cij=AiBj1×1的矩阵

      • B=P(i,j)B=P(i,j)B=P(i,j)

        • C=AB=(a11⋯aj⋯ai⋯a1n⋮⋮⋮⋮am1⋯aj⋯ai⋯amn)=(α1⋯αj⋯αi⋯αn)C=AB =\begin{pmatrix} a_{11} &\cdots &a_{j} &\cdots &a_{i} &\cdots &a_{1n} \\ \vdots & &\vdots & &\vdots & &\vdots \\ a_{m1} &\cdots &a_{j} &\cdots &a_{i} &\cdots &a_{mn} \\ \end{pmatrix}\\ =\begin{pmatrix} \alpha_{1}&\cdots&\alpha_{j}&\cdots&\alpha_{i}&\cdots&\alpha_{n} \end{pmatrix} C=AB=a11am1ajajaiaia1namn=(α1αjαiαn)

        • 相当于将矩阵A的第i列和第j列互换

      • B=P(i(k))B=P(i(k))B=P(i(k))

        • C=AB=(a11⋯kai⋯a1n⋮⋮⋮am1⋯kai⋯amn)=(α1⋯kαi⋯αn)C=AB=\begin{pmatrix} a_{11} &\cdots &ka_{i} &\cdots &a_{1n} \\ \vdots & &\vdots & &\vdots \\ a_{m1} &\cdots &ka_{i} &\cdots &a_{mn} \end{pmatrix} \\ =\begin{pmatrix} \alpha_{1}&\cdots &k\alpha_{i}&\cdots&\alpha_{n} \end{pmatrix} C=AB=a11am1kaikaia1namn=(α1kαiαn)

        • 相当于将矩阵A的第i列乘以k

      • B=P(i,j(k))B=P(i,j(k))B=P(i,j(k))

        • C=AB=(a11⋯ai+kaj⋯a1n⋮⋮⋮am1⋯ai+kaj⋯amn)=(α1⋯αi+kaj⋯αn)C=AB=\begin{pmatrix} a_{11} &\cdots &a_i+ka_{j} &\cdots &a_{1n} \\ \vdots & &\vdots & &\vdots \\ a_{m1} &\cdots &a_i+ka_{j} &\cdots &a_{mn} \end{pmatrix} \\ =\begin{pmatrix} \alpha_{1}&\cdots &\alpha_{i}+ka_j&\cdots&\alpha_{n} \end{pmatrix} C=AB=a11am1ai+kajai+kaja1namn=(α1αi+kajαn)

        • 相当于将矩阵A的第j列乘以k加到第i列

小结

  • 矩阵A和B等价的充分条件是存在初等矩阵P1,P2,⋯,PsP_1,P_2,\cdots,P_{s}P1,P2,,PsQ1,Q2,⋯,QtQ_1,Q_2,\cdots,Q_tQ1,Q2,,Qt使得

    • B=Ps⋯P2P1AQ1Q2⋯QtB=P_s\cdots{P_2P_1}AQ_{1}Q_{2}\cdots{Q_t}B=PsP2P1AQ1Q2Qt
  • 可逆矩阵和同阶单位矩阵等价

    • 任意矩阵可以化为(存在)与之等价的标准形矩阵:

      • PAQ=EPAQ=EPAQ=E

逆矩阵计算@初等变换法

  • n阶矩阵A可逆的充要条件是,A可以表示为若干个初等矩阵的乘积
    • 因为n阶矩阵A可逆(A与同阶单位阵E等价),则存在n阶初等矩阵序列P=P1P2⋯PsP=P_1P_2\cdots P_sP=P1P2PsQ=Q1Q2⋯QtQ=Q_1Q_2\cdots Q_tQ=Q1Q2Qt,使得
      • PAQ=EPAQ=EPAQ=E
      • 由于初等矩阵及其逆矩阵都是可逆矩阵(Pi,Qj(i=1,2⋯,s;j=1,2,⋯,t)P_i,Q_j(i=1,2\cdots,s;j=1,2,\cdots,t)Pi,Qj(i=1,2,s;j=1,2,,t))都可逆
      • 两侧同时左乘P−1=Ps−1⋯P2−1P1−1P^{-1}=P_s^{-1}\cdots{}{P_{2}^{-1}P_1^{-1}}P1=Ps1P21P11Q=Qt−1⋯Q2−1Q1−1Q=Q_{t}^{-1}\cdots{Q_{2}^{-1}Q_{1}^{-1}}Q=Qt1Q21Q11
      • A=PEQ=PQA=PEQ=PQA=PEQ=PQ
      • 为例更加统一,可以记W=PQ=P1P2⋯PsQ1Q2⋯Qt=W1W2⋯Wz,(z=s+t)W=PQ=P_1P_2\cdots P_sQ_1Q_2\cdots Q_t=W_1W_2\cdots W_{z},(z=s+t)W=PQ=P1P2PsQ1Q2Qt=W1W2Wz,(z=s+t)
      • 所以可逆矩阵A可以表示为若干个初等矩阵的乘积
  • 任意可逆矩阵都可以经过初等行变换化为单位矩阵
    • 设A为可逆矩阵
    • A可以表示为若干个初等矩阵的乘积,记为A=W1W2⋯WzA=W_1W_2\cdots W_{z}A=W1W2Wz
    • A−1=Ws−1⋯W2−1W1−1A^{-1}=W_s^{-1}\cdots{W_{2}^{-1}W_{1}^{-1}}A1=Ws1W21W11
    • 从而A−1A=Ws−1⋯W2−1W1−1A=EA^{-1}A=W_s^{-1}\cdots{W_{2}^{-1}W_{1}^{-1}}A=EA1A=Ws1W21W11A=E
    • 即可逆矩阵A左乘初等矩阵序列 Ws−1⋯W2−1W1−1W_s^{-1}\cdots{W_{2}^{-1}W_{1}^{-1}}Ws1W21W11得到单位矩阵
    • 对A做若干初等行变换(而不需要初等列变换)就可求出A−1A^{-1}A1
    • 利用这个原理可以求解可逆矩阵的逆

步骤

  • 对n阶可逆矩阵A,构造一个n×2nn\times{2n}n×2n的增广矩阵B=(A∣E)n×2nB=(A|E)_{n\times{2n}}B=(AE)n×2n

  • A可以表示为A=P1P2⋯PsA=P_{1}P_2\cdots{P_s}A=P1P2Ps

  • 易知对B=(A∣E)B=(A|E)B=(AE)左乘A−1A^{-1}A1,得到A−1B=(A−1A∣A−1E)=(E∣A−1)A^{-1}B=(A^{-1}A|A^{-1}E)=(E|A^{-1})A1B=(A1AA1E)=(EA1)

    • 由分块矩阵乘法公式可以得到上述结论(a)(bc)=(abac)(a)(b\ c)=(ab\ ac)(a)(b c)=(ab ac)

    • A−1A^{-1}A1是我们需要求解的东西,往往无法直接得到

    • 但是A−1=Ps−1⋯P2−1P1−1A^{-1}=P_s^{-1}\cdots{P_2^{-1}P_1^{-1}}A1=Ps1P21P11

      • Ps−1⋯P2−1P1−1B=(E∣A−1)P_s^{-1}\cdots{P_2^{-1}P_1^{-1}}B=(E|A^{-1})Ps1P21P11B=(EA1)
      • 也就是说通过观察矩阵A(或其增广阵B),对其进行一系列初等行变换来得到包含A−1A^{-1}A1的矩阵(E∣A−1)(E|A^{-1})(EA1),从该矩阵读出A−1A^{-1}A1

利用初等行变换直接求解简单矩阵方程

  • 对于基本的矩阵方程AX=BAX=BAX=B
    • A为n阶可逆方阵
    • B为n×mn\times{m}n×m矩阵
    • X则为n×mn\times{m}n×m矩阵
  • AX=BAX=BAX=B两边同时左乘A−1A^{-1}A1
    • X=A−1BX=A^{-1}BX=A1B
  • 构造C=(A∣B)C=(A|B)C=(AB)
    • 对两边同时左乘A−1A^{-1}A1,A−1C=A−1(A∣B)=(E∣A−1B)A^{-1}C=A^{-1}(A|B)=(E|A^{-1}B)A1C=A1(AB)=(EA1B)
    • 我们从(E∣A−1B)(E|A^{-1}B)(EA1B)可以直接读出X=A−1BX=A^{-1}BX=A1B(但是读不出A−1A^{-1}A1)

http://www.ppmy.cn/news/19048.html

相关文章

学习记录669@项目管理之项目合同管理

有效合同原则 有效合同应具备以下特点: (1)签订合同的当事人应当具有相应的民事权利能力和民事行为能力。 (2)意思表示真实。 (3)不违反法律或社会公共利益 与有效合同相对应,需要避免无效合同。无效合同通常需具备下列任一情形: (1)一方以欺诈、胁迫的手段订立合…

Flutter 基础-下

一、shared_preferences shared_preferences 是一个本地数据缓存库(类似 AsyncStorage) https://pub.dev/packages/shared_preferences 使用步骤 在 pubsepc.yaml 中添加 shared_preferences 依赖安装依赖(pub get | flutter packages get |…

JDBC快速入门,如何使用JDBC操作数据库?

文章目录1. 前言2. JDBC 概述2.1 概念2.2 优点3. JDBC 快速入门Java编程基础教程系列1. 前言 在 Java 开发中,使用 Java 语言操作数据库是非常重要的一部分,那么 Java 语言是如何操作数据库的呢?我们需要使用不同厂商的数据库时,…

sbt编程语言scala的构建工具配置及项目构建(附带网盘下载)

SBT简介 SBT 是 Scala 的构建工具,全称是 Simple Build Tool, 类似 Maven 或 Gradle。 Java可以用Maven快速构建项目,scala用SBT快速构建一个Scala项目。 sbt下载官网 百度网盘链接:https://pan.baidu.com/s/1eJkdWndZ0izcd3w…

Yolo系列理论

参考文章:AI菌的YOLO系列 目标检测-Yolo系列发展Anchors Base原理Anchors Free原理YOLO v1Yolov1网络结构Yolov1实现方法Yolov1损失函数Yolov1总结YOLO v2Yolov2网络结构Yolov2改进细节Yolov2总结YOLO V3Yolov3网络结构YOLO V4Yolov4网络结构YOLOv4 BackBone训练策略…

开篇点睛——Elasticsearch

在互联网当中我们的查询的信息主要包括文章、视频、图片、网站信息等各式各样的复杂海量信息。怎么才能快速、准确的检索到我们想要的信息呢? 传统意义上根据数据的格式,我们会将数据分为三个大类 结构化数据非结构化数据半结构化数据 接下来我们详细的了解一下这…

【目标检测论文解读复现NO.25】基于改进Yolov5的地铁隧道附属设施与衬砌表观病害检测方法

前言 此前出了目标改进算法专栏,但是对于应用于什么场景,需要什么改进方法对应与自己的应用场景有效果,并且多少改进点能发什么水平的文章,为解决大家的困惑,此系列文章旨在给大家解读最新目标检测算法论文&#xff0c…

【 uniapp - 黑马优购 | 登录与支付(2)】如何实现三秒后跳转和微信支付

个人名片: 🐼作者简介:一名大二在校生,讨厌编程🎋 🐻‍❄️个人主页🥇:小新爱学习. 🐼个人WeChat:见文末 🕊️系列专栏:🖼…