第 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 直接控制谱范数。