精校翻译 Ch.5 无独立性集中
第 5 章精校翻译:无独立性的集中

翻译状态说明

本页为第 5 章的精校翻译,覆盖正文、Notes、Exercises 5.1-5.32、原书 Figure 5.1-5.2,并为原文明确证明、正文隐藏验证和习题加入学习笔记证明跳转。

第 5 章 无独立性的集中

到目前为止,我们研究集中不等式时严重依赖随机变量之间的独立性。现在,我们转向另外一类方法:它们不依赖独立性。第 5.1 节以 Euclidean 球面为例,引入等周不等式方法;第 5.2 节再介绍其他度量测度空间上的集中现象。

第 5.3 节利用球面上的集中推出经典的 Johnson-Lindenstrauss 引理,这是高维数据降维中的核心结果。

第 5.4 节介绍矩阵集中不等式,重点是矩阵 Bernstein 不等式,它把经典 Bernstein 不等式推广到随机矩阵。第 5.5 与 5.6 节再把这一工具用于稀疏网络中的社群检测,以及更一般分布下的协方差估计。

不要跳过习题。本章习题会研究二值随机投影的降维版本(Exercise 5.14)、矩阵函数演算(Exercises 5.16-5.19)、多个矩阵集中不等式(Exercises 5.20-5.24),以及矩阵 sketching(Exercise 5.32)、社群检测(Exercise 5.25)等应用。

5.1 球面上 Lipschitz 函数的集中

给定随机向量 $X\in\mathbb R^n$ 和函数 $f:\mathbb R^n\to\mathbb R$,什么时候随机变量 $f(X)$ 会集中,即

$$ f(X)\approx \mathbb Ef(X) \qquad \text{with high probability}? $$

如果 $X$ 是正态向量且 $f$ 是线性函数,这很容易:$f(X)$ 仍是正态随机变量(Corollary 3.3.2),并且有很好的集中(Proposition 2.1.2)。

非线性函数 $f$ 呢?不能指望任意函数都有好集中;例如函数可以在极小区域内剧烈振荡。

查看学习笔记:为什么任意函数不可能都集中

不过,如果 $f$ 的振荡不太剧烈,我们就可以期待集中。为精确表达这一点,我们引入 Lipschitz 函数。它们正是用来排除极端振荡的。

5.1.1 Lipschitz 函数

Definition 5.1.1 Lipschitz 函数

设 $(X,d_X)$ 和 $(Y,d_Y)$ 是度量空间。如果存在 $L\in\mathbb R$,使得对所有 $u,v\in X$ 都有

$$ d_Y(f(u),f(v))\le L\,d_X(u,v), $$

则称函数 $f:X\to Y$ 是 Lipschitz 的。所有这类 $L$ 的下确界称为 $f$ 的 Lipschitz 范数,记为 $\|f\|_{\mathrm{Lip}}$。

换句话说,Lipschitz 函数不会把距离拉得太大。当 $\|f\|_{\mathrm{Lip}}\le1$ 时,它们是 contraction,因为它们只会缩短距离。Lipschitz 函数位于可微函数和一致连续函数之间:

$$ f\text{ 可微} \Longrightarrow f\text{ Lipschitz} \Longrightarrow f\text{ 一致连续}. $$

在 Exercise 5.1 中,你还会对 $f:\mathbb R^n\to\mathbb R$ 定量化第一个蕴含:

$$ \|f\|_{\mathrm{Lip}} \le \sup_{x\in\mathbb R^n}\|\nabla f(x)\|_2. $$

Example 5.1.2 向量、矩阵与范数给出的 Lipschitz 函数

(a) 固定 $\theta\in\mathbb R^n$。线性泛函

$$ f(x)=\langle x,\theta\rangle $$

的 Lipschitz 范数为 $\|f\|_{\mathrm{Lip}}=\|\theta\|_2$。

(b) 更一般地,对任意 $m\times n$ 矩阵 $A$,线性算子 $f(x)=Ax$ 的 Lipschitz 范数为 $\|A\|$。

(c) 对 $\mathbb R^n$ 上任意范数 $\|\cdot\|$,函数 $f(x)=\|x\|$ 的 Lipschitz 范数等于最小的 $L$,使得

$$ \|x\|\le L\|x\|_2 \qquad \text{for all }x\in\mathbb R^n. $$ 查看学习笔记:Example 5.1.2 验证

5.1.2 通过等周不等式得到集中

我们现在证明 Euclidean 球面

$$ S^{n-1}=\{x\in\mathbb R^n:\|x\|_2=1\} $$

上的任意 Lipschitz 函数都会集中。

Theorem 5.1.3 球面上 Lipschitz 函数的集中

设 $X$ 均匀分布在半径为 $\sqrt n$ 的 Euclidean 球面上,即 $X\sim\operatorname{Unif}(\sqrt n S^{n-1})$。那么,对任意 Lipschitz 函数 $f:\sqrt n S^{n-1}\to\mathbb R$,都有

$$ \|f(X)-\mathbb Ef(X)\|_{\psi_2} \le C\|f\|_{\mathrm{Lip}}. $$

等价地,对所有 $t\ge0$,

$$ \mathbb P\{|f(X)-\mathbb Ef(X)|\ge t\} \le 2\exp\left(-\frac{ct^2}{\|f\|_{\mathrm{Lip}}^2}\right). $$
阅读重点:这里没有坐标独立性。集中来自球面几何,而不是来自独立和。
查看学习笔记完整证明

我们已经在线性函数情形证明过 Theorem 5.1.3。Theorem 3.4.5 告诉我们,$X$ 是次高斯随机向量;按定义,$X$ 的任意线性函数都是次高斯随机变量。

为了完全证明 Theorem 5.1.3,需要说明任意 Lipschitz 函数至少像线性函数一样集中。我们不直接比较函数值,而是比较它们的 sublevel sets,即给定水平 $a$ 时球面上满足 $f(x)\le a$ 的区域。对线性函数而言,这些区域就是 spherical caps。比较一般集合和 spherical caps 的面积,需要一个重要几何原则:等周不等式。

Theorem 5.1.4 $\mathbb R^n$ 上的等周不等式

在 $\mathbb R^n$ 中,在所有给定体积的集合 $A$ 里,Euclidean 球具有最小表面积。更强地,对任意 $\varepsilon>0$,Euclidean 球也使 $A$ 的 $\varepsilon$-邻域体积最小,其中

$$ A_\varepsilon = \{x\in\mathbb R^n:\exists y\in A,\ \|x-y\|_2\le\varepsilon\} = A+\varepsilon B_2^n. $$

Theorem 5.1.4 的“更强”部分推出第一部分;令 $\varepsilon\to0$ 即可看出。Figure 5.1 给出了等周不等式的示意。

等周不等式中集合邻域体积比较
图 5.1 等周不等式:在所有给定体积的集合 $A$ 中,Euclidean 球使 $\varepsilon$-邻域 $A_\varepsilon$ 的体积最小。

球面 $S^{n-1}$ 上也有类似等周不等式,此时最小化者是 spherical caps,也就是某个点的邻域。用 $\sigma_{n-1}$ 表示球面 $S^{n-1}$ 上的归一化面积测度。

Theorem 5.1.5 球面上的等周不等式

设 $\varepsilon>0$。在所有具有给定面积 $\sigma_{n-1}(A)$ 的集合 $A\subset S^{n-1}$ 中,spherical caps 使邻域面积 $\sigma_{n-1}(A_\varepsilon)$ 最小,其中

$$ A_\varepsilon = \{x\in S^{n-1}:\exists y\in A,\ \|x-y\|_2\le\varepsilon\}. $$

本书不证明等周不等式(Theorem 5.1.4 与 5.1.5);本章 bibliographic notes 会给出若干已知证明的参考。

5.1.3 球面上集合的 blow-up

等周不等式推出一个非常反直觉的现象:如果集合 $A$ 至少覆盖球面一半面积,那么它的 $\varepsilon$-邻域 $A_\varepsilon$ 会覆盖球面的大部分。我们先陈述并证明这个 blow-up 现象,然后解释直觉。为了适配 Theorem 5.1.3,我们在半径为 $\sqrt n$ 的球面上工作。

Lemma 5.1.6 Blow-up

设 $A\subset\sqrt n S^{n-1}$,$\sigma$ 是该球面上的归一化面积测度。若 $\sigma(A)\ge1/2$,则对所有 $t\ge0$,

$$ \sigma(A_t)\ge1-2\exp(-ct^2). $$ 查看学习笔记完整证明
Proof 用半球作为等周比较对象

考虑第一坐标定义的半球

$$ H=\{x\in\sqrt n S^{n-1}:x_1\le0\}. $$

由假设 $\sigma(A)\ge1/2=\sigma(H)$,球面等周不等式给出

$$ \sigma(A_t)\ge\sigma(H_t). \tag{5.1} $$

$H_t$ 是半球 $H$ 的邻域。直接计算 spherical cap 面积可以完成证明,但更方便的是用 Theorem 3.4.5:若 $X\sim\operatorname{Unif}(\sqrt n S^{n-1})$,则 $X$ 是次高斯随机向量,且 $\|X\|_{\psi_2}\le C$。

由于 $\sigma$ 是球面均匀概率测度,

$$ \sigma(H_t)=\mathbb P\{X\in H_t\}. $$

邻域定义推出

$$ H_t\supset \{x\in\sqrt n S^{n-1}:x_1\le t/\sqrt2\}. \tag{5.2} $$ 查看学习笔记:为什么 (5.2) 成立

于是

$$ \sigma(H_t) \ge \mathbb P\{X_1\le t/\sqrt2\} \ge 1-2\exp(-ct^2). $$

最后一个不等式来自 $\|X_1\|_{\psi_2}\le\|X\|_{\psi_2}\le C$。结合 (5.1),引理得证。

Remark 5.1.7 更剧烈的 blow-up

Lemma 5.1.6 中面积为 $1/2$ 的数值不是本质的。它可以替换为任意常数,甚至可以替换为指数级小的量。Exercise 5.3 会让你验证这一点。

Remark 5.1.8 零一律

刚才看到的 blow-up 现象一开始可能很反直觉。一个指数级小的集合 $A$,为什么只在 $2t$ 的小扰动下就变成指数级大的集合 $A_{2t}$?这里 $t$ 可以远小于球面半径 $\sqrt n$。这正是高维空间的典型现象,类似概率论中的零一律:由许多随机变量共同影响的事件,概率往往接近 $0$ 或 $1$。

5.1.4 Theorem 5.1.3 的证明

不失一般性,假设 $\|f\|_{\mathrm{Lip}}=1$。

查看学习笔记:为什么可归一化 Lipschitz 范数

令 $M$ 是 $f(X)$ 的一个中位数,即

$$ \mathbb P\{f(X)\le M\}\ge\frac12, \qquad \mathbb P\{f(X)\ge M\}\ge\frac12. $$

考虑 sublevel set

$$ A=\{x\in\sqrt n S^{n-1}:f(x)\le M\}. $$

因为 $\mathbb P\{X\in A\}\ge1/2$,Lemma 5.1.6 给出

$$ \mathbb P\{X\in A_t\} \ge 1-2\exp(-ct^2). \tag{5.3} $$

另一方面,我们声称

$$ \mathbb P\{X\in A_t\} \le \mathbb P\{f(X)\le M+t\}. \tag{5.4} $$

确实,若 $X\in A_t$,则存在 $y\in A$ 使得 $\|X-y\|_2\le t$。由 $A$ 的定义,$f(y)\le M$。又因为 $\|f\|_{\mathrm{Lip}}=1$,

$$ f(X)\le f(y)+\|X-y\|_2\le M+t. $$

这就证明了 (5.4)。结合 (5.3) 和 (5.4),得到

$$ \mathbb P\{f(X)\le M+t\} \ge 1-2\exp(-ct^2). $$

对 $-f$ 重复同样论证,可得到 $f(X)\ge M-t$ 的概率下界。

查看学习笔记:补全对 $-f$ 的论证

合并两侧尾界,得到 $|f(X)-M|\le t$ 的类似概率界,从而

$$ \|f(X)-M\|_{\psi_2}\le C. $$

最后把中位数 $M$ 替换为均值 $\mathbb Ef(X)$,这可由中心化和 Lemma 2.7.8 得到。

查看学习笔记:从中位数改为均值

Theorem 5.1.3 得证。Exercise 5.7 会让你反向证明:集中现象与 blow-up 现象本质上等价。

5.2 其他度量测度空间上的集中

我们现在把球面上的集中推广到其他空间。Theorem 5.1.3 的证明依赖两个成分:

  1. 一个等周不等式。
  2. 最小化集合的 blow-up。

这两个成分并非球面独有;许多空间都有类似结构,因此也有类似集中结论。下面介绍两个关键例子:$\mathbb R^n$ 中的 Gaussian concentration,以及 Hamming cube 上的集中;之后再简要提到其他情形。

Remark 5.2.1 均值、中位数、$L^p$ 范数,都可以作为中心

集中意味着均值、中位数和 $L^p$ 范数彼此接近。因此,可以把 $\mathbb Ef(X)$ 替换为中位数(Exercise 5.6);若均值非负,也可以替换为任意 $p\ge1$ 的 $L^p$ 范数,不过常数可能依赖 $p$(Exercise 5.10)。

5.2.1 Gaussian concentration

Theorem 5.1.4 中 $\mathbb R^n$ 上的经典等周不等式不仅对体积成立,也对 $\mathbb R^n$ 上的 Gaussian measure 成立。对 Borel 集 $A\subset\mathbb R^n$,Gaussian measure 定义为

$$ \gamma_n(A) = \mathbb P\{X\in A\} = \frac1{(2\pi)^{n/2}}\int_A e^{-\|x\|_2^2/2}\,dx, $$

其中 $X\sim N(0,I_n)$。

Theorem 5.2.2 Gaussian 等周不等式

设 $\varepsilon>0$。在所有具有给定 Gaussian measure $\gamma_n(A)$ 的集合 $A\subset\mathbb R^n$ 中,half-spaces 使邻域 Gaussian measure $\gamma_n(A_\varepsilon)$ 最小。

用与球面相同的方法,可以推出下面的 Gaussian concentration inequality(见 Exercise 5.8)。

Theorem 5.2.3 Gaussian concentration

设 $X\sim N(0,I_n)$,并设 $f:\mathbb R^n\to\mathbb R$ 是关于 Euclidean 距离的 Lipschitz 函数。则

$$ \|f(X)-\mathbb Ef(X)\|_{\psi_2} \le C\|f\|_{\mathrm{Lip}}. \tag{5.5} $$ 查看学习笔记:由 Gaussian 等周推出集中
Example 5.2.4 两个熟悉特例

(a) 对线性函数 $f$,结论来自 $X\sim N(0,I_n)$ 是次高斯随机向量。

(b) 对 Euclidean 范数 $f(x)=\|x\|_2$,结论来自 norm concentration(Theorem 3.1.1)。

Exercise 5.9 会让你用 Gaussian concentration 证明 $n$ 个 Gaussian 最大值的集中。

5.2.2 Hamming cube

基于等周的集中方法也适用于 Hamming cube

$$ (\{0,1\}^n,d,\mathbb P), $$

其中 $d(x,y)$ 是归一化 Hamming 距离:

$$ d(x,y)=\frac1n|\{i:x_i\ne y_i\}|. $$

测度 $\mathbb P$ 是 cube 上的均匀概率测度:

$$ \mathbb P(A)=\frac{|A|}{2^n}, \qquad A\subset\{0,1\}^n. $$

Theorem 5.2.5 Hamming cube 上的集中

设 $X\sim\operatorname{Unif}\{0,1\}^n$。因此,$X$ 的坐标是独立的 $\operatorname{Ber}(1/2)$ 随机变量。则对任意 $f:\{0,1\}^n\to\mathbb R$,有

$$ \|f(X)-\mathbb Ef(X)\|_{\psi_2} \le \frac{C\|f\|_{\mathrm{Lip}}}{\sqrt n}. \tag{5.6} $$

这一结论来自 Hamming cube 上的等周不等式,其最小化者是 Hamming balls,即关于 Hamming 距离的单点邻域。

5.2.3 对称群

类似结论也适用于对称群 $S_n$,即 $n$ 个符号 $\{1,\ldots,n\}$ 的全部 $n!$ 个排列。把它看成度量测度空间

$$ (S_n,d,\mathbb P), $$

其中 $d(\pi,\rho)$ 是归一化 Hamming 距离:

$$ d(\pi,\rho)=\frac1n|\{i:\pi(i)\ne\rho(i)\}|, $$

而 $\mathbb P$ 是 $S_n$ 上的均匀概率测度:

$$ \mathbb P(A)=\frac{|A|}{n!}. $$

Theorem 5.2.6 对称群上的集中

设 $X\sim\operatorname{Unif}(S_n)$,并设 $f:S_n\to\mathbb R$。则集中不等式 (5.6) 成立。

5.2.4 正 Ricci 曲率 Riemannian 流形

Riemannian 流形提供了许多集中空间的例子。如果你不关注微分几何,可以跳过这一节剩余内容。

紧连通 Riemannian 流形 $(M,g)$ 带有 geodesic distance $d(x,y)$,即连接两点的最短曲线长度。它可看成度量测度空间

$$ (M,d,\mathbb P), $$

其中 $\mathbb P$ 是归一化 Riemannian 体积给出的均匀概率测度。令 $c(M)$ 表示 Ricci curvature tensor 在所有切向量上的下确界。若 $c(M)>0$,则可证明对任意 Lipschitz 函数 $f:M\to\mathbb R$,

$$ \|f(X)-\mathbb Ef(X)\|_{\psi_2} \le \frac{C\|f\|_{\mathrm{Lip}}}{\sqrt{c(M)}}. \tag{5.7} $$

例如 $c(S^{n-1})=n-1$,因此 (5.7) 给出单位球面集中不等式 (5.29) 的另一种证明。

5.2.5 特殊正交群

特殊正交群 $\operatorname{SO}(n)$ 由 $\mathbb R^n$ 中所有旋转组成,等价地说,是所有行列式为 $1$ 的 $n\times n$ 正交矩阵。把它看成度量测度空间

$$ (\operatorname{SO}(n),\|\cdot\|_F,\mathbb P), $$

距离由 Frobenius 范数给出,$\mathbb P$ 是均匀测度。

Theorem 5.2.7 特殊正交群上的集中

设随机正交矩阵 $X\sim\operatorname{Unif}(\operatorname{SO}(n))$,并设 $f:\operatorname{SO}(n)\to\mathbb R$。则集中不等式 (5.6) 成立。

Remark 5.2.8 Haar measure

生成 $X\sim\operatorname{Unif}(\operatorname{SO}(n))$ 的一种方法是:先取 $n\times n$ Gaussian 随机矩阵 $G$,其 entries 独立同分布为 $N(0,1)$;再计算 SVD $G=U\Sigma V^{\mathsf T}$,令 $X=UV^{\mathsf T}$。这样得到的 $X$ 均匀分布在 $\operatorname{SO}(n)$ 上。

$\operatorname{SO}(n)$ 上的均匀概率分布为

$$ \mu(A)=\mathbb P\{X\in A\}. $$

它是唯一的 rotation-invariant probability measure,称为 Haar measure。

查看学习笔记:旋转不变性检查

5.2.6 Grassmannian

Grassmann 流形 $G_{n,m}$ 由 $\mathbb R^n$ 中所有 $m$ 维子空间组成。当 $m=1$ 时,它可与球面 $S^{n-1}$ 识别起来,因此 Grassmannian 上的集中包含球面集中。把 $G_{n,m}$ 看成度量测度空间

$$ (G_{n,m},d,\mathbb P), $$

其中子空间 $E$ 与 $F$ 的距离定义为

$$ d(E,F)=\|P_E-P_F\|, $$

$P_E$ 和 $P_F$ 是对应正交投影。概率测度 $\mathbb P$ 仍是均匀的 Haar probability measure。随机子空间 $E\sim\operatorname{Unif}(G_{n,m})$ 可通过 $n\times m$ Gaussian 随机矩阵 $G$ 的像空间构造。

查看学习笔记:Grassmannian 构造的旋转不变性

Theorem 5.2.9 Grassmannian 上的集中

设随机子空间 $X\sim\operatorname{Unif}(G_{n,m})$,并设 $f:G_{n,m}\to\mathbb R$。则集中不等式 (5.6) 成立。

5.2.7 连续 cube 与 Euclidean ball

对单位 Euclidean cube $[0,1]^n$ 和 Euclidean ball $\sqrt n B_2^n$,也有类似的集中不等式。这可通过把 Gaussian measure push forward 到 cube 或 ball 上的均匀测度来证明。证明留给 Exercises 5.12、5.13。

Theorem 5.2.10 连续 cube 与 ball 上的集中

设 $T$ 是 cube $[0,1]^n$ 或 ball $\sqrt n B_2^n$。若 $X\sim\operatorname{Unif}(T)$,且 $f:T\to\mathbb R$ 是关于 Euclidean 距离的 Lipschitz 函数,则集中不等式 (5.5) 成立。

5.2.8 形如 $e^{-U(x)}$ 的密度

前一节的 push-forward 方法可用于 $\mathbb R^n$ 上许多其他分布。设随机向量 $X$ 有密度

$$ p(x)=e^{-U(x)} $$

其中 $U:\mathbb R^n\to\mathbb R$。例如,若 $X\sim N(0,I_n)$,则正态密度对应 $U(x)=\|x\|_2^2/2+c$,Gaussian concentration 成立。

如果一般函数 $U$ 的曲率至少像 $\|x\|_2^2$ 一样强,就应期待至少有 Gaussian concentration。$U$ 的曲率由 Hessian $\operatorname{Hess}U(x)$ 度量。

Theorem 5.2.11 强凸势能下的集中

设随机向量 $X$ 在 $\mathbb R^n$ 中的密度为 $p(x)=e^{-U(x)}$。若存在 $\kappa>0$,使得对所有 $x\in\mathbb R^n$ 都有

$$ \operatorname{Hess}U(x)\succeq \kappa I_n, $$

则任意 Lipschitz 函数 $f:\mathbb R^n\to\mathbb R$ 满足

$$ \|f(X)-\mathbb Ef(X)\|_{\psi_2} \le \frac{C\|f\|_{\mathrm{Lip}}}{\sqrt\kappa}. $$

5.2.9 独立有界坐标的随机向量

还有一个重要的部分推广:设 $X=(X_1,\ldots,X_n)$ 的坐标独立,且坐标分布任意但有界。通过缩放,可假设 $|X_i|\le1$。

Theorem 5.2.12 Talagrand concentration inequality

设 $X=(X_1,\ldots,X_n)$ 的坐标独立并满足 $|X_i|\le1$ a.s.。那么,对任意凸 Lipschitz 函数 $f:[-1,1]^n\to\mathbb R$,集中不等式 (5.5) 成立。

特别地,Talagrand concentration inequality 适用于 $\mathbb R^n$ 上任意范数。本书不证明该结果。

5.3 应用:Johnson-Lindenstrauss 引理

假设有 $N$ 个数据点位于 $\mathbb R^n$,而维度 $n$ 很大。能否在不严重损失数据几何结构的前提下降维?最简单的方法是把数据点投影到低维子空间

$$ E\subset\mathbb R^n, \qquad \dim(E)=m\ll n. $$

Figure 5.2 展示了这一想法。关键问题是:怎样选择子空间 $E$?维度 $m$ 可以多小?

Johnson-Lindenstrauss 随机投影降维示意图
图 5.2 Johnson-Lindenstrauss 引理通过随机投影到低维子空间来降低数据维度。

Johnson-Lindenstrauss 引理说明,只要把 $E$ 选为维度

$$ m\asymp \log N $$

的随机子空间,就能很好地保留数据几何。

Theorem 5.3.1 Johnson-Lindenstrauss 引理

设 $\mathcal X$ 是 $\mathbb R^n$ 中含有 $N$ 个点的集合,$\varepsilon>0$,并假设

$$ m\ge C\varepsilon^{-2}\log N. $$

令 $P$ 为 $\mathbb R^n$ 到随机 $m$ 维子空间 $E\sim\operatorname{Unif}(G_{n,m})$ 上的正交投影。则以至少

$$ 1-2\exp(-c\varepsilon^2m) $$

的概率,缩放投影 $Q=\sqrt{n/m}\,P$ 在 $\mathcal X$ 上是近似等距:

$$ (1-\varepsilon)\|x-y\|_2 \le \|Qx-Qy\|_2 \le (1+\varepsilon)\|x-y\|_2 \quad \text{for all }x,y\in\mathcal X. \tag{5.8} $$ 查看学习笔记完整证明

证明基于球面上 Lipschitz 函数的集中。先研究随机投影 $P$ 对固定向量 $x-y$ 的作用,再对所有 $N^2$ 个差向量做 union bound。

Lemma 5.3.2 随机投影

令 $P$ 是 $\mathbb R^n$ 到随机 $m$ 维子空间 $E\sim\operatorname{Unif}(G_{n,m})$ 上的投影。固定任意 $z\in\mathbb R^n$ 和 $\varepsilon>0$。则:

(a)

$$ \bigl(\mathbb E\|Pz\|_2^2\bigr)^{1/2} = \sqrt{\frac mn}\|z\|_2. $$

(b) 以至少 $1-2\exp(-c\varepsilon^2m)$ 的概率,

$$ (1-\varepsilon)\sqrt{\frac mn}\|z\|_2 \le \|Pz\|_2 \le (1+\varepsilon)\sqrt{\frac mn}\|z\|_2. $$ 查看学习笔记完整证明
Proof 固定向量上的随机投影

不失一般性设 $\|z\|_2=1$。换一个视角:随机 $m$ 维子空间 $E$ 可由固定坐标子空间 $\mathbb R^m$ 随机旋转得到。等价地,可以固定 $E=\mathbb R^m$,随机旋转向量 $z$;此时 $z$ 均匀分布在 $S^{n-1}$ 上。由旋转不变性,$\|Pz\|_2$ 的分布不变。

查看学习笔记:随机子空间和随机向量视角等价

(a) 此时 $P$ 投影到前 $m$ 个坐标,因此

$$ \mathbb E\|Pz\|_2^2 = \mathbb E\sum_{i=1}^m z_i^2 = m\mathbb Ez_1^2 = \frac mn. $$

最后一步来自 $\sum_{i=1}^n z_i^2=1$ 且各坐标同分布。

(b) 函数 $x\mapsto\|Px\|_2$ 在 $S^{n-1}$ 上 Lipschitz 范数至多为 $1$。

查看学习笔记:投影范数函数是 1-Lipschitz

由单位球面集中不等式 (5.30),

$$ \mathbb P\left\{ \left|\|Px\|_2-\sqrt{m/n}\right|\ge t \right\} \le 2\exp(-cnt^2). $$

这里用 Remark 5.2.1 把 $\mathbb E\|Px\|_2$ 换成 $(\mathbb E\|Px\|_2^2)^{1/2}$。取 $t=\varepsilon\sqrt{m/n}$,结论得证。

Proof of Theorem 5.3.1 对所有差向量做 union bound

考虑差集

$$ \mathcal X-\mathcal X=\{x-y:x,y\in\mathcal X\}. $$

我们希望以所需概率证明对所有 $z\in\mathcal X-\mathcal X$,

$$ (1-\varepsilon)\|z\|_2 \le \|Qz\|_2 \le (1+\varepsilon)\|z\|_2. $$

由于 $Q=\sqrt{n/m}\,P$,这等价于

$$ (1-\varepsilon)\sqrt{\frac mn}\|z\|_2 \le \|Pz\|_2 \le (1+\varepsilon)\sqrt{\frac mn}\|z\|_2. \tag{5.10} $$

对固定 $z$,Lemma 5.3.2 说明 (5.10) 失败的概率至多 $2\exp(-c\varepsilon^2m)$。对 $|\mathcal X-\mathcal X|\le N^2$ 个差向量做 union bound,可知 (5.10) 对所有差向量同时成立的概率至少为

$$ 1-N^2\cdot2\exp(-c\varepsilon^2m). $$

若 $m\ge C\varepsilon^{-2}\log N$,上式至少为 $1-2\exp(-c\varepsilon^2m/2)$。调整常数即得定理。

Remark 5.3.3 非自适应、与原维度无关

Johnson-Lindenstrauss 引理的一个突出特征是:降维映射是非自适应的,不依赖数据本身。注意,数据的 ambient dimension $n$ 在维度条件里没有出现。后续 Sections 9.2.4 和 9.7 会发展更高级的 Johnson-Lindenstrauss 版本。

Remark 5.3.4 最优性

Johnson-Lindenstrauss 引理把维度降到 $O(\log N)$。还能不能更低,例如 $o(\log N)$?Exercise 5.15 会说明不能:即使允许非线性映射,$\log N$ 量级也是最优的。

Exercise 5.14 会让你证明 Johnson-Lindenstrauss 引理的次高斯矩阵版本。

5.4 矩阵 Bernstein 不等式

这里,我们把独立随机变量和 $\sum X_i$ 的集中不等式推广到独立随机矩阵和。矩阵 Bernstein 不等式是 Theorem 2.9.5 的矩阵版本:把随机变量 $X_i$ 换成随机矩阵,把绝对值 $|\cdot|$ 换成算子范数 $\|\cdot\|$。注意,每个随机矩阵 $X_i$ 内部的 entries、rows 或 columns 不需要独立,这是非常一般的假设。

Theorem 5.4.1 Matrix Bernstein inequality

设 $X_1,\ldots,X_N$ 是独立、均值为零的 $n\times n$ 对称随机矩阵,并且对所有 $i$,

$$ \|X_i\|\le K \quad\text{a.s.} $$

则对所有 $t\ge0$,

$$ \mathbb P\left\{ \left\|\sum_{i=1}^N X_i\right\|\ge t \right\} \le 2n\exp\left( -\frac{t^2/2}{\sigma^2+Kt/3} \right), $$

其中

$$ \sigma^2 = \left\|\sum_{i=1}^N \mathbb E X_i^2\right\| $$

是矩阵和的 variance matrix 的算子范数。

等价地,右侧可写成次高斯与次指数混合尾:

$$ \mathbb P\left\{ \left\|\sum_{i=1}^N X_i\right\|\ge t \right\} \le 2n\exp\left[ -c\min\left(\frac{t^2}{\sigma^2},\frac tK\right) \right]. $$ 查看学习笔记完整证明

证明思路很简单:重复第 2.9 节的 MGF 论证,只是把标量替换为矩阵。大部分步骤都可工作,唯一的主要挑战是矩阵乘法不可交换。因此先准备矩阵函数演算。

5.4.1 矩阵函数演算

对 $n\times n$ 对称矩阵 $X$,取逆、平方等运算只作用在特征值上,特征向量保持不变。若谱分解为

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

$$ X^{-1}=\sum_i\lambda_i^{-1}u_iu_i^{\mathsf T}, \qquad X^2=\sum_i\lambda_i^2u_iu_i^{\mathsf T}, \qquad 2I_n-5X^3=\sum_i(2-5\lambda_i^3)u_iu_i^{\mathsf T}. \tag{5.11} $$

查看学习笔记:矩阵 Taylor 展开与 (5.11)

Definition 5.4.2 矩阵函数

设 $f:\mathbb R\to\mathbb R$,且 $X$ 是 $n\times n$ 对称矩阵,谱分解为

$$ X=\sum_{i=1}^n\lambda_i u_i u_i^{\mathsf T}. $$

定义

$$ f(X)=\sum_{i=1}^n f(\lambda_i)u_i u_i^{\mathsf T}. $$
Definition 5.4.3 Loewner order

若 $X$ 是对称半正定矩阵,记 $X\succeq0$。进一步,若 $X-Y\succeq0$,则记 $X\succeq Y$ 或 $Y\preceq X$。

这是偏序而不是全序,因为有些矩阵之间既没有 $X\succeq Y$,也没有 $Y\succeq X$。

Proposition 5.4.4 Loewner order 的简单性质

(a) 特征值单调性:$X\preceq Y$ 推出 $\lambda_i(X)\le\lambda_i(Y)$。

(b) Trace 单调性:若 $f:\mathbb R\to\mathbb R$ 弱增,则

$$ X\preceq Y \quad\Longrightarrow\quad \operatorname{tr}f(X)\le\operatorname{tr}f(Y). $$

(c) 算子范数:

$$ \|X\|\le a \quad\Longleftrightarrow\quad -aI_n\preceq X\preceq aI_n, \qquad a\ge0. \tag{5.12} $$

(d) 标量不等式升级为矩阵不等式:若对所有 $|x|\le a$ 都有 $f(x)\le g(x)$,则对所有 $\|X\|\le a$ 都有 $f(X)\preceq g(X)$。

查看学习笔记完整证明
Remark 5.4.5 算子范数是矩阵版“绝对值”

(5.12) 是标量事实 $|x|\le a\iff -a\le x\le a$ 的矩阵版本。这解释了为什么矩阵 Bernstein 不等式中自然出现算子范数。

Remark 5.4.6 矩阵单调性

能否把 Proposition 5.4.4(b) 的 trace monotonicity 升级为 matrix monotonicity,即

$$ X\preceq Y \quad\Longrightarrow\quad f(X)\preceq f(Y) \tag{5.13} $$

对所有弱增函数 $f$ 成立?若 $X$ 与 $Y$ 交换,答案是肯定的;一般情形是否定的(Exercise 5.17)。不过某些函数确实是矩阵单调的,例如 $1/x$ 和 $\log x$:

$$ 0\preceq X\preceq Y \quad\Longrightarrow\quad X^{-1}\succeq Y^{-1}\succeq0, \qquad \log X\preceq\log Y. $$ 查看学习笔记:$1/x$ 和 $\log x$ 的矩阵单调性

5.4.2 Trace inequalities

到目前为止,把标量概念推广到矩阵还比较顺利。但这并不总是成立。矩阵不可交换会使许多标量恒等式失效。例如标量恒等式 $e^{x+y}=e^xe^y$ 对矩阵一般不成立;Exercise 5.19 会让你找出 $n\times n$ 对称矩阵 $X,Y$,使得

$$ e^{X+Y}\ne e^Xe^Y. $$

这很麻烦,因为 $e^{x+y}=e^xe^y$ 正是标量 MGF 方法中分解和的关键。幸运的是,有一些 trace inequalities 可以替代缺失的恒等式。

Theorem 5.4.7 Golden-Thompson inequality

对任意 $n\times n$ 对称矩阵 $A,B$,

$$ \operatorname{tr}(e^{A+B}) \le \operatorname{tr}(e^Ae^B). $$

不过,Golden-Thompson 不等式不能直接推广到三个或更多矩阵。

Theorem 5.4.8 Lieb inequality

设 $H$ 是 $n\times n$ 对称矩阵。定义正定矩阵上的函数

$$ f(X)=\operatorname{tr}\exp(H+\log X). $$

则 $f$ 在正定 $n\times n$ 对称矩阵空间上是凹函数。

如果 $X$ 是随机矩阵,Lieb 与 Jensen 不等式推出

$$ \mathbb Ef(X)\le f(\mathbb EX). $$

查看学习笔记:为什么随机矩阵可用 Jensen

令 $X=e^Z$,得到下面的形式。

Lemma 5.4.9 随机矩阵版 Lieb inequality

设 $H$ 是固定的 $n\times n$ 对称矩阵,$Z$ 是随机 $n\times n$ 对称矩阵。则

$$ \mathbb E\operatorname{tr}\exp(H+Z) \le \operatorname{tr}\exp(H+\log\mathbb Ee^Z). $$

5.4.3 矩阵 Bernstein 不等式的证明

$$ S=\sum_{i=1}^N X_i. $$

为了控制 $\|S\|$,需要分别控制 $S$ 的最大特征值和最小特征值。记

$$ \lambda_{\max}(S)=\max_i\lambda_i(S). $$

$$ \|S\| = \max_i|\lambda_i(S)| = \max(\lambda_{\max}(S),\lambda_{\max}(-S)). \tag{5.14} $$

固定 $\lambda\ge0$。由 Markov 不等式,

$$ \mathbb P\{\lambda_{\max}(S)\ge t\} \le e^{-\lambda t}\mathbb E e^{\lambda\lambda_{\max}(S)}. \tag{5.15} $$

而 $e^{\lambda S}$ 的特征值是 $e^{\lambda\lambda_i(S)}$,所以

$$ \mathbb E e^{\lambda\lambda_{\max}(S)} \le \mathbb E\operatorname{tr}e^{\lambda S}. $$

用 Lemma 5.4.9 逐个分离 $\lambda X_N,\lambda X_{N-1},\ldots,\lambda X_1$,得到

$$ \mathbb E e^{\lambda\lambda_{\max}(S)} \le \operatorname{tr}\exp\left[ \sum_{i=1}^N\log\mathbb Ee^{\lambda X_i} \right]. \tag{5.16} $$

接下来需要控制单个矩阵的 MGF。

Lemma 5.4.10 矩阵 MGF 控制

设 $X$ 是均值为零的 $n\times n$ 对称随机矩阵,并满足 $\|X\|\le K$ a.s.。则当 $|\lambda|<3/K$ 时,

$$ \mathbb E\exp(\lambda X) \preceq \exp(g(\lambda)\mathbb EX^2), \qquad g(\lambda)=\frac{\lambda^2/2}{1-|\lambda|K/3}. $$ 查看学习笔记完整证明

把 Lemma 5.4.10 代入 (5.16),并使用 $\log x$ 的矩阵单调性与 trace monotonicity,得到

$$ \mathbb E e^{\lambda\lambda_{\max}(S)} \le \operatorname{tr}\exp(g(\lambda)Z), \qquad Z=\sum_{i=1}^N\mathbb EX_i^2. $$

由于 $Z\succeq0$,

$$ \operatorname{tr}\exp(g(\lambda)Z) \le n\exp(g(\lambda)\|Z\|) = n\exp(g(\lambda)\sigma^2). $$

代回 (5.15),得到

$$ \mathbb P\{\lambda_{\max}(S)\ge t\} \le n\exp[-\lambda t+g(\lambda)\sigma^2]. $$

$$ \lambda=\frac{t}{\sigma^2+Kt/3}, $$

可化简为

$$ \mathbb P\{\lambda_{\max}(S)\ge t\} \le n\exp\left(-\frac{t^2/2}{\sigma^2+Kt/3}\right). $$

对 $-S$ 重复论证,并用 (5.14) 合并,即得 Theorem 5.4.1。

查看学习笔记:补全 $-S$ 与双尾合并

Remark 5.4.11 矩阵 Bernstein:期望形式

矩阵 Bernstein 不等式给出高概率界。用尾积分公式可推出更简单但信息更少的期望界:

$$ \mathbb E\left\|\sum_{i=1}^N X_i\right\| \lesssim \left\|\sum_{i=1}^N\mathbb EX_i^2\right\|^{1/2} \sqrt{\log(2n)} + K\log(2n). \tag{5.17} $$ 查看学习笔记:由尾界推出期望界
Remark 5.4.12 对数代价

和标量情形相比,高维矩阵升级只多出一个对数因子。这是高维中的小代价,并且在一般情形下基本不可去掉;Exercise 5.28 给出例子。

5.4.4 矩阵 Hoeffding 与 Khintchine 不等式

Theorem 5.4.13 Matrix Hoeffding inequality

设 $\varepsilon_1,\ldots,\varepsilon_N$ 是独立 Rademacher 随机变量,$A_1,\ldots,A_N$ 是固定的对称 $n\times n$ 矩阵。则对任意 $t\ge0$,

$$ \mathbb P\left\{ \left\|\sum_{i=1}^N\varepsilon_iA_i\right\|\ge t \right\} \le 2n\exp\left(-\frac{t^2}{2\sigma^2}\right), $$

其中

$$ \sigma^2=\left\|\sum_{i=1}^N A_i^2\right\|. $$ 查看学习笔记完整证明
Theorem 5.4.14 Matrix Khintchine inequality

在 Theorem 5.4.13 的设定下,对每个 $p\in[1,\infty)$,

$$ \left( \mathbb E\left\|\sum_{i=1}^N\varepsilon_iA_i\right\|^p \right)^{1/p} \le C\sqrt{p+\log n} \left\|\sum_{i=1}^N A_i^2\right\|^{1/2}. $$ 查看学习笔记完整证明
Remark 5.4.15 非对称与矩形矩阵

矩阵集中不等式可用 Hermitian dilation 扩展到矩形矩阵。把每个 $m\times n$ 矩阵 $X_i$ 替换为对称块矩阵

$$ \mathcal D(X_i)= \begin{bmatrix} 0 & X_i\\ X_i^{\mathsf T} & 0 \end{bmatrix}, $$

再应用对称矩阵的集中不等式。Exercises 5.23、5.24 会给出矩形矩阵 Bernstein 与 Khintchine 版本。

5.5 应用:稀疏网络中的社群检测

第 4.5 节分析了 stochastic block model $G(n,p,q)$ 上的谱聚类,并说明当期望平均度 $\gtrsim\sqrt n$ 时算法有效。现在借助矩阵 Bernstein 不等式,我们将证明谱聚类可用于更稀疏的网络,期望平均度低至 $O(\log n)$。

Theorem 5.5.1 稀疏 stochastic block model 的谱聚类

设 $G\sim G(n,p,q)$,其中 $p=a/n$、$q=b/n$,且 $b<a<3b$。若

$$ (a-b)^2\ge Ca\log n, $$

则以至少 $0.99$ 的概率,谱聚类算法(Section 4.5.5)能以 $99\%$ 的准确率识别 $G$ 的两个社群,也就是误分类顶点数至多为 $0.01n$。

查看学习笔记完整证明

证明沿用 Section 4.5 的论证,只是把误差上界换成更锐利的矩阵 Bernstein 界。设邻接矩阵为 $A$,分解为

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

随机部分 $R$ 可写成独立、均值为零的对称矩阵和:

$$ R=\sum_{i\le j}Z_{ij}, $$

其中 $Z_{ij}$ 只保留 $(i,j)$ 与 $(j,i)$ 两个位置的随机噪声。因为 $A_{ij}\in\{0,1\}$,有 $|R_{ij}|\le1$,从而 $\|Z_{ij}\|\le1$。

查看学习笔记:为什么 $\|Z_{ij}\|\le1$

用矩阵 Bernstein 的期望形式和 Markov 不等式,可以以至少 $0.99$ 的概率得到

$$ \|R\| \lesssim \sigma\sqrt{\log n}+\log n, \qquad \sigma^2= \left\|\mathbb E\sum_{i\le j}Z_{ij}^2\right\|. \tag{5.18} $$

计算可得 $\sigma^2\le np=a$,因此

$$ \|R\| \lesssim \sqrt{a\log n}+\log n \lesssim \sqrt{a\log n}. \tag{5.19} $$

最后一步使用假设推出 $a\gtrsim\log n$。

查看学习笔记:为什么假设推出 $a\gtrsim\log n$

再对 $D$ 和 $A$ 应用 Davis-Kahan。第 4.5 节已经算出 $D$ 的第二特征值与其余谱的分离为

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

这里用到了 $a\le3b$。由 Davis-Kahan,

$$ \|\bar u_2(D)-\theta\bar u_2(A)\|_2 \le \frac{2\|R\|}{\delta} \le \frac{C_1\sqrt{a\log n}}{a-b} < \frac1{10}, $$

其中 $\theta\in\{-1,1\}$。若定理假设中的常数 $C$ 取得足够大,最后不等式成立。把单位向量还原为 $\pm1$ 系数向量,可得至少 $99\%$ 的顶点符号正确。

查看学习笔记:特征向量误差如何推出 99% 符号正确

Remark 5.5.2 稀疏性

Theorem 5.5.1 非平凡时,最稀疏图的期望平均度满足

$$ \frac{n(p+q)}2=\frac{a+b}{2}\asymp\log n. $$

这比第 4 章得到的 $O(\sqrt n)$ 稀疏度大幅改进。Exercise 5.25 会让你处理没有 loops 的 stochastic block model。

5.6 应用:一般分布的协方差估计

第 4.7 节中,我们学习了如何用 $O(n)$ 个样本估计 $\mathbb R^n$ 中次高斯分布的协方差矩阵。现在去掉次高斯假设,使方法适用于更广泛的分布,甚至离散分布;代价只是一个对数 oversampling 因子。

与第 4.7 节一样,我们用样本版本估计 second moment matrix

$$ \Sigma=\mathbb EXX^{\mathsf T}, \qquad \Sigma_m=\frac1m\sum_{i=1}^m X_iX_i^{\mathsf T}. $$

如果 $X$ 均值为零,则 $\Sigma$ 是 $X$ 的协方差矩阵,$\Sigma_m$ 是样本协方差矩阵。

Theorem 5.6.1 一般协方差估计

设 $X$ 是 $\mathbb R^n$ 中的随机向量,$n\ge2$。假设存在 $K\ge1$,使得

$$ \|X\|_2 \le K(\mathbb E\|X\|_2^2)^{1/2} \quad \text{a.s.} \tag{5.20} $$

那么对任意正整数 $m$,

$$ \mathbb E\|\Sigma_m-\Sigma\| \le C\left( \sqrt{\frac{K^2 n\log n}{m}} + \frac{K^2 n\log n}{m} \right)\|\Sigma\|. $$ 查看学习笔记完整证明

证明中用 Proposition 3.2.1(b) 得到 $\mathbb E\|X\|_2^2=\operatorname{tr}(\Sigma)$,所以 (5.20) 等价于

$$ \|X\|_2^2\le K^2\operatorname{tr}(\Sigma) \quad \text{a.s.} \tag{5.21} $$

对 i.i.d. 矩阵和 $\sum_i(X_iX_i^{\mathsf T}-\Sigma)$ 应用矩阵 Bernstein 的期望形式 (5.17),得到

$$ \mathbb E\|\Sigma_m-\Sigma\| \lesssim \frac1m(\sigma\sqrt{\log n}+M\log n), \tag{5.22} $$

其中

$$ \sigma^2 = m\left\|\mathbb E(XX^{\mathsf T}-\Sigma)^2\right\|, \qquad \|XX^{\mathsf T}-\Sigma\|\le M. $$

先估计 $\sigma^2$。展开平方:

$$ \mathbb E(XX^{\mathsf T}-\Sigma)^2 = \mathbb E(XX^{\mathsf T})^2-\Sigma^2 \preceq \mathbb E(XX^{\mathsf T})^2. \tag{5.23} $$

又由 (5.21),

$$ (XX^{\mathsf T})^2 = \|X\|_2^2XX^{\mathsf T} \preceq K^2\operatorname{tr}(\Sigma)\,XX^{\mathsf T}. $$

取期望得

$$ \mathbb E(XX^{\mathsf T})^2 \preceq K^2\operatorname{tr}(\Sigma)\Sigma, $$

所以

$$ \sigma^2 \le K^2m\operatorname{tr}(\Sigma)\|\Sigma\|. $$

再估计 $M$:

$$ \|XX^{\mathsf T}-\Sigma\| \le \|X\|_2^2+\|\Sigma\| \le K^2\operatorname{tr}(\Sigma)+\|\Sigma\| \le 2K^2\operatorname{tr}(\Sigma) =M. $$

把这些代入 (5.22),再用 $\operatorname{tr}(\Sigma)\le n\|\Sigma\|$ 即得 Theorem 5.6.1。

Remark 5.6.2 样本复杂度

Theorem 5.6.1 说明,对任意 $\varepsilon\in(0,1)$,只要

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

就能得到

$$ \mathbb E\|\Sigma_m-\Sigma\|\le\varepsilon\|\Sigma\|. \tag{5.24} $$

相比次高斯分布的 $m\asymp\varepsilon^{-2}n$,去掉次高斯假设只付出了一个小的对数 oversampling 因子。

Remark 5.6.3 低维分布与有效秩

证明末尾用了粗糙界 $\operatorname{tr}(\Sigma)\le n\|\Sigma\|$。若改用有效秩

$$ r(\Sigma)=\frac{\operatorname{tr}(\Sigma)}{\|\Sigma\|}, \tag{5.26} $$

则可得到更精细的界

$$ \mathbb E\|\Sigma_m-\Sigma\| \le C\left( \sqrt{\frac{K^2 r\log n}{m}} + \frac{K^2 r\log n}{m} \right)\|\Sigma\|. \tag{5.27} $$

因此,$m\asymp\varepsilon^{-2}r\log n$ 个样本足以估计协方差。由于 $r\le n$,它不会比 (5.25) 更差。

查看学习笔记:为什么 $r(\Sigma)\le n$
Remark 5.6.4 矩阵的有效秩与稳定秩

若 $\Sigma\succeq0$,则

$$ r(\Sigma) = \frac{\sum_{i=1}^n\lambda_i(\Sigma)} {\max_i\lambda_i(\Sigma)}. $$

它总是被实际 rank 控制,并且对近似低秩矩阵可以远小于实际 rank。相关概念是任意矩阵 $A$ 的 stable rank:

$$ s(A) = \frac{\|A\|_F^2}{\|A\|^2} = \frac{\sum_i s_i(A)^2}{\max_i s_i(A)^2} = r(A^{\mathsf T}A) = r(AA^{\mathsf T}). $$ 查看学习笔记:有效秩练习
Remark 5.6.5 高概率保证

上面给出的是期望界,但同一论证也给出高概率版本:对所有 $u\ge0$,以至少 $1-2e^{-u}$ 的概率,

$$ \|\Sigma_m-\Sigma\| \le C\left( \sqrt{\frac{K^2 r(\log n+u)}{m}} + \frac{K^2 r(\log n+u)}{m} \right)\|\Sigma\|, \tag{5.28} $$

其中 $r=\operatorname{tr}(\Sigma)/\|\Sigma\|\le n$。Exercise 5.26 会让你证明它。

查看学习笔记完整证明
Remark 5.6.6 有界性假设

(5.20) 的有界性假设看似强,但一般不能去掉:若 $X$ 各向同性但以很高概率为零,则样本很可能全是零,使协方差估计不可能。Exercise 5.27 会让你形式化这个论证。实践中通常用截断来保证这一条件,即丢弃少量范数最大的样本。

查看学习笔记完整证明

5.7 Notes

关于集中不等式,可参见本章 bibliographic notes 中列出的专著、教程和参考文献。第 5.1 节的等周方法源于 P. Levy。V. Milman 意识到 Levy 方法的威力后,推动了 concentration of measure principle 的广泛发展。

本章为了简洁省略了许多重要方法,包括 bounded differences、martingales、semigroups、transport、Poincare 与 log-Sobolev inequalities、hypercontractivity、Stein 方法和 Talagrand inequalities。Gaussian 等周不等式、Hamming cube 与对称群上的集中、正曲率 Riemannian 流形上的集中、连续 cube/ball 上的集中、强凸密度下的集中,以及 Talagrand concentration inequality,都有成熟文献证明。

Johnson-Lindenstrauss 引理的原始版本来自 Johnson 与 Lindenstrauss;条件 $m\gtrsim\varepsilon^{-2}\log N$ 是最优的。矩阵集中不等式部分沿用 Ahlswede-Winter 与 Tropp 发展出的思路;Golden-Thompson、Lieb inequality 和矩阵 Bernstein 的大量变体可在相关矩阵分析和随机矩阵文献中找到。矩阵 Bernstein 中出现的维度因子 $n$ 在某些版本中可被 intrinsic dimension 或 effective rank 替代。

本章还省略了 bounded differences inequality(McDiarmid inequality),它不仅适用于和,也适用于独立随机变量的一般函数,是 Hoeffding 不等式的推广。

Theorem 5.7.1 Bounded differences inequality

设 $X=(X_1,\ldots,X_N)$ 的坐标独立,$f$ 是可测函数。假设改变第 $i$ 个坐标最多使 $f$ 的值改变 $c_i>0$,即对所有 $i$ 和所有可能输入都有

$$ |f(x_1,\ldots,x_i,\ldots,x_N) - f(x_1,\ldots,x_i',\ldots,x_N)| \le c_i. $$

则对任意 $t>0$,

$$ \mathbb P\{f(X)-\mathbb Ef(X)\ge t\} \le \exp\left( -\frac{2t^2}{\sum_{i=1}^N c_i^2} \right). $$ 查看学习笔记完整证明

Exercises

下面保留原书习题内容,并为每题加入学习笔记证明入口。

Exercise 5.1连续、可微与 Lipschitz 函数

(a) 证明每个 Lipschitz 函数都是一致连续的。(b) 证明每个可微函数 $f:\mathbb R^n\to\mathbb R$ 若满足 $\sup_x\|\nabla f(x)\|_2<\infty$,则是 Lipschitz,并且 $\|f\|_{\mathrm{Lip}}\le\sup_x\|\nabla f(x)\|_2$。(c) 找一个非 Lipschitz 但一致连续的 $f:[-1,1]\to\mathbb R$。(d) 找一个不可微但 Lipschitz 的 $f:[-1,1]\to\mathbb R$。

查看学习笔记完整证明
Exercise 5.2验证 Example 5.1.2

检查 Example 5.1.2 中所有关于线性泛函、线性算子和范数函数的 Lipschitz 范数断言。

查看学习笔记完整证明
Exercise 5.3指数小集合的 blow-up

设 $A\subset\sqrt n S^{n-1}$ 且 $\sigma(A)>2\exp(-cs^2)$。(a) 证明 $\sigma(A_s)>1/2$。(b) 推出对任意 $t\ge s$,$\sigma(A_{2t})\ge1-2\exp(-ct^2)$。

查看学习笔记完整证明
Exercise 5.4Geodesic metric 下的集中

Theorem 5.1.3 是对球面上的 Euclidean 距离表述的。证明它对 geodesic metric 也成立,其中 geodesic metric 是两点之间最短球面弧长。

查看学习笔记完整证明
Exercise 5.5单位球面上的集中

由 Theorem 5.1.3 推出:若 $X\sim\operatorname{Unif}(S^{n-1})$,则单位球面上的 Lipschitz 函数满足

$$\|f(X)-\mathbb Ef(X)\|_{\psi_2}\le\frac{C\|f\|_{\mathrm{Lip}}}{\sqrt n}.\tag{5.29}$$

等价地,对任意 $t\ge0$,

$$\mathbb P\{|f(X)-\mathbb Ef(X)|\ge t\}\le2\exp\left(-\frac{cnt^2}{\|f\|_{\mathrm{Lip}}^2}\right).\tag{5.30}$$查看学习笔记完整证明
Exercise 5.6关于均值与中位数集中等价

设 $Z$ 的中位数为 $M$。证明

$$c\|Z-\mathbb EZ\|_{\psi_2}\le\|Z-M\|_{\psi_2}\le C\|Z-\mathbb EZ\|_{\psi_2}.$$查看学习笔记完整证明
Exercise 5.7集中等价于 blow-up

设随机向量 $X$ 取值于度量空间 $(T,d)$,且对所有 Lipschitz 函数 $f:T\to\mathbb R$ 都有 $\|f(X)-\mathbb Ef(X)\|_{\psi_2}\le K\|f\|_{\mathrm{Lip}}$。定义 $\sigma(A)=\mathbb P\{X\in A\}$。证明若 $\sigma(A)\ge1/2$,则对所有 $t\ge0$,

$$\sigma(A_t)\ge1-2\exp(-ct^2/K^2).$$查看学习笔记完整证明
Exercise 5.8Gaussian concentration

由 Gaussian 等周不等式 Theorem 5.2.2 推出 Gaussian concentration inequality Theorem 5.2.3。

查看学习笔记完整证明
Exercise 5.9最大值的集中

(a) 若 $X_1,\ldots,X_n$ 独立同分布为 $N(0,1)$,证明 $\|\max_iX_i-\mathbb E\max_iX_i\|_{\psi_2}\le C$。(b) 更一般地,若 $X_1,\ldots,X_n$ 联合正态,不必独立,证明

$$\|\max_iX_i-\mathbb E\max_iX_i\|_{\psi_2}\le C\max_i\sqrt{\operatorname{Var}(X_i)}.$$查看学习笔记完整证明
Exercise 5.10围绕 $L^p$ 范数的集中

设 $\mathbb EZ\ge0$ 且 $p\ge1$。证明

$$\|Z-\|Z\|_{L^p}\|_{\psi_2}\le C\sqrt p\,\|Z-\mathbb EZ\|_{\psi_2}.$$查看学习笔记完整证明
Exercise 5.11把 Gaussian 推送为均匀测度

令 $\Phi$ 为标准正态分布函数,$Z\sim N(0,I_n)$。检查

$$\phi(Z)=(\Phi(Z_1),\ldots,\Phi(Z_n))\sim\operatorname{Unif}([0,1]^n).$$查看学习笔记完整证明
Exercise 5.12连续 cube 上的集中

按如下方法证明 Theorem 5.2.10 中 $T=[0,1]^n$ 的情形:(a) 用 Exercise 5.11 表示 $X=\phi(Z)$,并用 Gaussian concentration 控制 $f(\phi(Z))$。(b) 证明 $\|\phi\|_{\mathrm{Lip}}$ 被绝对常数控制。

查看学习笔记完整证明
Exercise 5.13连续 ball 上的集中

按类似 Exercise 5.12 的策略证明 Theorem 5.2.10 中 $T=\sqrt n B_2^n$ 的情形:构造把 Gaussian measure 推送到 $\sqrt n B_2^n$ 均匀测度的映射 $\phi$,并检查其 Lipschitz 范数有界。

查看学习笔记完整证明
Exercise 5.14次高斯矩阵版 Johnson-Lindenstrauss

设 $A$ 是 $m\times n$ 随机矩阵,rows 独立、均值为零、次高斯、各向同性。证明 $Q=(1/\sqrt m)A$ 满足 Johnson-Lindenstrauss 引理。说明 $\pm1$ entries 的二值 Johnson-Lindenstrauss 是特例。

查看学习笔记完整证明
Exercise 5.15Johnson-Lindenstrauss 维度最优性

(a) 设 $z_1,\ldots,z_N\in\mathbb R^n$ 满足 $1<\|z_i-z_j\|_2\le2$,证明 $N\le5^n$。(b) 若 $n<\frac12\log N$,构造 $x_1,\ldots,x_N\in\mathbb R^N$,使不存在任何映射 $T:\mathbb R^N\to\mathbb R^n$ 对所有点对保持 $0.99$ 到 $1.01$ 的距离失真。

查看学习笔记完整证明
Exercise 5.16矩阵 Taylor series

验证矩阵函数定义与多项式、幂级数和矩阵指数 Taylor 展开一致。

查看学习笔记完整证明
Exercise 5.17矩阵单调性

(a) 若 $X,Y$ 交换,检查 $X\preceq Y$ 推出 $f(X)\preceq f(Y)$ 对任意弱增函数 $f$ 成立。(b) 给出非交换矩阵的反例。

查看学习笔记完整证明
Exercise 5.18$1/x$ 与 $\log x$ 的矩阵单调性

设 $0\preceq X\preceq Y$ 且 $X$ 可逆。证明 $Y$ 可逆且 $X^{-1}\succeq Y^{-1}\succeq0$;检查积分恒等式

$$\log x=\int_0^\infty\left(\frac1{1+t}-\frac1{x+t}\right)\,dt,$$

并推出 $\log X\preceq\log Y$。

查看学习笔记完整证明
Exercise 5.19矩阵指数

(a) 若对称矩阵 $X,Y$ 交换,证明 $e^{X+Y}=e^Xe^Y$。(b) 找出 $e^{X+Y}\ne e^Xe^Y$ 的例子。

查看学习笔记完整证明
Exercise 5.20矩阵 Bernstein 的期望界

从 Theorem 5.4.1 推出

$$\mathbb E\left\|\sum_{i=1}^N X_i\right\|\lesssim\left\|\sum_{i=1}^N\mathbb EX_i^2\right\|^{1/2}\sqrt{1+\log n}+K(1+\log n).$$查看学习笔记完整证明
Exercise 5.21Matrix Hoeffding inequality

仿照矩阵 Bernstein 的证明,证明 Theorem 5.4.13。

查看学习笔记完整证明
Exercise 5.22Matrix Khintchine inequality

从 Theorem 5.4.13 推出 Theorem 5.4.14。

查看学习笔记完整证明
Exercise 5.23矩形矩阵 Bernstein

设 $X_i$ 是独立、均值为零的 $m\times n$ 随机矩阵,且 $\|X_i\|\le K$ a.s.。证明矩形矩阵 Bernstein 不等式,维度因子为 $2(m+n)$,variance 参数由 $\sum\mathbb EX_i^{\mathsf T}X_i$ 与 $\sum\mathbb EX_iX_i^{\mathsf T}$ 控制。

查看学习笔记完整证明
Exercise 5.24矩形矩阵 Khintchine

对固定 $m\times n$ 矩阵 $A_i$ 和独立 Rademacher 变量 $\varepsilon_i$,证明

$$\left(\mathbb E\left\|\sum_i\varepsilon_iA_i\right\|^p\right)^{1/p}\le C\sigma\sqrt{p+\log(m+n)},$$

其中 $\sigma^2=\left\|\sum_iA_i^{\mathsf T}A_i\right\|+\left\|\sum_iA_iA_i^{\mathsf T}\right\|$。

查看学习笔记完整证明
Exercise 5.25无 loops 的 stochastic block model

修改 Definition 4.5.1,使 stochastic block model 不允许 loops,并证明 Theorem 5.5.1 的相应版本。

查看学习笔记完整证明
Exercise 5.26协方差估计的高概率版本

证明 Remark 5.6.5 中的高概率协方差估计结论。

查看学习笔记完整证明
Exercise 5.27有界性假设不能直接去掉

证明 Theorem 5.6.1 中的有界性假设 (5.20) 一般不能去掉。

查看学习笔记完整证明
Exercise 5.28协方差估计与矩阵 Bernstein 中的对数因子

(a) 构造一个分布,使得若 $m\not\gtrsim n\log n$,则 $\|\Sigma_m-\Sigma\|<\|\Sigma\|$ 以高概率失败。(b) 推出矩阵 Bernstein 期望界中的对数因子一般不能去掉。

查看学习笔记完整证明
Exercise 5.29有效秩

对 $\Sigma\succeq0$,$r(\Sigma)=\operatorname{tr}(\Sigma)/\|\Sigma\|$。(a) 证明 $1\le r(\Sigma)\le\operatorname{rank}(\Sigma)\le n$。(b) 说明这些不等式最优。(c) 说明有效秩是连续的。(d) 若 $X$ 取值于 $k$ 维子空间,证明 $\operatorname{rank}(\Sigma)\le k$,从而 $r(\Sigma)\le k$。

查看学习笔记完整证明
Exercise 5.30从 frames 采样

设 $(u_1,\ldots,u_M)$ 是 $\mathbb R^n$ 中等范数 Parseval frame。证明随机抽取 $m\gtrsim n\log n$ 个 frame elements 后,高概率仍形成近似 frame,并形式化“近似”的含义。

查看学习笔记完整证明
Exercise 5.31一般独立 rows 的随机矩阵

设 $A$ 是 $m\times n$ 随机矩阵,rows 独立各向同性,且 $\|A_i\|_2\le K\sqrt n$ a.s.。证明对 $t\ge1$,以至少 $1-2n^{-ct^2}$ 的概率,

$$\sqrt m-Kt\sqrt{n\log n}\le s_n(A)\le s_1(A)\le\sqrt m+Kt\sqrt{n\log n}.$$查看学习笔记完整证明
Exercise 5.32Matrix sketching

设 $A$ 是 tall $N\times n$ 矩阵,所有 rows 等范数。令 $B$ 为从 $A$ 中有放回均匀抽取 $m=O(n\log n)$ 行形成的矩阵。证明若 $m\ge Cn\log n$,则以至少 $0.9$ 的概率,

$$\max_{i=1,\ldots,n}\left|s_i(A)^2-\frac Nm s_i(B)^2\right|\le0.1s_1(A)^2.$$查看学习笔记完整证明

校对说明

  • concentration without independence 译为“无独立性的集中”,强调本章从独立和转向几何与矩阵工具。
  • isoperimetric inequality 译为“等周不等式”,blow-up 保留英文并解释为“邻域快速膨胀”。
  • matrix BernsteinLoewner orderHaar measureJohnson-Lindenstrauss 保留英文术语,方便与原书和后续章节对照。
  • OCR 中若干箭头、矩阵块和英文断字已按数学语义修正。
学习笔记 Ch.5 无独立性集中
第 5 章学习笔记:无独立性的集中

一句话定位

第 5 章把“集中来自独立和”的视角拓展为“集中来自空间几何与矩阵 MGF”:前半章用等周不等式解释 Lipschitz 函数为何集中,后半章用矩阵 Bernstein 把标量 Bernstein 升级到随机矩阵,并应用到稀疏社群检测和一般协方差估计。

本章导读

第 5 章的核心问题是:当没有坐标独立,或者研究对象本身是矩阵时,集中不等式还能从哪里来?本章给出两条路线。

章节 内容 在主线中的作用
5.1 球面上 Lipschitz 函数集中 用等周不等式替代独立性
5.2 其他度量测度空间 建立 concentration of measure 的空间视角
5.3 Johnson-Lindenstrauss 把球面集中用于随机降维
5.4 矩阵 Bernstein 把标量 MGF 方法升级到矩阵和
5.5 稀疏社群检测 用矩阵 Bernstein 控制稀疏图噪声
5.6 一般协方差估计 用矩阵 Bernstein 去掉次高斯假设
5.7 Notes 与 bounded differences 衔接更广的集中理论

本章读法是:前半章先理解“几何空间让所有 Lipschitz 函数集中”,后半章再理解“矩阵不可交换时如何做 Bernstein”。这两条线在第 5.3、5.5、5.6 节分别落到降维、网络和统计估计应用。

本页使用方式

第 5 章容易让初学者迷失,因为它同时出现几何测度、矩阵分析和数据应用。建议不要按术语堆叠来读,而是先判断自己卡在哪一种“非独立”上。

你现在卡在哪里 先看哪里 读完要形成的判断
不理解没有独立性为什么还能集中 5.1、Lemma 5.1.6、Theorem 5.1.3 集中可以来自等周不等式:大集合的邻域会迅速覆盖空间。
不知道 Lipschitz 条件在做什么 5.1.1、Example 5.1.2、Exercise 5.1-5.2 Lipschitz 是“函数不剧烈振荡”的可计算条件。
Johnson-Lindenstrauss 看起来像魔法 5.3、Lemma 5.3.2 固定差向量由球面集中控制,全部点对由 union bound 控制。
矩阵 Bernstein 推不动 5.4.1-5.4.3 标量 MGF 方法仍是骨架,Lieb inequality 负责处理不可交换性。
社群检测和协方差估计只是应用吗 5.5-5.6 它们展示同一个模板:写成信号/均值 + 随机矩阵噪声,再控制噪声范数。
准备做习题 Exercises 5.1-5.32 证明工作区 先分清题目训练的是几何集中、push-forward、JL、矩阵函数、矩阵集中还是应用建模。

本章主线

第 5 章的主线是从“独立和集中”升级到“空间与矩阵集中”。关键转折有两个:一是用等周不等式替代独立性;二是用 trace/Lieb 工具替代标量指数函数的可交换性。

推进层 要解决的问题 关键转折 后续用途
几何集中 没有独立坐标时为什么还会集中? 大集合的邻域快速 blow-up 球面、Gaussian、Hamming cube、群和流形
Lipschitz 观测 哪些函数能从空间集中获益? 函数只要不剧烈拉伸距离,就继承空间集中 随机投影、Gaussian 最大值、范数集中
随机降维 如何保持有限点集的所有距离? 固定差向量集中 + 对 $N^2$ 点对 union bound Johnson-Lindenstrauss 与 sketching
矩阵集中 随机矩阵和的范数如何集中? 用 Loewner order、trace inequalities 和 Lieb inequality 修复 MGF 方法 稀疏 SBM、协方差估计、矩阵 Hoeffding/Khintchine
统计应用 如何在更弱假设下估计结构? 控制噪声矩阵范数,再转成谱扰动或协方差误差 稀疏网络、一般分布 covariance、后续随机过程

本章学习路线

先抓住一个问题
集中不一定来自独立性。

第 5 章的前半部分说:高维空间的几何本身会迫使 Lipschitz 函数集中。后半部分说:即使对象是不可交换的矩阵,只要用对 trace 和谱序友好的 MGF 方法,仍能得到 Bernstein 型尾界。

初学者先抓三件事
  1. 等周不等式怎样推出 blow-up。
  2. 固定向量集中怎样推出 Johnson-Lindenstrauss。
  3. 标量 Bernstein 怎样变成矩阵 Bernstein。
1

先把集中读成几何事实

球面上,如果一个集合有一半面积,它的 $t$-邻域会覆盖几乎整个球面。Lipschitz 函数的 sublevel set 正好把这个几何事实转成函数值集中。

isoperimetryblow-upmedian
这一层要会问 哪个集合有一半测度?它的邻域对应函数值的哪一侧尾部?
2

再把几何集中用于随机投影

随机投影先固定一个差向量 $z$,用球面集中控制 $\|Pz\|_2$;再对所有 $x-y$ 做 union bound。这就是 Johnson-Lindenstrauss 的全部骨架。

random projectionJLunion bound
这一层要会问 固定 $z$ 的失败概率是否足够小,能不能支付 $N^2$ 个点对的成本?
3

最后进入矩阵 MGF

矩阵 Bernstein 仍然是 Markov + MGF + 优化参数,只是 $e^{A+B}=e^Ae^B$ 失效。Lieb inequality 让 trace exponential 可以逐个剥离独立项。

Loewner orderLiebmatrix Bernstein
这一层要会问 我是在控制最大特征值、最小特征值,还是两者合起来的算子范数?
等周blow-upLipschitz 集中随机投影矩阵 Bernstein稀疏图 / 协方差

分层阅读路线

层次 先抓什么 推荐入口 暂时怎么处理
第一遍:主线阅读 等周集中、Lipschitz 函数、JL、矩阵 Bernstein 本章主线、关键定理卡片 先分清“几何集中”和“矩阵集中”两条路线。
第二遍:证明精读 blow-up 到函数集中、球面投影 tail、Lieb inequality、矩阵 MGF Theorem 5.1.3、Lemma 5.3.2、Theorem 5.4.1 完整证明 重点理解不可交换矩阵为什么需要 trace/Lieb 工具。
第三遍:习题与应用 几何集中、push-forward、JL、矩阵函数、稀疏图和协方差估计 Exercises 5.1-5.32 把每题归入几何、矩阵或应用三类。
专题回看 concentration without independence、矩阵集中、随机降维 第 10 部分概率集中与随机矩阵路线 为稀疏网络和一般协方差估计补工具。

初学者补充:两条集中路线

Reading Pattern 几何集中路线
  1. 找到一个度量测度空间,例如球面、Gaussian 空间或 Hamming cube。
  2. 知道该空间有等周不等式:给定测度时,某类简单集合的邻域最小。
  3. 推出 blow-up:测度至少 $1/2$ 的集合,其邻域很快接近全空间。
  4. 对 Lipschitz 函数取中位数 sublevel set,把 blow-up 转成尾界。
Reading Pattern 矩阵集中路线
  1. 把矩阵和的范数拆成最大特征值和最小特征值。
  2. 对最大特征值用 Markov:$\mathbb P\{\lambda_{\max}(S)\ge t\}\le e^{-\lambda t}\mathbb Ee^{\lambda\lambda_{\max}(S)}$。
  3. 用 trace 上界最大特征值的指数。
  4. 用 Lieb inequality 逐个处理独立矩阵项。
  5. 用单项 MGF 控制,得到 Bernstein 型尾界。

核心对象与符号表

符号 / 对象 含义 本章用途
$\|f\|_{\mathrm{Lip}}$ Lipschitz 范数 控制函数继承空间集中的尺度
$A_t$ 集合 $A$ 的 $t$-邻域 blow-up 和等周不等式的核心对象
$\sigma$ 球面归一化面积测度 把面积读成概率
$\gamma_n$ $\mathbb R^n$ 上 Gaussian measure Gaussian 等周与集中
$G_{n,m}$ $m$ 维子空间组成的 Grassmannian 随机投影与 JL
$P_E$ 到子空间 $E$ 的正交投影 随机投影、Grassmannian 距离
$\preceq$ Loewner order 比较对称矩阵、矩阵函数和 MGF
$\operatorname{tr}e^X$ trace exponential 矩阵 MGF 方法的可处理替代品
$\sigma^2=\|\sum\mathbb EX_i^2\|$ 矩阵 variance 参数 矩阵 Bernstein 的二次尺度
$r(\Sigma)=\operatorname{tr}(\Sigma)/\|\Sigma\|$ 有效秩 协方差估计中替代维度 $n$

关键定理卡片

Core TheoremTheorem 5.1.3:球面 Lipschitz 集中

条件:$X\sim\operatorname{Unif}(\sqrt n S^{n-1})$,$f$ 是 Lipschitz 函数。

结论:$\|f(X)-\mathbb Ef(X)\|_{\psi_2}\le C\|f\|_{\mathrm{Lip}}$。

用途:把高维球面几何转成函数值集中,是 JL 引理的关键输入。

跳到完整证明
Core LemmaLemma 5.1.6:Blow-up

条件:$A\subset\sqrt n S^{n-1}$ 且 $\sigma(A)\ge1/2$。

结论:$\sigma(A_t)\ge1-2e^{-ct^2}$。

用途:这是从等周不等式到函数集中的中间桥梁。

跳到完整证明
Core TheoremTheorem 5.3.1:Johnson-Lindenstrauss

条件:$N$ 个点,随机投影维度 $m\gtrsim\varepsilon^{-2}\log N$。

结论:缩放随机投影以高概率保持所有点对距离。

用途:解释高维数据可以随机降维到对数维。

跳到完整证明
Core TheoremTheorem 5.4.1:矩阵 Bernstein

条件:$X_i$ 独立、均值零、对称,且 $\|X_i\|\le K$。

结论:$\|\sum_iX_i\|$ 满足 Bernstein 型混合尾。

用途:控制随机图噪声矩阵、样本协方差偏差和矩阵 sketching。

跳到完整证明
Application TheoremTheorem 5.6.1:一般协方差估计

条件:$\|X\|_2\le K(\mathbb E\|X\|_2^2)^{1/2}$ a.s.。

结论:$m\asymp K^2n\log n$ 量级样本可给出相对谱范数误差。

用途:把协方差估计从次高斯分布推广到更一般的有界分布。

跳到完整证明

关键定理完整证明

Complete ProofLemma 5.1.6:Blow-up
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明球面上测度至少为 $1/2$ 的集合,其 $t$-邻域以 $1-e^{-ct^2}$ 的速度覆盖球面。

证明思路

用等周不等式把一般集合 $A$ 的邻域和同面积半球的邻域比较;再用球面随机向量的一维坐标次高斯尾界估计半球邻域面积。

完整证明

$$ H=\{x\in\sqrt nS^{n-1}:x_1\le0\}. $$

因为 $\sigma(H)=1/2$ 且 $\sigma(A)\ge1/2$,球面等周不等式给出

$$ \sigma(A_t)\ge\sigma(H_t). $$

若 $X\sim\operatorname{Unif}(\sqrt nS^{n-1})$,则 $\sigma(H_t)=\mathbb P\{X\in H_t\}$。由几何包含

$$ H_t\supset\{x\in\sqrt nS^{n-1}:x_1\le t/\sqrt2\}, $$

可得

$$ \sigma(H_t)\ge\mathbb P\{X_1\le t/\sqrt2\}. $$

Theorem 3.4.5 告诉我们 $\|X_1\|_{\psi_2}\le C$。由于 $X_1$ 关于 $0$ 对称,

$$ \mathbb P\{X_1>t/\sqrt2\}\le2e^{-ct^2}. $$

因此

$$ \sigma(A_t)\ge\sigma(H_t)\ge1-2e^{-ct^2}. $$
Complete ProofTheorem 5.1.3:球面 Lipschitz 集中
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:用 blow-up 证明球面上 Lipschitz 函数围绕均值次高斯集中。

证明思路

先归一化 $\|f\|_{\mathrm{Lip}}=1$。取中位数 sublevel set $A=\{f\le M\}$,其测度至少 $1/2$。若 $X\in A_t$,Lipschitz 条件给出 $f(X)\le M+t$。对 $f$ 和 $-f$ 分别处理两侧尾部,最后把中位数换成均值。

完整证明

若 $\|f\|_{\mathrm{Lip}}=0$,则 $f$ 在球面上常数,结论成立。否则令 $g=f/\|f\|_{\mathrm{Lip}}$,证明 $g$ 的结论后再缩放回来。因此只需考虑 $\|f\|_{\mathrm{Lip}}=1$。

令 $M$ 为 $f(X)$ 的中位数,并设

$$ A=\{x:f(x)\le M\}. $$

则 $\sigma(A)\ge1/2$。由 Lemma 5.1.6,

$$ \mathbb P\{X\in A_t\}\ge1-2e^{-ct^2}. $$

若 $X\in A_t$,存在 $y\in A$ 使 $\|X-y\|_2\le t$,所以

$$ f(X)\le f(y)+\|X-y\|_2\le M+t. $$

于是

$$ \mathbb P\{f(X)>M+t\}\le2e^{-ct^2}. $$

对 $-f$ 重复同一论证。因为 $-M$ 是 $-f(X)$ 的中位数,得到

$$ \mathbb P\{f(X)<M-t\}\le2e^{-ct^2}. $$

合并两侧,

$$ \mathbb P\{|f(X)-M|>t\}\le4e^{-ct^2}, $$

这等价于 $\|f(X)-M\|_{\psi_2}\le C$。进一步,

$$ |\mathbb Ef(X)-M|\le\mathbb E|f(X)-M|\le C\|f(X)-M\|_{\psi_2}\le C. $$

由三角不等式,

$$ \|f(X)-\mathbb Ef(X)\|_{\psi_2} \le \|f(X)-M\|_{\psi_2}+|M-\mathbb Ef(X)| \le C. $$

缩放回一般 Lipschitz 范数,得到结论。

Complete ProofLemma 5.3.2:随机投影
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:固定一个向量 $z$,证明随机投影 $Pz$ 的长度集中在 $\sqrt{m/n}\|z\|_2$ 附近。

证明思路

用旋转不变性把“随机子空间 + 固定向量”变成“固定子空间 + 随机球面向量”。期望由坐标对称性计算;集中由单位球面 Lipschitz 集中给出。

完整证明

齐次性允许假设 $\|z\|_2=1$。由 Haar 不变性,随机子空间 $E$ 可写作 $U\mathbb R^m$,其中 $U$ 是均匀随机旋转。因此 $\|P_Ez\|_2$ 与 $\|P_{\mathbb R^m}U^{-1}z\|_2$ 同分布,而 $U^{-1}z$ 均匀分布在 $S^{n-1}$ 上。

于是只需令 $z\sim\operatorname{Unif}(S^{n-1})$,$P$ 为前 $m$ 个坐标投影。由坐标同分布,

$$ \mathbb E\|Pz\|_2^2 = \sum_{i=1}^m\mathbb Ez_i^2 = m\mathbb Ez_1^2. $$

又 $\sum_i z_i^2=1$,取期望得 $n\mathbb Ez_1^2=1$。所以 $\mathbb E\|Pz\|_2^2=m/n$。

函数 $h(x)=\|Px\|_2$ 满足

$$ |h(x)-h(y)|\le\|P(x-y)\|_2\le\|x-y\|_2, $$

故 Lipschitz 范数至多 $1$。由单位球面集中,

$$ \mathbb P\left\{\left|\|Pz\|_2-\mathbb E\|Pz\|_2\right|>u\right\} \le2e^{-cnu^2}. $$

Remark 5.2.1 允许把中心从 $\mathbb E\|Pz\|_2$ 换成 $(\mathbb E\|Pz\|_2^2)^{1/2}=\sqrt{m/n}$,只改变常数。取 $u=\varepsilon\sqrt{m/n}$ 即得 (b)。一般 $\|z\|_2$ 由齐次性恢复。

Complete ProofTheorem 5.3.1:Johnson-Lindenstrauss
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明随机 $m$ 维投影在 $m\gtrsim\varepsilon^{-2}\log N$ 时同时保持 $N$ 个点的全部点对距离。

证明思路

对每个差向量 $z=x-y$ 使用 Lemma 5.3.2,再对至多 $N^2$ 个差向量做 union bound。

完整证明

$$ \mathcal D=\mathcal X-\mathcal X=\{x-y:x,y\in\mathcal X\}. $$

对固定 $z\in\mathcal D$,Lemma 5.3.2 给出

$$ (1-\varepsilon)\sqrt{m/n}\|z\|_2 \le \|Pz\|_2 \le (1+\varepsilon)\sqrt{m/n}\|z\|_2 $$

失败概率至多 $2e^{-c\varepsilon^2m}$。因为 $|\mathcal D|\le N^2$,union bound 给出所有 $z\in\mathcal D$ 同时成功的概率至少为

$$ 1-2N^2e^{-c\varepsilon^2m}. $$

若 $m\ge C\varepsilon^{-2}\log N$,调整常数得

$$ 2N^2e^{-c\varepsilon^2m} \le 2e^{-c'\varepsilon^2m}. $$

最后乘以缩放因子 $Q=\sqrt{n/m}P$,即得到对所有 $x,y\in\mathcal X$ 的距离保持结论。

Complete ProofProposition 5.4.4:Loewner order 性质
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:验证 Loewner order 与特征值、trace、算子范数和矩阵函数的基本关系。

证明思路

始终把 $X\preceq Y$ 翻译为 $u^{\mathsf T}Xu\le u^{\mathsf T}Yu$,再用 min-max、谱分解和逐特征值函数演算。

完整证明

(a) 由 $X\preceq Y$,对所有单位向量 $u$ 有 $u^{\mathsf T}Xu\le u^{\mathsf T}Yu$。把这一点代入 Courant-Fischer min-max 公式,得到每个有序特征值均满足 $\lambda_i(X)\le\lambda_i(Y)$。

(b) $f(X)$ 的特征值是 $f(\lambda_i(X))$。由 (a) 和 $f$ 弱增,逐项求和得到 $\operatorname{tr}f(X)\le\operatorname{tr}f(Y)$。

(c) 若 $\|X\|\le a$,则所有特征值位于 $[-a,a]$,所以 $-aI\preceq X\preceq aI$。反过来若该 Loewner 夹逼成立,则所有 Rayleigh quotient 位于 $[-a,a]$,故所有特征值绝对值不超过 $a$,即 $\|X\|\le a$。

(d) 对 $h=g-f$,问题化为证明 $h(X)\succeq0$。若 $\|X\|\le a$,则 $X$ 的所有特征值 $\lambda_i$ 满足 $|\lambda_i|\le a$,于是 $h(\lambda_i)\ge0$。因此 $h(X)$ 的所有特征值非负,故 $h(X)\succeq0$。

Complete ProofLemma 5.4.10:矩阵 MGF 控制
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明有界零均值随机矩阵的矩阵 MGF 被二阶矩矩阵控制。

证明思路

先证明标量 Taylor 上界,再用 Proposition 5.4.4(d) 升级为矩阵不等式,最后取期望并用 $I+Z\preceq e^Z$。

完整证明

标量 Taylor 展开和 $p!\ge2\cdot3^{p-2}$ 给出:当 $|z|<3$ 时,

$$ e^z\le1+z+\frac{z^2/2}{1-|z|/3}. $$

令 $z=\lambda x$。若 $|x|\le K$ 且 $|\lambda|<3/K$,则

$$ e^{\lambda x}\le1+\lambda x+g(\lambda)x^2, \qquad g(\lambda)=\frac{\lambda^2/2}{1-|\lambda|K/3}. $$

对满足 $\|X\|\le K$ 的对称矩阵 $X$,由 Proposition 5.4.4(d),

$$ e^{\lambda X}\preceq I+\lambda X+g(\lambda)X^2. $$

取期望并用 $\mathbb EX=0$,得到

$$ \mathbb Ee^{\lambda X}\preceq I+g(\lambda)\mathbb EX^2. $$

由于标量不等式 $1+u\le e^u$ 可升级为矩阵不等式,$I+Z\preceq e^Z$ 对所有对称 $Z$ 成立。取 $Z=g(\lambda)\mathbb EX^2$,结论成立。

Complete ProofTheorem 5.4.1:矩阵 Bernstein
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明独立有界均值零对称随机矩阵和的算子范数满足 Bernstein 型尾界。

证明思路

对最大特征值做指数矩方法;用 trace exponential 控制最大特征值;用 Lieb inequality 递归剥离独立项;用 Lemma 5.4.10 控制单项矩阵 MGF;最后对 $S$ 和 $-S$ 合并。

完整证明

令 $S=\sum_iX_i$。由

$$ \|S\|=\max(\lambda_{\max}(S),\lambda_{\max}(-S)), $$

先控制 $\lambda_{\max}(S)$。对任意 $\lambda>0$,Markov 不等式给出

$$ \mathbb P\{\lambda_{\max}(S)\ge t\} \le e^{-\lambda t}\mathbb Ee^{\lambda\lambda_{\max}(S)}. $$

因为 $e^{\lambda S}$ 的最大特征值不超过其 trace,

$$ \mathbb Ee^{\lambda\lambda_{\max}(S)} \le \mathbb E\operatorname{tr}e^{\lambda S}. $$

对最后一项条件化并应用 Lemma 5.4.9,然后递归处理 $X_N,\ldots,X_1$,得到

$$ \mathbb E\operatorname{tr}e^{\lambda S} \le \operatorname{tr}\exp\left(\sum_i\log\mathbb Ee^{\lambda X_i}\right). $$

Lemma 5.4.10 给出

$$ \mathbb Ee^{\lambda X_i}\preceq\exp(g(\lambda)\mathbb EX_i^2). $$

由 $\log$ 的矩阵单调性和 trace monotonicity,

$$ \mathbb E\operatorname{tr}e^{\lambda S} \le \operatorname{tr}\exp(g(\lambda)Z), \qquad Z=\sum_i\mathbb EX_i^2. $$

由于 $Z\succeq0$,

$$ \operatorname{tr}\exp(g(\lambda)Z) \le n\exp(g(\lambda)\|Z\|) = n\exp(g(\lambda)\sigma^2). $$

所以

$$ \mathbb P\{\lambda_{\max}(S)\ge t\} \le n\exp[-\lambda t+g(\lambda)\sigma^2]. $$

取 $\lambda=t/(\sigma^2+Kt/3)$,它满足 $\lambda<3/K$,并化简得

$$ \mathbb P\{\lambda_{\max}(S)\ge t\} \le n\exp\left(-\frac{t^2/2}{\sigma^2+Kt/3}\right). $$

对 $-S$ 应用同一论证。由于 $-X_i$ 仍独立、均值零、范数不超过 $K$,且 $\mathbb E(-X_i)^2=\mathbb EX_i^2$,得到相同界。把两侧事件并合并,得到

$$ \mathbb P\{\|S\|\ge t\} \le 2n\exp\left(-\frac{t^2/2}{\sigma^2+Kt/3}\right). $$
Complete ProofTheorem 5.5.1:稀疏 SBM 谱聚类
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明在 $(a-b)^2\gtrsim a\log n$ 时,谱聚类可在稀疏 SBM 中达到 $99\%$ 准确率。

证明思路

把邻接矩阵写成 $A=D+R$。第 4 章已经知道 $D$ 的第二特征向量编码社群;本章用矩阵 Bernstein 证明 $\|R\|\lesssim\sqrt{a\log n}$,再用 Davis-Kahan 把谱范数误差转为特征向量误差。

完整证明

设 $D=\mathbb EA$,$R=A-D$。把 $R$ 分解为独立对称矩阵和 $R=\sum_{i\le j}Z_{ij}$,其中 $Z_{ij}$ 只保留 entry $(i,j)$ 和 $(j,i)$ 的中心化噪声。由于 $A_{ij}\in\{0,1\}$,$\|Z_{ij}\|\le1$。

用矩阵 Bernstein 的期望形式与 Markov 不等式,以概率至少 $0.99$ 得

$$ \|R\|\lesssim\sigma\sqrt{\log n}+\log n, \qquad \sigma^2=\left\|\mathbb E\sum_{i\le j}Z_{ij}^2\right\|. $$

直接计算 $Z_{ij}^2$ 是对角矩阵,且

$$ \sigma^2 = \max_i\sum_{j=1}^n\mathbb ER_{ij}^2 \le np = a. $$

由 $(a-b)^2\ge Ca\log n$ 和 $a-b<a$ 得 $a\ge C\log n$,因此

$$ \|R\|\lesssim\sqrt{a\log n}. $$

第 4.5 节给出 $D$ 的目标谱间隙

$$ \delta = \min(\lambda_2(D),\lambda_1(D)-\lambda_2(D)) = \frac{a-b}{2}, $$

这里使用 $a<3b$。Davis-Kahan 给出存在 $\theta\in\{-1,1\}$,使

$$ \|\bar u_2(D)-\theta\bar u_2(A)\|_2 \le \frac{2\|R\|}{\delta} \le \frac{C_1\sqrt{a\log n}}{a-b}. $$

若定理中的常数 $C$ 足够大,该误差小于 $0.1$。把单位向量乘以 $\sqrt n$,得到符号向量误差至多 $0.1\sqrt n$。每个符号错误坐标至少贡献 $1$ 的平方误差,因此错误坐标数至多 $0.01n$。这给出 $99\%$ 准确率。

Complete ProofTheorem 5.6.1:一般协方差估计
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:在只假设样本范数有界的情况下,用矩阵 Bernstein 控制样本协方差的谱范数误差。

证明思路

把 $\Sigma_m-\Sigma$ 写成均值零矩阵和。矩阵 Bernstein 的期望形式需要两个量:单项范数上界 $M$ 与 variance 参数 $\sigma^2$。二者都由 $\|X\|_2^2\le K^2\operatorname{tr}\Sigma$ 控制。

完整证明

由 Proposition 3.2.1(b),$\mathbb E\|X\|_2^2=\operatorname{tr}\Sigma$,所以假设给出

$$ \|X\|_2^2\le K^2\operatorname{tr}\Sigma. $$

令 $Y_i=X_iX_i^{\mathsf T}-\Sigma$。则 $\mathbb EY_i=0$,且

$$ \Sigma_m-\Sigma=\frac1m\sum_{i=1}^mY_i. $$

矩阵 Bernstein 期望形式给出

$$ \mathbb E\|\Sigma_m-\Sigma\| \lesssim \frac1m(\sigma\sqrt{\log n}+M\log n), $$

其中

$$ \sigma^2=\left\|\sum_{i=1}^m\mathbb EY_i^2\right\|, \qquad \|Y_i\|\le M. $$

先估计 variance。展开得

$$ \mathbb E(XX^{\mathsf T}-\Sigma)^2 = \mathbb E(XX^{\mathsf T})^2-\Sigma^2 \preceq \mathbb E(XX^{\mathsf T})^2. $$

$$ (XX^{\mathsf T})^2=\|X\|_2^2XX^{\mathsf T} \preceq K^2\operatorname{tr}(\Sigma)XX^{\mathsf T}. $$

取期望,

$$ \mathbb E(XX^{\mathsf T})^2 \preceq K^2\operatorname{tr}(\Sigma)\Sigma. $$

因此

$$ \sigma^2 \le K^2m\operatorname{tr}(\Sigma)\|\Sigma\|. $$

再估计单项范数:

$$ \|Y_i\| \le \|X_iX_i^{\mathsf T}\|+\|\Sigma\| = \|X_i\|_2^2+\|\Sigma\| \le 2K^2\operatorname{tr}(\Sigma). $$

代回并用 $\operatorname{tr}(\Sigma)\le n\|\Sigma\|$,得到

$$ \mathbb E\|\Sigma_m-\Sigma\| \le C\left( \sqrt{\frac{K^2n\log n}{m}} + \frac{K^2n\log n}{m} \right)\|\Sigma\|. $$
Complete ProofTheorem 5.7.1:Bounded differences inequality
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明独立输入函数在单坐标影响有界时满足 McDiarmid 上尾不等式。

证明思路

构造 Doob martingale $M_i=\mathbb E[f(X)\mid X_1,\ldots,X_i]$。单坐标有界差分保证 martingale increments 有界,然后应用 Azuma-Hoeffding。

完整证明

$$ M_i=\mathbb E[f(X)\mid X_1,\ldots,X_i], \qquad i=0,\ldots,N. $$

则 $M_0=\mathbb Ef(X)$,$M_N=f(X)$,且 $(M_i)$ 是 martingale。改变第 $i$ 个坐标最多改变 $f$ 的值 $c_i$,由条件期望的凸组合性质,martingale 差分满足

$$ |M_i-M_{i-1}|\le c_i. $$

Azuma-Hoeffding 不等式给出

$$ \mathbb P\{M_N-M_0\ge t\} \le \exp\left(-\frac{2t^2}{\sum_{i=1}^N c_i^2}\right). $$

这就是所需结论。

正文隐藏验证补全

Hidden Check5.1:为什么任意函数不可能都集中
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:说明若不限制函数振荡,$f(X)$ 可以完全不集中。

完整证明

取任意事件 $A$ 满足 $\mathbb P\{X\in A\}=1/2$,令 $f=\mathbf 1_A$。则 $f(X)$ 是 Bernoulli$(1/2)$ 随机变量,均值为 $1/2$,并且 $|f(X)-\mathbb Ef(X)|=1/2$ 以概率 $1$ 成立。若再把函数放大为 $Mf$,偏差规模变成 $M/2$。因此没有 Lipschitz 或有界振荡等限制时,不可能存在统一集中界。

Hidden Check5.1.6:为什么 (5.2) 成立
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明若 $x_1\le t/\sqrt2$,则 $x$ 到半球 $H=\{y:y_1\le0\}$ 的距离至多 $t$。

完整证明

若 $x_1\le0$,则 $x\in H$,距离为 $0$。下面设 $0<x_1=s\le t/\sqrt2$,球面半径为 $R=\sqrt n$。令 $x'=(x_2,\ldots,x_n)$,取

$$y=(0,\alpha x'),\qquad \alpha=\frac{R}{\|x'\|_2}=\frac{R}{\sqrt{R^2-s^2}}.$$

则 $y\in\sqrt nS^{n-1}$ 且 $y_1=0$,所以 $y\in H$。并且

$$\|x-y\|_2^2=s^2+\left(R-\sqrt{R^2-s^2}\right)^2\le2s^2\le t^2.$$

因此 $x\in H_t$。

Hidden Check5.1.4:为什么可假设 $\|f\|_{\mathrm{Lip}}=1$
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:解释 Lipschitz 范数归一化不会损失一般性。

完整证明

若 $\|f\|_{\mathrm{Lip}}=0$,则 $f$ 在连通球面上为常数,结论成立。若 $\|f\|_{\mathrm{Lip}}>0$,定义 $g=f/\|f\|_{\mathrm{Lip}}$,则 $\|g\|_{\mathrm{Lip}}=1$。若已知 $\|g(X)-\mathbb Eg(X)\|_{\psi_2}\le C$,两边乘以 $\|f\|_{\mathrm{Lip}}$ 即得 $\|f(X)-\mathbb Ef(X)\|_{\psi_2}\le C\|f\|_{\mathrm{Lip}}$。

Hidden Check5.1.4:补全对 $-f$ 的论证
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:从已证明的右尾界推出左尾界。

完整证明

若 $M$ 是 $f(X)$ 的中位数,则 $-M$ 是 $-f(X)$ 的中位数。对函数 $-f$ 应用同样论证,得到

$$\mathbb P\{-f(X)>-M+t\}\le2e^{-ct^2}.$$

该事件等价于 $f(X)<M-t$,因此

$$\mathbb P\{f(X)<M-t\}\le2e^{-ct^2}.$$
Hidden Check5.1.4:从中位数集中换到均值集中
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明若 $Z$ 围绕中位数 $M$ 次高斯集中,则 $Z$ 也围绕均值次高斯集中。

完整证明

设 $\|Z-M\|_{\psi_2}\le K$。由 Orlicz 范数控制一阶矩,$\mathbb E|Z-M|\le CK$。因此

$$|\mathbb EZ-M|\le\mathbb E|Z-M|\le CK.$$

再用三角不等式,

$$\|Z-\mathbb EZ\|_{\psi_2}\le\|Z-M\|_{\psi_2}+|M-\mathbb EZ|\le C'K.$$
Hidden Check5.2.8:Haar 构造的旋转不变性
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:说明 Gaussian 矩阵的 SVD 构造给出旋转不变的正交矩阵分布。

完整证明

若 $G$ 的 entries 独立 $N(0,1)$,则对任意固定正交矩阵 $W$,$WG$ 与 $G$ 同分布。若 $G=U\Sigma V^{\mathsf T}$,则 $WG=(WU)\Sigma V^{\mathsf T}$。因此由 $G$ 生成的正交因子 $UV^{\mathsf T}$ 与 $WUV^{\mathsf T}$ 同分布。这就是左旋转不变性;右旋转不变性同理来自 $GW$ 与 $G$ 同分布。故该分布是 Haar 型均匀分布。

Hidden Check5.2.6:Gaussian 矩阵像空间的旋转不变性
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明 $n\times m$ Gaussian 矩阵的列空间给出 Grassmannian 上的均匀随机子空间。

完整证明

设 $G$ 是 $n\times m$ 标准 Gaussian 矩阵。对任意固定正交矩阵 $U$,$UG$ 与 $G$ 同分布。于是 $\operatorname{im}(UG)=U(\operatorname{im}G)$ 与 $\operatorname{im}G$ 同分布。这说明列空间分布在任意旋转下不变,即为 $G_{n,m}$ 上的 Haar 概率测度。

Hidden Check5.3.2:随机子空间和随机向量视角等价
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:解释为什么可固定子空间、随机旋转向量。

完整证明

随机子空间可写为 $E=U E_0$,其中 $E_0=\mathbb R^m$,$U$ 是 Haar 正交矩阵。令 $P_E=UP_{E_0}U^{\mathsf T}$。则

$$\|P_Ez\|_2=\|UP_{E_0}U^{\mathsf T}z\|_2=\|P_{E_0}(U^{\mathsf T}z)\|_2.$$

由于 $U^{\mathsf T}z$ 均匀分布在半径 $\|z\|_2$ 的球面上,结论成立。

Hidden Check5.3.2:$x\mapsto\|Px\|_2$ 是 1-Lipschitz
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:验证投影长度函数的 Lipschitz 范数至多为 $1$。

完整证明

对任意 $x,y$,由反三角不等式和投影算子范数 $\|P\|\le1$,

$$|\|Px\|_2-\|Py\|_2|\le\|P(x-y)\|_2\le\|x-y\|_2.$$

因此该函数是 $1$-Lipschitz。

Hidden Check5.4.2:为什么 Jensen 可用于随机矩阵
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:说明 Lieb inequality 中的凹函数可对矩阵值随机变量使用 Jensen。

完整证明

正定矩阵空间是有限维实向量空间中凸集。若 $f$ 在该凸集上凹,则对有限支撑随机矩阵,Jensen 不等式由凹性和归纳直接得到。一般随机矩阵可用简单函数逼近并用连续性/单调收敛处理。因此

$$\mathbb Ef(X)\le f(\mathbb EX)$$

对这里的矩阵值随机变量成立。

Hidden Check5.4.3:补全 $-S$ 与算子范数合并
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:把最大特征值尾界升级为算子范数尾界。

完整证明

已证

$$\mathbb P\{\lambda_{\max}(S)\ge t\}\le n\exp\left(-\frac{t^2/2}{\sigma^2+Kt/3}\right).$$

对 $-S=\sum_i(-X_i)$ 应用同一结果,因为 $-X_i$ 仍独立均值零,$\|-X_i\|=\|X_i\|\le K$,且 $\mathbb E(-X_i)^2=\mathbb EX_i^2$。于是

$$\mathbb P\{\lambda_{\max}(-S)\ge t\}\le n\exp\left(-\frac{t^2/2}{\sigma^2+Kt/3}\right).$$

又 $\|S\|=\max(\lambda_{\max}(S),\lambda_{\max}(-S))$,union bound 给出最终的因子 $2n$。

Hidden Check5.5.1:为什么 $\|Z_{ij}\|\le1$
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:计算只含一对对称 entry 的矩阵算子范数。

完整证明

当 $i<j$ 时,$Z_{ij}=R_{ij}(e_ie_j^{\mathsf T}+e_je_i^{\mathsf T})$。它只在 $\operatorname{span}\{e_i,e_j\}$ 上非零,对应 $2\times2$ 块

$$\begin{bmatrix}0&R_{ij}\\R_{ij}&0\end{bmatrix},$$

特征值为 $\pm R_{ij}$,所以 $\|Z_{ij}\|=|R_{ij}|\le1$。当 $i=j$ 时,$Z_{ii}=R_{ii}e_ie_i^{\mathsf T}$,范数同样为 $|R_{ii}|\le1$。

Hidden Check5.5.1:为什么假设推出 $a\gtrsim\log n$
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:补全由 $(a-b)^2\ge Ca\log n$ 到 $a\gtrsim\log n$ 的代数步骤。

完整证明

因为 $b<a$,所以 $0<a-b<a$。由假设

$$Ca\log n\le(a-b)^2<a^2.$$

若 $a>0$,两边除以 $a$ 得 $a>C\log n$。这正是 $a\gtrsim\log n$。

Hidden Check5.5.1:特征向量误差推出 99% 符号正确
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:把 $\ell_2$ 特征向量误差转成误分类比例上界。

完整证明

令 $u=u_2(D)\in\{\pm1\}^n$ 是社群真值符号向量,$\hat u=\theta u_2(A)$。若 $\|u-\hat u\|_2\le0.1\sqrt n$,则设错误坐标集合为 $E=\{j:\operatorname{sign}\hat u_j\ne u_j\}$。对每个 $j\in E$,因为 $u_j=\pm1$ 且 $\hat u_j$ 符号相反,有 $|u_j-\hat u_j|\ge1$。所以

$$|E|\le\sum_{j\in E}|u_j-\hat u_j|^2\le\|u-\hat u\|_2^2\le0.01n.$$

故误分类顶点不超过 $1\%$。

Hidden Check5.6.3:为什么 $r(\Sigma)\le n$
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明正半定矩阵有效秩不超过维度。

完整证明

设 $\Sigma$ 的特征值为 $\lambda_1,\ldots,\lambda_n\ge0$,最大特征值为 $\lambda_{\max}=\|\Sigma\|$。若 $\Sigma\ne0$,则

$$r(\Sigma)=\frac{\sum_i\lambda_i}{\lambda_{\max}}\le\frac{n\lambda_{\max}}{\lambda_{\max}}=n.$$

若 $\Sigma=0$,有效秩通常不使用;也可约定为 $0$,仍不超过 $n$。

公式卡片

公式 读法
$\|f(X)-\mathbb Ef(X)\|_{\psi_2}\le C\|f\|_{\mathrm{Lip}}$ 空间集中传给 Lipschitz 观测
$\sigma(A_t)\ge1-2e^{-ct^2}$ 半空间/半球 blow-up 的概率形式
$m\gtrsim\varepsilon^{-2}\log N$ JL 的维度成本来自 $N^2$ 个点对的 union bound
$\mathbb P\{\|\sum_iX_i\|\ge t\}\le2n\exp[-(t^2/2)/(\sigma^2+Kt/3)]$ 矩阵 Bernstein,高维代价是 $n$
$\mathbb E\|\Sigma_m-\Sigma\|\lesssim(\sqrt{K^2n\log n/m}+K^2n\log n/m)\|\Sigma\|$ 一般协方差估计的样本复杂度
$r(\Sigma)=\operatorname{tr}\Sigma/\|\Sigma\|$ 用有效维度替代粗糙维度

Exercises 5.1-5.32 证明工作区

批次 建议目标
5.1-5.7 Lipschitz、blow-up、均值/中位数和几何集中基础
5.8-5.15 Gaussian concentration、push-forward、JL 与最优性
5.16-5.19 矩阵函数、矩阵单调性和矩阵指数
5.20-5.24 矩阵 Bernstein/Hoeffding/Khintchine 变体
5.25-5.32 稀疏 SBM、协方差估计、有效秩、frame、sketching
Exercise 5.1连续、可微与 Lipschitz
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:验证 Lipschitz、一致连续、可微有界梯度之间的基本关系,并给出反例。

证明思路

用 $\varepsilon$-$\delta$ 证明一致连续,用线段积分证明梯度上界推出 Lipschitz。

完整证明

(a) 若 $d_Y(f(u),f(v))\le Ld_X(u,v)$,给定 $\varepsilon>0$,取 $\delta=\varepsilon/L$ 即得一致连续;$L=0$ 时函数常值。

(b) 对 $x,y\in\mathbb R^n$,由微积分基本定理,

$$f(x)-f(y)=\int_0^1\langle\nabla f(y+t(x-y)),x-y\rangle\,dt.$$

取绝对值并用 Cauchy-Schwarz,得

$$|f(x)-f(y)|\le\sup_z\|\nabla f(z)\|_2\|x-y\|_2.$$

(c) $f(x)=\sqrt{|x|}$ 在 $[-1,1]$ 上一致连续但不是 Lipschitz,因为 $|f(h)-f(0)|/|h|=1/\sqrt{|h|}\to\infty$。(d) $f(x)=|x|$ Lipschitz,常数为 $1$,但在 $0$ 不可微。

Exercise 5.2Example 5.1.2 验证
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:计算线性泛函、矩阵算子和范数函数的 Lipschitz 范数。

证明思路

上界用 Cauchy-Schwarz 或算子范数定义;下界取达到或逼近上确界的方向。

完整证明

(a) $|\langle x-y,\theta\rangle|\le\|\theta\|_2\|x-y\|_2$,取 $x-y$ 平行于 $\theta$ 得等号,所以 Lipschitz 范数为 $\|\theta\|_2$。(b) $\|Ax-Ay\|_2\le\|A\|\|x-y\|_2$,而算子范数定义保证该常数最小。(c) 由三角不等式,$|\|x\|-\|y\||\le\|x-y\|$。若 $\|z\|\le L\|z\|_2$ 对所有 $z$ 成立,则 $f$ 是 $L$-Lipschitz;反过来取 $y=0$,任何 Lipschitz 常数 $L$ 都必须满足 $\|x\|\le L\|x\|_2$。

Exercise 5.3指数小集合的 blow-up
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:把 Lemma 5.1.6 的 $1/2$ 初始测度扩展到指数小测度。

证明思路

先反证 $A_s$ 的测度超过 $1/2$;再对 $A_s$ 应用 Lemma 5.1.6。

完整证明

(a) 反设 $\sigma(A_s)\le1/2$,令 $B=(A_s)^c$,则 $\sigma(B)\ge1/2$。由 Lemma 5.1.6,$\sigma(B_s)\ge1-2e^{-cs^2}$。但 $B_s\cap A=\varnothing$,否则 $A$ 中某点到 $B$ 距离至多 $s$,与 $B\subset(A_s)^c$ 矛盾。因此 $\sigma(A)\le2e^{-cs^2}$,矛盾。

(b) 由 (a),$\sigma(A_s)>1/2$。对 $A_s$ 应用 Lemma 5.1.6,得 $\sigma((A_s)_t)\ge1-2e^{-ct^2}$。且 $(A_s)_t\subset A_{s+t}\subset A_{2t}$,因为 $t\ge s$。结论成立。

Exercise 5.4Geodesic metric
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:说明球面集中也适用于 geodesic Lipschitz 函数。

证明思路

比较弧长距离和 chordal Euclidean 距离。

完整证明

单位球面上两点夹角为 $\theta\in[0,\pi]$ 时,geodesic 距离为 $\theta$,Euclidean 距离为 $2\sin(\theta/2)$。有 $\theta\le(\pi/2)2\sin(\theta/2)$。因此若 $f$ 关于 geodesic 距离 Lipschitz 常数为 $L$,则关于 Euclidean 距离 Lipschitz 常数至多 $(\pi/2)L$。把 Theorem 5.1.3 用于 Euclidean 距离即可,常数吸收 $\pi/2$。

Exercise 5.5单位球面集中
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:把半径 $\sqrt n$ 的球面集中缩放到单位球面。

证明思路

令 $Y=\sqrt n X$,把单位球面函数转成半径 $\sqrt n$ 球面函数。

完整证明

设 $X\sim\operatorname{Unif}(S^{n-1})$,令 $Y=\sqrt n X$,则 $Y\sim\operatorname{Unif}(\sqrt nS^{n-1})$。定义 $g(y)=f(y/\sqrt n)$。则

$$\|g\|_{\mathrm{Lip}}\le\frac1{\sqrt n}\|f\|_{\mathrm{Lip}}.$$

由 Theorem 5.1.3,

$$\|f(X)-\mathbb Ef(X)\|_{\psi_2}=\|g(Y)-\mathbb Eg(Y)\|_{\psi_2}\le\frac{C\|f\|_{\mathrm{Lip}}}{\sqrt n}.$$

尾界形式由 $\psi_2$ 范数定义推出。

Exercise 5.6均值和中位数等价
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明随机变量围绕均值集中和围绕中位数集中只差绝对常数。

证明思路

用 $\psi_2$ 控制一阶矩,证明均值和中位数距离被集中尺度控制。

完整证明

设 $K=\|Z-M\|_{\psi_2}$。则 $\mathbb E|Z-M|\le CK$,所以 $|\mathbb EZ-M|\le CK$。于是

$$\|Z-\mathbb EZ\|_{\psi_2}\le\|Z-M\|_{\psi_2}+|M-\mathbb EZ|\le C'K.$$

反过来设 $L=\|Z-\mathbb EZ\|_{\psi_2}$。若 $a=|M-\mathbb EZ|$,中位数性质给出 $\mathbb P\{|Z-\mathbb EZ|\ge a\}\ge1/2$,从次高斯尾界可得 $a\le CL$。因此

$$\|Z-M\|_{\psi_2}\le\|Z-\mathbb EZ\|_{\psi_2}+|\mathbb EZ-M|\le C'L.$$
Exercise 5.7集中等价于 blow-up
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:由所有 Lipschitz 函数集中推出集合 blow-up。

证明思路

取距离函数 $f(x)=d(x,A)$。它是 1-Lipschitz,且中位数为 $0$。

完整证明

令 $f(x)=d(x,A)$。则 $f$ 是 1-Lipschitz,并且因为 $\sigma(A)\ge1/2$,$0$ 是 $f(X)$ 的一个中位数。由假设和 Exercise 5.6,$\|f(X)\|_{\psi_2}\le CK$。于是

$$\mathbb P\{X\notin A_t\}=\mathbb P\{f(X)>t\}\le2\exp(-ct^2/K^2).$$

因此 $\sigma(A_t)\ge1-2\exp(-ct^2/K^2)$。

Exercise 5.8Gaussian concentration
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:由 Gaussian 等周不等式推出 Lipschitz 函数的 Gaussian concentration。

证明思路

完全仿照球面:用 half-space 作为等周最小化对象,先得 blow-up,再对中位数 sublevel set 应用。

完整证明

Gaussian 等周不等式说明,在给定 Gaussian measure 下,half-space 的邻域测度最小。若 $\gamma_n(A)\ge1/2$,与某个 measure 为 $1/2$ 的 half-space 比较,并用一维 Gaussian 尾界,得

$$\gamma_n(A_t)\ge1-2e^{-ct^2}.$$

对 Lipschitz 函数 $f$ 取中位数 $M$ 和集合 $A=\{f\le M\}$。若 $\|f\|_{\mathrm{Lip}}=1$,则 $X\in A_t$ 推出 $f(X)\le M+t$。对 $f$ 和 $-f$ 分别处理,再用 Exercise 5.6 把中位数换成均值,得到 Theorem 5.2.3。

Exercise 5.9Gaussian 最大值集中
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明 Gaussian 向量坐标最大值围绕均值次高斯集中。

证明思路

最大值函数是 Lipschitz 的;联合正态向量可写成标准 Gaussian 的线性像。

完整证明

(a) 函数 $g(x)=\max_i x_i$ 满足

$$|g(x)-g(y)|\le\max_i|x_i-y_i|\le\|x-y\|_2,$$

故 $1$-Lipschitz。对 $X\sim N(0,I_n)$ 用 Gaussian concentration 得结论。

(b) 联合正态向量可写为 $X=BG$,其中 $G\sim N(0,I_r)$。令 $b_i$ 是 $B$ 的第 $i$ 行,则 $\operatorname{Var}(X_i)=\|b_i\|_2^2$。函数 $h(g)=\max_i\langle b_i,g\rangle$ 的 Lipschitz 常数至多 $\max_i\|b_i\|_2$。用 Gaussian concentration 得结论。

Exercise 5.10$L^p$ 中心
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明中心可从均值换成 $L^p$ 范数,代价为 $\sqrt p$。

证明思路

用 subgaussian moment bound 控制 $\|Z-\mathbb EZ\|_{L^p}$,再用反三角不等式控制两个中心的距离。

完整证明

令 $\mu=\mathbb EZ$,$K=\|Z-\mu\|_{\psi_2}$。次高斯矩估计给出 $\|Z-\mu\|_{L^p}\le C\sqrt pK$。由反三角不等式,

$$\bigl|\|Z\|_{L^p}-|\mu|\bigr|\le\|Z-\mu\|_{L^p}\le C\sqrt pK.$$

由于 $\mu\ge0$,$|\mu|=\mu$。因此

$$\|Z-\|Z\|_{L^p}\|_{\psi_2}\le\|Z-\mu\|_{\psi_2}+|\mu-\|Z\|_{L^p}|\le C\sqrt pK.$$
Exercise 5.11Probability integral transform
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明标准正态 CDF 把标准正态变量推送为 $[0,1]$ 上均匀变量。

证明思路

一维用分布函数变换;多维用独立性。

完整证明

若 $U=\Phi(Z)$ 且 $Z\sim N(0,1)$,则对 $u\in[0,1]$,

$$\mathbb P\{U\le u\}=\mathbb P\{Z\le\Phi^{-1}(u)\}=u.$$

所以 $U\sim\operatorname{Unif}[0,1]$。因为 $Z_1,\ldots,Z_n$ 独立,坐标变换后 $\Phi(Z_i)$ 仍独立,故向量均匀分布在 $[0,1]^n$。

Exercise 5.12连续 cube 集中
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:用 Gaussian concentration 证明 $[0,1]^n$ 上的均匀集中。

证明思路

把 Gaussian 坐标用 CDF 推到 uniform 坐标,并检查推送映射 Lipschitz。

完整证明

令 $\phi(z)=(\Phi(z_1),\ldots,\Phi(z_n))$。由 Exercise 5.11,$\phi(Z)\sim\operatorname{Unif}([0,1]^n)$。又 $\Phi'(t)\le1/\sqrt{2\pi}$,所以

$$\|\phi(z)-\phi(w)\|_2\le(2\pi)^{-1/2}\|z-w\|_2.$$

若 $f$ 在 cube 上 Lipschitz,则 $f\circ\phi$ 在 $\mathbb R^n$ 上 Lipschitz,常数至多 $C\|f\|_{\mathrm{Lip}}$。对 $f(\phi(Z))$ 应用 Gaussian concentration 即得结论。

Exercise 5.13连续 ball 集中
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:说明可用 Lipschitz push-forward 把 Gaussian concentration 转成 ball 上均匀集中。

证明思路

使用径向传输:Gaussian 的方向本来均匀,只需把半径分布变成 ball 的均匀半径分布。

完整证明

把 $Z\sim N(0,I_n)$ 写成 $Z=RU$,其中 $U\sim\operatorname{Unif}(S^{n-1})$,$R=\|Z\|_2$,且二者独立。令 $F$ 为 $R$ 的分布函数。定义

$$\phi(0)=0,\qquad \phi(z)=\sqrt n\,F(\|z\|_2)^{1/n}\frac{z}{\|z\|_2}\quad(z\ne0).$$

因为 ball $\sqrt nB_2^n$ 的均匀半径分布满足 $\mathbb P\{\rho\le r\}=(r/\sqrt n)^n$,所以 $\phi(Z)$ 均匀分布在 $\sqrt nB_2^n$ 上。径向函数 $r\mapsto\sqrt nF(r)^{1/n}$ 是绝对常数 Lipschitz,并且不超过常数倍 $r$;结合径向映射的标准 Lipschitz 判别,$\phi$ 的 Lipschitz 范数被绝对常数控制。于是对 $f\circ\phi$ 应用 Gaussian concentration,得到 ball 上的集中。

Exercise 5.14次高斯 JL
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:用独立各向同性次高斯 rows 构造 Johnson-Lindenstrauss 映射。

证明思路

固定差向量后,$Az$ 的坐标是独立次高斯,Theorem 3.1.1 控制长度;再 union bound。

完整证明

固定 $z$。由于 rows $A_i$ 各向同性,$\mathbb E\langle A_i,z\rangle^2=\|z\|_2^2$;次高斯性保证 $\langle A_i,z\rangle$ 是次高斯。由 norm concentration,

$$\mathbb P\left\{\left|\frac1{\sqrt m}\|Az\|_2-\|z\|_2\right|>\varepsilon\|z\|_2\right\}\le2e^{-c\varepsilon^2m}.$$

对 $\mathcal X-\mathcal X$ 中至多 $N^2$ 个差向量做 union bound,若 $m\gtrsim\varepsilon^{-2}\log N$,得到 JL 结论。若 $A$ entries 为独立 Rademacher $\pm1$,rows 均值零、各向同性且次高斯,所以是特例。

Exercise 5.15JL 最优性
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明保持 $N$ 点距离的目标维度至少是 $\log N$ 量级。

证明思路

用 packing 体积比较限制低维空间中可容纳的分离点数量。

完整证明

(a) 以 $z_i$ 为中心、半径 $1/2$ 的球两两不交。又所有中心都在 $B(z_1,2)$ 中,所以这些小球都包含在 $B(z_1,5/2)$ 中。体积比较得

$$N(1/2)^n\le(5/2)^n,$$

即 $N\le5^n$。

(b) 取 $x_i=e_i\in\mathbb R^N$。若存在 $T$ 满足题中距离保持,则 $z_i=T(e_i)/(0.99\sqrt2)$ 在 $\mathbb R^n$ 中满足 $1<\|z_i-z_j\|_2\le2$。由 (a),$N\le5^n$。当 $n<\frac12\log N$ 时,$5^n<N$,矛盾。因此不存在这样的 $T$。

Exercise 5.16矩阵 Taylor series
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:验证矩阵函数定义与多项式、幂级数和指数展开一致。

证明思路

对谱分解逐特征值运算。

完整证明

设 $X=\sum_i\lambda_i u_iu_i^{\mathsf T}$。则 $X^k=\sum_i\lambda_i^k u_iu_i^{\mathsf T}$。因此对多项式 $f(x)=\sum_{k=0}^pa_kx^k$,

$$f(X)=\sum_i f(\lambda_i)u_iu_i^{\mathsf T}=\sum_{k=0}^pa_kX^k.$$

对收敛幂级数同理,矩阵级数在 Frobenius 范数中逐特征值绝对收敛,从而也在算子范数中收敛。取 $f(x)=e^x$ 即得 $e^X=I+X+X^2/2!+\cdots$。

Exercise 5.17矩阵单调性
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:说明交换矩阵可逐特征值比较,非交换时弱增函数未必矩阵单调。

证明思路

交换对称矩阵可同时正交对角化;反例用 $f(t)=t^2$。

完整证明

(a) 若 $X,Y$ 交换且对称,则可同时正交对角化:$X=U\operatorname{diag}(x_i)U^{\mathsf T}$,$Y=U\operatorname{diag}(y_i)U^{\mathsf T}$。$X\preceq Y$ 等价于 $x_i\le y_i$。若 $f$ 弱增,则 $f(x_i)\le f(y_i)$,故 $f(X)\preceq f(Y)$。

(b) 令

$$X=\begin{bmatrix}1&0\\0&0\end{bmatrix},\qquad C=\begin{bmatrix}1&1\\1&1\end{bmatrix},\qquad Y=X+C.$$

则 $C\succeq0$,所以 $X\preceq Y$,且 $X,Y\succeq0$。函数 $f(t)=t^2$ 在 $[0,\infty)$ 上弱增。但

$$Y^2-X^2=\begin{bmatrix}4&3\\3&2\end{bmatrix}$$

行列式为 $8-9=-1$,不是半正定。因此 $X^2\preceq Y^2$ 失败。

Exercise 5.18$1/x$ 与 $\log x$ 的矩阵单调性
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明逆函数反向矩阵单调,$\log$ 正向矩阵单调。

证明思路

逆函数用 congruence 归一化;$\log$ 用积分表示和逆函数单调性。

完整证明

由 $X$ 可逆且 $0\preceq X$,有 $X\succ0$。又 $Y\succeq X$,故 $Y\succ0$。令 $A=Y^{-1/2}XY^{-1/2}$,则 $0\prec A\preceq I$,所以 $A^{-1}\succeq I$。两边作 congruence 变换,得到 $X^{-1}\succeq Y^{-1}$,且 $Y^{-1}\succeq0$。

积分恒等式由

$$\int_0^M\left(\frac1{1+t}-\frac1{x+t}\right)dt=\log(1+M)-\log(1)-\log(x+M)+\log x$$

令 $M\to\infty$ 得到。若 $0\prec X\preceq Y$,则对每个 $t\ge0$,$X+tI\preceq Y+tI$,由逆函数反向单调,$(X+tI)^{-1}\succeq(Y+tI)^{-1}$。代入积分表示并积分,得到 $\log X\preceq\log Y$。

Exercise 5.19矩阵指数
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明交换时指数可拆,非交换时一般不可拆。

证明思路

交换时同时对角化;反例用两个不交换的 Pauli 型矩阵。

完整证明

(a) 若 $X,Y$ 交换且对称,则可同时正交对角化。逐特征值使用 $e^{x_i+y_i}=e^{x_i}e^{y_i}$,得到 $e^{X+Y}=e^Xe^Y$。

(b) 取

$$X=\begin{bmatrix}0&1\\1&0\end{bmatrix},\qquad Y=\begin{bmatrix}1&0\\0&-1\end{bmatrix}.$$

它们不交换。$e^{X+Y}$ 是对称矩阵;而

$$e^Xe^Y=(\cosh1\,I+\sinh1\,X)(\cosh1\,I+\sinh1\,Y)$$

由于 $XY\ne YX$,该乘积不是对称矩阵。因此 $e^{X+Y}\ne e^Xe^Y$。

Exercise 5.20矩阵 Bernstein 期望界
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:从矩阵 Bernstein 高概率尾界推出期望界。

证明思路

把尾界写成 $2n\exp[-c\min(t^2/\sigma^2,t/K)]$,选择阈值吸收维度因子,再积分。

完整证明

令 $Z=\|\sum_iX_i\|$。矩阵 Bernstein 给出

$$\mathbb P\{Z\ge t\}\le2n\exp[-c\min(t^2/\sigma^2,t/K)].$$

取 $u=\log(2n)$。该尾界等价推出

$$Z\lesssim \sigma\sqrt{u+s}+K(u+s)$$

的失败概率至多 $e^{-s}$。积分 $s\ge0$ 得

$$\mathbb EZ\lesssim\sigma\sqrt{1+\log n}+K(1+\log n).$$

这里 $\sigma^2=\|\sum_i\mathbb EX_i^2\|$。

Exercise 5.21Matrix Hoeffding
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明 Rademacher 矩阵和的 Hoeffding 型尾界。

证明思路

沿用矩阵 Bernstein 证明,只把单项 MGF 换成 $\mathbb Ee^{\lambda\varepsilon A}\preceq e^{\lambda^2A^2/2}$。

完整证明

对标量有 $\mathbb Ee^{\lambda\varepsilon x}=\cosh(\lambda x)\le e^{\lambda^2x^2/2}$。由函数演算升级为

$$\mathbb Ee^{\lambda\varepsilon A}\preceq e^{\lambda^2A^2/2}.$$

重复 Theorem 5.4.1 的 Lieb-MGF 证明,得到

$$\mathbb P\{\lambda_{\max}(S)\ge t\}\le n\exp\left(-\lambda t+\lambda^2\sigma^2/2\right),$$

其中 $S=\sum_i\varepsilon_iA_i$,$\sigma^2=\|\sum_iA_i^2\|$。优化 $\lambda=t/\sigma^2$ 得右尾。对 $-S$ 重复并合并,得到因子 $2n$。

Exercise 5.22Matrix Khintchine
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:由 Matrix Hoeffding 尾界推出 $L^p$ 范数界。

证明思路

对尾界积分,维度因子贡献 $\sqrt{\log n}$,$p$ 阶矩贡献 $\sqrt p$。

完整证明

令 $Z=\|\sum_i\varepsilon_iA_i\|$。Theorem 5.4.13 给出

$$\mathbb P\{Z\ge t\}\le2n\exp(-t^2/(2\sigma^2)).$$

等价地,$Z$ 的 subgaussian 尺度至多 $C\sigma$,但带有初始质量因子 $n$。用尾积分公式

$$\mathbb EZ^p=p\int_0^\infty t^{p-1}\mathbb P\{Z\ge t\}dt$$

并从 $t\asymp\sigma\sqrt{\log n}$ 开始积分,得到

$$\|Z\|_{L^p}\le C\sigma\sqrt{p+\log n}.$$
Exercise 5.23矩形矩阵 Bernstein
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:用 Hermitian dilation 把对称矩阵 Bernstein 推广到矩形矩阵。

证明思路

把 $X_i$ 嵌入对称块矩阵 $\mathcal D(X_i)$,其范数等于 $\|X_i\|$。

完整证明

定义

$$\mathcal D(X_i)=\begin{bmatrix}0&X_i\\X_i^{\mathsf T}&0\end{bmatrix}.$$

则 $\mathcal D(X_i)$ 是 $(m+n)\times(m+n)$ 对称矩阵,均值为零,且 $\|\mathcal D(X_i)\|=\|X_i\|\le K$。同时

$$\mathcal D(X_i)^2=\begin{bmatrix}X_iX_i^{\mathsf T}&0\\0&X_i^{\mathsf T}X_i\end{bmatrix}.$$

所以 variance 参数至多题中给出的 $\sigma^2$。对 $\sum_i\mathcal D(X_i)=\mathcal D(\sum_iX_i)$ 应用矩阵 Bernstein,并用 $\|\mathcal D(M)\|=\|M\|$,得到结论。

Exercise 5.24矩形矩阵 Khintchine
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明矩形 Rademacher 矩阵和的 Khintchine 型矩估计。

证明思路

对 $A_i$ 做 Hermitian dilation,然后应用对称矩阵 Khintchine。

完整证明

令 $\mathcal D(A_i)$ 如 Exercise 5.23。则

$$\left\|\sum_i\varepsilon_iA_i\right\|=\left\|\sum_i\varepsilon_i\mathcal D(A_i)\right\|,$$

并且

$$\left\|\sum_i\mathcal D(A_i)^2\right\|=\max\left(\left\|\sum_iA_iA_i^{\mathsf T}\right\|,\left\|\sum_iA_i^{\mathsf T}A_i\right\|\right)\le\sigma^2.$$

对维度 $m+n$ 的对称矩阵 Khintchine 应用 Theorem 5.4.14,即得

$$\left(\mathbb E\left\|\sum_i\varepsilon_iA_i\right\|^p\right)^{1/p}\le C\sigma\sqrt{p+\log(m+n)}.$$
Exercise 5.25无 loops 的 SBM
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明去掉 loops 不改变 Theorem 5.5.1 的谱聚类保证。

证明思路

无 loops 只删除对角项;随机噪声矩阵的 Bernstein 参数不增,信号矩阵只受有界对角扰动。

完整证明

定义无 loops 模型为 $A_{ii}=0$,$i<j$ 时按社群内/社群间概率 $p,q$ 独立连边,并令 $A_{ji}=A_{ij}$。随机部分分解为 $\sum_{i<j}Z_{ij}$,比含 loops 情况少了 $i=j$ 项,因此 $\|Z_{ij}\|\le1$ 且 $\sigma^2\le a$ 仍成立。

期望矩阵 $D$ 与含 loops 版本只差一个对角矩阵,其范数至多 $p\le a/n$,远小于谱间隙量级 $a-b$。把这一扰动并入 Davis-Kahan 的噪声项,常数调整后 Theorem 5.5.1 的条件仍给出 $99\%$ 准确率。

Exercise 5.26协方差估计高概率版
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明 Remark 5.6.5 的高概率界。

证明思路

在 Theorem 5.6.1 的证明中保留矩阵 Bernstein 的尾界,而不是积分成期望。

完整证明

沿 Theorem 5.6.1,$Y_i=X_iX_i^{\mathsf T}-\Sigma$ 满足

$$\sigma^2\le K^2m\operatorname{tr}(\Sigma)\|\Sigma\|,\qquad M\le2K^2\operatorname{tr}(\Sigma).$$

矩阵 Bernstein 给出以至少 $1-2e^{-u}$ 的概率,

$$\left\|\sum_iY_i\right\|\le C\left(\sigma\sqrt{\log n+u}+M(\log n+u)\right).$$

除以 $m$,并写 $r=\operatorname{tr}(\Sigma)/\|\Sigma\|$,得到

$$\|\Sigma_m-\Sigma\|\le C\left(\sqrt{\frac{K^2r(\log n+u)}m}+\frac{K^2r(\log n+u)}m\right)\|\Sigma\|.$$
Exercise 5.27有界性假设不能去掉
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:构造各向同性但大量为零的分布,说明没有有界性时样本可能完全看不到协方差。

证明思路

让 $X$ 以很高概率为 $0$,以很小概率取巨大值以维持协方差为 $I$。

完整证明

取任意小 $\delta>0$。令 $Y$ 是各向同性随机向量,例如 $Y$ 均匀取 $\pm\sqrt n e_j$,并定义

$$X=\begin{cases}0,&1-\delta,\\ \delta^{-1/2}Y,&\delta.\end{cases}$$

则 $\mathbb EXX^{\mathsf T}=\mathbb EYY^{\mathsf T}=I$,所以 $X$ 各向同性。但 $m$ 个样本全为零的概率是 $(1-\delta)^m$。若取 $\delta\ll1/m$,该概率接近 $1$;此时 $\Sigma_m=0$,而 $\Sigma=I$,误差为 $1$。因此没有 (5.20) 这种有界性或其他尾部条件时,不可能有统一协方差估计保证。

Exercise 5.28对数因子不可去掉
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:用 coupon collector 型例子说明一般协方差估计需要 $n\log n$ 样本。

证明思路

取 $X=\sqrt n\,\varepsilon e_J$。若某个坐标没被抽到,样本协方差在该方向为零,误差至少为 $1$。

完整证明

令 $J$ 均匀分布在 $\{1,\ldots,n\}$,$\varepsilon$ 为 Rademacher,$X=\sqrt n\,\varepsilon e_J$。则 $\mathbb EXX^{\mathsf T}=I$,且 $\|X\|_2=\sqrt n=(\mathbb E\|X\|_2^2)^{1/2}$,满足 Theorem 5.6.1 的有界性条件。

若 $m$ 个样本没有覆盖某个坐标 $k$,则 $\Sigma_m e_k=0$,而 $\Sigma e_k=e_k$,所以 $\|\Sigma_m-\Sigma\|\ge1$。coupon collector 说明若 $m\le c n\log n$,以高概率存在未出现的坐标。因此 $\|\Sigma_m-\Sigma\|<\|\Sigma\|$ 高概率失败。

若矩阵 Bernstein 期望界中的 $\log n$ 可一般去掉,Theorem 5.6.1 的证明会给出 $m\asymp n$ 的一般协方差估计,与上面的例子矛盾。因此对数因子一般不可去掉。

Exercise 5.29有效秩
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:理解有效秩与实际 rank、稳定性和低维支撑的关系。

证明思路

把有效秩写成特征值和除以最大特征值。

完整证明

(a) 若 $\Sigma\ne0$,设非零特征值为 $\lambda_i$,最大值为 $\lambda_{\max}$。则

$$1\le\frac{\sum_i\lambda_i}{\lambda_{\max}}\le\frac{\operatorname{rank}(\Sigma)\lambda_{\max}}{\lambda_{\max}}=\operatorname{rank}(\Sigma)\le n.$$

(b) 取 $\Sigma=\operatorname{diag}(1,\epsilon,\ldots,\epsilon)$,它满秩,但 $r(\Sigma)=1+(n-1)\epsilon\to1$。(c) 在 $\|\Sigma\|>0$ 的集合上,trace 和 operator norm 都连续,且分母非零,所以商连续。(d) 若 $X$ 取值于 $k$ 维子空间 $E$,则 $\Sigma$ 的像空间包含在 $E$ 中,因此 $\operatorname{rank}(\Sigma)\le k$,从而 $r(\Sigma)\le k$。

Exercise 5.30从 frames 采样
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明等范数 Parseval frame 随机采样后仍近似保持 frame operator。

证明思路

把 frame elements 看成一个有界各向同性分布,应用协方差估计。

完整证明

设 frame 为 $u_1,\ldots,u_M$,满足 $\sum_j u_ju_j^{\mathsf T}=I$,且等范数给出 $\|u_j\|_2^2=n/M$。令随机向量 $X=\sqrt M\,u_J$,其中 $J$ 均匀分布在 $\{1,\ldots,M\}$。则

$$\mathbb EXX^{\mathsf T}=\sum_j u_ju_j^{\mathsf T}=I,\qquad \|X\|_2=\sqrt n.$$

抽取 $m$ 个 frame elements 对应样本协方差

$$\Sigma_m=\frac Mm\sum_{\ell=1}^m u_{J_\ell}u_{J_\ell}^{\mathsf T}.$$

由 Theorem 5.6.1 的高概率形式,若 $m\gtrsim n\log n$,则 $\|\Sigma_m-I\|\le c$ 高概率成立。这就是近似 frame:其 frame operator 在 Loewner order 下满足 $(1-c)I\preceq\Sigma_m\preceq(1+c)I$。

Exercise 5.31一般独立 rows 的随机矩阵
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:在 row 有界但不要求次高斯时,证明随机矩阵奇异值仍有 $K\sqrt{n\log n}$ 级偏差控制。

证明思路

把 $\frac1mA^{\mathsf T}A-I$ 看成样本协方差偏差,并用高概率协方差估计控制。

完整证明

rows $A_i$ 各向同性意味着 $\mathbb EA_iA_i^{\mathsf T}=I_n$,所以

$$\frac1mA^{\mathsf T}A=\frac1m\sum_{i=1}^m A_iA_i^{\mathsf T}$$

是样本协方差。由 $\|A_i\|_2\le K\sqrt n$ 和 $\mathbb E\|A_i\|_2^2=n$,Theorem 5.6.1 的高概率版本给出

$$\left\|\frac1mA^{\mathsf T}A-I\right\|\le C Kt\sqrt{\frac{n\log n}{m}}$$

以至少 $1-2n^{-ct^2}$ 的概率成立;当右侧超过常数时,所需下界可能为负且上界可由 $\|A\|^2\le\sum_i\|A_i\|^2\le mK^2n$ 粗略吸收,调整常数后仍成立。

记 $\delta=CKt\sqrt{n\log n/m}$。由奇异值与 Gram 矩阵关系,

$$1-\delta\le s_n(A/\sqrt m)\le s_1(A/\sqrt m)\le1+\delta.$$

乘以 $\sqrt m$,得到题中不等式。

Exercise 5.32Matrix sketching
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明均匀抽样 rows 可近似保持 tall matrix 的奇异值平方。

证明思路

把随机 row 看成有界分布,样本协方差就是缩放后的 $B^{\mathsf T}B$。

完整证明

设 $a_1,\ldots,a_N$ 是 $A$ 的 rows,且 $\|a_j\|_2$ 相等。令随机向量 $X=\sqrt N\,a_J$,$J$ 均匀分布在 $\{1,\ldots,N\}$。则

$$\mathbb EXX^{\mathsf T}=A^{\mathsf T}A.$$

抽取 $m$ 次形成 $B$,样本二阶矩为

$$\Sigma_m=\frac1m\sum_{\ell=1}^mX_\ell X_\ell^{\mathsf T}=\frac Nm B^{\mathsf T}B.$$

等 row norm 给出有界性;有效秩不超过 $n$。由 Theorem 5.6.1 的高概率版本,若 $m\ge Cn\log n$,则以至少 $0.9$ 的概率,

$$\left\|\frac NmB^{\mathsf T}B-A^{\mathsf T}A\right\|\le0.1\|A^{\mathsf T}A\|=0.1s_1(A)^2.$$

Weyl 不等式把矩阵谱范数偏差转成所有特征值偏差:

$$\max_i\left|\frac Nm s_i(B)^2-s_i(A)^2\right|\le0.1s_1(A)^2.$$

易混点

易混点 正确读法
“无独立性集中”不是完全没有独立性 5.1-5.3 的球面集中不靠坐标独立;5.4 的矩阵 Bernstein 仍要求矩阵 summands 独立。
Lipschitz 集中不是所有函数集中 Lipschitz 条件排除高频振荡;没有它可任意放大坏事件指标函数。
JL 的随机性不依赖数据 随机投影先选好后同时适用于给定有限点集。
矩阵 Bernstein 不是 entrywise Bernstein 它控制算子范数,summand 内部 entries 不必独立。
有效秩不是 algebraic rank 它是特征值能量分散程度,近似低秩时可远小于实际 rank。

学习检查表

  • [ ] 能用一句话解释等周不等式、blow-up 和 Lipschitz 集中的关系。
  • [ ] 能复现 Theorem 5.1.3 的中位数 sublevel set 证明。
  • [ ] 能说明 Johnson-Lindenstrauss 中 $\log N$ 来自哪里。
  • [ ] 能解释矩阵 Bernstein 证明中 Lieb inequality 解决了什么问题。
  • [ ] 能计算 $Z_{ij}$ 的范数和稀疏 SBM 的 variance 参数。
  • [ ] 能复现 Theorem 5.6.1 中 $\sigma^2$ 与 $M$ 的估计。
  • [ ] 能说明协方差估计中为什么会出现 $n\log n$,以及有效秩如何替换 $n$。

后续衔接

第 5 章会在后续两条线上反复出现:一条是随机投影和 sketching,后面会在第 9 章继续深化;另一条是矩阵集中,后面会在随机过程、chaining 和随机矩阵在集合上的偏差中继续作为基础工具。复习时应把 Theorem 5.4.1 和 Theorem 5.6.1 与第 4 章的 net argument 区分清楚:第 4 章靠离散化所有方向,第 5 章靠矩阵 Laplace transform 直接控制谱范数。