M04 · 第 2 章 · 第一编 向量、矩阵与方程组

矩阵及其运算:从行列结构到正交变换

从矩阵的形状、索引和分块结构出发,以列的线性组合和映射复合推导矩阵乘法,再建立单位矩阵、逆矩阵与正交矩阵的基本性质。

报告页面错误
预备知识向量、坐标与线性组合线性组合

本章目标

  1. 用形状、行列索引和轴的语义准确描述矩阵。
  2. 计算矩阵加法、数乘、转置和分块运算,并检查尺寸相容性。
  3. 从列的线性组合与映射复合推导矩阵向量乘法和矩阵乘法。
  4. 解释单位矩阵与逆矩阵的作用,证明逆的唯一性并用双侧乘法核验候选逆。
  5. 判定低阶正交矩阵,并证明其保持内积、长度和夹角。
页面阅读位置0% · 仅保存在此浏览器
章节未开始
本册完成进度0/6 章 · 0%
本页目录

一张数表为何需要行列语义

矩阵看起来是一张长方形数表,但行和列从来不是纯排版。若四名学生各完成三次测验,可以约定每行对应一名学生、每列对应一次测验;若把约定反过来,同一批数字就具有不同的输入输出结构。计算前必须说明每条轴枚举什么对象,否则形状即使匹配,结果也可能没有合理含义。

本章使用大写字母表示矩阵、粗体小写字母表示列向量。若矩阵 AAmm 行、nn 列,写作 ARm×nA\in\mathbb R^{m\times n}。第 ii 行第 jj 列元素记为 aija_{ij}。索引范围始终是 1im1\le i\le m1jn1\le j\le n

矩阵、形状与相等

一个 m×nm\times n 实矩阵是按 mm 行、nn 列排列的实数数组

A=[a11a12a1na21a22a2nam1am2amn].A= \begin{bmatrix} a_{11}&a_{12}&\cdots&a_{1n}\\ a_{21}&a_{22}&\cdots&a_{2n}\\ \vdots&\vdots&\ddots&\vdots\\ a_{m1}&a_{m2}&\cdots&a_{mn} \end{bmatrix}.

两个矩阵相等,当且仅当形状相同,并且每个对应位置的元素相等。

矩阵的第 jj 列是 Rm\mathbb R^m 中的向量,第 ii 行则记录输出第 ii 个坐标如何读取全部输入坐标。方阵的行数与列数相等;对角矩阵在主对角线之外全为零;零矩阵必须连同形状说明,因为 2×32\times3 零矩阵与 3×23\times2 零矩阵不是同一对象。

加法数乘与转置

矩阵加法、数乘与转置

同形矩阵按对应元素相加,标量乘法也逐元素进行:

(A+B)ij=aij+bij,(cA)ij=caij.(A+B)_{ij}=a_{ij}+b_{ij}, \qquad (cA)_{ij}=ca_{ij}.

转置交换行列。若 ARm×nA\in\mathbb R^{m\times n},则 ATRn×mA^\mathsf T\in\mathbb R^{n\times m},并且

(AT)ij=aji.(A^\mathsf T)_{ij}=a_{ji}.

不同形状的矩阵不能相加,因为并非每个位置都有对应元素。转置也不是把图形顺时针旋转九十度;它交换两个索引的角色。由定义逐项可得

(A+B)T=AT+BT,(cA)T=cAT,(AT)T=A.(A+B)^\mathsf T=A^\mathsf T+B^\mathsf T, \qquad (cA)^\mathsf T=cA^\mathsf T, \qquad (A^\mathsf T)^\mathsf T=A.
从数据语义检查转置

四名学生的三次成绩排列为

X=[807291667578928895708184]R4×3.X= \begin{bmatrix} 80&72&91\\ 66&75&78\\ 92&88&95\\ 70&81&84 \end{bmatrix} \in\mathbb R^{4\times3}.

这里每行是一名学生,每列是一次测验。XTX^\mathsf T 的形状为 3×43\times4,其第一行 (80,66,92,70)(80,66,92,70) 汇总第一次测验的四个结果。转置没有改变任何成绩,只改变了读取数据的轴顺序。

对称部分与反对称部分

转置不仅交换索引,还能把方阵中“沿主对角线成对出现”的信息分开。若 A=ATA=A^\mathsf T,称 AA 为对称矩阵;若 AT=AA^\mathsf T=-A,称 AA 为反对称矩阵。实反对称矩阵的对角元必为零,因为 aii=aiia_{ii}=-a_{ii}

方阵的对称—反对称分解

每个实方阵 AA 都能唯一写成

A=S+K,A=S+K,

其中 SS 对称,KK 反对称。具体地,

S=A+AT2,K=AAT2.S=\frac{A+A^\mathsf T}{2}, \qquad K=\frac{A-A^\mathsf T}{2}.
证明

由转置的运算律,

ST=AT+A2=S,KT=ATA2=K,S^\mathsf T =\frac{A^\mathsf T+A}{2}=S, \qquad K^\mathsf T =\frac{A^\mathsf T-A}{2}=-K,

并且 S+K=AS+K=A,所以这种分解确实存在。若另有 A=S1+K1A=S_1+K_1,其中 S1S_1 对称、K1K_1 反对称,那么对等式取转置并与原式相加,得到 2S1=A+AT2S_1=A+A^\mathsf T;相减则得到 2K1=AAT2K_1=A-A^\mathsf T。因此 S1=SS_1=SK1=KK_1=K,分解唯一。

例如

A=[2514]=[2224]S+[0330]K.A=\begin{bmatrix}2&5\\-1&4\end{bmatrix} =\underbrace{\begin{bmatrix}2&2\\2&4\end{bmatrix}}_{S} +\underbrace{\begin{bmatrix}0&3\\-3&0\end{bmatrix}}_{K}.

对任意实列向量 x\mathbf x,标量 xTKx\mathbf x^\mathsf TK\mathbf x 取转置后等于它自身,也等于 xTKx-\mathbf x^\mathsf TK\mathbf x,故只能为零。因此二次表达式 xTAx\mathbf x^\mathsf TA\mathbf x 只读取 AA 的对称部分。以 x=(1,2)T\mathbf x=(1,2)^\mathsf T 为例, xTAx=26\mathbf x^\mathsf TA\mathbf x=26xTSx=26\mathbf x^\mathsf TS\mathbf x=26,而 xTKx=0\mathbf x^\mathsf TK\mathbf x=0

分块结构保留原有运算

大矩阵可沿相容的行列边界分成若干块。设

A=[A11A12A21A22],x=[x1x2].A= \begin{bmatrix} A_{11}&A_{12}\\ A_{21}&A_{22} \end{bmatrix}, \qquad \mathbf x= \begin{bmatrix} \mathbf x_1\\ \mathbf x_2 \end{bmatrix}.

只要块的尺寸相容,就有

Ax=[A11x1+A12x2A21x1+A22x2].A\mathbf x= \begin{bmatrix} A_{11}\mathbf x_1+A_{12}\mathbf x_2\\ A_{21}\mathbf x_1+A_{22}\mathbf x_2 \end{bmatrix}.

分块没有创造新乘法,只是把原来的求和按索引区间分组。块的位置相同不代表形状相同,块相加仍须逐项检查尺寸。

用分块分离两个变量组

A=[102011320],x=[425].A= \begin{bmatrix} 1&0&2\\ 0&1&-1\\ 3&2&0 \end{bmatrix}, \qquad \mathbf x= \begin{bmatrix} 4\\ -2\\ 5 \end{bmatrix}.

把前两个坐标作为 x1=(4,2)T\mathbf x_1=(4,-2)^\mathsf T,最后一个坐标作为 x2=5x_2=5。相应分块为

A11=[1001],A12=[21],A21=[32],A22=[0].A_{11}= \begin{bmatrix}1&0\\0&1\end{bmatrix}, \quad A_{12}= \begin{bmatrix}2\\-1\end{bmatrix}, \quad A_{21}=\begin{bmatrix}3&2\end{bmatrix}, \quad A_{22}=\begin{bmatrix}0\end{bmatrix}.

于是上方输出为 A11x1+A12x2=(14,7)TA_{11}\mathbf x_1+A_{12}x_2=(14,-7)^\mathsf T, 下方输出为 A21x1+A22x2=8A_{21}\mathbf x_1+A_{22}x_2=8。直接逐行计算同样得到 (14,7,8)T(14,-7,8)^\mathsf T

分块矩阵彼此相乘时,仍然执行“行乘列”,只是其中的元素换成了矩阵块。若

A=[A11A12A21A22],B=[B11B12B21B22],A=\begin{bmatrix}A_{11}&A_{12}\\A_{21}&A_{22}\end{bmatrix}, \qquad B=\begin{bmatrix}B_{11}&B_{12}\\B_{21}&B_{22}\end{bmatrix},

并且每个块乘积都有定义,则

AB=[A11B11+A12B21A11B12+A12B22A21B11+A22B21A21B12+A22B22].AB= \begin{bmatrix} A_{11}B_{11}+A_{12}B_{21}&A_{11}B_{12}+A_{12}B_{22}\\ A_{21}B_{11}+A_{22}B_{21}&A_{21}B_{12}+A_{22}B_{22} \end{bmatrix}.

这里的四个加号也要求两侧块同形。作为复算,取

A11=[1021],A12=[11],A21=[30],A22=[2],A_{11}=\begin{bmatrix}1&0\\2&1\end{bmatrix},\quad A_{12}=\begin{bmatrix}1\\-1\end{bmatrix},\quad A_{21}=\begin{bmatrix}3&0\end{bmatrix},\quad A_{22}=\begin{bmatrix}2\end{bmatrix},

以及

B11=[0110],B12=[21],B21=[11],B22=[3].B_{11}=\begin{bmatrix}0&1\\1&0\end{bmatrix},\quad B_{12}=\begin{bmatrix}2\\1\end{bmatrix},\quad B_{21}=\begin{bmatrix}1&-1\end{bmatrix},\quad B_{22}=\begin{bmatrix}3\end{bmatrix}.

四个结果块依次为

[1003],[52],[21],[12].\begin{bmatrix}1&0\\0&3\end{bmatrix}, \quad \begin{bmatrix}5\\2\end{bmatrix}, \quad \begin{bmatrix}2&1\end{bmatrix}, \quad \begin{bmatrix}12\end{bmatrix}.

拼回普通矩阵得到 AB=[1050322112]AB=\begin{bmatrix}1&0&5\\0&3&2\\2&1&12\end{bmatrix}。这一结果也可用九次行列配对逐项核对;分块的价值在于保留变量组之间的结构,而不是改变答案。

矩阵乘向量是列的线性组合

ARm×nA\in\mathbb R^{m\times n} 的列依次为 a1,,an\mathbf a_1,\ldots,\mathbf a_n,输入 x=(x1,,xn)T\mathbf x=(x_1,\ldots,x_n)^\mathsf T。定义

Ax=x1a1++xnan.A\mathbf x =x_1\mathbf a_1+\cdots+x_n\mathbf a_n.

输出属于 Rm\mathbb R^m。同一运算从第 ii 行看是

(Ax)i=j=1naijxj.(A\mathbf x)_i =\sum_{j=1}^{n}a_{ij}x_j.

“列的线性组合”和“每行读取输入”并非两种算法,而是同一乘法的两种解释。前者说明输出只能由矩阵列合成,后者便于逐坐标复算。

用列和行交叉核对

A=[121320],x=(2,1)T.A= \begin{bmatrix} 1&2\\ -1&3\\ 2&0 \end{bmatrix}, \qquad \mathbf x=(2,-1)^\mathsf T.

按列计算,

Ax=2(1,1,2)T(2,3,0)T=(0,5,4)T.A\mathbf x =2(1,-1,2)^\mathsf T-(2,3,0)^\mathsf T =(0,-5,4)^\mathsf T.

按行计算则得到 (12+2(1),12+3(1),22+0(1))T(1\cdot2+2\cdot(-1),-1\cdot2+3\cdot(-1),2\cdot2+0\cdot(-1))^\mathsf T, 仍为 (0,5,4)T(0,-5,4)^\mathsf T。输出有三个分量,也与 3×23\times2 矩阵乘二维向量的形状判断一致。

从左右两侧读取同一矩阵

右乘列向量 AxA\mathbf x 使用 x\mathbf x 的分量组合 AA 的各列;左乘行向量 yTA\mathbf y^\mathsf TA 则使用 y\mathbf y 的分量组合 AA 的各行。两种读取方式在标量 yTAx\mathbf y^\mathsf TA\mathbf x 中相遇:既可以先算 AxA\mathbf x 再与 y\mathbf y 做内积,也可以先算 yTA\mathbf y^\mathsf TA 再作用于 x\mathbf x

沿用上例的 A=[121320]A=\begin{bmatrix}1&2\\-1&3\\2&0\end{bmatrix}x=(2,1)T\mathbf x=(2,-1)^\mathsf T,再取 y=(1,2,1)T\mathbf y=(1,-2,1)^\mathsf T。按行组合得到

yTA=(1,2)2(1,3)+(2,0)=(5,4).\mathbf y^\mathsf TA =(1,2)-2(-1,3)+(2,0) =(5,-4).

因此 (yTA)x=52+(4)(1)=14(\mathbf y^\mathsf TA)\mathbf x=5\cdot2+(-4)(-1)=14。另一方面,上例已有 Ax=(0,5,4)TA\mathbf x=(0,-5,4)^\mathsf T,故 yT(Ax)=0+10+4=14\mathbf y^\mathsf T(A\mathbf x)=0+10+4=14。相等并非巧合,而是有限求和可以重新分组;这也是结合律在行向量、矩阵和列向量之间的具体表现。

矩阵乘法来自映射复合

矩阵乘法

ARm×nA\in\mathbb R^{m\times n}BRn×pB\in\mathbb R^{n\times p}。乘积 ABRm×pAB\in\mathbb R^{m\times p} 的元素定义为

(AB)ij=k=1naikbkj.(AB)_{ij}=\sum_{k=1}^{n}a_{ik}b_{kj}.

等价地,ABAB 的第 jj 列是 AA 乘以 BB 的第 jj 列。

BB 先把 pp 维输入变成 nn 维中间量,AA 再把中间量变成 mm 维输出,所以

A(Bx)=(AB)x.A(B\mathbf x)=(AB)\mathbf x.

内维 nn 对应被求和的中间坐标,外侧的 m,pm,p 成为结果形状。这个复合解释同时说明乘法顺序:在 ABxAB\mathbf x 中,右侧的 BB 先作用。

矩阵乘法对加法的分配律

只要形状相容,就有

A(B+C)=AB+AC,(A+B)C=AC+BC,A(B+C)=AB+AC, \qquad (A+B)C=AC+BC,

并且对标量 ccc(AB)=(cA)B=A(cB)c(AB)=(cA)B=A(cB)

证明

以第一式为例,第 (i,j)(i,j) 个元素为

kaik(bkj+ckj)=kaikbkj+kaikckj,\sum_k a_{ik}(b_{kj}+c_{kj}) =\sum_k a_{ik}b_{kj}+\sum_k a_{ik}c_{kj},

这正是 (AB+AC)ij(AB+AC)_{ij}。第二式同理;标量可以移入有限和中的任意因子,所以三种数乘写法也相等。

分配律允许把共同变换提取出来,也允许拆开输入贡献。例如 A(2u3v)=2Au3AvA(2\mathbf u-3\mathbf v)=2A\mathbf u-3A\mathbf v,说明矩阵作用与向量的线性组合相容。

矩阵乘法的结合律

只要各乘积形状相容,就有

(AB)C=A(BC).(AB)C=A(BC).
证明

ARm×nA\in\mathbb R^{m\times n}BRn×pB\in\mathbb R^{n\times p}CRp×qC\in\mathbb R^{p\times q}。比较第 (i,j)(i,j) 个元素:

((AB)C)ij==1p(AB)icj==1pk=1naikbkcj=k=1naik=1pbkcj=(A(BC))ij.\begin{aligned} ((AB)C)_{ij} &=\sum_{\ell=1}^{p}(AB)_{i\ell}c_{\ell j}\\ &=\sum_{\ell=1}^{p}\sum_{k=1}^{n}a_{ik}b_{k\ell}c_{\ell j}\\ &=\sum_{k=1}^{n}a_{ik}\sum_{\ell=1}^{p}b_{k\ell}c_{\ell j}\\ &=(A(BC))_{ij}. \end{aligned}

有限求和可以交换次序,因此所有对应元素相等。

结合律允许选择计算次序,却不允许交换矩阵次序。结合改变括号,交换则改变复合顺序。

缩放与剪切的顺序不同

A=[1101],B=[2001].A=\begin{bmatrix}1&1\\0&1\end{bmatrix}, \qquad B=\begin{bmatrix}2&0\\0&1\end{bmatrix}.

直接计算得

AB=[2101],BA=[2201].AB=\begin{bmatrix}2&1\\0&1\end{bmatrix}, \qquad BA=\begin{bmatrix}2&2\\0&1\end{bmatrix}.

e2=(0,1)T\mathbf e_2=(0,1)^\mathsf TABe2=(1,1)TAB\mathbf e_2=(1,1)^\mathsf TBAe2=(2,1)TBA\mathbf e_2=(2,1)^\mathsf T。先横向缩放再剪切,与先剪切再横向缩放产生不同结果,所以矩阵乘法一般不交换。

转置会反转乘法顺序:

(AB)T=BTAT.(AB)^\mathsf T=B^\mathsf TA^\mathsf T.

按元素展开可见,两边的第 (i,j)(i,j) 个元素都等于 kbkiajk\sum_k b_{ki}a_{jk}。反转顺序也符合复合关系:交换输入输出索引时,原先的最后一步先被读取。

括号不改结果,却会改变计算量

结合律保证两种加括号方式给出相同数学结果,但中间矩阵的形状可能相差悬殊。设 AR1000×2A\in\mathbb R^{1000\times2}BR2×1000B\in\mathbb R^{2\times1000}xR1000\mathbf x\in\mathbb R^{1000}。若先计算 ABAB,要形成一个 1000×10001000\times1000 矩阵:仅这一步就需要 100010002=20000001000\cdot1000\cdot2=2\,000\,000 次标量乘法,再乘 x\mathbf x 还需 10000001\,000\,000 次。

若按 A(Bx)A(B\mathbf x) 计算,BxB\mathbf x 只产生二维向量,需要 21000=20002\cdot1000=2\,000 次标量乘法;随后乘 AA 再需要 10002=20001000\cdot2=2\,000 次,总计 40004\,000 次。两条路径由结合律保证等价,但后者既少存储一个巨大中间量,也把此处的乘法次数从约三百万降到四千。实际计算中,先检查中间形状往往比盲目从左到右更重要。

矩阵幂记录重复复合

方阵的输出与输入处于同一坐标空间,因此可以重复作用。规定 A0=IA^0=I,并令 Ak+1=AkAA^{k+1}=A^kA。结合律保证 ArAs=Ar+sA^rA^s=A^{r+s};这里的幂表示重复矩阵乘法,绝不是把每个元素分别乘方。

例如

A=[1110],A2=[2111],A3=[3221].A=\begin{bmatrix}1&1\\1&0\end{bmatrix}, \qquad A^2=\begin{bmatrix}2&1\\1&1\end{bmatrix}, \qquad A^3=\begin{bmatrix}3&2\\2&1\end{bmatrix}.

从状态 s0=(1,0)T\mathbf s_0=(1,0)^\mathsf T 出发,连续三次作用得到

s1=(1,1)T,s2=(2,1)T,s3=(3,2)T.\mathbf s_1=(1,1)^\mathsf T, \qquad \mathbf s_2=(2,1)^\mathsf T, \qquad \mathbf s_3=(3,2)^\mathsf T.

每一步的第一分量都是上一步两个分量之和,第二分量复制上一步第一分量。矩阵幂由此把一条重复更新规则压缩成 sk=Aks0\mathbf s_k=A^k\mathbf s_0

单位矩阵与逆矩阵

单位矩阵与逆矩阵

nn 阶单位矩阵 InI_n 的主对角线元素为一,其余元素为零。它满足

ImA=A,AIn=AI_mA=A,\qquad AI_n=A

对每个 ARm×nA\in\mathbb R^{m\times n} 成立。

对方阵 ARn×nA\in\mathbb R^{n\times n},若存在方阵 BB 使 AB=BA=InAB=BA=I_n,则称 AA 可逆,BBAA 的逆,记为 A1A^{-1}

只有同阶方阵才可能拥有通常意义下的双侧逆。逆的定义要求两个方向都恢复单位作用;只检验一个乘积会遗漏形状或单侧恢复问题。

逆矩阵唯一

若方阵 AA 可逆,则其逆矩阵唯一。

证明

假设 BBCC 都是 AA 的逆。利用单位矩阵和结合律,

B=BI=B(AC)=(BA)C=IC=C.B=BI=B(AC)=(BA)C=IC=C.

所以任意两个候选逆相等。

用双侧乘法核验候选逆

A=[2111],B=[1112].A=\begin{bmatrix}2&1\\1&1\end{bmatrix}, \qquad B=\begin{bmatrix}1&-1\\-1&2\end{bmatrix}.

逐项相乘得到

AB=[212+2111+2]=I2,AB= \begin{bmatrix} 2-1&-2+2\\ 1-1&-1+2 \end{bmatrix} =I_2,

同时

BA=[21112+21+2]=I2.BA= \begin{bmatrix} 2-1&1-1\\ -2+2&-1+2 \end{bmatrix} =I_2.

因此 B=A1B=A^{-1}。这个例子只核验候选,不提供一般求逆算法;后续线性方程组章节会用系统方法判断可解性并求解。

AABB 都可逆,则

(AB)1=B1A1,(AB)^{-1}=B^{-1}A^{-1},

因为 (AB)(B1A1)=AIA1=I(AB)(B^{-1}A^{-1})=AIA^{-1}=I, 反向乘积也同样为 II。复合的撤销顺序必须与执行顺序相反。

并非每个方阵都可逆,而且无需先学一般判别公式也能识别某些失败。考虑

A=[1122],v=(1,1)T.A=\begin{bmatrix}1&1\\2&2\end{bmatrix}, \qquad \mathbf v=(1,-1)^\mathsf T.

直接计算得 Av=0A\mathbf v=\mathbf0,但 v0\mathbf v\ne\mathbf0。若 A1A^{-1} 存在,在等式两侧左乘它便会得到 v=A10=0\mathbf v=A^{-1}\mathbf0=\mathbf0,矛盾。因此 AA 不可逆。这里两列相同,使不同输入 (1,0)T(1,0)^\mathsf T(0,1)T(0,1)^\mathsf T 产生同一输出;一旦输入信息发生碰撞,单凭输出便不可能唯一撤销变换。

正交矩阵保持欧氏几何

正交矩阵

实方阵 QQ 若满足

QTQ=I,Q^\mathsf TQ=I,

则称为正交矩阵。此时 Q1=QTQ^{-1}=Q^\mathsf T

QTQ=IQ^\mathsf TQ=I 表示 QQ 的各列两两正交且长度为一。正交矩阵包括平面旋转和关于过原点直线的反射;一般缩放和剪切不是正交矩阵。

正交矩阵保持内积与长度

QQ 为正交矩阵,则对任意同维向量 x,y\mathbf x,\mathbf y

Qx,Qy=x,y.\langle Q\mathbf x,Q\mathbf y\rangle =\langle\mathbf x,\mathbf y\rangle.

因此 QQ 保持长度、距离、正交性和非零向量之间的夹角。

证明

使用标准欧氏内积和转置乘法,

Qx,Qy=(Qx)T(Qy)=xTQTQy=xTy.\langle Q\mathbf x,Q\mathbf y\rangle =(Q\mathbf x)^\mathsf T(Q\mathbf y) =\mathbf x^\mathsf TQ^\mathsf TQ\mathbf y =\mathbf x^\mathsf T\mathbf y.

y=x\mathbf y=\mathbf x 得长度平方保持;距离由 QxQy2=Q(xy)2\lVert Q\mathbf x-Q\mathbf y\rVert_2 =\lVert Q(\mathbf x-\mathbf y)\rVert_2 保持;夹角公式的分子和分母也都不变。

核对九十度旋转

Q=[0110].Q=\begin{bmatrix}0&-1\\1&0\end{bmatrix}.

QTQ=[0110][0110]=I2.Q^\mathsf TQ =\begin{bmatrix}0&1\\-1&0\end{bmatrix} \begin{bmatrix}0&-1\\1&0\end{bmatrix} =I_2.

x=(3,4)T\mathbf x=(3,4)^\mathsf TQx=(4,3)TQ\mathbf x=(-4,3)^\mathsf T。两者长度都为 55。 对 y=(1,0)T\mathbf y=(1,0)^\mathsf Tx,y=3\langle\mathbf x,\mathbf y\rangle=3,而 Qy=(0,1)TQ\mathbf y=(0,1)^\mathsf TQx,Qy=3\langle Q\mathbf x,Q\mathbf y\rangle=3,与定理一致。

平面旋转提供一整族正交矩阵:

Rθ=[cosθsinθsinθcosθ].R_\theta= \begin{bmatrix} \cos\theta&-\sin\theta\\ \sin\theta&\cos\theta \end{bmatrix}.

它的两列长度都为一,内积为零,所以 RθTRθ=I2R_\theta^\mathsf TR_\theta=I_2。对角矩阵 F=diag(1,1)F=\operatorname{diag}(1,-1) 则把第二坐标反号,表示关于第一坐标轴的反射,同样满足 FTF=I2F^\mathsf TF=I_2。二者都保持长度,但一个连续转动方向,另一个翻转方向;“保持长度”并不等于“什么都没有改变”。

正交矩阵对复合封闭

若同阶方阵 Q1,Q2Q_1,Q_2 都正交,则 Q1Q2Q_1Q_2 也正交。

证明

利用乘积转置会反序以及结合律,

(Q1Q2)T(Q1Q2)=Q2TQ1TQ1Q2=Q2TIQ2=I.(Q_1Q_2)^\mathsf T(Q_1Q_2) =Q_2^\mathsf TQ_1^\mathsf TQ_1Q_2 =Q_2^\mathsf TIQ_2 =I.

因此有限次旋转与反射的复合仍保持内积、长度、距离和夹角。乘积顺序仍不能随意交换,因为不同顺序可能对应不同的几何动作。

形状正确之后还要检查语义

面对矩阵表达式,可依次做三层检查。

第一层是形状:每个乘积的内维是否相等,输出形状是否符合问题。第二层是索引语义:一行代表样本、时刻还是输出坐标,一列代表特征、状态还是输入方向。第三层是单位:若一行把不同物理量线性组合,系数必须承担相应单位,使每一项能够相加。

例如 XRs×fX\in\mathbb R^{s\times f} 表示 ss 个样本、ff 个特征,权重 wRf\mathbf w\in\mathbb R^f,则 XwRsX\mathbf w\in\mathbb R^s 为每个样本产生一个输出。若误把样本放在列上,同一任务要改用转置或重新定义存储约定。形状匹配只是必要条件,不能替代对象与单位说明。

MIT OpenCourseWare 18.06SC 从方程组和列向量进入矩阵乘法、逆矩阵、转置与正交结构,并提供逐主题习题和解答。本章引用该课程支持这些线性代数对象的教学顺序,不把应用中的同形数表自动视为语义相同。

常见误区

矩阵乘法就是逐元素相乘

逐元素乘积要求两矩阵同形,结果仍同形;矩阵乘法要求左矩阵列数等于右矩阵行数,并沿中间索引求和。两种运算的定义、形状与用途均不同。

乘积为零就至少有一个零矩阵

实数满足零乘积性质,但矩阵一般不满足。例如

A=[1111],B=[1111]A=\begin{bmatrix}1&-1\\1&-1\end{bmatrix}, \qquad B=\begin{bmatrix}1&1\\1&1\end{bmatrix}

都不是零矩阵,却有 AB=0AB=0BB 的每个输出都是 (1,1)T(1,1)^\mathsf T 的倍数,而 A(1,1)T=0A(1,1)^\mathsf T=\mathbf0,所以复合后信息全部消失。这个现象的关键不在元素大小,而在前一步的输出方向恰好落入后一步会压成零的方向。类似地,从 AC=BCAC=BC 不能贸然约去 CC;只有已知 CC 可逆时,才能在等式右侧同乘 C1C^{-1} 得到 A=BA=B。矩阵消去必须由逆矩阵提供依据。

维度对上就一定可以解释

形状相容只说明代数乘法有定义。若输入坐标顺序、单位或样本轴约定错误,结果仍可能没有建模意义。

逆矩阵就是逐元素取倒数

逆矩阵必须通过矩阵乘法恢复单位矩阵。逐元素倒数既不处理非对角耦合,也可能在零元素处无定义。

保持长度的矩阵一定是单位矩阵

旋转和关于过原点直线的反射都保持长度,却通常不是单位矩阵。正确结构是正交矩阵。

练习:矩阵运算、逆与正交性

练习

AR2×3A\in\mathbb R^{2\times3}BR2×3B\in\mathbb R^{2\times3}CR3×2C\in\mathbb R^{3\times2}。 判断 A+BA+BA+CA+CAT+CA^\mathsf T+C 是否有定义,并写出有定义结果的形状。

查看提示
先写每个矩阵的行数和列数,再判断对应位置是否齐全。
查看解答

A+BA+B 有定义,形状为 2×32\times3A+CA+C 无定义,因为形状不同。 ATA^\mathsf T 的形状为 3×23\times2,所以 AT+CA^\mathsf T+C 有定义,形状为 3×23\times2

练习

A=[1201],B=[1031],x=(1,1)T.A=\begin{bmatrix}1&2\\0&1\end{bmatrix}, \quad B=\begin{bmatrix}1&0\\3&1\end{bmatrix}, \quad \mathbf x=(1,-1)^\mathsf T.

分别计算 A(Bx)A(B\mathbf x)(AB)x(AB)\mathbf x,并比较 ABABBABA

查看提示
先算右侧矩阵对向量的作用,再与乘积矩阵交叉核对。
查看解答

Bx=(1,2)TB\mathbf x=(1,2)^\mathsf T,所以 A(Bx)=(5,2)TA(B\mathbf x)=(5,2)^\mathsf T。另一方面

AB=[7231],AB=\begin{bmatrix}7&2\\3&1\end{bmatrix},

(AB)x=(5,2)T(AB)\mathbf x=(5,2)^\mathsf T。而

BA=[1237]AB.BA=\begin{bmatrix}1&2\\3&7\end{bmatrix}\ne AB.

结合律得到核验,非交换性也被具体展示。

练习

A=[201021113],x=(1,2,1)T.A= \begin{bmatrix} 2&0&1\\ 0&2&-1\\ 1&1&3 \end{bmatrix}, \qquad \mathbf x=(1,2,-1)^\mathsf T.

按前两维与最后一维分块计算 AxA\mathbf x,再逐行复算。

查看提示
把前两个坐标视为一个向量块,分别计算上下两个输出块。
查看解答

A11=2I2A_{11}=2I_2A12=(1,1)TA_{12}=(1,-1)^\mathsf TA21=(1,1)A_{21}=(1,1)A22=(3)A_{22}=(3)x1=(1,2)T\mathbf x_1=(1,2)^\mathsf Tx2=1x_2=-1。 上块为 2x1+A12x2=(1,5)T2\mathbf x_1+A_{12}x_2=(1,5)^\mathsf T, 下块为 A21x1+3x2=1+23=0A_{21}\mathbf x_1+3x_2=1+2-3=0。 所以 Ax=(1,5,0)TA\mathbf x=(1,5,0)^\mathsf T,逐行计算得到同一结果。

练习

A=[1101],B=[2001].A=\begin{bmatrix}1&1\\0&1\end{bmatrix}, \qquad B=\begin{bmatrix}2&0\\0&1\end{bmatrix}.

写出 A1A^{-1}B1B^{-1},并求 (AB)1(AB)^{-1}。用一次乘法核验结果。

查看提示
先验证候选逆的两个乘法方向,再使用复合逆公式。
查看解答

直接核验可得

A1=[1101],B1=[1/2001].A^{-1}=\begin{bmatrix}1&-1\\0&1\end{bmatrix}, \qquad B^{-1}=\begin{bmatrix}1/2&0\\0&1\end{bmatrix}.

因此

(AB)1=B1A1=[1/21/201].(AB)^{-1}=B^{-1}A^{-1} =\begin{bmatrix}1/2&-1/2\\0&1\end{bmatrix}.

AB=[2101]AB=\begin{bmatrix}2&1\\0&1\end{bmatrix}, 两者相乘得到 I2I_2,核验成立。

练习

判断

Q=15[3443]Q=\frac15 \begin{bmatrix} 3&-4\\ 4&3 \end{bmatrix}

是否为正交矩阵,并计算 x=(5,0)T\mathbf x=(5,0)^\mathsf T 的像及其长度。

查看提示
计算转置乘原矩阵;若得到单位矩阵,再比较变换前后的长度。
查看解答

矩阵两列的长度均为一,内积为 (3(4)+43)/25=0(3\cdot(-4)+4\cdot3)/25=0,所以 QTQ=I2Q^\mathsf TQ=I_2QQ 正交。并且

Qx=(3,4)T.Q\mathbf x=(3,4)^\mathsf T.

输入和输出长度都为 55,与正交矩阵保持长度的定理一致。

与后续知识的关系

  • 线性方程组 把多个约束统一写成矩阵方程,并用初等行变换分析解。
  • 线性变换 解释矩阵为何是选定基后的映射记录。
  • 行列式 给出方阵可逆性与有向体积缩放的标量判据。
  • 特征值与特征向量 寻找方阵作用下方向保持不变的向量。

已核实资源

课程 · 2011

MIT 18.06SC Linear Algebra

Gilbert Strang

提供线性代数核心概念的完整课程顺序、例题、习题和考试材料。

打开官方来源

MIT OpenCourseWare 18.06SC 从列向量、方程组和映射复合组织矩阵语言,并系统进入逆矩阵、转置、正交性与后续分解,可用于继续练习本章每类运算。

书籍 · 2019

Interactive Linear Algebra

Dan Margalit, Joseph Rabinoff

章节含几何解释、交互图、例题、练习提示和部分解答,可用于交叉核对本册六章的代数与几何主线。

打开官方来源

Georgia Tech 的 Interactive Linear Algebra 从矩阵方程、线性变换和矩阵代数解释乘法次序、逆矩阵与正交结构,适合交叉核对本章的代数和几何叙述。

下一章

下一章把约束写成 Ax=bA\mathbf x=\mathbf b,并通过保持解集不变的行操作判断无解、唯一解与无穷多解。进入该章前,应能熟练检查矩阵与向量形状,并用列的线性组合解释 AxA\mathbf x