Logo

矩阵运算

在介绍矩阵运算之前,我们有必要先熟悉矩阵元素的表示,这是理解矩阵加法、乘法等操作的基础。此外,我们要能够掌握常见的特殊矩阵类型,因为这些概念在后续的矩阵运算中将会频繁出现。

1. 矩阵元素与特殊矩阵类型

矩阵元素表示:矩阵被定义为一个 m×n ~m\times n~的数组,其中 m ~m~是矩阵的行数,n n~是列数。矩阵中的每个元素位于特定的行和列的位置,矩阵中第 i ~i~行第 j ~j~列的元素记为 aij ~a_{ij}~,并称为矩阵 (i,j) ~(i,j)~元素。

列向量表示:矩阵 Am×n ~\mathbf{A}_{m\times n}~的每一列可以看做是一个 m ~m~维向量,矩阵 A ~\mathbf{A}~可以写成由列向量 aj (0jn) ~\mathbf{a}_{j}~\scriptsize (0\leq j\leq n)~组成的形式,表示为 A=[a1a2an] ~\mathbf{A}=\begin{bmatrix}\mathbf{a}_{1}&\mathbf{a}_{2}&\dots& \mathbf{a}_{n}\end{bmatrix}~

对角线元素表示:对角线元素是矩阵中行号和列号相同的元素,表示为 a11,a22,a33 ~a_{11},a_{22},a_{33}\dots~,它们构成 A ~\mathbf{A}~的主对角线。

开通会员解锁全部动画

对角矩阵 (Diagonal Matrix) ~(\textbf{Diagonal Matrix})~:这是一种特殊的方阵(行数等于列数) ,除对角线元素外,其余元素均为零。对角矩阵在计算中可以简化运算,其乘法操作比一般矩阵更简单且高效。此外,对角矩阵在后续涉及的求逆矩阵和特征值分解等内容中将发挥十分重要的作用。

单位矩阵:对角线上全为 1,其余元素全为 0 的方阵。单位矩阵在矩阵运算中起到类似于数字 1 的作用,任何矩阵与单位矩阵相乘都会保持自身不变。

零矩阵:所有元素都为 0 的矩阵。零矩阵在矩阵加法中起到类似于数字 0 的作用。

开通会员解锁全部动画

2. 矩阵求和与标量乘法

只有相同维数的矩阵才能求和,我们说两个矩阵相等是指它们不仅维数相同,还意味着两个矩阵的对应元素必须完全相同。我们在 1.4 向量方程一节中,介绍过向量的加法,矩阵相加可以看做是多个向量相加。 矩阵加法、标量乘法规则如下:

请看下面示例:

开通会员解锁全部动画

矩阵的和与标量乘法遵循通常的代数规则,有如下定理:

3. 矩阵乘法

我们在 1.7 线性变换一节中提到过,矩阵 B ~\mathbf{B}~和向量 x ~\mathbf{x}~的乘积 Bx ~\mathbf{B}\mathbf{x}~可看作一个线性变换的过程,矩阵 B ~\mathbf{B}~把向量 x ~\mathbf{x}~变换成了向量 Bx ~\mathbf{B}\mathbf{x}~。如果继续对 Bx ~\mathbf{B}\mathbf{x}~施加变换 A ~\mathbf{A}~,就得到了向量 A(Bx) ~\mathbf{A}(\mathbf{B}\mathbf{x})~。在这个过程中,矩阵变换 T(x)=A(Bx) ~T(\mathbf{x})=\mathbf{A}(\mathbf{B}\mathbf{x})~就是一个复合变换。现在,我们希望只用一个矩阵 C ~\mathbf{C}~就能表示这个复合变换,那么这个 C ~\mathbf{C}~如何表示呢?
复合变换 T ~T~的标准矩阵 C ~\mathbf{C}~可以用[Ab1Ab2Abp] \begin{bmatrix}\mathbf{A}\mathbf{b}_{1}&\mathbf{A}\mathbf{b}_{2}&\dots & \mathbf{A}\mathbf{b}_{p}\end{bmatrix}~来表示,我们把它记为 AB ~\mathbf{A}\mathbf{B}~,也即:
T(x)=A(Bx)=(AB)x(1)T(\mathbf{x})=\mathbf{A}(\mathbf{B}\mathbf{x})=(\mathbf{A}\mathbf{B})\mathbf{x} \tag{1}
对矩阵乘积 AB ~\mathbf{A}\mathbf{B}~的定义如下:

根据矩阵乘积 AB ~\mathbf{A}\mathbf{B}~的定义,观察下面的计算过程:

开通会员解锁全部动画

为了理解矩阵乘法的几何意义,下面的动画展示了两次复合变换的过程,注意观察 A(Bx) ~\mathbf{A}(\mathbf{B}\mathbf{x})~ (AB)x ~(\mathbf{A}\mathbf{B})\mathbf{x}~的变换结果, 其中 B ~\mathbf{B}~是关于 x2=x1 ~x_2=-x_1~对称矩阵 A ~\mathbf{A}~水平剪切矩阵

开通会员解锁全部动画

4. 计算 AB ~\mathbf{A}\mathbf{B}~的行列法则

在实际计算矩阵乘积 AB ~\mathbf{A}\mathbf{B}~时,我们通常会使用更高效的行列法则 ( Row Column Rule ) ~(~\textbf{Row Column Rule}~)~,通过逐行逐列进行乘积求和来计算矩阵乘积的每个元素。

这个规则可以从计算Ax\mathbf{A}\mathbf{x}的行-向量规则推导出来。我们用这个规则重新计算上面示例中 AB ~\mathbf{A}\mathbf{B}~的乘积, 通过观察下面计算 AB ~\mathbf{A}\mathbf{B}~的动画过程来理解这个运算规则:

开通会员解锁全部动画

观察上面的计算过程,我们还可以对行列法则进行推广,当只需要矩阵乘积的某一行时,可以直接用矩阵 A ~\mathbf{A}~的该行与矩阵 B ~\mathbf{B}~相乘,有如下结论:
rowi(AB)=rowi(A)B\text{row}_{i}(\mathbf{A}\mathbf{B})=\text{row}_{i}(\mathbf{A})\cdot \mathbf{B}
对于上面的示例,如果我们只需要计算 AB ~\mathbf{A}\mathbf{B}~的第二行,那么计算过程如下:
row2(AB)=row2(A)B=[15][436123]=[1139]\begin{align*}\text{row}_{2}(\mathbf{A}\mathbf{B})&=\text{row}_{2}(\mathbf{A})\mathbf{B}\\[1ex] &=\begin{bmatrix} 1 & -5\end{bmatrix} \begin{bmatrix}4&3&6\\1&-2&3\end{bmatrix}\\[3ex] &=\begin{bmatrix}-1& 13& -9\end{bmatrix} \end{align*}

5. 矩阵乘法的性质

接下来我们要介绍的是矩阵乘法性质,如结合律、分配律以及单位矩阵的特性,这些性质的灵活运用能够显著提升运算效率,并帮助我们设计出更加高效的算法。

根据对矩阵乘积的定义来证明,设矩阵
C=[c1c2cp]\mathbf{C}=\begin{bmatrix}\mathbf{c}_{1}&\mathbf{c}_{2}&\dots&\mathbf{c}_{p}\end{bmatrix}
由矩阵乘法定义可得:
BC=[Bc1Bc2Bcp]\mathbf{B}\mathbf{C}=\begin{bmatrix} \mathbf{B}\mathbf{c}_{1}&\mathbf{B}\mathbf{c}_{2}&\dots&\mathbf{B}\mathbf{c}_{p} \end{bmatrix}
A(BC)=[A(Bc1)ABc2A(Bcp)]\mathbf{A}(\mathbf{B}\mathbf{C})=\begin{bmatrix} \mathbf{A}(\mathbf{B}\mathbf{c}_{1})&\mathbf{A}{\mathbf{B}\mathbf{c}_{2}}&\dots&\mathbf{A}(\mathbf{B}\mathbf{c}_{p}) \end{bmatrix}
 (1) ~(1)~式可知, AB ~\mathbf{A}\mathbf{B}~的定义使得对一切 x ~\mathbf{x}~ A(Bx)=(AB)x ~\mathbf{A}(\mathbf{B}\mathbf{x})=(\mathbf{A}\mathbf{B})\mathbf{x}~,所以
A(BC)=[(AB)c1(AB)c2(AB)cp]=(AB)C\mathbf{A}(\mathbf{B}\mathbf{C})=\begin{bmatrix}(\mathbf{A}\mathbf{B})\mathbf{c}_{1} & (\mathbf{A}\mathbf{B})\mathbf{c}_{2} & \dots & (\mathbf{A}\mathbf{B})\mathbf{c}_{p}\end{bmatrix} = (\mathbf{A}\mathbf{B})\mathbf{C}

设矩阵 A ~\mathbf{A}~ m×n ~m\times n~的矩阵,矩阵 B,C ~\mathbf{B},\mathbf{C}~ n×p ~n\times p~的矩阵, 根据矩阵乘法的行列法则,矩阵 A(B+C) ~\mathbf{A}(\mathbf{B}+\mathbf{C})~ (i,j) ~(i, j)~元素可以表示为:
(A(B+C))ij=ai1(b1j+c1j)++ain(bnj+cnj)=k=1naik(bkj+ckj)\begin{align*} (\mathbf{A}(\mathbf{B}+\mathbf{C}))_{ij}&=a_{i1}(b_{1j} + c_{1j}) + \cdots + a_{in}(b_{nj} + c_{nj})\\[2ex] &=\sum_{k=1}^{n} a_{ik}(b_{kj} + c_{kj})\\[2ex] \end{align*}
运用标量乘法分配律,再把求和公式进行拆分成两项:
(A(B+C))ij=k=1naikbkj+k=1naikckj\begin{align*}(\mathbf{A}(\mathbf{B}+\mathbf{C}))_{ij} &= \sum_{k=1}^{n} a_{ik}b_{kj} + \sum_{k=1}^{n} a_{ik}c_{kj} \end{align*}
右侧的表达式 AB+AC ~\mathbf{A}\mathbf{B}+\mathbf{A}\mathbf{C}~ (i,j) ~(i,j)~元素可以分别计算:
(AB+AC)ij=(AB)ij+(AC)ij=k=1naikbkj+k=1naikckj\begin{align*}(\mathbf{A}\mathbf{B}+\mathbf{A}\mathbf{C})_{ij} &=(\mathbf{A}\mathbf{B})_{ij} + (\mathbf{A}\mathbf{C})_{ij}\\[3ex] &=\sum_{k=1}^{n} a_{ik}b_{kj}+\sum_{k=1}^{n} a_{ik}c_{kj}\\[3ex] \end{align*}
结论:由于(A(B+C))ij(\mathbf{A}(\mathbf{B}+\mathbf{C}))_{ij} (AB+AC)ij ~(\mathbf{A}\mathbf{B}+\mathbf{A}\mathbf{C})_{ij}~对应元素完全相同,因此(A(B+C))=(AB+AC)(\mathbf{A}(\mathbf{B}+\mathbf{C}))=(\mathbf{A}\mathbf{B}+\mathbf{A}\mathbf{C})成立。

证明方法同乘法左分配律。

还是使用 AB ~\mathbf{A}\mathbf{B}~的行列法则来证明。设矩阵 A ~\mathbf{A}~ m×n ~m\times n~矩阵,矩阵 B ~B~ n×p ~n\times p~矩阵, r ~r~为实数。
根据 AB ~\mathbf{A}\mathbf{B}~的行列法则 AB ~\mathbf{A}\mathbf{B}~ (i,j) ~(i,j)~元素为:
(AB)ij=k=1naikbkj(\mathbf{A}\mathbf{B})_{ij}=\sum_{k=1}^{n} a_{ik}b_{kj}
因此, r(AB) ~r(\mathbf{A}\mathbf{B})~ (i,j) ~(i,j)~元素为:
(r(AB))ij=rk=1naikbkj(r(\mathbf{A}\mathbf{B}))_{ij}=r\cdot \sum_{k=1}^{n} a_{ik}b_{kj}
 (rA)B ~(r\mathbf{A})\mathbf{B}~ (i,j) ~(i,j)~元素为:
((rA)B)ij=k=1n(raik)bkj=rk=1naikbkj((r\mathbf{A})\mathbf{B})_{ij}=\sum_{k=1}^{n} (r\cdot a_{ik})b_{kj}=r\cdot \sum_{k=1}^{n} a_{ik}b_{kj}
由此可见, r(AB)ij=((rA)B)ij ~r(\mathbf{A}\mathbf{B})_{ij}=((r\mathbf{A})\mathbf{B})_{ij}~,因此 r(AB)=(rA)B~r(\mathbf{A}\mathbf{B})=(r\mathbf{A})\mathbf{B}。同理,可证 r(AB)=A(rB) ~r(\mathbf{A}\mathbf{B})=\mathbf{A}(r\mathbf{B})~

 Im ~\mathbf{I}_{m}~是一个 m×m ~m\times m~的单位矩阵,具有如下形式:
Im=[100010001]\mathbf{I}_m = \begin{bmatrix} 1 & 0 & \cdots & 0 \\ 0 & 1 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & 1 \end{bmatrix}
根据矩阵乘法行列法则, ImA ~\mathbf{I}_{m}\mathbf{A}~ (i,j) ~(i,j)~元素为:
(ImA)ij=k=1m(Im)ikakj(\mathbf{I}_{m}\mathbf{A})_{ij}=\sum_{k=1}^{m} (\mathbf{I_m})_{ik} \cdot a_{kj}
由于 Im ~\mathbf{I}_{m}~是单位矩阵,因此 (Im)ik ~(\mathbf{I}_{m})_{ik}~ i=k ~i=k~时为 1 ~1~,其他情况为 0 ~0~。因此,求和简化为:
(ImA)ij=1aij=aij(\mathbf{I}_{m}\mathbf{A})_{ij}=1\cdot a_{ij}=a_{ij}
这表明 ImA ~\mathbf{I}_{m}\mathbf{A}~的第 i ~i~行第 j ~j~列的元素等于 A ~\mathbf{A}~的UI应的元素。因此: Im=A ~\mathbf{I}_{m}=\mathbf{A}~。同理可证得, A=AIn ~\mathbf{A}=\mathbf{A}\mathbf{I}_{n}~

从前面对矩阵乘法的定义中,我们已经知道矩阵乘法不满足交换律。这一现象可以从几何角度来理解:对一个向量施加两个不同的矩阵变换时,如果变换的顺序不同,最终得到的结果可能会不同。例如: A ~\mathbf{A}~旋转矩阵 B ~\mathbf{B}~水平拉伸矩阵,请观察下面两种变换:

开通会员解锁全部动画

当然,交换变换顺序也有可能导致相同的结果( AB=BA ~\mathbf{A}\mathbf{B}=\mathbf{B}\mathbf{A}~)。此时,我们称 A ~\mathbf{A}~ B ~\mathbf{B}~乘积可交换(前提要满足可相乘的条件 Am×nBn×m=Bn×mAm×n ~\mathbf{A}_{m\times \textcolor{red}{n}}\mathbf{B}_{\textcolor{red}{n}\times m}=\mathbf{B}_{n\times \textcolor{red}{m}}\mathbf{A}_{\textcolor{red}{m}\times n}~)。定理 2(e) ~2(e)~中单位矩阵 In ~\mathbf{I}_{n}~与任何矩阵 An×n~\mathbf{A}_{n\times n}(方阵)相乘都是可交换的。此外,两个对角矩阵相乘也是可交换的,这是因为对角矩阵的非零元素只出现在主对角线上,而这些元素之间的乘法是独立进行的,不涉及到矩阵的其他元素。例如对角矩阵:A=[3002],B=[2001.5] \mathbf{A}=\begin{bmatrix}3 & 0 \\ 0&2\end{bmatrix},\mathbf{B}=\begin{bmatrix}2 & 0 \\ 0&1.5\end{bmatrix}~相乘可交换。还有,二维平面内的旋转矩阵(更高维度不满足)相乘也是可交换的,请观察下面的变换,矩阵 A ~\mathbf{A}~表示逆时针旋转 180 ~180^\circ~,矩阵 B ~\mathbf{B}~表示逆时针旋转 135 ~135^\circ~

开通会员解锁全部动画

6. 矩阵的幂

矩阵的幂是指将一个矩阵与自身相乘若干次。从矩阵变换的角度来看,矩阵的幂表示对一个向量重复施加同一个线性变换。例如,如果矩阵 A ~\mathbf{A}~表示一个变换,那么 A2 ~\mathbf{A}^2~表示将该变换施加两次,A3 \mathbf{A}^3~表示将该变换施加三次,以此类推。矩阵幂的结果可以用来分析在多次重复变换下系统或图形的行为,例如下面的矩阵幂可以对图像进行多次重复变换:

开通会员解锁全部动画

对矩阵 A ~\mathbf{A}~进行幂运算时,通常会先尝试对矩阵 A ~\mathbf{A}~进行对角化处理,将其表示为A=PDP1\mathbf{A}=\mathbf{P}\mathbf{D}\mathbf{P}^{-1},其中 D ~\mathbf{D}~是一个对角矩阵, P ~\mathbf{P}~是一个可逆矩阵。对角化的主要作用是简化运算,使得复杂的矩阵幂运算可以转化为更容易处理的形式,从而提高计算效率。需要注意的是,并非所有矩阵都可以对角化。

7. 转置矩阵

给定一个 m×n ~m\times n~矩阵 A ~\mathbf{A}~,矩阵 A ~\mathbf{A}~转置 (Transpose) ~(Transpose)~是一个 n×m ~n\times m~的矩阵,记作 A ~\mathbf{A}^\top~,其列由矩阵 A ~\mathbf{A}~的对应行构成。请观察下面对矩阵 A ~\mathbf{A}~的转置过程:

开通会员解锁全部动画

转置矩阵有以下的代数性质:

我们还是通过矩阵 AB ~\mathbf{A}\mathbf{B}~行列法则来证明。证明过程如下:

 A ~\mathbf{A}~是一个 m×n ~m\times n~的矩阵, B ~\mathbf{B}~是一个 n×p ~n\times p~的矩阵,那么它们的乘积 AB ~\mathbf{A}\mathbf{B}~是一个 m×p ~m\times p~的矩阵。矩阵 AB ~\mathbf{A}\mathbf{B}~ (i,j) ~(i,j)~元素为:
(AB)ij=k=1naikbkj(\mathbf{A}\mathbf{B})_{ij}=\sum_{k=1}^{n} a_{ik}b_{kj}
矩阵 AB ~\mathbf{A}\mathbf{B}~的转置 (AB) ~(\mathbf{A}\mathbf{B})^\top~是一个 p×m ~p\times m~的矩阵,它的 (i,j) ~(i,j)~元素为:
(AB)ij=(AB)ji=k=1najkbki(\mathbf{A}\mathbf{B})^\top_{ij} =(\mathbf{A}\mathbf{B})_{ji} =\sum_{k=1}^{n} a_{jk}b_{ki}
再求 (BA)ij ~(\mathbf{B}^\top\mathbf{A}^\top)_{ij}~
(BA)ij=k=1n(B)ik(A)kj=k=1nbkiajk\begin{align*}(\mathbf{B}^\top\mathbf{A}^\top)_{ij} &=\sum_{k=1}^{n}\textcolor{#2196F3}{(\mathbf{B}^\top)_{ik}} \textcolor{#ff8200}{(\mathbf{A}^\top)_{kj}} &=\sum_{k=1}^{n}\textcolor{#2196f3}{b_{ki}}\textcolor{#ff8200}{a_{jk}} \end{align*}
可得,
(AB)ij=(BA)ij(\mathbf{A}\mathbf{B})^\top_{ij}=(\mathbf{B}^\top\mathbf{A}^\top)_{ij}
由于 (AB) ~(\mathbf{A}\mathbf{B})^\top~ BA ~\mathbf{B}^\top\mathbf{A}^\top~的所有 (i,j) ~(i,j)~元素相同,因此:
(AB)=BA(\mathbf{A}\mathbf{B})^\top=\mathbf{B}^\top\mathbf{A}^\top

转置矩阵是线性代数中一个重要的工具,它几乎应用于所有的矩阵操作和运算。以下是几个常见的使用场景:在某些矩阵运算中,通过将矩阵的行和列互换(即计算矩阵的转置)可以简化问题的求解;在数据统计中,获取的数据集的行列格式有时并不符合特定算法的输入要求,因此也需要对矩阵进行转置操作以匹配算法的需求;在图像处理和计算机图形学中,转置矩阵常用于表示坐标变换,例如在二维或三维空间中执行旋转、缩放等操作时,矩阵的转置能够帮助转换坐标系或执行反向变换。请观察下面的动画,从几何角度去理解转置矩阵的作用:

开通会员解锁全部动画

转置矩阵一般化的几何解释:在几何上,矩阵的转置操作相当于将线性变换的基向量映射关系“翻转”,具体地说就是将原来变换中的行向量和列向量的角色互换。一般来说,矩阵 A ~\mathbf{A}^\top~对应的几何变换通常是矩阵 A ~\mathbf{A}~的逆向变换或与之对偶的变换。例如上面第一个示例中的 A ~\mathbf{A}~是一个旋转矩阵,那么 A ~\mathbf{A}^\top~则对应反方向的旋转。

第二个示例中缩放矩阵 A ~\mathbf{A}~和它的转置矩阵 A ~\mathbf{A}^\top~在几何上有相同的作用,这时我们称矩阵 A ~\mathbf{A}~对称的。下面是对称矩阵的定义:

结合前面对角矩阵的特性,我们可以判断出所有的对角矩阵都是对称矩阵。对称矩阵具有良好的数学性质(如特征值为实数、易于分解等),在许多工程问题上应用广泛。

8. 矩阵运算在模式识别中的应用

在人工智能和计算机视觉中,模式识别是关键任务之一。接下来要介绍的应用是一个简化的模式匹配问题,用于判断 2×2 ~2\times 2~颜色块是否符合特定棋盘格模式。通过矩阵运算(转置与乘法),计算机可自动识别模式,原理可推广至更复杂的图像处理任务。

开通会员解锁全部动画

首先,将给定的 2×2 ~2\times 2~的颜色块转换为 4×1 ~4\times 1~的向量 v ~\mathbf{v}~,其中蓝色块被编码为 1 ~1~,白色块编码为 0 ~0~。再设定一个用于专门用于检测 v ~\mathbf{v}~对应的棋盘格模式的矩阵 M ~\mathbf{M}~
M=[1001010000101001]\mathbf{M} = \begin{bmatrix} 1 & 0 & 0 & -1 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ -1 & 0 & 0 & 1 \end{bmatrix}
接下来,模式识别的思路是这样的:计算 xTMx ~\mathbf{x}^T\mathbf{M} \mathbf{x} ~,如果结果不为 0 ~0~,那么 xv ~\mathbf{x} \neq \mathbf{v}~匹配失败;如果结果为 0 ~0~,且 xx0 ~\mathbf{x}^\mathbf{x} \neq 0~,那么 x=v ~\mathbf{x} = \mathbf{v}~匹配成功。

开通会员解锁全部动画