《矩阵论》学习笔记 第0章 复习与引申

本文详细介绍了矩阵论的基本概念,包括矩阵乘法的性质,如非零因子、不可交换性、乘法消去律不成立以及代数恒等式。此外,还探讨了分块矩阵的操作,如上三角矩阵的乘积保持上三角特性,并证明了非齐次线性方程组的解的存在性。最后,讨论了向量组的极大无关组、矩阵的秩及其在满秩分解中的应用。

《矩阵论》学习笔记。学习视频为东南大学周建华教授的矩阵论课程,B站上可以搜到。

一、 矩阵乘法中应注意的问题

1. 矩阵乘法中存在非零因子

例1.
已知N=[01010⋱1⋱0]N=\left[\begin{matrix} 0 & 1 \\ & 0 & 1 \\ & & 0 & \ddots & 1 \\ & & & \ddots & 0 \end{matrix} \right]N=0101010NsN^sNs
解:Ns=[00⋯000⋯0⋮⋮⋱⋮00⋯0]N^s=\left[\begin{matrix} 0 & 0 & \cdots & 0 \\ 0 & 0 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & 0 \end{matrix} \right]Ns=000000000

2. 矩阵乘法不可交换

例2.
A=[d100d2]A=\left[\begin{matrix} d_1 & 0\\ 0 & d_2 \end{matrix} \right]A=[d100d2]其中,d1d_1d1d2d_2d2互异,求所有与AAA可交换的矩阵。
解:设B=[b11b12b21b22]B=\left[\begin{matrix} b_{11} & b_{12}\\ b_{21} & b_{22} \end{matrix} \right]B=[b11b21b12b22]
AB=[d1b11d1b12d2b21d2b22]AB=\left[\begin{matrix} d_1b_{11} & d_1b_{12}\\ d_2b_{21} & d_2b_{22} \end{matrix} \right]AB=[d1b11d2b21d1b12d2b22]
AB=BAAB=BAAB=BA ⇓\Downarrow
{d1b12=b12d2d2b21=b21d1\begin{cases} d_1b_{12}=b_{12}d_2 \\ d_2b_{21}=b_{21}d_1 \end{cases}{d1b12=b12d2d2b21=b21d1
⇓\Downarrow
b12=b21=0b_{12}=b_{21}=0b12=b21=0
只有对角矩阵相乘才可以交换。

3. 乘法消去律不成立

对给定的矩阵AAA,当AAA满足什么条件时,由AB=ACAB=ACAB=AC必可推出B=CB=CB=C?
AAA可逆。

4. 一些代数恒等式不再成立

如:
(A+B)(B+A)≠A2−B2(A+B)(B+A)\ne A^2-B^2(A+B)(B+A)=A2B2
AB=BAAB=BAAB=BA时,相应的二项式定理成立,即,
(A+B)m=Am+Cm1Am−1B+Cm2Am−2B2+⋯+Cmm−1ABm−1+Bm(A+B)^m=A^m+C_{m}^{1}A^{m-1}B+C_{m}^{2}A^{m-2}B^2+\cdots+C_{m}^{m-1}AB^{m-1}+B^m(A+B)m=Am+Cm1Am1B+Cm2Am2B2++Cmm1ABm1+Bm
例3.
计算下述矩阵的k次幂
A=[λ10⋯00λ1⋯0⋮⋮⋮⋱⋮000⋯1000⋯λ]A=\left[\begin{matrix} \lambda & 1 & 0 & \cdots & 0 \\ 0 & \lambda & 1 & \cdots & 0 \\ \vdots & \vdots & \vdots & \ddots & \vdots\\0 & 0 & 0 & \cdots & 1\\0&0&0&\cdots&\lambda \end{matrix} \right]A=λ0001λ000100001λ
解:A=λI+NA=\lambda I+NA=λI+Nλ\lambdaλ IIINNN可交换
Ak=(λI+N)k=(λI)k+Ck1(λI)k−1N+Ck2(λI)k−2N2+⋯+Ckk−1(λI)Nk−1+Nk=[λkCk1λk−1Ck2λk−2⋯Ckn−1λk−n+10λkCk1λk−1⋯Ckn−2λk−n+2⋮⋮⋮⋱⋮000⋯Ck2λk−2000⋯λk]A^k=(\lambda I+N)^k=(\lambda I)^k+C_{k}^{1}(\lambda I)^{k-1}N+C_{k}^{2}(\lambda I)^{k-2}N^2+\cdots+C_{k}^{k-1}(\lambda I)N^{k-1}+N^k \\ =\left[\begin{matrix} {\lambda}^k & C_{k}^{1}{\lambda}^{k-1} & C_{k}^{2}{\lambda}^{k-2} & \cdots & C_{k}^{n-1}{\lambda}^{k-n+1} \\ 0 & {\lambda}^k & C_{k}^{1}{\lambda}^{k-1} & \cdots & C_{k}^{n-2}{\lambda}^{k-n+2} \\ \vdots & \vdots & \vdots & \ddots & \vdots\\0 & 0 & 0 & \cdots & C_{k}^{2}{\lambda}^{k-2}\\0&0&0&\cdots&{\lambda}^k \end{matrix} \right]Ak=(λI+N)k=(λI)k+Ck1(λI)k1N+Ck2(λI)k2N2++Ckk1(λI)Nk1+Nk=λk000Ck1λk1λk00Ck2λk2Ck1λk100Ckn1λkn+1Ckn2λkn+2Ck2λk2λk
#二、分块矩阵
Aij=(aij)s×n,B=(bij)s×nA_{ij}=(a_{ij})_{s\times n},B=(b_{ij})_{s\times n}Aij=(aij)s×n,B=(bij)s×n

1. AAA,BBB均按按行进行分块

A=[α1α2⋮αs]B=[β1β2⋮βs]A=\left[\begin{matrix} \alpha_1 \\ \alpha_2 \\ \vdots \\ \alpha_s \end{matrix} \right] B=\left[\begin{matrix} \beta_1 \\ \beta_2 \\ \vdots \\ \beta_s \end{matrix} \right]A=α1α2αsB=β1β2βs
A+B=[α1+β1α2+β2⋮αs+βs]A+B=\left[\begin{matrix} \alpha_1+\beta_1 \\ \alpha_2+\beta_2 \\ \vdots \\ \alpha_s+\beta_s \end{matrix} \right]A+B=α1+β1α2+β2αs+βs
r(A+B)≤r(A)+r(B)r(A+B) \le r(A)+r(B)r(A+B)r(A)+r(B)

2. AAA按行分块,BBB不分块

AB=(α1,α2,⋯ ,αn)[b11b12⋯b1tb21b22⋯b2t⋮⋮⋱⋮bn1bn2⋯bnt]=(∑k=1nbk1αk,∑k=1nbk2αk,⋯ ,∑k=1nbktαk)AB=(\alpha_1,\alpha_2,\cdots,\alpha_n)\left[\begin{matrix}b_{11} & b_{12} & \cdots & b_{1t}\\b_{21} & b_{22} & \cdots & b_{2t}\\ \vdots & \vdots & \ddots & \vdots \\ b_{n1} & b_{n2} & \cdots & b_{nt}\end{matrix}\right] \\= (\sum_{k=1}^{n}b_{k1}\alpha_k,\sum_{k=1}^{n}b_{k2}\alpha_k,\cdots,\sum_{k=1}^{n}b_{kt}\alpha_k)AB=(α1,α2,,αn)b11b21bn1b12b22bn2b1tb2tbnt=(k=1nbk1αk,k=1nbk2αk,,k=1nbktαk)
从上式可知,ABABAB的列向量是AAA的列向量的线性组合,且组合系数刚好是BBB的相应的各个列的元素。因此
r(AB)≤r(A),r(B)r(AB)\le r(A),r(B)r(AB)r(A),r(B)

3. AAA不分块,BBB按行分块

AB=[a11a12⋯a1na21a22⋯a2n⋮⋮⋱⋮as1as2⋯asn][β1β2⋮βn]=[∑k=1na1kβk∑k=1na2kβk⋮∑k=1naskβk]AB=\left[\begin{matrix}a_{11}&a_{12}&\cdots&a_{1n}\\a_{21}&a_{22}&\cdots&a_{2n}\\ \vdots&\vdots&\ddots&\vdots\\a_{s1}&a_{s2}&\cdots&a_{sn}\end{matrix}\right]\left[\begin{matrix}\beta_1 \\ \beta_2 \\ \vdots \\ \beta_n\end{matrix}\right]\\=\left[\begin{matrix}\sum_{k=1}^{n}a_{1k}\beta_k\\\sum_{k=1}^{n}a_{2k}\beta_k\\\vdots\\\sum_{k=1}^{n}a_{sk}\beta_k\end{matrix}\right]AB=a11a21as1a12a22as2a1na2nasnβ1β2βn=k=1na1kβkk=1na2kβkk=1naskβk
从上式可知,ABABAB的行向量是BBB的行向量的线性组合,且组合系数刚好是AAA的相应的各个行的元素。因此
r(AB)≤r(A),r(B)r(AB)\le r(A),r(B)r(AB)r(A),r(B)

4. 将AAA看作一块,BBB按列分块

A(η1,η2,⋯ ,ηt)=(Aη1,Aη2,⋯ ,Aηt)A(\eta_1,\eta_2,\cdots,\eta_t)=(A\eta_1,A\eta_2,\cdots,A\eta_t)A(η1,η2,,ηt)=(Aη1,Aη2,,Aηt)
AB=0AB=0AB=0,则Aηj=0(j=1,2,⋯ ,t)A\eta_j=0(j=1,2,\cdots,t)Aηj=0(j=1,2,,t),ηj\eta_jηjAX=0AX=0AX=0的解向量,因此由n−rn-rnr个向量是线性无关的,rrrAAA的秩
r(η1,η2,⋯ ,ηt)≤n−rr(B)+r(A)≤nr(\eta_1,\eta_2,\cdots,\eta_t)\le n-r\\r(B)+r(A)\le nr(η1,η2,,ηt)nrr(B)+r(A)n

5. 将矩阵分成四块

例4.
证明:上三角矩阵的乘积仍是上三角矩阵
证:
A=[A1α0an×n],B=[B1β0bn×n]A=\left[\begin{array}{c|c}A_1&\alpha\\ \hline 0 & a_{n\times n}\end{array}\right] ,B=\left[\begin{array}{c|c}B_1 & \beta\\ \hline 0 & b_{n \times n} \end{array}\right]A=[A10αan×n],B=[B10βbn×n]
AB=[A1B1A1β+αbn×n0an×nbn×n]AB=\left[\begin{matrix} A_1B_1 & A_1\beta+\alpha b_{n\times n} \\ 0 & a_{n\times n}b_{n\times n}\end{matrix}\right]AB=[A1B10A1β+αbn×nan×nbn×n]
由上式可知,上三角矩阵的乘积仍是上三角矩阵。

三、 非齐次线性方程组

1. 线性方程组

对于非齐次线性方程组。
AX=bAX=bAX=b,其中A=(aij)s×n,b=(b1,b2,⋯ ,bs)TA=(a_{ij})_{s\times n},b=(b_1,b_2,\cdots,b_s)^TA=(aij)s×n,b=(b1,b2,,bs)T

  1. 有解⇔r(A)=r(A,b)\Leftrightarrow r(A)=r(A,b)r(A)=r(A,b)
  2. r(A)=r(A,b)=rr(A)=r(A,b)=rr(A)=r(A,b)=r则有唯一解⇔r=n\Leftrightarrow r=nr=n
  3. r(A)=r(A,b)=r<nr(A)=r(A,b)=r<nr(A)=r(A,b)=r<n,则通解中含有r−nr-nrn个自由未知量。

对于齐次线性方程组
AX=0AX=0AX=0,其中A=(aij)s×nA=(a_{ij})_{s\times n}A=(aij)s×n

  1. 有非零解当且仅当r(A)<n,n1,n2,⋯ ,ntr(A)<n,n_1,n_2,\cdots,n_tr(A)<nn1,n2,,nt(线性无关并且任一解都可由它们线性表示)是AX=0AX=0AX=0的解。
  2. r(A)<nr(A)<nr(A)<n,则其基础解系中含有n−rn-rnr个解向量。
  3. r(A)<rr(A)<rr(A)<r,则其任意n−rn-rnr个线性无关的解向量是其基础解系。

2. Gauss消元法

用初等变换将增广矩阵化成阶梯形矩阵。第一步,确定自由变量。第二步,用回代法找出通解。
例4.
AAAs×ns\times ns×n矩阵,bbbsss维列向量,证明

  1. r(A)=r(AHA)r(A)=r(A^HA)r(A)=r(AHA)
  2. 线性方程AHAX=AHbA^HAX=A^HbAHAX=AHb恒有解。

注:AHA^HAHAAA的共轭转置。

  1. 证:
    要证r(A)=r(AHA)r(A)=r(A^HA)r(A)=r(AHA),即证AX=0AX=0AX=0的基础解系S1S_1S1AHAX=0A^HAX=0AHAX=0的基础解系S2S_2S2相同。即S1∈S2S_1 \in S_2S1S2的同时S2∈S2S_2 \in S_2S2S2
    显然
    Aη=0⇒AHAηA\eta=0 \Rightarrow A^HA\etaAη=0AHAη
    AX=0AX=0AX=0的解必然是AHAX=0A^HAX=0AHAX=0的解。
    AHAη=0A^HA\eta=0AHAη=0
    等号两边同时乘以ηH\eta^HηH,得
    ηHAHAη=ηH0\eta^HA^HA\eta=\eta^H 0ηHAHAη=ηH0
    ηˉTAˉTAη=0\bar{\eta}^T\bar{A}^TA\eta=0ηˉTAˉTAη=0
    (Aˉηˉ)TAη=0({\bar{A}\bar{\eta}})^TA\eta=0(Aˉηˉ)TAη=0
    (Aη)HAη=0(A\eta)^HA\eta=0(Aη)HAη=0
    因为
    Aη=[k1k2⋮kn]A\eta=\left[\begin{matrix}k_1 \\ k_2 \\ \vdots \\ k_n \end{matrix}\right]Aη=k1k2kn
    所以
    (Aη)HAη=(kˉ1,kˉ2,⋯ ,kˉs)[k1k2⋮kn]=kˉ1k1+kˉ2k2+⋯+kˉsks=0(A\eta)^HA\eta=(\bar{k}_1,\bar{k}_2,\cdots,\bar{k}_s)\left[\begin{matrix}k_1 \\ k_2 \\ \vdots \\ k_n \end{matrix}\right]=\bar{k}_1k_1+\bar{k}_2k_2+\cdots+\bar{k}_sk_s=0(Aη)HAη=(kˉ1,kˉ2,,kˉs)k1k2kn=kˉ1k1+kˉ2k2++kˉsks=0
    因为
    kˉ1k1≥0,kˉ2k2≥0,⋯ ,kˉsks≥0\bar{k}_1k_1\ge 0,\bar{k}_2k_2\ge 0,\cdots,\bar{k}_sk_s\ge 0kˉ1k10,kˉ2k20,,kˉsks0
    因此
    k1=k2=⋯=ksk_1=k_2=\cdots=k_sk1=k2==ks
    ⇓\Downarrow
    Aη=0A\eta=0Aη=0
    2)证:
    要证AHAX=AHbA^HAX=A^HbAHAX=AHb恒有解,即证r(AHA)=r(AHA,AHb)r(A^HA)=r(A^HA,A^Hb)r(AHA)=r(AHA,AHb)
    显然
    r(AHA)≤r(AHA,AHb)r(A^HA)\le r(A^HA,A^Hb)r(AHA)r(AHA,AHb)
    r(AHA,AHb)=r(AH(A,b))≤r(AH)=r(A)r(A^HA,A^Hb)=r(A^H(A,b))\le r(A^H)=r(A)r(AHA,AHb)=r(AH(A,b))r(AH)=r(A)
    由第一问可知
    r(AHA)=r(A)r(A^HA)=r(A)r(AHA)=r(A)
    所以
    r(AHA,AHb)≤r(AHA)r(A^HA,A^Hb)\le r(A^HA)r(AHA,AHb)r(AHA)
    因此
    r(AHA)=r(AHA,AHb)r(A^HA)=r(A^HA,A^Hb)r(AHA)=r(AHA,AHb)
    AHAX=AHbA^HAX=A^HbAHAX=AHb恒有解

四、向量组的极大无关组和秩

若向量组α1,α2,⋯ ,αs\alpha_1,\alpha_2,\cdots,\alpha_sα1,α2,,αs的部分组αi1,αi2,⋯ ,αik\alpha_{i1},\alpha_{i2},\cdots,\alpha_{ik}αi1,αi2,,αik线性无关。且α1,α2,⋯ ,αs\alpha_1,\alpha_2,\cdots,\alpha_sα1,α2,,αs中每个向量都可由αi1,αi2,⋯ ,αir\alpha_{i1},\alpha_{i2},\cdots,\alpha_{ir}αi1,αi2,,αir线性表示,则称αi1,αi2,⋯ ,αir\alpha_{i1},\alpha_{i2},\cdots,\alpha_{ir}αi1,αi2,,αir是向量组α1,α2,⋯ ,αs\alpha_1,\alpha_2,\cdots,\alpha_sα1,α2,,αs的极大无关组。称rrr是向量组α1,α2,⋯ ,αs\alpha_1,\alpha_2,\cdots,\alpha_sα1,α2,,αs的秩。
若向量组α1,α2,⋯ ,αs\alpha_1,\alpha_2,\cdots,\alpha_sα1,α2,,αs的秩为rrr,则其中任意rrr个线性无关的向量是其极大无关组。

1. 矩阵的秩

矩阵A的秩=A中非零子式的最高阶数=A的行(列)向量组的秩

2. 有关矩阵的秩的不等式

  1. r(A+B)≤r(A)+r(B)r(A+B)\le r(A)+r(B)r(A+B)r(A)+r(B)
  2. r(AB)≤r(A),r(B)r(AB)\le r(A),r(B)r(AB)r(A),r(B)
  3. As×nBn×t=0r(A)+r(B)≤nA_{s\times n}B_{n\times t}=0 r(A)+r(B)\le nAs×nBn×t=0r(A)+r(B)n
  4. r(As×nBn×t)≥r(A)+r(B)−nr(A_{s\times n} B_{n\times t})\ge r(A)+r(B)-nr(As×nBn×t)r(A)+r(B)n
  5. M=[AC0B],r(M)≥r(A)+r(B)M=\left[\begin{matrix}A&C\\0&B\end{matrix}\right],r(M)\ge r(A)+r(B)M=[A0CB],r(M)r(A)+r(B)

例5.
AAA是可逆矩阵,证明r(AB)=r(B)r(AB)=r(B)r(AB)=r(B)
法一:
显然
r(AB)≤r(B)r(AB)\le r(B)r(AB)r(B)
同时
B=A−1ABB=A^{-1}ABB=A1AB
可得
r(B)≤r(A−1),r(AB)r(B)\le r(A^{-1}),r(AB)r(B)r(A1),r(AB)

r(B)≤r(AB)r(B)\le r(AB)r(B)r(AB)
所以
r(AB)=r(B)r(AB)=r(B)r(AB)=r(B)
法二:
∃\exists初等矩阵A1,A2,⋯ ,AsA_1,A_2,\cdots,A_sA1,A2,,As,使A=A1A2⋯AsA=A_1A_2{\cdots}A_sA=A1A2As
AB=(A1,A2,⋯ ,As)BAB=(A_1,A_2,\cdots,A_s)BAB=(A1,A2,,As)B
BBB经过一系列初等行变换变成A⇒r(AB)=r(B)A \Rightarrow r(AB)=r(B)Ar(AB)=r(B)

例6.
AAAnnn阶幂等矩阵,证明r(A)+r(I−A)=nr(A)+r(I-A)=nr(A)+r(IA)=n
证:
r(A)+r(I−A)≥r(A+I−A)=r(I)=nr(A)+r(I-A)\ge r(A+I-A)=r(I)=nr(A)+r(IA)r(A+IA)=r(I)=n
同时,由于AAA为幂等矩阵
A(I−A)=A−A2=0A(I-A)=A-A^2=0A(IA)=AA2=0
r(A)+r(I−A)≤nr(A)+r(I-A)\le nr(A)+r(IA)n
所以
r(A)+r(I−A)=nr(A)+r(I-A)=nr(A)+r(IA)=n

3. 矩阵的等价标准型

如果s×ns\times ns×n矩阵的秩为rrr,则AAA必与M=[Ir000]M=\left[\begin{matrix}I_r&0\\0&0\end{matrix}\right]M=[Ir000]等价
s×ns\times ns×n矩阵A,BA,BA,B等价⇔r(A)=r(B)⇔\Leftrightarrow r(A)=r(B) \Leftrightarrowr(A)=r(B)存在可逆矩阵P,QP,QP,Q使B=PAQB=PAQB=PAQ

4. 满秩分解

例7.
假设s×ns\times ns×n矩阵AAA的秩为rrr,证明存在s×rs\times rs×r矩阵BBBr×nr\times nr×n矩阵CCC,使得A=BCA=BCA=BC(矩阵的满秩分解)
证:
r(A)≤r(B),r(C)≤r⇒r(B)=r(C)=rr(A)\le r(B),r(C)\le r \Rightarrow r(B)=r(C)=rr(A)r(B),r(C)rr(B)=r(C)=r
由矩阵的等价标准型可知,∃s×n\exists s\times ns×n矩阵M=[Ir000]M=\left[\begin{matrix}I_r&0\\0&0\end{matrix}\right]M=[Ir000]使得AAAMMM等价,s×ns\times ns×n矩阵[Ir000]\left[\begin{matrix}I_r&0\\0&0\end{matrix}\right][Ir000]等于s×rs\times rs×r矩阵[Ir0]\left[\begin{matrix}I_r\\0\end{matrix}\right][Ir0]乘以r×nr\times nr×n矩阵(Ir0)(I_r 0)(Ir0)
A=P[Ir000]Q=P[Ir0](Ir0)QB=P[Ir0]C=(Ir0)QA=P\left[\begin{matrix}I_r&0\\0&0\end{matrix}\right]Q=P\left[\begin{matrix}I_r\\0\end{matrix}\right](I_r 0)Q \\B=P\left[\begin{matrix}I_r\\0\end{matrix}\right] C=(I_r 0)QA=P[Ir000]Q=P[Ir0](Ir0)QB=P[Ir0]C=(Ir0)Q
由此可知存在s×rs\times rs×r矩阵BBBr×nr\times nr×n矩阵CCC,使得A=BCA=BCA=BC
例8.
A=[10−1212310−1132−2−20−3−34−3]A=\left[\begin{matrix}1&0&-1&2&1\\2&3&1&0&-1\\1&3&2&-2&-2\\0&-3&-3&4&-3\end{matrix}\right]A=12100333112320241123的满秩分解
A→[10−1212310−1132−2−20−3−34−3]→[10−1010110−10001000000]A\rightarrow\left[\begin{matrix}1&0&-1&2&1\\2&3&1&0&-1\\1&3&2&-2&-2\\0&-3&-3&4&-3\end{matrix}\right]\rightarrow\left[\begin{matrix}1&0&-1&0&1\\0&1&1&0&-1\\0&0&0&1&0\\0&0&0&0&0\end{matrix}\right]A1210033311232024112310000100110000101100
A=[10223013−20−34][10−1010110−100010]A=\left[\begin{matrix}1&0&2\\2&3&0\\1&3&-2\\0&-3&4\end{matrix}\right]\left[\begin{matrix}1&0&-1&0&1\\0&1&1&0&-1\\0&0&0&1&0\end{matrix}\right]A=121003332024100010110001110

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值