第 4 章学习笔记:随机矩阵

一句话定位

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

本章导读

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

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

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

本页使用方式

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

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

本章主线

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

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

本章学习路线

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

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

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

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

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

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

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

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

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

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

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

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

分层阅读路线

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

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

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

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

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

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

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

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

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

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

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

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

样本协方差是

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

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

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

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

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

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

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

高斯混合模型可以写成

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

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

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

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

跳到 Gaussian mixture 完整证明

核心对象与符号表

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

关键定理卡片

Theorem 4.1.1 SVD

任意矩阵都可分解为

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

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

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

谱子空间扰动满足

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

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

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

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

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

以高概率成立。

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

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

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

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

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

证明模板:net argument

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

关键定理完整证明

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

证明思路

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

完整证明

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

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

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

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

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

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

因此作为算子,

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

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

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

证明思路

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

完整证明

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

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

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

于是

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

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

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

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

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

合并三式得到

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

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

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

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

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

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

最后直接计算

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

于是

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

证明思路

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

完整证明

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

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

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

$$ \|A\|\le2\max_{x\in\mathcal N,\ y\in\mathcal M}|\langle Ax,y\rangle|. $$

固定 $x,y$ 时,

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

是独立次高斯变量的线性组合。由独立次高斯和的范数估计,

$$ \|\langle Ax,y\rangle\|_{\psi_2}^2 \le C\sum_{i,j}x_j^2y_i^2\|A_{ij}\|_{\psi_2}^2 \le CK^2. $$

因此对所有 $u\ge0$,

$$ \mathbb P\{|\langle Ax,y\rangle|>u\}\le2e^{-cu^2/K^2}. $$

对 $\mathcal N\times\mathcal M$ 做并集界,得到

$$ \mathbb P\left\{ \max_{x\in\mathcal N,y\in\mathcal M} |\langle Ax,y\rangle|>u \right\} \le 2\cdot9^{m+n}e^{-cu^2/K^2}. $$

取 $u=C_0K(\sqrt m+\sqrt n+t)$,并把 $C_0$ 取为足够大的绝对常数,就有

$$ 2\cdot9^{m+n}e^{-cu^2/K^2} \le 2e^{-t^2}. $$

因此在概率至少 $1-2e^{-t^2}$ 的事件上,

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

这就是 Theorem 4.4.3。

Complete Proof Theorem 4.6.1:奇异值双侧界
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明各向同性次高斯行矩阵在所有方向上近似保持长度。

证明思路

对固定方向 $x$,$\|Ax\|_2^2$ 是独立次指数变量之和,Bernstein 给出集中。再对球面 net 做 union bound,并用逼近把 net 控制推广到所有方向。

完整证明

设 $A$ 的行 $A_i$ 独立、各向同性、次高斯,且 $K=\max_i\|A_i\|_{\psi_2}$。先证明更强的矩阵偏差界

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

以至少 $1-2e^{-t^2}$ 的概率成立。

取 $S^{n-1}$ 的 $1/4$-net $\mathcal N$,满足 $|\mathcal N|\le9^n$。由二次型 net 逼近引理,

$$ \left\| \frac1m A^{\mathsf T}A-I_n \right\| \le 2\max_{x\in\mathcal N} \left| \frac1m\|Ax\|_2^2-1 \right|. $$

下面固定 $x\in\mathcal N$。随机变量 $X_i=\langle A_i,x\rangle$ 独立,并由各向同性满足 $\mathbb E X_i^2=1$,且 $\|X_i\|_{\psi_2}\le K$。因此

$$ \|Ax\|_2^2=\sum_{i=1}^m\langle A_i,x\rangle^2 $$

且 $X_i^2-1$ 是独立中心化次指数变量,$\|X_i^2-1\|_{\psi_1}\le CK^2$。令

$$ \eta=K^2\max(\delta,\delta^2). $$

Bernstein 不等式给出

$$ \begin{aligned} \mathbb P\left\{ \left|\frac1m\|Ax\|_2^2-1\right| \ge \frac\eta2 \right\} &\le 2\exp\left[ -c_1m\min\left(\frac{\eta^2}{K^4},\frac{\eta}{K^2}\right) \right]\\ &= 2\exp(-c_1m\delta^2). \end{aligned} $$

等号处使用了 $\eta/K^2=\max(\delta,\delta^2)$,因此最小值为 $\delta^2$。又由 $\delta$ 的定义,$m\delta^2\ge C^2(n+t^2)$。对 $\mathcal N$ 做 union bound,

$$ \mathbb P\left\{ \max_{x\in\mathcal N} \left| \frac1m\|Ax\|_2^2-1 \right| \ge \frac\eta2 \right\} \le 9^n\cdot2e^{-c_1C^2(n+t^2)} \le 2e^{-t^2} $$

只要 $C$ 足够大。结合 net 逼近引理,得到上述矩阵偏差界。

最后把矩阵偏差转成奇异值界。矩阵 $m^{-1}A^{\mathsf T}A$ 的特征值是 $s_i(A)^2/m$。若

$$ \left\|\frac1m A^{\mathsf T}A-I_n\right\|\le\eta, $$

则所有 $i$ 满足

$$ 1-\eta\le \frac{s_i(A)^2}{m}\le1+\eta. $$

开方并使用 $\eta=K^2\max(\delta,\delta^2)$(下界为负时结论平凡),得到

$$ \sqrt m-CK^2(\sqrt n+t) \le s_n(A) \le s_1(A) \le \sqrt m+CK^2(\sqrt n+t). $$
Complete Proof Theorem 4.7.1:协方差估计
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:在一般协方差 $\Sigma$ 下,用白化和随机矩阵双侧界证明样本协方差的期望算子范数误差。

证明思路

先把 $X$ 白化为各向同性向量 $Z=\Sigma^{-1/2}X$。白化后的样本协方差偏差由 Theorem 4.6.1 的期望形式控制;再乘回 $\Sigma^{1/2}$,算子范数最多多出 $\|\Sigma\|$。

完整证明

设 $X_1,\ldots,X_m$ 是 $X$ 的独立样本,$\Sigma=\mathbb E XX^{\mathsf T}$,样本协方差为

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

先假设 $\Sigma$ 可逆。令 $Z=\Sigma^{-1/2}X$,$Z_i=\Sigma^{-1/2}X_i$。则 $Z_i$ 独立、各向同性,并且由假设

$$ \|\langle Z,x\rangle\|_{\psi_2} = \|\langle X,\Sigma^{-1/2}x\rangle\|_{\psi_2} \le K\|\langle X,\Sigma^{-1/2}x\rangle\|_{L^2} = K\|x\|_2, $$

所以 $\|Z_i\|_{\psi_2}\le K$。记

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

把 $Z_i^{\mathsf T}$ 作为行组成 $m\times n$ 矩阵 $A$,则 $R_m=m^{-1}A^{\mathsf T}A-I_n$。由 Theorem 4.6.1 的期望形式,

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

另一方面,$X_i=\Sigma^{1/2}Z_i$,所以

$$ \Sigma_m-\Sigma = \Sigma^{1/2}R_m\Sigma^{1/2}. $$

于是

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

取期望得到

$$ \mathbb E\|\Sigma_m-\Sigma\| \le CK^2\left(\sqrt{\frac nm}+\frac nm\right)\|\Sigma\|. $$

若 $\Sigma$ 不可逆,可在 $\operatorname{range}(\Sigma)$ 上做同样白化,或用 $\Sigma+\varepsilon I$ 近似后令 $\varepsilon\downarrow0$;结论不变。

章节阅读路径

小节 读法
4.1 重点掌握 SVD、算子范数、min-max、Weyl、Davis-Kahan。
4.2 把 covering/packing 当作“连续集合离散化”的工具。
4.3 通过 Hamming cube 理解 metric entropy。
4.4 学会随机矩阵范数上界的 net argument。
4.5 看懂“期望矩阵 + 噪声矩阵 + 扰动理论”的社区检测套路。
4.6 理解双侧奇异值界为何意味着近似等距。
4.7 把双侧界用于协方差估计和聚类。

易混点

易混点 正确读法
算子范数只是最大元素大小 错。算子范数是所有单位方向上的最大拉伸,通常比最大元素更全局。
固定 $x$ 的集中等于所有 $x$ 同时集中 不等于。必须支付 net size 的复杂度成本。
net argument 只是技术细节 它是第 4 章从标量/向量集中走向矩阵范数的核心桥梁。
Davis-Kahan 只需要噪声小 还需要信号谱间隙足够大;误差尺度是 $\|E\|/\text{gap}$。
协方差估计只是大数定律 高维中要控制所有方向,所以需要样本矩阵奇异值或算子范数界。

公式卡片

工具 公式 用途
SVD $A=\sum_i s_i u_i v_i^{\mathsf T}$ 分解矩阵作用
算子范数 $\|A\|=\sup_{x\in S^{n-1}}\|Ax\|_2$ 最大拉伸
Frobenius $\|A\|_F^2=\sum_i s_i^2$ 总能量
Net 大小 $\mathcal N(B_2^n,\varepsilon B_2^n)\le(3/\varepsilon)^n$ union bound 成本
随机矩阵范数 $\|A\|\lesssim\sqrt m+\sqrt n$ 次高斯矩阵
Davis-Kahan $\text{误差}\lesssim\|E\|/\text{gap}$ 谱方法稳定性

正文隐藏验证补全

Hidden Check 4.1.1:SVD 证明中 $m\ge n$ 用在哪里
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:解释原书在证明 SVD 时为什么先假设 $m\ge n$,以及该假设实际用于哪一步。

完整证明

证明从 $A^{\mathsf T}A$ 的谱分解得到 $n$ 个右奇异向量 $v_1,\ldots,v_n\in\mathbb R^n$。当 $s_i>0$ 时,左奇异向量必须定义为

$$ u_i=\frac{Av_i}{s_i}, $$

这些向量自动正交归一。问题出现在 $s_i=0$ 的方向:这时 $Av_i=0$,等式 $Av_i=s_i u_i$ 对任意单位向量 $u_i$ 都成立,但我们还要让全部 $u_1,\ldots,u_n$ 正交归一。

因此,原证明中真正用到 $m\ge n$ 的地方,是“任选 $u_i$,只要保持正交归一性”这一步。只有当 $\mathbb R^m$ 至少有 $n$ 维时,才一定能把已得到的左奇异向量补成 $n$ 个正交归一向量。

如果 $m<n$,不能在 $\mathbb R^m$ 中选出 $n$ 个正交归一向量。此时可以对 $A^{\mathsf T}$ 应用 $n\ge m$ 的情形,或直接从 $AA^{\mathsf T}$ 出发先构造 $m$ 个左奇异向量。最终得到的紧 SVD 始终只需要 $r=\min(m,n)$ 个方向。

Hidden Check 4.2.14:Hamming 距离确实是度量
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:验证 $d_H(x,y)=\#\{i:x(i)\ne y(i)\}$ 满足度量的三个条件。

完整证明

非负性来自 $d_H(x,y)$ 是一个集合的基数。并且 $d_H(x,y)=0$ 当且仅当没有坐标不同,也就是 $x=y$。

对称性也直接成立:$x(i)\ne y(i)$ 与 $y(i)\ne x(i)$ 是同一件事,所以 $d_H(x,y)=d_H(y,x)$。

最后证明三角不等式。对每个坐标 $i$,若 $x(i)\ne z(i)$,则不可能同时有 $x(i)=y(i)$ 且 $y(i)=z(i)$。因此指标函数逐坐标满足

$$ \mathbf 1_{\{x(i)\ne z(i)\}} \le \mathbf 1_{\{x(i)\ne y(i)\}} + \mathbf 1_{\{y(i)\ne z(i)\}}. $$

对 $i=1,\ldots,n$ 求和,得到

$$ d_H(x,z)\le d_H(x,y)+d_H(y,z). $$

所以 $d_H$ 是 $\{0,1\}^n$ 上的度量。

Hidden Check 4.6:由 (4.27) 推出奇异值界 (4.26)
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:补全正文中“利用 (4.15) 可以验证 (4.27) 推出 (4.26)”的代数步骤。

完整证明

把矩阵缩放为

$$ B=\frac1{\sqrt m}A. $$

那么

$$ B^{\mathsf T}B-I_n = \frac1mA^{\mathsf T}A-I_n, $$

且 $B$ 的奇异值是 $s_i(A)/\sqrt m$。由 (4.27),

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

可把 $K$ 替换为 $\max(K,1)$,这只会放宽结论。于是令 $\Delta=K^2\delta$,有

$$ K^2\max(\delta,\delta^2) \le \max(\Delta,\Delta^2). $$

现在应用 (4.15) 到矩阵 $B$,得到所有奇异值满足

$$ 1-\Delta \le s_n(B) \le s_1(B) \le 1+\Delta. $$

乘以 $\sqrt m$,并代入 $\Delta$,得到

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

若左端为负,该下界当然仍成立;这就是 (4.26)。

Exercises 4.1-4.51 证明工作区

批次 建议目标
4.1-4.22 线性代数、扰动理论、operator norms、cut norm、SDP。
4.23-4.33 nets、packing、covering、体积、Hamming cube。
4.34-4.44 net argument 和随机矩阵范数变体。
4.45-4.51 社区检测、协方差估计、低秩矩阵 covering、Gaussian mixture 聚类。

本区先把每道习题的证明入口固定下来,方便从译文、学习笔记和并排阅读页跳转。当前 Exercises 4.1-4.51 均已升级为“证明目标 + 证明思路 + 完整证明”。

Exercise 4.1 逆矩阵的 SVD
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:说明可逆性等价于所有奇异值非零,并写出 $A^{-1}$ 的 SVD 形式。

证明思路

SVD 把 $A$ 写成 $U\Sigma V^{\mathsf T}$。矩阵可逆当且仅当对角矩阵 $\Sigma$ 可逆;逆矩阵就是 $V\Sigma^{-1}U^{\mathsf T}$。

完整证明

由 SVD,$A=U\Sigma V^{\mathsf T}$,其中 $\Sigma=\operatorname{diag}(s_1,\ldots,s_n)$,$U,V$ 正交。由于正交矩阵可逆,$A$ 可逆当且仅当 $\Sigma$ 可逆,也就是当且仅当所有 $s_i\ne0$。

在这种情况下,

$$ A^{-1} = V\Sigma^{-1}U^{\mathsf T} = \sum_{i=1}^n s_i^{-1}v_i u_i^{\mathsf T}. $$
Exercise 4.2 算子范数的基本性质
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:验证算子范数的范数性质、转置不变性、次乘性和子矩阵单调性。

证明思路

始终回到定义 $\|A\|=\sup_{\|x\|_2=1}\|Ax\|_2$。转置不变性可用最大奇异值或双线性形式。

完整证明

(a) 非负性和齐次性直接来自 Euclidean 范数。若 $\|A\|=0$,则 $Ax=0$ 对所有单位向量成立,故 $A=0$。三角不等式为

$$ \|A+B\| = \sup_{\|x\|=1}\|(A+B)x\| \le \sup_{\|x\|=1}(\|Ax\|+\|Bx\|) \le \|A\|+\|B\|. $$

(b) $\|A\|=s_1(A)$,而 $A$ 与 $A^{\mathsf T}$ 有相同奇异值,所以 $\|A^{\mathsf T}\|=\|A\|$。

(c) 对任意单位向量 $x$,

$$ \|ABx\|\le\|A\|\|Bx\|\le\|A\|\|B\|. $$

取上确界得 $\|AB\|\le\|A\|\|B\|$。例子:$A=\operatorname{diag}(1,0)$,$B=\operatorname{diag}(0,1)$,则 $AB=0$ 但 $\|A\|\|B\|=1$。

(d) 子矩阵可写成 $PAQ$,其中 $P,Q$ 是坐标投影,且 $\|P\|,\|Q\|\le1$。所以 $\|PAQ\|\le\|A\|$。

Exercise 4.3 算子范数:简单例子
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:计算 rank-one 矩阵和对角矩阵的算子范数。

证明思路

rank-one 矩阵只把向量沿 $v$ 方向的分量送到 $u$ 方向;对角矩阵逐坐标缩放。

完整证明

(a) 对任意单位向量 $x$,

$$ uv^{\mathsf T}x=u\langle v,x\rangle, \qquad \|uv^{\mathsf T}x\|\le\|u\|\|v\|. $$

取 $x=v/\|v\|$ 达到等号,所以 $\|uv^{\mathsf T}\|=\|u\|\|v\|$。Frobenius 范数满足

$$ \|uv^{\mathsf T}\|_F^2 = \sum_{ij}u_i^2v_j^2 = \|u\|^2\|v\|^2. $$

(b) 若 $A=\operatorname{diag}(a_i)$,则

$$ \|Ax\|_2^2=\sum_i a_i^2x_i^2\le(\max_i|a_i|)^2\|x\|_2^2. $$

取 $x=e_j$,其中 $|a_j|=\max_i|a_i|$,得到等号。

Exercise 4.4 算子范数与 Frobenius 范数
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明 Frobenius 范数与算子范数、各向同性向量和左乘矩阵之间的基本关系。

证明思路

用奇异值表达 Frobenius 范数;用 trace 表达各向同性向量下的二阶矩;用 $\|Bz\|\le\|B\|\|z\|$ 控制左乘。

完整证明

(a) 若 $\operatorname{rank}(A)=r$,则只有 $r$ 个非零奇异值。因此

$$ \|A\|=s_1\le\left(\sum_i s_i^2\right)^{1/2}=\|A\|_F \le \sqrt r\,s_1=\sqrt r\,\|A\|. $$

rank-one 矩阵使左侧等号成立;$r$ 个相等非零奇异值使右侧等号成立。

(b) 若 $Z$ 各向同性,$\mathbb EZZ^{\mathsf T}=I$,则

$$ \mathbb E\|AZ\|_2^2 = \mathbb E\,Z^{\mathsf T}A^{\mathsf T}AZ = \operatorname{tr}(A^{\mathsf T}A\mathbb EZZ^{\mathsf T}) = \operatorname{tr}(A^{\mathsf T}A) = \|A\|_F^2. $$

(c) 设 $A_1,\ldots,A_n$ 是 $A$ 的列,或直接逐列看 $BA$。有

$$ \|BA\|_F^2 = \sum_j\|B A_j\|_2^2 \le \|B\|^2\sum_j\|A_j\|_2^2 = \|B\|^2\|A\|_F^2. $$

(d) 若 $A$ 对角,则 $BA$ 的第 $j$ 列是 $a_jB_j$,所以

$$ \|BA\|_F^2 = \sum_j a_j^2\|B_j\|_2^2 \le \|B\|_{1\to2}^2\sum_j a_j^2 = \|B\|_{1\to2}^2\|A\|_F^2. $$
Exercise 4.5 奇异值的一个界
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明第 $k$ 个奇异值由 Frobenius 范数控制。

证明思路

Frobenius 范数平方是所有奇异值平方之和,而前 $k$ 个奇异值都至少为 $s_k$。

完整证明

奇异值按降序排列,因此

$$ \|A\|_F^2 = \sum_i s_i(A)^2 \ge \sum_{i=1}^k s_i(A)^2 \ge k\,s_k(A)^2. $$

开平方得到 $s_k(A)\le \|A\|_F/\sqrt k$。

Exercise 4.6 Power method
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明 power method 的极限给出算子范数。

证明思路

在右奇异向量基下展开 Gaussian 初始向量。最高奇异值对应的分量几乎必然非零,反复作用会让最大奇异值主导。

完整证明

设 $A^{\mathsf T}A=V\operatorname{diag}(s_i^2)V^{\mathsf T}$,并写

$$ x=\sum_i\alpha_i v_i. $$

因为 $x$ 是标准 Gaussian,任意固定方向上的系数 $\alpha_i$ 几乎必然非零。于是

$$ (A^{\mathsf T}A)^kx = \sum_i \alpha_i s_i^{2k}v_i. $$

令 $s_1=\|A\|$。若 $s_1>0$,则

$$ \|(A^{\mathsf T}A)^kx\|_2 = s_1^{2k} \left(\sum_i\alpha_i^2(s_i/s_1)^{4k}\right)^{1/2}. $$

括号中的量在 $2k$ 次根下趋于 $1$,因为至少一个 top singular direction 的系数非零。因此

$$ \sqrt[2k]{\|(A^{\mathsf T}A)^kx\|_2}\to s_1=\|A\|. $$

若 $A=0$,结论平凡。

Exercise 4.7 算子范数与列范数
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:用列范数给出算子范数下界,并刻画等号情形。

证明思路

列向量就是 $Ae_i$。若某个列向量已经达到算子范数,则 $e_i$ 是 $A^{\mathsf T}A$ 的最大特征向量,迫使该列与其余列正交。

完整证明

(a) 对每个 $i$,

$$ \|A\|\ge\|Ae_i\|_2=\|A_i\|_2. $$

若列两两正交,则

$$ \|Ax\|_2^2 = \left\|\sum_i x_iA_i\right\|_2^2 = \sum_i x_i^2\|A_i\|_2^2 \le \max_i\|A_i\|_2^2\|x\|_2^2, $$

所以等号成立。

(b) 若 $\|A\|=\|A_i\|_2$,则 $e_i$ 使 Rayleigh quotient $x^{\mathsf T}A^{\mathsf T}Ax$ 取到最大值。因此 $e_i$ 是 $A^{\mathsf T}A$ 的特征向量。于是对 $j\ne i$,

$$ \langle A_i,A_j\rangle = e_i^{\mathsf T}A^{\mathsf T}Ae_j = 0. $$

(c) 对行应用同样结论于 $A^{\mathsf T}$。

Exercise 4.8 Schur bound
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明 Schur bound:算子范数由最大行绝对和与最大列绝对和控制。

证明思路

对每一行使用 Cauchy-Schwarz 的加权版本,再交换求和,把行和、列和分别提出来。

完整证明

$$ R=\max_i\sum_j|a_{ij}|,\qquad C=\max_j\sum_i|a_{ij}|. $$

对任意 $x$,

$$ |(Ax)_i|^2 = \left|\sum_j a_{ij}x_j\right|^2 \le \left(\sum_j|a_{ij}|\right) \left(\sum_j|a_{ij}||x_j|^2\right) \le R\sum_j|a_{ij}||x_j|^2. $$

对 $i$ 求和得到

$$ \|Ax\|_2^2 \le R\sum_j\left(\sum_i|a_{ij}|\right)|x_j|^2 \le RC\|x\|_2^2. $$

取单位向量上的上确界,$\|A\|\le\sqrt{RC}$。

Exercise 4.9 Walsh 矩阵
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明归一化 Walsh 矩阵是正交矩阵。

证明思路

用归纳法检查 $W_kW_k^{\mathsf T}=nI_n$,其中 $n$ 是矩阵阶数。

完整证明

当 $k=1$ 时,

$$ W_1W_1^{\mathsf T}=2I_2. $$

假设 $W_kW_k^{\mathsf T}=nI_n$。则

$$ W_{k+1}W_{k+1}^{\mathsf T} = \begin{bmatrix} W_k&W_k\\ W_k&-W_k \end{bmatrix} \begin{bmatrix} W_k^{\mathsf T}&W_k^{\mathsf T}\\ W_k^{\mathsf T}&-W_k^{\mathsf T} \end{bmatrix} = \begin{bmatrix} 2nI_n&0\\ 0&2nI_n \end{bmatrix}. $$

所以若 $W_{k+1}$ 阶数为 $2n$,则 $W_{k+1}W_{k+1}^{\mathsf T}=(2n)I_{2n}$。归纳完成,故 $W_k/\sqrt n$ 正交。

Exercise 4.10 具有 $\pm1$ 元素的矩阵
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明 $\pm1$ 矩阵算子范数介于 $\sqrt n$ 与 $n$ 之间,并说明两端可达。

证明思路

上界由 Frobenius 范数给出;下界由 Frobenius 范数和秩给出。全 $1$ 矩阵达到上界,Walsh 矩阵达到下界。

完整证明

因为 $A$ 有 $n^2$ 个元素且每个绝对值为 $1$,

$$ \|A\|_F=n. $$

上界由 $\|A\|\le\|A\|_F$ 得到:$\|A\|\le n$。下界由 Exercise 4.4(a) 得到:

$$ \|A\|_F\le\sqrt{\operatorname{rank}(A)}\,\|A\| \le \sqrt n\,\|A\|, $$

所以 $\|A\|\ge\sqrt n$。

全 $1$ 矩阵秩为 $1$,唯一非零奇异值为 $n$,达到上界。对 $n=2^k$,Walsh 矩阵满足 $WW^{\mathsf T}=nI_n$,所有奇异值都是 $\sqrt n$,因此达到下界。

Exercise 4.11 正交投影之差
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明任意两个正交投影之差的算子范数不超过 $1$。

证明思路

$P-Q$ 是自伴矩阵。对自伴矩阵,算子范数由 Rayleigh quotient 控制;而 $\langle(P-Q)x,x\rangle=\|Px\|^2-\|Qx\|^2$ 位于 $[-1,1]$。

完整证明

因为 $P,Q$ 是正交投影,$P=P^{\mathsf T}=P^2$,$Q=Q^{\mathsf T}=Q^2$,所以 $P-Q$ 自伴。对任意单位向量 $x$,

$$ \langle(P-Q)x,x\rangle = \langle Px,x\rangle-\langle Qx,x\rangle = \|Px\|_2^2-\|Qx\|_2^2. $$

两个平方范数都在 $[0,1]$ 中,因此该差属于 $[-1,1]$。由自伴矩阵的变分刻画,

$$ \|P-Q\| = \sup_{\|x\|_2=1}|\langle(P-Q)x,x\rangle| \le1. $$
Exercise 4.12 正交投影之差与乘积
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:用互补投影刻画两个正交投影之差的范数。

证明思路

先把 $P-Q$ 写成两个正交分量之和,再利用秩相等时 $PQ$ 与 $QP$ 的非零奇异值相同。

完整证明

先证明 (a)。恒等式

$$ P-Q=PQ_\perp-P_\perp Q $$

来自 $Q_\perp=I-Q$ 和 $P_\perp=I-P$。对任意单位向量 $x$,两个向量 $PQ_\perp x$ 与 $P_\perp Qx$ 分别落在 $\operatorname{range}(P)$ 与 $\operatorname{range}(P)^\perp$,所以正交。因此

$$ \|(P-Q)x\|_2^2 = \|PQ_\perp x\|_2^2+\|P_\perp Qx\|_2^2 \le M^2\bigl(\|Q_\perp x\|_2^2+\|Qx\|_2^2\bigr) = M^2, $$

其中 $M=\max(\|PQ_\perp\|,\|P_\perp Q\|)$。这给出 $\|P-Q\|\le M$。反向不等式也直接成立,因为

$$ PQ_\perp=(P-Q)Q_\perp, \qquad P_\perp Q=-P_\perp(P-Q), $$

而投影的算子范数不超过 $1$。故

$$ \|P-Q\| = \max(\|P_\perp Q\|,\|PQ_\perp\|). $$

下面证明 (b)。若 $\operatorname{rank}(P)>\operatorname{rank}(Q)$,则

$$ \dim\operatorname{range}(P)+\dim\ker(Q)>n, $$

所以存在单位向量 $x\in\operatorname{range}(P)\cap\ker(Q)$。于是 $(P-Q)x=x$,故 $\|P-Q\|\ge1$。结合 Exercise 4.11 的 $\|P-Q\|\le1$,得到 $\|P-Q\|=1$。若 $\operatorname{rank}(Q)>\operatorname{rank}(P)$,同理。

最后证明 (c)。令两个秩都为 $k$。有

$$ \|P_\perp Q\|^2 = \|QP_\perp Q\|, \qquad \|PQ_\perp\|^2 = \|P Q_\perp P\|. $$

在 $\operatorname{range}(Q)$ 上,$QPQ$ 的非零特征值等于算子 $PQ$ 的奇异值平方;在 $\operatorname{range}(P)$ 上,$PQP$ 的非零特征值也是同一批数。由于两个空间维数都为 $k$,零特征值的重数也相同。因此 $QP_\perp Q=Q-QPQ$ 与 $PQ_\perp P=P-PQP$ 有相同特征值,从而两个范数相等。结合 (a),得到

$$ \|P-Q\|=\|P_\perp Q\|=\|PQ_\perp\|. $$
Exercise 4.13 谱投影的 Davis-Kahan
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明 top $k$ 谱投影的 Davis-Kahan 型界。

证明思路

把 top $k$ 特征值看成一个谱簇。只要该谱簇与其余谱相隔 $\lambda_k-\lambda_{k+1}$,Davis-Kahan 定理直接控制两个谱投影。

完整证明

设 $\delta=\lambda_k-\lambda_{k+1}>0$。若 $2\|A-B\|\ge\delta$,右侧至少为 $1$,而 $\|P_A-P_B\|\le1$,结论平凡。

否则令 $\varepsilon=\|A-B\|<\delta/2$。由 Weyl 不等式,对 $j>k$ 有

$$ \mu_j\le\lambda_j+\varepsilon\le\lambda_{k+1}+\varepsilon. $$

而 $A$ 的 top $k$ 谱簇位于 $[\lambda_k,\infty)$ 中,所以该谱簇与 $B$ 的 bottom 谱簇 $\{\mu_j:j>k\}$ 的距离至少

$$ \lambda_k-(\lambda_{k+1}+\varepsilon) = \delta-\varepsilon \ge \delta/2. $$

把 Theorem 4.1.15 证明中的谱投影分离引理用于 $A$ 的 top $k$ 投影 $P_A$ 和 $B$ 的 bottom 投影 $I-P_B$,得到

$$ \|(I-P_B)P_A\| \le \frac{\varepsilon}{\delta/2} = \frac{2\|A-B\|}{\delta}. $$

由于 $P_A$ 与 $P_B$ 秩相同,Exercise 4.12(c) 给出 $\|P_A-P_B\|=\|(I-P_B)P_A\|$。因此

$$ \|P_A-P_B\| \le \frac{2\|A-B\|}{\lambda_k-\lambda_{k+1}}. $$
Exercise 4.14 Hermitian dilation
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明 Hermitian dilation 的非零特征值是原矩阵奇异值的正负号版本。

证明思路

把 SVD 中的左右奇异向量拼成分块向量,直接作用分块矩阵。

完整证明

若 $Av_i=s_iu_i$ 且 $A^{\mathsf T}u_i=s_iv_i$,则

$$ H\begin{bmatrix}u_i\\ v_i\end{bmatrix} = \begin{bmatrix}Av_i\\ A^{\mathsf T}u_i\end{bmatrix} = s_i\begin{bmatrix}u_i\\ v_i\end{bmatrix}, $$

并且

$$ H\begin{bmatrix}u_i\\ -v_i\end{bmatrix} = \begin{bmatrix}-Av_i\\ A^{\mathsf T}u_i\end{bmatrix} = -s_i\begin{bmatrix}u_i\\ -v_i\end{bmatrix}. $$

因此每个非零奇异值 $s_i$ 产生两个特征值 $\pm s_i$。剩余维度对应 $A$ 或 $A^{\mathsf T}$ 的核,特征值为 $0$。

Exercise 4.15 Wedin theorem
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:把矩形矩阵的奇异子空间扰动界化为 Hermitian dilation 的 Davis-Kahan 定理。

证明思路

构造 $H_A$ 和 $H_B$。它们的谱投影差由 Davis-Kahan 控制;左奇异子空间投影只是该分块谱投影的左上角压缩,所以范数不会增大。

完整证明

$$ H_A=\begin{bmatrix}0&A\\ A^{\mathsf T}&0\end{bmatrix}, \qquad H_B=\begin{bmatrix}0&B\\ B^{\mathsf T}&0\end{bmatrix}. $$

由 Exercise 4.14,$H_A$ 的正特征值是 $s_i(A)$,对应特征向量为 $[u_i;v_i]$ 的归一化版本。并且

$$ \|H_A-H_B\|=\|A-B\|. $$

令 $R_A$ 为 $H_A$ 对谱簇 $\{\pm s_1,\ldots,\pm s_k\}$ 的谱投影,并令 $R_B$ 为 $H_B$ 对对应谱簇的谱投影。这个投影同时包含正、负两组特征向量;对每个 $i\le k$,两个归一化特征向量 $[u_i;v_i]/\sqrt2$ 与 $[u_i;-v_i]/\sqrt2$ 的投影相加为

$$ \begin{bmatrix} u_i u_i^{\mathsf T} & 0\\ 0 & v_i v_i^{\mathsf T} \end{bmatrix}. $$

所以 $R_A$ 的左上角分块正是 $P_A$,右下角分块正是 top $k$ 右奇异向量投影。

若 $2\|A-B\|\ge s_k-s_{k+1}$,则右侧至少为 $1$,而两个正交投影之差的范数不超过 $1$,结论成立。否则 Weyl 不等式保证 $H_B$ 中对应的 $\pm t_1,\ldots,\pm t_k$ 谱簇仍与其余谱相隔至少 $(s_k-s_{k+1})/2$。对这个分离谱簇应用投影型 Davis-Kahan,得到

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

取左上角分块不会增加算子范数,于是 top $k$ 左奇异向量投影满足题中界。取右下角分块得到 top $k$ 右奇异向量投影的同样界。

单个第 $k$ 个奇异向量的版本作用于只含 $\pm s_k$ 的谱簇。若 $s_k$ 与相邻奇异值的间隔为

$$ \Delta_k=\min(s_{k-1}-s_k,\ s_k-s_{k+1}) $$

(端点处只保留存在的项),则投影误差至多 $2\|A-B\|/\Delta_k$;再用 Exercise 4.16 可把一维投影距离转为奇异向量在差一个全局符号意义下的 Euclidean 距离界。

Exercise 4.16 向量之间的夹角与距离
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:把小夹角转换成向量本身在差符号意义下的 Euclidean 距离。

证明思路

选取符号让内积非负。距离平方等于 $2-2|\langle u,v\rangle|$,再用 $\sin$ 控制。

完整证明

取 $\theta\in\{-1,1\}$ 使 $\langle u,\theta v\rangle=|\langle u,v\rangle|=\cos\angle(u,v)\ge0$。则

$$ \|u-\theta v\|_2^2 = 2-2\cos\angle(u,v). $$

对 $c\in[0,1]$ 有 $1-c\le1-c^2$。因此

$$ \|u-\theta v\|_2^2 \le 2(1-\cos^2\angle(u,v)) = 2\sin^2\angle(u,v) \le 2\varepsilon^2. $$

开平方得到结论。

Exercise 4.17 近似投影
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:把 isometry / approximate isometry 与 $AA^{\mathsf T}$ 接近投影联系起来。

证明思路

用 SVD。$A^{\mathsf T}A$ 和 $AA^{\mathsf T}$ 共享非零特征值;这些特征值就是奇异值平方。

完整证明

(a) 若 $A$ 是 isometry,则 $A^{\mathsf T}A=I_n$。矩阵 $AA^{\mathsf T}$ 对称,且

$$ (AA^{\mathsf T})^2=A(A^{\mathsf T}A)A^{\mathsf T}=AA^{\mathsf T}, $$

所以它是到 $\operatorname{range}(A)$ 的正交投影。反过来,若 $AA^{\mathsf T}$ 是秩为 $n$ 的正交投影,则其非零特征值全为 $1$。这些正是 $A^{\mathsf T}A$ 的特征值,所以 $A^{\mathsf T}A=I_n$,即 $A$ 是 isometry。

(b) 写完整的 $n$ 个右奇异方向版本:$A=U\Sigma V^{\mathsf T}$,其中 $U$ 是 $m\times n$ 的列正交矩阵;若 $A$ 有零奇异值,就把对应左奇异向量任意补成正交列。令 $P=UU^{\mathsf T}$,这是秩为 $n$ 的正交投影。则

$$ AA^{\mathsf T}=U\Sigma^2U^{\mathsf T}, \qquad AA^{\mathsf T}-P=U(\Sigma^2-I)U^{\mathsf T}. $$

因此

$$ \|AA^{\mathsf T}-P\| = \|\Sigma^2-I\| = \|A^{\mathsf T}A-I\|. $$

这正是 approximate isometry 的二次型表述,故两者等价。

Exercise 4.18 $p\to q$ 范数
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:建立 $p\to q$ operator norm 的双线性表达、范数性质和对偶关系。

证明思路

用 $\ell^q$ 与 $\ell^{q'}$ 的对偶性:$\|z\|_q=\max_{\|y\|_{q'}=1}|y^{\mathsf T}z|$。

完整证明

(a) 对任意 $x$,

$$ \|Ax\|_q = \max_{\|y\|_{q'}=1}|y^{\mathsf T}Ax|. $$

再对 $\|x\|_p=1$ 最大化,得到 (4.31)。

(b) 非负性、齐次性和三角不等式都由 $\|Ax\|_q/\|x\|_p$ 的定义直接推出;若范数为 $0$,则 $Ax=0$ 对所有 $x$ 成立,所以 $A=0$。

(c) 由双线性表达,

$$ \|A^{\mathsf T}\|_{p\to q} = \max_{\|x\|_p=1,\ \|y\|_{q'}=1} |y^{\mathsf T}A^{\mathsf T}x| = \max_{\|y\|_{q'}=1,\ \|x\|_p=1} |x^{\mathsf T}Ay| = \|A\|_{q'\to p'}. $$
Exercise 4.19 $1\to\infty$ 和 $1\to2$ 范数
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:给出 $1\to\infty$、$1\to2$、$2\to\infty$ 范数的元素级表达。

证明思路

$\ell^1$ 单位球的极点是 $\pm e_j$;$\ell^\infty$ 输出看行,$\ell^2$ 输出看列。

完整证明

(a) 对任意 $\|x\|_1\le1$,

$$ |(Ax)_i| \le \sum_j |A_{ij}||x_j| \le \max_j|A_{ij}|. $$

再对 $i$ 取最大,得 $\|A\|_{1\to\infty}\le\max_{ij}|A_{ij}|$。取 $x=\pm e_j$ 达到最大元素,故等号成立。

(b) 设 $A_j$ 为第 $j$ 列。若 $\|x\|_1\le1$,则

$$ \|Ax\|_2 = \left\|\sum_jx_jA_j\right\|_2 \le \sum_j|x_j|\|A_j\|_2 \le \max_j\|A_j\|_2. $$

取 $x=e_j$ 得等号。因此 $\|A\|_{1\to2}=\max_j\|A_{:j}\|_2$。对 $2\to\infty$,第 $i$ 个坐标是行向量 $A_{i:}$ 与 $x$ 的内积,Cauchy-Schwarz 给出上界,取 $x$ 平行于最大行得到等号。

Exercise 4.20 $\infty\to1$ 范数
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明 $\infty\to1$ 范数可在符号向量上取到,并给出 rank-one 对偶表述和二次型反例。

证明思路

线性函数在盒子上最大值取在顶点。rank-one 表述把 $y^{\mathsf T}Ax$ 写成矩阵内积。二次型反例来自非凸二次函数在盒子内部达到更大绝对值。

完整证明

(a) 由 $\ell^1$ 与 $\ell^\infty$ 对偶,

$$ \|A\|_{\infty\to1} = \max_{\|x\|_\infty\le1,\ \|y\|_\infty\le1}|y^{\mathsf T}Ax|. $$

该双线性函数在两个盒子上最大,固定一边时另一边的最大值都在顶点取得,所以可取 $x,y\in\{-1,1\}$。

(b) 对秩一矩阵 $Z=yx^{\mathsf T}$,有 $\langle A,Z\rangle=y^{\mathsf T}Ax$。若 $\|x\|_\infty,\|y\|_\infty\le1$,则 $|Z_{ij}|\le1$。反过来,任意秩一 $Z=uv^{\mathsf T}$ 且 $|Z_{ij}|\le1$,可重新缩放为 $Z=yx^{\mathsf T}$,其中 $\|x\|_\infty,\|y\|_\infty\le1$。于是得到所需上确界。

(c) 取

$$ A= \begin{bmatrix} -1&0\\ 0&2 \end{bmatrix}. $$

在顶点 $x\in\{-1,1\}^2$ 上,$x^{\mathsf T}Ax=1$,所以最大绝对值为 $1$。但在 $x=(0,1)$ 处,$x^{\mathsf T}Ax=2$。因此二次型版本不能只检查符号顶点。

Exercise 4.21 Cut norm
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明 cut norm 与 $\infty\to1$ 范数在常数因子内等价。

证明思路

cut norm 是用 $0/1$ 指示向量选择子矩阵;$\infty\to1$ 范数是用 $\pm1$ 符号向量。两者可通过把符号向量分解为正负集合互相转换。

完整证明

对任意 $I,J$,令 $y=\mathbf 1_I$、$x=\mathbf 1_J$,它们的 $\ell^\infty$ 范数不超过 $1$。因此

$$ \left|\sum_{i\in I,j\in J}A_{ij}\right| = |y^{\mathsf T}Ax| \le \|A\|_{\infty\to1}. $$

所以 $\|A\|_{\mathrm{cut}}\le\|A\|_{\infty\to1}$。

反过来,由 Exercise 4.20,只需考虑 $x\in\{-1,1\}^n$、$y\in\{-1,1\}^m$。令 $I_\pm=\{i:y_i=\pm1\}$,$J_\pm=\{j:x_j=\pm1\}$。则 $y^{\mathsf T}Ax$ 是四个矩形和的带符号组合:

$$ S(I_+,J_+)-S(I_+,J_-)-S(I_-,J_+)+S(I_-,J_-). $$

每个矩形和的绝对值都不超过 cut norm,因此

$$ |y^{\mathsf T}Ax|\le4\|A\|_{\mathrm{cut}}. $$

对符号向量取最大,即得 $\|A\|_{\infty\to1}\le4\|A\|_{\mathrm{cut}}$。

Exercise 4.22 SDP 松弛
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:说明若干难算矩阵范数可由 SDP 在常数因子内近似。

证明思路

核心是把符号或盒约束优化松弛为向量 Gram 矩阵优化。Grothendieck 不等式保证 $\infty\to1$ 的 SDP 松弛只损失常数;其他范数通过等价关系或平方化归约到同类 SDP。

完整证明

对 $\|A\|_{\infty\to1}$,Exercise 4.20 把它写成符号双线性优化。把符号 $x_j,y_i$ 放松为单位向量 $u_j,v_i$,目标变为

$$ \max \sum_{ij}A_{ij}\langle v_i,u_j\rangle, $$

这是 Gram 矩阵变量上的 SDP。Grothendieck 不等式说明这个 SDP 值与原值只差绝对常数因子。

cut norm 与 $\infty\to1$ 范数由 Exercise 4.21 在常数因子内等价,因此也可常数近似。

对于 $\|A\|_{\infty\to2}$,有

$$ \|A\|_{\infty\to2}^2 = \max_{x\in\{-1,1\}^n}x^{\mathsf T}A^{\mathsf T}Ax $$

在常数因子意义下可由正半定二次 Grothendieck/SDP 松弛近似。最后 $\|A\|_{2\to1}=\|A^{\mathsf T}\|_{\infty\to2}$,所以同样可由 SDP 常数近似。

Exercise 4.23 $\varepsilon$-net 的传递性
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明 net 的覆盖误差满足三角不等式式的传递。

证明思路

对 $K$ 中任意点,先用 $\mathcal M$ 近似,再用 $\mathcal N$ 近似 $\mathcal M$ 中的点。两次误差相加。

完整证明

任取 $x\in K$。因为 $\mathcal M$ 是 $K$ 的 $\delta$-net,存在 $m\in\mathcal M$ 使

$$ d(x,m)\le\delta. $$

又因为 $\mathcal N$ 是 $\mathcal M$ 的 $\varepsilon$-net,存在 $y\in\mathcal N$ 使

$$ d(m,y)\le\varepsilon. $$

由三角不等式,

$$ d(x,y)\le d(x,m)+d(m,y)\le\delta+\varepsilon. $$

这对任意 $x\in K$ 成立,所以 $\mathcal N$ 是 $K$ 的 $(\varepsilon+\delta)$-net。

Exercise 4.24 Packing balls
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:说明 separated set 与半径 $\varepsilon/2$ balls 两两不交的逆命题在一般 metric space 中失败,但在 normed space 中成立。

证明思路

一般 metric space 可能没有两个点之间的“中点”。normed space 有中点;若两中心距离不超过 $\varepsilon$,中点会同时落入两个闭球。

完整证明

(a) 取 $T=\{a,b\}$,定义 $d(a,b)=0.75\varepsilon$。以 $a,b$ 为中心的闭 $\varepsilon/2$-balls 分别为 $\{a\}$ 和 $\{b\}$,因为 $0.75\varepsilon>\varepsilon/2$,所以它们不交。但 $d(a,b)<\varepsilon$,因此 $\{a,b\}$ 不是 $\varepsilon$-separated。这给出一般 metric space 中的反例。

(b) 设 $T$ 是 normed space 的子集,且以 $\mathcal N$ 中点为中心的闭 $\varepsilon/2$-balls 两两不交。若存在不同 $x,y\in\mathcal N$ 满足 $\|x-y\|\le\varepsilon$,则中点 $z=(x+y)/2$ 满足

$$ \|z-x\|=\|z-y\|=\frac12\|x-y\|\le\frac\varepsilon2. $$

于是 $z$ 同时落入两个闭 $\varepsilon/2$-balls,矛盾。因此任意不同中心距离都大于 $\varepsilon$,即 $\mathcal N$ 是 $\varepsilon$-separated。

Exercise 4.25 External covering numbers
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明允许球心在 $K$ 外部不会本质改变 covering number,只会改变半径常数。

证明思路

内部 covering 是 external covering 的特殊情形,所以外部覆盖数不比内部覆盖数大。反向时,把每个外部中心移动到它所覆盖到的一个 $K$ 中点上,半径从 $\varepsilon/2$ 增大到 $\varepsilon$。

完整证明

第一部分直接成立:普通 covering 要求中心在 $K$ 中,是 external covering 的特殊情形,所以

$$ \overline{\mathcal N}(K,d,\varepsilon)\le \mathcal N(K,d,\varepsilon). $$

下面证明第二个不等式。设 $y_1,\ldots,y_M$ 是 $K$ 的 external $\varepsilon/2$-net,其中

$$ M=\overline{\mathcal N}(K,d,\varepsilon/2). $$

丢掉那些半径 $\varepsilon/2$ 球与 $K$ 没有交集的中心。对每个剩余的 $y_i$,选取一个 $x_i\in K$ 满足 $d(x_i,y_i)\le\varepsilon/2$。若 $z\in K$ 被 $y_i$ 覆盖,即 $d(z,y_i)\le\varepsilon/2$,则

$$ d(z,x_i)\le d(z,y_i)+d(y_i,x_i)\le\varepsilon. $$

所以这些 $x_i\in K$ 构成 $K$ 的普通 $\varepsilon$-net,且数量不超过 $M$。因此

$$ \mathcal N(K,d,\varepsilon)\le\overline{\mathcal N}(K,d,\varepsilon/2). $$
Exercise 4.26 Covering number 的单调性
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:说明 covering number 关于半径单调,但关于集合包含只满足带半径损失的近似单调性。

证明思路

反例来自“内部球心”的限制:大集合可能包含一个中心点,可以同时覆盖小集合的两端;小集合自身没有这个中心点。近似单调性用 Exercise 4.25 的移动中心技巧。

完整证明

关于 $\varepsilon$ 的单调性是直接的:半径越大,每个球覆盖越多,所以 $\varepsilon_1\le\varepsilon_2$ 时

$$ \mathcal N(K,d,\varepsilon_2)\le \mathcal N(K,d,\varepsilon_1). $$

(a) 取实线中的 $K=\{-1,0,1\}$,$L=\{-1,1\}$,并令 $\varepsilon=1.1$。集合 $K$ 可由一个以 $0\in K$ 为中心的 $\varepsilon$-ball 覆盖,所以 $\mathcal N(K,d,\varepsilon)=1$。但 $L$ 的中心必须在 $L$ 中;以 $-1$ 或 $1$ 为中心的半径 $1.1$ 球都不能覆盖另一个点,因为两点距离为 $2$。因此 $\mathcal N(L,d,\varepsilon)=2$,这与普通单调性相反。

(b) 设 $L\subset K$。取 $K$ 的一个 $\varepsilon/2$-net,中心为 $x_1,\ldots,x_M\in K$,其中 $M=\mathcal N(K,d,\varepsilon/2)$。只保留那些 $\varepsilon/2$-球与 $L$ 相交的中心;对每个保留中心 $x_i$,选一个 $y_i\in L$ 满足 $d(x_i,y_i)\le\varepsilon/2$。若 $z\in L$ 被 $x_i$ 覆盖,则

$$ d(z,y_i)\le d(z,x_i)+d(x_i,y_i)\le\varepsilon. $$

于是这些 $y_i\in L$ 构成 $L$ 的 $\varepsilon$-net,数量不超过 $M$。所以

$$ \mathcal N(L,d,\varepsilon)\le\mathcal N(K,d,\varepsilon/2). $$
Exercise 4.27 Euclidean ball 的体积:几何论证
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:用 simplex、$\ell^1$ ball 和 cube 的几何关系给出 Euclidean ball 体积的粗界。

证明思路

先算 canonical simplex 的体积,再把 $\ell^1$ ball 分成 $2^n$ 个 simplex。Euclidean ball 夹在一个 cube 和缩放的 $\ell^1$ ball 之间。

完整证明

(a) 对 $\Delta_n=\{x_i\ge0,\ \sum_i x_i\le1\}$,按最后一个坐标积分:

$$ \operatorname{Vol}(\Delta_n) = \int_0^1 \operatorname{Vol}\{x_1,\ldots,x_{n-1}\ge0,\ \sum_{i=1}^{n-1}x_i\le1-t\}\,dt. $$

由归纳假设,该截面体积为 $(1-t)^{n-1}/(n-1)!$。因此

$$ \operatorname{Vol}(\Delta_n) = \frac1{(n-1)!}\int_0^1(1-t)^{n-1}\,dt = \frac1{n!}. $$

(b) $B_1^n$ 被 $2^n$ 个正交象限分割;每个象限中的部分都是一个 canonical simplex。因此

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

由 $n!\ge(n/e)^n$,得到 $\operatorname{Vol}(B_1^n)\le(2e/n)^n$。

(c) 若 $\|x\|_\infty\le1/\sqrt n$,则 $\|x\|_2\le1$,所以

$$ [-1/\sqrt n,1/\sqrt n]^n\subset B_2^n. $$

这给出下界 $\operatorname{Vol}(B_2^n)\ge(2/\sqrt n)^n$。另一方面,由 Cauchy-Schwarz,$\|x\|_1\le\sqrt n\|x\|_2$,所以

$$ B_2^n\subset \sqrt n\,B_1^n. $$

因此

$$ \operatorname{Vol}(B_2^n) \le n^{n/2}\operatorname{Vol}(B_1^n) \le n^{n/2}\left(\frac{2e}{n}\right)^n = \left(\frac{2e}{\sqrt n}\right)^n. $$
Exercise 4.28 Euclidean ball 的体积:概率论证
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:用 Gaussian 概率给出单位 Euclidean ball 体积的渐近尖锐上界。

证明思路

把 Gaussian 落入半径 $\sqrt n$ 球的概率用密度下界控制。概率不超过 $1$,于是反推出球体积上界。

完整证明

设 $g\sim N(0,I_n)$。在球 $\sqrt n B_2^n$ 上,Gaussian 密度满足

$$ (2\pi)^{-n/2}e^{-\|x\|_2^2/2} \ge (2\pi)^{-n/2}e^{-n/2}. $$

因此

$$ 1 \ge \mathbb P\{g\in\sqrt n B_2^n\} \ge (2\pi)^{-n/2}e^{-n/2}\operatorname{Vol}(\sqrt n B_2^n). $$

而 $\operatorname{Vol}(\sqrt n B_2^n)=n^{n/2}\operatorname{Vol}(B_2^n)$,所以

$$ \operatorname{Vol}(B_2^n) \le \left(\sqrt{\frac{2\pi e}{n}}\right)^n. $$
Exercise 4.29 Euclidean ball 的体积:解析论证
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明一般范数的径向积分公式,并由此精确计算 Euclidean ball 体积和单位体积半径渐近。

证明思路

范数球满足 $\operatorname{Vol}(tB)=t^n\operatorname{Vol}(B)$。把积分按层集写成 Stieltjes 积分,再分部积分。代入 Gaussian kernel 后使用 Gamma 函数和 Stirling 公式。

完整证明

(a) 因为 $\{x:\|x\|\le t\}=tB$,其体积为 $t^n\operatorname{Vol}(B)$。对递减可微且趋于 $0$ 的 $f$,层集积分给出

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

分部积分,边界项为零,于是

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

(b) 对 Euclidean norm 取 $f(t)=e^{-t^2/2}$。左侧为

$$ \int_{\mathbb R^n}e^{-\|x\|_2^2/2}\,dx=(2\pi)^{n/2}. $$

右侧中的积分为

$$ -\int_0^\infty t^n f'(t)\,dt = \int_0^\infty t^{n+1}e^{-t^2/2}\,dt = 2^{n/2}\Gamma(n/2+1). $$

所以

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

(c) 体积为 $1$ 的 Euclidean ball 半径 $R_n$ 满足

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

由 Stirling 公式 $\Gamma(n/2+1)=\sqrt{\pi n}(n/(2e))^{n/2}(1+o(1))$,得到

$$ \operatorname{Vol}(B_2^n)^{1/n} = \sqrt{\frac{2\pi e}{n}}\,(1+o(1)). $$

因此

$$ R_n=\operatorname{Vol}(B_2^n)^{-1/n} = \sqrt{\frac{n}{2\pi e}}\,(1+o(1)). $$
Exercise 4.30 $\ell^p$ ball 的体积
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:给出 $\ell^p$ 单位球体积的粗界和 Gamma 函数精确公式。

证明思路

粗界来自 cube 与缩放 $\ell^1$ ball 的夹逼。精确公式沿用 Exercise 4.29 的径向积分,取 $f(t)=e^{-t^p}$。

完整证明

(a) 若 $\|x\|_\infty\le n^{-1/p}$,则 $\|x\|_p\le1$,所以

$$ [-n^{-1/p},n^{-1/p}]^n\subset B_p^n. $$

这给出下界 $(2/n^{1/p})^n$。另一方面,Hölder 不等式给出 $\|x\|_1\le n^{1-1/p}\|x\|_p$,因此

$$ B_p^n\subset n^{1-1/p}B_1^n. $$

结合 $\operatorname{Vol}(B_1^n)\le(2e/n)^n$,得到

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

(b) 在 Exercise 4.29 的公式中取 $\|\cdot\|=\|\cdot\|_p$ 和 $f(t)=e^{-t^p}$。左侧分解为坐标积分:

$$ \int_{\mathbb R^n}e^{-\|x\|_p^p}\,dx = \left(2\int_0^\infty e^{-s^p}\,ds\right)^n = \bigl(2\Gamma(1+1/p)\bigr)^n. $$

右侧的径向积分为

$$ -\int_0^\infty t^n f'(t)\,dt = p\int_0^\infty t^{n+p-1}e^{-t^p}\,dt = \Gamma(n/p+1). $$

因此

$$ \operatorname{Vol}(B_p^n) = \frac{\bigl(2\Gamma(1/p+1)\bigr)^n}{\Gamma(n/p+1)}. $$
Exercise 4.31 Lattice 是 $\varepsilon$-net
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明缩放整数 lattice 在单位 Euclidean ball 中给出显式 $\varepsilon$-net,并估计其大小。

证明思路

对每个坐标向零取整,保证格点仍在单位球内,同时每个坐标误差不超过 $\varepsilon/\sqrt n$。大小估计用格点小立方体和 Euclidean ball 的体积上界。

完整证明

令格距 $h=\varepsilon/\sqrt n$。给定 $x\in B_2^n$,定义 $y_i$ 为与 $x_i$ 同号且满足

$$ |y_i|=h\left\lfloor\frac{|x_i|}{h}\right\rfloor. $$

则 $y\in h\mathbb Z^n$,并且 $|y_i|\le |x_i|$,所以 $\|y\|_2\le\|x\|_2\le1$,即 $y\in\mathcal N$。同时 $|x_i-y_i|\le h$,故

$$ \|x-y\|_2\le \sqrt n\,h=\varepsilon. $$

因此 $\mathcal N$ 是 $B_2^n$ 的 $\varepsilon$-net。这个构造也给出了快速近似算法:逐坐标向零取整即可。

下面估计大小。以每个 $y\in\mathcal N$ 为角放置半开小立方体 $y+[0,h)^n$;这些小立方体互不相交,并且每个点到 $y$ 的 Euclidean 距离至多 $\sqrt n h=\varepsilon$。因此它们都包含在 $(1+\varepsilon)B_2^n$ 中。于是

$$ |\mathcal N|h^n \le \operatorname{Vol}((1+\varepsilon)B_2^n). $$

用 Exercise 4.28 的体积上界,

$$ |\mathcal N| \le \left(\frac{\sqrt n}{\varepsilon}\right)^n (1+\varepsilon)^n \left(\sqrt{\frac{2\pi e}{n}}\right)^n = \left(\frac{\sqrt{2\pi e}(1+\varepsilon)}{\varepsilon}\right)^n. $$

对通常的 net 范围 $0<\varepsilon\le1$,右侧被 $e^n(2/\varepsilon+1)^n$ 控制;若 $\varepsilon>1$,该界更平凡成立。故得到题中估计。

Exercise 4.32 Hamming cube 的 covering 与 packing numbers
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明 Proposition 4.2.15 中 Hamming cube 的 covering 与 packing number 界。

证明思路

把 Euclidean 体积换成 Hamming ball 的基数。半径 $m$ 的 Hamming ball 大小为 $\sum_{k=0}^m\binom nk$。

完整证明

在 $\{0,1\}^n$ 中,任意中心的半径 $m$ Hamming ball 的基数为

$$ V_m=\sum_{k=0}^m\binom nk, $$

因为可以选择至多 $m$ 个坐标翻转。

若 $M$ 个半径 $m$ 球覆盖整个 cube,则 $M V_m\ge2^n$。因此

$$ \mathcal N(K,d_H,m)\ge\frac{2^n}{V_m}. $$

另一方面,取极大的 $m$-separated 子集。极大性保证它是一个 $m$-net,所以

$$ \mathcal N(K,d_H,m)\le\mathcal P(K,d_H,m). $$

最后,若 $\mathcal P$ 是一个 $m$-separated 集,则以其点为中心、半径 $\lfloor m/2\rfloor$ 的 Hamming balls 两两不交。每个这样的球大小为

$$ V_{\lfloor m/2\rfloor}=\sum_{k=0}^{\lfloor m/2\rfloor}\binom nk. $$

这些不交球都包含在 cube 中,因此

$$ |\mathcal P|V_{\lfloor m/2\rfloor}\le2^n. $$

取最大的 packing,即得

$$ \mathcal P(K,d_H,m) \le \frac{2^n}{\sum_{k=0}^{\lfloor m/2\rfloor}\binom nk}. $$
Exercise 4.33 Error correcting codes 的限制
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明能纠正 $r$ 个错误的 $n$ 位编码必须有至少 $r\log_2(n/r)$ 个冗余位。

证明思路

若 code 能纠正 $r$ 个错误,则不同 codewords 的半径 $r$ Hamming balls 必须两两不交。比较这些球的总基数与整个 Hamming cube 的基数。

完整证明

一个把 $k$-bit 字符串编码为 $n$-bit 字符串的 code 有 $2^k$ 个 codewords。若它能纠正 $r$ 个错误,则以不同 codeword 为中心、半径 $r$ 的 Hamming balls 必须两两不交;否则某个接收词会同时位于两个 codeword 的 $r$ 邻域内,无法唯一解码。

每个半径 $r$ Hamming ball 的大小为

$$ \sum_{j=0}^r\binom nj \ge \binom nr. $$

不交性给出 Hamming bound:

$$ 2^k\sum_{j=0}^r\binom nj\le2^n. $$

因此

$$ 2^{n-k}\ge \binom nr. $$

使用标准组合估计 $\binom nr\ge(n/r)^r$,得到

$$ n-k \ge \log_2\binom nr \ge r\log_2\left(\frac nr\right). $$

这就是所需限制。

Exercise 4.34 $\varepsilon$-net expansion
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明任意单位向量可按 $\varepsilon$-net 点作收敛展开,并用它重证 Lemma 4.4.1。

证明思路

反复用 net 点逼近当前残差方向。每一步残差长度最多乘以 $\varepsilon$,因此系数按 $\varepsilon^k$ 衰减。

完整证明

(a) 令 $r_0=x$。若 $r_k\ne0$,从 net 中取 $x_k\in\mathcal N$ 使

$$ \left\|\frac{r_k}{\|r_k\|_2}-x_k\right\|_2\le\varepsilon, \qquad \lambda_k=\|r_k\|_2, $$

并定义 $r_{k+1}=r_k-\lambda_kx_k$。则

$$ \|r_{k+1}\|_2 = \lambda_k \left\|\frac{r_k}{\|r_k\|_2}-x_k\right\|_2 \le \varepsilon\lambda_k. $$

由 $\lambda_0=1$ 归纳得 $0\le\lambda_k\le\varepsilon^k$,且 $\|r_k\|_2\to0$。因此

$$ x=\sum_{k=0}^{\infty}\lambda_kx_k. $$

(b) 令 $M=\sup_{z\in\mathcal N}\|Az\|_2$。对任意 $x\in S^{n-1}$ 使用 (a),

$$ \|Ax\|_2 \le \sum_{k=0}^{\infty}\lambda_k\|Ax_k\|_2 \le M\sum_{k=0}^{\infty}\varepsilon^k = \frac{M}{1-\varepsilon}. $$

对 $x$ 取上确界得到 Lemma 4.4.1 的上界;下界来自 $\mathcal N\subset S^{n-1}$,即在 net 上取上确界不会超过在整个球面上取上确界。

Exercise 4.35 在 $\varepsilon$-net 上计算范数
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:用球面 $\varepsilon$-net 近似计算向量范数。

证明思路

范数是 $\sup_{u\in S^{n-1}}\langle x,u\rangle$。用 net 点逼近最大方向 $u=x/\|x\|_2$。

完整证明

下界直接由 Cauchy-Schwarz 给出:

$$ \sup_{y\in\mathcal N}\langle x,y\rangle\le\|x\|_2. $$

若 $x=0$,上界两边都为 $0$。若 $x\ne0$,令 $u=x/\|x\|_2$。取 $y\in\mathcal N$ 使 $\|u-y\|_2\le\varepsilon$。则

$$ \|x\|_2 = \langle x,u\rangle = \langle x,y\rangle+\langle x,u-y\rangle \le \sup_{z\in\mathcal N}\langle x,z\rangle + \varepsilon\|x\|_2. $$

移项得

$$ \|x\|_2 \le \frac1{1-\varepsilon} \sup_{z\in\mathcal N}\langle x,z\rangle. $$
Exercise 4.36 在 $\varepsilon$-net 上最大化二次型
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明双线性型版本的 net 近似 Lemma 4.4.2。

证明思路

取几乎达到 $\|A\|$ 的方向 $x,y$,分别用 net 点 $x_0,y_0$ 逼近。两个逼近误差各贡献至多 $\varepsilon\|A\|$。

完整证明

$$ M=\sup_{x\in\mathcal N,\ y\in\mathcal M}|\langle Ax,y\rangle|. $$

下界 $M\le\|A\|$ 来自 $\mathcal N\subset S^{n-1}$、$\mathcal M\subset S^{m-1}$。为证上界,先假设 $0<\varepsilon<1/2$;否则分母 $1-2\varepsilon$ 不为正,该形式的估计没有内容。取 $x\in S^{n-1}$、$y\in S^{m-1}$ 使 $|\langle Ax,y\rangle|$ 任意接近 $\|A\|$。选择 $x_0\in\mathcal N$、$y_0\in\mathcal M$ 满足

$$ \|x-x_0\|_2\le\varepsilon, \qquad \|y-y_0\|_2\le\varepsilon. $$

$$ \begin{aligned} |\langle Ax,y\rangle| &\le |\langle Ax_0,y_0\rangle| + |\langle A(x-x_0),y\rangle| + |\langle Ax_0,y-y_0\rangle|\\ &\le M+2\varepsilon\|A\|. \end{aligned} $$

取上确界并整理,得到

$$ \|A\|\le\frac{M}{1-2\varepsilon}. $$

若 $A$ 对称且只用同一个 net,则从 $|\langle Ax,x\rangle|$ 的变分刻画出发,同样用 $x_0$ 逼近 $x$,误差仍至多 $2\varepsilon\|A\|$,得到相同结论。

Exercise 4.37 $\varepsilon$-net 上的范数偏差
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:把 net 上的长度偏差推广到整个球面。

证明思路

当 net 偏差已经大于中心量级时,Lemma 4.4.1 直接给出粗控制;否则先控制平方范数对应的二次型,再除回长度偏差。

完整证明

$$ \Delta=\sup_{x\in\mathcal N}\bigl|\|Ax\|_2-\mu\bigr|. $$

若 $\mu<0$,则 $|\|Ax\|_2-\mu|=\|Ax\|_2+|\mu|$。由 Lemma 4.4.1,

$$ \sup_{x\in S^{n-1}}\|Ax\|_2 \le \frac{1}{1-\varepsilon} \sup_{x\in\mathcal N}\|Ax\|_2 \le \frac{\Delta}{1-\varepsilon}, $$

且 $|\mu|\le\Delta$,所以结论成立。下面只考虑 $\mu\ge0$。

若 $\Delta\ge\mu$,则 $\sup_{\mathcal N}\|Ax\|_2\le\mu+\Delta\le2\Delta$。再次由 Lemma 4.4.1,

$$ \sup_{x\in S^{n-1}}\|Ax\|_2 \le \frac{2\Delta}{1-\varepsilon}, $$

于是 $\sup_{S^{n-1}}|\|Ax\|_2-\mu|\le C\Delta/(1-2\varepsilon)$。

剩下情形是 $0\le\Delta<\mu$。设

$$ B=A^{\mathsf T}A-\mu^2I_n. $$

对 $x\in\mathcal N$,

$$ |x^{\mathsf T}Bx| = \bigl|\|Ax\|_2^2-\mu^2\bigr| = \bigl|\|Ax\|_2-\mu\bigr|(\|Ax\|_2+\mu) \le \Delta(2\mu+\Delta) \le 3\mu\Delta. $$

由 Lemma 4.4.2 的对称情形,

$$ \|B\| \le \frac{3\mu\Delta}{1-2\varepsilon}. $$

因此对任意 $x\in S^{n-1}$,

$$ \bigl|\|Ax\|_2-\mu\bigr| = \frac{\bigl|\|Ax\|_2^2-\mu^2\bigr|}{\|Ax\|_2+\mu} \le \frac{\|B\|}{\mu} \le \frac{3\Delta}{1-2\varepsilon}. $$

合并三种情形即可得到题设不等式,其中 $C$ 为绝对常数。

Exercise 4.38 $\varepsilon$-net 的 convex hull
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明球面 $\varepsilon$-net 的凸包包含 $(1-\varepsilon)B_2^n$。

证明思路

对任意方向 $u$,net 中有点 $z$ 与 $u$ 接近,因此 $\langle u,z\rangle\ge1-\varepsilon$。这说明凸包在每个方向上的支撑函数至少为 $1-\varepsilon$。

完整证明

令 $C=\operatorname{conv}(\mathcal N)$。对任意 $u\in S^{n-1}$,取 $z\in\mathcal N$ 使 $\|u-z\|_2\le\varepsilon$。由

$$ \|u-z\|_2^2 = 2-2\langle u,z\rangle $$

可得 $\langle u,z\rangle\ge1-\varepsilon^2/2\ge1-\varepsilon$。因此

$$ \sup_{w\in C}\langle u,w\rangle\ge1-\varepsilon \qquad\text{for all }u\in S^{n-1}. $$

若存在 $y\in(1-\varepsilon)B_2^n$ 但 $y\notin C$,由超平面分离定理,存在 $u\in S^{n-1}$ 使

$$ \langle u,y\rangle> \sup_{w\in C}\langle u,w\rangle. $$

但 $\langle u,y\rangle\le\|y\|_2\le1-\varepsilon$,与上式的支撑函数下界矛盾。因此 $(1-\varepsilon)B_2^n\subset C$。

Exercise 4.39 随机点形成 $\varepsilon$-net
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明指数多个标准 Gaussian 点经 $\sqrt n$ 缩放后,以高概率在固定半径球内形成 $\varepsilon$-net。

证明思路

先取 $RB_2^n$ 的一个确定性细 net。对每个 net 中心,Gaussian 缩放点落入它附近的概率至少为 $e^{-C(R,\varepsilon)n}$。指数多个独立点加 union bound 即覆盖所有中心。

完整证明

令 $Z_i=g_i/\sqrt n$。取 $RB_2^n$ 的一个 $(\varepsilon/4)$-net $\mathcal M$,满足

$$ |\mathcal M|\le\left(\frac{CR}{\varepsilon}\right)^n. $$

对每个 $y\in\mathcal M$,把它稍微向内收缩为 $y'=(1-\varepsilon/(8R))y$。若 $\|Z_i-y'\|_2\le\varepsilon/8$,则 $Z_i\in RB_2^n$,且 $\|Z_i-y\|_2\le\varepsilon/4$。

Gaussian 密度在固定半径球内的指数阶下界给出

$$ \mathbb P\{\|Z_i-y'\|_2\le\varepsilon/8\} \ge \exp[-C(R,\varepsilon)n] =:p. $$

这里可以直接验证:在球 $B(y',\varepsilon/8)$ 中所有点范数不超过 $R$,所以 $Z_i$ 的密度至少为 $(n/2\pi)^{n/2}e^{-nR^2/2}$;再用 $\operatorname{Vol}(B_2^n)\ge(c/\sqrt n)^n$,得到上述指数型下界。

因此某个固定 $y\in\mathcal M$ 没有被任何样本点命中的概率至多 $(1-p)^N\le e^{-Np}$。由 union bound,所有 net 中心都被命中的失败概率至多

$$ |\mathcal M|e^{-Np}. $$

只要 $N\ge\exp(C'(R,\varepsilon)n)$,上式至多为 $e^{-cn}$。在好事件上,每个 $x\in RB_2^n$ 先找 $y\in\mathcal M$ 使 $\|x-y\|_2\le\varepsilon/4$,再找样本点 $Z_i\in RB_2^n$ 使 $\|Z_i-y\|_2\le\varepsilon/4$,于是 $\|x-Z_i\|_2\le\varepsilon$。这正是所需 net 性质。

Exercise 4.40 过多随机点不处于 convex position
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明当 Gaussian 点数指数级过大时,它们高概率不全是凸包顶点。

证明思路

把点分成两批。第一批足够多,缩放后在 $2B_2^n$ 内形成细 net,从而其凸包包含一个固定小球;第二批中高概率至少有一个点落入这个小球,于是该点在其他点的凸包内。

完整证明

缩放不改变 convex position,所以考虑 $Z_i=g_i/\sqrt n$。把样本分成两组,各有指数多个点。由 Exercise 4.39,若第一组大小至少为 $e^{C_1n}$,则以至少 $1-e^{-c_1n}$ 的概率,第一组中落在 $2B_2^n$ 内的点构成 $2B_2^n$ 的一个足够小的 net。

特别地,它包含球面 $S^{n-1}$ 的一个 $\varepsilon$-net。由 Exercise 4.38,第一组点的凸包包含 $(1-\varepsilon)B_2^n$。

另一方面,对单个 $Z_i$,事件 $\|Z_i\|_2\le1-\varepsilon$ 有概率至少 $e^{-C_2n}$。如果第二组大小至少为 $e^{C_3n}$ 且 $C_3>C_2$,则第二组至少有一个点落在 $(1-\varepsilon)B_2^n$ 内的概率至少为 $1-e^{-c_2n}$。

在两件好事件同时发生时,第二组中的这个点属于第一组点的凸包,因此它不是整个点集凸包的极点。于是这些 Gaussian 点不处于 convex position。选择 $C$ 足够大并合并失败概率,得到至少 $1-e^{-cn}$ 的结论。

Exercise 4.41 随机矩阵 operator norm 的期望
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:把 Theorem 4.4.3 和 Theorem 4.6.1 的高概率界积分成期望界。

证明思路

使用尾积分公式:若 $Z\le a+bt$ 以概率至少 $1-2e^{-t^2}$ 成立,则 $\mathbb EZ\le a+Cb$。

完整证明

(a) Theorem 4.4.3 给出对所有 $t\ge0$,

$$ \mathbb P\left\{ \|A\|>CK(\sqrt m+\sqrt n+t) \right\} \le 2e^{-t^2}. $$

由尾积分公式,

$$ \mathbb E\|A\| \le CK(\sqrt m+\sqrt n) + CK\int_0^\infty 2e^{-t^2}\,dt \le C'K(\sqrt m+\sqrt n). $$

(b) Theorem 4.6.1 证明中的矩阵形式给出

$$ \left\| \frac1mA^{\mathsf T}A-I_n \right\| \le CK^2\left( \sqrt{\frac nm}+\frac{t}{\sqrt m} + \frac nm+\frac{t^2}{m} \right) $$

以至少 $1-2e^{-t^2}$ 的概率成立。再次积分尾部,得到

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

同理,把 Theorem 4.6.1 的奇异值双侧尾界积分,得到

$$ \sqrt m-CK^2\sqrt n \le \mathbb E\,s_n(A) \le \mathbb E\,s_1(A) \le \sqrt m+CK^2\sqrt n. $$
Exercise 4.42 随机矩阵范数的下界
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明随机矩阵算子范数至少有 $\sqrt m+\sqrt n$ 的量级。

证明思路

算子范数至少不小于任意一列和任意一行的 Euclidean 范数。分别对第一列和第一行使用范数集中下界。

完整证明

对第一列 $Ae_1$,它是 $m$ 维独立、方差为一、次高斯坐标向量。范数集中下界给出

$$ \|Ae_1\|_2 \ge \sqrt m-t $$

以至少 $1-\exp(-ct^2/K^4)$ 的概率成立。同理,对第一行 $A^{\mathsf T}e_1$,

$$ \|A^{\mathsf T}e_1\|_2 \ge \sqrt n-t $$

也以同样概率成立。两者同时成立的概率至少为 $1-2\exp(-ct^2/K^4)$。

在该事件上,

$$ \|A\| \ge \max\{\|Ae_1\|_2,\|A^{\mathsf T}e_1\|_2\} \ge \max\{\sqrt m-t,\sqrt n-t\}. $$

由于 $\max(a,b)\ge(a+b)/2$,得到

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

把 $2t$ 重新记为 $t$ 并调整常数,即得题设形式。

Exercise 4.43 Subgaussian 矩阵上界:放松独立性
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:说明 Theorem 4.4.3 的 net proof 只需要固定双线性型是次高斯,不一定需要所有元素独立。

证明思路

对固定 $x,y$,$\langle Ax,y\rangle$ 是矩阵随机向量 $A$ 在 Frobenius 单位方向 $yx^{\mathsf T}$ 上的投影。只要这个投影次高斯,后续 net 和 union bound 完全相同。

完整证明

(a) 把 $A$ 看成 $\mathbb R^{mn}$ 中的随机向量。对固定 $x\in S^{n-1}$、$y\in S^{m-1}$,

$$ \langle Ax,y\rangle = \langle A,yx^{\mathsf T}\rangle_F, \qquad \|yx^{\mathsf T}\|_F=1. $$

因此

$$ \|\langle Ax,y\rangle\|_{\psi_2} \le \|A\|_{\psi_2}. $$

对两个球面的 $1/4$-net 使用 union bound,并用 Lemma 4.4.2 从 net 回到整个球面,得到 Theorem 4.4.3 的同阶结论,其中 $K=\|A\|_{\psi_2}$。

(b) 若 $A$ 的行 $A_i$ 独立、均值为零且次高斯,则

$$ \langle Ax,y\rangle = \sum_{i=1}^m y_i\langle A_i,x\rangle. $$

这些项独立、均值为零,并满足 $\|\langle A_i,x\rangle\|_{\psi_2}\le K$,其中 $K=\max_i\|A_i\|_{\psi_2}$。独立次高斯和的估计给出 $\|\langle Ax,y\rangle\|_{\psi_2}\le CK$。于是 (a) 中的 net argument 仍然适用。列独立情形对 $A^{\mathsf T}$ 使用同样论证。

Exercise 4.44 随机矩阵的一些 $p\to q$ 范数
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:估计独立次高斯随机矩阵的若干 $p\to q$ 范数,并说明 Gaussian 情形下量级最优。

证明思路

$1\to\infty$ 范数就是最大元素;$1\to2$ 范数就是最大列范数。分别用最大次高斯变量和列范数集中。

完整证明

(a) 由定义,

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

对 $mn$ 个次高斯变量使用 union bound,得到以高概率

$$ \max_{i,j}|A_{ij}| \le CK\sqrt{\log(mn)} \le CK(\sqrt{\log m}+\sqrt{\log n}). $$

(b) 对任意矩阵,

$$ \|A\|_{1\to2} = \max_{j\le n}\|Ae_j\|_2, \qquad \|A^{\mathsf T}\|_{2\to\infty} = \max_{j\le n}\|Ae_j\|_2, $$

所以二者相等。每一列是 $m$ 维独立、方差为一、次高斯坐标向量,范数集中给出

$$ \|Ae_j\|_2 \le \sqrt m+CK^2t $$

以概率至少 $1-2e^{-t^2}$ 成立。取 $t=C\sqrt{\log n}$ 并对列作 union bound,得到

$$ \|A\|_{1\to2} \le \sqrt m+CK^2\sqrt{\log n} $$

以高概率成立。

(c) 若 $A_{ij}$ 为独立标准正态,则最大元素满足

$$ \max_{i,j}|A_{ij}| \ge c\sqrt{\log(mn)} \ge c'(\sqrt{\log m}+\sqrt{\log n}) $$

以高概率成立,证明 (a) 的量级最优。另一方面,

$$ \|A\|_{1\to2} = \max_j\|Ae_j\|_2 \ge \|Ae_1\|_2 \ge c\sqrt m $$

以高概率成立;同时 $\|A\|_{1\to2}\ge\max_j|A_{1j}|\ge c\sqrt{\log n}$。合并两个下界,得到

$$ \|A\|_{1\to2} \ge c(\sqrt m+\sqrt{\log n}) $$

以高概率成立。

Exercise 4.45 Community detection
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:去掉 Theorem 4.5.2 中 $q$ 不太小的假设,仍得到 $C/(p-q)^2$ 个误分类顶点。

证明思路

原证明中 $q$ 出现在 $\lambda_1(D)$ 与 $\lambda_2(D)$ 的间隔。改用中心化或投影到 $\mathbf1^\perp$ 的谱聚类,只追踪 community 方向,此时谱间隙为 $(p-q)n/2$。

完整证明

设 $\mathbf1=(1,\ldots,1)$,令 $P$ 是到 $\mathbf1^\perp$ 的正交投影。对邻接矩阵 $A=D+R$,考虑中心化矩阵

$$ \widetilde A=PAP. $$

在 $\mathbf1^\perp$ 上,信号部分 $\widetilde D=PDP$ 只有一个非零特征方向,即 community 向量

$$ u=(1,\ldots,1,-1,\ldots,-1), $$

对应特征值

$$ \lambda=\frac{p-q}{2}n. $$

因此在 $\mathbf1^\perp$ 上的相关谱间隙是 $\lambda$,不再涉及 $q$。

噪声部分为 $\widetilde R=PRP$,且 $\|\widetilde R\|\le\|R\|$。由于 Bernoulli 噪声有界、中心化且次高斯,Corollary 4.4.7 给出

$$ \|\widetilde R\|\le C\sqrt n $$

以高概率成立。Davis-Kahan 不等式给出单位化特征向量的误差

$$ \sin\angle(\bar u,\bar v) \le \frac{C\sqrt n}{(p-q)n} = \frac{C}{(p-q)\sqrt n}, $$

其中 $\bar v$ 是 $\widetilde A$ 在 $\mathbf1^\perp$ 上 top eigenvector。乘以 $\sqrt n$ 后,存在全局符号 $\theta$ 使

$$ \|u-\theta\sqrt n\,\bar v\|_2 \le \frac{C}{p-q}. $$

每个误分类坐标至少贡献一个常数量级的平方误差,因此误分类数至多

$$ \frac{C}{(p-q)^2}. $$

算法就是:计算 $PAP$ 在 $\mathbf1^\perp$ 上的 top eigenvector,并按其坐标符号分类。

Exercise 4.46 随机矩阵双边界的另一个证明
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:用固定方向的范数集中加 net argument 重新证明次高斯矩阵奇异值双侧界。

证明思路

固定 $x$ 时,$Ax$ 是有独立次高斯坐标的向量,Theorem 3.1.1 控制 $\|Ax\|_2$。再在球面 net 上 union bound,并用 Exercise 4.37 推回所有方向。

完整证明

固定 $x\in S^{n-1}$。矩阵行 $A_i$ 独立、各向同性、次高斯,因此随机变量 $\langle A_i,x\rangle$ 独立、均值为零、方差为一且次高斯范数由 $K$ 控制。Theorem 3.1.1 给出

$$ \mathbb P\left\{ \bigl|\|Ax\|_2-\sqrt m\bigr|>CK^2t \right\} \le 2e^{-t^2}. $$

取 $S^{n-1}$ 的 $1/4$-net $\mathcal N$,满足 $|\mathcal N|\le9^n$。令 $t=C(\sqrt n+u)$,对 net 点作 union bound,可得以至少 $1-2e^{-u^2}$ 的概率,

$$ \sup_{x\in\mathcal N} \bigl|\|Ax\|_2-\sqrt m\bigr| \le CK^2(\sqrt n+u). $$

由 Exercise 4.37,把 net 上的偏差推广到整个球面:

$$ \sup_{x\in S^{n-1}} \bigl|\|Ax\|_2-\sqrt m\bigr| \le CK^2(\sqrt n+u). $$

最后注意

$$ s_1(A)=\sup_{x\in S^{n-1}}\|Ax\|_2, \qquad s_n(A)=\inf_{x\in S^{n-1}}\|Ax\|_2. $$

于是得到 Theorem 4.6.1 的双侧奇异值界。

Exercise 4.47 随机矩阵的中间 singular values
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:从全维最小奇异值界推出第 $k$ 个奇异值的下界。

证明思路

把 $A$ 限制到任意固定 $k$ 维子空间,例如前 $k$ 个坐标。限制矩阵是 $m\times k$ 次高斯矩阵;其最小奇异值给出原矩阵第 $k$ 个奇异值的下界。

完整证明

令 $E\subset\mathbb R^n$ 为前 $k$ 个坐标张成的子空间,令 $A_E$ 表示 $A$ 在 $E$ 上的限制,即 $A$ 的前 $k$ 列。它是 $m\times k$ 随机矩阵,行仍独立、各向同性且次高斯范数由 $K$ 控制。

由 Theorem 4.6.1 应用于 $A_E$,以至少 $1-2e^{-t^2}$ 的概率,

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

由 Courant-Fischer 变分刻画,

$$ s_k(A) = \max_{\dim F=k} \min_{x\in F\cap S^{n-1}}\|Ax\|_2 \ge \min_{x\in E\cap S^{n-1}}\|Ax\|_2 = s_k(A_E). $$

合并即得

$$ s_k(A) \ge \sqrt m-CK^2(\sqrt k+t). $$
Exercise 4.48 相对误差下的 covariance estimation
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:把 covariance estimation 改写成白化后的相对误差估计。

证明思路

令 $Y=\Sigma^{-1/2}X$。则 $Y$ 各向同性,样本协方差为 $\Sigma^{-1/2}\Sigma_m\Sigma^{-1/2}$。相对误差正是这个白化样本协方差与 $I$ 的算子范数距离。

完整证明

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

则 $\mathbb EYY^{\mathsf T}=I_n$,并且

$$ \frac1m\sum_{i=1}^mY_iY_i^{\mathsf T} = \Sigma^{-1/2}\Sigma_m\Sigma^{-1/2} =:\widetilde\Sigma_m. $$

条件 (4.29) 在白化后给出 $\|Y\|_{\psi_2}\le K$。由 Theorem 4.7.1 应用于 $Y$,

$$ \mathbb E\|\widetilde\Sigma_m-I\| \le CK^2\left(\sqrt{\frac nm}+\frac nm\right). $$

另一方面,对任意 $v\ne0$,令 $w=\Sigma^{1/2}v$,则

$$ \frac{v^{\mathsf T}\Sigma_m v}{v^{\mathsf T}\Sigma v} = \frac{w^{\mathsf T}\widetilde\Sigma_m w}{\|w\|_2^2}. $$

因此

$$ \sup_{v\in\mathbb R^n} \left| \frac{v^{\mathsf T}\Sigma_m v}{v^{\mathsf T}\Sigma v}-1 \right| = \|\widetilde\Sigma_m-I\|. $$

取期望即得所需相对误差界。

Exercise 4.49 High-probability covariance estimation
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:验证 covariance estimation 的高概率版本。

证明思路

Theorem 4.7.1 的证明本质上对白化样本矩阵应用 Theorem 4.6.1。保留其中的尾参数并令 $t=\sqrt u$ 即得高概率形式。

完整证明

先假设 $\Sigma$ 可逆,令 $Y=\Sigma^{-1/2}X$。如 Exercise 4.48,$Y$ 各向同性且次高斯参数由 $K$ 控制。Theorem 4.6.1 的矩阵形式给出:对任意 $u\ge0$,以至少 $1-2e^{-u}$ 的概率,

$$ \left\| \frac1m\sum_{i=1}^mY_iY_i^{\mathsf T}-I \right\| \le CK^2\left( \sqrt{\frac{n+u}{m}} + \frac{n+u}{m} \right). $$

乘回 $\Sigma^{1/2}$,得到 Loewner 意义下

$$ -\delta\Sigma \preceq \Sigma_m-\Sigma \preceq \delta\Sigma, $$

其中 $\delta$ 为右侧。于是

$$ \|\Sigma_m-\Sigma\| \le \delta\|\Sigma\|. $$

若 $\Sigma$ 不可逆,可在 $\Sigma+\eta I$ 上做上述论证后令 $\eta\downarrow0$,或限制到 $\Sigma$ 的值域。于是得到 Remark 4.7.3 的高概率界。

Exercise 4.50 低秩矩阵的 covering numbers
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明秩 $r$、Frobenius 范数至多 $1$ 的矩阵类覆盖数按约 $(m+n)r$ 维增长。

证明思路

上界直接离散化 SVD 的三个部分:左奇异向量、右奇异向量、奇异值。下界不必使用流形理论;只要把一个全维矩阵球嵌入低秩矩阵类即可。

完整证明

(a) 注意 $\|U-V\|_{1\to2}$ 就是列差的最大 Euclidean 范数:

$$ \|U-V\|_{1\to2} = \max_{j\le r}\|U_{:j}-V_{:j}\|_2. $$

取单位球面 $S^{m-1}$ 的一个 $(\varepsilon/2)$-net $\mathcal N$,满足 $|\mathcal N|\le(C/\varepsilon)^m$。对 $U\in O_{m,r}$ 的每一列 $U_{:j}$,选择 $u_j\in\mathcal N$ 使 $\|U_{:j}-u_j\|_2\le\varepsilon/2$。所有 $r$ 元组 $(u_1,\ldots,u_r)$ 构成 $O_{m,r}$ 的一个外部 $(\varepsilon/2)$-net;这些近似矩阵本身不必仍在 $O_{m,r}$ 中。由 Exercise 4.25,可把外部中心移回 $O_{m,r}$,并得到内部 $\varepsilon$-net,数量不增加。于是

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

(b) 任意 $A\in M_{m,n,r}$ 可写成 SVD

$$ A=U\operatorname{diag}(s)V^{\mathsf T}, \qquad U\in O_{m,r},\quad V\in O_{n,r},\quad \|s\|_2\le1. $$

分别取 $O_{m,r}$、$O_{n,r}$ 和 Euclidean 单位球 $B_2^r$ 的 $\varepsilon/3$-nets。若 $U,V,s$ 分别由 $U_0,V_0,s_0$ 逼近,则

$$ \|U\operatorname{diag}(s)V^{\mathsf T} - U_0\operatorname{diag}(s_0)V_0^{\mathsf T}\|_F \le \|U-U_0\|_{1\to2}\|s\|_2 + \|s-s_0\|_2 + \|V-V_0\|_{1\to2}\|s_0\|_2 \le \varepsilon. $$

三个 net 的大小相乘,得到

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

(c) 先假设 $m\ge n$。考虑所有只在前 $r$ 列可能非零的矩阵

$$ A=[B\ 0], \qquad B\in\mathbb R^{m\times r}, \qquad \|B\|_F\le1. $$

其中满列秩的 $B$ 给出 rank 正好为 $r$ 的矩阵,且满列秩矩阵在 $\mathbb R^{mr}$ 的单位 Frobenius 球中占满测度。若 $M_{m,n,r}$ 可由 $N$ 个 Frobenius 半径 $\varepsilon$ 的球覆盖,则这些球也覆盖上述满秩子集。忽略测度为零的秩亏集合并比较 $\mathbb R^{mr}$ 中体积,得到

$$ N\cdot \operatorname{Vol}(\varepsilon B_2^{mr}) \ge \operatorname{Vol}(B_2^{mr}), $$

所以 $N\ge(c/\varepsilon)^{mr}$(当 $\varepsilon$ 大于固定常数时右侧至多为 $1$,结论平凡)。由于 $m\ge n$,有 $mr\ge (m+n)r/2$,故

$$ \mathcal N(M_{m,n,r},\|\cdot\|_F,\varepsilon) \ge \left(\frac c\varepsilon\right)^{(m+n)r/2}. $$

若 $n>m$,则改为只让前 $r$ 行可能非零,得到 $N\ge(c/\varepsilon)^{nr}$,同样推出所需下界。这说明上界中的指数阶在常数因子意义下正确。

Exercise 4.51 Gaussian mixture model 的 spectral clustering
状态:完整证明已按当前笔记标准给出闭合推导。
证明目标:证明样本协方差的 top eigenvector 可以学习 Gaussian mixture model 的分离方向,从而实现谱聚类。

证明思路

总体协方差是 $I+\mu\mu^{\mathsf T}$,top eigenvector 是 $\mu$ 方向。协方差估计给出样本 top eigenvector 接近该方向;投影到近似方向后,Gaussian 噪声不能改变大多数点的符号。

完整证明

(a) 写 $X=\theta\mu+g$,其中 $\theta$ 是 Rademacher 变量,$g\sim N(0,I_n)$。由于 $\theta$ 与 $g$ 独立且均值为零,

$$ \Sigma=\mathbb EXX^{\mathsf T} = \mu\mu^{\mathsf T}+I_n. $$

因此 top eigenvector 为

$$ u=\frac{\mu}{\|\mu\|_2}, $$

对应特征值 $1+\|\mu\|_2^2$;其余方向特征值为 $1$,谱间隙为 $\|\mu\|_2^2$。

(b) 若使用真实方向 $u$ 分类,则

$$ \langle X,u\rangle = \theta\|\mu\|_2+\langle g,u\rangle. $$

令 $M=\|\mu\|_2$。若真实标签为 $\theta_i$,则

$$ \theta_i\langle X_i,u\rangle = M+\theta_i\langle g_i,u\rangle. $$

因此真实方向误分类概率为 $\mathbb P\{N(0,1)\le-M\}$。更强地,真实方向 margin 小于 $M/2$ 的概率为 $\mathbb P\{N(0,1)\le-M/2\}$。当 $M\ge C$ 且 $C$ 足够大时,这个概率小于 $0.001$。对 $m$ 个样本用 Chernoff 或 Hoeffding 不等式,可得至少以概率 $0.997$,满足 $\theta_i\langle X_i,u\rangle<M/2$ 的样本比例至多 $0.3\%$。

(c) 该 mixture 分布满足 Theorem 4.7.1 的相对次高斯条件,且参数为绝对常数:对任意 $x$,

$$ \|\langle X,x\rangle\|_{\psi_2} \le C\left(|\langle\mu,x\rangle|+\|x\|_2\right) \le C\|\langle X,x\rangle\|_{L^2}. $$

若 $m\ge Cn$ 且常数 $C$ 足够大,covariance estimation 的高概率版本给出

$$ \|\Sigma_m-\Sigma\| \le c\|\mu\|_2^2 $$

以概率至少 $0.995$ 成立,其中 $c>0$ 可任意小。

(d) 由 Davis-Kahan 不等式和谱间隙 $\|\mu\|_2^2$,样本 top eigenvector $v$ 满足,在差一个全局符号后,

$$ \|v-u\|_2\le \eta, $$

其中 $\eta>0$ 可通过增大样本量常数取得足够小。同时在同一好事件上,

$$ \|\Sigma_m\| \le \|\Sigma\|+\|\Sigma_m-\Sigma\| \le C(1+M^2). $$

(e) 现在处理关键的依赖性问题:$v$ 是由样本算出的,不能把 $\langle X_i,v-u\rangle$ 当作固定方向上的独立噪声。我们改用确定性计数。记 $h=v-u$。若某个样本满足

$$ \theta_i\langle X_i,u\rangle\ge M/2 \quad\text{且}\quad |\langle X_i,h\rangle|<M/2, $$

则它按 $v$ 的符号也会被正确分类,因为

$$ \theta_i\langle X_i,v\rangle = \theta_i\langle X_i,u\rangle + \theta_i\langle X_i,h\rangle > 0. $$

因此额外错误只能来自 $|\langle X_i,h\rangle|\ge M/2$ 的样本。由 Markov 型计数,

$$ \#\{i:|\langle X_i,h\rangle|\ge M/2\} \le \frac{4}{M^2} \sum_{i=1}^m\langle X_i,h\rangle^2 = \frac{4m}{M^2}h^{\mathsf T}\Sigma_m h. $$

在 (d) 的好事件上,

$$ h^{\mathsf T}\Sigma_m h \le \|\Sigma_m\|\|h\|_2^2 \le C(1+M^2)\eta^2. $$

由于 $M\ge C\ge1$,额外错误比例至多 $C'\eta^2$。选择样本量常数使 $\eta$ 足够小,可令它小于 $0.3\%$。再加上 (b) 中真实方向 margin 不足的样本,误分类比例至多 $1\%$。合并两个好事件,概率至少 $0.99$。这证明了 Theorem 4.7.5。

学习检查表

  • [ ] 能从 $A^{\mathsf T}A$ 推出 SVD。
  • [ ] 能解释算子范数、最大奇异值和双线性形式之间的等价关系。
  • [ ] 能写出 net argument 的四步证明模板。
  • [ ] 能解释为什么随机矩阵范数量级是 $\sqrt m+\sqrt n$。
  • [ ] 能把社区检测拆成“信号矩阵、噪声矩阵、扰动定理”。
  • [ ] 能说明协方差估计为什么需要控制样本矩阵的奇异值。
  • [ ] 能从译文、学习笔记或并排页跳转到对应证明位置。

后续衔接

第 4 章当前已经完成结构化学习页、关键定理证明区、51 个习题完整证明和 proof link 锚点。后续若继续精修,应优先复核应用题中的常数条件,以及低秩矩阵 covering number 下界的体积论证细节。