1.矩阵是什么?#
在正式讲解矩阵之前,先来讨论一下什么是矩阵:
数学上,一个 的矩阵(Matrix)是一个有 行(row) 列(column)元素的矩形阵列。矩阵里的元素可以是数字或符号甚至是函数。以数学符号而言,就是:
Ms×n(K):={数域K上s×n的阵列}当 时,该矩阵又叫做 级矩阵,也记作 。当我们说两个矩阵相等,实际上是说他们的 行数相等,列数相等,以及对应的元素相等。
定义元素都是 的矩阵叫做一个零矩阵,记作 或 。
2.矩阵的运算#
2.1.矩阵的基础运算#
我们考虑两个矩阵能否进行加法操作,对于两个形状相等的矩阵,我们很自然的定义定义:
a11⋮as1......a1n⋮asn+b11⋮bs1......b1n⋮bsn=a11+b11⋮as1+bs1......a1n+b1n⋮asn+bsn同理,对于数量乘法,有:
k⋅a11⋮as1......a1n⋮asn=ka11⋮kas1......ka1n⋮kasn易验证, 组成的集合满足线性空间条件的八条性质,故矩阵 可以构成数域 一个线性空间。
2.2.矩阵的乘法#
看一个例子:过定点 的平面上,有一向量绕 点旋转,从 点转到 点,大小都是 ,经过 角,我们记住这个旋转操作叫做 。
屏幕截图 2026-07-31 223835.png设 与 轴正向夹角为 ,则:
{x′=rcosαcosθ−rsinαsinθ=xcosθ−ysinθy′=rsinαcosθ+rcosαsinθ=xsinθ+ycosθ我们把得到的 的系数排成一张表,将旋转 角的操作 对应到这个矩阵上。变成:
θ⟼(cosθsinθ−sinθcosθ)=A而倘若向量继续绕 点旋转,转到 点,经过 角,这个旋转操作叫做 。则其也有对应的:
τ⟼(cosφsinφ−sinφcosφ)=B而将两者合起来看,考虑转角为 的映射乘积 。则有:
στ⟼(cos(θ+φ)sin(θ+φ)−sin(θ+φ)cos(θ+φ))=A⋅B展开就有:
A⋅B=(cosθcosφ−sinθsinφsinsθcosφ+cosθsinφ−sinθcosφ−cosθsinφcosθcosφ−sinθsinφ)很自然地,我们就推出来了矩阵乘法的操作。不难看出,他们大抵是按照这样的方式相乘的:
Pasted image 20260731230511.png由此引出矩阵乘法的定义:
A=a11⋮as1......a1n⋮asn,B=b11⋮bn1......b1m⋮bnm我们规定矩阵的乘法为:
A⋅B=a11b11+...+a1nbn1⋮as1b11+...+asnbn1......a11b1m+...+a1nbnm⋮as1b1m+...+asnbnms×m矩阵乘法有如下要点:
①:左矩阵的列数 = 右矩阵的行数
②:对于乘积矩阵里的元素,有:
AB(i;j)=ai1bij+ai2b2j+...+ainbnj=k=1∑naikbkj,i=1,...,sj=1,...,n③:矩阵 的行数 = 的行数 , 的列数 = 的列数
来看几个例子:
1.若 , 则 :
2.对于行向量 : ( 个) 可以看成一个 的矩阵 ,同理对于列向量 : ( 个) 也可以看成一个 的矩阵,他们之间做乘法,有:
(1,1,...,1)⋅11⋮1=(n)=n11⋮1⋅(1,1,...,1)=1⋮1⋯⋱⋯1⋮1=:J显然,我们能看到矩阵的乘法不满足交换律,两个矩阵对调位置,结果有时候会有质的变化,那么矩阵的乘法都满足什么运算法则呢?
我们上节知道了映射的乘法满足一部分运算法则,故我们猜测矩阵乘法满足以下规律:
①: (结合律)
②: (左分配律)
③: (右分配律)
④:
对于一个特殊的 矩阵:它的主对角线元素都是一,其余元素都是零,并定义为:
I=10⋮001⋮0⋯⋯⋱⋯00⋮1我们称这类矩阵叫做 单位矩阵 。
单位矩阵和任何矩阵相乘都得这个矩阵本身,即:
A⋅I=I⋅A=A⑤:
2.2.5.矩阵的幂#
我们已经定义了矩阵的乘法,多个相等的矩阵连续做乘积便构成了矩阵的幂,若矩阵为方阵,定义:
Am:=A⋅A⋅A⋅...⋅A(m个,其中m∈N∗)当 时, 。且有:
(Ak)l=Akl,Ak⋅Al=Ak+l这与代数的幂基本一致。
2.3.矩阵转置的运算#
对于矩阵的转置,常见的有以下运算:设矩阵 ,有:
①:
②:
③:
证明:*由乘法定义,有:
(AB)′(i;j)=AB(j;i)=k=1∑najkbki且:
(B′A′)(i;j)=k=1∑nB′(i;k)A′(k;j)=k=1∑nbkiajk因此二者相等,证毕。*
了解了以上运算,对于一个普通的线性方程组:
⎩⎨⎧a11x1+a12x2+...+a1nxn=b1...an1x1+an2x2+...+annxn=bs(1)能写成:
a11⋮as1⋯⋱⋯a1n⋮asn⋅x1⋮xn=b1⋮bn的形式,我们用 分别代表这三个矩阵,则,原本形式复杂的线性方程组就可以表示成:
AX=β相应的齐次线性方程组就是:
AX=0这也意味着我们也可以用矩阵的理论去解线性方程组。
2.3.乘积的秩#
我们回过头来再看矩阵乘积:
A⋅B=a11⋮as1......a1n⋮asnb11⋮bn1......b1m⋮bnm我们把A矩阵的每一个列向量从左到右记作: ,现在我们考虑这个 的第 列是多少,经过计算应该能看出第 列为:
a11b1j+...+a1nbnj⋮as1b1j+...+asnbnj=b1ja11⋮as1+...+bnja1n⋮asn=b1jα1+...+bnjαn故:
A⋅B=(α1,α2,,...,αn)b11⋮bn1......b1m⋮bnm=(b11α1+...+bn1αn,...,b1mα1+...+bnmαn)观察可以发现我们实际上就是在用这个行向量 去乘了矩阵 的每一列.而且 的列向量组可以由 的列向量组线性表出。于是:
rank(A⋅B)≤rank(A)而 ,所以:
rank(A⋅B)≤rank(B)我们总结一下:
3.特殊矩阵#
3.1.基本矩阵/单元素矩阵#
只有一个元素为一,位于第 行第 列,而其余元素都是 的矩阵,叫做一个基本矩阵/单元素矩阵,记作:
Eij:=00⋮001⋮0⋯⋯⋱⋯00⋮0由矩阵的加法,我们能知道任何一个矩阵 :
A=a11⋮as1......a1n⋮asn都可以表示成诸如:
a11E11+...+a1nE1n+...+asnEsn的形式,而且显然可知道 线性无关。所以以矩阵 中每一个元素作为基本矩阵元素的位置构成的这一堆矩阵,就是矩阵 的一个基,且
dimMs×n(K)=s×n特别地,
dimMn(K)=n23.2.对角矩阵#
除主对角线之外的元素皆为 的方形矩阵叫对角矩阵(英语:diagonal matrix)。对角线上的元素可以为 或其他值。记作
D=d1d2⋱dn数域 K 上 n 级对角矩阵组成的集合是 的一个子空间 ,对于一般的对角矩阵,有:
d1⋱ds⋅γ1⋮γs=d1γ1⋮dsγs其中 是某个矩阵的行向量组。还有:
(α1,...,αn)⋅d1⋱dn=(d1α1,...,dnαn)两个对角矩阵相乘结果仍旧是一个对角矩阵,具体有:
d1⋱dn⋅c1⋱cn=d1c1⋱dncn3.2.5.数量矩阵#
主对角线元素都相等的对角矩阵叫做数量矩阵。又叫做纯量阵。记作:
k⋅I:=kk⋱k取满数域 上所有数的数量矩阵构成的集合也是 的一个子空间 。
对于任意一个 级矩阵 :
A⋅(k⋅I)=k⋅(A⋅I)=k⋅A(k⋅I)⋅A=k⋅(I⋅A)=k⋅A因此:
(k⋅I)⋅A=A⋅(k⋅I)也就是说,虽然矩阵与矩阵之间的乘法并不满足交换律,但是 数量矩阵与矩阵是的确满足交换率 的。
3.3.三角矩阵#
三角矩阵(英语:triangular matrix)因其非零系数的排列呈三角形状而得名。三角矩阵分上三角矩阵和下三角矩阵两种。
3.3.1.上三角矩阵#
主对角线左下方的系数全部为零的n级矩阵叫做n级上三角矩阵,常记作:
U=a110⋮0a12a22⋮0⋯⋯⋱⋯a1na2n⋮ann同理,他组成的集合也是 的一个子空间 。对于上三角矩阵的乘法:
a110⋮0a12a22⋮0⋯⋯⋱⋯a1na2n⋮annb110⋮0b12b22⋮0⋯⋯⋱⋯b1nb2n⋮bnn=a11b110⋮0∗a22b22⋮0⋯⋯⋱⋯∗∗⋮annbnn不细考虑其他元素,上三角矩阵相乘结果还是上三角矩阵,其中主对角线元素为两上三角矩阵元素乘积。
3.3.2.下三角矩阵#
有上三角矩阵,类似的也有下三角矩阵:
L=a11a21⋮am10a22⋮am2⋯⋯⋱⋯00⋮amn同理,他组成的集合也是 的一个子空间 。对于下三角矩阵的乘法:
a11a21⋮am10a22⋮am2⋯⋯⋱⋯00⋮amnb11b21⋮bm10b22⋮bm2⋯⋯⋱⋯00⋮bmn=a11b11∗⋮∗0a22b22⋮∗⋯⋯⋱⋯00⋮annbnn不细考虑其他元素,下三角矩阵相乘结果还是下三角矩阵,其中主对角线元素为两下三角矩阵元素乘积。
3.4.初等矩阵#
线性代数中,初等矩阵(又称为基本矩阵)是一个与单位矩阵只有微小区别的矩阵。常记作 P 具体来说,一个 阶单位矩阵 经过一次初等行变换或一次初等列变换所得矩阵称为 阶初等矩阵。
同理,他组成的集合也是 的一个子空间 。
初等矩阵具体有以下几类:
考虑单位矩阵 :
I=1⋱1⋱1⋱1①:我们将其中第 行的 倍加到第 行,(或者把第 列的 倍加到第 列)即记作:
P(j,i(k))=1⋱1⋮k⋱⋯1⋱1②:我们将其中的第 行 和 第 行(或第 列 和 第 列)互换,即记作
P(i,j)=1⋱01⋱10⋱1③:我们将其中的第 行(或第 列)乘上一个常数 ,即记作
P(i(c))=1⋱c⋱1⋱1不难发现,初等矩阵即是将我们已知的 3 种初等变换应用于一单位矩阵的结果。
接下来考虑初等矩阵左乘普通矩阵会发生什么:我们将矩阵 以行向量组的方式表出,即:
P(j,i(k))⋅γ1⋮γi⋮γj⋮γs=1⋱1⋮k⋱⋯1⋱1⋅γ1⋮γi⋮γj⋮γs=γ1⋮γi⋮kγi+γj⋮γs也就是说,将一个矩阵 左乘一个①型初等矩阵,就可以把矩阵 的第 行的 倍加到第 行上,相当于做了对应的初等行变换,效果是一样的。,当然其余型同理。
看完左乘,我们继续看右乘:我们将矩阵 以列向量组的方式表出,即:
(α1,...,αi,...,αj,...,αn)⋅1⋱1⋮k⋱⋯1⋱1=(α1,...,αi+kαj,...,αj,...,αn)也就是说,将一个矩阵 右乘一个①型初等矩阵,就可以把矩阵 的第 列的 倍加到第 列上,相当于做了对应的初等列变换,效果是一样的。,当然其余型同理。
于是我们可以总结:将矩阵 A 左乘一包含某初等变换的初等矩阵,其自身将做对应初等矩阵的初等 ==行将矩阵 A 左乘一包含某初等变换的初等矩阵,其自身将做对应初等矩阵的初等 ==行 变换;若右乘一包含某初等变换的初等矩阵,其自身将做对应初等矩阵的初等 列 变换。*
Notice:对于初等矩阵,有:
P(j,i(k))−1=P(j,i(−k)),P(i,j)−1=P(i,j),P(i(c))−1=P(i(c−1))逆矩阵的定义请看下方第四节.
3.5 对称矩阵#
对于数域 上的 级矩阵 :若有
AT=A则称这个矩阵叫做一个对称矩阵,根据转置的定义,显然应该有:。化成矩阵形式,即是形如:
a11a12⋮a1na12a22⋮a2n⋯⋯⋱⋯a1na2n⋮ann的矩阵,显然能看出来对称矩阵中的右上至左下方向元素以主对角线为轴进行对称。
同理,他组成的集合也是 的一个子空间 。
对称矩阵具有如下性质:
①:对于任何 级方形矩阵 , 是对称矩阵。
证: ,
②:对角矩阵都是对称矩阵。
③:当且仅当两者的乘法可交换(即 )时,两个对称矩阵的积()是对称矩阵。
3.5.5 斜对称矩阵#
对于数域 上的 级矩阵 :若有
AT=−A则称这个矩阵叫做一个斜对称矩阵/反对称矩阵,根据转置的定义,显然应该有:。特别地,若 , ,化成矩阵形式,即是形如:
0−a12⋮−a1na120⋮−a2n⋯⋯⋱⋯a1na2n⋮0的矩阵,显然能看出来在非偶数域中,斜对称矩阵中的主对角线元素皆为 。
同理,他组成的集合也是 的一个子空间 。
4.矩阵的逆#
4.1 可逆矩阵#
类比我们之前所述的映射的逆,我们定义矩阵的逆为:
设 ,如果存在 ,使得:
A⋅B=B⋅A=I我们就称 是一个可逆矩阵, 称为 的逆矩阵,常记作。即有:
A⋅A−1=A−1⋅A=I任意矩阵的逆矩阵是唯一的。
由逆矩阵的定义,有: 。注意,只有方阵(n×n的矩阵) 才可能有逆矩阵。若方阵 的逆矩阵存在,则称 为非奇异方阵或可逆方阵。
4.2 求可逆矩阵的伴随矩阵法#
设一矩阵 可逆, 则有: ,所以有 ,即 。反过来说,设 ,我们要找到一个矩阵,使得它能等于单位矩阵,则我们肯定有:
a11a21⋮an1a12a22⋮an2⋯⋯⋱⋯a1na2n⋮ann⋅A11A12⋮A1nA21A22⋮A2n⋯⋯⋱⋯An1An2⋮Ann=∣A∣0⋮00∣A∣⋮0⋯⋯⋱⋯00⋮∣A∣其中: 是元素 的代数余子式,由于 ,所以有上面的代数关系。我们把找到的这个矩阵称作矩阵 的 伴随矩阵 ,记作 或 。则显然,我们有如下关系:
A⋅A∗=∣A∣⋅I类似可得:
A∗⋅A=∣A∣⋅I由于 ,因此:
A(∣A∣1A∗)=I=(∣A∣1A∗)A从而 A 可逆,且
A−1=∣A∣A∗我们可以总结为以下结论:
数域 上的 级矩阵 可逆 满秩,且当 A 可逆时:
A−1=∣A∣A∗设 ,如果我们现在只知道 ,由先前证明可知 ,故 满秩,所以 可逆,同理, 也可逆,由于 ,两边同时左乘 ,有
B=A−1同理:也能得到:
A=B−14.2.5 伴随矩阵补充#
在线性代数中,一个方形矩阵的伴随矩阵(英语:adjugate matrix)是一个类似于逆矩阵的概念。如果矩阵可逆,那么它的逆矩阵和它的伴随矩阵之间只差一个系数。然而,伴随矩阵对不可逆的矩阵也有定义,并且不需要用到除法。
A的伴随矩阵是一个 的矩阵(记作 ),使得其第 行第 列的元素是 A 关于第 j 行第 i 列的代数余子式。 简言之,伴随矩阵就是把原来余子矩阵C每一列的代数余子式横着写.
最简单的例子,有:
A=(acbd),A∗=(d−c−ba)4.3 可逆矩阵的性质#
若 都是可逆矩阵,则:
①:
②:
③:
推广:
④: 若 可逆,则 可逆,且
⑤:
⑥: 级矩阵 可逆,则其一定可以通过初等行变换化成单位矩阵
4.4 求可逆矩阵的初等变换法#
由刚刚描述的性质 ⑥ ,对于 级矩阵 ,这两个矩阵可以只经由初等行变换,或者只经由初等列变换,变为单位矩阵。
因为对矩阵 施以初等行变换(初等列变换)就相当于在 的左边(右边)乘以相应的初等矩阵,即
(Pt...P2P1)A=I所以有:
A=(Pt...P2P1)−1=P1−1P2−1...Pt−1即: 级矩阵 可逆 等于一些初等矩阵的乘积,证明如上
又设 , 是 级的可逆矩阵,由于:,因此:
rank(PA)=rank[(P1P2...Pn)A]=rank(A)设 是 级可逆矩阵,有:
rank(AQ)=rank[(AQ)T]=rank(QTAT)=rank(AT)=rank(A)即:用可逆矩阵左(右)乘矩阵 ,不改变矩阵 的秩,证明如上
若 可逆,则
(Pt...P2P1)A=IA−1=Pt...P2P1I这两个式子就意味着:我们可以同时对 和 施以相同的初等行变换(初等列变换)。这样,当矩阵 被变为 时, 就被变为 的逆阵。
这是我们最常用的求逆矩阵的最基本的方法,称为初等变换法。
5.分块矩阵#
一个分块矩阵或是分段矩阵就是将矩阵分割出较小的矩形矩阵,这些较小的矩阵就称为区块。
换个方式来说,就是以较小的矩阵组合成一个矩阵。分块矩阵的分割原则是以水平线和垂直线进行划分。
设一个矩阵 ,我们将其按行分为 个,按列分为 个,则就有:
a11a21⋮am1a12a22⋮am2⋯⋯⋱⋯a1na2n⋮amn=A11⋮Au1⋯⋱⋯A1t⋮Aut其中:
A11=a11⋮as11⋯⋱⋯a1n1⋮as1n1,...若我们再设矩阵 ,则:
AB=A11⋮Au1⋯⋱⋯A1t⋮AutnB11⋮Bt1⋯⋱⋯B1v⋮Btv每个分块都像普通矩阵乘法那样运算,在乘完之后分成了许多的小矩阵乘法,继续算每个小矩阵乘法的结果加起来就是对应分块矩阵乘完对应位置的结果。
5.1.分块矩阵的初等行变换#
①:把第一行左乘可逆矩阵 P 加到第二行
(A1A3A2A4)把第一行左乘可逆矩阵P加到第二行→(A1PA1+A3A2PA2+A4)②:互换两个块行位置。
③:用一个可逆矩阵 P 左乘某一个块行。
类似的,也有初等列变换,不过需要将左乘变为右乘。
对于我们之前讨论的初等矩阵,我们可以尝试将它分块,一边更好地表达:
考虑单位矩阵 :
I=1⋱1⋱1⋱1将它化作分块矩阵:
(IOOI)①:将其中第 块行的 倍加到第 块行,(或者把第 块列的 倍加到第 块列)可以化为:
(IPOI)②:将其中的第 块行 和 第 块行(或第 块列 和 第 块列)互换,可以化为
(OIIO)③:将其中的第 块行(或第 块列)乘上一个常数 ,可以化为
(POOI)这些化为分块矩阵的初等矩阵叫分块可逆矩阵。且不难看出它们都是可逆矩阵,它们可以用于分块矩阵的初等行变换,这也间接解释了为什么需要左乘。根据先前我们所讲的知识,很容易也能知道:
分块矩阵的初等行(列)变换不改变矩阵的秩
6. 柯西-比内公式#
6.1 引入和证明#
∣AB∣=0若 ,那么 等于 的所有 阶子式以及 的相应 阶子式的乘积之和,即
∣AB∣=1≤v1<...<vs≤n∑A(1,...,sv1,...,vs)B(v1,...,vs1,...,s)证明:
情景①:若 ,,于是 不满秩,从而 。
情景②:若 ,考虑:
(AIOB)把其的第二块行乘以 -A 加到第一块行,变成:
(OI−ABB)从而:
(IsOOIn)(AIOB)=(OI−ABB)情景②.①: 时
若 ,则,从而 ;
若 ,则 可逆,
A=P1P2...Pm→AB=(P1P2...Pm)B其中:
∣P(j,i(k))B∣=∣B∣,∣P(i,j)B∣=−∣B∣,∣P(i(c))B∣=c∣B∣特别地,当 时,从上述式子得:
∣P(j,i(k))∣=∣I∣=1,∣P(i,j)∣=−∣I∣=−1,∣P(i(c))∣=c∣I∣=c于是:
∣P(j,i(k))B∣=∣P(j,i(k))∣⋅∣B∣,∣P(i,j)B∣=∣P(i,j)∣⋅∣B∣,∣P(i(c))B∣=∣P(i(c))∣⋅∣B∣于是:
∣AB∣=∣P1∣⋅∣P2∣⋅...⋅∣Pm∣⋅∣B∣特别地,当 时:
∣A∣=∣AI∣=∣P1∣⋅∣P2∣⋅...⋅∣Pm∣所以:
∣AB∣=∣A∣⋅∣B∣情景②.②: 时:依旧考虑:
IsOOInAIOB=OI−ABB前者为 1,从而:
AIOB=OI−ABB将后面那个行列式按第一块行(设第一块行占了 行,剩余占了 行)展开:
OI−ABB=∣−AB∣⋅(−1)(1+...+s)+[(n+1)+...+(n+s)]⋅∣I∣=(−1)s+sn∣AB∣再考虑前面那个行列式(设第一块行占了 行, 列,剩余占了 行, 列):
AIOB=a11⋮as11⋮0⋯⋯⋯⋯a1n⋮asn0⋮10⋮0b11⋮bn1⋯⋯⋯⋯0⋮0b1s⋮bns设单位子矩阵的列向量记作:我们按前 行展开:
1≤v1<...<vs≤n∑A(1,...,sv1,...,vs)(−1)(1+...+s)+(v1+...+vs)∣(εμ1,...,εμn−s,B)∣其中令 且 。化简得:
1≤v1<...<vs≤n∑A(1,...,sv1,...,vs)(−1)(1+...+s)+(v1+...+vs)∣In−s∣B(v1,...,vs1,...,s)(−1)(u1+...+un−s)+[1+...+(n−s)]=1≤v1<...<vs≤n∑A(1,...,sv1,...,vs)B(v1,...,vs1,...,s)(−1)s2+n2−ns+n将刚刚两个式子带入我们化成的这几个式子,先看 项:
(−1)s+sn(−1)s2+n2−ns+n=(−1)s2+n2−ns+n−s−sn=(−1)s2+n2+n−s=(−1)s(s−1)+n(n+1)=1刚好为一,消掉了,再把其余项组合:
∣AB∣=1≤v1<...<vs≤n∑A(1,...,sv1,...,vs)B(v1,...,vs1,...,s)证毕。
6.2 正式介绍#
线性代数中,柯西-比内公式(Cauchy–Binet formula)将行列式的可乘性(两个方块矩阵的行列式等于两个行列式的乘积)推广到非方块矩阵。
这条公式的意义在于将两个不规则的矩阵乘积变成规则方阵的乘积之和,在理论上与其它公式证明中具有巨大意义。