· Matrix
Fragment 1 矩阵概念
/Define/
将 mn 个数 aij (1≤i≤m, 1≤j≤n) 排成一个矩阵,称为一个 m×n 阶矩阵:
A=a11a21⋮am1a12a22⋮am2⋯⋯⋱⋯a1na2n⋮amn=(aij)m×n
aij 称为第 (i,j) 元素。若 aij∈R(或 C),则称 A 为实(复)矩阵;
若 ∀i,j, aij=0,则称 A 为零矩阵,记作 Om×n
若 m=n,则称 A 为 n 阶方阵.
设 A=(aij)n×n 为方阵,则 a11,a22,…,ann 连成的线称为 A 的主对角线.
(一般对称方阵才称主对角线)
若 ∀i=j, aij=0,则称 A 为对角矩阵.
A=a110⋮00a22⋮0⋯⋯⋱⋯00⋮ann
简记为diag {a11,⋯,ann}
进一步,若a11=a22=⋯=ann=1 ,这样的矩阵称为单位阵.记作:
In=10⋮001⋮0⋯⋯⋱⋯00⋮1
单位阵在矩阵乘法里起的作用与 1 在数字乘法里起的作用一样.
在有了主对角线之后,我们可以定义上下三角矩阵
/Define/
若 aij=0 (∀i>j),则称 A 为上三角阵;
若 aij=0 (∀i<j),则称 A 为下三角阵.
那我们如何判断两个矩阵是否相等?
/Define/
设 A=(aij)m×n, B=(bij)s×t
A=B⟺defm=s,n=t,aij=bij(∀1≤i≤m,∀1≤j≤n)
两矩阵相等,当且仅当它们一模一样 (阶数相等,元素相同)
/example/
[0000]2×2=[000000]2×3,[1000]=[1000]
接下来接着给定义。 我们将空间解析几何中,三维的行/列向量,利用矩阵的定义推广到一般的n维的行/列向量。
/Define/
设 A=(aij)m×n,则 [ai1,ai2,…,ain] 称为 A 的第 i 个行向量, a1ja2j⋮amj 称为 A 的第 j 个列向量
称 1×n 矩阵 (a1,a2,…,an) 为 n 维行向量,称 n×1 矩阵 a1a2⋮an 为 n 维列向量
接下来我们看矩阵和行列式有什么联系:
/Claim/
记 Mn(R/C) 为所有 n 阶实方阵 (复方阵) 全体构成的集合
映射(函数):
det:Mn(R)⟶RA⟼∣A∣
每个 n 阶方阵对应过去都有一个数值(行列式),这样构造成一个映射(函数).
这里就诞生了两个问题:
/question/
(1) n 阶行列式的值在多大程度上反映 n 阶方阵的性质? (2) 映射 det 具有怎样的性质?
我们需要先定义行列式运算,才能给出性质与意义.
Fragment 2 矩阵运算
下面会给出矩阵的加减法、数乘、乘法、乘方、转置、共轭.
· 矩阵加减
/Define/
设 A=(aij)m×n, B=(bij)m×n
A+B=def(aij+bij)m×n
A−B=def(aij−bij)m×n
Am×n±Om×n=Am×n
/example/
[1−12301]+[1−1−2−310]=[2−20011]
/property/
(1) 交换律: A+B=B+A
(2) 结合律: (A+B)+C=A+(B+C)
(3) 零矩阵: A+0=A (加上同阶零矩阵)
(4) 负矩阵: A+(−A)=0 (5) A−B=A+(−B)
· 矩阵数乘
/Define/
设 A=(aij)m×n, C 为常数
C 与 A 的数乘为(scalar product) C⋅A=def(C⋅aij)m×n
负矩阵 −A=def(−1)⋅A=(−aij)m×n
/example/
2⋅[1324]=[2648]
/property/
(1) 数的分配律: (c+d)⋅A=c⋅A+d⋅A
(2) 矩阵的分配律: C⋅(A+B)=C⋅A+C⋅B
(3) 数乘结合律: (c⋅d)⋅A=C⋅(d⋅A)
(4) 数乘单位元: 1⋅A=A
(5) 数乘的零元: 0⋅Am×n=Om×n
· 矩阵乘法
/Define/
设 A=(aij)m×k, B=(bij)k×n, 定义不同于 A 与 B 的点积 C 是 m×n 矩阵, 其元素:
Cij=ai1b1j+ai2b2j+⋯+aikbkj=r=1∑kairbrj
不同于矩阵数乘。
”定义一看, 来得毫无征兆“
为什么不像前面一样对应元素相加?
——历史上有过定义, 但用处不大。
(1).
只有 A 的列数 = B 的行数, AB 才有意义
AB 的行数 = A 的行数, AB 列数 = B 的列数
Am×kBk×n⟶(AB)m×n
(2).
AB 的第 (i,j) 元素 = A 的第 i 行与 B 的第 j 列 对应元素乘积之和
(3).
矩阵乘法一般不满足交换律 → 一般来说, 即使 AB, BA 都有意义, 但 AB=BA
/example/
[120110]2×311−101012−11−103×4=[03010411]
/example/
A=[111]B=111AB=[3]=AB
试反驳: A, B 不同阶 用同阶矩阵举例
/example/
A=[0010]B=[0001]AB=[0010]=BA=[0000]
结论:矩阵乘法一般来说, 不满足交换律
/property/
(1) 结合律: (AB)C=A(BC)
(2) 分配律: (A+B)C=AC+BC ,A⋅(B+C)=AB+AC
(3) 与数乘的相容性 c⋅(AB)=(c⋅A)B=A(c⋅B)
(4) 单位元: 设 Am×n, 则 ImA=A=AIn
然后我们来康康性质(1)的证明(很有精神地打字)
/proof/
A=(aij)m×n, B=(bij)n×p, C=(cij)p×q
对于 (AB)C,先考虑 AB 的第 (i,j) 元素 ∑k=1naikbkj
再考虑 (AB)C 的第 (i,j) 元素
r=1∑p(k=1∑naikbkr)crj=r=1∑pk=1∑naikbkrcrj
对于 A(BC) ,先考虑 BC 的第 (k,i) 元素 ∑r=1pbkrcri ⇒ 第 (i,j) 元素
k=1∑naik(r=1∑pbkrcri)=k=1∑nr=1∑paikbkrcri
⇒(AB)⋅C=A(BC) 成立
PS:以后写 3 个矩阵相乘时:不用打括号
· 矩阵乘方
/Define/
乘方是一种特殊的乘法
∀A1,A2,…,Am ,乘积写为 A1A2⋯Am
设 A 为 n 阶方阵,则 An 有定义 Ak=defA⋅A⋯A (k 个)
PS:如果一个矩阵有平方,则它一定是方阵
“那方阵有什么性质?”
/property/
(1) Ar⋅As=Ar+s (r,s∈Z+)
(2) (Ar)s=Ars 先前乘法的性质都满足。
关键点:
(1) 设 A,B 为 n 阶方阵,一般地,(AB)r=ArBr (无交换律)
关于(1),这里有一个例子:
/example/
A=[1324],B=[0010]
PS:后面的很多例子都和 [0010] 有关
经过计算,我们不难发现:AB=BA ,(AB)2=0,A2B2=0
(2) 若 AB=BA,则 (AB)r=ArBr,进一步有
(A+B)m=Am+Cm1Am−1B+⋯+Cmm−1ABm−1+Bm
(若 AB=BA ,中学里所有关于数的公式都能用)
特别地: A⋅In=In⋅A=A
A(CIn)=(CIn)A=CA
(cIn+A)m=CmIn+Cm1cm−1A+⋯+Cmm−1cAm−1+Am
此时二项式展开不需要任何条件.
(3) 矩阵乘法的律一般不成立
/example/
A=B=[0010],C=[0000]A=B=0,AB=0
· 矩阵转置
/Define/
设 A=(aij)m×n, A 的转置,其中 AT 是 n×m 阶矩阵, AT=(bij)n×m
其中 bij=aji, 即 AT 的第 i 行是 A 的第 i 列, AT 的第 j 列是 A 的第 j 行.
把原来 A 的行变成列, 即转置.
接下来给出对称阵与反对称阵的定义。
/Define/
若 AT=A,则称 A 为对称阵。A=[1324]
若 AT=−A,则称 A 为反对称阵。A=[0−110]
接下来给出矩阵转置的性质:
/property/
(1) ATT=A
(2) (A+B)T=AT+BT
(3) (c⋅A)T=c⋅AT,c 是一个数
(4) (A⋅B)T=BT⋅AT
证明过程在此略去.
· 矩阵共轭
/Define/
设 A=(aij)m×n 复矩阵,复转置 A=(aij)m×n a+bi→a−bi
对复矩阵每个元素取共轭,即矩阵共轭.
/example/
A=1+i2i−4−i0−13i1+2i−4ii,A=1−i−2i−4+i0−1−3i1−2i4i−i
然后是共轭的性质:
/property/
(1) A=A (2) A+B=A+B (3) c⋅A=c⋅A (4) AB=A⋅B (5) A′=A′
到这里矩阵计算的内容基本结束了
Fragment 3 方阵的逆阵
上节引入矩阵定义、性质,最重要:乘法;定理的含义。
在第四章,我们会了解矩阵的几何意义是线性映射,矩阵乘法是线性映射的复合
矩阵乘法的应用可以是解决线性方程组:
(∗)⎩⎨⎧a11x1+a12x2+⋯+a1nxn=b1a21x1+a22x2+⋯+a2nxn=b2⋯am1x1+am2x2+⋯+amnxn=bm
每次写线性方程组会很烦,但是我们引入矩阵乘法后会变得简洁
我们让所有系数构成一个矩阵(线性方程组的系数矩阵):
A=a11a21⋮am1a12a22⋮am2⋯⋯⋱⋯a1na2n⋮amn
我们把所有的未定元拼成列向量,把常数项拼成列向量:
x=x1x2⋮xnb=b1b2⋮bm
自然地,根据矩阵乘法 Ax=b⟺(∗)
写成矩阵能否借此来解线性方程组?(第三章)
· 定义
我们学过 b=0, ba=a⋅b−1
那矩阵能否求逆: A∼A−1
自然地想:Ax=b⇒x=A−1b
那我们给出矩阵逆阵的定义:
/Define/
设 A 为 n 阶方阵,若存在 n 阶方阵 B,使得
AB=BA=In
(n 阶单位阵),则称 B 为 A 的逆阵,记为 A−1
若矩阵 A 存在逆阵,则称为非奇异阵或可逆阵 (nonsingular , invertible)若矩阵 A 不存在逆阵,则称为奇异阵 (singular).
由式 AB=BA=In,In 迫使 AB 定为 n 阶方阵。
/remark/
(1) 只有对方阵才有逆阵的定义,当 m=n 时,没有逆阵的定义
(2) 非零方阵不一定有逆阵:
/example/
A=[1010]B=[b11b21b12b22]
AB=[b11+b210b12+b220]=[1001]
⇒A 不可逆
当然,我们也可以自然引出问题:一个方阵有逆的充分必要条件是什么? (本节课主要定理)
(3) 一般而言 A−1B=BA−1 矩阵乘法一般来说不可交换。
· 性质
/property/
(i) 若 A 可逆,则逆阵必唯一。
( A 有无逆阵不确定,但若存在必唯一)
/proof/
设 B, C 均为 A 的逆阵,即 AB=BA=In,AC=CA=In
B=BIn=B(AC)=(BA)C=InC=C
此处主要是矩阵乘法的结合律。
(iii) 设 A, B 可逆,则 AB 也可逆,(AB)−1=B−1A−1
/proof/
(AB)(B−1A−1)=A(BB−1)A−1=AInA−1=AA−1=In
(B−1A−1)(AB)=B−1(A−1A)B=B−1InB=B−1B=In
推广:
设 A1, ⋯, Am 均为 n 阶可逆阵 (m⩾2)
则 A1, ⋯, Am 可逆,(A1⋯Am)−1=Am−1⋯A1−1
是可逆阵!
(A1⋯Am)(Am−1⋯A1−1)=In,(Am−1⋯A1−1)(A1⋯Am)=In
(有限个同阶可逆矩阵乘积仍是可逆阵)
(iii) 设 A 可逆,c=0 常数,则 cA 仍可逆(c⋅A)−1=c−1A−1
/proof/
(cA)(c−1A−1)=(cc−1)AA−1
(c−1A−1)(cA)=(c−1c)A−1A
(iiii) 设 A 可逆,则 A′ 也可逆且 (A′)−1=(A−1)′
/proof/
A′(A−1)′=(A−1A)′=In′=In
(A−1)′A′=(A′A−1)′=In′=In
(v) 设 A 可逆,则 (A−1)−1=A
/proof/
A−1A=A⋅A−1=In=(A+I)−1=A
(vi) 对可逆阵而言,乘法消去律成立
设 A 为可逆:
{AB=AC⟹B=CBA=CA⟹B=C
/proof/
A−1AB=A−1AC⟹InB=InC⟹B=C
两边同时左乘 A−1 即 B=C
BAA−1=CAA−1⟹BIn=CIn⟹B=C
/remark/
先前说过,矩阵乘法还要有个特性
两个不等于0的矩阵相等,结果不一定是零矩阵
(vii) 整性对可逆阵成立,即A 可逆,
B=0⟹AB=0
C=0⟹CA=0
/proof/
设 AB=0=A⋅0⟹B=0 , 矛盾
利用 若 AB=0 且 A 可逆,则 B=0
· 伴随矩阵
介绍完可逆矩阵性质,我们自然提出一个问题,怎么判断和求解矩阵逆阵?
于是就引入伴随矩阵:
/Define/
设 A=(aij)n×n,Aij 是 aij 在 ∣A∣ 中的代数余子式,则下列矩阵称为 A 的伴随矩阵,记为 A∗
A∗=A11A12⋮A1nA21A22⋮A2n⋯⋯⋱⋯An1An2⋮Ann
注意:A 不一定有逆阵,但伴随阵一定存在。
通过代数余子式定义,( n−1 )阶行列式都可以计算,不依赖是否可逆
/example/
A=[1010],A∗=[00−11]
PS:此处计算涉及转置与 (−1)i+j ,易错
引理:设 A 为 n 阶方阵,则 AA∗=A∗A=∣A∣In。
/proof/
a11a21⋮an1a12a22⋮an2⋯⋯⋱⋯a1na2n⋮annA11A12⋮A1nA21A22⋮A2n⋯⋯⋱⋯An1An2⋮Ann=∣A∣0⋮00∣A∣⋮0⋯⋯⋱⋯00⋮∣A∣=∣A∣In.
上式矩阵中第 (i,j) 个元素如下:
ai1Aj1+ai2Aj2+⋯+ainAjn=δij⋅∣A∣
根据第一章内容,我们容易证明上式
定理:设 A 为 n 阶方阵,若 ∣A∣=0,则 A 可逆,且 A−1=∣A∣1A∗。
/proof/
AA∗=A∗A=∣A∣In。又 A(∣A∣1A∗)=(∣A∣1A∗)A=In。
定理: 设 A,B 为 n 阶方阵,则 ∣AB∣=∣A∣∣B∣。
/proof/
设 A=(aij)n×n,B=(bij)n×n,构造 2n 阶方阵 C=[AOOB]。
∣C∣=a11⋮an1−1⋮0⋯⋱⋯⋯⋱⋯a1n⋮ann0⋮−10⋮0b11⋮bn1⋯⋱⋯⋯⋱⋯0⋮0b1n⋮bnn=0⋮00⋮0⋯⋱⋯⋯⋱⋯a1n⋮ann0⋮0a11b11⋮an1bn1b11⋮bn1⋯⋱⋯⋯⋱⋯a1nbnn⋮an1bnnb1n⋮bnn.
用 aij 乘以 C 的第 i+j 行加到第 j 行上,得
∣C∣=0⋮0−1⋮0⋯⋱⋯⋯⋱⋯0⋮00⋮−1Σa1rbr1⋮Σanrbr1b11⋮bn1⋯⋱⋯⋯⋱⋯Σa1rbrn⋮Σanrbrnb1n⋮bnn=O−InABB.
用 Laplace 定理展开,按前 n 行展开
左边 =∣A∣∣B∣,右边 =∣AB∣。
∣AB∣=∣A∣∣B∣.
定理6:设 A 为 n 阶方阵,则 A 可逆 ⇔∣A∣=0。
/proof/
充分性定理已经证明,只需证明必要性:
A 可逆,由定义 ∃ n 阶方阵 B 使 AB=BA=In。
1=∣In∣=∣BA∣=∣B∣∣A∣, 从而 ∣A∣=0.
推论:
- 设 A1,…,Am 为 n 阶方阵,若 ∃i 使 i 是奇数,则 A1⋯Am 是奇异矩阵。
证明:∣A1⋯Am∣=∣A1∣⋯∣Am∣=0。由以上知 A1⋯Am 奇异。
-
设 A 可逆,则 ∣A−1∣=∣A∣−1。
-
设 A,B 为 n 阶方阵,若 AB=In 并且 BA=In,则 B=A−1。
/proof/
只证 AB=In 情形
1=∣In∣=∣AB∣=∣A∣∣B∣⇒∣A∣=0⇒A−1 exist
B=InB=(A−1A)B=A−1(AB)=A−1A=In
/example/
A=(a1+b1)n(a2+b1)n⋮(an+b1)n(a1+b2)n(a2+b2)n⋮(an+b2)n⋯⋯⋱⋯(a1+bn)n(a2+bn)n⋮(an+bn)n(ai+bj)n=k=0∑nCknaikbjn−k
/Claim/
有时间再补充
/example/
A=xyzwy−xwzz−wx−ywz−y−x
/Claim/
AA′=u0000u0000u0000uu=x2+y2+z2+w2
∣A∣2=∣AA′∣=u4∣A∣=±u=±(x2+y2+z2+w2)2
令 x=1, y=z=w=0,
∣A∣=1⇒∣A∣=(x2+y2+z2+w2)2
Fragment 4 矩阵初等变换
· 矩阵的初等变换
行列式性质:
设 A,B 为 n 阶方阵,c 为常数。
-
∣A+B∣=∣A∣+∣B∣
-
∣cA∣=cn∣A∣
-
∣AB∣=∣A∣∣B∣
-
∣AAT∣=∣A∣2
-
若 A 可逆,则 ∣A−1∣=∣A∣1
-
设 n≥2,则 ∣A∗∣=∣A∣n−1
/proof/
先证 A 为可逆阵,A−1=ATA∗,A∗=∣A∣A−1
∣A∗∣=∣A∣∣A−1∣=∣A∣∣ATA∗∣=∣A∣∣AT∣∣A∗∣=∣A∣∣A∣−1∣A∗∣=∣A∣n−1
若 A 为奇异阵,即 ∣A∣=0,证 ∣A∗∣=0 即可。
反证:设 A∗=0,则 A∗ 为非奇异阵。
AA∗=A∗A=∣A∣In=0⇒A>0⇒A∗=0 矛盾。
/example/ 用Gauss消元法求解
PS: 系数矩阵+常数列为增广矩阵。
给出如下定义:
定义1 矩阵的初等变换
第一类:对换矩阵的两行(列)
第二类:矩阵某一行(列)乘上非零常数
第三类:矩阵某一行乘常数加到另一行上
定义2
若矩阵 A 通过若干的初等变换变为矩阵 B,则称 A 相抵于 B,记 A∼B
定义3
设 A=(aij)m×n 相抵于 (Ir000) 称为 A 的相抵标准型。
注意:r≤min{m,n},r 不依赖于初等变换选取,r 唯一确定。
/proof/
对 min{m,n} 进行归纳。若 min{m,n}=0,则归纳过程结束。
设 min{m,n}<k 归纳成立,下证 min{m,n}=k 情形。
若 A=0 则结论成立,下设 A=0,不妨设 aij=0,i≤r, j≤s.
此时 aij 变到了 (1,1) 位置,以下不妨设 aij=1.
a11⋮am1⋯⋯a1n⋮amn→10⋮00b22⋮bm2⋯⋯⋱⋯0b2n⋮bmn
此时 B 是 (m−1)×(n−1) 矩阵,
min{(m−1),(n−1)}=min{m,n}−1=k−1.
由归纳假设,B 可变换为 [Ir000],A∼[Ir000].
/example/ (1,2,3)∼(1,0,0)
定义4
设 A=(aij)m×n,对任给 m,
ki={+∞minj≥1{j∣aij=0}if ai1=0if ai1=0
给出阶梯点定义:
A 为阶梯形矩阵 ⟺def 存在 r,使得 k1<k2<⋯<kr,kr+1=⋯=km=+∞。
定义5:设 A=(aij)m×n,则通过初等变换可将 A 化为阶梯形矩阵。
/proof/
对列数 n 进行归纳,n=0 表示归纳过程已结束。
设列数 <n 时结论成立,现证 n 列情况。
Case 1
A 的第一列元素为 0,A=0⋮A1,A1 为 m×(n−1) 矩阵。
由归纳假设,A1⟶初等变换B1(阶梯形)。
则 A⟶初等变换0⋮0⋯B10 仍为阶梯形。
Case 2
A 的第一列不全为 0,由行变换可将非零元移到第一行。
不妨设 a11=0,
A=a11a21⋮an1a12⋯∗a1n→a110⋮0a12⋯A2a1n
A2 是 (m−1)×(n−1) 矩阵,由归纳假设,
A2⟶初等变换B2(阶梯形)。
则 A⟶初等变换a11⋮0a12⋯B2a1n
· 初等矩阵
定义6
对单位矩阵 In 实施三类初等行变换后得到的矩阵称为三类初等矩阵。
Pij=1⋱0110⋱1
Pi(c)=1⋮0⋮0⋯⋯⋯0⋮c⋮0⋯⋯⋯0⋮0⋮1
Tij(c)=1⋱1c1⋱1
定理7
初等矩阵行列变换等价于左(右)乘对应初等矩阵,左行右列
CjiA=PijA⇒A=Pi(c)ACjiiA=Tij(c)A
引理8
- 初等矩阵都是可逆矩阵,且逆矩阵为同类初等矩阵。
Pij−1=PijPi(c)−1=Pi(−c)Tij(c)−1=Tij(−c)
- 矩阵实施第三类初等行变换后 ∣A∣ 不变。
∣Pij∣=−1∣Pi(c)∣=c∣Tij(c)∣=1
- 非奇异矩阵经初等变换后仍为非奇异矩阵。奇异矩阵 ∼ 奇异矩阵
/proof/
定理7 ⇒ A实施初等变换等价于 P1⋯PrAQ1⋯Qs
其中 Pi,Qj 都是非奇异矩阵,从而都是非奇异矩阵。
- A非异 ⇒ B非异
- A奇异 ⇒ B奇异
relation
集合A :给定子集 R⊂A×A={(a,b)∣a,b∈A},若 (a,b)∈R 则称 a R b
若又满足以下三条性质,则称为等价关系:
- 自反性:a R a
- 对称性:a R b 且 b R a
- 传递性:若 a R b 且 b R c,则 a R c
定理9
矩阵的相抵关系是等价关系即:
- 自反性 A∼A
- 对称性 A∼B⇔B∼A
- 传递性 A∼B 且 B∼C⇒A∼C
/proof/
-
A=In⋅A⇒A∼A
-
由 A∼B 可知,存在初等矩阵 P1⋯Pr 和 Q1⋯Qs 使得
B=P1⋯PrAQ1⋯Qs⇒P1⋯PrBQ1⋯Qs−1=A⇒B∼A
-
B=P1⋯PrAQ1⋯Qs
C=R1⋯RtBT1⋯Tl
其中 Pi,Qj,Rk,Tl 都是初等矩阵。
⇒C=R1⋯RtP1⋯PrAQ1⋯QsT1⋯Tl⇒A∼C
下面研究非异阵的一些性质:
定理1
设 A 为 n 阶方阵,则下列结论等价:
- A 为非奇异矩阵
- A 的相抵标准型为 In
- A 通过初等行变换或列变换能变为 In
- A 是若干个初等矩阵的乘积
/proof/
(1) ⇒ (2) 设 A 相抵标准型为 (Ir000),由 A 的非奇异性和初等矩阵性质可知 r=n。
(2) ⇒ (3) 设 A∼In 即 ∃P1⋯Pr,Q1⋯Qs 使得 P1⋯PrAQ1⋯Qs=In
P1⋯PrA=Ini.e.Q1⋯QsP1⋯PrA=In
(3) ⇒ (4)
设 A 通过行变换变为 In,即存在非奇异矩阵 P1,P2,…,Pr 使得 Pr⋯P1A=In。
⇒A=(Pr⋯P1)−1=P1−1⋯Pr−1
即 A 是若干个初等矩阵的乘积。
(4) ⇒ (1)
设 A=P1⋯Pr,Pi 为初等矩阵。则 Pi 为可逆矩阵 ⇒A 为可逆矩阵。
推论2
设 A=(aij)m×n,则存在 m 阶非奇异矩阵 P 和 n 阶非奇异矩阵 Q 使得 PAQ=(Ir000)。
/proof/
存在 m 阶非奇异矩阵 P 和 n 阶非奇异矩阵 Q 使得 P⋯P1AQ1⋯Qs=(Ir000)。
令 P=Pr⋯P1,Q=Q1⋯Qs。
· 逆矩阵求解
设 A 为 n 阶非奇异矩阵,则存在初等矩阵 P1 使得 P1⋯PrA=In⇒A−1=Pr⋯P1。
可以推出求逆矩阵方法:
[A:In]→[In:A−1]
/proof/
(A:In)→(PA:P1In)→(Pr⋯P1A:Pr⋯P1In)→(In:A−1)
列变换求逆矩阵:
[AIn]
矩阵方程
AX=BX=A−1B(A∣B)→(In∣A−1B)
XA=BX=BA−1[AB]→[InBA−1]
Fragment 5 分块矩阵
· 定义
/Define/
定义1
分块矩阵: A=(aij)r×s ,分块行、分块列。
A=10−−00−−011−−00−−0∣∣−−∣∣−−∣00−−11−−000−−01−−0∣∣−−∣∣−−∣00−−00−−1=A11000A22000A33
定义2
设 A=(aij)r×s , B=(bij)k×l 为分块矩阵。
则 A=B 当且仅当 r=k , s=l , Aij=Bij ( ∀1≤i≤r , ∀1≤j≤s)
定义3
设 A=(aij)r×s , B=(bij)r×s 且 A , B 有相同的分块方式,则 A±B=(aij±bij)r×s 。
定义4
设 A=(aij)r×s , c 为常数,则 cA=(c⋅aij)r×s 。
定义5
设 A=(aij)r×s , B=(bij)s×k 且 A 的列分块与 B 的行分块方式相同。
定义 C=AB=(cij)r×k ,其中 cij=ai1b1j+ai2b2j+⋯+aisbsj 。
定义6
设 A=(Aij)r×s 则 A 的转置 AT 是 s×r 分块矩阵。
AT=A11′A12′⋮A1s′A21′A22′⋮A2s′⋯⋯⋱⋯Ar1′Ar2′⋮Ars′
定义7
设 A=(Aij)r×s,复分块矩阵 A=(Aij)r×s
我们给出如下例子
/example/
A=A1A2⋱Ak=diag(A1,A2,⋯,Ak)B=B1B2⋱Bk
Ai 与 Bi 为同阶方阵,(1≤i≤k) , AB=diag(A1B1,A2B2,⋯,AkBk)
推论:若 A1,A2,⋯,Ak 为非零矩阵,则 A 为非零矩阵。
/proof/
∣A∣=∣A1∣∣A2∣⋯∣Ak∣=0A−1=diag(A1−1,A2−1,⋯,Ak−1)
此处再给出一个例子:
设 $A = (a_{ij}){m \times n} , B = (b{ij})_{n \times p}$
A 行分块:$A = \begin{bmatrix}
\alpha_1 \
\alpha_2 \
\vdots \
\alpha_m
\end{bmatrix} , B $ 列分块:B=[β1β2⋯βp]
AB=α1α2⋮αm[β1β2⋯βp]
AB=A(β1,β2,⋯,βp)=(Aβ1,Aβ2,⋯,Aβp)
下面我们接着给出定义:
/Define/
定义8
设 A=(Aij)r×s 为分块矩阵,以下称为三类分块初等变换。
- 第一类:对换 A 的两个分块行(分块列)。
- 第二类:用非零矩阵 M 左乘 A 的某一分块行(右乘 A 的某一分块列)。
- 第三类:将一个矩阵 M 左乘 A 某一分块行加到另一个分块行上(右乘列)。
(1).Pij=Im1⋱0⋮Imi⋯⋯Imj⋮0⋱Imk(2).Pi(M)=Im1⋱M⋱Imk
(3).Pij(M)=Im1⋱Imi⋮M⋱⋯Imj⋱Imk
接着给出两个定理:
/Theorem/
定理9
(1). 分块初等阵是非零阵:
PijT=Pij,Pij(M)T=Pij(M),Tij(M)T=Tij(−M)
(2).
∣Pij∣=1,l=m+n+(m+n),i=j,m∣Pij(M)∣=∣M∣,∣Tij(M)∣=1
(3).分块初等行变换等价于左乘对应的分块初等阵。
(4). 第三类分块初等变换不改变行列式的值。
∣Mij(M)⋅A∣=∣Tij(M)⋅A∣=∣Tij(M)∣⋅∣A∣=∣A∣
定理10
设 M=[ACBD]
(1). 若 A 可逆,则 ∣M∣=∣A∣⋅∣D−CA−1B∣
(2). 若 D 可逆,则 ∣M∣=∣D∣⋅∣A−BD−1C∣
(3). 若 A 和 D 同时可逆,则 ∣M∣=∣A∣⋅∣D−CA−1B∣=∣D∣⋅∣A−BD−1C∣
针对定理10给出证明:
/proof/
(1).
M=[ACBD]→[A0BD−CA−1B]
根据定理9的(4),有
∣M∣=A0BD−CA−1B=∣A∣⋅∣D−CA−1B∣
(2). 同理
接下来给出两道例题:
/example/
A=a12an+1a1⋮an+1a1a1an+1a22⋮an+1a2⋯⋯⋱⋯a1an+1a2an+1⋮an+12
/solution/
M=−In+(aij+1)n×n=−In+a1a2⋮an11⋮1[a11a21⋯⋯an1]=−In−a1a2⋮an11⋮1(−I2)−1[a11a21⋯⋯an1]
构造 [ACBD] , A=−I , B′=C=a1a2⋮an ,D=−In
∣M∣=∣−I2∣⋅∣M∣=∣−In∣⋅−In−(a11a21⋯⋯an1)(−In)−1a1a2⋮an11⋮1
=(−1)n∑ai∑ai2−1n−1∑ai=(−1)n[(∑ai2−1)(n−1)−(∑ai)2]
Q.E.D.
/example/ 求 [ADBD] 的逆矩阵,其中 A,D 可逆。
/solution/
[AOBDImOOIn]→[AOODImO−BD−1In]→[ImOOInA−1OA−1BD−1D−1]
Fragment 6 Cauthy-Binet 公式