精校翻译 Ch.4 随机矩阵
第 4 章精校翻译:随机矩阵

翻译状态说明

本页为第 4 章的精校翻译,已覆盖正文、Notes、Exercises 4.1-4.51、明确 Proof 卡片和原书 Figure 4.1-4.7。

第 4 章 随机矩阵

本章开始研究随机矩阵的非渐近理论。第 4.1 节快速回顾线性代数,重点不是重复 SVD 等基础事实,而是为后续随机矩阵分析准备奇异值、算子范数、低秩逼近和扰动理论。

第 4.2 节引入 $\varepsilon$-net、covering number 与 packing number。第 4.3 节把这些几何概念连接到纠错码。第 4.4 与第 4.6 节发展 $\varepsilon$-net argument,用于证明次高斯随机矩阵算子范数上界以及奇异值双侧界。随后应用到社区检测、协方差估计和聚类。

4.1 线性代数快速回顾

先回顾一些关于矩阵的基本材料。你很可能已经在线性代数课程中学过其中不少内容;熟悉的部分可以跳过,但后面的一些内容可能是新的。

4.1.1 奇异值分解

谱分解 (3.7) 是强大的工具,但它只适用于对称矩阵。奇异值分解(singular value decomposition, SVD)把这一思想推广到所有矩阵。

Theorem 4.1.1 奇异值分解

任意实 $m\times n$ 矩阵 $A$ 都可以写成

$$ A = \sum_{i=1}^r s_i u_i v_i^{\mathsf T}, \qquad r=\min(m,n). \tag{4.1} $$

这里 $s_i$ 是非负数,称为 $A$ 的奇异值;$u_i\in\mathbb R^m$ 是正交归一向量,称为 $A$ 的左奇异向量;$v_i\in\mathbb R^n$ 是正交归一向量,称为 $A$ 的右奇异向量。

查看学习笔记完整证明
Proof 从 $A^{\mathsf T}A$ 的谱分解得到 SVD

不失一般性,可假设 $m\ge n$。由于 $A^{\mathsf T}A$ 是 $n\times n$ 对称正半定矩阵,谱定理告诉我们它有实非负特征值 $s_1^2,\ldots,s_n^2$ 和正交归一特征向量 $v_1,\ldots,v_n\in\mathbb R^n$,满足

$$ A^{\mathsf T}Av_i=s_i^2v_i. $$

向量 $Av_i$ 两两正交,因为

$$ \langle Av_i,Av_j\rangle = \langle A^{\mathsf T}Av_i,v_j\rangle = s_i^2\langle v_i,v_j\rangle = s_i^2\delta_{ij}, \qquad i,j=1,\ldots,n. \tag{4.2} $$

因此存在正交归一向量 $u_1,\ldots,u_n\in\mathbb R^m$,使得

$$ Av_i=s_iu_i, \qquad i=1,\ldots,n. \tag{4.3} $$

当 $s_i\ne0$ 时,(4.3) 唯一定义 $u_i$,而 (4.2) 保证它们正交归一。当 $s_i=0$ 时,由 (4.2)(取 $j=i$)得 $Av_i=0$,所以 (4.3) 平凡成立;此时可任选 $u_i$,只要保持正交归一性。

由于 $v_1,\ldots,v_n$ 构成 $\mathbb R^n$ 的正交归一基,

$$ I_n=\sum_{i=1}^n v_iv_i^{\mathsf T}. $$

左乘 $A$ 并使用 (4.3),得到所需 SVD:

$$ A = \sum_{i=1}^n (Av_i)v_i^{\mathsf T} = \sum_{i=1}^n s_i u_i v_i^{\mathsf T}. $$

可以思考一下:这个论证在哪一步用到了 $m\ge n$?

查看学习笔记:$m\ge n$ 用在哪里
Remark 4.1.2 矩阵做了什么

从 (4.3) 可见,SVD 给出了矩阵的几何图像:$A$ 先把每个正交方向 $v_i$ 按 $s_i$ 拉伸,然后把空间“旋转”,把正交归一基 $(v_i)$ 映到 $(u_i)$。

Remark 4.1.3 矩阵形式的 SVD

为方便起见,通常对 $i>r$ 令 $s_i=0$,并把 $s_i$ 按非增顺序排列。也可把 $(u_i)$ 和 $(v_i)$ 分别扩展为 $\mathbb R^m$ 和 $\mathbb R^n$ 的正交归一基。这样 SVD (4.1) 可写成矩阵形式

$$ A=U\Sigma V^{\mathsf T}, \tag{4.4} $$

其中 $U$ 是以左奇异向量 $u_i$ 为列的 $m\times m$ 正交矩阵,$V$ 是以右奇异向量 $v_i$ 为列的 $n\times n$ 正交矩阵,$\Sigma$ 是以奇异值 $s_i$ 为对角元素的 $m\times n$ 对角矩阵。

类似地,对称 $n\times n$ 矩阵 $A$ 的谱分解也可写成

$$ A = \sum_{i=1}^n \lambda_i u_i u_i^{\mathsf T} = U\Lambda U^{\mathsf T}, $$

其中 $U$ 是以特征向量 $u_i$ 为列的 $n\times n$ 正交矩阵,$\Lambda$ 是以特征值 $\lambda_i$ 为对角元素的 $n\times n$ 对角矩阵。

Remark 4.1.4 谱分解与 SVD

谱分解与奇异值分解紧密相关。由于

$$ AA^{\mathsf T} = \sum_{i=1}^r s_i^2 u_i u_i^{\mathsf T}, \qquad A^{\mathsf T}A = \sum_{i=1}^r s_i^2 v_i v_i^{\mathsf T}, $$

左奇异向量 $u_i$ 是 $AA^{\mathsf T}$ 的特征向量,右奇异向量 $v_i$ 是 $A^{\mathsf T}A$ 的特征向量,而 $A$ 的奇异值 $s_i$ 是 $AA^{\mathsf T}$ 和 $A^{\mathsf T}A$ 的特征值平方根:

$$ s_i(A) = \sqrt{\lambda_i(AA^{\mathsf T})} = \sqrt{\lambda_i(A^{\mathsf T}A)}. \tag{4.5} $$
Example 4.1.5 正交投影

考虑 $\mathbb R^n$ 中到 $k$ 维子空间 $E$ 上的正交投影 $P$。若 $u_1,\ldots,u_k$ 是 $E$ 的正交归一基,则向量 $x$ 到 $E$ 上的投影为

$$ Px = \sum_{i=1}^k \langle u_i,x\rangle u_i. $$

因此可把 $P$ 写成谱分解

$$ P = \sum_{i=1}^k u_i u_i^{\mathsf T} = UU^{\mathsf T}, $$

其中 $U$ 是以 $u_i$ 为正交归一列的 $n\times k$ 矩阵。特别地,$P$ 是对称矩阵,特征值为 $k$ 个 $1$ 和 $n-k$ 个 $0$。

4.1.2 Min-max 定理

我们在 Proposition 3.2.2 中见过一种基于优化的特征值描述。下面是另一种非常有用的描述。

Theorem 4.1.6 特征值的 Min-max 定理

$n\times n$ 对称矩阵 $A$ 的第 $k$ 大特征值可写成

$$ \lambda_k(A) = \max_{\dim E=k} \min_{x\in S(E)} x^{\mathsf T}Ax = \min_{\dim E=n-k+1} \max_{x\in S(E)} x^{\mathsf T}Ax, \tag{4.6} $$

其中最大值和最小值分别取遍所有维数为 $k$ 和 $n-k+1$ 的子空间 $E$,而 $S(E)$ 表示 $E$ 中的 Euclidean 单位球面,即 $E$ 中所有单位向量的集合。

Proof 在大方向与小方向之间夹住 Rayleigh quotient

只证明第一个等式。为证明 $\lambda_k=\lambda_k(A)$ 的一个方向,需要找到一个 $k$ 维子空间 $E$,使得

$$ x^{\mathsf T}Ax\ge \lambda_k \qquad \text{for all }x\in S(E). $$

如何找到这样的 $E$?取谱分解

$$ A=\sum_{i=1}^n \lambda_i u_i u_i^{\mathsf T}, $$

并选择 $A$ “最大”的子空间 $E=\operatorname{span}(u_1,\ldots,u_k)$。任意 $x\in S(E)$ 可写成 $x=\sum_{i=1}^k a_i u_i$。由 $u_i$ 的正交归一性和 $\lambda_i$ 的单调性,

$$ x^{\mathsf T}Ax = \sum_{i=1}^k \lambda_i a_i^2 \ge \lambda_k\sum_{i=1}^k a_i^2 = \lambda_k. $$

反过来,要在任意 $k$ 维子空间 $E$ 中找到一个 $x\in S(E)$,使得 $x^{\mathsf T}Ax\le\lambda_k$。看 $A$ “最小”的方向所在的子空间

$$ F=\operatorname{span}(u_k,\ldots,u_n). $$

由于 $\dim(E)+\dim(F)=n+1$,交集 $E\cap F$ 非平凡,因此存在单位向量 $x\in E\cap F$。写成 $x=\sum_{i=k}^n a_i u_i$,得到

$$ x^{\mathsf T}Ax = \sum_{i=k}^n \lambda_i a_i^2 \le \lambda_k\sum_{i=k}^n a_i^2 = \lambda_k. $$

这就得到 (4.6) 中第一个等式。第二个等式可通过把第一个等式应用于 $-A$ 并按相反顺序排列特征值得到。

把 Theorem 4.1.6 应用于 $A^{\mathsf T}A$,并使用 (4.5),立刻得到下面的奇异值版本。

Corollary 4.1.7 奇异值的 Min-max 定理

设 $A$ 是 $m\times n$ 矩阵,奇异值满足 $s_1\ge s_2\ge\cdots\ge s_n\ge0$。那么

$$ s_k(A) = \max_{\dim E=k} \min_{x\in S(E)} \|Ax\|_2 = \min_{\dim E=n-k+1} \max_{x\in S(E)} \|Ax\|_2, $$

其中最大值和最小值分别取遍所有维数为 $k$ 和 $n-k+1$ 的子空间 $E$,$S(E)$ 表示 $E$ 中的 Euclidean 单位球面。

4.1.3 Frobenius 范数与算子范数

所有 $m\times n$ 矩阵构成的线性空间有几种经典范数。最简单的是 Frobenius 范数,也称 Hilbert-Schmidt 范数;它就是 $\mathbb R^{m\times n}$ 上的 $\ell^2$ 范数:

$$ \|A\|_F := \left( \sum_{i=1}^m\sum_{j=1}^n A_{ij}^2 \right)^{1/2}. $$

简言之,Frobenius 范数就是把矩阵向量化之后的 Euclidean 范数。类似地,矩阵内积就是 $\mathbb R^{m\times n}$ 上通常的点积:

$$ \langle A,B\rangle = \sum_{i=1}^m\sum_{j=1}^n A_{ij}B_{ij} = \operatorname{tr}(A^{\mathsf T}B). \tag{4.7} $$

矩阵与自身的内积给出范数平方:

$$ \|A\|_F^2 = \langle A,A\rangle = \operatorname{tr}(A^{\mathsf T}A). \tag{4.8} $$

另一个关键矩阵范数是算子范数,也称 spectral norm。它度量线性变换 $A$ 最多能把向量拉伸多少。

Definition 4.1.8 算子范数

$m\times n$ 矩阵 $A$ 的算子范数,是使得

$$ \|Ax\|_2 \le K\|x\|_2 \qquad \text{for all }x\in\mathbb R^n $$

成立的最小数 $K$。等价地,

$$ \|A\| = \max_{x\ne0}\frac{\|Ax\|_2}{\|x\|_2} = \max_{\|x\|_2\le1}\|Ax\|_2 = \max_{\|x\|_2=1}\|Ax\|_2 = \max_{\|x\|_2=\|y\|_2=1} |y^{\mathsf T}Ax|. \tag{4.9} $$

前三个等式由缩放得到。最后一个等式来自 Euclidean 范数的对偶公式 (1.6):

$$ \|Ax\|_2 = \max_{\|y\|_2=1}\langle Ax,y\rangle. $$

在 (4.9) 中保留或去掉绝对值都得到同一个数;原因是如果 $y$ 可取,$-y$ 也可取。你将在 Exercise 4.2 中验证 $\|A\|$ 确实是一个范数。

Remark 4.1.9 其他算子范数

可以把 Definition 4.1.8 中的 $\ell^2$ 范数替换为其他范数,从而得到更一般的算子范数概念。你将在 Exercises 4.18-4.22 中探索这一点。

4.1.4 矩阵范数与谱

Lemma 4.1.10 正交不变性

Frobenius 范数和 spectral norm 都是正交不变的。也就是说,对任意矩阵 $A$ 以及维度匹配的正交矩阵 $Q$ 和 $R$,有

$$ \|QAR\|_F=\|A\|_F \qquad\text{and}\qquad \|QAR\|=\|A\|. $$
Proof 迹的循环性质与单位球面不变性

第一部分由 (4.8) 和 trace 的循环性质得到:

$$ \begin{aligned} \|QAR\|_F^2 &= \operatorname{tr} (R^{\mathsf T}A^{\mathsf T}Q^{\mathsf T}QAR)\\ &= \operatorname{tr} (R^{\mathsf T}A^{\mathsf T}AR)\\ &= \operatorname{tr} (RR^{\mathsf T}A^{\mathsf T}A)\\ &= \operatorname{tr}(A^{\mathsf T}A) = \|A\|_F^2. \end{aligned} $$

第二部分中,(4.9) 表明 $\|QAR\|$ 可通过在单位向量 $x,y$ 上最大化双线性型

$$ y^{\mathsf T}QARx = (Q^{\mathsf T}y)^{\mathsf T}A(Rx) $$

得到。由于 $Q$ 和 $R$ 是正交矩阵,$Q^{\mathsf T}y$ 与 $Rx$ 仍遍历所有单位向量,因此最大值正是 $\|A\|$。

Lemma 4.1.11 由奇异值表示矩阵范数

对任意 $m\times n$ 矩阵 $A$,若奇异值满足 $s_1(A)\ge\cdots\ge s_n(A)$,则

$$ \|A\|_F = \left( \sum_{i=1}^n s_i(A)^2 \right)^{1/2}, \qquad \|A\|=s_1(A). $$
Proof 用 SVD 和正交不变性化到对角矩阵

对 Frobenius 范数部分,使用 SVD (4.4) 与正交不变性(Lemma 4.1.10):

$$ \|A\|_F = \|U\Sigma V^{\mathsf T}\|_F = \|\Sigma\|_F. $$

而 $\Sigma$ 唯一可能非零的元素正是 $s_1,\ldots,s_n$。算子范数部分直接由奇异值 min-max 定理(Corollary 4.1.7)取 $k=1$ 得到。

Remark 4.1.12 对称矩阵

如果 $A$ 是对称矩阵,特征值为 $\lambda_k(A)$,那么

$$ \|A\| = \max_k|\lambda_k(A)| = \max_{\|x\|_2=1}|x^{\mathsf T}Ax|. \tag{4.10} $$

因此在算子范数定义 (4.9) 中可以取 $x=y$。第一个等式来自 Lemma 4.1.11,因为对称矩阵的奇异值就是 $|\lambda_k(A)|$。Min-max 定理(Theorem 4.1.6)给出

$$ |\lambda_k(A)| \le \max_{\|x\|_2=1}|x^{\mathsf T}Ax|, $$

这证明了 (4.10) 中的一个方向;另一个方向由在 (4.9) 中取 $x=y$ 立即得到。

可以通过 Exercises 4.2-4.10 进一步熟悉 Frobenius 范数和算子范数。

4.1.5 低秩逼近

假设我们想用一个给定秩 $k$ 的矩阵 $B$ 逼近矩阵 $A$。应该如何构造最佳的 $B$?以算子范数度量时,最小逼近误差是多少?答案是:截断 $A$ 的 SVD,逼近误差正是 $A$ 的第 $k+1$ 个奇异值。

Theorem 4.1.13 Eckart-Young-Mirsky 定理

设 $A$ 是 $m\times n$ 矩阵,其 SVD 为

$$ A = \sum_{i=1}^n s_i u_i v_i^{\mathsf T}. $$

那么对任意 $1\le k\le n$,有

$$ \min_{\operatorname{rank}(B)=k} \|A-B\| = s_{k+1}, $$

其中需要时约定 $s_{n+1}=0$。最小值在

$$ B = \sum_{i=1}^k s_i u_i v_i^{\mathsf T} $$

处达到。

Proof 核空间下界与截断 SVD 上界

设 $B$ 是任意秩为 $k$ 的 $m\times n$ 矩阵,其核 $E=\ker(B)$ 的维数为 $n-k$。由奇异值 min-max 定理(Corollary 4.1.7)对 $k+1$ 使用,可得

$$ \begin{aligned} \|A-B\| &\ge \max_{x\in S(E)} \|(A-B)x\|_2\\ &= \max_{x\in S(E)} \|Ax\|_2\\ &\ge s_{k+1}(A). \end{aligned} $$

反方向,取

$$ B = \sum_{i=1}^k s_i u_i v_i^{\mathsf T}, $$

$$ A-B = \sum_{i=k+1}^n s_i u_i v_i^{\mathsf T}. $$

这个矩阵的最大奇异值是 $s_{k+1}$,由 Lemma 4.1.11,它也就是算子范数。

4.1.6 扰动理论

扰动理论研究矩阵扰动时特征值和特征向量如何变化。对特征值而言,min-max 定理(Theorem 4.1.6 和 Corollary 4.1.7)的一个应用给出如下结果。

Lemma 4.1.14 Weyl 不等式

同维对称矩阵 $A$ 和 $B$ 的第 $k$ 大特征值满足

$$ |\lambda_k(A)-\lambda_k(B)| \le \|A-B\|. $$

类似地,一般矩形矩阵的第 $k$ 大奇异值满足

$$ |s_k(A)-s_k(B)| \le \|A-B\|. $$

特征向量也有类似结果,但必须小心追踪扰动前后的同一个特征向量。如果特征值太接近,一个小扰动可能交换它们的位置,导致特征向量误差很大,因为对应特征向量正交、相距很远。为避免这种情形,我们假设 $A$ 的特征值彼此分离得足够开。

Theorem 4.1.15 Davis-Kahan 不等式

考虑两个对称矩阵 $A$ 和 $B$,其谱分解分别为

$$ A=\sum_{i=1}^n\lambda_i u_i u_i^{\mathsf T}, \qquad B=\sum_{i=1}^n\mu_i v_i v_i^{\mathsf T}, $$

其中特征值按非增顺序排列。假设 $A$ 的第 $k$ 大特征值与其余特征值 $\delta$-分离:

$$ \min_{i:\,i\ne k}|\lambda_k-\lambda_i| = \delta>0. $$

那么特征向量 $u_k$ 和 $v_k$ 之间的夹角(取 $0$ 到 $\pi/2$ 之间的数)满足

$$ \sin\angle(u_k,v_k) \le \frac{2\|A-B\|}{\delta}. $$ 查看学习笔记完整证明

下面从一个更强的结果推出它。这个更强结果可同时处理多个特征向量。我们会用到 spectral projection,也就是投影到某些特征向量张成空间上的正交投影。

Lemma 4.1.16 谱投影的 Davis-Kahan 不等式

考虑两个对称矩阵 $A$ 和 $B$,谱分解分别为

$$ A=\sum_{i=1}^n\lambda_i u_i u_i^{\mathsf T}, \qquad B=\sum_{j=1}^n\mu_j v_j v_j^{\mathsf T}. $$

设 $I,J$ 是 $\mathbb R$ 中两个 $\delta$-分离的子集,且 $I$ 是区间。则谱投影

$$ P=\sum_{i:\lambda_i\in I}u_i u_i^{\mathsf T}, \qquad Q=\sum_{j:\mu_j\in J}v_j v_j^{\mathsf T} $$

满足

$$ \|QP\| \le \frac{\|A-B\|}{\delta}. $$
Proof 通过 $H=B-A$ 控制两个谱子空间的相互作用

不失一般性,假设区间 $I$ 有界且闭。给 $A$ 和 $B$ 加上相同的单位矩阵倍数,可把 $I$ 居中为 $I=[-r,r]$,使得当 $\lambda_i\in I$ 时 $|\lambda_i|\le r$,而当 $\mu_j\in J$ 时 $|\mu_j|\ge r+\delta$。令 $H:=B-A$。核心是研究 $P$ 和 $Q$ 如何通过 $H$ 相互作用:

$$ \|H\| \ge \|QHP\| = \|QBP-QAP\| \ge \|QBP\|-\|QAP\|. \tag{4.11} $$

谱投影 $Q$ 与 $B$ 交换,因此

$$ \|QBP\| = \|BQP\| \ge (r+\delta)\|QP\|. \tag{4.12} $$

最后一个不等式的原因是:$Q$ 的像空间由满足 $|\mu_j|\ge r+\delta$ 的正交向量 $v_j$ 张成,而 $B$ 把每个这样的 $v_j$ 映为 $\mu_jv_j$,至少按 $r+\delta$ 倍拉伸。因此对任意 $x$,有 $\|B QP x\|_2\ge (r+\delta)\|QP x\|_2$;再对单位向量 $x$ 取上确界即得 (4.12)。

另一方面,

$$ AP=PAP = \sum_{i:\lambda_i\in I} \lambda_i u_i u_i^{\mathsf T}, $$

所以

$$ \|QAP\| = \|QPAP\| \le \|QP\|\cdot\|AP\| \le r\|QP\|, \tag{4.13} $$

因为 $\|AP\|=\max_{\lambda_i\in I}|\lambda_i|\le r$。把 (4.12) 与 (4.13) 代入 (4.11),得到 $\|H\|\ge\delta\|QP\|$,证明完成。

Proof of Theorem 4.1.15 由谱投影版本推出单个特征向量版本

令 $\varepsilon:=\|A-B\|$。可假设 $\varepsilon\le\delta/2$,否则结论平凡。由 Weyl 不等式(Lemma 4.1.14),对每个 $j$ 有 $|\lambda_j-\mu_j|\le\varepsilon$,因此

$$ \min_{j:\,j\ne k}|\lambda_k-\mu_j| \ge \min_{j:\,j\ne k}|\lambda_k-\lambda_j|-\varepsilon = \delta-\varepsilon \ge \delta/2. $$

把 Lemma 4.1.16 应用于 $\delta/2$-分离的集合 $I=\{\lambda_k\}$ 和 $J=\{\mu_j:j\ne k\}$,得到

$$ \|QP\|\le \frac{2\varepsilon}{\delta}. $$

最后,$P$ 和 $I_n-Q$ 分别是到 $u_k$ 与 $v_k$ 方向上的正交投影,直接计算得到

$$ \|QP\| = \|Qu_k\|_2 = \sin\angle(u_k,v_k). $$

定理得证。

Exercises 4.11-4.16 会继续探索 Davis-Kahan 不等式的其他版本。

4.1.7 Isometries

矩阵 $A$ 的极端奇异值 $s_1$ 和 $s_n$ 有重要几何意义。由 min-max 定理(Corollary 4.1.7),它们可表示为

$$ s_1(A) = \max_{\|x\|_2=1}\|Ax\|_2, \qquad s_n(A) = \min_{\|x\|_2=1}\|Ax\|_2. \tag{4.14} $$

把这个不等式应用到 $x-y$ 而不是 $x$,得到对所有 $x,y\in\mathbb R^n$,

$$ s_n(A)\|x-y\|_2 \le \|Ax-Ay\|_2 \le s_1(A)\|x-y\|_2. $$

所以,极端奇异值给出了线性映射 $A$ 扭曲空间的上下限。条件数

$$ \kappa(A)=\frac{s_1(A)}{s_n(A)} $$

度量最坏情形下的扭曲因子,是数值算法中的关键量。

如果矩阵精确保距离,也就是

$$ \|Ax\|_2=\|x\|_2 \qquad \text{for all }x\in\mathbb R^n, $$

则称它为 isometry。基本例子是 $m\times m$ 正交矩阵 $U$。更一般地,任取 $U$ 的 $n$ 列构成 $m\times n$ 矩阵,它也是 isometry,给出从 $\mathbb R^n$ 到 $\mathbb R^m$ 的 isometric embedding。事实上,所有 isometry 都是这种形式:对 $m\times n$ 矩阵 $A$ 且 $m\ge n$,下面性质等价:

(a) $A$ 的列正交归一,即 $A^{\mathsf T}A=I_n$;

(b) $A$ 是 isometry;

(c) $A$ 的所有奇异值都等于 $1$。

下面证明一个更强的近似版本;它在处理随机矩阵时会很有用。

Lemma 4.1.17 近似等距

设 $A$ 是 $m\times n$ 矩阵,$m\ge n$,并设 $\varepsilon\ge0$。下面性质等价:

(a) $\|A^{\mathsf T}A-I_n\|\le\varepsilon$。

(b) 对任意 $x\in\mathbb R^n$,

$$ (1-\varepsilon)\|x\|_2^2 \le \|Ax\|_2^2 \le (1+\varepsilon)\|x\|_2^2. $$

(c)

$$ 1-\varepsilon \le s_n(A)^2 \le s_1(A)^2 \le 1+\varepsilon. $$
Proof 把近似等距写成二次型和奇异值条件

先证明 (a) $\Leftrightarrow$ (b)。由缩放,可在 (b) 中假设 $\|x\|_2=1$。由 (4.10),

$$ \begin{aligned} \|A^{\mathsf T}A-I_n\| &= \max_{\|x\|_2=1} \left| x^{\mathsf T}(A^{\mathsf T}A-I_n)x \right|\\ &= \max_{\|x\|_2=1} \left| \|Ax\|_2^2-1 \right|. \end{aligned} $$

这个量由 $\varepsilon$ 控制,正等价于 (b) 对所有单位向量成立。

(b) $\Leftrightarrow$ (c) 直接来自 (4.14)。

Remark 4.1.18 更方便的奇异值版本

下面是 Lemma 4.1.17 中关键蕴含 (a) $\Rightarrow$ (c) 的一个更方便版本。对任意 $z,\delta\ge0$,有

$$ |z^2-1| \le \max(\delta,\delta^2) \quad\Longrightarrow\quad |z-1|\le\delta. $$

因此,令 $\varepsilon=\max(\delta,\delta^2)$,可得

$$ \|A^{\mathsf T}A-I_n\| \le \max(\delta,\delta^2) \quad\Longrightarrow\quad 1-\delta \le s_n(A) \le s_1(A) \le 1+\delta. \tag{4.15} $$

Exercise 4.17 会为 Lemma 4.1.17 增加另一个方便的等价性质。

4.2 Nets, covering and packing

我们将发展一个简单但强大的方法:$\varepsilon$-net argument,并展示它在随机矩阵分析中的作用。本节先回顾 $\varepsilon$-net 的概念;你可能在实分析中见过它。随后把它与 covering、packing、entropy、volume 和 coding 等基本概念联系起来。

Definition 4.2.1 $\varepsilon$-net

设 $(T,d)$ 是度量空间,$K\subset T$,$\varepsilon>0$。子集 $\mathcal N\subset K$ 称为 $K$ 的 $\varepsilon$-net,如果 $K$ 中每个点都在 $\mathcal N$ 中某点的 $\varepsilon$ 距离内,也就是

$$ \forall x\in K\ \exists x_0\in\mathcal N: d(x,x_0)\le\varepsilon. $$

等价地,如果以 $\mathcal N$ 中点为中心、半径为 $\varepsilon$ 的球覆盖了 $K$,那么 $\mathcal N$ 就是 $K$ 的 $\varepsilon$-net;见图 4.1(a)。

如果这个一般性让人感到抽象,请记住一个关键例子:令 $T=\mathbb R^n$,并取 Euclidean 距离

$$ d(x,y)=\|x-y\|_2, \qquad x,y\in\mathbb R^n. \tag{4.16} $$

在这个情形中,我们用圆球覆盖子集 $K\subset\mathbb R^n$,如图 4.1(a)。在 Corollary 0.0.3 中我们已经见过这种覆盖的一个例子,当时 $K$ 是一个 polytope。

覆盖与装填示意
图 4.1 覆盖与装填。(a) 用六个 $\varepsilon$-球覆盖多边形 $K$,说明 $\mathcal N(K,\varepsilon)\le6$。(b) $\mathcal P(K,\varepsilon)\ge6$ 表示 $K$ 中存在六个两两 $\varepsilon$-分离的点;以这些点为中心的 $\varepsilon/2$-球两两不交。
Definition 4.2.2 覆盖数

$K$ 的 $\varepsilon$-net 的最小基数称为 $K$ 的覆盖数,记为 $\mathcal N(K,d,\varepsilon)$。等价地,$\mathcal N(K,d,\varepsilon)$ 是用中心在 $K$ 中、半径为 $\varepsilon$ 的闭球覆盖 $K$ 所需的最少球数。

Remark 4.2.3 紧性

实分析中的一个重要结果说:完备度量空间 $(T,d)$ 的子集 $K$ 是 precompact 的,也就是其闭包紧,当且仅当

$$ \mathcal N(K,d,\varepsilon)<\infty \qquad \text{for every }\varepsilon>0. $$

因此,可以把覆盖数 $\mathcal N(K,d,\varepsilon)$ 看作 $K$ 的紧性的定量度量。

与 covering 密切相关的是 packing。

Definition 4.2.4 装填数

度量空间 $(T,d)$ 的子集 $\mathcal N$ 称为 $\varepsilon$-separated,如果

$$ d(x,y)>\varepsilon \qquad \text{for any distinct points }x,y\in\mathcal N. $$

给定集合 $K\subset T$ 中 $\varepsilon$-separated 子集的最大可能基数,称为 $K$ 的装填数,记为 $\mathcal P(K,d,\varepsilon)$。

Remark 4.2.5 把球装入 $K$

如果 $\mathcal N$ 是 $\varepsilon$-separated 的,则由三角不等式,以 $\mathcal N$ 中点为中心的闭 $\varepsilon/2$-球两两不交。因此,总可以像图 4.1(b) 那样,向 $K$ 中“装入”至少 $\mathcal P(K,d,\varepsilon)$ 个两两不交的 $\varepsilon/2$-球。

Lemma 4.2.6 由 separated set 构造 net

设 $\mathcal N$ 是 $K$ 的一个极大 $\varepsilon$-separated 子集。那么 $\mathcal N$ 是 $K$ 的一个 $\varepsilon$-net。

Proof 极大性排除远点

任取 $x\in K$。我们要证明存在 $x_0\in\mathcal N$,使得 $d(x,x_0)\le\varepsilon$。如果 $x\in\mathcal N$,取 $x_0=x$ 即可。若 $x\notin\mathcal N$,由极大性,$\mathcal N\cup\{x\}$ 不是 $\varepsilon$-separated 的。这正意味着存在 $x_0\in\mathcal N$,使得 $d(x,x_0)\le\varepsilon$。

Remark 4.2.7 构造 net

Lemma 4.2.6 给出了构造给定集合 $K$ 的 $\varepsilon$-net 的迭代算法。先任取一点 $x_1\in K$,再取与 $x_1$ 距离大于 $\varepsilon$ 的 $x_2\in K$,再取与 $x_1,x_2$ 距离都大于 $\varepsilon$ 的 $x_3$,依此类推。如果 $K$ 紧,这个过程会在有限步后停止,并给出 $K$ 的一个 $\varepsilon$-net。

覆盖数和装填数本质上等价。

Lemma 4.2.8 覆盖数与装填数的等价

对任意集合 $K\subset T$ 和任意 $\varepsilon>0$,有

$$ \mathcal P(K,d,2\varepsilon) \le \mathcal N(K,d,\varepsilon) \le \mathcal P(K,d,\varepsilon). $$
Proof 极大 separated set 与抽屉原理

上界来自 Lemma 4.2.6:取 $K$ 的一个极大 $\varepsilon$-separated 子集,它既是 $\varepsilon$-net,其基数也不超过 $\mathcal P(K,d,\varepsilon)$。

为证明下界,取 $K$ 中任意 $2\varepsilon$-separated 子集 $\mathcal P=\{x_i\}$,以及任意 $\varepsilon$-net $\mathcal N=\{y_j\}$。按定义,每个 $x_i$ 都落在某个以 $y_j$ 为中心的 $\varepsilon$-球中。任意一个闭 $\varepsilon$-球不可能包含两个 $2\varepsilon$-separated 点,所以每个以 $y_j$ 为中心的球至多包含一个 $x_i$。抽屉原理给出 $|\mathcal P|\le|\mathcal N|$。由于 $\mathcal P$ 和 $\mathcal N$ 任意,下界成立。

Exercises 4.23-4.26 会帮助你练习 covering 和 packing。

4.2.1 覆盖数与体积

现在研究最重要例子中的覆盖数:$T=\mathbb R^n$,度量为 Euclidean 距离

$$ d(x,y)=\|x-y\|_2, $$

如 (4.16) 所示。为简化记号,当度量已明确时,我们常省略它,写作

$$ \mathcal N(K,\varepsilon) = \mathcal N(K,d,\varepsilon). $$

如果覆盖数度量 $K$ 的大小,那么它与最经典的大小度量,也就是 $\mathbb R^n$ 中 $K$ 的体积,有什么关系?两者不可能完全等价,因为“扁平”集合体积为零但覆盖数非零。尽管如此,有一个有用的部分等价,而且常常相当尖锐。它基于 $\mathbb R^n$ 中集合的 Minkowski 和。

Definition 4.2.9 Minkowski 和

设 $A$ 和 $B$ 是 $\mathbb R^n$ 的子集。Minkowski 和 $A+B$ 定义为

$$ A+B := \{a+b:\ a\in A,\ b\in B\}. $$

图 4.2 展示了平面上两个集合的 Minkowski 和。

正方形和圆的 Minkowski 和
图 4.2 正方形与圆的 Minkowski 和是一个带圆角的正方形。
Proposition 4.2.10 覆盖数与体积

设 $K\subset\mathbb R^n$ 且 $\varepsilon>0$。那么

$$ \frac{\operatorname{Vol}(K)} {\operatorname{Vol}(\varepsilon B_2^n)} \le \mathcal N(K,\varepsilon) \le \mathcal P(K,\varepsilon) \le \frac{ \operatorname{Vol}\bigl(K+(\varepsilon/2)B_2^n\bigr) } { \operatorname{Vol}\bigl((\varepsilon/2)B_2^n\bigr) }. $$

这里 $B_2^n$ 表示 $\mathbb R^n$ 中的 Euclidean 单位球,因此 $\varepsilon B_2^n$ 是半径为 $\varepsilon$ 的 Euclidean 球。

Proof 覆盖体积下界与装填体积上界

中间不等式来自 Lemma 4.2.8,所以只需证明左右两端。

下界。 令 $N:=\mathcal N(K,\varepsilon)$。那么 $K$ 可由 $N$ 个半径为 $\varepsilon$ 的球覆盖。比较体积得到

$$ \operatorname{Vol}(K) \le N\cdot\operatorname{Vol}(\varepsilon B_2^n), $$

这证明了下界。

上界。 令 $N:=\mathcal P(K,\varepsilon)$。可以找到 $N$ 个两两不交的闭 $\varepsilon/2$-球,其中心 $x_i$ 属于 $K$(Remark 4.2.5)。这些球未必完全包含在 $K$ 中,但它们都包含在膨胀集合 $K+(\varepsilon/2)B_2^n$ 中。比较体积得到

$$ N\cdot \operatorname{Vol}\bigl((\varepsilon/2)B_2^n\bigr) \le \operatorname{Vol}\bigl(K+(\varepsilon/2)B_2^n\bigr), $$

从而得到命题中的上界。

体积界的一个重要结论是:覆盖数以及装填数通常随维数 $n$ 指数增长。

Corollary 4.2.11 欧氏球的覆盖数

对任意 $\varepsilon>0$,Euclidean 单位球 $B_2^n$ 的覆盖数满足

$$ \left(\frac{1}{\varepsilon}\right)^n \le \mathcal N(B_2^n,\varepsilon) \le \left(\frac{2}{\varepsilon}+1\right)^n. \tag{4.17} $$

同样的上界也适用于 Euclidean 单位球面 $S^{n-1}$。

Proof 体积缩放

下界直接来自 Proposition 4.2.10,因为 $\mathbb R^n$ 中体积按如下方式缩放:

$$ \operatorname{Vol}(\varepsilon B_2^n) = \varepsilon^n\operatorname{Vol}(B_2^n). $$

上界同样来自 Proposition 4.2.10:

$$ \mathcal N(B_2^n,\varepsilon) \le \frac{ \operatorname{Vol}\bigl((1+\varepsilon/2)B_2^n\bigr) } { \operatorname{Vol}\bigl((\varepsilon/2)B_2^n\bigr) } = \frac{(1+\varepsilon/2)^n}{(\varepsilon/2)^n} = \left(\frac{2}{\varepsilon}+1\right)^n. $$

球面的上界也可用同样方式证明。

为简化 (4.17),注意在非平凡范围 $\varepsilon\in(0,1]$ 中,有

$$ \left(\frac{1}{\varepsilon}\right)^n \le \mathcal N(B_2^n,\varepsilon) \le \left(\frac{3}{\varepsilon}\right)^n. \tag{4.18} $$

在平凡范围 $\varepsilon>1$ 中,一个 $\varepsilon$-球就可以覆盖单位球,所以 $\mathcal N(B_2^n,\varepsilon)=1$。

Remark 4.2.12 球的体积

Corollary 4.2.11 的证明使用了 Euclidean 球的体积,但巧妙地避免了实际计算这个体积。可以在 Exercises 4.27-4.29 中尝试用三种方式计算它:几何、概率和分析方式;并在 Exercise 4.30 中把它推广到 $\ell^p$ 球。

Remark 4.2.13 如何构造 net

Remark 4.2.7 描述了构造 $\varepsilon$-net 的一般迭代算法。但对 Euclidean 球而言,也可以直接使用缩放后的整数格点(Exercise 4.31),甚至使用随机点(Exercise 4.39)。

刚才的体积论证很灵活,可用于许多其他场景。下面是一个重要例子。

Definition 4.2.14 Hamming 距离

Hamming 立方体 $\{0,1\}^n$ 由所有长度为 $n$ 的二进制字符串构成。为了把它变成度量空间,定义 Hamming 距离 $d_H(x,y)$ 为两个字符串 $x$ 和 $y$ 不同的比特数:

$$ d_H(x,y) := \#\{i:\ x(i)\ne y(i)\}, \qquad x,y\in\{0,1\}^n. $$

可以验证这确实是一个度量。

查看学习笔记:Hamming 距离为何是度量
Proposition 4.2.15 Hamming 立方体的覆盖数和装填数

Hamming 立方体 $K=\{0,1\}^n$ 的覆盖数和装填数满足:对任意整数 $m\in[0,n]$,

$$ \frac{2^n} {\sum_{k=0}^m \binom{n}{k}} \le \mathcal N(K,d_H,m) \le \mathcal P(K,d_H,m) \le \frac{2^n} {\sum_{k=0}^{\lfloor m/2\rfloor}\binom{n}{k}}. $$

证明只是对体积论证做一个小修改:把体积替换为基数。请在 Exercise 4.32 中尝试证明。为进一步练习体积论证,可在 Exercise 4.50 中计算低秩矩阵的覆盖数。

4.3 应用:纠错码

Covering 和 packing 论证经常出现在编码理论中。这里给出两个例子,把 covering/packing numbers 与复杂度和纠错联系起来。

4.3.1 Metric entropy 与复杂度

直观上,覆盖数和装填数度量集合 $K$ 的复杂度。覆盖数的对数 $\log_2\mathcal N(K,\varepsilon)$ 常称为 $K$ 的 metric entropy。下面会看到,metric entropy 等价于用给定精度编码 $K$ 中点所需的比特数。

Proposition 4.3.1 Metric entropy 与编码

设 $(T,d)$ 是度量空间,$K\subset T$。令 $\mathcal C(K,d,\varepsilon)$ 表示用度量 $d$ 中精度 $\varepsilon$ 指定 $K$ 中每个点所需的最少比特数。那么

$$ \log_2\mathcal N(K,d,\varepsilon) \le \mathcal C(K,d,\varepsilon) \le \left\lceil \log_2\mathcal N(K,d,\varepsilon/2) \right\rceil. $$
Proof 编码划分与 net 编码

下界。 假设 $\mathcal C(K,d,\varepsilon)\le N$。这表示存在一个把点 $x\in K$ 编码为长度为 $N$ 的比特串的映射,并且能以精度 $\varepsilon$ 指定每个点。这个编码把定义域 $K$ 分成至多 $2^N$ 个子集,每个子集由映到同一个比特串的点组成;见图 4.3。每个子集的直径至多为 $\varepsilon$,因此可用一个中心在 $K$ 中、半径为 $\varepsilon$ 的球覆盖。于是 $K$ 可由至多 $2^N$ 个半径为 $\varepsilon$ 的球覆盖,即 $\mathcal N(K,d,\varepsilon)\le2^N$。取对数得到下界。

上界。 假设 $\log_2\mathcal N(K,d,\varepsilon/2)\le N$,其中 $N$ 是整数。这表示 $K$ 有一个 $(\varepsilon/2)$-net $\mathcal N$,其点数至多为 $2^N$。对每个 $x\in K$,指定一个最近的 $x_0\in\mathcal N$。由于 $\mathcal N$ 中至多有 $2^N$ 个点,$N$ 个比特足以指定 $x_0$。

如果 $x$ 和 $y$ 被编码为同一个 $x_0$,则由三角不等式,

$$ d(x,y) \le d(x,x_0)+d(y,x_0) \le \frac{\varepsilon}{2} + \frac{\varepsilon}{2} = \varepsilon. $$

因此该编码以精度 $\varepsilon$ 表示 $K$ 中点,从而 $\mathcal C(K,d,\varepsilon)\le N$。

用 N 比特编码把集合分割成至多 2^N 个子集
图 4.3 用 $N$ 比特串编码 $K$ 中的点,会把 $K$ 分割成至多 $2^N$ 个子集。

4.3.2 纠错码

假设 Alice 想给 Bob 发送一条有 $k$ 个字母的消息,例如 “fill the glass”。再假设攻击者可以篡改 Alice 的消息,最多改变 $r$ 个字母。如果 $r=2$,Bob 可能收到类似 “bill the class” 的消息。有没有办法保护 Alice 和 Bob 之间的通信信道,使它可以纠正这种对抗性错误?

一种常见方法是使用冗余。Alice 把她的 $k$ 字母消息编码成更长的 $n$ 字母消息($n>k$),希望额外信息能帮助 Bob 在最多 $r$ 个错误出现时恢复原消息。

Example 4.3.2 重复码

Alice 可以简单地把消息重复若干次发送给 Bob。例如把 “fill the glass” 重复多遍。Bob 可以使用多数解码:检查每个位置收到的多个副本,选择出现次数最多的字母。如果原始消息 $x$ 被重复 $2r+1$ 次,那么即使 $E(x)$ 中有 $r$ 个字母被篡改,多数解码也能正确恢复 $x$。

多数解码的问题是效率低。它要用

$$ n=(2r+1)k \tag{4.19} $$

个字母来编码一个 $k$ 字母消息。马上会看到,存在所需 $n$ 小得多的纠错码。

先形式化纠错码的概念:它是一种把 $k$ 字母字符串编码为 $n$ 字母字符串,并可纠正 $r$ 个错误的编码。为简单起见,我们使用二进制 alphabet,也就是只有 $0$ 和 $1$ 两个字母,而不是英文字母。

Definition 4.3.3 纠错码

一个把 $k$-bit 字符串编码为 $n$-bit 字符串、并能纠正 $r$ 个错误的纠错码,由编码映射

$$ E:\{0,1\}^k\to\{0,1\}^n $$

和解码映射

$$ D:\{0,1\}^n\to\{0,1\}^k $$

组成,并满足:对任意 word $x\in\{0,1\}^k$,以及任意与 $E(x)$ 至多相差 $r$ 个 bit 的字符串 $y\in\{0,1\}^n$,都有

$$ D(y)=x. $$

现在把纠错与 Hamming 立方体 $\{0,1\}^n$ 的装填数联系起来;这里使用 Definition 4.2.14 中引入的 Hamming 度量。

Lemma 4.3.4 纠错与 packing

假设正整数 $k,n,r$ 满足

$$ \log_2 \mathcal P(\{0,1\}^n,d_H,2r) \ge k. $$

那么存在一个纠错码,它把 $k$-bit 字符串编码为 $n$-bit 字符串,并能纠正 $r$ 个错误。

Proof 把 codewords 选成 Hamming cube 中的 packing

由假设,存在子集 $\mathcal N\subset\{0,1\}^n$,满足 $|\mathcal N|=2^k$,并且以 $\mathcal N$ 中点为中心、半径为 $r$ 的闭球两两不交。令编码器 $E:\{0,1\}^k\to\mathcal N$ 是任意一一映射,并令 $D:\{0,1\}^n\to\{0,1\}^k$ 是最近邻解码器。

如果 $y\in\{0,1\}^n$ 与 $E(x)$ 至多相差 $r$ 个 bit,那么 $y$ 落在以 $E(x)$ 为中心、半径为 $r$ 的闭球中。由于这些球两两不交,$y$ 严格更接近 $E(x)$,而不是任何其他 codeword $E(x')$。因此最近邻解码会正确解码 $y$,即 $D(y)=x$。

把 Proposition 4.2.15 中 Hamming 立方体的 packing number 界代入 Lemma 4.3.4,得到下面的保证。

Theorem 4.3.5 纠错码保证

假设正整数 $k,n,r$ 满足

$$ n-k \ge 2r\log_2\left(\frac{en}{2r}\right). $$

那么存在一个纠错码,把 $k$-bit 字符串编码为 $n$-bit 字符串,并能纠正 $r$ 个错误。

Proof Hamming cube 的体积界给出足够多 codewords

先用 Lemma 4.2.8 从 packing 转到 covering,再使用 Proposition 4.2.15 中的覆盖数界,并用 Exercise 0.6 的估计化简,得到

$$ \begin{aligned} \mathcal P(\{0,1\}^n,d_H,2r) &\ge \mathcal N(\{0,1\}^n,d_H,2r)\\ &\ge \frac{2^n}{\sum_{i=0}^{2r}\binom{n}{i}}\\ &\ge 2^n\left(\frac{2r}{en}\right)^{2r}. \end{aligned} $$

由定理假设,最后这个量至少为 $2^k$。应用 Lemma 4.3.4 即完成证明。

Remark 4.3.6 额外比特数几乎随错误数线性增长

Theorem 4.3.5 表明,为纠正 $r$ 个错误,只需要 $n-k$ 近似按 $r$ 线性增长,忽略一个对数因子。这比重复码 (4.19) 高效得多,并且是最优的;见 Exercise 4.33。

4.4 次高斯随机矩阵的上界

现在可以进入随机矩阵的非渐近理论。随机矩阵理论研究具有随机元素的 $m\times n$ 矩阵 $A$。核心问题围绕奇异值、特征值(如果 $A$ 对称)和特征向量的分布。

Theorem 4.4.3 将给出第一个关于独立次高斯元素随机矩阵的算子范数(即最大奇异值)的界。它既不是最尖锐的,也不是最一般的结果;我们会在 Sections 4.6 和 6.4 中改进和推广它。

在此之前,先看看 $\varepsilon$-net 如何帮助计算矩阵的算子范数。

4.4.1 在 $\varepsilon$-net 上计算范数

Definition 4.1.8 中引入的 $m\times n$ 矩阵 $A$ 的算子范数为

$$ \|A\| = \max_{x\in S^{n-1}}\|Ax\|_2. $$

为了估计 $\|A\|$,需要在整个球面 $S^{n-1}$ 上一致控制 $\|Ax\|_2$,这有时很困难。下面说明,其实只需控制球面的一个 $\varepsilon$-net(Euclidean 度量下)即可。

Lemma 4.4.1 在 net 上计算算子范数

设 $A$ 是 $m\times n$ 矩阵,$\varepsilon\in[0,1)$。那么对球面 $S^{n-1}$ 的任意 $\varepsilon$-net $\mathcal N$,有

$$ \sup_{x\in\mathcal N}\|Ax\|_2 \le \|A\| \le \frac{1}{1-\varepsilon} \sup_{x\in\mathcal N}\|Ax\|_2. $$
Proof 用 net 点逼近最大方向

下界显然成立,因为 $\mathcal N\subset S^{n-1}$。为证明上界,固定 $x\in S^{n-1}$ 使得 $\|A\|=\|Ax\|_2$,并选择 $x_0\in\mathcal N$ 逼近 $x$,使得 $\|x-x_0\|_2\le\varepsilon$。由算子范数定义,

$$ \|Ax-Ax_0\|_2 = \|A(x-x_0)\|_2 \le \|A\|\|x-x_0\|_2 \le \varepsilon\|A\|. $$

由三角不等式,

$$ \|Ax_0\|_2 \ge \|Ax\|_2-\|Ax-Ax_0\|_2 \ge (1-\varepsilon)\|A\|. $$

除以 $1-\varepsilon$ 即得上界。

Lemma 4.4.1 很灵活。下面是一个有用版本。由 (4.9),$m\times n$ 矩阵 $A$ 的算子范数也可通过最大化双线性型得到:

$$ \|A\| = \max_{x\in S^{n-1},\,y\in S^{m-1}} |\langle Ax,y\rangle|. $$

如果 $A$ 对称,还可取 $x=y$;见 (4.10)。我们可以把两个球面都替换成各自的 net。

Lemma 4.4.2 在 net 上最大化二次型

设 $A$ 是 $m\times n$ 矩阵,$\varepsilon\in[0,1/2)$。那么对球面 $S^{n-1}$ 的任意 $\varepsilon$-net $\mathcal N$ 和球面 $S^{m-1}$ 的任意 $\varepsilon$-net $\mathcal M$,有

$$ \sup_{x\in\mathcal N,\,y\in\mathcal M} |\langle Ax,y\rangle| \le \|A\| \le \frac{1}{1-2\varepsilon} \sup_{x\in\mathcal N,\,y\in\mathcal M} |\langle Ax,y\rangle|. $$

此外,如果 $m=n$、$A$ 对称且 $\mathcal N=\mathcal M$,则可取 $x=y$。

这个结论可通过微调 Lemma 4.4.1 的证明得到;见 Exercise 4.36。另一种方法见 Exercise 4.34,更多练习见 Exercise 4.37。

4.4.2 次高斯随机矩阵的范数

现在给出第一个随机矩阵结果。它说明,具有独立次高斯元素的 $m\times n$ 随机矩阵 $A$ 的算子范数满足

$$ \|A\|\lesssim \sqrt m+\sqrt n \qquad \text{with high probability}. $$

Theorem 4.4.3 次高斯矩阵的范数

设 $A$ 是 $m\times n$ 随机矩阵,其元素 $A_{ij}$ 独立、均值为零且次高斯。那么对任意 $t>0$,有

$$ \|A\| \le CK(\sqrt m+\sqrt n+t) $$

以至少 $1-2\exp(-t^2)$ 的概率成立。这里

$$ K=\max_{i,j}\|A_{ij}\|_{\psi_2}. $$ 查看学习笔记完整证明
Proof $\varepsilon$-net argument

这个证明是 $\varepsilon$-net argument 的典型例子。我们需要对单位球面上所有 $x,y$ 控制 $\langle Ax,y\rangle$。为此,先用 net 离散化球面(逼近步骤),再对来自 net 的固定 $x,y$ 建立集中界(集中步骤),最后对 net 中所有点做 union bound。

Step 1:逼近。 取 $\varepsilon=1/4$。由 Corollary 4.2.11,可取球面 $S^{n-1}$ 的 $\varepsilon$-net $\mathcal N$ 和球面 $S^{m-1}$ 的 $\varepsilon$-net $\mathcal M$,满足

$$ |\mathcal N|\le 9^n, \qquad |\mathcal M|\le 9^m. \tag{4.20} $$

由 Lemma 4.4.2,

$$ \|A\| \le 2 \max_{x\in\mathcal N,\ y\in\mathcal M} |\langle Ax,y\rangle|. \tag{4.21} $$

Step 2:集中。 固定 $x\in\mathcal N$ 和 $y\in\mathcal M$。双线性型

$$ \langle Ax,y\rangle = \sum_{i=1}^m\sum_{j=1}^n A_{ij}x_jy_i $$

是独立次高斯随机变量之和。由 Proposition 2.7.1,

$$ \begin{aligned} \|\langle Ax,y\rangle\|_{\psi_2}^2 &\le C\sum_{i=1}^m\sum_{j=1}^n \|A_{ij}x_jy_i\|_{\psi_2}^2\\ &\le CK^2 \sum_{i=1}^m\sum_{j=1}^n x_j^2y_i^2\\ &= CK^2 \left(\sum_{j=1}^n x_j^2\right) \left(\sum_{i=1}^m y_i^2\right) = CK^2. \end{aligned} $$

由 Proposition 2.6.6(i),这可写成尾界

$$ \mathbb P\{|\langle Ax,y\rangle|\ge u\} \le 2\exp(-cu^2/K^2), \qquad u\ge0. \tag{4.22} $$

Step 3:Union bound。 现在通过 union bound 解除 $x,y$ 固定的限制。事件

$$ \max_{x\in\mathcal N,\,y\in\mathcal M} |\langle Ax,y\rangle| \ge u $$

表示存在 $x\in\mathcal N$ 和 $y\in\mathcal M$ 使得 $|\langle Ax,y\rangle|\ge u$。因此

$$ \begin{aligned} &\mathbb P\left\{ \max_{x\in\mathcal N,\,y\in\mathcal M} |\langle Ax,y\rangle| \ge u \right\}\\ &\qquad\le \sum_{x\in\mathcal N,\,y\in\mathcal M} \mathbb P\{|\langle Ax,y\rangle|\ge u\}. \end{aligned} $$

利用尾界 (4.22) 和 net 大小估计 (4.20),上式不超过

$$ 9^{n+m}\cdot 2\exp(-cu^2/K^2). \tag{4.23} $$

$$ u=CK(\sqrt n+\sqrt m+t). \tag{4.24} $$

则 $u^2\ge C^2K^2(n+m+t^2)$。若常数 $C$ 足够大,则指数足够大,例如 $cu^2/K^2\ge3(n+m)+t^2$。于是

$$ \mathbb P\left\{ \max_{x\in\mathcal N,\,y\in\mathcal M} |\langle Ax,y\rangle| \ge u \right\} \le 2\exp(-t^2). $$

最后结合 (4.21),得到

$$ \mathbb P\{\|A\|\ge2u\} \le 2\exp(-t^2). $$

再代入 (4.24),并把常数 $2$ 吸收到 $C$ 中,证明完成。

Remark 4.4.4 期望界

像 Theorem 4.4.3 这样的高概率界,通常可以通过 integrated tail formula(Lemma 1.6.1)转化为更简单但信息较少的期望界。尝试 Exercise 4.41,可得到

$$ \mathbb E\|A\| \le CK(\sqrt m+\sqrt n). $$
Remark 4.4.5 最优性

Theorem 4.4.3 通常是紧的,因为矩阵的算子范数至少不小于任意一行或一列的 Euclidean 范数(Exercise 4.7)。例如,如果 $A$ 的元素是 Rademacher 随机变量,那么它的列范数为 $\sqrt m$,行范数为 $\sqrt n$,因此

$$ \|A\| \ge \max(\sqrt m,\sqrt n) \ge \frac{1}{2}(\sqrt m+\sqrt n) $$

以概率 $1$ 成立。完全一般的下界见 Exercise 4.42。

Remark 4.4.6 放松独立性

Theorem 4.4.3 中的独立性假设可以放松:只需要 $A$ 的行(或列)独立,行内元素可相依;见 Exercise 4.43。

4.4.3 对称矩阵

Theorem 4.4.3 很容易推广到对称矩阵,给出高概率界

$$ \|A\|\lesssim \sqrt n. $$

Corollary 4.4.7 具有次高斯元素的对称矩阵范数

设 $A$ 是 $n\times n$ 对称随机矩阵,其对角线及其上方的元素 $A_{ij}$ 独立、均值为零且次高斯。那么对任意 $t>0$,

$$ \|A\| \le CK(\sqrt n+t) $$

以至少 $1-4\exp(-t^2)$ 的概率成立。这里 $K=\max_{i,j}\|A_{ij}\|_{\psi_2}$。

Proof 分解为上下三角部分

把 $A$ 分成上三角部分 $A^+$ 和下三角部分 $A^-$。对角线可以放到任意一边;为具体起见,把它放入 $A^+$。于是

$$ A=A^++A^-. $$

Theorem 4.4.3 可分别应用于 $A^+$ 和 $A^-$。由 union bound,以至少 $1-4\exp(-t^2)$ 的概率,这两个界同时成立:

$$ \|A^+\|\le CK(\sqrt n+t), \qquad \|A^-\|\le CK(\sqrt n+t). $$

再由三角不等式 $\|A\|\le\|A^+\|+\|A^-\|$,证明完成。

关于随机矩阵范数的更多练习见 Exercise 4.44。

4.5 应用:网络中的社区检测

随机矩阵理论有许多应用。这里给出一个网络分析中的例子。

真实世界的网络常有 communities,也就是紧密连接的节点簇。准确且高效地识别它们,是 community detection problem 的主要挑战。

4.5.1 Stochastic Block Model

我们在一个简单的两社区概率网络模型中研究 community detection。它是 Section 2.5 中介绍的 Erdős-Rényi 随机图模型的直接推广。

Definition 4.5.1 Stochastic block model

把 $n$ 个顶点分成两个大小均为 $n/2$ 的组,也就是两个 communities。构造随机图 $G$:每对顶点独立连接;如果两个顶点在同一个 community 中,则以概率 $p$ 连接;如果它们在不同 communities 中,则以概率 $q$ 连接。自环也按这个规则包含在模型中。这个随机图模型称为 stochastic block model,记为 $G(n,p,q)$。

当 $p=q$ 时,得到一个带自环的 Erdős-Rényi 模型 $G(n,p)$。但当 $p>q$ 时,community 内部的边比 community 之间的边更可能出现,于是形成 community structure;见图 4.4。

随机块模型生成的随机图
图 4.4 随机块模型 $G(n,p,q)$ 生成的图,参数为 $n=200$、$p=1/20$、$q=1/200$。

4.5.2 期望邻接矩阵包含关键信息

图 $G$ 可方便地由邻接矩阵 $A$ 表示;邻接矩阵已在 Definition 3.6.2 中引入。对随机图 $G\sim G(n,p,q)$,邻接矩阵 $A$ 是随机矩阵,因此可以用本章前面发展的工具分析它。

把 $A$ 分成确定性部分和随机部分:

$$ A=D+R, \qquad D=\mathbb EA. $$

把 $D$ 看作 signal(有信息的部分),把 $R$ 看作 noise。

为了看出为什么 $D$ 有信息,计算它的特征结构。元素 $A_{ij}$ 服从 Bernoulli 分布:根据顶点 $i$ 和 $j$ 是否属于同一个 community,分别是 $\operatorname{Ber}(p)$ 或 $\operatorname{Ber}(q)$。因此 $D$ 的元素根据 community 关系取值为 $p$ 或 $q$。例如,如果按 community 排列顶点,当 $n=4$ 时,

$$ D=\mathbb EA = \begin{bmatrix} p&p&q&q\\ p&p&q&q\\ q&q&p&p\\ q&q&p&p \end{bmatrix}. $$

更一般地,$D$ 由四个 $n/2\times n/2$ 常数块组成。它的秩为 $2$,非零特征值和对应特征向量为

$$ \lambda_1(D) = \left(\frac{p+q}{2}\right)n, \qquad u_1(D) = (1,\ldots,1,1,\ldots,1)^{\mathsf T}, $$

以及

$$ \lambda_2(D) = \left(\frac{p-q}{2}\right)n, \qquad u_2(D) = (1,\ldots,1,-1,\ldots,-1)^{\mathsf T}. $$

关键对象是第二个特征向量 $u_2(D)$,它包含了全部 community structure 信息。如果知道 $u_2(D)$,就可以根据它各坐标的符号识别两个 communities。

4.5.3 实际邻接矩阵是良好近似

但是我们不知道期望邻接矩阵 $D=\mathbb EA$,因此无法直接访问 $u_2(D)$。我们知道的是实际邻接矩阵 $A=D+R$,它是 $D$ 的噪声版本。signal $D$ 的水平为

$$ \|D\| = \lambda_1(D) \asymp n, $$

而 noise $R$ 的水平可用 Corollary 4.4.7 估计:

$$ \|R\| \le C\sqrt n \qquad \text{with probability at least }1-4e^{-n}. \tag{4.25} $$

因此,当 $n$ 很大时,noise $R$ 远小于 signal $D$。这表示 $A$ 接近 $D$,所以可以用 $A$ 代替 $D$ 来提取 community 信息。下面用 Section 4.1.6 中的矩阵扰动理论证明这一点。

4.5.4 扰动理论

对 $D$ 和 $A$ 应用 Davis-Kahan 不等式(Theorem 4.1.15),并关注第二大特征值。需要检查 $\lambda_2(D)$ 与 $D$ 的其余谱,也就是与 $0$ 和 $\lambda_1(D)$,分离得足够开。距离为

$$ \delta = \min\left(\lambda_2(D),\lambda_1(D)-\lambda_2(D)\right) = \min\left(\frac{p-q}{2},q\right)n =: \mu n. $$

回忆 (4.25) 中对 $R=A-D$ 的界,Davis-Kahan 不等式给出 $D$ 和 $A$ 的单位特征向量之间的夹角界。这里用上横线表示单位化后的向量:

$$ \sin\angle\left(\bar u_2(D),\bar u_2(A)\right) \le \frac{2\|R\|}{\delta} \lesssim \frac{\sqrt n}{\mu n} \lesssim \frac{1}{\mu\sqrt n}. $$

如果两个单位向量之间夹角的正弦很小,那么它们在差一个符号的意义下接近;见 Exercise 4.16。因此存在 $\theta\in\{-1,1\}$,使得

$$ \|\bar u_2(D)-\theta\bar u_2(A)\|_2 \lesssim \frac{1}{\mu\sqrt n}. $$

我们已经算出 $D$ 的特征向量 $u_2(D)$,但它不是单位向量;因为坐标均为 $\pm1$,所以范数为 $\sqrt n$。两边乘以 $\sqrt n$,得到

$$ \|u_2(D)-\theta u_2(A)\|_2 \lesssim \frac{1}{\mu}. $$

这意味着 $u_2(D)$ 和 $\theta u_2(A)$ 的大多数坐标符号必须一致。事实上,把上式平方写成

$$ \sum_{j=1}^n \left| u_2(D)_j-\theta u_2(A)_j \right|^2 \lesssim \frac{1}{\mu^2}, $$

并注意 $u_2(D)$ 的所有坐标都是 $\pm1$,可知每一个符号不一致的坐标至少贡献 $1$ 到这项和。因此,符号不一致的数量为 $\lesssim1/\mu^2$。

4.5.5 Spectral Clustering

总结一下,我们可以用向量 $u_2(A)$ 准确估计 $u_2=u_2(D)$;后者坐标为 $\pm1$,并识别两个 communities。这种 community detection 方法通常称为 spectral clustering。

Algorithm Spectral Clustering

Input: 图 $G$。

Output: 将 $G$ 的顶点分成两个 communities。

  1. 计算图的邻接矩阵 $A$。
  2. 计算 $A$ 第二大特征值对应的特征向量 $v_2(A)$。
  3. 根据 $v_2(A)$ 的坐标符号,把顶点分成两个 communities;零坐标可任意分配。

我们已经证明了如下结果。

Theorem 4.5.2 随机块模型的谱聚类

设 $G\sim G(n,p,q)$,并令 $\min(q,p-q)=\mu>0$。那么以至少 $1-4e^{-n}$ 的概率,spectral clustering algorithm 能识别 $G$ 的 communities,且误分类顶点数至多为 $C/\mu^2$。

总结来说,只要随机图足够稠密($q\ge\text{const}$)并且 community 内外连边概率分离良好($p-q\ge\text{const}$),spectral clustering algorithm 就能正确分类除 $O(1)$ 个顶点以外的所有顶点。条件 $q\ge\text{const}$ 并非本质;可通过 Exercise 4.45 尝试去掉它。

Remark 4.5.3 稀疏性

Theorem 4.5.2 非平凡,也就是 $C/\mu^2\le n$,所允许的最稀疏图具有期望平均度

$$ \frac{n(p+q)}{2} \asymp \sqrt n. $$

借助更多工具,我们将在 Section 5.5 处理稀疏得多的图。

4.6 次高斯矩阵的双侧界

回到 Theorem 4.4.3。它给出了 $m\times n$ 次高斯随机矩阵 $A$ 的奇异值上界:

$$ s_1(A)=\|A\| \le C(\sqrt m+\sqrt n) $$

以高概率成立。现在证明对 $A$ 整个谱更尖锐的双侧界:

$$ \sqrt m-C\sqrt n \le s_i(A) \le \sqrt m+C\sqrt n. $$

换句话说,我们将证明一个高矩阵 $(1/\sqrt m)A$(当 $m\gg n$)是近似等距映射;见 Section 4.1.7。

Theorem 4.6.1 次高斯矩阵奇异值的双侧界

设 $A$ 是 $m\times n$ 随机矩阵,其行 $A_i$ 独立、均值为零、次高斯且各向同性。那么对任意 $t\ge0$,有

$$ \sqrt m - CK^2(\sqrt n+t) \le s_n(A) \le s_1(A) \le \sqrt m + CK^2(\sqrt n+t) \tag{4.26} $$

以至少 $1-2\exp(-t^2)$ 的概率成立。这里 $K=\max_i\|A_i\|_{\psi_2}$。

查看学习笔记完整证明

我们将证明比 (4.26) 稍强的结论:

$$ \left\| \frac{1}{m}A^{\mathsf T}A-I_n \right\| \le K^2\max(\delta,\delta^2), \qquad \delta = C\left( \sqrt{\frac{n}{m}} + \frac{t}{\sqrt m} \right). \tag{4.27} $$

利用 (4.15) 可以验证 (4.27) 确实推出 (4.26)。

查看学习笔记:从矩阵偏差推出奇异值界

Proof net argument + Bernstein concentration

像 Theorem 4.4.3 一样,我们用 $\varepsilon$-net argument 证明 (4.27),但这次使用 Bernstein 集中不等式,而不是 Hoeffding 型界。

Step 1:逼近。 由 Corollary 4.2.11,可取单位球面 $S^{n-1}$ 的一个 $1/4$-net $\mathcal N$,满足

$$ |\mathcal N|\le9^n. $$

由 Lemma 4.4.2,可在 $\mathcal N$ 上估计 (4.27) 中的算子范数:

$$ \begin{aligned} \left\| \frac{1}{m}A^{\mathsf T}A-I_n \right\| &\le 2\max_{x\in\mathcal N} \left| \left\langle \left( \frac{1}{m}A^{\mathsf T}A-I_n \right)x,x \right\rangle \right|\\ &= 2\max_{x\in\mathcal N} \left| \frac{1}{m}\|Ax\|_2^2-1 \right|. \end{aligned} $$

因此,为证明 (4.27),只需证明以所需概率有

$$ \max_{x\in\mathcal N} \left| \frac{1}{m}\|Ax\|_2^2-1 \right| \le \frac{\eta}{2}, \qquad \eta:=K^2\max(\delta,\delta^2). $$

Step 2:集中。 固定 $x\in\mathcal N$,把 $\|Ax\|_2^2$ 写成独立随机变量之和:

$$ \|Ax\|_2^2 = \sum_{i=1}^m \langle A_i,x\rangle^2 =: \sum_{i=1}^m X_i^2. \tag{4.28} $$

由假设,行 $A_i$ 是独立、各向同性、次高斯的随机向量,且 $\|A_i\|_{\psi_2}\le K$。因此 $X_i=\langle A_i,x\rangle$ 是独立次高斯随机变量,满足 $\mathbb E X_i^2=1$ 且 $\|X_i\|_{\psi_2}\le K$。于是 $X_i^2-1$ 是独立、均值为零的次指数随机变量,并且

$$ \|X_i^2-1\|_{\psi_1} \le CK^2. $$

使用 Bernstein 不等式(Corollary 2.9.2),得到

$$ \begin{aligned} &\mathbb P\left\{ \left| \frac{1}{m}\|Ax\|_2^2-1 \right| \ge \frac{\eta}{2} \right\}\\ &\qquad= \mathbb P\left\{ \left| \frac{1}{m}\sum_{i=1}^m X_i^2-1 \right| \ge \frac{\eta}{2} \right\}\\ &\qquad\le 2\exp\left[ -c_1\min\left( \frac{\eta^2}{K^4}, \frac{\eta}{K^2} \right)m \right]\\ &\qquad= 2\exp(-c_1\delta^2m)\\ &\qquad\le 2\exp\left[-c_1C^2(n+t^2)\right]. \end{aligned} $$

倒数第二步使用了 $\eta/K^2=\max(\delta,\delta^2)$,因而括号中的最小值等于 $\delta^2$。最后一步来自 $\delta$ 的定义和 $(a+b)^2\ge a^2+b^2$。

Step 3:Union bound。 现在对 $x\in\mathcal N$ 做 union bound。由 $|\mathcal N|\le9^n$,

$$ \begin{aligned} &\mathbb P\left\{ \max_{x\in\mathcal N} \left| \frac{1}{m}\|Ax\|_2^2-1 \right| \ge \frac{\eta}{2} \right\}\\ &\qquad\le 9^n\cdot 2\exp\left[-c_1C^2(n+t^2)\right] \le 2\exp(-t^2), \end{aligned} $$

只要 (4.27) 中的绝对常数 $C$ 足够大。由 Step 1,这完成了 (4.27) 以及定理的证明。

Remark 4.6.2 期望形式

如 Remark 4.4.4 所述,高概率界可转化为期望界。做 Exercise 4.41 可得到 Theorem 4.6.1 的如下期望形式:

$$ \mathbb E \left\| \frac{1}{m}A^{\mathsf T}A-I_n \right\| \le CK^2 \left( \sqrt{\frac{n}{m}} + \frac{n}{m} \right). $$

Exercise 4.46 会给出 Theorem 4.6.1 的另一种证明。

4.7 应用:协方差估计与聚类

假设我们想分析一些高维数据,给定为从 $\mathbb R^n$ 中某个未知分布抽样得到的点 $X_1,\ldots,X_m$。探索这类数据的基本工具是主成分分析(PCA),我们在 Section 3.2.2 中已经接触过它。

PCA 把数据协方差矩阵的 top eigenvectors 作为 principal components。虽然我们不知道底层分布的协方差矩阵,也就是 population covariance matrix,但可以用样本 $X_1,\ldots,X_m$ 近似估计它。随后 Davis-Kahan theorem 4.1.15 可帮助估计底层分布的 principal components。

如何从数据估计协方差矩阵?令 $X$ 表示来自未知分布的随机向量。为简单起见,假设 $X$ 均值为零,并把协方差矩阵记为

$$ \Sigma=\mathbb E XX^{\mathsf T}. $$

为估计 $\Sigma$,使用由样本 $X_1,\ldots,X_m$ 计算得到的样本协方差矩阵

$$ \Sigma_m = \frac{1}{m} \sum_{i=1}^m X_iX_i^{\mathsf T}. $$

这里做的只是把总体平均替换为样本平均。由于 $X_i$ 与 $X$ 同分布,该估计是无偏的:

$$ \mathbb E\Sigma_m=\Sigma. $$

对 $\Sigma$ 的每个元素应用大数定律(Theorem 1.7.1),可得当样本量 $m\to\infty$ 时,

$$ \Sigma_m\to\Sigma \qquad \text{almost surely}. $$

这引出一个定量问题:样本量 $m$ 需要多大,才能保证

$$ \Sigma_m\approx\Sigma \qquad \text{with high probability?} $$

从维数角度看,至少需要 $m\gtrsim n$ 个样本点。下面说明 $m\asymp n$ 就足够。

Theorem 4.7.1 协方差估计

设 $X$ 是 $\mathbb R^n$ 中的次高斯随机向量。更具体地,假设存在 $K\ge1$,使得对任意 $x\in\mathbb R^n$,

$$ \|\langle X,x\rangle\|_{\psi_2} \le K\|\langle X,x\rangle\|_{L^2}. \tag{4.29} $$

那么对每个正整数 $m$,

$$ \mathbb E\|\Sigma_m-\Sigma\| \le CK^2 \left( \sqrt{\frac{n}{m}} + \frac{n}{m} \right) \|\Sigma\|. $$ 查看学习笔记完整证明
Proof 白化后应用次高斯矩阵双侧界

先把随机向量 $X,X_1,\ldots,X_m$ 放到 isotropic position。为简单起见,假设 $\Sigma$ 可逆;这个条件可像 Exercise 3.10 那样去掉。令

$$ Z=\Sigma^{-1/2}X, \qquad Z_i=\Sigma^{-1/2}X_i. $$

则 $Z,Z_1,\ldots,Z_m$ 是独立且各向同性的随机向量,并满足

$$ X=\Sigma^{1/2}Z, \qquad X_i=\Sigma^{1/2}Z_i. $$

假设 (4.29) 推出

$$ \|Z\|_{\psi_2}\le K, \qquad \|Z_i\|_{\psi_2}\le K. $$

于是

$$ \|\Sigma_m-\Sigma\| = \|\Sigma^{1/2}R_m\Sigma^{1/2}\| \le \|R_m\|\|\Sigma\|, \tag{4.30} $$

其中

$$ R_m := \frac{1}{m} \sum_{i=1}^m Z_iZ_i^{\mathsf T} - I_n. $$

考虑以 $Z_i^{\mathsf T}$ 为行的 $m\times n$ 随机矩阵 $A$。那么

$$ \frac{1}{m}A^{\mathsf T}A-I_n = \frac{1}{m} \sum_{i=1}^m Z_iZ_i^{\mathsf T} - I_n = R_m. $$

应用 Theorem 4.6.1 的期望形式(Remark 4.6.2),得到

$$ \mathbb E\|R_m\| \le CK^2 \left( \sqrt{\frac{n}{m}} + \frac{n}{m} \right). $$

代入 (4.30) 即得结论。

Remark 4.7.2 样本复杂度

Theorem 4.7.1 表明,对任意 $\varepsilon\in(0,1)$,只要样本量满足

$$ m\asymp \varepsilon^{-2}n, $$

就可以用小的相对误差估计协方差矩阵:

$$ \mathbb E\|\Sigma_m-\Sigma\| \le \varepsilon\|\Sigma\|. $$

因此,如果样本量 $m$ 与维数 $n$ 成正比,样本协方差矩阵就能很好估计总体协方差矩阵。

Remark 4.7.3 高概率界

上述论证也给出高概率界:对任意 $u\ge0$,以至少 $1-2e^{-u}$ 的概率有

$$ \|\Sigma_m-\Sigma\| \le CK^2 \left( \sqrt{\frac{n+u}{m}} + \frac{n+u}{m} \right) \|\Sigma\|. $$

请在 Exercise 4.49 中验证这一点。

我们将在 Sections 5.6 和 9.2.3 中重新讨论协方差估计,处理 heavy-tailed 和近似低维分布。目前可尝试 Exercise 4.48,以小相对误差捕捉所有一维边缘。

4.7.1 应用:点集聚类

下面用一个聚类应用说明 Theorem 4.7.1。它遵循与 Section 4.5 类似的思想,但这次是在 $\mathbb R^n$ 的点集中寻找 clusters,而不是在网络中寻找 communities。虽然 cluster 的概念没有严格定义,常识上同一 cluster 中的点应比不同 cluster 中的点更接近。

和网络一样,先建立一个简单的 $\mathbb R^n$ 中双簇点集概率模型。

Definition 4.7.4 Gaussian mixture model

按如下方式在 $\mathbb R^n$ 中生成 $m$ 个随机点:抛一枚公平硬币;若为正面,则从 $N(\mu,I_n)$ 抽取一个点;若为反面,则从 $N(-\mu,I_n)$ 抽取一个点。这个分布称为均值为 $\pm\mu$ 的 Gaussian mixture model。

等价地,考虑随机向量

$$ X=\theta\mu+g, $$

其中 $\theta$ 是 Rademacher 随机变量,$g\sim N(0,I_n)$,并且 $\theta$ 与 $g$ 独立。抽取与 $X$ 同分布的独立样本 $X_1,\ldots,X_m$,该样本即服从 Gaussian mixture model;见图 4.5。

Gaussian mixture model 的二维点云
图 4.5 从均值为 $-\mu$ 与 $\mu$ 的 Gaussian mixture model 中抽取的 $m=3000$ 个点;两个大黑点为均值,$\mu=(-1.6,0)$。

给定来自 Gaussian mixture model 的 $m$ 个样本点,目标是判断哪些点属于哪个 cluster。为此,可以使用 Section 4.5 中网络 spectral clustering algorithm 的一个变体。

为什么谱方法可能有效?注意 $X$ 的分布并非各向同性,而是在 $\mu$ 方向上被拉伸,也就是图 4.5 中的水平方向。因此,可以通过计算数据的第一主成分,即样本协方差矩阵最大特征值对应的特征向量,近似找到 $\mu$。然后把数据点投影到这个主成分上,并根据它们位于原点哪一侧进行分类。算法如下。

Algorithm Spectral Clustering for Point Sets

Input: 点 $X_1,\ldots,X_m\in\mathbb R^n$。

Output: 将点划分为两个 clusters。

  1. 计算样本协方差矩阵 $\Sigma_m=\frac{1}{m}\sum_{i=1}^m X_iX_i^{\mathsf T}$。
  2. 计算 top eigenvector $v=v_1(\Sigma_m)$。
  3. 根据 $\langle X_i,v\rangle$ 的符号,把点 $X_i$ 分到两个 clusters。
Theorem 4.7.5 Gaussian mixture model 的谱聚类

从均值为 $\mu$ 和 $-\mu$ 的 Gaussian mixture model 中抽取点 $X_1,\ldots,X_m\in\mathbb R^n$。若 $m\ge Cn$ 且 $\|\mu\|_2\ge C$,则以至少 $0.99$ 的概率,spectral clustering algorithm 能识别 communities,且误分类点数至多为 $1\%$。

现在轮到你了:请在 Exercise 4.51 中证明 Theorem 4.7.5。

值得注意的是,即使 cluster separation $\|\mu\|_2$ 远小于 cluster 直径 $\asymp\sqrt n$,仍然可以实现准确分类。

4.8 Notes

Min-max theorem(Theorem 4.1.6)也称为 Courant-Fischer-Weyl min-max principle。

Davis-Kahan theorem(Theorem 4.1.15 和 Lemma 4.1.16 的一种形式)最初由 [95] 证明,现在已成为数值分析和统计学中不可或缺的工具。这个定理有大量扩展、变体和不同证明,尤其可参见 [346, 352]、[41, Section VII.3]、[306, Chapter V]。对随机扰动而言,Davis-Kahan 还可以改进,见 [343, 120, 262, 263, 345]。

Section 4.2 引入了 covering numbers、packing numbers 和 metric entropy。进一步内容可见 [21, Chapter 4] 和 [272];带有应用视角的阐述可见 [280]。

Section 4.3.2 涵盖了纠错码的一些基本结果。书籍 [334] 提供了更系统的介绍。Theorem 4.3.5 是纠错码 rate 的经典 Gilbert-Varshamov bound 的简化版本;其反向结果(Exercise 4.33)是 Hamming bound 的简化版本。

在 Section 4.5 中,我们给出了随机矩阵理论在网络中的一个应用。关于 network analysis 这一跨学科领域的系统介绍,可参见书籍 [257]。Stochastic block models(Definition 4.5.1)由 [163] 引入。Stochastic block models 中的 community detection problem 吸引了大量研究;可参见书籍 [257]、综述 [1, 126],以及包括 [229, 353, 254, 153, 2, 48, 86, 206, 151, 175, 120] 在内的许多论文。

在 Section 4.7 中,我们讨论了 covariance estimation problem;更一般的结果会在 Sections 5.6 和 9.2.3 中出现。这个主题在高维统计中已有大量研究,参见例如 [340, 284, 190, 69, 211, 82, 320, 238, 4, 261]。

在 Section 4.7.1 中,我们给出了 Gaussian mixture models 聚类的一个应用。这个问题在统计学和计算机科学中都有研究,参见 [250, Chapter 6] 和 [181, 251, 33, 165, 19, 146, 212]。

Power method(Exercise 4.6)也称为 von Mises iteration,可用于寻找 top eigenvector。Schur bound,也称 Schur test(Exercise 4.8),最初由 I. Schur [298, p. 6] 证明。

Walsh matrices(Exercise 4.9)是 Hadamard matrices 的一个特殊情形,也就是具有 $\pm1$ 元素的正交矩阵。该练习为所有 $2^k$ 阶构造 Hadamard 矩阵,而人们猜想对所有 $4k$ 阶都存在 Hadamard 矩阵。

连接正交投影之差的范数与乘积的范数的等式(Exercise 4.12),有时称为 Krein-Krasnoselskii-Milman formula,归功于 [194]。Wedin [347] 的工作显式证明了这个公式。[264] 给出了一个几何证明,[239, p. 454] 给出了一个线性代数证明;Exercise 4.12 的提示指向后者。

Wedin theorem(Exercise 4.15 的一种形式)由 [346] 建立。

Hermitian dilation(Exercise 4.14)是一个简单但强大的技巧,有许多应用 [306, Chapter I, Section 4]。

Cut norm(Exercise 4.21)在理论计算机科学 [16]、组合学 [173]、数值逼近 [130] 等领域中很重要。Cut norm 与 $\infty\to1$ 范数的等价性(Exercise 4.21),以及计算 cut norm 的 SDP 松弛(Exercise 4.22),来自 N. Alon 和 A. Naor 的开创性论文 [16]。

$\varepsilon$-net expansion 的思想(Exercises 4.34、4.35)可追溯到 [246]。

Exercise 4.48 讨论 direction-aware covariance estimation;关于这个问题的更深结果,见 [3]。

Exercise 4.50 中低秩矩阵覆盖数的界归功于 E. Candes 和 Y. Plan [71]。

Exercises

Exercise 4.1 逆矩阵的 SVD

设 $A$ 是 $n\times n$ 矩阵,奇异值分解为

$$ A=\sum_{i=1}^n s_i u_i v_i^{\mathsf T}. $$

证明 $A$ 可逆当且仅当所有奇异值 $s_i$ 都非零。在这种情况下,验证

$$ A^{-1} = \sum_{i=1}^n s_i^{-1}v_i u_i^{\mathsf T}. $$ 查看学习笔记证明路线
Exercise 4.2 算子范数的基本性质

(a) 证明 Definition 4.1.8 中引入的算子范数 $\|A\|$ 确实是 $m\times n$ 矩阵空间上的范数。

(b) 对任意矩阵,证明

$$ \|A^{\mathsf T}\|=\|A\|. $$

(c) 对任意维度匹配的矩阵,证明

$$ \|AB\|\le \|A\|\cdot\|B\|. $$

找一个例子,使左侧为零但右侧不为零。

(d) 证明任意子矩阵的算子范数都不超过原矩阵的算子范数。

查看学习笔记证明路线
Exercise 4.3 算子范数:简单例子

算子范数通常很难用矩阵元素表示,但下面是两个例外。

(a) Rank-one:对任意 $u\in\mathbb R^m$ 和 $v\in\mathbb R^n$,验证

$$ \|uv^{\mathsf T}\| = \|uv^{\mathsf T}\|_F = \|u\|_2\|v\|_2. $$

(b) Diagonal:对任意对角矩阵 $A$,其对角元素为 $a_1,\ldots,a_n$,验证

$$ \|A\|=\max_i|a_i|. $$ 查看学习笔记证明路线
Exercise 4.4 算子范数与 Frobenius 范数

设 $A$ 是 $m\times n$ 矩阵。

(a) 如果 $A$ 的秩为 $r$,证明

$$ \|A\|\le \|A\|_F\le \sqrt r\,\|A\|. $$

说明对任意 $m,n$ 和 $1\le r\le\min(m,n)$,两个界都可以达到。

(b) 如果 $Z$ 是 $\mathbb R^n$ 中的各向同性随机向量,证明

$$ \mathbb E\|AZ\|_2^2=\|A\|_F^2. $$

(c) 如果 $B$ 是 $k\times m$ 矩阵,证明

$$ \|BA\|_F\le \|B\|\|A\|_F. $$

(d) 如果 $A$ 是对角矩阵,把 (c) 改进为

$$ \|BA\|_F \le \|B\|_{1\to2}\|A\|_F, $$

其中 $\|B\|_{1\to2}$ 表示 $B$ 的列的最大 Euclidean 范数。我们将在 Exercises 4.18、4.19 中进一步讨论 $1\to2$ 范数和一般 $p\to q$ 范数。

查看学习笔记证明路线
Exercise 4.5 奇异值的一个界

证明任意矩阵 $A$ 的奇异值 $s_i(A)$ 满足

$$ s_k(A) \le \frac{\|A\|_F}{\sqrt k}, \qquad k=1,2,\ldots. $$ 查看学习笔记证明路线
Exercise 4.6 Power method

下面是一种计算友好的方法:不显式计算谱,也能近似矩阵的算子范数。设 $A$ 是 $m\times n$ 矩阵,$x\sim N(0,I_n)$。证明以概率 $1$,

$$ \sqrt[2k]{\|(A^{\mathsf T}A)^k x\|_2} \to \|A\| \qquad \text{as }k\to\infty. $$ 查看学习笔记证明路线
Exercise 4.7 算子范数与列范数

矩阵的算子范数很难用元素表示,但本题和下一题会证明一些有用界。设 $A$ 是 $m\times n$ 矩阵,列为 $A_1,\ldots,A_n$。

(a) 证明

$$ \|A\| \ge \max_i\|A_i\|_2, $$

当 $A_i$ 两两正交时取等号。

(b) 证明如果对某个 $i$ 有 $\|A\|=\|A_i\|_2$,则 $A_i$ 与所有其他列正交。

(c) 说明对行也有同样结论。

查看学习笔记证明路线
Exercise 4.8 Schur bound

设 $A$ 是 $m\times n$ 矩阵,行为 $A_{i:}$,列为 $A_{:j}$。证明

$$ \|A\| \le \sqrt{ \max_i\|A_{i:}\|_1 \cdot \max_j\|A_{:j}\|_1 }. $$ 查看学习笔记证明路线
Exercise 4.9 Walsh 矩阵

构造“最分散”的正交矩阵,也就是所有元素绝对值都相同的正交矩阵。从 $2\times2$ 矩阵

$$ W_1= \begin{bmatrix} 1&1\\ 1&-1 \end{bmatrix} $$

开始,迭代定义分块矩阵

$$ W_{k+1} = \begin{bmatrix} W_k&W_k\\ W_k&-W_k \end{bmatrix}, \qquad k=1,2,\ldots. $$

因此 $W_2$ 是 $4\times4$ 矩阵,$W_3$ 是 $8\times8$ 矩阵,依此类推。这些具有 $\pm1$ 元素的矩阵称为 Walsh 矩阵。证明对每个 $k$,$\frac{1}{\sqrt n}W_k$ 是正交矩阵,其中 $n$ 是 $W_k$ 的阶数。

查看学习笔记证明路线
Exercise 4.10 具有 $\pm1$ 元素的矩阵

证明任意 $n\times n$、元素为 $\pm1$ 的矩阵 $A$ 都满足

$$ \sqrt n\le\|A\|\le n, $$

并说明两个界都可在无穷多个 $n$ 上达到。

查看学习笔记证明路线
Exercise 4.11 正交投影之差

证明 $\mathbb R^n$ 中任意两个正交投影 $P$ 和 $Q$ 都满足 $\|P-Q\|\le1$。如果想挑战一下,可以做下一题得到更精细的界。

查看学习笔记证明路线
Exercise 4.12 正交投影之差与乘积

对 $\mathbb R^n$ 中任意两个正交投影 $P$ 和 $Q$,证明如下结论。

(a)

$$ \|P-Q\| = \max\bigl(\|P_{\perp}Q\|,\|PQ_{\perp}\|\bigr), $$

其中 $P_{\perp}=I_n-P$ 且 $Q_{\perp}=I_n-Q$。

(b) 如果 $P$ 和 $Q$ 的秩不同,则 $\|P-Q\|=1$。

(c) 如果 $P$ 和 $Q$ 的秩相同,则

$$ \|P-Q\| = \|P_{\perp}Q\| = \|PQ_{\perp}\|. $$ 查看学习笔记证明路线
Exercise 4.13 谱投影的 Davis-Kahan

证明 Davis-Kahan theorem(Theorem 4.1.15)关于 top $k$ eigenvectors 投影的一个版本。考虑两个对称矩阵 $A$ 和 $B$,谱分解为

$$ A=\sum_i\lambda_i u_i u_i^{\mathsf T}, \qquad B=\sum_i\mu_i v_i v_i^{\mathsf T}, $$

其中特征值按非增顺序排列。证明谱投影

$$ P_A=\sum_{i=1}^k u_i u_i^{\mathsf T}, \qquad P_B=\sum_{i=1}^k v_i v_i^{\mathsf T} $$

满足

$$ \|P_A-P_B\| \le \frac{2\|A-B\|}{\lambda_k-\lambda_{k+1}}. $$ 查看学习笔记证明路线
Exercise 4.14 Hermitian dilation

把任意 $m\times n$ 矩阵 $A$ 变成对称矩阵的两种常见方式是:(1) 使用 $A^{\mathsf T}A$ 或 $AA^{\mathsf T}$;(2) 取它的 Hermitian dilation,即 $(m+n)\times(m+n)$ 分块矩阵

$$ H= \begin{bmatrix} 0&A\\ A^{\mathsf T}&0 \end{bmatrix}. $$

Hermitian dilation 有几个优点:它是线性变换,并且保持稀疏性。设

$$ A=\sum_i s_i u_i v_i^{\mathsf T} $$

是奇异值分解。证明 $H$ 的非零特征值只有 $\pm s_i$,对应特征向量为

$$ \begin{bmatrix} u_i\\ \pm v_i \end{bmatrix}. $$ 查看学习笔记证明路线
Exercise 4.15 Wedin theorem

Davis-Kahan inequality(Theorem 4.1.15)适用于对称矩阵。下面证明一般矩形矩阵的版本。设 $A$ 和 $B$ 是 $m\times n$ 矩阵,奇异值分解为

$$ A=\sum_i s_i u_i v_i^{\mathsf T}, \qquad B=\sum_i t_i w_i z_i^{\mathsf T}, $$

其中奇异值按非增顺序排列。证明 top $k$ 左奇异向量上的投影

$$ P_A=\sum_{i=1}^k u_i u_i^{\mathsf T}, \qquad P_B=\sum_{i=1}^k w_i w_i^{\mathsf T} $$

满足

$$ \|P_A-P_B\| \le \frac{2\|A-B\|}{s_k-s_{k+1}}. $$

给出 top $k$ 右奇异向量投影的类似界。并证明第 $k$ 个奇异向量的 Davis-Kahan 型版本。

查看学习笔记证明路线
Exercise 4.16 向量之间的夹角与距离

简化 Davis-Kahan inequality(Theorem 4.1.15)的结论。设两个单位向量 $u,v\in\mathbb R^n$ 的夹角(取 $0$ 到 $\pi/2$ 之间的数)很小,即

$$ \sin\angle(u,v)\le\varepsilon \qquad \text{for some }\varepsilon>0. $$

证明 $u$ 与 $v$ 在差一个符号意义下接近:

$$ \|u-\theta v\|_2 \le \sqrt2\,\varepsilon \qquad \text{for some }\theta\in\{-1,1\}. $$ 查看学习笔记证明路线
Exercise 4.17 近似投影

为 Lemma 4.1.17 增加一个等价性质。设 $A$ 是 $m\times n$ 矩阵,$m\ge n$,$\varepsilon\ge0$。

(a) 证明 $A$ 是 isometry,当且仅当 $AA^{\mathsf T}$ 是 $\mathbb R^m$ 中到某个 $n$ 维子空间上的正交投影。

(b) 更一般地,证明 $A$ 是 $\varepsilon$-approximate isometry(即 Lemma 4.1.17 中三个等价性质成立),当且仅当存在 $\mathbb R^m$ 中秩为 $n$ 的正交投影 $P$,使得

$$ \|AA^{\mathsf T}-P\|\le\varepsilon. $$ 查看学习笔记证明路线
Exercise 4.18 $p\to q$ 范数

(4.9) 中定义的 $m\times n$ 矩阵 $A$ 的算子范数,度量 $A$ 对向量的拉伸。这里使用的是 $\ell^2$ 范数;也可以使用任意两个范数,例如 $p,q\in[1,\infty]$ 的 $\ell^p$ 和 $\ell^q$。定义 $A$ 的 $\ell^p\to\ell^q$ operator norm,简称 $p\to q$ 范数,为

$$ \|A\|_{p\to q} = \max_{x\ne0} \frac{\|Ax\|_q}{\|x\|_p} = \max_{\|x\|_p=\|y\|_{q'}=1} |y^{\mathsf T}Ax|, \tag{4.31} $$

其中 $q'$ 是 (1.5) 中定义的 $q$ 的共轭指数。

(a) 把 (4.9) 中的等式推广到 $p\to q$ 范数,并说明理由。

(b) 验证这确实定义了 $m\times n$ 矩阵空间上的一个范数。

(c) Duality:证明

$$ \|A^{\mathsf T}\|_{p\to q} = \|A\|_{q'\to p'}, $$

其中 $p'$ 和 $q'$ 是 $p$ 和 $q$ 的共轭指数。

查看学习笔记证明路线
Exercise 4.19 $1\to\infty$ 和 $1\to2$ 范数

设 $A$ 是 $m\times n$ 矩阵。对大多数 $p,q$,很难用元素表达 $\|A\|_{p\to q}$,但下面是一些例外。

(a) 证明 $1\to\infty$ 范数等于元素绝对值最大值:

$$ \|A\|_{1\to\infty} = \max_{i,j}|A_{ij}|. $$

(b) 证明 $1\to2$ 范数等于列的最大 Euclidean 范数,而 $2\to\infty$ 范数等于行的最大 Euclidean 范数:

$$ \|A\|_{1\to2} = \max_{j\le n}\|A_{:j}\|_2, \qquad \|A\|_{2\to\infty} = \max_{i\le m}\|A_{i:}\|_2. $$ 查看学习笔记证明路线
Exercise 4.20 $\infty\to1$ 范数

这个范数我们已经见过:它隐含在 Grothendieck 不等式(Theorem 3.5.1)的假设中。设 $A$ 是 $m\times n$ 矩阵。

(a) 证明

$$ \|A\|_{\infty\to1} = \max_{\substack{x\in\{-1,1\}^n\\ y\in\{-1,1\}^m}} |y^{\mathsf T}Ax|. \tag{4.32} $$

(b) Duality:证明

$$ \|A\|_{\infty\to1} = \sup_{\|Z\|_{\infty}\le1,\ \operatorname{rank}(Z)=1} |\langle A,Z\rangle|, $$

其中上确界取遍所有秩一 $m\times n$ 矩阵 $Z$,其所有元素绝对值不超过 $1$;矩阵内积定义在 (4.7) 中。

(c) 说明 (4.32) 的恒等式不能推广到二次型:找出 $n\times n$ 对称矩阵 $A$,使得

$$ \max_{x\in[-1,1]^n} |x^{\mathsf T}Ax| \ne \max_{x\in\{-1,1\}^n} |x^{\mathsf T}Ax|. $$ 查看学习笔记证明路线
Exercise 4.21 Cut norm

这个范数在理论计算机科学和图论中很重要。为了求 $m\times n$ 矩阵 $A$ 的 cut norm,我们对 $A$ 的每个子矩阵求元素和,取绝对值,再对所有子矩阵最大化。形式化地,定义

$$ \|A\|_{\mathrm{cut}} = \max_{I,J} \left| \sum_{i\in I,\ j\in J} A_{ij} \right|, $$

其中最大值取遍所有 $I\subset\{1,\ldots,m\}$ 和 $J\subset\{1,\ldots,n\}$。证明 cut norm 与 $\infty\to1$ 范数在常数因子内等价:

$$ \|A\|_{\mathrm{cut}} \le \|A\|_{\infty\to1} \le 4\|A\|_{\mathrm{cut}}. $$ 查看学习笔记证明路线
Exercise 4.22 SDP 松弛

设 $A$ 是 $m\times n$ 矩阵。说明 $\|A\|_{\infty\to1}$、cut norm、$\|A\|_{\infty\to2}$ 和 $\|A\|_{2\to1}$ 都可以通过求解一个半定规划,在绝对常数因子内近似。

查看学习笔记证明路线
Exercise 4.23 $\varepsilon$-net 的传递性

若 $\mathcal N$ 是 $\mathcal M$ 的 $\varepsilon$-net,而 $\mathcal M$ 是 $K$ 的 $\delta$-net,证明 $\mathcal N$ 是 $K$ 的 $(\varepsilon+\delta)$-net。

查看学习笔记证明路线
Exercise 4.24 Packing balls

Remark 4.2.5 中我们指出:若 $\mathcal N$ 是 metric space $(T,d)$ 中的 $\varepsilon$-separated 子集,则以 $\mathcal N$ 中点为中心的闭 $\varepsilon/2$-balls 两两不交。

(a) 说明其逆命题在一般 metric space $(T,d)$ 中是假的。

(b) 证明其逆命题在任意 normed space $(T,d)$ 中是真的。

查看学习笔记证明路线
Exercise 4.25 External covering numbers

在定义 $K$ 的 covering number 时,我们要求球心 $x_i$ 位于 $K$ 中。若放宽这一点,external covering number $\overline{\mathcal N}(K,d,\varepsilon)$ 允许球心位于 $K$ 外。证明通常 covering number 与 external covering number 本质等价:

$$ \overline{\mathcal N}(K,d,\varepsilon) \le \mathcal N(K,d,\varepsilon) \le \overline{\mathcal N}(K,d,\varepsilon/2). $$ 查看学习笔记证明路线
Exercise 4.26 Covering number 的单调性

Covering number 关于 $\varepsilon$ 是单调的(为什么?),但它不一定关于集合 $K$ 单调。

(a) 举例说明一般情形下并不总有

$$ L\subset K \quad\Longrightarrow\quad \mathcal N(L,d,\varepsilon) \le \mathcal N(K,d,\varepsilon). $$

(b) 证明一个近似单调性版本:

$$ L\subset K \quad\Longrightarrow\quad \mathcal N(L,d,\varepsilon) \le \mathcal N(K,d,\varepsilon/2). $$ 查看学习笔记证明路线
Exercise 4.27 Euclidean ball 的体积:几何论证

在本题和接下来的两题中,我们用三种方式计算 $\mathbb R^n$ 中单位 Euclidean ball 的体积:几何方式、概率方式和解析方式。先从几何方式开始。

三维和二维 canonical simplices
图 4.6 $\mathbb R^3$ 与 $\mathbb R^2$ 中的 canonical simplices。

(a) Canonical simplex $\Delta_n$ 如 Figure 4.6 阴影部分所示,由 $\mathbb R^n$ 中所有非负坐标且坐标和至多为 $1$ 的点组成。证明

$$ \operatorname{Vol}(\Delta_n)=\frac1{n!}. $$

(b) 回忆 $\mathbb R^n$ 中单位 $\ell^1$ ball 为

$$ B_1^n:=\{x\in\mathbb R^n:\|x\|_1\le1\}. $$

像 Figure 4.6 那样把它分割为 simplices,证明

$$ \operatorname{Vol}(B_1^n) = \frac{2^n}{n!} \le \left(\frac{2e}{n}\right)^n. $$

(c) 像 Figure 4.6 那样,把单位 Euclidean ball $B_2^n$ 夹在适当缩放的 cube 与 $\ell^1$ ball 之间,从而推出

$$ \left(\frac2{\sqrt n}\right)^n \le \operatorname{Vol}(B_2^n) \le \left(\frac{2e}{\sqrt n}\right)^n. \tag{4.33} $$

因此,高维单位球的体积是指数级小的。虽然这些界仍有松弛,但对多数用途已经足够。

查看学习笔记证明路线
Exercise 4.28 Euclidean ball 的体积:概率论证

现在改进 (4.33) 中的上界。利用 small ball probability(Exercise 3.7)推出

$$ \operatorname{Vol}(B_2^n) \le \left(\sqrt{\frac{2\pi e}{n}}\right)^n. $$

这个界是渐近尖锐的;下一题会看到这一点。

查看学习笔记证明路线
Exercise 4.29 Euclidean ball 的体积:解析论证

最后精确计算 Euclidean ball 的体积。

(a) 设 $\|\cdot\|$ 是 $\mathbb R^n$ 上的一个范数。考虑该 normed space 的单位球

$$ B:=\{x\in\mathbb R^n:\|x\|\le1\}. $$

令 $f:\mathbb R_+\to\mathbb R_+$ 为递减、可微函数,并满足 $\lim_{x\to\infty}f(x)=0$。验证恒等式

$$ \int_{\mathbb R^n} f(\|x\|)\,dx = -\operatorname{Vol}(B) \int_0^\infty t^n f'(t)\,dt. $$

(b) 代入 $f(t)=e^{-t^2/2}$,推出

$$ \operatorname{Vol}(B_2^n) = \frac{\pi^{n/2}}{\Gamma(n/2+1)}. $$

(c) 记 $R_n$ 为体积等于 $1$ 的 Euclidean ball 的半径。证明

$$ R_n = \sqrt{\frac{n}{2\pi e}}\,(1+o(1)) \quad\text{as } n\to\infty. $$

$R_n$ 很大这一事实说明单位球 $B_2^n$ 的体积很小。

查看学习笔记证明路线
Exercise 4.30 $\ell^p$ ball 的体积

令 $p\in[1,\infty]$。

(a) 重复 Exercise 4.27 中的几何论证,推出下面的界:

$$ \left(\frac2{n^{1/p}}\right)^n \le \operatorname{Vol}(B_p^n) \le \left(\frac{2e}{n^{1/p}}\right)^n. $$

(b) 使用 Exercise 4.29 中的解析论证,精确计算体积:

$$ \operatorname{Vol}(B_p^n) = \frac{\bigl(2\Gamma(1/p+1)\bigr)^n} {\Gamma(n/p+1)}. $$ 查看学习笔记证明路线
Exercise 4.31 Lattice 是 $\varepsilon$-net

$\varepsilon$-net 的构造(Remark 4.2.7)并不高效,因为 nets 往往指数级大且难以存储。换一种方式:令 $\mathcal N$ 为缩放整数 lattice

$$ \frac{\varepsilon}{\sqrt n}\mathbb Z^n $$

中落在单位球 $B_2^n$ 内的所有点。验证 $\mathcal N$ 在 Euclidean metric 下构成 $B_2^n$ 的 $\varepsilon$-net,并且

$$ |\mathcal N| \le e^n\left(\frac2\varepsilon+1\right)^n. $$

这个界几乎匹配 Corollary 4.2.11 中的界。另请说明如何快速把单位球中的任意给定向量近似为 $\mathcal N$ 中的向量。

查看学习笔记证明路线
Exercise 4.32 Hamming cube 的 covering 与 packing numbers

通过调整 volumetric method,证明 Proposition 4.2.15。

查看学习笔记证明路线
Exercise 4.33 Error correcting codes 的限制

证明 Theorem 4.3.5 的一个逆向结果。对任意把 $k$-bit 字符串编码为 $n$-bit 字符串且能纠正 $r$ 个错误的 error correcting code,证明

$$ n-k \ge r\log_2\left(\frac nr\right). $$ 查看学习笔记证明路线
Exercise 4.34 $\varepsilon$-net expansion

$\varepsilon$-nets 可以帮助近似向量。下面说明如何用它们给出精确表示。

(a) 令 $\mathcal N$ 是 $\mathbb R^n$ 的单位球面 $S^{n-1}$ 的一个 $\varepsilon$-net。证明任意向量 $x\in S^{n-1}$ 都可以写成收敛级数

$$ x=\sum_{k=0}^\infty \lambda_k x_k \quad\text{for some coefficients }0\le\lambda_k\le\varepsilon^k, $$

其中 $x_k\in\mathcal N$。

(b) 使用 $\varepsilon$-net expansion 给出 Lemma 4.4.1 的另一个证明。

查看学习笔记证明路线
Exercise 4.35 在 $\varepsilon$-net 上计算范数

令 $x\in\mathbb R^n$,并令 $\mathcal N$ 是球面 $S^{n-1}$ 的一个 $\varepsilon$-net。证明

$$ \sup_{y\in\mathcal N}\langle x,y\rangle \le \|x\|_2 \le \frac1{1-\varepsilon} \sup_{y\in\mathcal N}\langle x,y\rangle. $$ 查看学习笔记证明路线
Exercise 4.36 在 $\varepsilon$-net 上最大化二次型

通过修改 Lemma 4.4.1 的证明,证明 Lemma 4.4.2。

查看学习笔记证明路线
Exercise 4.37 $\varepsilon$-net 上的范数偏差

令 $A$ 是 $m\times n$ 矩阵,$\mu\in\mathbb R$ 且 $\varepsilon\in[0,1/2)$。证明对球面 $S^{n-1}$ 的任意 $\varepsilon$-net $\mathcal N$,都有

$$ \sup_{x\in S^{n-1}} \bigl|\|Ax\|_2-\mu\bigr| \le \frac{C}{1-2\varepsilon} \sup_{x\in\mathcal N} \bigl|\|Ax\|_2-\mu\bigr|. $$ 查看学习笔记证明路线
Exercise 4.38 $\varepsilon$-net 的 convex hull

证明球面的 $\varepsilon$-net 会“吸收”一个略小的球,见 Figure 4.7。令 $\mathcal N$ 是单位球面 $S^{n-1}$ 的一个 $\varepsilon$-net,其中 $\varepsilon\in(0,1)$。证明

$$ (1-\varepsilon)B_2^n \subset \operatorname{conv}(\mathcal N). $$

这里一如往常,$B_2^n$ 表示 $\mathbb R^n$ 中的单位 Euclidean ball,$\operatorname{conv}(\cdot)$ 表示 Appetizer 中介绍的 convex hull。

球面 net 的凸包包含稍小球
图 4.7 球面 $\varepsilon$-net 的 convex hull 包含一个略小的虚线球,见 Exercise 4.38。
查看学习笔记证明路线
Exercise 4.39 随机点形成 $\varepsilon$-net

令 $g_1,\ldots,g_N$ 为独立的 $N(0,I_n)$ 随机向量。对任意 $R,\varepsilon>0$,证明只要 $N\ge e^{C(R,\varepsilon)n}$,则以至少 $1-e^{-cn}$ 的概率,集合

$$ \left\{ \frac{g_1}{\sqrt n},\ldots,\frac{g_N}{\sqrt n} \right\} \cap RB_2^n $$

构成 $RB_2^n$ 的一个 $\varepsilon$-net(即半径为 $R$、中心在原点的球的 $\varepsilon$-net)。这里 $C(R,\varepsilon)$ 只允许依赖于 $R$ 与 $\varepsilon$。

查看学习笔记证明路线
Exercise 4.40 过多随机点不处于 convex position

Exercise 3.23 中我们看到,如果 $N\le e^{cn}$,则独立 Gaussian 随机向量 $g_1,\ldots,g_N\sim N(0,I_n)$ 以高概率处于 convex position。证明其逆向结论:若 $N\ge e^{Cn}$,则这些向量以至少 $1-e^{-cn}$ 的概率不处于 convex position。

查看学习笔记证明路线
Exercise 4.41 随机矩阵 operator norm 的期望

(a) 从 Theorem 4.4.3 推出

$$ \mathbb E\|A\| \le CK(\sqrt m+\sqrt n). $$

(b) 从 Theorem 4.6.1 推出

$$ \mathbb E \left\| \frac1m A^{\mathsf T}A-I_n \right\| \le CK^2 \left(\sqrt{\frac nm}+\frac nm\right), $$

并且

$$ \sqrt m-CK^2\sqrt n \le \mathbb E\,s_n(A) \le \mathbb E\,s_1(A) \le \sqrt m+CK^2\sqrt n. $$ 查看学习笔记证明路线
Exercise 4.42 随机矩阵范数的下界

证明 Theorem 4.4.3 的匹配下界。令 $A$ 是 $m\times n$ 随机矩阵,其元素 $A_{ij}$ 独立、subgaussian,并满足 $\mathbb E A_{ij}^2=1$。证明对任意 $t>0$,都有

$$ \|A\| \ge \frac12(\sqrt m+\sqrt n-t) $$

以至少 $1-2\exp(-ct^2/K^4)$ 的概率成立。这里

$$ K=\max_{i,j}\|A_{ij}\|_{\psi_2}. $$ 查看学习笔记证明路线
Exercise 4.43 Subgaussian 矩阵上界:放松独立性

(a) 如果 $A$ 是一个 subgaussian matrix(即 $A$ 是 $\mathbb R^{m\times n}$ 中的 subgaussian random vector,见 Definition 3.4.1),证明 Theorem 4.4.3 在不作任何独立性假设时仍成立,其中

$$ K=\|A\|_{\psi_2}. $$

(b) 特别地,如果 $A$ 有独立、均值为零、subgaussian 的行(或列)$A_i$,证明 Theorem 4.4.3 仍成立,其中

$$ K=\max_i\|A_i\|_{\psi_2}. $$ 查看学习笔记证明路线
Exercise 4.44 随机矩阵的一些 $p\to q$ 范数

Theorem 4.4.3 研究的是随机矩阵的 $2\to2$ operator norm。那么其他 $p\to q$ operator norm 呢?我们在 Exercise 4.18 中定义了它们,并在 Exercises 4.19、4.20、4.21 中给出了一些例子。令 $A$ 是 $m\times n$ 矩阵,其元素 $A_{ij}$ 独立、均值为零、方差为一且 subgaussian。记

$$ K=\max_{i,j}\|A_{ij}\|_{\psi_2}. $$

(a) 证明

$$ \|A\|_{1\to\infty} \le CK(\sqrt{\log m}+\sqrt{\log n}). $$

(b) 证明

$$ \|A\|_{1\to2} = \|A^{\mathsf T}\|_{2\to\infty} \le \sqrt m+CK^2\sqrt{\log n}. $$

之后借助更高级的工具,我们将能够处理所有 $p\to q$ 范数(见 Exercise 8.41)。

(c) 检查 (a) 和 (b) 中的界本质上是最优的。如果 $A$ 的元素为独立标准正态随机变量,验证

$$ \|A\|_{1\to\infty} \ge c(\sqrt{\log m}+\sqrt{\log n}) \quad\text{and}\quad \|A\|_{1\to2} \ge c(\sqrt m+\sqrt{\log n}). $$ 查看学习笔记证明路线
Exercise 4.45 Community detection

Community detection 的 spectral clustering algorithm 保证(Theorem 4.5.2)不仅假设 $p-q$ 不太小,也假设跨社区边的概率 $q$ 不太小。现在去掉后一个假设。设计一个 spectral clustering algorithm 的版本,使误分类顶点数至多为

$$ \frac{C}{(p-q)^2}. $$ 查看学习笔记证明路线
Exercise 4.46 随机矩阵双边界的另一个证明

给出 Theorem 4.6.1 的一个更简单证明:使用 Theorem 3.1.1 得到 $\|Ax\|_2$ 的 concentration bound,再用 Exercise 4.37 把问题化为 net 上的 union bound。

查看学习笔记证明路线
Exercise 4.47 随机矩阵的中间 singular values

从 Theorem 4.6.1 推出如下关于 intermediate singular values 的界。对任意固定 $1\le k\le n$ 和 $t\ge0$,都有

$$ s_k(A) \ge \sqrt m-CK^2(\sqrt k+t) $$

以至少 $1-2\exp(-t^2)$ 的概率成立。

查看学习笔记证明路线
Exercise 4.48 相对误差下的 covariance estimation

给出 covariance estimation(Theorem 4.7.1)的一个更敏感版本:使用相对误差而非绝对误差。假设 $\Sigma$ 可逆,证明

$$ \mathbb E\sup_{v\in\mathbb R^n} \left| \frac{v^{\mathsf T}\Sigma_m v} {v^{\mathsf T}\Sigma v} -1 \right| \le CK^2 \left( \sqrt{\frac nm} + \frac nm \right). $$

换言之,我们得到了 $X$ 的一维 marginals 的方差的 uniform relative approximation(回忆 (3.6))。

查看学习笔记证明路线
Exercise 4.49 High-probability covariance estimation

验证 Remark 4.7.3 中提到的 covariance estimation 的 high-probability 保证。

查看学习笔记证明路线
Exercise 4.50 低秩矩阵的 covering numbers

秩为 $r$ 的 $m\times n$ 矩阵可以通过它的 SVD 用 $(m+n+1)r$ 个参数描述:左 singular vectors 需要 $mr$ 个参数,右 singular vectors 需要 $nr$ 个参数,singular values 需要 $r$ 个参数。由于 covering numbers 通常随维度指数增长(见 Corollary 4.2.11),我们可以猜测集合

$$ M_{m,n,r} = \{\,m\times n\text{ matrices }A \text{ of rank }r \text{ and }\|A\|_F\le1\,\} $$

的 covering numbers 按 $(m+n+1)r$ 指数级增长。事实确实如此。

(a) 考虑由所有列正交归一的 $m\times r$ 矩阵组成的集合 $O_{m,r}$。证明对每个 $\varepsilon>0$,都有

$$ \mathcal N(O_{m,r},\|\cdot\|_{1\to2},\varepsilon) \le \left(\frac C\varepsilon\right)^{mr}. $$

(b) 证明对每个 $\varepsilon>0$,都有

$$ \mathcal N(M_{m,n,r},\|\cdot\|_F,\varepsilon) \le \left(\frac C\varepsilon\right)^{(m+n+1)r}. $$

(c) 反过来,证明对每个 $\varepsilon>0$,都有

$$ \mathcal N(M_{m,n,r},\|\cdot\|_F,\varepsilon) \ge \left(\frac c\varepsilon\right)^{(m+n)r/2}. $$ 查看学习笔记证明路线
Exercise 4.51 Gaussian mixture model 的 spectral clustering

证明 Theorem 4.7.5,也就是 spectral clustering 可以学习 Gaussian mixture model。步骤如下。

(a) 计算 $X$ 的 covariance matrix $\Sigma$,并注意 top eigenvector $u$ 与 $\mu$ 共线。

(b) 推出 $\langle X_i,u\rangle$ 的符号可以正确分类大多数点 $X_i$。

(c) 使用 covariance estimation 的结果,说明 sample covariance matrix $\Sigma_m$ 近似 $\Sigma$。

(d) 使用 Davis-Kahan inequality(Theorem 4.1.15),推出 top eigenvector $v=v_1(\Sigma_m)$ 近似 $u$。

(e) 使用 (b),推出 $\langle X_i,v\rangle$ 的符号可以正确分类大多数点 $X_i$。

查看学习笔记证明路线

校对说明

  • 本页已经迁移 Figure 4.1-4.7。
  • 本页已覆盖 Exercises 4.1-4.51,并统一渲染为 exercise 卡片。
  • 本页的公式优先采用 MathJax 兼容写法,避免 OCR 中常见的不兼容粗体命令、断裂定界符和错误空格。
学习笔记 Ch.4 随机矩阵
第 4 章学习笔记:随机矩阵

一句话定位

第 4 章把前面学到的高维概率工具用于矩阵:用 net argument 把无穷多个方向化为有限集合,再用次高斯集中和 union bound 控制算子范数、奇异值、谱聚类和协方差估计。

本章导读

第 4 章的核心问题是:如何从固定方向上的随机集中,升级到对所有方向同时成立的矩阵范数、奇异值和协方差估计?章节安排体现了很清楚的“理论 -> 应用”路线。

章节 内容 在主线中的作用
4.1 SVD、算子范数、扰动理论 建立随机矩阵结论的线性代数语言
4.2 Nets、covering、packing 把连续方向集合离散化
4.3 纠错码 用 Hamming cube 展示 metric entropy 的信息论意义
4.4 次高斯随机矩阵范数上界 得到第一条核心矩阵控制工具
4.5 社区检测 把谱范数界用于“信号 + 噪声 + 扰动”问题
4.6 奇异值双侧界 从上界升级到所有方向近似等距
4.7 协方差估计与聚类 把奇异值界用于统计学习应用

本章读法是:先学会用线性代数描述矩阵,再用 net argument 控制所有方向,最后把谱范数和奇异值界转成数据应用中的误差保证。

本页使用方式

第 4 章的核心障碍不是“矩阵很多”,而是“要同时控制无穷多个方向”。初学者应先把线性代数对象读顺,再学 net argument 如何把无穷检查变成有限检查。

你现在卡在哪里 先看哪里 本章要建立的判断
SVD、算子范数、奇异值不稳 4.1、核心对象与符号表 能把 $\|A\|$ 读成“最大拉伸”,把 $s_{\min}(A)$ 读成“最小保真”。
不理解为什么固定方向容易、所有方向难 本章学习路线、4.2-4.4 固定 $x$ 用标量集中;所有 $x$ 需要 $\varepsilon$-net + union bound。
net argument 每次都看丢 关键定理卡片、关键定理完整证明 能说出三步:建 net、在 net 上控制、用近似引理推回整个球面。
随机矩阵定理太多 公式卡片、易混点 先分清目标:控制 $\|A\|$、控制奇异值、控制协方差偏差、控制低秩扰动。
应用部分没有直觉 stochastic block model、协方差估计、PCA 相关卡片 把每个应用翻译成“信号矩阵 + 噪声矩阵,证明噪声谱范数不太大”。
准备做题 Exercises 4.1-4.51 证明工作区 先判断题目训练的是覆盖数、谱范数、奇异值、矩阵偏差还是应用建模。

本章主线

第 4 章的主线是把第 2 章的标量集中和第 3 章的随机向量语言进一步升级为矩阵级控制。最重要的转折是:固定方向上的估计不等于对所有方向同时成立,必须借助 net argument 或类似离散化机制。

推进层 要解决的问题 关键转折 后续用途
线性代数语言 矩阵误差应该怎样度量? 用 SVD、算子范数、奇异值和谱扰动描述最大影响 所有随机矩阵定理的表达语言
离散化 怎样从固定方向控制所有方向? 用 $\varepsilon$-net 把球面上的无穷检查变成有限检查 算子范数上界、奇异值界、协方差估计
矩阵集中 随机矩阵的范数通常多大? 固定 $x,y$ 后退回标量集中,再 union bound over net 控制噪声矩阵的基本工具
近似等距 随机矩阵是否保持所有向量长度? 控制 $s_{\min}(A)$ 和 $s_{\max}(A)$,把矩阵看成几何嵌入 Johnson-Lindenstrauss、压缩感知、协方差估计
应用 谱范数界怎样变成统计结论? 写成“信号 + 噪声”,再用谱扰动读出结构 社区检测、PCA、聚类

本章学习路线

先抓住一个问题
随机矩阵会不会把所有向量都差不多保持住?

第 4 章的多数定理都在回答这个问题:矩阵 $A$ 作用在一个固定向量上不难分析,难点是要同时控制所有方向。学习时先围绕“固定方向 → 所有方向 → 数据任务”这条线走。

初学者先不急着全背
  1. 先会读 $\|A\|$、$s_i(A)$ 和谱间隙。
  2. 再理解 net argument 为什么能处理“所有方向”。
  3. 最后看它如何变成社区检测、协方差估计和聚类。
1

先把矩阵看成“拉伸机器”

第 4.1 节只需要先建立三个词的直觉:SVD 告诉你矩阵沿哪些方向拉伸,operator norm 是最大拉伸,Davis-Kahan 告诉你小扰动会不会改变重要方向。

SVD $s_1(A)=\|A\|$ Weyl Davis-Kahan
这一层要会问 $A$ 最多把向量拉长多少?如果 $A$ 加一点噪声,主方向还能不能信?
2

再学会把“所有方向”变成“有限检查”

第 4.2-4.4 节的核心不是覆盖数本身,而是这个套路:先固定 $x,y$ 做集中,再把球面换成有限的 $\varepsilon$-net,最后用近似引理推回整个球面。

$\varepsilon$-net covering number union bound net argument
这一层要会问 固定方向的概率界,付出多少 net 大小成本后,能同时覆盖所有方向?
3

最后把范数界翻译成数据结论

第 4.5-4.7 节反复使用同一个模板:把观测矩阵写成“信号 + 噪声”,先控制噪声矩阵范数,再用谱扰动把矩阵误差转成结构恢复误差。

$A=D+R$ 噪声范数 谱间隙 聚类 / PCA
这一层要会问 信号矩阵的谱结构是什么?噪声相对谱间隙够不够小?
固定方向 集中不等式 $\varepsilon$-net 算子范数 / 奇异值 谱方法应用

分层阅读路线

层次 先抓什么 推荐入口 暂时怎么处理
第一遍:主线阅读 SVD、operator norm、net argument、随机矩阵范数和奇异值 本章主线、关键定理卡片 先会说清“固定方向 -> net -> 所有方向”。
第二遍:证明精读 net 近似引理、最小奇异值、Davis-Kahan、Wedin、低秩覆盖 关键定理完整证明、正文隐藏验证补全 把线性代数条件和概率界的连接补完整。
第三遍:习题与应用 covering number、谱范数、SBM、协方差估计、PCA/聚类 Exercises 4.1-4.51 先识别题目是线性代数、离散化、矩阵集中还是应用建模。
专题回看 随机矩阵、谱方法、协方差估计、PCA 第 5/6/9 章相关路线 为矩阵 Bernstein、Hanson-Wright 和 matrix deviation 做准备。

初学者补充:工具链与应用读法

下面把本章容易断开的几个技术环节展开:线性代数工具、社区检测流程、奇异值双边界、协方差估计和高斯混合聚类。这里不再重复章节顺序,而是帮助初学者把每个环节读成可操作的证明或算法模板。

Tool Chain 4.1 线性代数工具链
工具 本章用途 初学者读法
SVD 把矩阵作用拆成正交方向上的拉伸 先记 $A=\sum_i s_i u_iv_i^{\mathsf T}$,随机矩阵问题常常就是估计 $s_i(A)$
Min-max 把特征值和奇异值写成优化问题 最大方向、最小方向和谱间隙都从这里来
算子范数 度量噪声矩阵最大影响 读成“最坏方向上的放大倍数”
低秩逼近 解释 PCA、谱聚类和信号矩阵 信号常在少数主方向上,噪声用范数控制
Davis-Kahan / Wedin 把矩阵误差转成特征向量或奇异向量误差 永远同时检查噪声范数和谱间隙
近似等距 连接奇异值双侧界和协方差估计 所有方向长度几乎保持,等价于所有奇异值都接近同一尺度
Application Pattern 4.5 社区检测的谱方法流程

随机块模型不是独立应用,它是本章“信号 + 噪声 + 扰动”的第一个完整例子。

  1. 把邻接矩阵写成 $A=D+R$,其中 $D=\mathbb EA$ 是信号,$R=A-D$ 是噪声。
  2. 先读 $D$ 的谱结构:社团方向对应一个显著特征向量。
  3. 用随机矩阵范数界证明 $\|R\|$ 不大。
  4. 用 Davis-Kahan 把 $\|R\|$ 和谱间隙转成特征向量误差。
  5. 最后按经验特征向量的符号分类节点。
判断句:只要噪声范数小于信号谱间隙,谱聚类就能把社团方向读出来。
跳到社区检测补充证明
Proof Chain 4.6 次高斯矩阵双边界怎么读

Theorem 4.6.1 的目标不是只给出上界,而是同时控制所有奇异值:

$$ \sqrt m-CK^2(\sqrt n+t) \le s_n(A) \le s_1(A) \le \sqrt m+CK^2(\sqrt n+t). $$

这句话的几何意思是:当 $m$ 足够大时,$\frac1{\sqrt m}A$ 在所有方向上近似保持长度。

证明路线 控制对象 适合解决的问题
原证明 $\left\|\frac1mA^{\mathsf T}A-I_n\right\|$ 协方差估计、近似等距、奇异值双边界
Exercise 4.46 的证明 $\sup_{x\in S^{n-1}}\left|\|Ax\|_2-\sqrt m\right|$ 用 Theorem 3.1.1 的长度集中快速得到双边界

原证明更强,因为它直接控制样本协方差矩阵;练习路线更直观,因为它先固定一个方向 $x$,把 $Ax$ 看成随机向量,再用 net argument 推到所有方向。

跳到 Theorem 4.6.1 完整证明 跳到 Exercise 4.46 另一种证明
Why Sample Size 4.7 为什么协方差估计需要 $m\gtrsim n$

样本协方差是

$$ \Sigma_m=\frac1m\sum_{i=1}^m X_iX_i^{\mathsf T}. $$

每个 $X_iX_i^{\mathsf T}$ 的秩至多为 $1$,所以

$$ \operatorname{rank}(\Sigma_m)\le m. $$

如果总体协方差是满秩的,例如 $\Sigma=I_n$,但 $m<n$,那么 $\Sigma_m$ 必然有非零向量 $v$ 落在核空间中,使得 $\Sigma_m v=0$。这时

$$ \|(\Sigma_m-I_n)v\|_2=\|v\|_2, $$

因此算子范数误差不可能小。这个秩原因说明,想在所有方向上估计协方差,样本量至少要和维数同阶。

读法:大数定律只保证每个固定元素收敛;高维协方差估计要同时控制所有方向,所以需要随机矩阵奇异值工具。
跳到 Theorem 4.7.1 完整证明 跳到高概率形式
Algorithm Card 4.7.1 高斯混合模型的谱聚类

高斯混合模型可以写成

$$ X=\theta\mu+g, \qquad \theta\in\{-1,1\}, \qquad g\sim N(0,I_n). $$

总体协方差为 $I_n+\mu\mu^{\mathsf T}$,所以最大特征向量就是 $\mu$ 的方向。谱聚类算法读法如下:

  1. 输入样本 $X_1,\ldots,X_m$。
  2. 计算样本协方差 $\Sigma_m=\frac1m\sum_iX_iX_i^{\mathsf T}$。
  3. 取 $\Sigma_m$ 的最大特征向量 $v$。
  4. 按 $\langle X_i,v\rangle$ 的符号给样本分簇。

证明路线和 4.5 完全平行:协方差估计说明 $\Sigma_m$ 接近 $\Sigma$,Davis-Kahan 说明 $v$ 接近 $\mu/\|\mu\|_2$,最后一维投影的 Gaussian 尾界控制误分类比例。

跳到 Gaussian mixture 完整证明

核心对象与符号表

符号 / 对象 在原书中的角色 学习时要抓住的意思
$A\in\mathbb R^{m\times n}$ 随机矩阵或数据矩阵。 研究它如何拉伸所有向量。
$s_i(A)$ 奇异值。 描述矩阵沿主方向的拉伸;最大奇异值是算子范数。
$\|A\|$ 算子范数。 最大拉伸,也是噪声矩阵最常用的度量。
$\|A\|_F$ Frobenius 范数。 所有奇异值平方和的平方根,衡量总能量。
$\varepsilon$-net 有限近似集合。 把球面上的无限方向变成有限检查。
covering number 覆盖一个集合所需点数。 union bound 的复杂度成本。
net argument 固定方向、net、union bound、逼近。 第 4 章证明随机矩阵范数的核心套路。
Davis-Kahan 谱子空间扰动界。 把矩阵噪声范数转成特征向量/子空间误差。
$A=D+R$ 信号加噪声分解。 社区检测、谱聚类和 PCA 估计的统一模板。
sample covariance $\frac1N\sum_iX_iX_i^{\mathsf T}$。 用随机矩阵奇异值控制协方差估计误差。

关键定理卡片

Theorem 4.1.1 SVD

任意矩阵都可分解为

$$ A=\sum_i s_i u_i v_i^{\mathsf T}. $$
阅读重点:随机矩阵的核心问题通常就是控制这些 $s_i$。
查看完整证明
Definition 4.1.8 算子范数

算子范数是矩阵对单位向量的最大拉伸:

$$ \|A\|=\sup_{x\in S^{n-1}}\|Ax\|_2. $$
记忆方式:最大奇异值 $s_1(A)$。
Theorem 4.1.15 Davis-Kahan

谱子空间扰动满足

$$ \text{误差}\lesssim\frac{\|E\|}{\text{谱间隙}}. $$
使用场景:社区检测、PCA、谱聚类;先控噪声矩阵范数,再转成特征向量误差。
查看完整证明
Definition 4.2.1 $\varepsilon$-net

用有限点集近似连续集合。它把“对所有方向成立”变成“对有限多个方向成立 + 近似误差”。

Theorem 4.4.3 次高斯随机矩阵范数上界

若 $A$ 是 $m\times n$ 独立次高斯元素矩阵,则

$$ \|A\|\lesssim K^2(\sqrt m+\sqrt n) $$

以高概率成立。

查看完整证明
Theorem 4.6.1 奇异值双侧界

次高斯随机矩阵不仅最大奇异值可控,最小奇异值也可控;这说明 $A$ 在所有方向上近似等距。

查看完整证明
Theorem 4.7.1 协方差估计

样本协方差的误差可由样本矩阵的奇异值界控制。

使用场景:PCA、统计学习、聚类和协方差矩阵估计。
查看完整证明

证明模板:net argument

Proof Pattern 从固定方向到所有方向
  1. 固定方向 $x,y$,用集中不等式控制 $\langle Ax,y\rangle$ 或 $\|Ax\|_2$。
  2. 构造球面的 $\varepsilon$-net,大小约为 $(C/\varepsilon)^n$。
  3. 对 net 中所有点使用 union bound。
  4. 用近似引理把 net 上的控制扩展到整个球面。

关键定理完整证明

Complete Proof Theorem 4.1.1:奇异值分解
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明任意矩阵都可写成正交方向上的拉伸和。

证明思路

对 $A^{\mathsf T}A$ 做谱分解,右奇异向量来自其特征向量,奇异值是特征值平方根,左奇异向量由 $u_i=Av_i/s_i$ 给出。

完整证明

设 $A$ 是 $m\times n$ 矩阵。矩阵 $A^{\mathsf T}A$ 对称半正定,所以存在 $\mathbb R^n$ 的正交基 $v_i$ 和特征值 $\lambda_i\ge0$,使得

$$ A^{\mathsf T}Av_i=\lambda_i v_i. $$

令 $s_i=\sqrt{\lambda_i}$,并设正奇异值个数为 $r$。当 $1\le i\le r$ 时定义 $u_i=Av_i/s_i$,则

$$ \langle u_i,u_j\rangle = \frac{\langle Av_i,Av_j\rangle}{s_is_j} = \frac{\langle v_i,A^{\mathsf T}Av_j\rangle}{s_is_j} = \delta_{ij}. $$

于是 $Av_i=s_i u_i$。若 $i>r$,则 $s_i=0$ 且 $\|Av_i\|_2^2=\langle A^{\mathsf T}Av_i,v_i\rangle=0$,所以 $Av_i=0$。对任意 $x=\sum_i\alpha_iv_i$,

$$ Ax=\sum_i\alpha_iAv_i=\sum_i s_i u_i\langle v_i,x\rangle. $$

因此作为算子,

$$ A=\sum_{i=1}^r s_i u_iv_i^{\mathsf T}, $$

这就是紧 SVD。若需要完整 SVD,可把 $u_1,\ldots,u_r$ 扩充为 $\mathbb R^m$ 的正交基,并把其余奇异值补为 $0$。

Complete Proof Theorem 4.1.15:Davis-Kahan 的读法
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明单个特征向量的扰动由噪声范数和原矩阵谱间隙控制。

证明思路

先用 Weyl 不等式保证扰动后目标特征值仍与其他谱分离;再用谱投影分离引理,把投影差控制为 $\|A-B\|$ 与谱间隙之比。

完整证明

先证明需要的谱投影分离引理。设 $A=\sum_i\lambda_i u_iu_i^{\mathsf T}$,$B=\sum_j\mu_jv_jv_j^{\mathsf T}$。令

$$ P=\sum_{\lambda_i\in I}u_iu_i^{\mathsf T}, \qquad Q=\sum_{\mu_j\in J}v_jv_j^{\mathsf T}, $$

其中 $I,J\subset\mathbb R$ 相距至少 $\delta$,且 $I$ 是区间。令 $H=B-A$。给 $A,B$ 同时减去某个标量矩阵不改变 $H$、$P$、$Q$,所以可把 $I$ 写成 $[-r,r]$,并使得 $\mu_j\in J$ 时 $|\mu_j|\ge r+\delta$。

于是

$$ \|H\|\ge\|QHP\| = \|QBP-QAP\| \ge \|QBP\|-\|QAP\|. $$

因为 $Q$ 与 $B$ 交换,且 $B$ 在 $Q$ 的像空间上至少拉伸 $r+\delta$ 倍,

$$ \|QBP\|=\|BQP\|\ge(r+\delta)\|QP\|. $$

另一方面,$AP=PAP$,且 $\|AP\|\le r$,因此

$$ \|QAP\|=\|QPAP\|\le r\|QP\|. $$

合并三式得到

$$ \|QP\|\le\frac{\|A-B\|}{\delta}. $$

现在回到单个特征向量。令 $\varepsilon=\|A-B\|$。若 $\varepsilon>\delta/2$,则右侧 $2\varepsilon/\delta>1$,而 $\sin\angle(u_k,v_k)\le1$,结论平凡。下面设 $\varepsilon\le\delta/2$。

由 Weyl 不等式,$|\lambda_j-\mu_j|\le\varepsilon$。因此对所有 $j\ne k$,

$$ |\lambda_k-\mu_j| \ge |\lambda_k-\lambda_j|-|\lambda_j-\mu_j| \ge \delta-\varepsilon \ge \delta/2. $$

把上面的谱投影引理用于 $I=\{\lambda_k\}$ 和 $J=\{\mu_j:j\ne k\}$。此时 $P=u_ku_k^{\mathsf T}$,而 $Q$ 是 $v_k^\perp$ 上的投影,故

$$ \|QP\|\le\frac{2\varepsilon}{\delta}. $$

最后直接计算

$$ \|QP\|=\|Qu_k\|_2=\sin\angle(u_k,v_k). $$

于是

$$ \sin\angle(u_k,v_k) \le \frac{2\|A-B\|}{\delta}. $$
Complete Proof Theorem 4.4.3:次高斯随机矩阵范数上界
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明独立次高斯元素矩阵的算子范数量级为 $\sqrt m+\sqrt n$。

证明思路

算子范数是对所有 $x,y$ 的双线性型上确界。先在固定 $x,y$ 上使用次高斯集中;再取两个球面的 $1/4$-net,对 net 上所有点 union bound;最后用 net 逼近引理回到整个球面。

完整证明

取 $S^{n-1}$ 的 $1/4$-net $\mathcal N$ 和 $S^{m-1}$ 的 $1/4$-net $\mathcal M$,满足

$$ |\mathcal N|\le9^n,\qquad |\mathcal M|\le9^m. $$

由 net 版本的算子范数估计,

$$ \|A\|\le2\max_{x\in\mathcal N,\ y\in\mathcal M}|\langle Ax,y\rangle|. $$

固定 $x,y$ 时,

$$ \langle Ax,y\rangle = \sum_{i=1}^m\sum_{j=1}^n A_{ij}x_jy_i $$

是独立次高斯变量的线性组合。由独立次高斯和的范数估计,

$$ \|\langle Ax,y\rangle\|_{\psi_2}^2 \le C\sum_{i,j}x_j^2y_i^2\|A_{ij}\|_{\psi_2}^2 \le CK^2. $$

因此对所有 $u\ge0$,

$$ \mathbb P\{|\langle Ax,y\rangle|>u\}\le2e^{-cu^2/K^2}. $$

对 $\mathcal N\times\mathcal M$ 做并集界,得到

$$ \mathbb P\left\{ \max_{x\in\mathcal N,y\in\mathcal M} |\langle Ax,y\rangle|>u \right\} \le 2\cdot9^{m+n}e^{-cu^2/K^2}. $$

取 $u=C_0K(\sqrt m+\sqrt n+t)$,并把 $C_0$ 取为足够大的绝对常数,就有

$$ 2\cdot9^{m+n}e^{-cu^2/K^2} \le 2e^{-t^2}. $$

因此在概率至少 $1-2e^{-t^2}$ 的事件上,

$$ \|A\| \le 2u \le CK(\sqrt m+\sqrt n+t). $$

这就是 Theorem 4.4.3。

Complete Proof Theorem 4.6.1:奇异值双侧界
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明各向同性次高斯行矩阵在所有方向上近似保持长度。

证明思路

对固定方向 $x$,$\|Ax\|_2^2$ 是独立次指数变量之和,Bernstein 给出集中。再对球面 net 做 union bound,并用逼近把 net 控制推广到所有方向。

完整证明

设 $A$ 的行 $A_i$ 独立、各向同性、次高斯,且 $K=\max_i\|A_i\|_{\psi_2}$。先证明更强的矩阵偏差界

$$ \left\| \frac1m A^{\mathsf T}A-I_n \right\| \le K^2\max(\delta,\delta^2), \qquad \delta=C\left(\sqrt{\frac nm}+\frac{t}{\sqrt m}\right) $$

以至少 $1-2e^{-t^2}$ 的概率成立。

取 $S^{n-1}$ 的 $1/4$-net $\mathcal N$,满足 $|\mathcal N|\le9^n$。由二次型 net 逼近引理,

$$ \left\| \frac1m A^{\mathsf T}A-I_n \right\| \le 2\max_{x\in\mathcal N} \left| \frac1m\|Ax\|_2^2-1 \right|. $$

下面固定 $x\in\mathcal N$。随机变量 $X_i=\langle A_i,x\rangle$ 独立,并由各向同性满足 $\mathbb E X_i^2=1$,且 $\|X_i\|_{\psi_2}\le K$。因此

$$ \|Ax\|_2^2=\sum_{i=1}^m\langle A_i,x\rangle^2 $$

且 $X_i^2-1$ 是独立中心化次指数变量,$\|X_i^2-1\|_{\psi_1}\le CK^2$。令

$$ \eta=K^2\max(\delta,\delta^2). $$

Bernstein 不等式给出

$$ \begin{aligned} \mathbb P\left\{ \left|\frac1m\|Ax\|_2^2-1\right| \ge \frac\eta2 \right\} &\le 2\exp\left[ -c_1m\min\left(\frac{\eta^2}{K^4},\frac{\eta}{K^2}\right) \right]\\ &= 2\exp(-c_1m\delta^2). \end{aligned} $$

等号处使用了 $\eta/K^2=\max(\delta,\delta^2)$,因此最小值为 $\delta^2$。又由 $\delta$ 的定义,$m\delta^2\ge C^2(n+t^2)$。对 $\mathcal N$ 做 union bound,

$$ \mathbb P\left\{ \max_{x\in\mathcal N} \left| \frac1m\|Ax\|_2^2-1 \right| \ge \frac\eta2 \right\} \le 9^n\cdot2e^{-c_1C^2(n+t^2)} \le 2e^{-t^2} $$

只要 $C$ 足够大。结合 net 逼近引理,得到上述矩阵偏差界。

最后把矩阵偏差转成奇异值界。矩阵 $m^{-1}A^{\mathsf T}A$ 的特征值是 $s_i(A)^2/m$。若

$$ \left\|\frac1m A^{\mathsf T}A-I_n\right\|\le\eta, $$

则所有 $i$ 满足

$$ 1-\eta\le \frac{s_i(A)^2}{m}\le1+\eta. $$

开方并使用 $\eta=K^2\max(\delta,\delta^2)$(下界为负时结论平凡),得到

$$ \sqrt m-CK^2(\sqrt n+t) \le s_n(A) \le s_1(A) \le \sqrt m+CK^2(\sqrt n+t). $$
Complete Proof Theorem 4.7.1:协方差估计
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:在一般协方差 $\Sigma$ 下,用白化和随机矩阵双侧界证明样本协方差的期望算子范数误差。

证明思路

先把 $X$ 白化为各向同性向量 $Z=\Sigma^{-1/2}X$。白化后的样本协方差偏差由 Theorem 4.6.1 的期望形式控制;再乘回 $\Sigma^{1/2}$,算子范数最多多出 $\|\Sigma\|$。

完整证明

设 $X_1,\ldots,X_m$ 是 $X$ 的独立样本,$\Sigma=\mathbb E XX^{\mathsf T}$,样本协方差为

$$ \Sigma_m=\frac1m\sum_{i=1}^m X_iX_i^{\mathsf T}. $$

先假设 $\Sigma$ 可逆。令 $Z=\Sigma^{-1/2}X$,$Z_i=\Sigma^{-1/2}X_i$。则 $Z_i$ 独立、各向同性,并且由假设

$$ \|\langle Z,x\rangle\|_{\psi_2} = \|\langle X,\Sigma^{-1/2}x\rangle\|_{\psi_2} \le K\|\langle X,\Sigma^{-1/2}x\rangle\|_{L^2} = K\|x\|_2, $$

所以 $\|Z_i\|_{\psi_2}\le K$。记

$$ R_m=\frac1m\sum_{i=1}^m Z_iZ_i^{\mathsf T}-I_n. $$

把 $Z_i^{\mathsf T}$ 作为行组成 $m\times n$ 矩阵 $A$,则 $R_m=m^{-1}A^{\mathsf T}A-I_n$。由 Theorem 4.6.1 的期望形式,

$$ \mathbb E\|R_m\| \le CK^2\left(\sqrt{\frac nm}+\frac nm\right). $$

另一方面,$X_i=\Sigma^{1/2}Z_i$,所以

$$ \Sigma_m-\Sigma = \Sigma^{1/2}R_m\Sigma^{1/2}. $$

于是

$$ \|\Sigma_m-\Sigma\| \le \|\Sigma^{1/2}\|^2\|R_m\| = \|\Sigma\|\,\|R_m\|. $$

取期望得到

$$ \mathbb E\|\Sigma_m-\Sigma\| \le CK^2\left(\sqrt{\frac nm}+\frac nm\right)\|\Sigma\|. $$

若 $\Sigma$ 不可逆,可在 $\operatorname{range}(\Sigma)$ 上做同样白化,或用 $\Sigma+\varepsilon I$ 近似后令 $\varepsilon\downarrow0$;结论不变。

章节阅读路径

小节 读法
4.1 重点掌握 SVD、算子范数、min-max、Weyl、Davis-Kahan。
4.2 把 covering/packing 当作“连续集合离散化”的工具。
4.3 通过 Hamming cube 理解 metric entropy。
4.4 学会随机矩阵范数上界的 net argument。
4.5 看懂“期望矩阵 + 噪声矩阵 + 扰动理论”的社区检测套路。
4.6 理解双侧奇异值界为何意味着近似等距。
4.7 把双侧界用于协方差估计和聚类。

易混点

易混点 正确读法
算子范数只是最大元素大小 错。算子范数是所有单位方向上的最大拉伸,通常比最大元素更全局。
固定 $x$ 的集中等于所有 $x$ 同时集中 不等于。必须支付 net size 的复杂度成本。
net argument 只是技术细节 它是第 4 章从标量/向量集中走向矩阵范数的核心桥梁。
Davis-Kahan 只需要噪声小 还需要信号谱间隙足够大;误差尺度是 $\|E\|/\text{gap}$。
协方差估计只是大数定律 高维中要控制所有方向,所以需要样本矩阵奇异值或算子范数界。

公式卡片

工具 公式 用途
SVD $A=\sum_i s_i u_i v_i^{\mathsf T}$ 分解矩阵作用
算子范数 $\|A\|=\sup_{x\in S^{n-1}}\|Ax\|_2$ 最大拉伸
Frobenius $\|A\|_F^2=\sum_i s_i^2$ 总能量
Net 大小 $\mathcal N(B_2^n,\varepsilon B_2^n)\le(3/\varepsilon)^n$ union bound 成本
随机矩阵范数 $\|A\|\lesssim\sqrt m+\sqrt n$ 次高斯矩阵
Davis-Kahan $\text{误差}\lesssim\|E\|/\text{gap}$ 谱方法稳定性

正文隐藏验证补全

Hidden Check 4.1.1:SVD 证明中 $m\ge n$ 用在哪里
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:解释原书在证明 SVD 时为什么先假设 $m\ge n$,以及该假设实际用于哪一步。

完整证明

证明从 $A^{\mathsf T}A$ 的谱分解得到 $n$ 个右奇异向量 $v_1,\ldots,v_n\in\mathbb R^n$。当 $s_i>0$ 时,左奇异向量必须定义为

$$ u_i=\frac{Av_i}{s_i}, $$

这些向量自动正交归一。问题出现在 $s_i=0$ 的方向:这时 $Av_i=0$,等式 $Av_i=s_i u_i$ 对任意单位向量 $u_i$ 都成立,但我们还要让全部 $u_1,\ldots,u_n$ 正交归一。

因此,原证明中真正用到 $m\ge n$ 的地方,是“任选 $u_i$,只要保持正交归一性”这一步。只有当 $\mathbb R^m$ 至少有 $n$ 维时,才一定能把已得到的左奇异向量补成 $n$ 个正交归一向量。

如果 $m<n$,不能在 $\mathbb R^m$ 中选出 $n$ 个正交归一向量。此时可以对 $A^{\mathsf T}$ 应用 $n\ge m$ 的情形,或直接从 $AA^{\mathsf T}$ 出发先构造 $m$ 个左奇异向量。最终得到的紧 SVD 始终只需要 $r=\min(m,n)$ 个方向。

Hidden Check 4.2.14:Hamming 距离确实是度量
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:验证 $d_H(x,y)=\#\{i:x(i)\ne y(i)\}$ 满足度量的三个条件。

完整证明

非负性来自 $d_H(x,y)$ 是一个集合的基数。并且 $d_H(x,y)=0$ 当且仅当没有坐标不同,也就是 $x=y$。

对称性也直接成立:$x(i)\ne y(i)$ 与 $y(i)\ne x(i)$ 是同一件事,所以 $d_H(x,y)=d_H(y,x)$。

最后证明三角不等式。对每个坐标 $i$,若 $x(i)\ne z(i)$,则不可能同时有 $x(i)=y(i)$ 且 $y(i)=z(i)$。因此指标函数逐坐标满足

$$ \mathbf 1_{\{x(i)\ne z(i)\}} \le \mathbf 1_{\{x(i)\ne y(i)\}} + \mathbf 1_{\{y(i)\ne z(i)\}}. $$

对 $i=1,\ldots,n$ 求和,得到

$$ d_H(x,z)\le d_H(x,y)+d_H(y,z). $$

所以 $d_H$ 是 $\{0,1\}^n$ 上的度量。

Hidden Check 4.6:由 (4.27) 推出奇异值界 (4.26)
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:补全正文中“利用 (4.15) 可以验证 (4.27) 推出 (4.26)”的代数步骤。

完整证明

把矩阵缩放为

$$ B=\frac1{\sqrt m}A. $$

那么

$$ B^{\mathsf T}B-I_n = \frac1mA^{\mathsf T}A-I_n, $$

且 $B$ 的奇异值是 $s_i(A)/\sqrt m$。由 (4.27),

$$ \|B^{\mathsf T}B-I_n\| \le K^2\max(\delta,\delta^2), \qquad \delta=C\left(\sqrt{\frac nm}+\frac{t}{\sqrt m}\right). $$

可把 $K$ 替换为 $\max(K,1)$,这只会放宽结论。于是令 $\Delta=K^2\delta$,有

$$ K^2\max(\delta,\delta^2) \le \max(\Delta,\Delta^2). $$

现在应用 (4.15) 到矩阵 $B$,得到所有奇异值满足

$$ 1-\Delta \le s_n(B) \le s_1(B) \le 1+\Delta. $$

乘以 $\sqrt m$,并代入 $\Delta$,得到

$$ \sqrt m-CK^2(\sqrt n+t) \le s_n(A) \le s_1(A) \le \sqrt m+CK^2(\sqrt n+t). $$

若左端为负,该下界当然仍成立;这就是 (4.26)。

Exercises 4.1-4.51 证明工作区

批次 建议目标
4.1-4.22 线性代数、扰动理论、operator norms、cut norm、SDP。
4.23-4.33 nets、packing、covering、体积、Hamming cube。
4.34-4.44 net argument 和随机矩阵范数变体。
4.45-4.51 社区检测、协方差估计、低秩矩阵 covering、Gaussian mixture 聚类。

本区先把每道习题的证明入口固定下来,方便从译文、学习笔记和并排阅读页跳转。当前 Exercises 4.1-4.51 均已升级为“证明目标 + 证明思路 + 完整证明”。

Exercise 4.1 逆矩阵的 SVD
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:说明可逆性等价于所有奇异值非零,并写出 $A^{-1}$ 的 SVD 形式。

证明思路

SVD 把 $A$ 写成 $U\Sigma V^{\mathsf T}$。矩阵可逆当且仅当对角矩阵 $\Sigma$ 可逆;逆矩阵就是 $V\Sigma^{-1}U^{\mathsf T}$。

完整证明

由 SVD,$A=U\Sigma V^{\mathsf T}$,其中 $\Sigma=\operatorname{diag}(s_1,\ldots,s_n)$,$U,V$ 正交。由于正交矩阵可逆,$A$ 可逆当且仅当 $\Sigma$ 可逆,也就是当且仅当所有 $s_i\ne0$。

在这种情况下,

$$ A^{-1} = V\Sigma^{-1}U^{\mathsf T} = \sum_{i=1}^n s_i^{-1}v_i u_i^{\mathsf T}. $$
Exercise 4.2 算子范数的基本性质
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:验证算子范数的范数性质、转置不变性、次乘性和子矩阵单调性。

证明思路

始终回到定义 $\|A\|=\sup_{\|x\|_2=1}\|Ax\|_2$。转置不变性可用最大奇异值或双线性形式。

完整证明

(a) 非负性和齐次性直接来自 Euclidean 范数。若 $\|A\|=0$,则 $Ax=0$ 对所有单位向量成立,故 $A=0$。三角不等式为

$$ \|A+B\| = \sup_{\|x\|=1}\|(A+B)x\| \le \sup_{\|x\|=1}(\|Ax\|+\|Bx\|) \le \|A\|+\|B\|. $$

(b) $\|A\|=s_1(A)$,而 $A$ 与 $A^{\mathsf T}$ 有相同奇异值,所以 $\|A^{\mathsf T}\|=\|A\|$。

(c) 对任意单位向量 $x$,

$$ \|ABx\|\le\|A\|\|Bx\|\le\|A\|\|B\|. $$

取上确界得 $\|AB\|\le\|A\|\|B\|$。例子:$A=\operatorname{diag}(1,0)$,$B=\operatorname{diag}(0,1)$,则 $AB=0$ 但 $\|A\|\|B\|=1$。

(d) 子矩阵可写成 $PAQ$,其中 $P,Q$ 是坐标投影,且 $\|P\|,\|Q\|\le1$。所以 $\|PAQ\|\le\|A\|$。

Exercise 4.3 算子范数:简单例子
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:计算 rank-one 矩阵和对角矩阵的算子范数。

证明思路

rank-one 矩阵只把向量沿 $v$ 方向的分量送到 $u$ 方向;对角矩阵逐坐标缩放。

完整证明

(a) 对任意单位向量 $x$,

$$ uv^{\mathsf T}x=u\langle v,x\rangle, \qquad \|uv^{\mathsf T}x\|\le\|u\|\|v\|. $$

取 $x=v/\|v\|$ 达到等号,所以 $\|uv^{\mathsf T}\|=\|u\|\|v\|$。Frobenius 范数满足

$$ \|uv^{\mathsf T}\|_F^2 = \sum_{ij}u_i^2v_j^2 = \|u\|^2\|v\|^2. $$

(b) 若 $A=\operatorname{diag}(a_i)$,则

$$ \|Ax\|_2^2=\sum_i a_i^2x_i^2\le(\max_i|a_i|)^2\|x\|_2^2. $$

取 $x=e_j$,其中 $|a_j|=\max_i|a_i|$,得到等号。

Exercise 4.4 算子范数与 Frobenius 范数
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明 Frobenius 范数与算子范数、各向同性向量和左乘矩阵之间的基本关系。

证明思路

用奇异值表达 Frobenius 范数;用 trace 表达各向同性向量下的二阶矩;用 $\|Bz\|\le\|B\|\|z\|$ 控制左乘。

完整证明

(a) 若 $\operatorname{rank}(A)=r$,则只有 $r$ 个非零奇异值。因此

$$ \|A\|=s_1\le\left(\sum_i s_i^2\right)^{1/2}=\|A\|_F \le \sqrt r\,s_1=\sqrt r\,\|A\|. $$

rank-one 矩阵使左侧等号成立;$r$ 个相等非零奇异值使右侧等号成立。

(b) 若 $Z$ 各向同性,$\mathbb EZZ^{\mathsf T}=I$,则

$$ \mathbb E\|AZ\|_2^2 = \mathbb E\,Z^{\mathsf T}A^{\mathsf T}AZ = \operatorname{tr}(A^{\mathsf T}A\mathbb EZZ^{\mathsf T}) = \operatorname{tr}(A^{\mathsf T}A) = \|A\|_F^2. $$

(c) 设 $A_1,\ldots,A_n$ 是 $A$ 的列,或直接逐列看 $BA$。有

$$ \|BA\|_F^2 = \sum_j\|B A_j\|_2^2 \le \|B\|^2\sum_j\|A_j\|_2^2 = \|B\|^2\|A\|_F^2. $$

(d) 若 $A$ 对角,则 $BA$ 的第 $j$ 列是 $a_jB_j$,所以

$$ \|BA\|_F^2 = \sum_j a_j^2\|B_j\|_2^2 \le \|B\|_{1\to2}^2\sum_j a_j^2 = \|B\|_{1\to2}^2\|A\|_F^2. $$
Exercise 4.5 奇异值的一个界
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明第 $k$ 个奇异值由 Frobenius 范数控制。

证明思路

Frobenius 范数平方是所有奇异值平方之和,而前 $k$ 个奇异值都至少为 $s_k$。

完整证明

奇异值按降序排列,因此

$$ \|A\|_F^2 = \sum_i s_i(A)^2 \ge \sum_{i=1}^k s_i(A)^2 \ge k\,s_k(A)^2. $$

开平方得到 $s_k(A)\le \|A\|_F/\sqrt k$。

Exercise 4.6 Power method
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明 power method 的极限给出算子范数。

证明思路

在右奇异向量基下展开 Gaussian 初始向量。最高奇异值对应的分量几乎必然非零,反复作用会让最大奇异值主导。

完整证明

设 $A^{\mathsf T}A=V\operatorname{diag}(s_i^2)V^{\mathsf T}$,并写

$$ x=\sum_i\alpha_i v_i. $$

因为 $x$ 是标准 Gaussian,任意固定方向上的系数 $\alpha_i$ 几乎必然非零。于是

$$ (A^{\mathsf T}A)^kx = \sum_i \alpha_i s_i^{2k}v_i. $$

令 $s_1=\|A\|$。若 $s_1>0$,则

$$ \|(A^{\mathsf T}A)^kx\|_2 = s_1^{2k} \left(\sum_i\alpha_i^2(s_i/s_1)^{4k}\right)^{1/2}. $$

括号中的量在 $2k$ 次根下趋于 $1$,因为至少一个 top singular direction 的系数非零。因此

$$ \sqrt[2k]{\|(A^{\mathsf T}A)^kx\|_2}\to s_1=\|A\|. $$

若 $A=0$,结论平凡。

Exercise 4.7 算子范数与列范数
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:用列范数给出算子范数下界,并刻画等号情形。

证明思路

列向量就是 $Ae_i$。若某个列向量已经达到算子范数,则 $e_i$ 是 $A^{\mathsf T}A$ 的最大特征向量,迫使该列与其余列正交。

完整证明

(a) 对每个 $i$,

$$ \|A\|\ge\|Ae_i\|_2=\|A_i\|_2. $$

若列两两正交,则

$$ \|Ax\|_2^2 = \left\|\sum_i x_iA_i\right\|_2^2 = \sum_i x_i^2\|A_i\|_2^2 \le \max_i\|A_i\|_2^2\|x\|_2^2, $$

所以等号成立。

(b) 若 $\|A\|=\|A_i\|_2$,则 $e_i$ 使 Rayleigh quotient $x^{\mathsf T}A^{\mathsf T}Ax$ 取到最大值。因此 $e_i$ 是 $A^{\mathsf T}A$ 的特征向量。于是对 $j\ne i$,

$$ \langle A_i,A_j\rangle = e_i^{\mathsf T}A^{\mathsf T}Ae_j = 0. $$

(c) 对行应用同样结论于 $A^{\mathsf T}$。

Exercise 4.8 Schur bound
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明 Schur bound:算子范数由最大行绝对和与最大列绝对和控制。

证明思路

对每一行使用 Cauchy-Schwarz 的加权版本,再交换求和,把行和、列和分别提出来。

完整证明

$$ R=\max_i\sum_j|a_{ij}|,\qquad C=\max_j\sum_i|a_{ij}|. $$

对任意 $x$,

$$ |(Ax)_i|^2 = \left|\sum_j a_{ij}x_j\right|^2 \le \left(\sum_j|a_{ij}|\right) \left(\sum_j|a_{ij}||x_j|^2\right) \le R\sum_j|a_{ij}||x_j|^2. $$

对 $i$ 求和得到

$$ \|Ax\|_2^2 \le R\sum_j\left(\sum_i|a_{ij}|\right)|x_j|^2 \le RC\|x\|_2^2. $$

取单位向量上的上确界,$\|A\|\le\sqrt{RC}$。

Exercise 4.9 Walsh 矩阵
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明归一化 Walsh 矩阵是正交矩阵。

证明思路

用归纳法检查 $W_kW_k^{\mathsf T}=nI_n$,其中 $n$ 是矩阵阶数。

完整证明

当 $k=1$ 时,

$$ W_1W_1^{\mathsf T}=2I_2. $$

假设 $W_kW_k^{\mathsf T}=nI_n$。则

$$ W_{k+1}W_{k+1}^{\mathsf T} = \begin{bmatrix} W_k&W_k\\ W_k&-W_k \end{bmatrix} \begin{bmatrix} W_k^{\mathsf T}&W_k^{\mathsf T}\\ W_k^{\mathsf T}&-W_k^{\mathsf T} \end{bmatrix} = \begin{bmatrix} 2nI_n&0\\ 0&2nI_n \end{bmatrix}. $$

所以若 $W_{k+1}$ 阶数为 $2n$,则 $W_{k+1}W_{k+1}^{\mathsf T}=(2n)I_{2n}$。归纳完成,故 $W_k/\sqrt n$ 正交。

Exercise 4.10 具有 $\pm1$ 元素的矩阵
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明 $\pm1$ 矩阵算子范数介于 $\sqrt n$ 与 $n$ 之间,并说明两端可达。

证明思路

上界由 Frobenius 范数给出;下界由 Frobenius 范数和秩给出。全 $1$ 矩阵达到上界,Walsh 矩阵达到下界。

完整证明

因为 $A$ 有 $n^2$ 个元素且每个绝对值为 $1$,

$$ \|A\|_F=n. $$

上界由 $\|A\|\le\|A\|_F$ 得到:$\|A\|\le n$。下界由 Exercise 4.4(a) 得到:

$$ \|A\|_F\le\sqrt{\operatorname{rank}(A)}\,\|A\| \le \sqrt n\,\|A\|, $$

所以 $\|A\|\ge\sqrt n$。

全 $1$ 矩阵秩为 $1$,唯一非零奇异值为 $n$,达到上界。对 $n=2^k$,Walsh 矩阵满足 $WW^{\mathsf T}=nI_n$,所有奇异值都是 $\sqrt n$,因此达到下界。

Exercise 4.11 正交投影之差
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明任意两个正交投影之差的算子范数不超过 $1$。

证明思路

$P-Q$ 是自伴矩阵。对自伴矩阵,算子范数由 Rayleigh quotient 控制;而 $\langle(P-Q)x,x\rangle=\|Px\|^2-\|Qx\|^2$ 位于 $[-1,1]$。

完整证明

因为 $P,Q$ 是正交投影,$P=P^{\mathsf T}=P^2$,$Q=Q^{\mathsf T}=Q^2$,所以 $P-Q$ 自伴。对任意单位向量 $x$,

$$ \langle(P-Q)x,x\rangle = \langle Px,x\rangle-\langle Qx,x\rangle = \|Px\|_2^2-\|Qx\|_2^2. $$

两个平方范数都在 $[0,1]$ 中,因此该差属于 $[-1,1]$。由自伴矩阵的变分刻画,

$$ \|P-Q\| = \sup_{\|x\|_2=1}|\langle(P-Q)x,x\rangle| \le1. $$
Exercise 4.12 正交投影之差与乘积
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:用互补投影刻画两个正交投影之差的范数。

证明思路

先把 $P-Q$ 写成两个正交分量之和,再利用秩相等时 $PQ$ 与 $QP$ 的非零奇异值相同。

完整证明

先证明 (a)。恒等式

$$ P-Q=PQ_\perp-P_\perp Q $$

来自 $Q_\perp=I-Q$ 和 $P_\perp=I-P$。对任意单位向量 $x$,两个向量 $PQ_\perp x$ 与 $P_\perp Qx$ 分别落在 $\operatorname{range}(P)$ 与 $\operatorname{range}(P)^\perp$,所以正交。因此

$$ \|(P-Q)x\|_2^2 = \|PQ_\perp x\|_2^2+\|P_\perp Qx\|_2^2 \le M^2\bigl(\|Q_\perp x\|_2^2+\|Qx\|_2^2\bigr) = M^2, $$

其中 $M=\max(\|PQ_\perp\|,\|P_\perp Q\|)$。这给出 $\|P-Q\|\le M$。反向不等式也直接成立,因为

$$ PQ_\perp=(P-Q)Q_\perp, \qquad P_\perp Q=-P_\perp(P-Q), $$

而投影的算子范数不超过 $1$。故

$$ \|P-Q\| = \max(\|P_\perp Q\|,\|PQ_\perp\|). $$

下面证明 (b)。若 $\operatorname{rank}(P)>\operatorname{rank}(Q)$,则

$$ \dim\operatorname{range}(P)+\dim\ker(Q)>n, $$

所以存在单位向量 $x\in\operatorname{range}(P)\cap\ker(Q)$。于是 $(P-Q)x=x$,故 $\|P-Q\|\ge1$。结合 Exercise 4.11 的 $\|P-Q\|\le1$,得到 $\|P-Q\|=1$。若 $\operatorname{rank}(Q)>\operatorname{rank}(P)$,同理。

最后证明 (c)。令两个秩都为 $k$。有

$$ \|P_\perp Q\|^2 = \|QP_\perp Q\|, \qquad \|PQ_\perp\|^2 = \|P Q_\perp P\|. $$

在 $\operatorname{range}(Q)$ 上,$QPQ$ 的非零特征值等于算子 $PQ$ 的奇异值平方;在 $\operatorname{range}(P)$ 上,$PQP$ 的非零特征值也是同一批数。由于两个空间维数都为 $k$,零特征值的重数也相同。因此 $QP_\perp Q=Q-QPQ$ 与 $PQ_\perp P=P-PQP$ 有相同特征值,从而两个范数相等。结合 (a),得到

$$ \|P-Q\|=\|P_\perp Q\|=\|PQ_\perp\|. $$
Exercise 4.13 谱投影的 Davis-Kahan
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明 top $k$ 谱投影的 Davis-Kahan 型界。

证明思路

把 top $k$ 特征值看成一个谱簇。只要该谱簇与其余谱相隔 $\lambda_k-\lambda_{k+1}$,Davis-Kahan 定理直接控制两个谱投影。

完整证明

设 $\delta=\lambda_k-\lambda_{k+1}>0$。若 $2\|A-B\|\ge\delta$,右侧至少为 $1$,而 $\|P_A-P_B\|\le1$,结论平凡。

否则令 $\varepsilon=\|A-B\|<\delta/2$。由 Weyl 不等式,对 $j>k$ 有

$$ \mu_j\le\lambda_j+\varepsilon\le\lambda_{k+1}+\varepsilon. $$

而 $A$ 的 top $k$ 谱簇位于 $[\lambda_k,\infty)$ 中,所以该谱簇与 $B$ 的 bottom 谱簇 $\{\mu_j:j>k\}$ 的距离至少

$$ \lambda_k-(\lambda_{k+1}+\varepsilon) = \delta-\varepsilon \ge \delta/2. $$

把 Theorem 4.1.15 证明中的谱投影分离引理用于 $A$ 的 top $k$ 投影 $P_A$ 和 $B$ 的 bottom 投影 $I-P_B$,得到

$$ \|(I-P_B)P_A\| \le \frac{\varepsilon}{\delta/2} = \frac{2\|A-B\|}{\delta}. $$

由于 $P_A$ 与 $P_B$ 秩相同,Exercise 4.12(c) 给出 $\|P_A-P_B\|=\|(I-P_B)P_A\|$。因此

$$ \|P_A-P_B\| \le \frac{2\|A-B\|}{\lambda_k-\lambda_{k+1}}. $$
Exercise 4.14 Hermitian dilation
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明 Hermitian dilation 的非零特征值是原矩阵奇异值的正负号版本。

证明思路

把 SVD 中的左右奇异向量拼成分块向量,直接作用分块矩阵。

完整证明

若 $Av_i=s_iu_i$ 且 $A^{\mathsf T}u_i=s_iv_i$,则

$$ H\begin{bmatrix}u_i\\ v_i\end{bmatrix} = \begin{bmatrix}Av_i\\ A^{\mathsf T}u_i\end{bmatrix} = s_i\begin{bmatrix}u_i\\ v_i\end{bmatrix}, $$

并且

$$ H\begin{bmatrix}u_i\\ -v_i\end{bmatrix} = \begin{bmatrix}-Av_i\\ A^{\mathsf T}u_i\end{bmatrix} = -s_i\begin{bmatrix}u_i\\ -v_i\end{bmatrix}. $$

因此每个非零奇异值 $s_i$ 产生两个特征值 $\pm s_i$。剩余维度对应 $A$ 或 $A^{\mathsf T}$ 的核,特征值为 $0$。

Exercise 4.15 Wedin theorem
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:把矩形矩阵的奇异子空间扰动界化为 Hermitian dilation 的 Davis-Kahan 定理。

证明思路

构造 $H_A$ 和 $H_B$。它们的谱投影差由 Davis-Kahan 控制;左奇异子空间投影只是该分块谱投影的左上角压缩,所以范数不会增大。

完整证明

$$ H_A=\begin{bmatrix}0&A\\ A^{\mathsf T}&0\end{bmatrix}, \qquad H_B=\begin{bmatrix}0&B\\ B^{\mathsf T}&0\end{bmatrix}. $$

由 Exercise 4.14,$H_A$ 的正特征值是 $s_i(A)$,对应特征向量为 $[u_i;v_i]$ 的归一化版本。并且

$$ \|H_A-H_B\|=\|A-B\|. $$

令 $R_A$ 为 $H_A$ 对谱簇 $\{\pm s_1,\ldots,\pm s_k\}$ 的谱投影,并令 $R_B$ 为 $H_B$ 对对应谱簇的谱投影。这个投影同时包含正、负两组特征向量;对每个 $i\le k$,两个归一化特征向量 $[u_i;v_i]/\sqrt2$ 与 $[u_i;-v_i]/\sqrt2$ 的投影相加为

$$ \begin{bmatrix} u_i u_i^{\mathsf T} & 0\\ 0 & v_i v_i^{\mathsf T} \end{bmatrix}. $$

所以 $R_A$ 的左上角分块正是 $P_A$,右下角分块正是 top $k$ 右奇异向量投影。

若 $2\|A-B\|\ge s_k-s_{k+1}$,则右侧至少为 $1$,而两个正交投影之差的范数不超过 $1$,结论成立。否则 Weyl 不等式保证 $H_B$ 中对应的 $\pm t_1,\ldots,\pm t_k$ 谱簇仍与其余谱相隔至少 $(s_k-s_{k+1})/2$。对这个分离谱簇应用投影型 Davis-Kahan,得到

$$ \|R_A-R_B\| \le \frac{2\|A-B\|}{s_k-s_{k+1}}. $$

取左上角分块不会增加算子范数,于是 top $k$ 左奇异向量投影满足题中界。取右下角分块得到 top $k$ 右奇异向量投影的同样界。

单个第 $k$ 个奇异向量的版本作用于只含 $\pm s_k$ 的谱簇。若 $s_k$ 与相邻奇异值的间隔为

$$ \Delta_k=\min(s_{k-1}-s_k,\ s_k-s_{k+1}) $$

(端点处只保留存在的项),则投影误差至多 $2\|A-B\|/\Delta_k$;再用 Exercise 4.16 可把一维投影距离转为奇异向量在差一个全局符号意义下的 Euclidean 距离界。

Exercise 4.16 向量之间的夹角与距离
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:把小夹角转换成向量本身在差符号意义下的 Euclidean 距离。

证明思路

选取符号让内积非负。距离平方等于 $2-2|\langle u,v\rangle|$,再用 $\sin$ 控制。

完整证明

取 $\theta\in\{-1,1\}$ 使 $\langle u,\theta v\rangle=|\langle u,v\rangle|=\cos\angle(u,v)\ge0$。则

$$ \|u-\theta v\|_2^2 = 2-2\cos\angle(u,v). $$

对 $c\in[0,1]$ 有 $1-c\le1-c^2$。因此

$$ \|u-\theta v\|_2^2 \le 2(1-\cos^2\angle(u,v)) = 2\sin^2\angle(u,v) \le 2\varepsilon^2. $$

开平方得到结论。

Exercise 4.17 近似投影
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:把 isometry / approximate isometry 与 $AA^{\mathsf T}$ 接近投影联系起来。

证明思路

用 SVD。$A^{\mathsf T}A$ 和 $AA^{\mathsf T}$ 共享非零特征值;这些特征值就是奇异值平方。

完整证明

(a) 若 $A$ 是 isometry,则 $A^{\mathsf T}A=I_n$。矩阵 $AA^{\mathsf T}$ 对称,且

$$ (AA^{\mathsf T})^2=A(A^{\mathsf T}A)A^{\mathsf T}=AA^{\mathsf T}, $$

所以它是到 $\operatorname{range}(A)$ 的正交投影。反过来,若 $AA^{\mathsf T}$ 是秩为 $n$ 的正交投影,则其非零特征值全为 $1$。这些正是 $A^{\mathsf T}A$ 的特征值,所以 $A^{\mathsf T}A=I_n$,即 $A$ 是 isometry。

(b) 写完整的 $n$ 个右奇异方向版本:$A=U\Sigma V^{\mathsf T}$,其中 $U$ 是 $m\times n$ 的列正交矩阵;若 $A$ 有零奇异值,就把对应左奇异向量任意补成正交列。令 $P=UU^{\mathsf T}$,这是秩为 $n$ 的正交投影。则

$$ AA^{\mathsf T}=U\Sigma^2U^{\mathsf T}, \qquad AA^{\mathsf T}-P=U(\Sigma^2-I)U^{\mathsf T}. $$

因此

$$ \|AA^{\mathsf T}-P\| = \|\Sigma^2-I\| = \|A^{\mathsf T}A-I\|. $$

这正是 approximate isometry 的二次型表述,故两者等价。

Exercise 4.18 $p\to q$ 范数
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:建立 $p\to q$ operator norm 的双线性表达、范数性质和对偶关系。

证明思路

用 $\ell^q$ 与 $\ell^{q'}$ 的对偶性:$\|z\|_q=\max_{\|y\|_{q'}=1}|y^{\mathsf T}z|$。

完整证明

(a) 对任意 $x$,

$$ \|Ax\|_q = \max_{\|y\|_{q'}=1}|y^{\mathsf T}Ax|. $$

再对 $\|x\|_p=1$ 最大化,得到 (4.31)。

(b) 非负性、齐次性和三角不等式都由 $\|Ax\|_q/\|x\|_p$ 的定义直接推出;若范数为 $0$,则 $Ax=0$ 对所有 $x$ 成立,所以 $A=0$。

(c) 由双线性表达,

$$ \|A^{\mathsf T}\|_{p\to q} = \max_{\|x\|_p=1,\ \|y\|_{q'}=1} |y^{\mathsf T}A^{\mathsf T}x| = \max_{\|y\|_{q'}=1,\ \|x\|_p=1} |x^{\mathsf T}Ay| = \|A\|_{q'\to p'}. $$
Exercise 4.19 $1\to\infty$ 和 $1\to2$ 范数
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:给出 $1\to\infty$、$1\to2$、$2\to\infty$ 范数的元素级表达。

证明思路

$\ell^1$ 单位球的极点是 $\pm e_j$;$\ell^\infty$ 输出看行,$\ell^2$ 输出看列。

完整证明

(a) 对任意 $\|x\|_1\le1$,

$$ |(Ax)_i| \le \sum_j |A_{ij}||x_j| \le \max_j|A_{ij}|. $$

再对 $i$ 取最大,得 $\|A\|_{1\to\infty}\le\max_{ij}|A_{ij}|$。取 $x=\pm e_j$ 达到最大元素,故等号成立。

(b) 设 $A_j$ 为第 $j$ 列。若 $\|x\|_1\le1$,则

$$ \|Ax\|_2 = \left\|\sum_jx_jA_j\right\|_2 \le \sum_j|x_j|\|A_j\|_2 \le \max_j\|A_j\|_2. $$

取 $x=e_j$ 得等号。因此 $\|A\|_{1\to2}=\max_j\|A_{:j}\|_2$。对 $2\to\infty$,第 $i$ 个坐标是行向量 $A_{i:}$ 与 $x$ 的内积,Cauchy-Schwarz 给出上界,取 $x$ 平行于最大行得到等号。

Exercise 4.20 $\infty\to1$ 范数
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明 $\infty\to1$ 范数可在符号向量上取到,并给出 rank-one 对偶表述和二次型反例。

证明思路

线性函数在盒子上最大值取在顶点。rank-one 表述把 $y^{\mathsf T}Ax$ 写成矩阵内积。二次型反例来自非凸二次函数在盒子内部达到更大绝对值。

完整证明

(a) 由 $\ell^1$ 与 $\ell^\infty$ 对偶,

$$ \|A\|_{\infty\to1} = \max_{\|x\|_\infty\le1,\ \|y\|_\infty\le1}|y^{\mathsf T}Ax|. $$

该双线性函数在两个盒子上最大,固定一边时另一边的最大值都在顶点取得,所以可取 $x,y\in\{-1,1\}$。

(b) 对秩一矩阵 $Z=yx^{\mathsf T}$,有 $\langle A,Z\rangle=y^{\mathsf T}Ax$。若 $\|x\|_\infty,\|y\|_\infty\le1$,则 $|Z_{ij}|\le1$。反过来,任意秩一 $Z=uv^{\mathsf T}$ 且 $|Z_{ij}|\le1$,可重新缩放为 $Z=yx^{\mathsf T}$,其中 $\|x\|_\infty,\|y\|_\infty\le1$。于是得到所需上确界。

(c) 取

$$ A= \begin{bmatrix} -1&0\\ 0&2 \end{bmatrix}. $$

在顶点 $x\in\{-1,1\}^2$ 上,$x^{\mathsf T}Ax=1$,所以最大绝对值为 $1$。但在 $x=(0,1)$ 处,$x^{\mathsf T}Ax=2$。因此二次型版本不能只检查符号顶点。

Exercise 4.21 Cut norm
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明 cut norm 与 $\infty\to1$ 范数在常数因子内等价。

证明思路

cut norm 是用 $0/1$ 指示向量选择子矩阵;$\infty\to1$ 范数是用 $\pm1$ 符号向量。两者可通过把符号向量分解为正负集合互相转换。

完整证明

对任意 $I,J$,令 $y=\mathbf 1_I$、$x=\mathbf 1_J$,它们的 $\ell^\infty$ 范数不超过 $1$。因此

$$ \left|\sum_{i\in I,j\in J}A_{ij}\right| = |y^{\mathsf T}Ax| \le \|A\|_{\infty\to1}. $$

所以 $\|A\|_{\mathrm{cut}}\le\|A\|_{\infty\to1}$。

反过来,由 Exercise 4.20,只需考虑 $x\in\{-1,1\}^n$、$y\in\{-1,1\}^m$。令 $I_\pm=\{i:y_i=\pm1\}$,$J_\pm=\{j:x_j=\pm1\}$。则 $y^{\mathsf T}Ax$ 是四个矩形和的带符号组合:

$$ S(I_+,J_+)-S(I_+,J_-)-S(I_-,J_+)+S(I_-,J_-). $$

每个矩形和的绝对值都不超过 cut norm,因此

$$ |y^{\mathsf T}Ax|\le4\|A\|_{\mathrm{cut}}. $$

对符号向量取最大,即得 $\|A\|_{\infty\to1}\le4\|A\|_{\mathrm{cut}}$。

Exercise 4.22 SDP 松弛
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:说明若干难算矩阵范数可由 SDP 在常数因子内近似。

证明思路

核心是把符号或盒约束优化松弛为向量 Gram 矩阵优化。Grothendieck 不等式保证 $\infty\to1$ 的 SDP 松弛只损失常数;其他范数通过等价关系或平方化归约到同类 SDP。

完整证明

对 $\|A\|_{\infty\to1}$,Exercise 4.20 把它写成符号双线性优化。把符号 $x_j,y_i$ 放松为单位向量 $u_j,v_i$,目标变为

$$ \max \sum_{ij}A_{ij}\langle v_i,u_j\rangle, $$

这是 Gram 矩阵变量上的 SDP。Grothendieck 不等式说明这个 SDP 值与原值只差绝对常数因子。

cut norm 与 $\infty\to1$ 范数由 Exercise 4.21 在常数因子内等价,因此也可常数近似。

对于 $\|A\|_{\infty\to2}$,有

$$ \|A\|_{\infty\to2}^2 = \max_{x\in\{-1,1\}^n}x^{\mathsf T}A^{\mathsf T}Ax $$

在常数因子意义下可由正半定二次 Grothendieck/SDP 松弛近似。最后 $\|A\|_{2\to1}=\|A^{\mathsf T}\|_{\infty\to2}$,所以同样可由 SDP 常数近似。

Exercise 4.23 $\varepsilon$-net 的传递性
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明 net 的覆盖误差满足三角不等式式的传递。

证明思路

对 $K$ 中任意点,先用 $\mathcal M$ 近似,再用 $\mathcal N$ 近似 $\mathcal M$ 中的点。两次误差相加。

完整证明

任取 $x\in K$。因为 $\mathcal M$ 是 $K$ 的 $\delta$-net,存在 $m\in\mathcal M$ 使

$$ d(x,m)\le\delta. $$

又因为 $\mathcal N$ 是 $\mathcal M$ 的 $\varepsilon$-net,存在 $y\in\mathcal N$ 使

$$ d(m,y)\le\varepsilon. $$

由三角不等式,

$$ d(x,y)\le d(x,m)+d(m,y)\le\delta+\varepsilon. $$

这对任意 $x\in K$ 成立,所以 $\mathcal N$ 是 $K$ 的 $(\varepsilon+\delta)$-net。

Exercise 4.24 Packing balls
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:说明 separated set 与半径 $\varepsilon/2$ balls 两两不交的逆命题在一般 metric space 中失败,但在 normed space 中成立。

证明思路

一般 metric space 可能没有两个点之间的“中点”。normed space 有中点;若两中心距离不超过 $\varepsilon$,中点会同时落入两个闭球。

完整证明

(a) 取 $T=\{a,b\}$,定义 $d(a,b)=0.75\varepsilon$。以 $a,b$ 为中心的闭 $\varepsilon/2$-balls 分别为 $\{a\}$ 和 $\{b\}$,因为 $0.75\varepsilon>\varepsilon/2$,所以它们不交。但 $d(a,b)<\varepsilon$,因此 $\{a,b\}$ 不是 $\varepsilon$-separated。这给出一般 metric space 中的反例。

(b) 设 $T$ 是 normed space 的子集,且以 $\mathcal N$ 中点为中心的闭 $\varepsilon/2$-balls 两两不交。若存在不同 $x,y\in\mathcal N$ 满足 $\|x-y\|\le\varepsilon$,则中点 $z=(x+y)/2$ 满足

$$ \|z-x\|=\|z-y\|=\frac12\|x-y\|\le\frac\varepsilon2. $$

于是 $z$ 同时落入两个闭 $\varepsilon/2$-balls,矛盾。因此任意不同中心距离都大于 $\varepsilon$,即 $\mathcal N$ 是 $\varepsilon$-separated。

Exercise 4.25 External covering numbers
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明允许球心在 $K$ 外部不会本质改变 covering number,只会改变半径常数。

证明思路

内部 covering 是 external covering 的特殊情形,所以外部覆盖数不比内部覆盖数大。反向时,把每个外部中心移动到它所覆盖到的一个 $K$ 中点上,半径从 $\varepsilon/2$ 增大到 $\varepsilon$。

完整证明

第一部分直接成立:普通 covering 要求中心在 $K$ 中,是 external covering 的特殊情形,所以

$$ \overline{\mathcal N}(K,d,\varepsilon)\le \mathcal N(K,d,\varepsilon). $$

下面证明第二个不等式。设 $y_1,\ldots,y_M$ 是 $K$ 的 external $\varepsilon/2$-net,其中

$$ M=\overline{\mathcal N}(K,d,\varepsilon/2). $$

丢掉那些半径 $\varepsilon/2$ 球与 $K$ 没有交集的中心。对每个剩余的 $y_i$,选取一个 $x_i\in K$ 满足 $d(x_i,y_i)\le\varepsilon/2$。若 $z\in K$ 被 $y_i$ 覆盖,即 $d(z,y_i)\le\varepsilon/2$,则

$$ d(z,x_i)\le d(z,y_i)+d(y_i,x_i)\le\varepsilon. $$

所以这些 $x_i\in K$ 构成 $K$ 的普通 $\varepsilon$-net,且数量不超过 $M$。因此

$$ \mathcal N(K,d,\varepsilon)\le\overline{\mathcal N}(K,d,\varepsilon/2). $$
Exercise 4.26 Covering number 的单调性
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:说明 covering number 关于半径单调,但关于集合包含只满足带半径损失的近似单调性。

证明思路

反例来自“内部球心”的限制:大集合可能包含一个中心点,可以同时覆盖小集合的两端;小集合自身没有这个中心点。近似单调性用 Exercise 4.25 的移动中心技巧。

完整证明

关于 $\varepsilon$ 的单调性是直接的:半径越大,每个球覆盖越多,所以 $\varepsilon_1\le\varepsilon_2$ 时

$$ \mathcal N(K,d,\varepsilon_2)\le \mathcal N(K,d,\varepsilon_1). $$

(a) 取实线中的 $K=\{-1,0,1\}$,$L=\{-1,1\}$,并令 $\varepsilon=1.1$。集合 $K$ 可由一个以 $0\in K$ 为中心的 $\varepsilon$-ball 覆盖,所以 $\mathcal N(K,d,\varepsilon)=1$。但 $L$ 的中心必须在 $L$ 中;以 $-1$ 或 $1$ 为中心的半径 $1.1$ 球都不能覆盖另一个点,因为两点距离为 $2$。因此 $\mathcal N(L,d,\varepsilon)=2$,这与普通单调性相反。

(b) 设 $L\subset K$。取 $K$ 的一个 $\varepsilon/2$-net,中心为 $x_1,\ldots,x_M\in K$,其中 $M=\mathcal N(K,d,\varepsilon/2)$。只保留那些 $\varepsilon/2$-球与 $L$ 相交的中心;对每个保留中心 $x_i$,选一个 $y_i\in L$ 满足 $d(x_i,y_i)\le\varepsilon/2$。若 $z\in L$ 被 $x_i$ 覆盖,则

$$ d(z,y_i)\le d(z,x_i)+d(x_i,y_i)\le\varepsilon. $$

于是这些 $y_i\in L$ 构成 $L$ 的 $\varepsilon$-net,数量不超过 $M$。所以

$$ \mathcal N(L,d,\varepsilon)\le\mathcal N(K,d,\varepsilon/2). $$
Exercise 4.27 Euclidean ball 的体积:几何论证
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:用 simplex、$\ell^1$ ball 和 cube 的几何关系给出 Euclidean ball 体积的粗界。

证明思路

先算 canonical simplex 的体积,再把 $\ell^1$ ball 分成 $2^n$ 个 simplex。Euclidean ball 夹在一个 cube 和缩放的 $\ell^1$ ball 之间。

完整证明

(a) 对 $\Delta_n=\{x_i\ge0,\ \sum_i x_i\le1\}$,按最后一个坐标积分:

$$ \operatorname{Vol}(\Delta_n) = \int_0^1 \operatorname{Vol}\{x_1,\ldots,x_{n-1}\ge0,\ \sum_{i=1}^{n-1}x_i\le1-t\}\,dt. $$

由归纳假设,该截面体积为 $(1-t)^{n-1}/(n-1)!$。因此

$$ \operatorname{Vol}(\Delta_n) = \frac1{(n-1)!}\int_0^1(1-t)^{n-1}\,dt = \frac1{n!}. $$

(b) $B_1^n$ 被 $2^n$ 个正交象限分割;每个象限中的部分都是一个 canonical simplex。因此

$$ \operatorname{Vol}(B_1^n)=\frac{2^n}{n!}. $$

由 $n!\ge(n/e)^n$,得到 $\operatorname{Vol}(B_1^n)\le(2e/n)^n$。

(c) 若 $\|x\|_\infty\le1/\sqrt n$,则 $\|x\|_2\le1$,所以

$$ [-1/\sqrt n,1/\sqrt n]^n\subset B_2^n. $$

这给出下界 $\operatorname{Vol}(B_2^n)\ge(2/\sqrt n)^n$。另一方面,由 Cauchy-Schwarz,$\|x\|_1\le\sqrt n\|x\|_2$,所以

$$ B_2^n\subset \sqrt n\,B_1^n. $$

因此

$$ \operatorname{Vol}(B_2^n) \le n^{n/2}\operatorname{Vol}(B_1^n) \le n^{n/2}\left(\frac{2e}{n}\right)^n = \left(\frac{2e}{\sqrt n}\right)^n. $$
Exercise 4.28 Euclidean ball 的体积:概率论证
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:用 Gaussian 概率给出单位 Euclidean ball 体积的渐近尖锐上界。

证明思路

把 Gaussian 落入半径 $\sqrt n$ 球的概率用密度下界控制。概率不超过 $1$,于是反推出球体积上界。

完整证明

设 $g\sim N(0,I_n)$。在球 $\sqrt n B_2^n$ 上,Gaussian 密度满足

$$ (2\pi)^{-n/2}e^{-\|x\|_2^2/2} \ge (2\pi)^{-n/2}e^{-n/2}. $$

因此

$$ 1 \ge \mathbb P\{g\in\sqrt n B_2^n\} \ge (2\pi)^{-n/2}e^{-n/2}\operatorname{Vol}(\sqrt n B_2^n). $$

而 $\operatorname{Vol}(\sqrt n B_2^n)=n^{n/2}\operatorname{Vol}(B_2^n)$,所以

$$ \operatorname{Vol}(B_2^n) \le \left(\sqrt{\frac{2\pi e}{n}}\right)^n. $$
Exercise 4.29 Euclidean ball 的体积:解析论证
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明一般范数的径向积分公式,并由此精确计算 Euclidean ball 体积和单位体积半径渐近。

证明思路

范数球满足 $\operatorname{Vol}(tB)=t^n\operatorname{Vol}(B)$。把积分按层集写成 Stieltjes 积分,再分部积分。代入 Gaussian kernel 后使用 Gamma 函数和 Stirling 公式。

完整证明

(a) 因为 $\{x:\|x\|\le t\}=tB$,其体积为 $t^n\operatorname{Vol}(B)$。对递减可微且趋于 $0$ 的 $f$,层集积分给出

$$ \int_{\mathbb R^n}f(\|x\|)\,dx = \int_0^\infty f(t)\,d\bigl(t^n\operatorname{Vol}(B)\bigr). $$

分部积分,边界项为零,于是

$$ \int_{\mathbb R^n}f(\|x\|)\,dx = -\operatorname{Vol}(B)\int_0^\infty t^n f'(t)\,dt. $$

(b) 对 Euclidean norm 取 $f(t)=e^{-t^2/2}$。左侧为

$$ \int_{\mathbb R^n}e^{-\|x\|_2^2/2}\,dx=(2\pi)^{n/2}. $$

右侧中的积分为

$$ -\int_0^\infty t^n f'(t)\,dt = \int_0^\infty t^{n+1}e^{-t^2/2}\,dt = 2^{n/2}\Gamma(n/2+1). $$

所以

$$ \operatorname{Vol}(B_2^n) = \frac{(2\pi)^{n/2}}{2^{n/2}\Gamma(n/2+1)} = \frac{\pi^{n/2}}{\Gamma(n/2+1)}. $$

(c) 体积为 $1$ 的 Euclidean ball 半径 $R_n$ 满足

$$ R_n^n\operatorname{Vol}(B_2^n)=1. $$

由 Stirling 公式 $\Gamma(n/2+1)=\sqrt{\pi n}(n/(2e))^{n/2}(1+o(1))$,得到

$$ \operatorname{Vol}(B_2^n)^{1/n} = \sqrt{\frac{2\pi e}{n}}\,(1+o(1)). $$

因此

$$ R_n=\operatorname{Vol}(B_2^n)^{-1/n} = \sqrt{\frac{n}{2\pi e}}\,(1+o(1)). $$
Exercise 4.30 $\ell^p$ ball 的体积
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:给出 $\ell^p$ 单位球体积的粗界和 Gamma 函数精确公式。

证明思路

粗界来自 cube 与缩放 $\ell^1$ ball 的夹逼。精确公式沿用 Exercise 4.29 的径向积分,取 $f(t)=e^{-t^p}$。

完整证明

(a) 若 $\|x\|_\infty\le n^{-1/p}$,则 $\|x\|_p\le1$,所以

$$ [-n^{-1/p},n^{-1/p}]^n\subset B_p^n. $$

这给出下界 $(2/n^{1/p})^n$。另一方面,Hölder 不等式给出 $\|x\|_1\le n^{1-1/p}\|x\|_p$,因此

$$ B_p^n\subset n^{1-1/p}B_1^n. $$

结合 $\operatorname{Vol}(B_1^n)\le(2e/n)^n$,得到

$$ \operatorname{Vol}(B_p^n) \le n^{n(1-1/p)}\left(\frac{2e}{n}\right)^n = \left(\frac{2e}{n^{1/p}}\right)^n. $$

(b) 在 Exercise 4.29 的公式中取 $\|\cdot\|=\|\cdot\|_p$ 和 $f(t)=e^{-t^p}$。左侧分解为坐标积分:

$$ \int_{\mathbb R^n}e^{-\|x\|_p^p}\,dx = \left(2\int_0^\infty e^{-s^p}\,ds\right)^n = \bigl(2\Gamma(1+1/p)\bigr)^n. $$

右侧的径向积分为

$$ -\int_0^\infty t^n f'(t)\,dt = p\int_0^\infty t^{n+p-1}e^{-t^p}\,dt = \Gamma(n/p+1). $$

因此

$$ \operatorname{Vol}(B_p^n) = \frac{\bigl(2\Gamma(1/p+1)\bigr)^n}{\Gamma(n/p+1)}. $$
Exercise 4.31 Lattice 是 $\varepsilon$-net
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明缩放整数 lattice 在单位 Euclidean ball 中给出显式 $\varepsilon$-net,并估计其大小。

证明思路

对每个坐标向零取整,保证格点仍在单位球内,同时每个坐标误差不超过 $\varepsilon/\sqrt n$。大小估计用格点小立方体和 Euclidean ball 的体积上界。

完整证明

令格距 $h=\varepsilon/\sqrt n$。给定 $x\in B_2^n$,定义 $y_i$ 为与 $x_i$ 同号且满足

$$ |y_i|=h\left\lfloor\frac{|x_i|}{h}\right\rfloor. $$

则 $y\in h\mathbb Z^n$,并且 $|y_i|\le |x_i|$,所以 $\|y\|_2\le\|x\|_2\le1$,即 $y\in\mathcal N$。同时 $|x_i-y_i|\le h$,故

$$ \|x-y\|_2\le \sqrt n\,h=\varepsilon. $$

因此 $\mathcal N$ 是 $B_2^n$ 的 $\varepsilon$-net。这个构造也给出了快速近似算法:逐坐标向零取整即可。

下面估计大小。以每个 $y\in\mathcal N$ 为角放置半开小立方体 $y+[0,h)^n$;这些小立方体互不相交,并且每个点到 $y$ 的 Euclidean 距离至多 $\sqrt n h=\varepsilon$。因此它们都包含在 $(1+\varepsilon)B_2^n$ 中。于是

$$ |\mathcal N|h^n \le \operatorname{Vol}((1+\varepsilon)B_2^n). $$

用 Exercise 4.28 的体积上界,

$$ |\mathcal N| \le \left(\frac{\sqrt n}{\varepsilon}\right)^n (1+\varepsilon)^n \left(\sqrt{\frac{2\pi e}{n}}\right)^n = \left(\frac{\sqrt{2\pi e}(1+\varepsilon)}{\varepsilon}\right)^n. $$

对通常的 net 范围 $0<\varepsilon\le1$,右侧被 $e^n(2/\varepsilon+1)^n$ 控制;若 $\varepsilon>1$,该界更平凡成立。故得到题中估计。

Exercise 4.32 Hamming cube 的 covering 与 packing numbers
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明 Proposition 4.2.15 中 Hamming cube 的 covering 与 packing number 界。

证明思路

把 Euclidean 体积换成 Hamming ball 的基数。半径 $m$ 的 Hamming ball 大小为 $\sum_{k=0}^m\binom nk$。

完整证明

在 $\{0,1\}^n$ 中,任意中心的半径 $m$ Hamming ball 的基数为

$$ V_m=\sum_{k=0}^m\binom nk, $$

因为可以选择至多 $m$ 个坐标翻转。

若 $M$ 个半径 $m$ 球覆盖整个 cube,则 $M V_m\ge2^n$。因此

$$ \mathcal N(K,d_H,m)\ge\frac{2^n}{V_m}. $$

另一方面,取极大的 $m$-separated 子集。极大性保证它是一个 $m$-net,所以

$$ \mathcal N(K,d_H,m)\le\mathcal P(K,d_H,m). $$

最后,若 $\mathcal P$ 是一个 $m$-separated 集,则以其点为中心、半径 $\lfloor m/2\rfloor$ 的 Hamming balls 两两不交。每个这样的球大小为

$$ V_{\lfloor m/2\rfloor}=\sum_{k=0}^{\lfloor m/2\rfloor}\binom nk. $$

这些不交球都包含在 cube 中,因此

$$ |\mathcal P|V_{\lfloor m/2\rfloor}\le2^n. $$

取最大的 packing,即得

$$ \mathcal P(K,d_H,m) \le \frac{2^n}{\sum_{k=0}^{\lfloor m/2\rfloor}\binom nk}. $$
Exercise 4.33 Error correcting codes 的限制
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明能纠正 $r$ 个错误的 $n$ 位编码必须有至少 $r\log_2(n/r)$ 个冗余位。

证明思路

若 code 能纠正 $r$ 个错误,则不同 codewords 的半径 $r$ Hamming balls 必须两两不交。比较这些球的总基数与整个 Hamming cube 的基数。

完整证明

一个把 $k$-bit 字符串编码为 $n$-bit 字符串的 code 有 $2^k$ 个 codewords。若它能纠正 $r$ 个错误,则以不同 codeword 为中心、半径 $r$ 的 Hamming balls 必须两两不交;否则某个接收词会同时位于两个 codeword 的 $r$ 邻域内,无法唯一解码。

每个半径 $r$ Hamming ball 的大小为

$$ \sum_{j=0}^r\binom nj \ge \binom nr. $$

不交性给出 Hamming bound:

$$ 2^k\sum_{j=0}^r\binom nj\le2^n. $$

因此

$$ 2^{n-k}\ge \binom nr. $$

使用标准组合估计 $\binom nr\ge(n/r)^r$,得到

$$ n-k \ge \log_2\binom nr \ge r\log_2\left(\frac nr\right). $$

这就是所需限制。

Exercise 4.34 $\varepsilon$-net expansion
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明任意单位向量可按 $\varepsilon$-net 点作收敛展开,并用它重证 Lemma 4.4.1。

证明思路

反复用 net 点逼近当前残差方向。每一步残差长度最多乘以 $\varepsilon$,因此系数按 $\varepsilon^k$ 衰减。

完整证明

(a) 令 $r_0=x$。若 $r_k\ne0$,从 net 中取 $x_k\in\mathcal N$ 使

$$ \left\|\frac{r_k}{\|r_k\|_2}-x_k\right\|_2\le\varepsilon, \qquad \lambda_k=\|r_k\|_2, $$

并定义 $r_{k+1}=r_k-\lambda_kx_k$。则

$$ \|r_{k+1}\|_2 = \lambda_k \left\|\frac{r_k}{\|r_k\|_2}-x_k\right\|_2 \le \varepsilon\lambda_k. $$

由 $\lambda_0=1$ 归纳得 $0\le\lambda_k\le\varepsilon^k$,且 $\|r_k\|_2\to0$。因此

$$ x=\sum_{k=0}^{\infty}\lambda_kx_k. $$

(b) 令 $M=\sup_{z\in\mathcal N}\|Az\|_2$。对任意 $x\in S^{n-1}$ 使用 (a),

$$ \|Ax\|_2 \le \sum_{k=0}^{\infty}\lambda_k\|Ax_k\|_2 \le M\sum_{k=0}^{\infty}\varepsilon^k = \frac{M}{1-\varepsilon}. $$

对 $x$ 取上确界得到 Lemma 4.4.1 的上界;下界来自 $\mathcal N\subset S^{n-1}$,即在 net 上取上确界不会超过在整个球面上取上确界。

Exercise 4.35 在 $\varepsilon$-net 上计算范数
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:用球面 $\varepsilon$-net 近似计算向量范数。

证明思路

范数是 $\sup_{u\in S^{n-1}}\langle x,u\rangle$。用 net 点逼近最大方向 $u=x/\|x\|_2$。

完整证明

下界直接由 Cauchy-Schwarz 给出:

$$ \sup_{y\in\mathcal N}\langle x,y\rangle\le\|x\|_2. $$

若 $x=0$,上界两边都为 $0$。若 $x\ne0$,令 $u=x/\|x\|_2$。取 $y\in\mathcal N$ 使 $\|u-y\|_2\le\varepsilon$。则

$$ \|x\|_2 = \langle x,u\rangle = \langle x,y\rangle+\langle x,u-y\rangle \le \sup_{z\in\mathcal N}\langle x,z\rangle + \varepsilon\|x\|_2. $$

移项得

$$ \|x\|_2 \le \frac1{1-\varepsilon} \sup_{z\in\mathcal N}\langle x,z\rangle. $$
Exercise 4.36 在 $\varepsilon$-net 上最大化二次型
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明双线性型版本的 net 近似 Lemma 4.4.2。

证明思路

取几乎达到 $\|A\|$ 的方向 $x,y$,分别用 net 点 $x_0,y_0$ 逼近。两个逼近误差各贡献至多 $\varepsilon\|A\|$。

完整证明

$$ M=\sup_{x\in\mathcal N,\ y\in\mathcal M}|\langle Ax,y\rangle|. $$

下界 $M\le\|A\|$ 来自 $\mathcal N\subset S^{n-1}$、$\mathcal M\subset S^{m-1}$。为证上界,先假设 $0<\varepsilon<1/2$;否则分母 $1-2\varepsilon$ 不为正,该形式的估计没有内容。取 $x\in S^{n-1}$、$y\in S^{m-1}$ 使 $|\langle Ax,y\rangle|$ 任意接近 $\|A\|$。选择 $x_0\in\mathcal N$、$y_0\in\mathcal M$ 满足

$$ \|x-x_0\|_2\le\varepsilon, \qquad \|y-y_0\|_2\le\varepsilon. $$

$$ \begin{aligned} |\langle Ax,y\rangle| &\le |\langle Ax_0,y_0\rangle| + |\langle A(x-x_0),y\rangle| + |\langle Ax_0,y-y_0\rangle|\\ &\le M+2\varepsilon\|A\|. \end{aligned} $$

取上确界并整理,得到

$$ \|A\|\le\frac{M}{1-2\varepsilon}. $$

若 $A$ 对称且只用同一个 net,则从 $|\langle Ax,x\rangle|$ 的变分刻画出发,同样用 $x_0$ 逼近 $x$,误差仍至多 $2\varepsilon\|A\|$,得到相同结论。

Exercise 4.37 $\varepsilon$-net 上的范数偏差
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:把 net 上的长度偏差推广到整个球面。

证明思路

当 net 偏差已经大于中心量级时,Lemma 4.4.1 直接给出粗控制;否则先控制平方范数对应的二次型,再除回长度偏差。

完整证明

$$ \Delta=\sup_{x\in\mathcal N}\bigl|\|Ax\|_2-\mu\bigr|. $$

若 $\mu<0$,则 $|\|Ax\|_2-\mu|=\|Ax\|_2+|\mu|$。由 Lemma 4.4.1,

$$ \sup_{x\in S^{n-1}}\|Ax\|_2 \le \frac{1}{1-\varepsilon} \sup_{x\in\mathcal N}\|Ax\|_2 \le \frac{\Delta}{1-\varepsilon}, $$

且 $|\mu|\le\Delta$,所以结论成立。下面只考虑 $\mu\ge0$。

若 $\Delta\ge\mu$,则 $\sup_{\mathcal N}\|Ax\|_2\le\mu+\Delta\le2\Delta$。再次由 Lemma 4.4.1,

$$ \sup_{x\in S^{n-1}}\|Ax\|_2 \le \frac{2\Delta}{1-\varepsilon}, $$

于是 $\sup_{S^{n-1}}|\|Ax\|_2-\mu|\le C\Delta/(1-2\varepsilon)$。

剩下情形是 $0\le\Delta<\mu$。设

$$ B=A^{\mathsf T}A-\mu^2I_n. $$

对 $x\in\mathcal N$,

$$ |x^{\mathsf T}Bx| = \bigl|\|Ax\|_2^2-\mu^2\bigr| = \bigl|\|Ax\|_2-\mu\bigr|(\|Ax\|_2+\mu) \le \Delta(2\mu+\Delta) \le 3\mu\Delta. $$

由 Lemma 4.4.2 的对称情形,

$$ \|B\| \le \frac{3\mu\Delta}{1-2\varepsilon}. $$

因此对任意 $x\in S^{n-1}$,

$$ \bigl|\|Ax\|_2-\mu\bigr| = \frac{\bigl|\|Ax\|_2^2-\mu^2\bigr|}{\|Ax\|_2+\mu} \le \frac{\|B\|}{\mu} \le \frac{3\Delta}{1-2\varepsilon}. $$

合并三种情形即可得到题设不等式,其中 $C$ 为绝对常数。

Exercise 4.38 $\varepsilon$-net 的 convex hull
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明球面 $\varepsilon$-net 的凸包包含 $(1-\varepsilon)B_2^n$。

证明思路

对任意方向 $u$,net 中有点 $z$ 与 $u$ 接近,因此 $\langle u,z\rangle\ge1-\varepsilon$。这说明凸包在每个方向上的支撑函数至少为 $1-\varepsilon$。

完整证明

令 $C=\operatorname{conv}(\mathcal N)$。对任意 $u\in S^{n-1}$,取 $z\in\mathcal N$ 使 $\|u-z\|_2\le\varepsilon$。由

$$ \|u-z\|_2^2 = 2-2\langle u,z\rangle $$

可得 $\langle u,z\rangle\ge1-\varepsilon^2/2\ge1-\varepsilon$。因此

$$ \sup_{w\in C}\langle u,w\rangle\ge1-\varepsilon \qquad\text{for all }u\in S^{n-1}. $$

若存在 $y\in(1-\varepsilon)B_2^n$ 但 $y\notin C$,由超平面分离定理,存在 $u\in S^{n-1}$ 使

$$ \langle u,y\rangle> \sup_{w\in C}\langle u,w\rangle. $$

但 $\langle u,y\rangle\le\|y\|_2\le1-\varepsilon$,与上式的支撑函数下界矛盾。因此 $(1-\varepsilon)B_2^n\subset C$。

Exercise 4.39 随机点形成 $\varepsilon$-net
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明指数多个标准 Gaussian 点经 $\sqrt n$ 缩放后,以高概率在固定半径球内形成 $\varepsilon$-net。

证明思路

先取 $RB_2^n$ 的一个确定性细 net。对每个 net 中心,Gaussian 缩放点落入它附近的概率至少为 $e^{-C(R,\varepsilon)n}$。指数多个独立点加 union bound 即覆盖所有中心。

完整证明

令 $Z_i=g_i/\sqrt n$。取 $RB_2^n$ 的一个 $(\varepsilon/4)$-net $\mathcal M$,满足

$$ |\mathcal M|\le\left(\frac{CR}{\varepsilon}\right)^n. $$

对每个 $y\in\mathcal M$,把它稍微向内收缩为 $y'=(1-\varepsilon/(8R))y$。若 $\|Z_i-y'\|_2\le\varepsilon/8$,则 $Z_i\in RB_2^n$,且 $\|Z_i-y\|_2\le\varepsilon/4$。

Gaussian 密度在固定半径球内的指数阶下界给出

$$ \mathbb P\{\|Z_i-y'\|_2\le\varepsilon/8\} \ge \exp[-C(R,\varepsilon)n] =:p. $$

这里可以直接验证:在球 $B(y',\varepsilon/8)$ 中所有点范数不超过 $R$,所以 $Z_i$ 的密度至少为 $(n/2\pi)^{n/2}e^{-nR^2/2}$;再用 $\operatorname{Vol}(B_2^n)\ge(c/\sqrt n)^n$,得到上述指数型下界。

因此某个固定 $y\in\mathcal M$ 没有被任何样本点命中的概率至多 $(1-p)^N\le e^{-Np}$。由 union bound,所有 net 中心都被命中的失败概率至多

$$ |\mathcal M|e^{-Np}. $$

只要 $N\ge\exp(C'(R,\varepsilon)n)$,上式至多为 $e^{-cn}$。在好事件上,每个 $x\in RB_2^n$ 先找 $y\in\mathcal M$ 使 $\|x-y\|_2\le\varepsilon/4$,再找样本点 $Z_i\in RB_2^n$ 使 $\|Z_i-y\|_2\le\varepsilon/4$,于是 $\|x-Z_i\|_2\le\varepsilon$。这正是所需 net 性质。

Exercise 4.40 过多随机点不处于 convex position
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明当 Gaussian 点数指数级过大时,它们高概率不全是凸包顶点。

证明思路

把点分成两批。第一批足够多,缩放后在 $2B_2^n$ 内形成细 net,从而其凸包包含一个固定小球;第二批中高概率至少有一个点落入这个小球,于是该点在其他点的凸包内。

完整证明

缩放不改变 convex position,所以考虑 $Z_i=g_i/\sqrt n$。把样本分成两组,各有指数多个点。由 Exercise 4.39,若第一组大小至少为 $e^{C_1n}$,则以至少 $1-e^{-c_1n}$ 的概率,第一组中落在 $2B_2^n$ 内的点构成 $2B_2^n$ 的一个足够小的 net。

特别地,它包含球面 $S^{n-1}$ 的一个 $\varepsilon$-net。由 Exercise 4.38,第一组点的凸包包含 $(1-\varepsilon)B_2^n$。

另一方面,对单个 $Z_i$,事件 $\|Z_i\|_2\le1-\varepsilon$ 有概率至少 $e^{-C_2n}$。如果第二组大小至少为 $e^{C_3n}$ 且 $C_3>C_2$,则第二组至少有一个点落在 $(1-\varepsilon)B_2^n$ 内的概率至少为 $1-e^{-c_2n}$。

在两件好事件同时发生时,第二组中的这个点属于第一组点的凸包,因此它不是整个点集凸包的极点。于是这些 Gaussian 点不处于 convex position。选择 $C$ 足够大并合并失败概率,得到至少 $1-e^{-cn}$ 的结论。

Exercise 4.41 随机矩阵 operator norm 的期望
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:把 Theorem 4.4.3 和 Theorem 4.6.1 的高概率界积分成期望界。

证明思路

使用尾积分公式:若 $Z\le a+bt$ 以概率至少 $1-2e^{-t^2}$ 成立,则 $\mathbb EZ\le a+Cb$。

完整证明

(a) Theorem 4.4.3 给出对所有 $t\ge0$,

$$ \mathbb P\left\{ \|A\|>CK(\sqrt m+\sqrt n+t) \right\} \le 2e^{-t^2}. $$

由尾积分公式,

$$ \mathbb E\|A\| \le CK(\sqrt m+\sqrt n) + CK\int_0^\infty 2e^{-t^2}\,dt \le C'K(\sqrt m+\sqrt n). $$

(b) Theorem 4.6.1 证明中的矩阵形式给出

$$ \left\| \frac1mA^{\mathsf T}A-I_n \right\| \le CK^2\left( \sqrt{\frac nm}+\frac{t}{\sqrt m} + \frac nm+\frac{t^2}{m} \right) $$

以至少 $1-2e^{-t^2}$ 的概率成立。再次积分尾部,得到

$$ \mathbb E \left\| \frac1mA^{\mathsf T}A-I_n \right\| \le CK^2\left( \sqrt{\frac nm}+\frac nm \right). $$

同理,把 Theorem 4.6.1 的奇异值双侧尾界积分,得到

$$ \sqrt m-CK^2\sqrt n \le \mathbb E\,s_n(A) \le \mathbb E\,s_1(A) \le \sqrt m+CK^2\sqrt n. $$
Exercise 4.42 随机矩阵范数的下界
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明随机矩阵算子范数至少有 $\sqrt m+\sqrt n$ 的量级。

证明思路

算子范数至少不小于任意一列和任意一行的 Euclidean 范数。分别对第一列和第一行使用范数集中下界。

完整证明

对第一列 $Ae_1$,它是 $m$ 维独立、方差为一、次高斯坐标向量。范数集中下界给出

$$ \|Ae_1\|_2 \ge \sqrt m-t $$

以至少 $1-\exp(-ct^2/K^4)$ 的概率成立。同理,对第一行 $A^{\mathsf T}e_1$,

$$ \|A^{\mathsf T}e_1\|_2 \ge \sqrt n-t $$

也以同样概率成立。两者同时成立的概率至少为 $1-2\exp(-ct^2/K^4)$。

在该事件上,

$$ \|A\| \ge \max\{\|Ae_1\|_2,\|A^{\mathsf T}e_1\|_2\} \ge \max\{\sqrt m-t,\sqrt n-t\}. $$

由于 $\max(a,b)\ge(a+b)/2$,得到

$$ \|A\| \ge \frac12(\sqrt m+\sqrt n-2t). $$

把 $2t$ 重新记为 $t$ 并调整常数,即得题设形式。

Exercise 4.43 Subgaussian 矩阵上界:放松独立性
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:说明 Theorem 4.4.3 的 net proof 只需要固定双线性型是次高斯,不一定需要所有元素独立。

证明思路

对固定 $x,y$,$\langle Ax,y\rangle$ 是矩阵随机向量 $A$ 在 Frobenius 单位方向 $yx^{\mathsf T}$ 上的投影。只要这个投影次高斯,后续 net 和 union bound 完全相同。

完整证明

(a) 把 $A$ 看成 $\mathbb R^{mn}$ 中的随机向量。对固定 $x\in S^{n-1}$、$y\in S^{m-1}$,

$$ \langle Ax,y\rangle = \langle A,yx^{\mathsf T}\rangle_F, \qquad \|yx^{\mathsf T}\|_F=1. $$

因此

$$ \|\langle Ax,y\rangle\|_{\psi_2} \le \|A\|_{\psi_2}. $$

对两个球面的 $1/4$-net 使用 union bound,并用 Lemma 4.4.2 从 net 回到整个球面,得到 Theorem 4.4.3 的同阶结论,其中 $K=\|A\|_{\psi_2}$。

(b) 若 $A$ 的行 $A_i$ 独立、均值为零且次高斯,则

$$ \langle Ax,y\rangle = \sum_{i=1}^m y_i\langle A_i,x\rangle. $$

这些项独立、均值为零,并满足 $\|\langle A_i,x\rangle\|_{\psi_2}\le K$,其中 $K=\max_i\|A_i\|_{\psi_2}$。独立次高斯和的估计给出 $\|\langle Ax,y\rangle\|_{\psi_2}\le CK$。于是 (a) 中的 net argument 仍然适用。列独立情形对 $A^{\mathsf T}$ 使用同样论证。

Exercise 4.44 随机矩阵的一些 $p\to q$ 范数
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:估计独立次高斯随机矩阵的若干 $p\to q$ 范数,并说明 Gaussian 情形下量级最优。

证明思路

$1\to\infty$ 范数就是最大元素;$1\to2$ 范数就是最大列范数。分别用最大次高斯变量和列范数集中。

完整证明

(a) 由定义,

$$ \|A\|_{1\to\infty} = \max_{i,j}|A_{ij}|. $$

对 $mn$ 个次高斯变量使用 union bound,得到以高概率

$$ \max_{i,j}|A_{ij}| \le CK\sqrt{\log(mn)} \le CK(\sqrt{\log m}+\sqrt{\log n}). $$

(b) 对任意矩阵,

$$ \|A\|_{1\to2} = \max_{j\le n}\|Ae_j\|_2, \qquad \|A^{\mathsf T}\|_{2\to\infty} = \max_{j\le n}\|Ae_j\|_2, $$

所以二者相等。每一列是 $m$ 维独立、方差为一、次高斯坐标向量,范数集中给出

$$ \|Ae_j\|_2 \le \sqrt m+CK^2t $$

以概率至少 $1-2e^{-t^2}$ 成立。取 $t=C\sqrt{\log n}$ 并对列作 union bound,得到

$$ \|A\|_{1\to2} \le \sqrt m+CK^2\sqrt{\log n} $$

以高概率成立。

(c) 若 $A_{ij}$ 为独立标准正态,则最大元素满足

$$ \max_{i,j}|A_{ij}| \ge c\sqrt{\log(mn)} \ge c'(\sqrt{\log m}+\sqrt{\log n}) $$

以高概率成立,证明 (a) 的量级最优。另一方面,

$$ \|A\|_{1\to2} = \max_j\|Ae_j\|_2 \ge \|Ae_1\|_2 \ge c\sqrt m $$

以高概率成立;同时 $\|A\|_{1\to2}\ge\max_j|A_{1j}|\ge c\sqrt{\log n}$。合并两个下界,得到

$$ \|A\|_{1\to2} \ge c(\sqrt m+\sqrt{\log n}) $$

以高概率成立。

Exercise 4.45 Community detection
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:去掉 Theorem 4.5.2 中 $q$ 不太小的假设,仍得到 $C/(p-q)^2$ 个误分类顶点。

证明思路

原证明中 $q$ 出现在 $\lambda_1(D)$ 与 $\lambda_2(D)$ 的间隔。改用中心化或投影到 $\mathbf1^\perp$ 的谱聚类,只追踪 community 方向,此时谱间隙为 $(p-q)n/2$。

完整证明

设 $\mathbf1=(1,\ldots,1)$,令 $P$ 是到 $\mathbf1^\perp$ 的正交投影。对邻接矩阵 $A=D+R$,考虑中心化矩阵

$$ \widetilde A=PAP. $$

在 $\mathbf1^\perp$ 上,信号部分 $\widetilde D=PDP$ 只有一个非零特征方向,即 community 向量

$$ u=(1,\ldots,1,-1,\ldots,-1), $$

对应特征值

$$ \lambda=\frac{p-q}{2}n. $$

因此在 $\mathbf1^\perp$ 上的相关谱间隙是 $\lambda$,不再涉及 $q$。

噪声部分为 $\widetilde R=PRP$,且 $\|\widetilde R\|\le\|R\|$。由于 Bernoulli 噪声有界、中心化且次高斯,Corollary 4.4.7 给出

$$ \|\widetilde R\|\le C\sqrt n $$

以高概率成立。Davis-Kahan 不等式给出单位化特征向量的误差

$$ \sin\angle(\bar u,\bar v) \le \frac{C\sqrt n}{(p-q)n} = \frac{C}{(p-q)\sqrt n}, $$

其中 $\bar v$ 是 $\widetilde A$ 在 $\mathbf1^\perp$ 上 top eigenvector。乘以 $\sqrt n$ 后,存在全局符号 $\theta$ 使

$$ \|u-\theta\sqrt n\,\bar v\|_2 \le \frac{C}{p-q}. $$

每个误分类坐标至少贡献一个常数量级的平方误差,因此误分类数至多

$$ \frac{C}{(p-q)^2}. $$

算法就是:计算 $PAP$ 在 $\mathbf1^\perp$ 上的 top eigenvector,并按其坐标符号分类。

Exercise 4.46 随机矩阵双边界的另一个证明
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:用固定方向的范数集中加 net argument 重新证明次高斯矩阵奇异值双侧界。

证明思路

固定 $x$ 时,$Ax$ 是有独立次高斯坐标的向量,Theorem 3.1.1 控制 $\|Ax\|_2$。再在球面 net 上 union bound,并用 Exercise 4.37 推回所有方向。

完整证明

固定 $x\in S^{n-1}$。矩阵行 $A_i$ 独立、各向同性、次高斯,因此随机变量 $\langle A_i,x\rangle$ 独立、均值为零、方差为一且次高斯范数由 $K$ 控制。Theorem 3.1.1 给出

$$ \mathbb P\left\{ \bigl|\|Ax\|_2-\sqrt m\bigr|>CK^2t \right\} \le 2e^{-t^2}. $$

取 $S^{n-1}$ 的 $1/4$-net $\mathcal N$,满足 $|\mathcal N|\le9^n$。令 $t=C(\sqrt n+u)$,对 net 点作 union bound,可得以至少 $1-2e^{-u^2}$ 的概率,

$$ \sup_{x\in\mathcal N} \bigl|\|Ax\|_2-\sqrt m\bigr| \le CK^2(\sqrt n+u). $$

由 Exercise 4.37,把 net 上的偏差推广到整个球面:

$$ \sup_{x\in S^{n-1}} \bigl|\|Ax\|_2-\sqrt m\bigr| \le CK^2(\sqrt n+u). $$

最后注意

$$ s_1(A)=\sup_{x\in S^{n-1}}\|Ax\|_2, \qquad s_n(A)=\inf_{x\in S^{n-1}}\|Ax\|_2. $$

于是得到 Theorem 4.6.1 的双侧奇异值界。

Exercise 4.47 随机矩阵的中间 singular values
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:从全维最小奇异值界推出第 $k$ 个奇异值的下界。

证明思路

把 $A$ 限制到任意固定 $k$ 维子空间,例如前 $k$ 个坐标。限制矩阵是 $m\times k$ 次高斯矩阵;其最小奇异值给出原矩阵第 $k$ 个奇异值的下界。

完整证明

令 $E\subset\mathbb R^n$ 为前 $k$ 个坐标张成的子空间,令 $A_E$ 表示 $A$ 在 $E$ 上的限制,即 $A$ 的前 $k$ 列。它是 $m\times k$ 随机矩阵,行仍独立、各向同性且次高斯范数由 $K$ 控制。

由 Theorem 4.6.1 应用于 $A_E$,以至少 $1-2e^{-t^2}$ 的概率,

$$ s_k(A_E) \ge \sqrt m-CK^2(\sqrt k+t). $$

由 Courant-Fischer 变分刻画,

$$ s_k(A) = \max_{\dim F=k} \min_{x\in F\cap S^{n-1}}\|Ax\|_2 \ge \min_{x\in E\cap S^{n-1}}\|Ax\|_2 = s_k(A_E). $$

合并即得

$$ s_k(A) \ge \sqrt m-CK^2(\sqrt k+t). $$
Exercise 4.48 相对误差下的 covariance estimation
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:把 covariance estimation 改写成白化后的相对误差估计。

证明思路

令 $Y=\Sigma^{-1/2}X$。则 $Y$ 各向同性,样本协方差为 $\Sigma^{-1/2}\Sigma_m\Sigma^{-1/2}$。相对误差正是这个白化样本协方差与 $I$ 的算子范数距离。

完整证明

$$ Y=\Sigma^{-1/2}X, \qquad Y_i=\Sigma^{-1/2}X_i. $$

则 $\mathbb EYY^{\mathsf T}=I_n$,并且

$$ \frac1m\sum_{i=1}^mY_iY_i^{\mathsf T} = \Sigma^{-1/2}\Sigma_m\Sigma^{-1/2} =:\widetilde\Sigma_m. $$

条件 (4.29) 在白化后给出 $\|Y\|_{\psi_2}\le K$。由 Theorem 4.7.1 应用于 $Y$,

$$ \mathbb E\|\widetilde\Sigma_m-I\| \le CK^2\left(\sqrt{\frac nm}+\frac nm\right). $$

另一方面,对任意 $v\ne0$,令 $w=\Sigma^{1/2}v$,则

$$ \frac{v^{\mathsf T}\Sigma_m v}{v^{\mathsf T}\Sigma v} = \frac{w^{\mathsf T}\widetilde\Sigma_m w}{\|w\|_2^2}. $$

因此

$$ \sup_{v\in\mathbb R^n} \left| \frac{v^{\mathsf T}\Sigma_m v}{v^{\mathsf T}\Sigma v}-1 \right| = \|\widetilde\Sigma_m-I\|. $$

取期望即得所需相对误差界。

Exercise 4.49 High-probability covariance estimation
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:验证 covariance estimation 的高概率版本。

证明思路

Theorem 4.7.1 的证明本质上对白化样本矩阵应用 Theorem 4.6.1。保留其中的尾参数并令 $t=\sqrt u$ 即得高概率形式。

完整证明

先假设 $\Sigma$ 可逆,令 $Y=\Sigma^{-1/2}X$。如 Exercise 4.48,$Y$ 各向同性且次高斯参数由 $K$ 控制。Theorem 4.6.1 的矩阵形式给出:对任意 $u\ge0$,以至少 $1-2e^{-u}$ 的概率,

$$ \left\| \frac1m\sum_{i=1}^mY_iY_i^{\mathsf T}-I \right\| \le CK^2\left( \sqrt{\frac{n+u}{m}} + \frac{n+u}{m} \right). $$

乘回 $\Sigma^{1/2}$,得到 Loewner 意义下

$$ -\delta\Sigma \preceq \Sigma_m-\Sigma \preceq \delta\Sigma, $$

其中 $\delta$ 为右侧。于是

$$ \|\Sigma_m-\Sigma\| \le \delta\|\Sigma\|. $$

若 $\Sigma$ 不可逆,可在 $\Sigma+\eta I$ 上做上述论证后令 $\eta\downarrow0$,或限制到 $\Sigma$ 的值域。于是得到 Remark 4.7.3 的高概率界。

Exercise 4.50 低秩矩阵的 covering numbers
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明秩 $r$、Frobenius 范数至多 $1$ 的矩阵类覆盖数按约 $(m+n)r$ 维增长。

证明思路

上界直接离散化 SVD 的三个部分:左奇异向量、右奇异向量、奇异值。下界不必使用流形理论;只要把一个全维矩阵球嵌入低秩矩阵类即可。

完整证明

(a) 注意 $\|U-V\|_{1\to2}$ 就是列差的最大 Euclidean 范数:

$$ \|U-V\|_{1\to2} = \max_{j\le r}\|U_{:j}-V_{:j}\|_2. $$

取单位球面 $S^{m-1}$ 的一个 $(\varepsilon/2)$-net $\mathcal N$,满足 $|\mathcal N|\le(C/\varepsilon)^m$。对 $U\in O_{m,r}$ 的每一列 $U_{:j}$,选择 $u_j\in\mathcal N$ 使 $\|U_{:j}-u_j\|_2\le\varepsilon/2$。所有 $r$ 元组 $(u_1,\ldots,u_r)$ 构成 $O_{m,r}$ 的一个外部 $(\varepsilon/2)$-net;这些近似矩阵本身不必仍在 $O_{m,r}$ 中。由 Exercise 4.25,可把外部中心移回 $O_{m,r}$,并得到内部 $\varepsilon$-net,数量不增加。于是

$$ \mathcal N(O_{m,r},\|\cdot\|_{1\to2},\varepsilon) \le \left(\frac C\varepsilon\right)^{mr}. $$

(b) 任意 $A\in M_{m,n,r}$ 可写成 SVD

$$ A=U\operatorname{diag}(s)V^{\mathsf T}, \qquad U\in O_{m,r},\quad V\in O_{n,r},\quad \|s\|_2\le1. $$

分别取 $O_{m,r}$、$O_{n,r}$ 和 Euclidean 单位球 $B_2^r$ 的 $\varepsilon/3$-nets。若 $U,V,s$ 分别由 $U_0,V_0,s_0$ 逼近,则

$$ \|U\operatorname{diag}(s)V^{\mathsf T} - U_0\operatorname{diag}(s_0)V_0^{\mathsf T}\|_F \le \|U-U_0\|_{1\to2}\|s\|_2 + \|s-s_0\|_2 + \|V-V_0\|_{1\to2}\|s_0\|_2 \le \varepsilon. $$

三个 net 的大小相乘,得到

$$ \mathcal N(M_{m,n,r},\|\cdot\|_F,\varepsilon) \le \left(\frac C\varepsilon\right)^{mr} \left(\frac C\varepsilon\right)^{nr} \left(\frac C\varepsilon\right)^r = \left(\frac C\varepsilon\right)^{(m+n+1)r}. $$

(c) 先假设 $m\ge n$。考虑所有只在前 $r$ 列可能非零的矩阵

$$ A=[B\ 0], \qquad B\in\mathbb R^{m\times r}, \qquad \|B\|_F\le1. $$

其中满列秩的 $B$ 给出 rank 正好为 $r$ 的矩阵,且满列秩矩阵在 $\mathbb R^{mr}$ 的单位 Frobenius 球中占满测度。若 $M_{m,n,r}$ 可由 $N$ 个 Frobenius 半径 $\varepsilon$ 的球覆盖,则这些球也覆盖上述满秩子集。忽略测度为零的秩亏集合并比较 $\mathbb R^{mr}$ 中体积,得到

$$ N\cdot \operatorname{Vol}(\varepsilon B_2^{mr}) \ge \operatorname{Vol}(B_2^{mr}), $$

所以 $N\ge(c/\varepsilon)^{mr}$(当 $\varepsilon$ 大于固定常数时右侧至多为 $1$,结论平凡)。由于 $m\ge n$,有 $mr\ge (m+n)r/2$,故

$$ \mathcal N(M_{m,n,r},\|\cdot\|_F,\varepsilon) \ge \left(\frac c\varepsilon\right)^{(m+n)r/2}. $$

若 $n>m$,则改为只让前 $r$ 行可能非零,得到 $N\ge(c/\varepsilon)^{nr}$,同样推出所需下界。这说明上界中的指数阶在常数因子意义下正确。

Exercise 4.51 Gaussian mixture model 的 spectral clustering
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明样本协方差的 top eigenvector 可以学习 Gaussian mixture model 的分离方向,从而实现谱聚类。

证明思路

总体协方差是 $I+\mu\mu^{\mathsf T}$,top eigenvector 是 $\mu$ 方向。协方差估计给出样本 top eigenvector 接近该方向;投影到近似方向后,Gaussian 噪声不能改变大多数点的符号。

完整证明

(a) 写 $X=\theta\mu+g$,其中 $\theta$ 是 Rademacher 变量,$g\sim N(0,I_n)$。由于 $\theta$ 与 $g$ 独立且均值为零,

$$ \Sigma=\mathbb EXX^{\mathsf T} = \mu\mu^{\mathsf T}+I_n. $$

因此 top eigenvector 为

$$ u=\frac{\mu}{\|\mu\|_2}, $$

对应特征值 $1+\|\mu\|_2^2$;其余方向特征值为 $1$,谱间隙为 $\|\mu\|_2^2$。

(b) 若使用真实方向 $u$ 分类,则

$$ \langle X,u\rangle = \theta\|\mu\|_2+\langle g,u\rangle. $$

令 $M=\|\mu\|_2$。若真实标签为 $\theta_i$,则

$$ \theta_i\langle X_i,u\rangle = M+\theta_i\langle g_i,u\rangle. $$

因此真实方向误分类概率为 $\mathbb P\{N(0,1)\le-M\}$。更强地,真实方向 margin 小于 $M/2$ 的概率为 $\mathbb P\{N(0,1)\le-M/2\}$。当 $M\ge C$ 且 $C$ 足够大时,这个概率小于 $0.001$。对 $m$ 个样本用 Chernoff 或 Hoeffding 不等式,可得至少以概率 $0.997$,满足 $\theta_i\langle X_i,u\rangle<M/2$ 的样本比例至多 $0.3\%$。

(c) 该 mixture 分布满足 Theorem 4.7.1 的相对次高斯条件,且参数为绝对常数:对任意 $x$,

$$ \|\langle X,x\rangle\|_{\psi_2} \le C\left(|\langle\mu,x\rangle|+\|x\|_2\right) \le C\|\langle X,x\rangle\|_{L^2}. $$

若 $m\ge Cn$ 且常数 $C$ 足够大,covariance estimation 的高概率版本给出

$$ \|\Sigma_m-\Sigma\| \le c\|\mu\|_2^2 $$

以概率至少 $0.995$ 成立,其中 $c>0$ 可任意小。

(d) 由 Davis-Kahan 不等式和谱间隙 $\|\mu\|_2^2$,样本 top eigenvector $v$ 满足,在差一个全局符号后,

$$ \|v-u\|_2\le \eta, $$

其中 $\eta>0$ 可通过增大样本量常数取得足够小。同时在同一好事件上,

$$ \|\Sigma_m\| \le \|\Sigma\|+\|\Sigma_m-\Sigma\| \le C(1+M^2). $$

(e) 现在处理关键的依赖性问题:$v$ 是由样本算出的,不能把 $\langle X_i,v-u\rangle$ 当作固定方向上的独立噪声。我们改用确定性计数。记 $h=v-u$。若某个样本满足

$$ \theta_i\langle X_i,u\rangle\ge M/2 \quad\text{且}\quad |\langle X_i,h\rangle|<M/2, $$

则它按 $v$ 的符号也会被正确分类,因为

$$ \theta_i\langle X_i,v\rangle = \theta_i\langle X_i,u\rangle + \theta_i\langle X_i,h\rangle > 0. $$

因此额外错误只能来自 $|\langle X_i,h\rangle|\ge M/2$ 的样本。由 Markov 型计数,

$$ \#\{i:|\langle X_i,h\rangle|\ge M/2\} \le \frac{4}{M^2} \sum_{i=1}^m\langle X_i,h\rangle^2 = \frac{4m}{M^2}h^{\mathsf T}\Sigma_m h. $$

在 (d) 的好事件上,

$$ h^{\mathsf T}\Sigma_m h \le \|\Sigma_m\|\|h\|_2^2 \le C(1+M^2)\eta^2. $$

由于 $M\ge C\ge1$,额外错误比例至多 $C'\eta^2$。选择样本量常数使 $\eta$ 足够小,可令它小于 $0.3\%$。再加上 (b) 中真实方向 margin 不足的样本,误分类比例至多 $1\%$。合并两个好事件,概率至少 $0.99$。这证明了 Theorem 4.7.5。

学习检查表

  • [ ] 能从 $A^{\mathsf T}A$ 推出 SVD。
  • [ ] 能解释算子范数、最大奇异值和双线性形式之间的等价关系。
  • [ ] 能写出 net argument 的四步证明模板。
  • [ ] 能解释为什么随机矩阵范数量级是 $\sqrt m+\sqrt n$。
  • [ ] 能把社区检测拆成“信号矩阵、噪声矩阵、扰动定理”。
  • [ ] 能说明协方差估计为什么需要控制样本矩阵的奇异值。
  • [ ] 能从译文、学习笔记或并排页跳转到对应证明位置。

后续衔接

第 4 章当前已经完成结构化学习页、关键定理证明区、51 个习题完整证明和 proof link 锚点。后续若继续精修,应优先复核应用题中的常数条件,以及低秩矩阵 covering number 下界的体积论证细节。