本节属“补充篇 · 数值线性代数工具箱” :2.x 是补充篇内部编号,不是课堂第 2 章的节号。课堂第 2 章是线性方程组的解法 。
第一章把「向量有多长」公理化为范数,但长度不是几何的全部:要谈角度与「垂直」,还需要一个能同时给出长度与角度的装置——内积 。有了内积才有正交、标准正交基,把它们按列排成方阵就得到正交矩阵 :「不改变长度与角度的坐标变换」,也是 QR 分解、SVD 与全部稳定算法的基石。本节三步:内积公理化(2.1.1);正交与标准正交基(2.1.2);刻画正交矩阵(2.1.3)。
在 R n \mathbb{R}^n R n 里,我们早就会算 x T y = ∑ i x i y i \boldsymbol{x}^{\mathrm T}\boldsymbol{y}=\sum_ix_iy_i x T y = ∑ i x i y i ,并用它定义长度与夹角。但数值计算中有时需要别的 内积:加权最小二乘要给不同观测不同的权重,有限元要在函数空间里做同一套几何,统计里要按协方差结构度量距离。这些内积都给不出 x T y \boldsymbol{x}^{\mathrm T}\boldsymbol{y} x T y 这个具体公式,却都服从同一套规则。把规则抽出来,就得到内积空间。
直觉 · 内积是「长度 + 角度」的共同来源
内积要同时管两件事:让 ⟨ x , x ⟩ \langle\boldsymbol{x},\boldsymbol{x}\rangle ⟨ x , x ⟩ 当长度的平方 (非负,且只有零向量给零);让 ⟨ x , y ⟩ \langle\boldsymbol{x},\boldsymbol{y}\rangle ⟨ x , y ⟩ 当投影 / 角度 的来源(对称 + 对加法、数乘线性)。满足这几条,「长度、距离、垂直、投影、勾股定理」就自动长出来——同一套几何,换一把尺子 。
定义 2.1.1 · 内积与内积空间
设 V V V 是实数域 R \mathbb{R} R 上的线性空间。若对任意 x , y ∈ V \boldsymbol{x},\boldsymbol{y}\in V x , y ∈ V 都指定一个实数,记作 ⟨ x , y ⟩ \langle\boldsymbol{x},\boldsymbol{y}\rangle ⟨ x , y ⟩ ,满足:
对称性 :⟨ x , y ⟩ = ⟨ y , x ⟩ \langle\boldsymbol{x},\boldsymbol{y}\rangle=\langle\boldsymbol{y},\boldsymbol{x}\rangle ⟨ x , y ⟩ = ⟨ y , x ⟩ ;
对第一变元的线性 (对第二变元由对称性同理):⟨ k x , y ⟩ = k ⟨ x , y ⟩ \langle k\boldsymbol{x},\boldsymbol{y}\rangle=k\langle\boldsymbol{x},\boldsymbol{y}\rangle ⟨ k x , y ⟩ = k ⟨ x , y ⟩ ;⟨ x + z , y ⟩ = ⟨ x , y ⟩ + ⟨ z , y ⟩ \langle\boldsymbol{x}+\boldsymbol{z},\boldsymbol{y}\rangle=\langle\boldsymbol{x},\boldsymbol{y}\rangle+\langle\boldsymbol{z},\boldsymbol{y}\rangle ⟨ x + z , y ⟩ = ⟨ x , y ⟩ + ⟨ z , y ⟩ ;
正定性 :⟨ x , x ⟩ ≥ 0 \langle\boldsymbol{x},\boldsymbol{x}\rangle\ge0 ⟨ x , x ⟩ ≥ 0 ,且 ⟨ x , x ⟩ = 0 ⟺ x = 0 \langle\boldsymbol{x},\boldsymbol{x}\rangle=0\iff\boldsymbol{x}=\boldsymbol{0} ⟨ x , x ⟩ = 0 ⟺ x = 0 ;
则称 ⟨ ⋅ , ⋅ ⟩ \langle\cdot,\cdot\rangle ⟨ ⋅ , ⋅ ⟩ 为 V V V 上的内积 (inner product),( V , ⟨ ⋅ , ⋅ ⟩ ) (V,\langle\cdot,\cdot\rangle) ( V , ⟨ ⋅ , ⋅ ⟩) 称为内积空间 (inner product space)。由内积可诱导出长度
∥ x ∥ = ⟨ x , x ⟩ , (2.1) \lVert\boldsymbol{x}\rVert=\sqrt{\langle\boldsymbol{x},\boldsymbol{x}\rangle}, \tag{2.1} ∥ x ∥ = ⟨ x , x ⟩ , ( 2.1 )
以及距离 d ( x , y ) = ∥ x − y ∥ d(\boldsymbol{x},\boldsymbol{y})=\lVert\boldsymbol{x}-\boldsymbol{y}\rVert d ( x , y ) = ∥ x − y ∥ 。
在 R n \mathbb{R}^n R n 中最常用的内积 是 ⟨ x , y ⟩ = x T y = ∑ i x i y i \langle\boldsymbol{x},\boldsymbol{y}\rangle=\boldsymbol{x}^{\mathrm T}\boldsymbol{y}=\sum_ix_iy_i ⟨ x , y ⟩ = x T y = ∑ i x i y i ,称为 Euclid 内积 ,它诱导的范数就是 ∥ ⋅ ∥ 2 \lVert\cdot\rVert_2 ∥ ⋅ ∥ 2 。但满足内积公理的远不止它一个:
例 2.1.1(新增) · 加权内积是一族合法的内积
设 D = diag ( d 1 , … , d n ) D=\operatorname{diag}(d_1,\dots,d_n) D = diag ( d 1 , … , d n ) ,其中所有 d i > 0 d_i\gt0 d i > 0 。定义
⟨ x , y ⟩ D = x T D y = ∑ i = 1 n d i x i y i . \langle\boldsymbol{x},\boldsymbol{y}\rangle_D=\boldsymbol{x}^{\mathrm T}D\boldsymbol{y}=\sum_{i=1}^n d_ix_iy_i . ⟨ x , y ⟩ D = x T D y = i = 1 ∑ n d i x i y i .
验证它是 R n \mathbb{R}^n R n 上的内积,并写出诱导范数。
解
对称性: x T D y = ∑ i d i x i y i = y T D x \boldsymbol{x}^{\mathrm T}D\boldsymbol{y}=\sum_id_ix_iy_i=\boldsymbol{y}^{\mathrm T}D\boldsymbol{x} x T D y = ∑ i d i x i y i = y T D x (D D D 是对角阵,故 D T = D D^{\mathrm T}=D D T = D ;求和号内只是纯量相乘,次序可以交换)。线性: 内积是 x \boldsymbol{x} x 的线性函数(x \boldsymbol{x} x 的每个分量都是齐次线性组合)。正定性: x T D x = ∑ i d i x i 2 ≥ 0 \boldsymbol{x}^{\mathrm T}D\boldsymbol{x}=\sum_id_ix_i^2\ge0 x T D x = ∑ i d i x i 2 ≥ 0 ;等于 0 0 0 当且仅当每个 d i x i 2 = 0 d_ix_i^2=0 d i x i 2 = 0 ,由 d i > 0 d_i\gt0 d i > 0 得 x i = 0 x_i=0 x i = 0 (i = 1 , … , n i=1,\dots,n i = 1 , … , n )。
诱导范数 ∥ x ∥ D = ∑ i d i x i 2 \lVert\boldsymbol{x}\rVert_D=\sqrt{\sum_id_ix_i^2} ∥ x ∥ D = ∑ i d i x i 2 ——给每个坐标轴不同的「刻度」 :一旦某个 d i > 1 d_i\gt1 d i > 1 ,对应的坐标方向就更「贵」。更一般地,只要 D D D 对称正定,x T D y \boldsymbol{x}^{\mathrm T}D\boldsymbol{y} x T D y 就是内积(证明:令 D = L D 0 L T D=LD_0L^{\mathrm T} D = L D 0 L T ,则 x T D y = ( L T x ) T D 0 ( L T y ) \boldsymbol{x}^{\mathrm T}D\boldsymbol{y}=(L^{\mathrm T}\boldsymbol{x})^{\mathrm T}D_0(L^{\mathrm T}\boldsymbol{y}) x T D y = ( L T x ) T D 0 ( L T y ) ,化归对角情形)。
用途: 2.2.5 的加权最小二乘、加权内积下的 Rayleigh 商(2.5.4 )、第五章 的预条件都建立在它上面。
∎
内积诱导的范数自动 满足第一章的三条范数公理——前提是 Cauchy–Schwarz 成立,而它在这里是可以证明 的结论(第一章只是把它作为引理引入):
定理 2.1.1(新增) · 内积诱导范数的三个基本性质
设 ⟨ ⋅ , ⋅ ⟩ \langle\cdot,\cdot\rangle ⟨ ⋅ , ⋅ ⟩ 是 V V V 上的内积,∥ x ∥ = ⟨ x , x ⟩ \lVert\boldsymbol{x}\rVert=\sqrt{\langle\boldsymbol{x},\boldsymbol{x}\rangle} ∥ x ∥ = ⟨ x , x ⟩ ,则对任意 x , y ∈ V \boldsymbol{x},\boldsymbol{y}\in V x , y ∈ V :
(i) Cauchy–Schwarz: ∣ ⟨ x , y ⟩ ∣ ≤ ∥ x ∥ ∥ y ∥ , (ii) 三角不等式: ∥ x + y ∥ ≤ ∥ x ∥ + ∥ y ∥ , \text{(i) Cauchy–Schwarz:}\quad\lvert\langle\boldsymbol{x},\boldsymbol{y}\rangle\rvert\le\lVert\boldsymbol{x}\rVert\,\lVert\boldsymbol{y}\rVert,\qquad \text{(ii) 三角不等式:}\quad\lVert\boldsymbol{x}+\boldsymbol{y}\rVert\le\lVert\boldsymbol{x}\rVert+\lVert\boldsymbol{y}\rVert, (i) Cauchy–Schwarz : ∣⟨ x , y ⟩∣ ≤ ∥ x ∥ ∥ y ∥ , (ii) 三角不等式: ∥ x + y ∥ ≤ ∥ x ∥ + ∥ y ∥ ,
且 ∥ ⋅ ∥ \lVert\cdot\rVert ∥ ⋅ ∥ 的正定性、齐次性成立——即内积诱导的范数确实是范数 。
证明
(i) y = 0 \boldsymbol{y}=\boldsymbol{0} y = 0 时两端为零。设 y ≠ 0 \boldsymbol{y}\neq\boldsymbol{0} y = 0 ,对任意实数 t t t ,
0 ≤ ⟨ x − t y , x − t y ⟩ = ∥ x ∥ 2 − 2 t ⟨ x , y ⟩ + t 2 ∥ y ∥ 2 . 0\le\langle\boldsymbol{x}-t\boldsymbol{y},\boldsymbol{x}-t\boldsymbol{y}\rangle=\lVert\boldsymbol{x}\rVert^2-2t\langle\boldsymbol{x},\boldsymbol{y}\rangle+t^2\lVert\boldsymbol{y}\rVert^2 . 0 ≤ ⟨ x − t y , x − t y ⟩ = ∥ x ∥ 2 − 2 t ⟨ x , y ⟩ + t 2 ∥ y ∥ 2 .
这是关于 t t t 的首项系数 ∥ y ∥ 2 > 0 \lVert\boldsymbol{y}\rVert^2\gt0 ∥ y ∥ 2 > 0 的二次式,恒非负,故判别式 ≤ 0 \le0 ≤ 0 :
4 ⟨ x , y ⟩ 2 − 4 ∥ x ∥ 2 ∥ y ∥ 2 ≤ 0 ⟹ ∣ ⟨ x , y ⟩ ∣ ≤ ∥ x ∥ ∥ y ∥ . 4\langle\boldsymbol{x},\boldsymbol{y}\rangle^2-4\lVert\boldsymbol{x}\rVert^2\lVert\boldsymbol{y}\rVert^2\le0\ \Longrightarrow\ \lvert\langle\boldsymbol{x},\boldsymbol{y}\rangle\rvert\le\lVert\boldsymbol{x}\rVert\lVert\boldsymbol{y}\rVert . 4 ⟨ x , y ⟩ 2 − 4 ∥ x ∥ 2 ∥ y ∥ 2 ≤ 0 ⟹ ∣⟨ x , y ⟩∣ ≤ ∥ x ∥ ∥ y ∥ .
(也可以直接取极小点 t = ⟨ x , y ⟩ / ∥ y ∥ 2 t=\langle\boldsymbol{x},\boldsymbol{y}\rangle/\lVert\boldsymbol{y}\rVert^2 t = ⟨ x , y ⟩ / ∥ y ∥ 2 代入——与 引理 1.3.1 的证法相同,只是这里对一般 内积也成立。)
(ii) 展开后放缩:
∥ x + y ∥ 2 = ∥ x ∥ 2 + 2 ⟨ x , y ⟩ + ∥ y ∥ 2 ≤ ∥ x ∥ 2 + 2 ∥ x ∥ ∥ y ∥ + ∥ y ∥ 2 = ( ∥ x ∥ + ∥ y ∥ ) 2 . \lVert\boldsymbol{x}+\boldsymbol{y}\rVert^2=\lVert\boldsymbol{x}\rVert^2+2\langle\boldsymbol{x},\boldsymbol{y}\rangle+\lVert\boldsymbol{y}\rVert^2\le\lVert\boldsymbol{x}\rVert^2+2\lVert\boldsymbol{x}\rVert\lVert\boldsymbol{y}\rVert+\lVert\boldsymbol{y}\rVert^2=\bigl(\lVert\boldsymbol{x}\rVert+\lVert\boldsymbol{y}\rVert\bigr)^2 . ∥ x + y ∥ 2 = ∥ x ∥ 2 + 2 ⟨ x , y ⟩ + ∥ y ∥ 2 ≤ ∥ x ∥ 2 + 2 ∥ x ∥ ∥ y ∥ + ∥ y ∥ 2 = ( ∥ x ∥ + ∥ y ∥ ) 2 .
开方即得。正定性 与齐次性 直接来自内积公理:∥ k x ∥ 2 = ⟨ k x , k x ⟩ = k 2 ∥ x ∥ 2 \lVert k\boldsymbol{x}\rVert^2=\langle k\boldsymbol{x},k\boldsymbol{x}\rangle=k^2\lVert\boldsymbol{x}\rVert^2 ∥ k x ∥ 2 = ⟨ k x , k x ⟩ = k 2 ∥ x ∥ 2 (先对第一变元、再对第二变元用线性)。
∎
不是每个范数都能由内积诱导:∥ ⋅ ∥ 1 , ∥ ⋅ ∥ ∞ \lVert\cdot\rVert_1,\lVert\cdot\rVert_\infty ∥ ⋅ ∥ 1 , ∥ ⋅ ∥ ∞ 就不行。判别工具是下面两条恒等式 ——一个范数究竟来自内积还是另有出处,会在它身上留下可检验的痕迹。
直觉 · 平行四边形恒等式在说什么
内积空间里的范数满足初中就学过的平行四边形法则 :两条对角线的平方和 = = = 四条边的平方和。这条性质之所以关键,是因为内积可以用范数反算 回来(极化恒等式)——所以「范数来自内积」和「范数满足平行四边形恒等式」是同一件事。∥ ⋅ ∥ 1 \lVert\cdot\rVert_1 ∥ ⋅ ∥ 1 、∥ ⋅ ∥ ∞ \lVert\cdot\rVert_\infty ∥ ⋅ ∥ ∞ 在 x = ( 1 , 0 ) , y = ( 0 , 1 ) \boldsymbol{x}=(1,0),\boldsymbol{y}=(0,1) x = ( 1 , 0 ) , y = ( 0 , 1 ) 上立刻露馅:它们的「单位球」有棱角,不是椭圆。
定理 2.1.2(新增) · 平行四边形恒等式与极化恒等式
在内积空间中,对任意 x , y \boldsymbol{x},\boldsymbol{y} x , y :
平行四边形恒等式: ∥ x + y ∥ 2 + ∥ x − y ∥ 2 = 2 ∥ x ∥ 2 + 2 ∥ y ∥ 2 ; (2.2) \text{平行四边形恒等式:}\quad\lVert\boldsymbol{x}+\boldsymbol{y}\rVert^2+\lVert\boldsymbol{x}-\boldsymbol{y}\rVert^2=2\lVert\boldsymbol{x}\rVert^2+2\lVert\boldsymbol{y}\rVert^2 ; \tag{2.2} 平行四边形恒等式: ∥ x + y ∥ 2 + ∥ x − y ∥ 2 = 2 ∥ x ∥ 2 + 2 ∥ y ∥ 2 ; ( 2.2 )
极化恒等式: ⟨ x , y ⟩ = 1 4 ( ∥ x + y ∥ 2 − ∥ x − y ∥ 2 ) . (2.3) \text{极化恒等式:}\quad\langle\boldsymbol{x},\boldsymbol{y}\rangle=\frac{1}{4}\Bigl(\lVert\boldsymbol{x}+\boldsymbol{y}\rVert^2-\lVert\boldsymbol{x}-\boldsymbol{y}\rVert^2\Bigr). \tag{2.3} 极化恒等式: ⟨ x , y ⟩ = 4 1 ( ∥ x + y ∥ 2 − ∥ x − y ∥ 2 ) . ( 2.3 )
反之,若 ∥ ⋅ ∥ \lVert\cdot\rVert ∥ ⋅ ∥ 是 R n \mathbb{R}^n R n 上的范数且满足(2.2),则由(2.3)定义的函数是一个内积,而它诱导出的范数恰好就是原来的 ∥ ⋅ ∥ \lVert\cdot\rVert ∥ ⋅ ∥ 。
证明
(2.2) 把两个平方都展开:
∥ x + y ∥ 2 = ∥ x ∥ 2 + 2 ⟨ x , y ⟩ + ∥ y ∥ 2 , ∥ x − y ∥ 2 = ∥ x ∥ 2 − 2 ⟨ x , y ⟩ + ∥ y ∥ 2 . \lVert\boldsymbol{x}+\boldsymbol{y}\rVert^2=\lVert\boldsymbol{x}\rVert^2+2\langle\boldsymbol{x},\boldsymbol{y}\rangle+\lVert\boldsymbol{y}\rVert^2,\qquad \lVert\boldsymbol{x}-\boldsymbol{y}\rVert^2=\lVert\boldsymbol{x}\rVert^2-2\langle\boldsymbol{x},\boldsymbol{y}\rangle+\lVert\boldsymbol{y}\rVert^2 . ∥ x + y ∥ 2 = ∥ x ∥ 2 + 2 ⟨ x , y ⟩ + ∥ y ∥ 2 , ∥ x − y ∥ 2 = ∥ x ∥ 2 − 2 ⟨ x , y ⟩ + ∥ y ∥ 2 .
两式相加即得(2.2);相减即得(2.3)。
逆方向 (Jordan–von Neumann 定理的有限维版):由(2.3)定义 ⟨ ⋅ , ⋅ ⟩ \langle\cdot,\cdot\rangle ⟨ ⋅ , ⋅ ⟩ ,对称性显然(交换 x , y \boldsymbol{x},\boldsymbol{y} x , y 只对调两个平方项),正定性来自范数的正定性(y = x \boldsymbol{y}=\boldsymbol{x} y = x 代入(2.3)即得诱导范数还原:⟨ x , x ⟩ = ∥ x ∥ 2 \langle\boldsymbol{x},\boldsymbol{x}\rangle=\lVert\boldsymbol{x}\rVert^2 ⟨ x , x ⟩ = ∥ x ∥ 2 )。
线性性是真正的难关:加法要靠一条关键恒等式 ⟨ x + z , y ⟩ + ⟨ x − z , y ⟩ = 2 ⟨ x , y ⟩ \langle\boldsymbol{x}+\boldsymbol{z},\boldsymbol{y}\rangle+\langle\boldsymbol{x}-\boldsymbol{z},\boldsymbol{y}\rangle=2\langle\boldsymbol{x},\boldsymbol{y}\rangle ⟨ x + z , y ⟩ + ⟨ x − z , y ⟩ = 2 ⟨ x , y ⟩ (把(2.3)代入、四个平方和用(2.2)两两配对);齐次性则按「整数归纳 → 有理数 → 实数(用范数连续性)」三段推进。完整证明见附录 A.2.1 。
∎
例 2.1.2(新增) · 1-范数与 ∞-范数不来自任何内积
取 x = ( 1 , 0 ) T , y = ( 0 , 1 ) T \boldsymbol{x}=(1,0)^{\mathrm T},\boldsymbol{y}=(0,1)^{\mathrm T} x = ( 1 , 0 ) T , y = ( 0 , 1 ) T ,验证 ∥ ⋅ ∥ 1 \lVert\cdot\rVert_1 ∥ ⋅ ∥ 1 与 ∥ ⋅ ∥ ∞ \lVert\cdot\rVert_\infty ∥ ⋅ ∥ ∞ 都不满足平行四边形恒等式(2.2)。
解
x + y = ( 1 , 1 ) T \boldsymbol{x}+\boldsymbol{y}=(1,1)^{\mathrm T} x + y = ( 1 , 1 ) T ,x − y = ( 1 , − 1 ) T \boldsymbol{x}-\boldsymbol{y}=(1,-1)^{\mathrm T} x − y = ( 1 , − 1 ) T 。
∞-范数: ∥ x + y ∥ ∞ 2 = 1 \lVert\boldsymbol{x}+\boldsymbol{y}\rVert_\infty^2=1 ∥ x + y ∥ ∞ 2 = 1 ,∥ x − y ∥ ∞ 2 = 1 \lVert\boldsymbol{x}-\boldsymbol{y}\rVert_\infty^2=1 ∥ x − y ∥ ∞ 2 = 1 ,左端 = 2 =2 = 2 ;右端 = 2 × 1 + 2 × 1 = 4 =2\times1+2\times1=4 = 2 × 1 + 2 × 1 = 4 。2 ≠ 4 2\neq4 2 = 4 ,失败。
1-范数: ∥ x + y ∥ 1 2 = 4 \lVert\boldsymbol{x}+\boldsymbol{y}\rVert_1^2=4 ∥ x + y ∥ 1 2 = 4 ,∥ x − y ∥ 1 2 = 4 \lVert\boldsymbol{x}-\boldsymbol{y}\rVert_1^2=4 ∥ x − y ∥ 1 2 = 4 ,左端 = 8 =8 = 8 ;右端 = 2 × 1 + 2 × 1 = 4 =2\times1+2\times1=4 = 2 × 1 + 2 × 1 = 4 。8 ≠ 4 8\neq4 8 = 4 ,失败。
几何读法: 内积空间里,任意两个向量张成的「平面切片」上的单位球必是椭圆(没有棱角)。∥ ⋅ ∥ 1 \lVert\cdot\rVert_1 ∥ ⋅ ∥ 1 的单位球是菱形、∥ ⋅ ∥ ∞ \lVert\cdot\rVert_\infty ∥ ⋅ ∥ ∞ 是正方形——有棱角,就不是内积空间 。这也解释了为什么最小二乘、投影、SVD 这些工具天然属于 2-范数:只有 2-范数来自内积。
∎
课堂对照(2026 秋·庞刚)
例 2.1.2 的「单位球有棱角就不是内积空间」,在第 2 讲精编 里有对应的板书几何直观:∥ x ∥ 2 = 1 \lVert\boldsymbol{x}\rVert_2=1 ∥ x ∥ 2 = 1 是球面、∥ x ∥ 1 \lVert\boldsymbol{x}\rVert_1 ∥ x ∥ 1 (三维)是正八面体、∥ x ∥ ∞ \lVert\boldsymbol{x}\rVert_\infty ∥ x ∥ ∞ 是立方体(盒子)——有没有棱角,一画便知。
思考 · 2-范数为什么在数值分析里「特权」这么大
(2.1)与(2.2)划出一条明确的分界:只有 2-范数拥有内积结构 ,因此正交、投影、正交变换、奇异值这一整套几何工具也只对 2-范数开放。1-范数与 ∞-范数在误差界 与稀疏性 (1.3.6 、第 12 章 )上有独特价值,但「正交化」「投影」「SVD」这些词只对 2-范数有意义。这不是偏好,而是数学结构决定的。
直觉 · 正交基 = 坐标轴互相垂直,且长度都是 1
在一组「互相垂直的坐标轴」上写坐标是最舒服的:求坐标只要投影 ,也就是做一次内积;各坐标之间互不干扰;勾股定理照样成立。本节就把这些直觉逐条证明成定理。先记住一句话:在标准正交基下,「求解」退化成「做内积」 ——而内积是 O ( n ) O(n) O ( n ) 的廉价操作,所有快速算法都从这里起步。
定义 2.1.2 · 正交、标准正交组与标准正交基
设 ⟨ ⋅ , ⋅ ⟩ \langle\cdot,\cdot\rangle ⟨ ⋅ , ⋅ ⟩ 是内积空间 V V V 上的内积。
若 ⟨ x , y ⟩ = 0 \langle\boldsymbol{x},\boldsymbol{y}\rangle=0 ⟨ x , y ⟩ = 0 ,称 x \boldsymbol{x} x 与 y \boldsymbol{y} y 正交 (orthogonal),记 x ⊥ y \boldsymbol{x}\perp\boldsymbol{y} x ⊥ y ;约定零向量与任何向量正交;
若非零向量组 v 1 , … , v k \boldsymbol{v}_1,\dots,\boldsymbol{v}_k v 1 , … , v k 两两正交(i ≠ j ⇒ v i ⊥ v j i\neq j\Rightarrow\boldsymbol{v}_i\perp\boldsymbol{v}_j i = j ⇒ v i ⊥ v j ),称为正交组 (orthogonal set);
若正交组中每个向量的范数都是 1 1 1 ,称为标准正交组 (orthonormal set);
若 V V V 的标准正交组含 dim V \dim V dim V 个向量,则它是 V V V 的标准正交基 (orthonormal basis)。
定理 2.1.3(新增) · 正交组必线性无关
内积空间中任一不含零向量的正交组 { v 1 , … , v k } \{\boldsymbol{v}_1,\dots,\boldsymbol{v}_k\} { v 1 , … , v k } 线性无关。
证明
设 ∑ i = 1 k c i v i = 0 \sum_{i=1}^{k}c_i\boldsymbol{v}_i=\boldsymbol{0} ∑ i = 1 k c i v i = 0 。对任意固定的 j j j ,两边与 v j \boldsymbol{v}_j v j 作内积:
0 = ⟨ ∑ i c i v i , v j ⟩ = ∑ i c i ⟨ v i , v j ⟩ = c j ⟨ v j , v j ⟩ = c j ∥ v j ∥ 2 . 0=\Bigl\langle\sum_ic_i\boldsymbol{v}_i,\ \boldsymbol{v}_j\Bigr\rangle=\sum_ic_i\langle\boldsymbol{v}_i,\boldsymbol{v}_j\rangle=c_j\langle\boldsymbol{v}_j,\boldsymbol{v}_j\rangle=c_j\lVert\boldsymbol{v}_j\rVert^2 . 0 = ⟨ i ∑ c i v i , v j ⟩ = i ∑ c i ⟨ v i , v j ⟩ = c j ⟨ v j , v j ⟩ = c j ∥ v j ∥ 2 .
由 v j ≠ 0 \boldsymbol{v}_j\neq\boldsymbol{0} v j = 0 得 ∥ v j ∥ 2 > 0 \lVert\boldsymbol{v}_j\rVert^2\gt0 ∥ v j ∥ 2 > 0 ,故 c j = 0 c_j=0 c j = 0 。j j j 任意,故所有系数为零,线性无关。
证明里透出的「正交投影」味道: 用内积把第 j j j 个系数单独「筛」出来,其余各项一概不受影响——这正是正交基下「坐标 = 内积」的雏形。
∎
定理 2.1.4(新增) · 勾股定理与标准正交基下的坐标
设 { q 1 , … , q k } \{\boldsymbol{q}_1,\dots,\boldsymbol{q}_k\} { q 1 , … , q k } 是标准正交组。
勾股定理: 若 x = ∑ i α i q i \boldsymbol{x}=\sum_i\alpha_i\boldsymbol{q}_i x = ∑ i α i q i ,y = ∑ i β i q i \boldsymbol{y}=\sum_i\beta_i\boldsymbol{q}_i y = ∑ i β i q i 且 x ⊥ y \boldsymbol{x}\perp\boldsymbol{y} x ⊥ y ,则 ∥ x + y ∥ 2 = ∥ x ∥ 2 + ∥ y ∥ 2 \lVert\boldsymbol{x}+\boldsymbol{y}\rVert^2=\lVert\boldsymbol{x}\rVert^2+\lVert\boldsymbol{y}\rVert^2 ∥ x + y ∥ 2 = ∥ x ∥ 2 + ∥ y ∥ 2 ;
坐标公式: α i = ⟨ x , q i ⟩ \alpha_i=\langle\boldsymbol{x},\boldsymbol{q}_i\rangle α i = ⟨ x , q i ⟩ ,即 x = ∑ i = 1 k ⟨ x , q i ⟩ q i \boldsymbol{x}=\sum_{i=1}^{k}\langle\boldsymbol{x},\boldsymbol{q}_i\rangle\boldsymbol{q}_i x = ∑ i = 1 k ⟨ x , q i ⟩ q i ;
Parseval 恒等式: ∥ x ∥ 2 = ∑ i = 1 k ∣ ⟨ x , q i ⟩ ∣ 2 \lVert\boldsymbol{x}\rVert^2=\sum_{i=1}^{k}\lvert\langle\boldsymbol{x},\boldsymbol{q}_i\rangle\rvert^2 ∥ x ∥ 2 = ∑ i = 1 k ∣⟨ x , q i ⟩ ∣ 2 。
证明
(1) 展开:∥ x + y ∥ 2 = ∥ x ∥ 2 + 2 ⟨ x , y ⟩ + ∥ y ∥ 2 = ∥ x ∥ 2 + ∥ y ∥ 2 \lVert\boldsymbol{x}+\boldsymbol{y}\rVert^2=\lVert\boldsymbol{x}\rVert^2+2\langle\boldsymbol{x},\boldsymbol{y}\rangle+\lVert\boldsymbol{y}\rVert^2=\lVert\boldsymbol{x}\rVert^2+\lVert\boldsymbol{y}\rVert^2 ∥ x + y ∥ 2 = ∥ x ∥ 2 + 2 ⟨ x , y ⟩ + ∥ y ∥ 2 = ∥ x ∥ 2 + ∥ y ∥ 2 (用 ⟨ x , y ⟩ = 0 \langle\boldsymbol{x},\boldsymbol{y}\rangle=0 ⟨ x , y ⟩ = 0 )。一般地,对任意 标准正交组,若把 x \boldsymbol{x} x 分解为「在张成方向上的部分」加上「剩下的部分」,两部分必正交——见 定理 2.2.3 。
(2) 设 x = ∑ i α i q i \boldsymbol{x}=\sum_i\alpha_i\boldsymbol{q}_i x = ∑ i α i q i ,两边与 q j \boldsymbol{q}_j q j 作内积并用标准正交性 ⟨ q i , q j ⟩ = δ i j \langle\boldsymbol{q}_i,\boldsymbol{q}_j\rangle=\delta_{ij} ⟨ q i , q j ⟩ = δ ij :
⟨ x , q j ⟩ = ∑ i α i ⟨ q i , q j ⟩ = α j . \langle\boldsymbol{x},\boldsymbol{q}_j\rangle=\sum_i\alpha_i\langle\boldsymbol{q}_i,\boldsymbol{q}_j\rangle=\alpha_j . ⟨ x , q j ⟩ = i ∑ α i ⟨ q i , q j ⟩ = α j .
(3) 用(2)的展开直接算:
∥ x ∥ 2 = ⟨ ∑ i ⟨ x , q i ⟩ q i , ∑ j ⟨ x , q j ⟩ q j ⟩ = ∑ i ∣ ⟨ x , q i ⟩ ∣ 2 . \lVert\boldsymbol{x}\rVert^2=\Bigl\langle\sum_i\langle\boldsymbol{x},\boldsymbol{q}_i\rangle\boldsymbol{q}_i,\ \sum_j\langle\boldsymbol{x},\boldsymbol{q}_j\rangle\boldsymbol{q}_j\Bigr\rangle=\sum_i\lvert\langle\boldsymbol{x},\boldsymbol{q}_i\rangle\rvert^2 . ∥ x ∥ 2 = ⟨ i ∑ ⟨ x , q i ⟩ q i , j ∑ ⟨ x , q j ⟩ q j ⟩ = i ∑ ∣⟨ x , q i ⟩ ∣ 2 .
∎
(2)是最值钱的公式 :标准正交基下求坐标不需要解方程组,每个坐标就是一次内积(O ( n 3 ) O(n^3) O ( n 3 ) 的求解退化为 O ( n 2 ) O(n^2) O ( n 2 ) 的内积)。2.3 节的 QR 分解,本质就是制造 这样一组基,把求解变廉价。把标准正交基按列排成方阵会得到什么?答案是正交矩阵 ——2.1.3 的主题。
例 2.1.3(新增) · 标准正交基下的坐标与 Parseval 核对
(1)验证
q 1 = 1 2 ( 1 1 0 ) , q 2 = 1 2 ( 1 − 1 0 ) , q 3 = ( 0 0 1 ) \boldsymbol{q}_1=\frac{1}{\sqrt2}\begin{pmatrix}1\\1\\0\end{pmatrix},\qquad \boldsymbol{q}_2=\frac{1}{\sqrt2}\begin{pmatrix}1\\-1\\0\end{pmatrix},\qquad \boldsymbol{q}_3=\begin{pmatrix}0\\0\\1\end{pmatrix} q 1 = 2 1 1 1 0 , q 2 = 2 1 1 − 1 0 , q 3 = 0 0 1
是 R 3 \mathbb{R}^3 R 3 的标准正交基,并把 x = ( 6 , 2 , 1 ) T \boldsymbol{x}=(6,2,1)^{\mathrm T} x = ( 6 , 2 , 1 ) T 用它们表示、核对 Parseval 恒等式。
(2)下面这组向量看起来 也像标准正交基,请检验它是否真的成立:
r 1 = 1 2 ( 1 1 0 ) , r 2 = 1 3 ( 1 − 1 1 ) , r 3 = 1 6 ( 1 1 − 2 ) . \boldsymbol{r}_1=\frac{1}{\sqrt2}\begin{pmatrix}1\\1\\0\end{pmatrix},\qquad \boldsymbol{r}_2=\frac{1}{\sqrt3}\begin{pmatrix}1\\-1\\1\end{pmatrix},\qquad \boldsymbol{r}_3=\frac{1}{\sqrt6}\begin{pmatrix}1\\1\\-2\end{pmatrix}. r 1 = 2 1 1 1 0 , r 2 = 3 1 1 − 1 1 , r 3 = 6 1 1 1 − 2 .
解
(1)验证正交性。 逐个算两两内积(分母暂时抛开,只看分子):
对
分子 q i T q j \boldsymbol{q}_i^{\mathrm T}\boldsymbol{q}_j q i T q j
结论
q 1 , q 2 \boldsymbol{q}_1,\boldsymbol{q}_2 q 1 , q 2
( 1 ) ( 1 ) + ( 1 ) ( − 1 ) + ( 0 ) ( 0 ) = 0 (1)(1)+(1)(-1)+(0)(0)=0 ( 1 ) ( 1 ) + ( 1 ) ( − 1 ) + ( 0 ) ( 0 ) = 0
正交 ✓
q 1 , q 3 \boldsymbol{q}_1,\boldsymbol{q}_3 q 1 , q 3
( 1 ) ( 0 ) + ( 1 ) ( 0 ) + ( 0 ) ( 1 ) = 0 (1)(0)+(1)(0)+(0)(1)=0 ( 1 ) ( 0 ) + ( 1 ) ( 0 ) + ( 0 ) ( 1 ) = 0
正交 ✓
q 2 , q 3 \boldsymbol{q}_2,\boldsymbol{q}_3 q 2 , q 3
( 1 ) ( 0 ) + ( − 1 ) ( 0 ) + ( 0 ) ( 1 ) = 0 (1)(0)+(-1)(0)+(0)(1)=0 ( 1 ) ( 0 ) + ( − 1 ) ( 0 ) + ( 0 ) ( 1 ) = 0
正交 ✓
范数:∥ q 1 ∥ 2 2 = 1 + 1 2 = 1 \lVert\boldsymbol{q}_1\rVert_2^2=\frac{1+1}{2}=1 ∥ q 1 ∥ 2 2 = 2 1 + 1 = 1 ,∥ q 2 ∥ 2 2 = 1 + 1 2 = 1 \lVert\boldsymbol{q}_2\rVert_2^2=\frac{1+1}{2}=1 ∥ q 2 ∥ 2 2 = 2 1 + 1 = 1 ,∥ q 3 ∥ 2 2 = 1 \lVert\boldsymbol{q}_3\rVert_2^2=1 ∥ q 3 ∥ 2 2 = 1 。由 定理 2.1.3 ,三个两两正交的非零向量线性无关,故是 R 3 \mathbb{R}^3 R 3 的标准正交基。
求坐标(用 定理 2.1.4 的坐标公式,只做内积):
α 1 = q 1 T x = 6 + 2 2 = 4 2 , α 2 = q 2 T x = 6 − 2 2 = 2 2 , α 3 = q 3 T x = 1. \alpha_1=\boldsymbol{q}_1^{\mathrm T}\boldsymbol{x}=\frac{6+2}{\sqrt2}=4\sqrt2,\qquad \alpha_2=\boldsymbol{q}_2^{\mathrm T}\boldsymbol{x}=\frac{6-2}{\sqrt2}=2\sqrt2,\qquad \alpha_3=\boldsymbol{q}_3^{\mathrm T}\boldsymbol{x}=1 . α 1 = q 1 T x = 2 6 + 2 = 4 2 , α 2 = q 2 T x = 2 6 − 2 = 2 2 , α 3 = q 3 T x = 1.
x = 4 2 q 1 + 2 2 q 2 + 1 ⋅ q 3 . \boldsymbol{x}=4\sqrt2\,\boldsymbol{q}_1+2\sqrt2\,\boldsymbol{q}_2+1\cdot\boldsymbol{q}_3 . x = 4 2 q 1 + 2 2 q 2 + 1 ⋅ q 3 .
重建核对: 4 2 q 1 = ( 4 , 4 , 0 ) T 4\sqrt2\,\boldsymbol{q}_1=(4,4,0)^{\mathrm T} 4 2 q 1 = ( 4 , 4 , 0 ) T ,2 2 q 2 = ( 2 , − 2 , 0 ) T 2\sqrt2\,\boldsymbol{q}_2=(2,-2,0)^{\mathrm T} 2 2 q 2 = ( 2 , − 2 , 0 ) T ,1 ⋅ q 3 = ( 0 , 0 , 1 ) T 1\cdot\boldsymbol{q}_3=(0,0,1)^{\mathrm T} 1 ⋅ q 3 = ( 0 , 0 , 1 ) T ,三者相加 = ( 6 , 2 , 1 ) T = x =(6,2,1)^{\mathrm T}=\boldsymbol{x} = ( 6 , 2 , 1 ) T = x ✓。
Parseval 核对: α 1 2 + α 2 2 + α 3 2 = 32 + 8 + 1 = 41 \alpha_1^2+\alpha_2^2+\alpha_3^2=32+8+1=41 α 1 2 + α 2 2 + α 3 2 = 32 + 8 + 1 = 41 ,而 ∥ x ∥ 2 2 = 36 + 4 + 1 = 41 \lVert\boldsymbol{x}\rVert_2^2=36+4+1=41 ∥ x ∥ 2 2 = 36 + 4 + 1 = 41 ✓。Parseval 恒等式在本例中精确成立 ——这也从另一个角度印证了这组基确实是标准正交的。
(2)r \boldsymbol{r} r 组不成立。 范数都对(1 + 1 2 = 1 \frac{1+1}{2}=1 2 1 + 1 = 1 、1 + 1 + 1 3 = 1 \frac{1+1+1}{3}=1 3 1 + 1 + 1 = 1 、1 + 1 + 4 6 = 1 \frac{1+1+4}{6}=1 6 1 + 1 + 4 = 1 ),三对两两内积里 r 1 T r 2 = 0 \boldsymbol{r}_1^{\mathrm T}\boldsymbol{r}_2=0 r 1 T r 2 = 0 也确实成立,可是
r 1 T r 3 = 1 ⋅ 1 + 1 ⋅ 1 + 0 ⋅ ( − 2 ) 12 = 2 12 ≈ 0.577 ≠ 0 , r 2 T r 3 = 1 − 1 − 2 18 = − 2 18 ≈ − 0.471 ≠ 0. \boldsymbol{r}_1^{\mathrm T}\boldsymbol{r}_3=\frac{1\cdot1+1\cdot1+0\cdot(-2)}{\sqrt{12}}=\frac{2}{\sqrt{12}}\approx0.577\neq0,\qquad \boldsymbol{r}_2^{\mathrm T}\boldsymbol{r}_3=\frac{1-1-2}{\sqrt{18}}=-\frac{2}{\sqrt{18}}\approx-0.471\neq0 . r 1 T r 3 = 12 1 ⋅ 1 + 1 ⋅ 1 + 0 ⋅ ( − 2 ) = 12 2 ≈ 0.577 = 0 , r 2 T r 3 = 18 1 − 1 − 2 = − 18 2 ≈ − 0.471 = 0.
两对不正交,故 { r 1 , r 2 , r 3 } \{\boldsymbol{r}_1,\boldsymbol{r}_2,\boldsymbol{r}_3\} { r 1 , r 2 , r 3 } 不是 标准正交基。若把它误当标准正交基套用坐标公式,拼回去得到
( x T r 1 ) r 1 + ( x T r 2 ) r 2 + ( x T r 3 ) r 3 = 4 2 r 1 + 5 3 r 2 + 6 r 3 = ( 20 3 , 10 3 , − 1 3 ) ≠ x , (\boldsymbol{x}^{\mathrm T}\boldsymbol{r}_1)\boldsymbol{r}_1+(\boldsymbol{x}^{\mathrm T}\boldsymbol{r}_2)\boldsymbol{r}_2+(\boldsymbol{x}^{\mathrm T}\boldsymbol{r}_3)\boldsymbol{r}_3 =4\sqrt2\,\boldsymbol{r}_1+\frac{5}{\sqrt3}\,\boldsymbol{r}_2+\sqrt6\,\boldsymbol{r}_3 =\Bigl(\frac{20}{3},\ \frac{10}{3},\ -\frac13\Bigr)\neq\boldsymbol{x}, ( x T r 1 ) r 1 + ( x T r 2 ) r 2 + ( x T r 3 ) r 3 = 4 2 r 1 + 3 5 r 2 + 6 r 3 = ( 3 20 , 3 10 , − 3 1 ) = x ,
且「Parseval」给出 32 + 25 3 + 6 ≈ 46.3 ≠ 41 32+\frac{25}{3}+6\approx46.3\neq41 32 + 3 25 + 6 ≈ 46.3 = 41 ——两处同时失败 ,因为正交性正是这两个公式的唯一 依据。
(2)的修复: 正确的第三个方向应满足 v ⊥ r 1 , v ⊥ r 2 \boldsymbol{v}\perp\boldsymbol{r}_1,\boldsymbol{v}\perp\boldsymbol{r}_2 v ⊥ r 1 , v ⊥ r 2 。解 v T r 1 = 0 , v T r 2 = 0 \boldsymbol{v}^{\mathrm T}\boldsymbol{r}_1=0,\ \boldsymbol{v}^{\mathrm T}\boldsymbol{r}_2=0 v T r 1 = 0 , v T r 2 = 0 得 v ∝ ( 1 , − 1 , − 2 ) T \boldsymbol{v}\propto(1,-1,-2)^{\mathrm T} v ∝ ( 1 , − 1 , − 2 ) T (复核:( 1 − 1 + 0 ) / 2 = 0 (1-1+0)/\sqrt2=0 ( 1 − 1 + 0 ) / 2 = 0 ,( 1 + 1 − 2 ) / 3 = 0 (1+1-2)/\sqrt3=0 ( 1 + 1 − 2 ) / 3 = 0 ✓),单位化后
q 3 ′ = 1 6 ( 1 − 1 − 2 ) , { r 1 , r 2 , q 3 ′ } 才是标准正交基 . \boldsymbol{q}_3'=\frac{1}{\sqrt6}\begin{pmatrix}1\\-1\\-2\end{pmatrix},\qquad \{\boldsymbol{r}_1,\boldsymbol{r}_2,\boldsymbol{q}_3'\}\ \text{才是标准正交基}. q 3 ′ = 6 1 1 − 1 − 2 , { r 1 , r 2 , q 3 ′ } 才是标准正交基 .
此时坐标变为 α 1 = 4 2 , α 2 = 5 3 , α 3 ′ = 6 − 2 − 2 6 = 2 3 \alpha_1=4\sqrt2,\ \alpha_2=\frac{5}{\sqrt3},\ \alpha_3'=\frac{6-2-2}{\sqrt6}=\sqrt{\frac23} α 1 = 4 2 , α 2 = 3 5 , α 3 ′ = 6 6 − 2 − 2 = 3 2 ,Parseval 给出 32 + 25 3 + 2 3 = 41 = ∥ x ∥ 2 2 32+\frac{25}{3}+\frac{2}{3}=41=\lVert\boldsymbol{x}\rVert_2^2 32 + 3 25 + 3 2 = 41 = ∥ x ∥ 2 2 ✓。
本题的收获: 判定一组向量是否标准正交,必须逐个验证全部两两内积 ——只看范数、或只看一部分内积,都会漏掉那对不正交的向量,后果是坐标公式与 Parseval 同时静默失效 (见 2.1.4 常见错误 1)。上面「修一个方向」的做法,正是 2.3.1 的 Gram–Schmidt 正交化。
∎
直觉 · 正交矩阵 = 不改变长度和角度的变换
正交矩阵,就是把标准正交基按列排成的方阵 ——它的列就是一组互相垂直、长度归一的新坐标轴。换个角度看:把矩阵 Q Q Q 看成对空间的一次变换 (x ↦ Q x \boldsymbol{x}\mapsto Q\boldsymbol{x} x ↦ Q x )。Q Q Q 一旦是正交矩阵,这次变换就只做两件事:旋转 或镜像 (也可能两者复合)——刚体被搬来搬去,但不拉伸、不压缩、不歪斜 。因此:
长度不变(∥ Q x ∥ = ∥ x ∥ \lVert Q\boldsymbol{x}\rVert=\lVert\boldsymbol{x}\rVert ∥ Q x ∥ = ∥ x ∥ ),距离与夹角也不变;
逆变换极其便宜:Q − 1 = Q T Q^{-1}=Q^{\mathrm T} Q − 1 = Q T (转置不用真的算,只是换个存储顺序);
误差不被放大 :∥ Q x ∥ = ∥ x ∥ \lVert Q\boldsymbol{x}\rVert=\lVert\boldsymbol{x}\rVert ∥ Q x ∥ = ∥ x ∥ 意味着浮点误差经过正交变换后仍然只受机器精度支配。
前两条让正交矩阵成为最理想的坐标系 (QR 分解、SVD、特征值算法的全部中间步骤都在反复使用正交变换),第三条则是数值稳定性理论的核心(2.3.5 会用到)。
定义 2.1.3 · 正交矩阵
设 Q ∈ R n × n Q\in\mathbb{R}^{n\times n} Q ∈ R n × n 。若
Q T Q = I , (2.4) Q^{\mathrm T}Q=I, \tag{2.4} Q T Q = I , ( 2.4 )
则称 Q Q Q 为正交矩阵 (orthogonal matrix)。对长方阵 Q ∈ R m × n Q\in\mathbb{R}^{m\times n} Q ∈ R m × n (m ≥ n m\ge n m ≥ n ),若满足 Q T Q = I n Q^{\mathrm T}Q=I_n Q T Q = I n ,则称其列向量组标准正交,Q Q Q 为列正交矩阵 ——QR 分解里的 Q Q Q 正是这一类。
定理 2.1.5(新增) · 正交矩阵的等价刻画
对 Q ∈ R n × n Q\in\mathbb{R}^{n\times n} Q ∈ R n × n ,下列命题互相等价:
#
命题
一句话含义
1
Q T Q = I Q^{\mathrm T}Q=I Q T Q = I
定义
2
Q Q T = I QQ^{\mathrm T}=I Q Q T = I
左逆即右逆
3
Q − 1 = Q T Q^{-1}=Q^{\mathrm T} Q − 1 = Q T
逆就是转置
4
Q Q Q 的列向量组是 R n \mathbb{R}^n R n 的标准正交基
列正交且范数 1
5
Q Q Q 的行向量组是 R n \mathbb{R}^n R n 的标准正交基
行同样正交
6
⟨ Q x , Q y ⟩ = ⟨ x , y ⟩ \langle Q\boldsymbol{x},Q\boldsymbol{y}\rangle=\langle\boldsymbol{x},\boldsymbol{y}\rangle ⟨ Q x , Q y ⟩ = ⟨ x , y ⟩ 对一切 x , y \boldsymbol{x},\boldsymbol{y} x , y
保内积
7
∥ Q x ∥ 2 = ∥ x ∥ 2 \lVert Q\boldsymbol{x}\rVert_2=\lVert\boldsymbol{x}\rVert_2 ∥ Q x ∥ 2 = ∥ x ∥ 2 对一切 x \boldsymbol{x} x
保长度
8
∥ Q x − Q y ∥ 2 = ∥ x − y ∥ 2 \lVert Q\boldsymbol{x}-Q\boldsymbol{y}\rVert_2=\lVert\boldsymbol{x}-\boldsymbol{y}\rVert_2 ∥ Q x − Q y ∥ 2 = ∥ x − y ∥ 2
保距离
9
∣ det Q ∣ = 1 \lvert\det Q\rvert=1 ∣ det Q ∣ = 1
体积不变(方向可翻)
10
Q Q Q 的奇异值全为 1 1 1 ;∥ Q ∥ 2 = c o n d 2 ( Q ) = 1 \lVert Q\rVert_2=\mathrm{cond}_2(Q)=1 ∥ Q ∥ 2 = cond 2 ( Q ) = 1
奇异值视角
11
Q Q Q 的(复)特征值模全为 1 1 1
特征值视角
证明(挑六个关键的,其余留作[习题 B1](exercises.md))
1 ⟹ 3 ⟹ 2: 由 Q T Q = I Q^{\mathrm T}Q=I Q T Q = I 知 Q Q Q 可逆(det ( Q T Q ) = ( det Q ) 2 = 1 ≠ 0 \det(Q^{\mathrm T}Q)=(\det Q)^2=1\neq0 det ( Q T Q ) = ( det Q ) 2 = 1 = 0 ),于是 Q − 1 = Q − 1 I = Q − 1 ( Q T Q ) = ( Q − 1 Q ) Q T = Q T Q^{-1}=Q^{-1}I=Q^{-1}(Q^{\mathrm T}Q)=(Q^{-1}Q)Q^{\mathrm T}=Q^{\mathrm T} Q − 1 = Q − 1 I = Q − 1 ( Q T Q ) = ( Q − 1 Q ) Q T = Q T ,这就是 3;再由 Q Q T = Q Q − 1 = I QQ^{\mathrm T}=Q Q^{-1}=I Q Q T = Q Q − 1 = I 得 2。
1 ⟺ 4: 记 Q = [ q 1 , … , q n ] Q=[\boldsymbol{q}_1,\dots,\boldsymbol{q}_n] Q = [ q 1 , … , q n ] ,则 ( Q T Q ) i j = q i T q j = δ i j (Q^{\mathrm T}Q)_{ij}=\boldsymbol{q}_i^{\mathrm T}\boldsymbol{q}_j=\delta_{ij} ( Q T Q ) ij = q i T q j = δ ij ,正是「两两正交且范数为 1」。
3 ⟹ 6: ⟨ Q x , Q y ⟩ = ( Q x ) T ( Q y ) = x T Q T Q y = x T y = ⟨ x , y ⟩ \langle Q\boldsymbol{x},Q\boldsymbol{y}\rangle=(Q\boldsymbol{x})^{\mathrm T}(Q\boldsymbol{y})=\boldsymbol{x}^{\mathrm T}Q^{\mathrm T}Q\boldsymbol{y}=\boldsymbol{x}^{\mathrm T}\boldsymbol{y}=\langle\boldsymbol{x},\boldsymbol{y}\rangle ⟨ Q x , Q y ⟩ = ( Q x ) T ( Q y ) = x T Q T Q y = x T y = ⟨ x , y ⟩ 。
6 ⟹ 7: 取 y = x \boldsymbol{y}=\boldsymbol{x} y = x :∥ Q x ∥ 2 2 = ⟨ Q x , Q x ⟩ = ⟨ x , x ⟩ = ∥ x ∥ 2 2 \lVert Q\boldsymbol{x}\rVert_2^2=\langle Q\boldsymbol{x},Q\boldsymbol{x}\rangle=\langle\boldsymbol{x},\boldsymbol{x}\rangle=\lVert\boldsymbol{x}\rVert_2^2 ∥ Q x ∥ 2 2 = ⟨ Q x , Q x ⟩ = ⟨ x , x ⟩ = ∥ x ∥ 2 2 。
7 ⟹ 1: 由 7,x T Q T Q x = x T x \boldsymbol{x}^{\mathrm T}Q^{\mathrm T}Q\boldsymbol{x}=\boldsymbol{x}^{\mathrm T}\boldsymbol{x} x T Q T Q x = x T x 对一切 x \boldsymbol{x} x ,即 x T ( Q T Q − I ) x = 0 \boldsymbol{x}^{\mathrm T}(Q^{\mathrm T}Q-I)\boldsymbol{x}=0 x T ( Q T Q − I ) x = 0 。记 S = Q T Q − I S=Q^{\mathrm T}Q-I S = Q T Q − I ,它是对称阵。取 x \boldsymbol{x} x 为 e i \boldsymbol{e}_i e i 得 s i i = 0 s_{ii}=0 s ii = 0 ;取 x = e i + e j \boldsymbol{x}=\boldsymbol{e}_i+\boldsymbol{e}_j x = e i + e j 得 2 s i j = 0 2s_{ij}=0 2 s ij = 0 。故 S = O S=O S = O ,即 Q T Q = I Q^{\mathrm T}Q=I Q T Q = I 。
9: 1 = det I = det ( Q T Q ) = det ( Q T ) ⋅ det Q = ( det Q ) 2 1=\det I=\det(Q^{\mathrm T}Q)=\det(Q^{\mathrm T})\cdot\det Q=(\det Q)^2 1 = det I = det ( Q T Q ) = det ( Q T ) ⋅ det Q = ( det Q ) 2 ,故 ∣ det Q ∣ = 1 \lvert\det Q\rvert=1 ∣ det Q ∣ = 1 。
10: Q T Q = I Q^{\mathrm T}Q=I Q T Q = I 的特征值全是 1 1 1 ,故奇异值全为 1 1 1 ,∥ Q ∥ 2 = σ max = 1 \lVert Q\rVert_2=\sigma_{\max}=1 ∥ Q ∥ 2 = σ m a x = 1 、c o n d 2 ( Q ) = σ 1 / σ n = 1 \mathrm{cond}_2(Q)=\sigma_1/\sigma_n=1 cond 2 ( Q ) = σ 1 / σ n = 1 。(用 定理 1.3.5 与 2.6.2 的奇异值刻画。)
11: 设 Q v = λ v Q\boldsymbol{v}=\lambda\boldsymbol{v} Q v = λ v (v ≠ 0 \boldsymbol{v}\neq\boldsymbol{0} v = 0 ),取 2-范数并用 7:∥ v ∥ 2 = ∥ Q v ∥ 2 = ∣ λ ∣ ∥ v ∥ 2 \lVert\boldsymbol{v}\rVert_2=\lVert Q\boldsymbol{v}\rVert_2=\lvert\lambda\rvert\lVert\boldsymbol{v}\rVert_2 ∥ v ∥ 2 = ∥ Q v ∥ 2 = ∣ λ ∣ ∥ v ∥ 2 ,故 ∣ λ ∣ = 1 \lvert\lambda\rvert=1 ∣ λ ∣ = 1 。
∎
关于第 9 条的几何含义: ∣ det Q ∣ = 1 \lvert\det Q\rvert=1 ∣ det Q ∣ = 1 说明正交变换保持体积 (第一章的 Jacobi 行列式出现在换元积分里,而正交变换的 Jacobi 行列式恰好是 ± 1 \pm1 ± 1 )。再看符号:det Q = + 1 \det Q=+1 det Q = + 1 对应旋转 (保定向),det Q = − 1 \det Q=-1 det Q = − 1 对应镜像 / 反射 (反定向)——这正是两类正交变换的本质区别,也是 2.3.3 的 Householder 反射(det = − 1 \det=-1 det = − 1 )与 2.3.4 的 Givens 旋转(det = + 1 \det=+1 det = + 1 )的分界线。
例 2.1.4(新增) · 二维的旋转与反射
验证平面旋转矩阵与反射矩阵都是正交矩阵,并计算各自的行列式:
R θ = ( cos θ − sin θ sin θ cos θ ) , F = ( cos 2 θ sin 2 θ sin 2 θ − cos 2 θ ) . R_\theta=\begin{pmatrix}\cos\theta&-\sin\theta\\ \sin\theta&\cos\theta\end{pmatrix},\qquad F=\begin{pmatrix}\cos2\theta&\sin2\theta\\ \sin2\theta&-\cos2\theta\end{pmatrix}. R θ = ( cos θ sin θ − sin θ cos θ ) , F = ( cos 2 θ sin 2 θ sin 2 θ − cos 2 θ ) .
解
旋转: R θ T R θ = ( cos θ sin θ − sin θ cos θ ) ( cos θ − sin θ sin θ cos θ ) = ( cos 2 θ + sin 2 θ 0 0 cos 2 θ + sin 2 θ ) = I R_\theta^{\mathrm T}R_\theta=\begin{pmatrix}\cos\theta&\sin\theta\\-\sin\theta&\cos\theta\end{pmatrix}\begin{pmatrix}\cos\theta&-\sin\theta\\ \sin\theta&\cos\theta\end{pmatrix}=\begin{pmatrix}\cos^2\theta+\sin^2\theta&0\\0&\cos^2\theta+\sin^2\theta\end{pmatrix}=I R θ T R θ = ( cos θ − sin θ sin θ cos θ ) ( cos θ sin θ − sin θ cos θ ) = ( cos 2 θ + sin 2 θ 0 0 cos 2 θ + sin 2 θ ) = I 。故正交。det R θ = cos 2 θ + sin 2 θ = 1 > 0 \det R_\theta=\cos^2\theta+\sin^2\theta=1\gt0 det R θ = cos 2 θ + sin 2 θ = 1 > 0 :旋转保定向 。
反射: 记 c = cos 2 θ , s = sin 2 θ c=\cos2\theta,s=\sin2\theta c = cos 2 θ , s = sin 2 θ 。F T F = ( c s s − c ) 2 = ( c 2 + s 2 c s − s c c s − s c s 2 + c 2 ) = I F^{\mathrm T}F=\begin{pmatrix}c&s\\s&-c\end{pmatrix}^2=\begin{pmatrix}c^2+s^2&cs-sc\\cs-sc&s^2+c^2\end{pmatrix}=I F T F = ( c s s − c ) 2 = ( c 2 + s 2 cs − sc cs − sc s 2 + c 2 ) = I 。det F = − c 2 − s 2 = − 1 < 0 \det F=-c^2-s^2=-1\lt0 det F = − c 2 − s 2 = − 1 < 0 :反射反定向 。
几何验证(数值): 取 θ = 3 0 ∘ \theta=30^{\circ} θ = 3 0 ∘ ,x = ( 1 , 0 ) T \boldsymbol{x}=(1,0)^{\mathrm T} x = ( 1 , 0 ) T 。
R 3 0 ∘ x = ( cos 3 0 ∘ , sin 3 0 ∘ ) T ≈ ( 0.8660 , 0.5 ) T R_{30^{\circ}}\boldsymbol{x}=(\cos30^{\circ},\sin30^{\circ})^{\mathrm T}\approx(0.8660,0.5)^{\mathrm T} R 3 0 ∘ x = ( cos 3 0 ∘ , sin 3 0 ∘ ) T ≈ ( 0.8660 , 0.5 ) T ,长度 1 1 1 不变,方向转了 3 0 ∘ 30^{\circ} 3 0 ∘ ;
F x = ( cos 6 0 ∘ , sin 6 0 ∘ ) T ≈ ( 0.5 , 0.8660 ) T F\boldsymbol{x}=(\cos60^{\circ},\sin60^{\circ})^{\mathrm T}\approx(0.5,0.8660)^{\mathrm T} F x = ( cos 6 0 ∘ , sin 6 0 ∘ ) T ≈ ( 0.5 , 0.8660 ) T ,长度不变,方向关于 3 0 ∘ 30^{\circ} 3 0 ∘ 方向镜像(入射角 = 反射角)。
任何正交矩阵都是「若干旋转 × 若干反射」的复合 (这是正交群的几何分类);此外 Q Q Q 的复特征值成对出现在单位圆上(e ± i θ e^{\pm i\theta} e ± i θ ,对应旋转角 θ \theta θ ),实特征值只能是 ± 1 \pm1 ± 1 。
∎
下面这条定理的证明只有两行,却是整个数值线性代数的支柱之一 ——它说明,正交变换不改变矩阵的「大小」。
定理 2.1.6(新增) · 正交变换保持 2-范数与 Frobenius 范数
设 Q ∈ R m × m Q\in\mathbb{R}^{m\times m} Q ∈ R m × m 、P ∈ R n × n P\in\mathbb{R}^{n\times n} P ∈ R n × n 为正交矩阵,A ∈ R m × n A\in\mathbb{R}^{m\times n} A ∈ R m × n ,则
∥ Q A ∥ 2 = ∥ A ∥ 2 , ∥ Q A ∥ F = ∥ A ∥ F , ∥ A P ∥ F = ∥ A ∥ F . \lVert QA\rVert_2=\lVert A\rVert_2,\qquad \lVert QA\rVert_F=\lVert A\rVert_F,\qquad \lVert AP\rVert_F=\lVert A\rVert_F . ∥ Q A ∥ 2 = ∥ A ∥ 2 , ∥ Q A ∥ F = ∥ A ∥ F , ∥ A P ∥ F = ∥ A ∥ F .
证明
Frobenius:
∥ Q A ∥ F 2 = tr ( ( Q A ) T Q A ) = tr ( A T Q T Q A ) = tr ( A T A ) = ∥ A ∥ F 2 . \lVert QA\rVert_F^2=\operatorname{tr}\bigl((QA)^{\mathrm T}QA\bigr)=\operatorname{tr}(A^{\mathrm T}Q^{\mathrm T}QA)=\operatorname{tr}(A^{\mathrm T}A)=\lVert A\rVert_F^2 . ∥ Q A ∥ F 2 = tr ( ( Q A ) T Q A ) = tr ( A T Q T Q A ) = tr ( A T A ) = ∥ A ∥ F 2 .
(用了 tr ( X Y ) = tr ( Y X ) \operatorname{tr}(XY)=\operatorname{tr}(YX) tr ( X Y ) = tr ( Y X ) 与 Q T Q = I Q^{\mathrm T}Q=I Q T Q = I 。)同理 ∥ A P ∥ F 2 = tr ( P T A T A P ) \lVert AP\rVert_F^2=\operatorname{tr}(P^{\mathrm T}A^{\mathrm T}AP) ∥ A P ∥ F 2 = tr ( P T A T A P ) ——注意 tr ( P Y ) = tr ( Y P ) \operatorname{tr}(PY)=\operatorname{tr}(YP) tr ( P Y ) = tr ( Y P ) 且这里 P T = P − 1 P^{\mathrm T}=P^{-1} P T = P − 1 ,故 tr ( P T A T A P ) = tr ( A T A P P T ) = tr ( A T A ) \operatorname{tr}(P^{\mathrm T}A^{\mathrm T}AP)=\operatorname{tr}(A^{\mathrm T}AP P^{\mathrm T})=\operatorname{tr}(A^{\mathrm T}A) tr ( P T A T A P ) = tr ( A T A P P T ) = tr ( A T A ) 。
谱范数: 用 定理 1.3.5 的 ∥ A ∥ 2 = λ max ( A T A ) \lVert A\rVert_2=\sqrt{\lambda_{\max}(A^{\mathrm T}A)} ∥ A ∥ 2 = λ m a x ( A T A ) 。令 B = A T A B=A^{\mathrm T}A B = A T A 与 B ′ = ( Q A ) T ( Q A ) = A T Q T Q A = A T A = B B'=(QA)^{\mathrm T}(QA)=A^{\mathrm T}Q^{\mathrm T}QA=A^{\mathrm T}A=B B ′ = ( Q A ) T ( Q A ) = A T Q T Q A = A T A = B ——完全相同 ,故特征值相同、∥ Q A ∥ 2 = ∥ A ∥ 2 \lVert QA\rVert_2=\lVert A\rVert_2 ∥ Q A ∥ 2 = ∥ A ∥ 2 。
注意: 上文只直接验证了 ∥ ⋅ ∥ 2 \lVert\cdot\rVert_2 ∥ ⋅ ∥ 2 在左乘 正交矩阵时不变;右乘同样成立 (∥ A P ∥ 2 = ∥ A ∥ 2 \lVert AP\rVert_2=\lVert A\rVert_2 ∥ A P ∥ 2 = ∥ A ∥ 2 ,因为 ( A P ) T ( A P ) = P T A T A P (AP)^{\mathrm T}(AP)=P^{\mathrm T}A^{\mathrm T}AP ( A P ) T ( A P ) = P T A T A P 与 A T A A^{\mathrm T}A A T A 相似);而 ∥ ⋅ ∥ F \lVert\cdot\rVert_F ∥ ⋅ ∥ F 左右乘都不变。这就是「正交变换不改变矩阵的大小」的精确含义 ——它是 QR 分解与 SVD 中「用正交变换把矩阵化简、又不失真」的依据。
∎
思考 · 正交变换为什么能「不放大误差」
用浮点数执行一次正交变换,误差来自两处:变换本身的舍入(相对误差 ∼ u \sim u ∼ u ,见 1.2.5 )与放大因子 (由 定理 2.1.6 ,放大因子 = 1 =1 = 1 )。因此
∥ Q A ~ − Q A ∥ F ≈ ∥ Q A ∥ F ⋅ c u = ∥ A ∥ F ⋅ c u \lVert \widetilde{Q A}-QA\rVert_F\approx\lVert QA\rVert_F\cdot c\,u=\lVert A\rVert_F\cdot c\,u ∥ Q A − Q A ∥ F ≈ ∥ Q A ∥ F ⋅ c u = ∥ A ∥ F ⋅ c u
——误差与矩阵大小同比例,不随变换次数累积 。对比一下:若换成非正交的变换矩阵 S S S (如一般的消元算子),放大因子是 c o n d ( S ) \mathrm{cond}(S) cond ( S ) ,c o n d ≫ 1 \mathrm{cond}\gg1 cond ≫ 1 时误差就会被放大。这正是 2.3 节的正交化(Householder / Givens)与 2.6 节的 SVD 都坚持用正交变换的原因 ,也是 QR 方法解最小二乘比正态方程稳定的原因(2.3.5 )。
本节公式速查
概念
公式 / 结论
内积公理
对称、线性、正定(定义 2.1.1 )
诱导范数
∥ x ∥ = ⟨ x , x ⟩ \lVert\boldsymbol{x}\rVert=\sqrt{\langle\boldsymbol{x},\boldsymbol{x}\rangle} ∥ x ∥ = ⟨ x , x ⟩
Cauchy–Schwarz
∣ ⟨ x , y ⟩ ∣ ≤ ∥ x ∥ ∥ y ∥ \lvert\langle\boldsymbol{x},\boldsymbol{y}\rangle\rvert\le\lVert\boldsymbol{x}\rVert\lVert\boldsymbol{y}\rVert ∣⟨ x , y ⟩∣ ≤ ∥ x ∥ ∥ y ∥
平行四边形
∥ x + y ∥ 2 + ∥ x − y ∥ 2 = 2 ∥ x ∥ 2 + 2 ∥ y ∥ 2 \lVert\boldsymbol{x}+\boldsymbol{y}\rVert^2+\lVert\boldsymbol{x}-\boldsymbol{y}\rVert^2=2\lVert\boldsymbol{x}\rVert^2+2\lVert\boldsymbol{y}\rVert^2 ∥ x + y ∥ 2 + ∥ x − y ∥ 2 = 2 ∥ x ∥ 2 + 2 ∥ y ∥ 2
极化恒等式
⟨ x , y ⟩ = 1 4 ( ∥ x + y ∥ 2 − ∥ x − y ∥ 2 ) \langle\boldsymbol{x},\boldsymbol{y}\rangle=\frac14\bigl(\lVert\boldsymbol{x}+\boldsymbol{y}\rVert^2-\lVert\boldsymbol{x}-\boldsymbol{y}\rVert^2\bigr) ⟨ x , y ⟩ = 4 1 ( ∥ x + y ∥ 2 − ∥ x − y ∥ 2 )
正交组
两两内积为 0 0 0 的非零向量组 ⟹ 线性无关
正交基坐标
x = ∑ i ⟨ x , q i ⟩ q i \boldsymbol{x}=\sum_i\langle\boldsymbol{x},\boldsymbol{q}_i\rangle\boldsymbol{q}_i x = ∑ i ⟨ x , q i ⟩ q i
Parseval
∥ x ∥ 2 = ∑ i ∣ ⟨ x , q i ⟩ ∣ 2 \lVert\boldsymbol{x}\rVert^2=\sum_i\lvert\langle\boldsymbol{x},\boldsymbol{q}_i\rangle\rvert^2 ∥ x ∥ 2 = ∑ i ∣⟨ x , q i ⟩ ∣ 2
正交矩阵
Q T Q = I Q^{\mathrm T}Q=I Q T Q = I ;Q − 1 = Q T Q^{-1}=Q^{\mathrm T} Q − 1 = Q T ;∣ det Q ∣ = 1 \lvert\det Q\rvert=1 ∣ det Q ∣ = 1 ;c o n d 2 ( Q ) = 1 \mathrm{cond}_2(Q)=1 cond 2 ( Q ) = 1
正交不变性
∥ Q A ∥ F = ∥ A ∥ F \lVert QA\rVert_F=\lVert A\rVert_F ∥ Q A ∥ F = ∥ A ∥ F ,∥ Q A ∥ 2 = ∥ A ∥ 2 \lVert QA\rVert_2=\lVert A\rVert_2 ∥ Q A ∥ 2 = ∥ A ∥ 2
常见错误清单 · 自查五条
「范数都是 1 1 1 」就当作标准正交。 范数为 1 1 1 只保证「标准」,「正交」必须逐个验证内积为零 ——例 2.1.3 的三个向量范数全对,但 r 2 T r 3 = − 2 / 18 ≠ 0 \boldsymbol{r}_2^{\mathrm T}\boldsymbol{r}_3=-2/\sqrt{18}\neq0 r 2 T r 3 = − 2/ 18 = 0 ,于是坐标公式与 Parseval 同时静默失效 。判定流程:先算全部两两内积,再算范数。
把 ∥ ⋅ ∥ 1 \lVert\cdot\rVert_1 ∥ ⋅ ∥ 1 或 ∥ ⋅ ∥ ∞ \lVert\cdot\rVert_\infty ∥ ⋅ ∥ ∞ 当成内积范数。 它们不满足平行四边形恒等式(例 2.1.2 ),因此没有 正交、投影、SVD 这套工具。见到「正交」二字,默认内积是 Euclid 内积。
以为「正交矩阵的逆要算」 。Q − 1 = Q T Q^{-1}=Q^{\mathrm T} Q − 1 = Q T ,转置是 O ( n 2 ) O(n^2) O ( n 2 ) 的数据搬移,不需要 O ( n 3 ) O(n^3) O ( n 3 ) 的求逆。代码里出现 np.linalg.inv(Q),基本可以断定是错的。
混淆「列正交」与「正交矩阵」。 长方阵 Q ∈ R m × n Q\in\mathbb{R}^{m\times n} Q ∈ R m × n (m > n m\gt n m > n )满足 Q T Q = I n Q^{\mathrm T}Q=I_n Q T Q = I n ,但 Q Q T ≠ I m QQ^{\mathrm T}\neq I_m Q Q T = I m ——它是一个 m m m 阶投影矩阵(2.2.4 )。只有方阵才有 Q − 1 = Q T Q^{-1}=Q^{\mathrm T} Q − 1 = Q T 。
以为内积只有 x T y \boldsymbol{x}^{\mathrm T}\boldsymbol{y} x T y 一种。 加权内积 x T D y \boldsymbol{x}^{\mathrm T}D\boldsymbol{y} x T D y (D D D 对称正定)同样合法且常用;不同的内积给出不同的 正交性与投影——「投影」必须说明是在哪个内积下。切换内积时,x T y = 0 \boldsymbol{x}^{\mathrm T}\boldsymbol{y}=0 x T y = 0 不再是正交的判据。
课程预告 · 本节之后
有了正交,就能定义投影 (把向量「垂直地」拍到一个子空间上),进而得到最小二乘与四个基本子空间的正交补结构——2.2 ;
有了正交矩阵,就能用正交变换 把矩阵化为三角形(QR 分解),并证明它的误差不放大——2.3 ;
有了内积,就可以定义奇异值 (2.6)与 Rayleigh 商 (2.5.4),把「矩阵有多大」「多敏感」讲成几何。