精校翻译 Ch.9 矩阵偏差
第 9 章精校翻译:集合上的随机矩阵偏差

翻译状态说明

本页为第 9 章的结构化精校翻译,覆盖正文、Notes、Exercises 9.1-9.43、原书 Figure 9.1-9.10,并为正式证明、正文隐藏验证和习题加入学习笔记证明跳转。

第 9 章 Deviations of Random Matrices on Sets

固定一个向量 $x\in\mathbb R^n$ 时,若 $A$ 是 $m\times n$ 随机矩阵,行独立、isotropic 且 subgaussian,则 norm concentration 告诉我们

$$ \|Ax\|_2\approx \sqrt m\|x\|_2. \tag{9.1} $$

本章问一个更强的问题:这个近似能否同时对集合 $T\subset\mathbb R^n$ 中所有 $x$ 成立?答案由 $T$ 的 Gaussian complexity 控制。第 9.1 节给出 matrix deviation inequality;后续章节把它依次应用到 random projections、covariance estimation、Johnson-Lindenstrauss、random sections、linear inverse problems、sparse recovery、general norm deviations、two-sided Chevet 和 Dvoretzky-Milman theorem。

9.1 Matrix deviation inequality

Theorem 9.1.1Matrix deviation inequality

令 $A$ 为 $m\times n$ 随机矩阵,行 $A_i$ 独立、isotropic 且 subgaussian。则对任意 $T\subset\mathbb R^n$,

$$\mathbb E\sup_{x\in T}\left|\|Ax\|_2-\sqrt m\|x\|_2\right|\le CK^2\gamma(T),$$

其中 $\gamma(T)$ 是 Gaussian complexity,$K=\max_i\|A_i\|_{\psi_2}$。

查看学习笔记:Theorem 9.1.1 完整证明

证明计划是把

$$ Z_x=\|Ax\|_2-\sqrt m\|x\|_2 \tag{9.2} $$

看成由 $x$ 索引的随机过程,并用第 8 章 Talagrand comparison 的几何形式。关键是验证该过程具有 subgaussian increments。

Theorem 9.1.2Subgaussian increments

在 Theorem 9.1.1 的假设下,过程 (9.2) 满足

$$\|Z_x-Z_y\|_{\psi_2}\le CK^2\|x-y\|_2,\qquad x,y\in\mathbb R^n.$$ 查看学习笔记:Theorem 9.1.2 完整证明

Theorem 9.1.2 的证明分四步:先处理 $x$ 为单位向量、$y=0$;再对单位向量比较 squared norms;接着回到原始 norm;最后用齐次性和投影到球面处理一般 $x,y$。

Reverse triangle geometry
Figure 9.1:当 $\bar y=y/\|y\|_2$ 时,可近似反向使用三角不等式:$\|x-\bar y\|_2+\|\bar y-y\|_2\le\sqrt2\|x-y\|_2$。

查看学习笔记:Figure 9.1 反向三角不等式验证

Remark 9.1.3-9.1.5三种偏差版本

Theorem 9.1.1 可通过 centering 得到围绕 $\mathbb E\|Ax\|_2$ 的版本;通过高概率 Talagrand comparison 得到 tail bound;也可推出 quadratic process $\|Ax\|_2^2-m\|x\|_2^2$ 的偏差界。

查看 Exercise 9.2 查看 Exercise 9.3

9.2 Random matrices, covariance estimation, and Johnson-Lindenstrauss

9.2.1 Singular values of random matrices

取 $T=S^{n-1}$,matrix deviation inequality 给出

$$ \sup_{x\in S^{n-1}}\left|\|Ax\|_2-\sqrt m\right| \lesssim K^2\sqrt n, $$

这等价于随机矩阵奇异值集中在 $\sqrt m$ 附近。

9.2.2 Random projections of sets

Proposition 9.2.1Sizes of random projections of sets

令 $P=A/\sqrt n$ 为 subgaussian projection。对有界集合 $T\subset\mathbb R^n$,投影后直径由 $\sqrt{m/n}\operatorname{diam}(T)$ 与 spherical width $w_s(T)$ 共同控制。

查看学习笔记:Proposition 9.2.1 完整证明

该命题把第 7.6 节的 random projection phase transition 推广到 subgaussian matrices。

9.2.3 Covariance estimation for low-dimensional distributions

Theorem 9.2.2Covariance estimation for low-dimensional distributions

若 $X$ 是 subgaussian random vector,$\Sigma=\mathbb EXX^{\mathsf T}$,$\Sigma_m=m^{-1}\sum_{i=1}^mX_iX_i^{\mathsf T}$,且 $r=\operatorname{tr}(\Sigma)/\|\Sigma\|$,则

$$\mathbb E\|\Sigma_m-\Sigma\|\le CK^4\left(\sqrt{\frac rm}+\frac rm\right)\|\Sigma\|.$$ 查看学习笔记:Theorem 9.2.2 完整证明

证明把分布写成 $X=\Sigma^{1/2}Z$,其中 $Z$ isotropic。然后把 sample covariance operator norm 写成 ellipsoid $T=\Sigma^{1/2}S^{n-1}$ 上的 quadratic matrix deviation。

9.2.4 Johnson-Lindenstrauss lemma for infinite sets

Lemma 9.2.4Additive Johnson-Lindenstrauss lemma

若 $\mathcal X\subset\mathbb R^n$ 有界,$Q=A/\sqrt m$,则以高概率

$$\left|\|Qx-Qy\|_2-\|x-y\|_2\right|\le \delta,\qquad x,y\in\mathcal X,$$

其中 $\delta=CK^2w(\mathcal X)/\sqrt m$。

查看学习笔记:Lemma 9.2.4 完整证明

有限集的经典 JL lemma 是该结果的相对误差版本;无限集只能期待 additive error,尺度由 Gaussian width 控制。

9.3 Random sections: the $M^*$ bound and escape theorem

9.3.1 The $M^*$ bound

Theorem 9.3.1$M^*$ bound

令 $E=\ker A$,其中 $A$ 的行独立、isotropic 且 subgaussian。则任意有界 $T\subset\mathbb R^n$ 满足

$$\mathbb E\operatorname{diam}(T\cap E)\le \frac{CK^2w(T)}{\sqrt m}.$$ 查看学习笔记:Theorem 9.3.1 完整证明

证明直接对 $T-T$ 应用 matrix deviation inequality;当 $x,y\in T\cap\ker A$ 时,$\|Ax-Ay\|_2=0$,于是可反推出 $\|x-y\|_2$ 的上界。

Random section of cross-polytope Random section of convex set
Figure 9.2:随机子空间切割 convex set;对 cross-polytope,随机切片通常错开尖刺并穿过主体。

9.3.2 The escape theorem

Theorem 9.3.4Escape theorem

若 $T\subset S^{n-1}$ 且

$$m\ge CK^4w(T)^2,$$

则 $E=\ker A$ 以概率至少 $1-2e^{-cm/K^4}$ 满足 $T\cap E=\varnothing$。

查看学习笔记:Theorem 9.3.4 完整证明
Escape theorem
Figure 9.3:escape theorem 描述随机子空间何时避开球面子集。

9.4 Application: high-dimensional linear models

线性观测模型写成

$$ y=Ax+w, $$

其中 $A$ 已知,$w$ 是噪声,目标是恢复高维向量 $x$。

High-dimensional linear model
Figure 9.4:高维线性模型 $y=Ax+w$。
Audio sampling recovery
Figure 9.5:audio sampling 中从少量随机采样恢复信号。

9.4.1 Constrained recovery

Theorem 9.4.4Constrained recovery

若 $x\in T$ 且 $y=Ax$,用约束程序在 $T$ 中寻找满足 $Ax'=y$ 的解 $\hat x$,则

$$\mathbb E\|\hat x-x\|_2\lesssim \frac{K^2w(T)}{\sqrt m}.$$ 查看学习笔记:Theorem 9.4.4 完整证明
Constrained high-dimensional linear problem
Figure 9.6:约束恢复问题的几何图像。

9.4.2 Example: sparse recovery

Corollary 9.4.8Sparse recovery

若未知向量 $x$ 是 $s$-sparse,$\|x\|_2\le1$,则 $\ell^1$ 约束恢复满足

$$\mathbb E\|\hat x-x\|_2\lesssim K^2\sqrt{\frac{s\log n}{m}}.$$ 查看学习笔记:Corollary 9.4.8 完整证明

更紧的 convexification 可把 $\log n$ 改进为 $\log(en/s)$;对应练习见 Exercises 9.25-9.27。

9.4.3 Example: low-rank recovery

Corollary 9.4.11Low-rank matrix recovery

若 $d\times d$ 矩阵 $X$ 的 rank 至多 $r$ 且 $\|X\|_F\le1$,Gaussian linear measurements 下 nuclear-norm constrained recovery 满足

$$\mathbb E\|\hat X-X\|_F\lesssim \sqrt{\frac{rd}{m}}.$$ 查看学习笔记:Corollary 9.4.11 完整证明

9.5 Application: exact sparse recovery

9.5.1 Exact recovery based on the escape theorem

Theorem 9.5.1Exact sparse recovery

若 $m\gtrsim K^4s\log(en/s)$,则随机矩阵 $A$ 以高概率能通过 basis pursuit 精确恢复所有 $s$-sparse 向量。

查看学习笔记:Theorem 9.5.1 完整证明

几何上,exact recovery 等价于随机 null space 与 $\ell^1$ ball 在 $x$ 处的 tangent cone 只在原点相交。

Exact sparse recovery
Figure 9.7:exact sparse recovery 的 tangent cone 几何。
Lemma 9.5.2-9.5.3Support 与 approximate sparsity

若 $\hat x$ 是 $\ell^1$ 最小化解,误差 $h=\hat x-x$ 在 $x$ 的 support 上不会太小;进一步,$h$ 可被控制在 approximate sparse set 中。

Lemma 9.5.2 Lemma 9.5.3

9.5.2 Restricted isometries

Definition 9.5.5RIP

矩阵 $A$ 对 $s$-sparse vectors 满足 restricted isometry property,若它在所有 $s$-sparse 向量上近似保持 Euclidean norm。

Theorem 9.5.6RIP implies exact recovery

若 $A$ 满足足够阶数和常数的 RIP,则每个 $s$-sparse vector 都可由 $\ell^1$ minimization 精确恢复。

查看学习笔记:Theorem 9.5.6 完整证明
Theorem 9.5.7Random matrices satisfy RIP

若 $m\gtrsim K^4s\log(en/s)$,则 isotropic subgaussian random matrices 以高概率满足 $s$-sparse 集合上的 RIP。

查看学习笔记:Theorem 9.5.7 完整证明

9.6 Deviations of random matrices for general norms

Definition 9.6.1Positive homogeneous and subadditive functions

这里允许 $f$ 不是非负 norm,但要求 positive homogeneous 与 subadditive,并满足 $f(x)\le b\|x\|_2$。

Theorem 9.6.3General matrix deviation inequality

若 $A$ 是 i.i.d. Gaussian matrix,$f$ positive homogeneous、subadditive 且 $f(x)\le b\|x\|_2$,则

$$\mathbb E\sup_{x\in T}|f(Ax)-\mathbb Ef(Ax)|\le Cb\gamma(T).$$ 查看学习笔记:Theorem 9.6.3 完整证明
Theorem 9.6.4Subgaussian increments for general norms

过程 $Z_x=f(Ax)-\mathbb Ef(Ax)$ 满足

$$\|Z_x-Z_y\|_{\psi_2}\le Cb\|x-y\|_2.$$ 查看学习笔记:Theorem 9.6.4 完整证明
Orthogonal u and v from x and y
Figure 9.8:由 $x,y$ 构造正交向量 $u=(x+y)/2$ 和 $v=(x-y)/2$。

9.7 Two-sided Chevet inequality and Dvoretzky-Milman theorem

9.7.1 Two-sided Chevet inequality

Theorem 9.7.1Two-sided Chevet inequality

若 $A$ 为 i.i.d. Gaussian matrix,$T\subset\mathbb R^n$、$S\subset\mathbb R^m$ 有界,则

$$\mathbb E\sup_{x\in T}\left|\sup_{y\in S}\langle Ax,y\rangle-w(S)\|x\|_2\right| \le C\gamma(T)\operatorname{rad}(S).$$ 查看学习笔记:Theorem 9.7.1 完整证明

9.7.2 Dvoretzky-Milman Theorem

Theorem 9.7.2Dvoretzky-Milman theorem

若 $A$ 为 $m\times n$ i.i.d. Gaussian matrix,$T\subset\mathbb R^n$ 有界,则以概率至少 $0.99$,

$$r_-B_2^m\subset \operatorname{conv}(AT)\subset r_+B_2^m,$$

其中 $r_\pm=w(T)\pm C\sqrt m\,\operatorname{rad}(T)$。

查看学习笔记:Theorem 9.7.2 完整证明
Dvoretzky-Milman random projection
Figure 9.9:8 维 cube 与 $10^4$ 个 Gaussian points 投影到平面的近圆形效果。
Remark 9.7.3-9.7.5Effective dimension 与 projection phase transition

若 $m\lesssim d(T)\asymp w(T)^2/\operatorname{rad}(T)^2$,随机投影后的凸包接近 Euclidean ball。与第 7.6 和第 9.2.2 节结合,可看到随机投影在 effective dimension 附近发生相变。

9.8 Notes

Matrix deviation inequality 是本章的主引擎,它把第 8 章的 Talagrand comparison 与随机矩阵问题结合起来。$M^*$ bound、escape theorem、compressed sensing、RIP、general norm deviations、two-sided Chevet 和 Dvoretzky-Milman theorem 都可看作同一原则的不同投影。

Exercises

Exercises 9.1-9.6Matrix deviation 的基础变体

这些题补反向三角不等式、centered deviation、quadratic deviation、anisotropic deviation、quadratic empirical process 与 random projection 版本。

9.1 9.2 9.3 9.4 9.5 9.6
Exercises 9.7-9.16随机投影、$M^*$ 与 escape theorem

这些题补 projection tail、Grassmannian projection、covariance high-probability、JL 推导、additive JL 必要性、affine/high-probability $M^*$、$\ell^p$ ball slice、escape tightness 与 sticker 版本。

9.7 9.8 9.9 9.10 9.11 9.12 9.13 9.14 9.15 9.16
Exercises 9.17-9.29高维线性模型、稀疏恢复与低秩恢复

这些题补 constrained recovery 的 MSE、norm minimization、noisy recovery、Lasso、well-posed sparse recovery、非凸 $\ell^p$、approximate sparse recovery、convexification、Gaussian width、Garnaev-Gluskin 与低秩恢复扩展。

9.17 9.18 9.19 9.20 9.21 9.22 9.23 9.24 9.25 9.26 9.27 9.28 9.29
Exercises 9.30-9.43Exact recovery、RIP、general norms 与 Dvoretzky

这些题补 exact recovery 几何、noisy exact recovery、nullspace property、random projection RIP、subadditivity、anisotropic/general high-probability deviation、general norm JL、$\ell^1/\ell^\infty$ embeddings、duality、high-probability Dvoretzky、Gaussian cloud 与 true projection version。

9.30 9.31 9.32 9.33 9.34 9.35 9.36 9.37 9.38 9.39 9.40 9.41 9.42 9.43
lp balls for p<1
Figure 9.10:$\ell^p$ unit balls 在 $0<p<1$ 时非凸,因此 $\|\cdot\|_p$ 不是 norm。
学习笔记 Ch.9 矩阵偏差
第 9 章学习笔记:集合上的随机矩阵偏差

一句话定位

第 9 章把第 8 章的 chaining / Talagrand comparison 变成随机矩阵的统一偏差工具:先证明 matrix deviation inequality,也就是 $\|Ax\|_2$ 在整个集合 $T$ 上同时接近 $\sqrt m\|x\|_2$,再把它用于随机投影、协方差估计、线性反问题、稀疏恢复、RIP、一般范数偏差和 Dvoretzky-Milman theorem

本章导读

本章不是一串孤立应用,而是一条“一个主定理,多种几何后果”的线。Matrix deviation inequality 负责把随机矩阵在集合上的最大偏差控制在 Gaussian complexity 量级;当集合换成 $T-T$、sphere subset、$\ell^1$ ball、low-rank model、support function 时,就分别得到 $M^*$ boundescape theorem、compressed sensing、low-rank recovery、two-sided Chevet 与 Dvoretzky-Milman。

章节 内容 在主线中的作用
9.1 Matrix deviation inequality 本章主引擎
9.2 Random projections、covariance、JL 解释偏差定理如何变成嵌入与估计
9.3 $M^*$ bound 与 escape theorem 用 null space 几何控制随机截面
9.4 High-dimensional linear models 把几何偏差转成恢复误差
9.5 Exact sparse recovery 与 RIP 从 approximate recovery 走到 exact recovery
9.6 General norm deviations 从 Euclidean norm 推广到一般 subadditive homogeneous function
9.7 Two-sided Chevet、Dvoretzky 用 support function 描述随机投影后的近圆性

本页使用方式

你现在卡在哪里 先看哪里 读完应形成的判断
不知道本章主线 9.1 和本章主线表 所有应用都来自 matrix deviation 的不同索引集合。
Theorem 9.1.2 证明太长 关键定理证明 核心是 squared norm 差值 + Bernstein + 齐次化。
$M^*$ 和 escape theorem 的关系 9.3 一个控制随机截面直径,一个控制随机子空间避开集合。
稀疏恢复为什么与 tangent cone 有关 9.5 exact recovery 等价于 null space 不碰 tangent cone 的球面部分。
Dvoretzky 为什么是 Chevet 的后果 9.7 support function 接近 Euclidean norm 等价于 convex hull 接近 Euclidean ball。

本章主线

推进层 要解决的问题 关键转折 后续用途
集合上矩阵偏差 $\|Ax\|_2$ 是否同时近似 $\sqrt m\|x\|_2$? $Z_x$ 具有 subgaussian increments 所有后续应用
投影与估计 随机矩阵如何保距离、估协方差? 对差集、ellipsoid 用 deviation JL、covariance estimation
随机截面 $\ker A$ 与集合怎样相交? 对 $T-T$ 或球面子集用 deviation $M^*$、escape theorem
线性反问题 约束恢复误差多大? 误差落在 $T-T\cap\ker A$ sparse/low-rank recovery
Exact recovery 何时唯一精确恢复? null space 避开 tangent cone basis pursuit、RIP
一般范数 Euclidean norm 以外怎么办? support function + Gaussian concentration two-sided Chevet、Dvoretzky

分层阅读路线

先抓住一个主引擎
第 9 章几乎所有结论都在问同一件事:随机矩阵在某个集合上偏离多少?

读法不是把应用逐个背下来,而是每遇到一个应用就问:这里的索引集合是什么?它的 Gaussian width 是多少?偏差界如何转成几何或恢复结论?

第一遍先抓三条线
  1. Matrix deviation 是统一工具。
  2. $M^*$ 与 escape theorem 是 null space 几何。
  3. 恢复问题本质是误差集合与 $\ker A$ 的相交问题。
Matrix deviationRandom sectionsRecovery errorExact recoveryGeneral normsDvoretzky
层次 先掌握什么 关键入口 暂时怎么处理
第一遍:主线阅读 Matrix deviation、随机截面、线性反问题、exact sparse recovery 的逻辑链 Theorem 9.1.1、Theorem 9.3.1、Theorem 9.4.4、Theorem 9.5.1 先看懂索引集合是什么、Gaussian width 是多少、误差集合怎样进入 $\ker A$。
第二遍:证明精读 Matrix deviation 的增量证明、$M^*$、escape、sharp sparse / low-rank width、general norm deviation Theorem 9.1.2、Theorem 9.3.4、Exercises 9.25-9.29、Theorem 9.6.3 常数和 tail 细节放在这一遍集中处理。
第三遍:习题与应用 把偏差定理迁移到 covariance、JL、recovery、RIP、Dvoretzky Exercises 9.1-9.43,尤其 9.17-9.31、9.37-9.43 按基础验证、核心证明、高价值挑战分层做题。
专题回看 compressed sensing、low-rank recovery、conic geometry、Dvoretzky-Milman 第 10 部分深入阅读路线与 References 适合在第 7、8 章 Gaussian width、Chevet 和 chaining 读熟后再进入。
练习层级 建议题目 训练目的
基础验证 9.1-9.2、9.6-9.8、9.11、9.16、9.22-9.23 检查投影、null space、norm 和稀疏性定义能否直接使用。
核心证明 9.3-9.5、9.9-9.15、9.17-9.21、9.24-9.27 把 matrix deviation、$M^*$、escape 和 recovery 主线串起来。
高价值挑战 9.28、9.31、9.37-9.43 适合第二遍证明精读或专题回看:分别对应 Garnaev-Gluskin、带噪恢复、一般范数 JL 与 Dvoretzky。

核心对象与符号表

符号 含义 初学者读法
$\gamma(T)$ Gaussian complexity $\mathbb E\sup_{x\in T}|\langle g,x\rangle|$ 控制集合上绝对过程大小
$w(T)$ Gaussian width $\mathbb E\sup_{x\in T}\langle g,x\rangle$ 控制非对称方向平均宽度
$\operatorname{rad}(T)$ $\sup_{x\in T}\|x\|_2$ 半径
$d(T)$ effective dimension $w(T)^2/\operatorname{diam}(T)^2$ 或对应半径版本
$E=\ker A$ 随机 null space 约束恢复中的不可见方向
tangent cone 从 $x$ 沿 feasible set 出发的误差方向 exact recovery 要避开的集合
RIP sparse vectors 上的近等距性质 compressed sensing 的 uniform guarantee
support function $h_T(y)=\sup_{x\in T}\langle x,y\rangle$ convex body 的对偶描述

关键定理卡片

定理 输入 输出 证明入口
Matrix deviation isotropic subgaussian rows $\sup_T|\|Ax\|-\sqrt m\|x\||$ 证明
Subgaussian increments $Z_x=\|Ax\|-\sqrt m\|x\|$ $\psi_2$ increment 证明
Covariance estimation effective rank $r$ $\|\Sigma_m-\Sigma\|$ 证明
$M^*$ bound $T$ 与 $\ker A$ random section diameter 证明
Escape theorem $T\subset S^{n-1}$ $\ker A$ misses $T$ 证明
Constrained recovery $x\in T$ error $\lesssim w(T)/\sqrt m$ 证明
Exact sparse recovery $m\gtrsim s\log(en/s)$ basis pursuit exact 证明
General deviation Gaussian matrix + subadditive $f$ $f(Ax)$ uniform deviation 证明
Dvoretzky-Milman Gaussian projection convex hull near ball 证明

关键定理完整证明

Theorem 9.1.1Matrix deviation inequality
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:证明 $\mathbb E\sup_{x\in T}|\|Ax\|_2-\sqrt m\|x\|_2|\le CK^2\gamma(T)$。

完整证明:定义 $Z_x=\|Ax\|_2-\sqrt m\|x\|_2$。Theorem 9.1.2 给出 $\|Z_x-Z_y\|_{\psi_2}\le CK^2\|x-y\|_2$,并且 $Z_0=0$。对过程 $Z_x$ 和 $-Z_x$ 分别应用 Talagrand comparison 的几何形式,得到

$$\mathbb E\sup_{x\in T}Z_x\le CK^2w(T),\qquad \mathbb E\sup_{x\in T}(-Z_x)\le CK^2w(-T).$$

二者合并等价于绝对值上确界,用 Gaussian complexity $\gamma(T)=\mathbb E\sup_{x\in T}|\langle g,x\rangle|$ 表示,得到结论。

Theorem 9.1.2Subgaussian increments
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:证明 $Z_x=\|Ax\|_2-\sqrt m\|x\|_2$ 的 $\psi_2$ 增量。

完整证明:第一步,若 $\|x\|_2=1$ 且 $y=0$,则 $Ax$ 的坐标独立、isotropic、subgaussian,norm concentration 给出 $\|\|Ax\|_2-\sqrt m\|_{\psi_2}\le CK^2$。第二步,若 $x,y$ 都是单位向量,展开

$$\|Ax\|_2^2-\|Ay\|_2^2=\sum_i\langle A_i,x+y\rangle\langle A_i,x-y\rangle.$$

除以 $\|x-y\|_2$ 后得到独立均值零 subexponential 和,Bernstein inequality 给出该平方差除以 $\|x-y\|_2$ 的 tail 尺度为 $\sqrt m$。第三步,用

$$|\|Ax\|_2-\|Ay\|_2|=\frac{|\|Ax\|_2^2-\|Ay\|_2^2|}{\|Ax\|_2+\|Ay\|_2}$$

并把事件分成 $\|Ax\|_2\ge\sqrt m/2$ 与补事件;补事件由第一步控制,从而得到单位球面上的 $\psi_2$ 增量。第四步,对一般 $x,y$ 通过缩放令 $\|x\|_2=1\le\|y\|_2$,设 $\bar y=y/\|y\|_2$。三角不等式、齐次性和第一步给出

$$\|Z_x-Z_y\|_{\psi_2}\le CK^2(\|x-\bar y\|_2+\|\bar y-y\|_2).$$

Figure 9.1 的几何验证给出括号不超过 $\sqrt2\|x-y\|_2$,结论成立。

Proposition 9.2.1Sizes of random projections
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:由 matrix deviation 控制 $\operatorname{diam}(PT)$。

完整证明:令 $P=A/\sqrt n$。对差集 $T-T$ 应用 Theorem 9.1.1,得到

$$\mathbb E\sup_{z\in T-T}\left|\|Az\|_2-\sqrt m\|z\|_2\right|\le CK^2\gamma(T-T).$$

除以 $\sqrt n$,并用 $\gamma(T-T)\le2\gamma(T)\asymp2\sqrt n\,w_s(T)$,得到

$$\mathbb E\operatorname{diam}(PT)\le\sqrt{m/n}\operatorname{diam}(T)+CK^2w_s(T).$$

下界或双边形式同样由偏差界应用于实现直径的差向量得到。

Theorem 9.2.2Covariance estimation
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:证明 effective rank 控制 subgaussian covariance estimation。

完整证明:写 $X=\Sigma^{1/2}Z$,其中 $Z$ isotropic subgaussian。令 $A$ 的行为 $Z_i$,则

$$\|\Sigma_m-\Sigma\|=\frac1m\sup_{x\in T}\left|\|Ax\|_2^2-m\|x\|_2^2\right|,\qquad T=\Sigma^{1/2}S^{n-1}.$$

Exercise 9.3 的 quadratic deviation 给出

$$\mathbb E\|\Sigma_m-\Sigma\|\le\frac{CK^4\gamma(T)^2+CK^2\sqrt m\,\operatorname{rad}(T)\gamma(T)}{m}.$$

对 ellipsoid,$\operatorname{rad}(T)=\|\Sigma\|^{1/2}$,$\gamma(T)\le(\operatorname{tr}\Sigma)^{1/2}$。代入并用 $r=\operatorname{tr}(\Sigma)/\|\Sigma\|$,得到 $CK^4(\sqrt{r/m}+r/m)\|\Sigma\|$。

Lemma 9.2.4Additive Johnson-Lindenstrauss
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:证明无限集合的 additive JL bound。

完整证明:对差集 $T=\mathcal X-\mathcal X$ 应用 Theorem 9.1.1 的高概率版本,得到

$$\sup_{z\in T}\left|\|Az\|_2-\sqrt m\|z\|_2\right|\le CK^2\gamma(T).$$

由 Lemma 7.5.11,$\gamma(\mathcal X-\mathcal X)\le Cw(\mathcal X)$。令 $Q=A/\sqrt m$,两边除以 $\sqrt m$,即得对所有 $x,y\in\mathcal X$ 的 additive error $\delta=CK^2w(\mathcal X)/\sqrt m$。

Theorem 9.3.1$M^*$ bound
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:控制随机截面 $T\cap\ker A$ 的直径。

完整证明:对差集 $T-T$ 应用 Theorem 9.1.1:

$$\mathbb E\sup_{x,y\in T}\left|\|A(x-y)\|_2-\sqrt m\|x-y\|_2\right|\le CK^2\gamma(T-T)\le CK^2w(T).$$

若 $x,y\in T\cap\ker A$,则 $A(x-y)=0$,所以上式内部等于 $\sqrt m\|x-y\|_2$。因此

$$\mathbb E\operatorname{diam}(T\cap\ker A)\le CK^2w(T)/\sqrt m.$$
Theorem 9.3.4Escape theorem
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:证明 $\ker A$ 以高概率避开 $T\subset S^{n-1}$。

完整证明:高概率 matrix deviation 给出

$$\sup_{x\in T}\left|\|Ax\|_2-\sqrt m\right|\le C_1K^2(w(T)+u)$$

概率至少 $1-2e^{-u^2}$。取 $u=\sqrt m/(2C_1K^2)$。若存在 $x\in T\cap\ker A$,则左边至少 $\sqrt m$,从而 $\sqrt m\le C_1K^2w(T)+\sqrt m/2$。当 $m\ge CK^4w(T)^2$ 且 $C$ 充分大时矛盾。因此该事件上 $T\cap\ker A=\varnothing$,概率为 $1-2e^{-cm/K^4}$。

Theorem 9.4.4Constrained recovery
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:证明约束恢复误差由 $w(T)/\sqrt m$ 控制。

完整证明:真实 $x$ 与恢复 $\hat x$ 都在 $T$ 中,且满足 $A\hat x=Ax$。因此误差 $h=\hat x-x$ 属于 $(T-T)\cap\ker A$。由 $M^*$ bound 应用于 $T$,

$$\mathbb E\|h\|_2\le\mathbb E\operatorname{diam}(T\cap(x+\ker A))\le CK^2w(T)/\sqrt m,$$

仿射版本可由对 $T-T$ 的同一证明得到。因此结论成立。

Corollary 9.4.8Sparse recovery
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:由 constrained recovery 推出 sparse recovery。

完整证明:若 $x$ 是 $s$-sparse 且 $\|x\|_2\le1$,则 $\|x\|_1\le\sqrt s$,所以 $x\in T=\sqrt s B_1^n\cap B_2^n$。用 $\ell^1$ constrained program 等价于在该 convex feasible set 中恢复。由 Theorem 9.4.4,误差由 $w(T)/\sqrt m$ 控制。粗略估计 $w(\sqrt sB_1^n\cap B_2^n)\le\sqrt s\,w(B_1^n)\lesssim\sqrt{s\log n}$,得到

$$\mathbb E\|\hat x-x\|_2\lesssim K^2\sqrt{s\log n/m}.$$
Corollary 9.4.11Low-rank recovery
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:证明 nuclear norm constrained recovery 的误差界。

完整证明:若 $\operatorname{rank}(X)\le r$ 且 $\|X\|_F\le1$,则 $\|X\|_*\le\sqrt r$。取 feasible set $T=\sqrt r B_*\cap B_F$。Theorem 9.4.4 的矩阵版给出误差由 $w(T)/\sqrt m$ 控制。Gaussian width 满足 $w(T)\lesssim\sqrt{rd}$,因为 nuclear norm ball 的 polar 是 operator norm ball,而 Gaussian matrix 的 operator norm 期望为 $O(\sqrt d)$。因此 $\mathbb E\|\hat X-X\|_F\lesssim\sqrt{rd/m}$。

Theorem 9.5.1Exact sparse recovery
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:用 escape theorem 证明 basis pursuit 精确恢复。

完整证明:设 $x$ 为 $s$-sparse,$\hat x$ 为 $\ell^1$ minimization 解,$h=\hat x-x$。若 $h\ne0$,则 $Ah=0$ 且 $\|x+h\|_1\le\|x\|_1$。Lemma 9.5.2 和 Lemma 9.5.3 说明归一化误差 $h/\|h\|_2$ 属于一个 approximate sparse spherical set $S$,且 $w(S)^2\lesssim s\log(en/s)$。若 $m\ge CK^4s\log(en/s)$,escape theorem 给出 $\ker A$ 与 $S$ 不相交。由于 $h/\|h\|_2\in\ker A\cap S$ 会矛盾,故 $h=0$,精确恢复成立。

Lemma 9.5.2Error is heavier on support
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:证明误差在 support 外的 $\ell^1$ 质量受 support 内控制。

完整证明:令 $S=\operatorname{supp}(x)$。由 $\ell^1$ 最小化,$\|x+h\|_1\le\|x\|_1$。分解 support 得

$$\|x_S+h_S\|_1+\|h_{S^c}\|_1\le\|x_S\|_1.$$

三角不等式给出 $\|x_S+h_S\|_1\ge\|x_S\|_1-\|h_S\|_1$,代入得到 $\|h_{S^c}\|_1\le\|h_S\|_1$。

Lemma 9.5.3Error is approximately sparse
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:由 support inequality 推出误差属于 approximate sparse set。

完整证明:由 Lemma 9.5.2,$\|h\|_1=\|h_S\|_1+\|h_{S^c}\|_1\le2\|h_S\|_1$。Cauchy-Schwarz 给出 $\|h_S\|_1\le\sqrt s\|h_S\|_2\le\sqrt s\|h\|_2$。因此 $\|h\|_1\le2\sqrt s\|h\|_2$。若归一化 $\|h\|_2=1$,则 $h\in2\sqrt sB_1^n\cap S^{n-1}$,这就是 approximate sparse spherical set。

Theorem 9.5.6RIP implies exact recovery
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:由 RIP 推出 $\ell^1$ exact recovery。

完整证明:令 $h=\hat x-x\in\ker A$,$S=\operatorname{supp}(x)$。按系数大小把 $S^c$ 分块,每块大小约为 $\lambda s$。Cone inequality 给出 $\|h_{S^c}\|_1\le\|h_S\|_1$,从而 tail blocks 的 $\ell^2$ 总量被 $\|h_S\|_2$ 控制。RIP 对 $S$ 与首个 tail block 的并集给出 $\|Ah_{S\cup T_1}\|_2$ 的下界,对剩余 tail blocks 给出上界。由于 $Ah=0$,两边比较并利用 $\lambda>(\beta/\alpha)^2$,得到 $\|h_S\|_2=0$,进而 $h=0$。

Theorem 9.5.7Random matrices satisfy RIP
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:证明 subgaussian matrices 在 sparse set 上满足 RIP。

完整证明:令 $T=S_{n,s}=\{x:\|x\|_0\le s,\|x\|_2=1\}$。Matrix deviation 的高概率版本给出 $\sup_{x\in T}|\|Ax\|_2-\sqrt m|$ 由 $K^2(w(T)+u)$ 控制。Gaussian width 满足 $w(T)\lesssim\sqrt{s\log(en/s)}$。若 $m\ge CK^4s\log(en/s)$,取 $u\asymp\sqrt m/K^2$ 的小常数倍即可使偏差小于 $\varepsilon\sqrt m$,从而得到所有 sparse vectors 上的 norm preservation,即 RIP。

Theorem 9.6.3General matrix deviation
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:证明 $f(Ax)$ 的 uniform centered deviation。

完整证明:定义 $Z_x=f(Ax)-\mathbb Ef(Ax)$。Theorem 9.6.4 证明 $\|Z_x-Z_y\|_{\psi_2}\le Cb\|x-y\|_2$。并且 $Z_0=0$。对 $Z_x$ 和 $-Z_x$ 用 Talagrand comparison 的几何形式,得到 $\mathbb E\sup_{x\in T}|Z_x|\le Cb\gamma(T)$。

Theorem 9.6.4General norm increments
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:证明 $Z_x=f(Ax)-\mathbb Ef(Ax)$ 的 $\psi_2$ 增量。

完整证明:先设 $b=1$ 且 $\|x\|_2=\|y\|_2=1$。令 $u=(x+y)/2$、$v=(x-y)/2$,则 $u\perp v$,所以 Gaussian vectors $Au$ 与 $Av$ 独立。条件在 $a=Au$ 上,$Ax=a+Av$ 与 $Ay=a-Av$。函数 $z\mapsto f(a+\|v\|_2z)$ 是 $\|v\|_2$-Lipschitz,因为 subadditivity 给出 $f(r)-f(s)\le f(r-s)\le\|r-s\|_2$。Gaussian concentration 给出 $f(a+Av)$ 与其条件期望的 $\psi_2$ 范数至多 $C\|v\|_2$;$f(a-Av)$ 同样成立且条件期望相同。相减后得到 $\|f(Ax)-f(Ay)\|_{\psi_2}\le C\|x-y\|_2$。一般 $x,y$ 按 Theorem 9.1.2 的齐次化步骤处理;恢复 $b$ 得结论。

Theorem 9.7.1Two-sided Chevet
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:由 general matrix deviation 推出 two-sided Chevet。

完整证明:取 support function $f(z)=\sup_{y\in S}\langle z,y\rangle$。它 positive homogeneous、subadditive,且 $f(z)\le\operatorname{rad}(S)\|z\|_2$。Theorem 9.6.3 给出

$$\mathbb E\sup_{x\in T}|f(Ax)-\mathbb Ef(Ax)|\le C\gamma(T)\operatorname{rad}(S).$$

由于 $Ax\sim\|x\|_2g$,$\mathbb Ef(Ax)=\|x\|_2\mathbb E\sup_{y\in S}\langle g,y\rangle=\|x\|_2w(S)$。代入即得。

Theorem 9.7.2Dvoretzky-Milman
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:证明 $\operatorname{conv}(AT)$ 夹在两个 Euclidean balls 之间。

完整证明:把 Theorem 9.7.1 应用于 $A^{\mathsf T}$,并取 $S=S^{m-1}$,得到

$$\mathbb E\sup_{y\in S^{m-1}}\left|\sup_{x\in T}\langle Ax,y\rangle-w(T)\right|\le C\sqrt m\,\operatorname{rad}(T).$$

Markov inequality 给出概率至少 $0.99$ 的同阶界。于是对所有单位 $y$,support function $h_{\operatorname{conv}(AT)}(y)$ 位于 $w(T)\pm C\sqrt m\operatorname{rad}(T)$。Support function 与 convex body inclusion 的对偶关系给出

$$r_-B_2^m\subset\operatorname{conv}(AT)\subset r_+B_2^m,$$

其中 $r_\pm=w(T)\pm C\sqrt m\operatorname{rad}(T)$。

正文隐藏验证补全

Hidden Check9.1:近似反向三角不等式
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:证明 $\|x-\bar y\|_2+\|\bar y-y\|_2\le\sqrt2\|x-y\|_2$。

完整证明:设 $\|x\|_2=1$,$r=\|y\|_2\ge1$,$\bar y=y/r$,并令 $\theta$ 为 $x$ 与 $\bar y$ 的夹角。则

$$\|x-y\|_2^2=\|x-r\bar y\|_2^2=(r-1)^2+2r(1-\cos\theta),$$

而 $\|x-\bar y\|_2^2=2(1-\cos\theta)$、$\|\bar y-y\|_2=r-1$。由 Cauchy-Schwarz,

$$(\|x-\bar y\|_2+\|\bar y-y\|_2)^2\le2(\|x-\bar y\|_2^2+\|\bar y-y\|_2^2)\le2\|x-y\|_2^2,$$

因为 $r\ge1$ 使 $2r(1-\cos\theta)\ge2(1-\cos\theta)$。

Hidden Check9.2:ellipsoid 的 radius 与 complexity
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:计算 $T=\Sigma^{1/2}S^{n-1}$ 的两个量。

完整证明:半径为 $\sup_{\|x\|=1}\|\Sigma^{1/2}x\|_2=\|\Sigma\|^{1/2}$。Gaussian complexity 满足

$$\gamma(T)=\mathbb E\sup_{\|x\|=1}|\langle g,\Sigma^{1/2}x\rangle|=\mathbb E\|\Sigma^{1/2}g\|_2\le(\mathbb E g^{\mathsf T}\Sigma g)^{1/2}=(\operatorname{tr}\Sigma)^{1/2}.$$
Hidden Check9.5:tangent cone 与 exact recovery
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:说明 exact recovery 等价于 null space 避开 tangent cone。

完整证明:若存在非零 $h\in\ker A$ 使 $x+h$ 仍在 $\ell^1$ ball 中,则 $A(x+h)=Ax$ 且 $\ell^1$ minimization 至少有另一个可行解,精确性失败。反过来,若精确性失败,则 $\hat x\ne x$ 且 $h=\hat x-x\in\ker A$,并且 $\|x+h\|_1\le\|x\|_1$,所以 $h$ 属于从 $x$ 指向 $\ell^1$ ball 的 tangent cone。归一化后就是 cone 的 spherical part 与 null space 相交。

Hidden Check9.6:$u=(x+y)/2$ 与 $v=(x-y)/2$ 正交
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:证明单位向量 $x,y$ 给出的 $u,v$ 正交。

完整证明:计算内积:

$$\langle u,v\rangle=\frac14\langle x+y,x-y\rangle=\frac14(\|x\|_2^2-\|y\|_2^2)=0.$$

因此 Gaussian matrix 作用后 $Au$ 与 $Av$ 是 independent Gaussian vectors。

Hidden Check9.7:support function 与球包含
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:说明 $h_K(y)$ 夹住 Euclidean norm 等价于 $K$ 夹在 Euclidean balls 中。

完整证明:若 $K\subset r_+B_2$,则 $h_K(y)=\sup_{x\in K}\langle x,y\rangle\le r_+\|y\|_2$。若 $r_-B_2\subset K$,则 $h_K(y)\ge\sup_{\|x\|\le r_-}\langle x,y\rangle=r_-\|y\|_2$。反向使用 separation theorem:若某点 $z\in K$ 满足 $\|z\|>r_+$,取 $y=z$ 得 $h_K(y)>r_+\|y\|$;若 $r_-B_2$ 中有点不在 $K$,可用 separating hyperplane 找到某个 $y$ 使 $h_K(y)<r_-\|y\|$。

Exercises 完整证明

Exercise 9.1Reverse triangle inequality
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:证明 Figure 9.1 的几何不等式。

完整证明:见 [隐藏验证](#proof-check-9-1-reverse-triangle)。核心是把 $\|x-y\|^2$ 写成 radial part $(\|y\|-1)^2$ 与 angular part,再用 Cauchy-Schwarz 合并两段路径长度。

Exercise 9.2Deviations from the mean
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:由 Theorem 9.1.1 推出 centered version。

完整证明:对任意 $x$,

$$|\|Ax\|_2-\mathbb E\|Ax\|_2|\le |\|Ax\|_2-\sqrt m\|x\|_2|+\mathbb E|\|Ax\|_2-\sqrt m\|x\|_2|.$$

对 $x\in T$ 取 supremum 和期望,第二项由 Jensen 和 Theorem 9.1.1 控制,得到同阶 $CK^2\gamma(T)$。

Exercise 9.3Quadratic matrix deviation
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:推出 $\|Ax\|_2^2-m\|x\|_2^2$ 的上界。

完整证明:令 $\Delta_x=\|Ax\|_2-\sqrt m\|x\|_2$。则

$$|\|Ax\|_2^2-m\|x\|_2^2|=|\Delta_x|(\|Ax\|_2+\sqrt m\|x\|_2).$$

上确界由 $\sup|\Delta_x|^2+2\sqrt m\,\operatorname{rad}(T)\sup|\Delta_x|$ 控制。对期望用 Theorem 9.1.1 及高概率版本积分得到 $\mathbb E\sup|\Delta_x|^2\le CK^4\gamma(T)^2$,并得到题设 bound。

Exercise 9.4Anisotropic matrix deviation
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:去掉 isotropic 假设。

完整证明:写 $B_i=\Sigma^{1/2}A_i$,其中 $A_i$ 在 $\Sigma$ 的 range 上 isotropic。则 $\|Bx\|_2=\|A(\Sigma^{1/2}x)\|_2$,而 deterministic target 为 $\sqrt m\|\Sigma^{1/2}x\|_2$。把 isotropic matrix deviation 应用于集合 $\Sigma^{1/2}T$,得到

$$\mathbb E\sup_{x\in T}\left|\|Bx\|_2-\sqrt m\|\Sigma^{1/2}x\|_2\right|\le CK^2\gamma(\Sigma^{1/2}T).$$
Exercise 9.5Quadratic empirical process
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:把 matrix deviation 推广到函数类的 $L^2$ norm。

完整证明:把每个 $f\in\mathcal F$ 映到向量 $(f(X_1),\dots,f(X_m))$。目标偏差为 empirical $L^2$ norm 与 population $L^2$ norm 的差。假设给出过程在 metric $d(f,g)$ 下的 subgaussian increments;按 Theorem 9.1.2 的证明,norm deviation process 具有 $CK^2d(f,g)/\sqrt m$ 尺度的 increments。Talagrand comparison 给出

$$\mathbb E\sup_f\left|\left(m^{-1}\sum_if(X_i)^2\right)^{1/2}-(\mathbb Ef(X)^2)^{1/2}\right|\le CK^2\gamma_2(\mathcal F,d)/\sqrt m.$$

线性函数类 $f_x(z)=\langle z,x\rangle$ 且 $x\in T\subset S^{n-1}$ 时回到 Theorem 9.1.1。

Exercise 9.6Deviation of random projections
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:证明 Haar 随机正交投影的 matrix deviation 版本:

$$\mathbb E\sup_{x\in T}\left|\|Px\|_2-\sqrt{\frac mn}\|x\|_2\right|\le \frac{C\gamma(T)}{\sqrt n}.$$

完整证明:先处理单位球面上的增量。令 $P=U^\mathsf T P_mU$,其中 $U$ 为 Haar 正交矩阵,$P_m$ 是前 $m$ 个坐标的正交投影。对固定单位向量 $x,y$,函数

$$F(U)=\|P_mUx\|_2-\|P_mUy\|_2$$

在正交群上是 $C\|x-y\|_2$-Lipschitz。由正交群上的 concentration inequality,

$$\|F-\mathbb EF\|_{\psi_2}\le \frac{C\|x-y\|_2}{\sqrt n}.$$

这一步也可按原书提示直接验证:由旋转不变性令 $x=e_1$、$y=\sqrt{1-\varepsilon^2}e_1+\varepsilon e_2$,展开 $\|Px\|_2^2-\|Py\|_2^2$,主随机项由投影矩阵元素 $P_{12}$ 控制,而 $P_{12}$ 在 Haar 投影下具有 $O(1/\sqrt n)$ 的 subgaussian 尺度。

$$Z_x=\|Px\|_2-\mathbb E\|Px\|_2.$$

上一步给出 $\|Z_x-Z_y\|_{\psi_2}\le C\|x-y\|_2/\sqrt n$。对过程 $Z_x$ 和 $-Z_x$ 用 Talagrand comparison 的几何形式,得到

$$\mathbb E\sup_{x\in T}|Z_x-Z_0|\le \frac{C\gamma(T)}{\sqrt n}.$$

最后处理均值偏差。对固定 $x$,$\|Px\|_2/\|x\|_2$ 与随机球面点前 $m$ 个坐标长度同分布,其平方均值为 $m/n$,且方差尺度为 $O(1/n)$;因此

$$\left|\mathbb E\|Px\|_2-\sqrt{\frac mn}\|x\|_2\right|\le \frac{C\|x\|_2}{\sqrt n}.$$

这里的 $\gamma(T)$ 是本书定义的 Gaussian complexity $\mathbb E\sup_{x\in T}|\langle g,x\rangle|$,因此 $\gamma(T)\ge c\sup_{x\in T}\|x\|_2$。所以均值项也被 $C\gamma(T)/\sqrt n$ 控制。合并随机中心化项与均值项,即得结论。

Exercise 9.7Projection size tail
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:证明 Proposition 9.2.1 的高概率版本。

完整证明:对 $T-T$ 使用 matrix deviation 高概率形式:

$$\sup_{z\in T-T}\left|\|Az\|_2-\sqrt m\|z\|_2\right|\le CK^2[w(T-T)+u\operatorname{diam}(T)].$$

令 $P=A/\sqrt n$,取 $u=c\varepsilon\sqrt m/K^2$,并用 $w(T-T)\le2\sqrt n\,w_s(T)$。整理得到题设直径上界,概率为 $1-\exp(-c\varepsilon^2m/K^4)$。

Exercise 9.8Grassmannian projection version
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:把 Proposition 9.2.1 改成真正随机子空间投影。

完整证明:用 Exercise 9.6 的 random projection deviation 代替 subgaussian matrix deviation。对差集 $T-T$ 应用该结果,得到 $\operatorname{diam}(PT)$ 与 $\sqrt{m/n}\operatorname{diam}(T)$ 的偏差由 $\gamma(T-T)/\sqrt n\asymp w_s(T)$ 控制。

Exercise 9.9Covariance estimation high probability
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:证明 Theorem 9.2.2 的 tail version。

完整证明:在 Theorem 9.2.2 的 proof 中,把 quadratic deviation 的 expectation version 换成 high-probability version。对 ellipsoid $T$,tail 参数 $u$ 把 $\gamma(T)^2$ 替换为 $\gamma(T)^2+u\operatorname{rad}(T)^2$,把 $\sqrt m\operatorname{rad}(T)\gamma(T)$ 替换为 $\sqrt m\operatorname{rad}(T)(\gamma(T)+\sqrt u\,\operatorname{rad}(T))$。代入 $\operatorname{rad}(T)^2=\|\Sigma\|$ 与 $\gamma(T)^2\le r\|\Sigma\|$,得到 $CK^4(\sqrt{(r+u)/m}+(r+u)/m)\|\Sigma\|$。

Exercise 9.10Matrix deviation to JL
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:用 matrix deviation 推出有限集 JL。

完整证明:对 normalized differences $T=\{(x-y)/\|x-y\|:x,y\in\mathcal X,x\ne y\}$ 应用 matrix deviation 高概率形式。该集合大小至多 $N^2$,所以 $\gamma(T)\le C\sqrt{\log N}$。若 $m\ge C K^4\varepsilon^{-2}\log N$,则所有差向量满足 $\left|\|A z\|_2/\sqrt m-1\right|\le\varepsilon$。这等价于 $Q=A/\sqrt m$ 是 $\varepsilon$-isometry。

Exercise 9.11Additive JL must be additive
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:说明 infinite set 不能期待 uniform relative error。

完整证明:若 $\mathcal X$ 含有一条线段或有聚点,则存在任意接近的 $x,y$。对任意非等距线性映射 $Q$,在某些极小差向量上相对误差不会因距离缩小而改善。若要求所有非零差向量的 relative error 很小,就要求 $Q$ 在 $\mathcal X-\mathcal X$ 的 span 上近等距;当该 span 维度超过 $m$ 时不可能。因此无限集自然得到 additive error。

Exercise 9.12Affine $M^*$ bound
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:证明所有 affine sections 的 $M^*$ bound。

完整证明:若 $x,y\in T\cap(z+\ker A)$,则 $x-y\in T-T$ 且 $A(x-y)=0$。因此任意 affine section 的直径都被 $\sup\{\|h\|_2:h\in T-T,Ah=0\}$ 控制。对 $T-T$ 重复 Theorem 9.3.1 证明,即得 $\mathbb E\max_z\operatorname{diam}(T\cap(z+\ker A))\le CK^2w(T)/\sqrt m$。

Exercise 9.13High-probability $M^*$
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:证明 $M^*$ bound 的 tail version。

完整证明:在 affine proof 中使用 matrix deviation 高概率形式作用于 $T-T$。得到概率至少 $1-2e^{-u^2}$ 的界

$$\operatorname{diam}(T\cap(z+\ker A))\le \frac{CK^2(w(T)+u\operatorname{diam}(T))}{\sqrt m}$$

对所有 $z$ 同时成立。

Exercise 9.14Slicing $\ell^p$ balls
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:若 $E$ 是随机 $k$ 维子空间且 $1\le k\le0.99n$,证明

$$\mathbb E\operatorname{diam}(B_p^n\cap E)\asymp_p n^{1/2-1/p},\qquad 1<p<\infty.$$

完整证明:先证上界。设 $m=n-k$ 为余维,则 $m\ge0.01n$。对 $T=B_p^n$ 使用 $M^*$ bound,

$$\mathbb E\operatorname{diam}(B_p^n\cap E)\le \frac{Cw(B_p^n)}{\sqrt m}.$$

由对偶性,$w(B_p^n)=\mathbb E\|g\|_{p'}$,其中 $1/p+1/p'=1$。Gaussian $p'$-范数的矩估计给出 $\mathbb E\|g\|_{p'}\asymp_p n^{1/p'}$。因此

$$\mathbb E\operatorname{diam}(B_p^n\cap E)\le C_p n^{1/p'-1/2}=C_p n^{1/2-1/p}.$$

再证下界。随机子空间 $E$ 中取一条 Haar 随机直线 $L=\operatorname{span}(u)$,其方向 $u$ 在 $S^{n-1}$ 上均匀。因此

$$\operatorname{diam}(B_p^n\cap E)\ge\operatorname{diam}(B_p^n\cap L)=\frac{2}{\|u\|_p}.$$

写 $u=g/\|g\|_2$,则

$$\frac{2}{\|u\|_p}=2\frac{\|g\|_2}{\|g\|_p}.$$

以正概率同时有 $\|g\|_2\ge c\sqrt n$ 且 $\|g\|_p\le C_p n^{1/p}$,所以

$$\mathbb E\operatorname{diam}(B_p^n\cap E)\ge c_p n^{1/2-1/p}.$$

上下界合并即得。几何解释也随之清楚:当 $p\le2$ 时,$B_p^n$ 的最大 Euclidean 内球半径是 $n^{1/2-1/p}$;当 $p\ge2$ 时,$B_p^n$ 的 Euclidean 外接半径是同一尺度。随机大维截面看到的正是这个 Euclidean 尺度。

Exercise 9.15Tightness of escape theorem
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:说明 escape theorem 中 $m\gtrsim w(T)^2$ 的采样阈值一般不能改进。

完整证明:固定一个 $d$ 维线性子空间 $F\subset\mathbb R^n$,令 $T=S^{n-1}\cap F$。则

$$w(T)=\mathbb E\sup_{x\in S^{n-1}\cap F}\langle g,x\rangle=\mathbb E\|P_Fg\|_2\asymp\sqrt d.$$

令 $E=\ker A$ 为随机余维 $m$ 子空间。无论 $E$ 是否随机,都有线性代数维数下界

$$\dim(E\cap F)\ge \dim E+\dim F-n=(n-m)+d-n=d-m.$$

如果 $m<d$,则 $E\cap F$ 含有非零向量,归一化后得到 $T\cap E\ne\varnothing$。也就是说,当 $m<w(T)^2$ 的常数倍时,escape 结论必然失败。取 $d$ 与 $m$ 同阶即可说明 Theorem 9.3.4 的 $w(T)^2$ 阈值在一般情形下是 sharp 的。

Exercise 9.16Sticker on the soccer ball
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:证明存在旋转使 $UT$ 避开有限点集 $\mathcal X$。

完整证明:随机取 Haar rotation $U$。对固定 $x\in\mathcal X$,$U^{-1}x$ 在球面均匀分布,所以 $\mathbb P\{x\in UT\}=\sigma_{n-1}(T)<1/N$。Union bound 给出 $\mathbb P\{UT\cap\mathcal X\ne\varnothing\}<1$。因此存在一个旋转满足 $UT\cap\mathcal X=\varnothing$。

Exercise 9.17Constrained recovery MSE
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:把 Theorem 9.4.4 的一阶误差界升级为 mean squared error。

完整证明:设 $h=\hat x-x$。和 Theorem 9.4.4 一样,$h$ 属于某个 affine section 的差集,因此 high-probability $M^*$ bound 给出:对所有 $u\ge0$,以概率至少 $1-2e^{-u^2}$,

$$\|h\|_2\le \frac{CK^2}{\sqrt m}\bigl(w(T)+u\,\operatorname{diam}(T)\bigr).$$

令 $a=CK^2w(T)/\sqrt m$,$b=CK^2\operatorname{diam}(T)/\sqrt m$。上式等价于 $\mathbb P\{\|h\|_2>a+bu\}\le2e^{-u^2}$。由 tail integration,

$$\mathbb E\|h\|_2^2\le C(a^2+b^2)\le \frac{CK^4}{m}\bigl(w(T)^2+\operatorname{diam}(T)^2\bigr).$$

若 $T$ 已按半径或直径为常数归一化,这就是对应的 MSE 形式;一般尺度由同一个公式保留。

Exercise 9.18Recovery by norm minimization
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:把 constrained recovery 写成 norm minimization。

完整证明:若 $T$ 是 norm unit ball,且真实 $x$ 满足 $\|x\|_T\le1$,则 program $\min\|x'\|_T$ subject to $Ax'=y$ 的解 $\hat x$ 满足 $\|\hat x\|_T\le\|x\|_T\le1$,所以 $\hat x,x\in T$ 且 $A\hat x=Ax$。直接应用 Theorem 9.4.4。

Exercise 9.19Noisy constrained optimization
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:证明 noisy constrained recovery bound。

完整证明:最优性给出 $\|A\hat x-y\|_2\le\|Ax-y\|_2=\|w\|_2$。因此 $\|A(\hat x-x)\|_2\le2\|w\|_2$。Matrix deviation 在 $T-T$ 上给出 $\|Ah\|_2\ge\sqrt m\|h\|_2-CK^2w(T)$。合并并取期望得 $\mathbb E\|h\|_2\lesssim(K^2w(T)+\|w\|_2)/\sqrt m$。

Exercise 9.20Unconstrained optimization
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:证明 penalized least squares 在 $\lambda\asymp\|w\|_2^2/\|x\|_T$ 时满足

$$\mathbb E\|\hat x-x\|_2\lesssim\frac{K^2w(T)\|x\|_T+\|w\|_2}{\sqrt m}.$$

完整证明:记 $R=\|x\|_T$,$h=\hat x-x$。若 $R=0$,结论退化为纯噪声情形,按同样论证去掉 $Rw(T)$ 项即可。下面设 $R>0$,并取 $\lambda=c\|w\|_2^2/R$ 到常数因子。由最优性,

$$\|A\hat x-y\|_2^2+\lambda\|\hat x\|_T\le \|Ax-y\|_2^2+\lambda\|x\|_T=\|w\|_2^2+\lambda R.$$

因此 $\|\hat x\|_T\le C R$,且 $\|A\hat x-y\|_2\le C\|w\|_2$。由 $y=Ax+w$,

$$\|Ah\|_2\le \|A\hat x-y\|_2+\|w\|_2\le C\|w\|_2.$$

同时 $\|h\|_T\le\|\hat x\|_T+\|x\|_T\le CR$,所以 $h\in CR\,T$,其中 $T$ 是该 norm 的单位球。对集合 $CR\,T$ 使用 matrix deviation,得到期望意义下的 uniform lower bound

$$\mathbb E\sup_{z\in CR\,T}\left|\|Az\|_2-\sqrt m\|z\|_2\right|\le CK^2 R w(T).$$

于是

$$\sqrt m\,\mathbb E\|h\|_2\le \mathbb E\|Ah\|_2+CK^2Rw(T)\le C\|w\|_2+CK^2Rw(T).$$

两边除以 $\sqrt m$ 即得结论。

Exercise 9.21Lasso
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:把 Exercise 9.20 专门化到 $\ell^1$ norm。

完整证明:若 $x$ 为 $s$-sparse,则 $\|x\|_1\le\sqrt s\|x\|_2$。Exercise 9.20 中 $T=B_1^n$,但有效集合缩放到 $\|x\|_1B_1^n$。Gaussian width 为 $\|x\|_1w(B_1^n)\lesssim\sqrt{s\log n}\|x\|_2$。归一化 $\|x\|_2\le1$ 后得到 $\mathbb E\|\hat x-x\|_2\lesssim(K^2\sqrt{s\log n}+\|w\|_2)/\sqrt m$。

Exercise 9.22Sparse recovery is well posed
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:证明 general position 下唯一性与已知 support 的求解。

完整证明:若两个 $s$-sparse 解 $x,x'$ 满足 $Ax=Ax'$,则 $h=x-x'$ 是 $2s$-sparse 且在 $\ker A$ 中。General position 假设任意 $2s$ 列线性无关;当 $m\ge2s$ 时,这迫使 $h=0$。若 support 已知为 $S$,只需解 $A_Sx_S=y$ 的最小二乘或线性方程;$A_S$ 满列秩时解唯一且可高效计算。

Exercise 9.23$\ell^p$ for $p<1$
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:说明 $\ell^0$ 与 $0<p<1$ 不是 norm,并证明极限。

完整证明:$\|2x\|_0=\|x\|_0$,不满足 norm 的齐次性。对 $0<p<1$,unit ball 非凸,例如 $(1,0)$ 与 $(0,1)$ 在 ball 中,但中点的 $p$-quasi norm 可超过 $1$,三角不等式失败。最后,若 $x_i\ne0$,则 $|x_i|^p\to1$;若 $x_i=0$,则 $|x_i|^p=0$。有限求和后 $\lim_{p\to0+}\|x\|_p^p=\#\{i:x_i\ne0\}=\|x\|_0$。

Exercise 9.24Approximate sparse recovery
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:把 sparse recovery 从 exactly sparse 推广到 approximately sparse。

完整证明:先看 exactly sparse。若 $x$ 是 $s$-sparse,则 $\|x\|_1\le\sqrt s\|x\|_2$。$\ell^1$ minimization 的解 $\hat x$ 满足 $\|\hat x\|_1\le\|x\|_1$,所以 $x,\hat x$ 都属于 $\sqrt s\|x\|_2B_1^n$。对该集合应用 constrained recovery,并用 $w(B_1^n)\lesssim\sqrt{\log n}$,得到

$$\mathbb E\|\hat x-x\|_2\le CK^2\sqrt{\frac{s\log n}{m}}\,\|x\|_2.$$

更精细地用 Exercise 9.26 的 truncated sparse hull,可把 $\log n$ 改成 $\log(en/s)$。

对 approximately sparse 的版本,令 $x_s$ 为 $x$ 的最佳 $s$-term approximation,$r=x-x_s$。把观测写成

$$y=Ax=A x_s+Ar,$$

即把 $Ar$ 视为噪声。使用带噪约束版本:在 $\|z\|_1\le\|x_s\|_1$ 中最小化 residual,得到

$$\mathbb E\|\hat x-x_s\|_2\le C\left(K^2\sqrt{\frac{s\log(en/s)}{m}}\|x_s\|_2+\frac{\mathbb E\|Ar\|_2}{\sqrt m}\right).$$

由于行各向同性,$\mathbb E\|Ar\|_2\le\sqrt m\|r\|_2$。对最佳 $s$-term tail,按坐标大小分块有 $\|r\|_2\le \|r\|_1/\sqrt s$。因此

$$\mathbb E\|\hat x-x\|_2\le C K^2\sqrt{\frac{s\log(en/s)}{m}}\|x\|_2+C\frac{\|x-x_s\|_1}{\sqrt s}.$$

这就是 approximately sparse recovery:第一项是估计 $s$-sparse 主体的随机测量误差,第二项是不可避免的稀疏近似误差。

Exercise 9.25Convexifying sparse vectors
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:证明 $\operatorname{conv}(S_{n,s})\subset T_{n,s}\subset2\operatorname{conv}(S_{n,s})$。

完整证明:第一包含来自 convexity:$S_{n,s}$ 中向量满足 $\|x\|_1\le\sqrt s$ 且 $\|x\|_2\le1$,所以其 convex hull 在 $T_{n,s}$ 中。反向,取 $x\in T_{n,s}$,按坐标绝对值每 $s$ 个一组分块 $I_j$。有 $\sum_{j\ge1}\|x_{I_j}\|_2\le\|x_{I_1}\|_2+s^{-1/2}\sum_{j\ge2}\|x_{I_{j-1}}\|_1\le1+s^{-1/2}\|x\|_1\le2$。每个非零 $x_{I_j}/\|x_{I_j}\|_2$ 属于 $S_{n,s}$,所以 $x$ 是这些点的非负组合,总系数至多 $2$,得到第二包含。

Exercise 9.26Log improvement
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:改进 sparse recovery 的 logarithmic factor。

完整证明:由 Exercise 9.25,$w(T_{n,s})\le2w(S_{n,s})$。对 $S_{n,s}$,按 support union bound:每个 support 上 width 为 $\sqrt s$,support 数为 $\binom ns$,最大值估计给出 $w(S_{n,s})\le C\sqrt{s\log(en/s)}$。代入 constrained recovery 得 $m\gtrsim s\log(en/s)$ 的改进界。

Exercise 9.27Gaussian width of sparse vectors
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:证明 sparse set width 的上下界。

完整证明:上界见 Exercise 9.26。下界取 $\binom ns$ 个 support 中的 packing,或取前 $s$ 大 Gaussian 坐标:$\sup_{x\in S_{n,s}}\langle g,x\rangle=(\sum_{i=1}^s g_{(i)}^2)^{1/2}$,其中 $g_{(i)}$ 为绝对值降序。order statistics 给出期望 $\gtrsim\sqrt{s\log(en/s)}$。convexified set 与 $S_{n,s}$ 同阶由 Exercise 9.25 得到。

Exercise 9.28Garnaev-Gluskin theorem
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:把 $M^*$ bound 对 $B_1^n$ 随机切片给出的粗略 $\sqrt{\log n/m}$ 改进为

$$\mathbb E\operatorname{diam}(B_1^n\cap E)\lesssim \sqrt{\frac{\log(en/m)}{m}},$$

其中 $E=\ker A$ 是随机余维 $m$ 子空间。

完整证明:记 $T=B_1^n$,并对 $0<\rho\le1$ 定义截断 cross-polytope

$$T_\rho=B_1^n\cap \rho B_2^n.$$

第一步是把全局半径问题化为截断半径问题。若某个 $\delta\le\rho$ 满足 $\operatorname{rad}(T_\rho\cap E)\le\delta$,则 $\operatorname{rad}(T\cap E)\le\delta$。证明如下:若存在 $x\in T\cap E$ 且 $\|x\|_2>\delta$,分两种情况。若 $\|x\|_2\le\rho$,则 $x\in T_\rho\cap E$,与截断半径界矛盾;若 $\|x\|_2>\rho$,令 $y=\rho x/\|x\|_2$,则 $y\in E$、$\|y\|_2=\rho>\delta$,且 $\|y\|_1=(\rho/\|x\|_2)\|x\|_1\le1$,所以 $y\in T_\rho\cap E$,仍矛盾。

第二步估计 $T_\rho$ 的 Gaussian width。设 $s=\lceil\rho^{-2}\rceil$。由于

$$T_\rho=B_1^n\cap\rho B_2^n\subset \rho\bigl(\sqrt{s}B_1^n\cap B_2^n\bigr)=\rho T_{n,s},$$

Exercise 9.26 给出

$$w(T_\rho)\le \rho\,w(T_{n,s})\lesssim \rho\sqrt{s\log(en/s)} \lesssim \sqrt{\log(en\rho^2)},$$

这里最后一步用 $s\asymp\rho^{-2}$,并把 ceiling 带来的常数吸收进绝对常数中。

第三步使用高概率 $M^*$ bound。对 $T_\rho$ 有 $\operatorname{rad}(T_\rho)\le\rho$,因此 Exercise 9.13 型 tail 形式给出

$$\operatorname{rad}(T_\rho\cap E)\le C\frac{w(T_\rho)+u\rho}{\sqrt m}$$

以概率至少 $1-2e^{-u^2}$ 成立。选择 $u=c\sqrt m\,\rho$。只要

$$m\rho^2\ge C_0\log(en\rho^2),$$

并取常数 $C_0$ 足够大,就有

$$C\frac{w(T_\rho)+u\rho}{\sqrt m}\le \rho/2.$$

于是事件上 $\operatorname{rad}(T_\rho\cap E)\le\rho/2$,由第一步推出 $\operatorname{rad}(T\cap E)\le\rho/2$。

最后选择参数。令

$$L=\log(en/m),\qquad \rho^2=A\,\frac{L}{m},$$

其中 $A$ 是足够大的绝对常数。若 $\rho>1$,则右侧目标已经是常数阶,结论由 $\operatorname{diam}(B_1^n)\le2$ 平凡成立。下面假设 $\rho\le1$。此时

$$\log(en\rho^2)=\log\!\left(\frac{en}{m}\,A L\right)\le C L,$$

所以 $m\rho^2=AL$ 可以压过 $\log(en\rho^2)$。由上一步,除概率至多 $2e^{-c m\rho^2}\le2e^{-cAL}$ 的坏事件外,$\operatorname{rad}(T\cap E)\le\rho/2$。坏事件上只用 $\operatorname{rad}(T\cap E)\le1$,得到

$$\mathbb E\operatorname{rad}(T\cap E)\le \rho/2+2e^{-cAL}\lesssim \sqrt{\frac{\log(en/m)}{m}}.$$

乘以 $2$ 把半径换成直径,即得

$$\mathbb E\operatorname{diam}(B_1^n\cap E)\lesssim \sqrt{\frac{\log(en/m)}{m}}.$$

这个证明的核心是截断、稀疏凸包宽度估计与 $M^*$ tail bound 的自洽参数选择:选 $\rho$ 时必须让截断切片先落回 $\rho B_2^n$ 内,再把这个半径传回整个 $B_1^n$。

Exercise 9.29Low-rank recovery extensions
状态:完整证明已按当前笔记标准给出闭合推导。

证明目标:把 Corollary 9.4.11 扩展到 nuclear norm minimization、approximately low-rank 和矩形矩阵。

完整证明:(a) 若真实矩阵 $X$ 满足 $\operatorname{rank}(X)\le r$ 且 $\|X\|_F\le1$,则 $\|X\|_*\le\sqrt r$。nuclear norm minimization 的解 $\hat X$ 满足 $\|\hat X\|_*\le\|X\|_*\le\sqrt r$。因此 $X,\hat X$ 都落在 $\sqrt r B_*$ 中,且 $A(\hat X-X)=0$。对 $T=\sqrt r B_*$ 应用 constrained recovery,

$$\mathbb E\|\hat X-X\|_F\le \frac{Cw(\sqrt r B_*)}{\sqrt m}.$$

由于 nuclear norm ball 的 polar 是 operator norm ball,

$$w(\sqrt r B_*)=\sqrt r\,\mathbb E\|G\|\lesssim\sqrt{rd}$$

在 $d\times d$ 情形成立,故得到 $\mathbb E\|\hat X-X\|_F\lesssim\sqrt{rd/m}$。

(b) 若 $X$ 不是 rank-$r$,令 $X_r$ 为最佳 rank-$r$ approximation,$R=X-X_r$。把观测写成 $y_i=\langle A_i,X_r\rangle+\langle A_i,R\rangle$,把尾部看成噪声。由带噪 constrained recovery,

$$\mathbb E\|\hat X-X_r\|_F\lesssim \sqrt{\frac{rd}{m}}+\frac{\mathbb E\|A(R)\|_2}{\sqrt m}.$$

各向同性给出 $\mathbb E\|A(R)\|_2\le\sqrt m\|R\|_F$,而奇异值尾部满足 $\|R\|_F\le\|R\|_*/\sqrt r$。于是误差增加

$$\frac{\|X-X_r\|_*}{\sqrt r}.$$

(c) 对 $d_1\times d_2$ 矩阵,仍有 $w(\sqrt r B_*)=\sqrt r\,\mathbb E\|G\|$,而 Gaussian 矩阵 operator norm 满足 $\mathbb E\|G\|\lesssim\sqrt{d_1}+\sqrt{d_2}$。因此

$$w(\sqrt r B_*)\lesssim\sqrt{r(d_1+d_2)},$$

代入 constrained recovery 得

$$\mathbb E\|\hat X-X\|_F\lesssim\sqrt{\frac{r(d_1+d_2)}{m}},$$

approximately low-rank 的矩形版本同理再加 $\|X-X_r\|_*/\sqrt r$。

Exercise 9.30Geometry of exact sparse recovery
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:解释 Theorem 9.5.1 的 tangent cone 图像。

完整证明:见 [隐藏验证](#proof-check-9-5-tangent-cone)。证明说明 basis pursuit 成功当且仅当 $\ker A$ 不包含任何从 $x$ 进入 $\ell^1$ ball 的非零方向;归一化后,就是 $\ker A$ 避开 tangent cone 的 spherical part。

Exercise 9.31Noisy measurements
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:把 Theorem 9.5.1 的 noiseless exact sparse recovery 改成 noisy measurements $y=Ax+w$ 的稳定恢复。

完整证明:设已知噪声上界 $\eta\ge\|w\|_2$,考虑 basis pursuit denoising:

$$\min_{z\in\mathbb R^n}\|z\|_1\quad\text{subject to}\quad \|Az-y\|_2\le\eta.$$

真实 $x$ 可行,所以解 $\hat x$ 满足 $\|\hat x\|_1\le\|x\|_1$。令 $h=\hat x-x$,若 $S=\operatorname{supp}(x)$ 且 $|S|=s$,与 Lemma 9.5.2 同样的三角不等式给出 cone constraint

$$\|h_{S^c}\|_1\le \|h_S\|_1.$$

因此 $\|h\|_1\le2\sqrt s\|h\|_2$,也就是说 $h/\|h\|_2$ 落在 approximate sparse spherical set

$$T_s=\{u\in S^{n-1}:\|u\|_1\le2\sqrt s\}.$$

另一方面,$\hat x$ 和 $x$ 都满足噪声约束,所以

$$\|Ah\|_2\le\|A\hat x-y\|_2+\|Ax-y\|_2\le2\eta.$$

当 $m\ge CK^4s\log(en/s)$ 时,matrix deviation 或 RIP 证明给出 uniform lower bound:对所有 $u\in T_s$,$\|Au\|_2\ge c\sqrt m$,概率至少 $1-2e^{-cm/K^4}$。若 $h\ne0$,代入 $u=h/\|h\|_2$ 得

$$c\sqrt m\|h\|_2\le\|Ah\|_2\le2\eta,$$

从而

$$\|\hat x-x\|_2\le C\frac{\eta}{\sqrt m}.$$

若 $x$ 只是 approximately sparse,把 $S$ 取为最大 $s$ 个坐标集合,则 cone constraint 变为 $\|h_{S^c}\|_1\le\|h_S\|_1+2\|x_{S^c}\|_1$,同一论证额外产生 tail 项 $C\|x_{S^c}\|_1/\sqrt s$。

Exercise 9.32Nullspace property
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:证明 NSP 等价于所有 $s$-sparse exact recovery,并证明随机矩阵满足。

完整证明:若 NSP 成立,任取 $s$-sparse $x$ 和非零 $h\in\ker A$,令 $S=\operatorname{supp}(x)$。则 $\|x+h\|_1\ge\|x\|_1-\|h_S\|_1+\|h_{S^c}\|_1>\|x\|_1$,故 $x$ 是唯一 minimizer。反向,若 NSP 失败,存在 $h$ 与 $S$ 使 $\|h_S\|_1\ge\|h_{S^c}\|_1$,取 $x=-h_S$,则 $x+h=h_{S^c}$ 也是可行且 $\ell^1$ 不大,唯一恢复失败。随机矩阵部分由 Theorem 9.5.1 对所有 supports 的 uniform statement 给出。

Exercise 9.33Random projections satisfy RIP
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:证明 Grassmannian projections 的 sparse RIP 与 exact recovery。

完整证明:对 sparse set $S_{n,s}$ 使用 Exercise 9.6 的 random projection deviation。若 $m\gtrsim s\log(en/s)$,则所有 $s$-sparse unit vectors 满足 $\|Px\|_2\approx\sqrt{m/n}$。归一化 $\sqrt{n/m}P$ 后得到 RIP。再应用 Theorem 9.5.6 或 nullspace property,得到 random projections 下的 exact recovery。

Exercise 9.34Subadditivity
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:证明 $f(x)-f(y)\le f(x-y)$。

完整证明:由 subadditivity,$f(x)=f((x-y)+y)\le f(x-y)+f(y)$。移项即得 $f(x)-f(y)\le f(x-y)$。

Exercise 9.35Anisotropic general deviation
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:推广 Theorem 9.6.3 到 $N(0,\Sigma)$ rows。

完整证明:写 $A=G\Sigma^{1/2}$,其中 $G$ 为 standard Gaussian matrix。于是 $f(Ax)=f(G(\Sigma^{1/2}x))$。对集合 $\Sigma^{1/2}T$ 应用 Theorem 9.6.3,得到 $Cb\gamma(\Sigma^{1/2}T)$。

Exercise 9.36General deviation high probability
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:证明 Theorem 9.6.3 的 tail version。

完整证明:Theorem 9.6.4 给出 increments,generic chaining high-probability 或 Talagrand comparison tail 给出

$$\sup_{x\in T}|f(Ax)-\mathbb Ef(Ax)|\le Cb(\gamma(T)+u\operatorname{rad}(T))$$

概率至少 $1-2e^{-u^2}$。

Exercise 9.37JL for general norms
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:用 Theorem 9.6.3 得到一般 norm 的 JL。

完整证明:令 $f=\|\cdot\|$,假设 $f(z)\le b\|z\|_2$。对 normalized differences $T$ 应用 general deviation:$f(Az)$ 同时接近 $\mathbb Ef(Az)=\|z\|_2\mathbb Ef(g)$。若 $T$ 来自 $N$ 点集,则 $\gamma(T)\le C\sqrt{\log N}$。当 $m$ 或 norm 参数使 $b\sqrt{\log N}$ 小于 $\varepsilon\mathbb Ef(g)$ 时,得到一般 norm 中的 JL embedding。

Exercise 9.38JL into $\ell^1$
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:证明 Gaussian map 到 $\ell^1$ 的 JL。

完整证明:取 $f(z)=\|z\|_1$,则 $b=\sqrt m$,且 $\mathbb E\|g\|_1=m\sqrt{2/\pi}$。对 normalized differences 应用 Theorem 9.6.3 的高概率版本,偏差为 $C\sqrt m\sqrt{\log N}$。除以 $m\sqrt{2/\pi}$ 后相对误差为 $C\sqrt{\log N/m}$。若 $m\ge C(\varepsilon)\log N$,则 $Q=\sqrt{\pi/2}\,A/m$ 满足题设双边界。

Exercise 9.39JL into $\ell^\infty$
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:证明 Gaussian map 到 $\ell^\infty$ 的 embedding。

完整证明:取 $f(z)=\|z\|_\infty$,则 $b=1$,且 $\mathbb E\|g\|_\infty\asymp\sqrt{\log m}$。对 $N^2$ 个 normalized differences 应用 general deviation,偏差为 $C\sqrt{\log N}$。要使其小于 $\varepsilon\sqrt{\log m}$,需 $\log m\ge C(\varepsilon)\log N$,即 $m\ge N^{C(\varepsilon)}$。取 $Q=C(\log m)^{-1/2}A$ 得题设结论。

Exercise 9.40Duality
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:证明 support function 与 Euclidean ball inclusion 等价。

完整证明:见 [隐藏验证](#proof-check-9-7-support-duality)。关键是 $h_{\operatorname{conv}(V)}=h_V$,以及 closed convex body 可由其所有 supporting halfspaces 表示。

Exercise 9.41High-probability Dvoretzky-Milman
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:给 Theorem 9.7.2 加入 tail 参数,得到随机 Gaussian projection 后凸包夹在两个 Euclidean balls 之间的高概率版本。

完整证明:令 $A$ 为 $m\times n$ Gaussian matrix。对任意 $y\in S^{m-1}$,support function 满足

$$h_{AT}(y)=\sup_{x\in T}\langle Ax,y\rangle=\sup_{x\in T}\langle A^\top y,x\rangle.$$

由于 $A^\top y\sim N(0,I_n)$,有 $\mathbb Eh_{AT}(y)=w(T)$。Theorem 9.7.1 的高概率版本应用到索引集 $S^{m-1}$ 给出:对所有 $u\ge0$,以概率至少 $1-2e^{-u^2}$,

$$\sup_{y\in S^{m-1}}|h_{AT}(y)-w(T)|\le C\operatorname{rad}(T)(\sqrt m+u).$$

$$r_-=w(T)-C\operatorname{rad}(T)(\sqrt m+u),\qquad r_+=w(T)+C\operatorname{rad}(T)(\sqrt m+u).$$

如果 $r_-<0$,内含球半径按 $0$ 理解。对每个 $y\in S^{m-1}$,上式等价于

$$r_-\le h_{AT}(y)\le r_+.$$

support function 的对偶刻画说明 $h_{\operatorname{conv}(AT)}=h_{AT}$,并且对闭凸集 $K$,$K\subset r_+B_2^m$ 等价于 $h_K(y)\le r_+$ 对所有 $y\in S^{m-1}$,而 $r_-B_2^m\subset K$ 等价于 $h_K(y)\ge r_-$ 对所有 $y$。因此

$$r_-B_2^m\subset \operatorname{conv}(AT)\subset r_+B_2^m$$

以概率至少 $1-2e^{-u^2}$ 成立。这就是 Theorem 9.7.2 的 high-probability 形式。

Exercise 9.42Gaussian cloud is nearly round
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:证明 $n$ 个 Gaussian points 的 convex hull 近似 Euclidean ball。

完整证明:令 $T=\{e_1,\dots,e_n\}\subset\mathbb R^n$,Gaussian matrix $A$ 的列就是 $g_1,\dots,g_n$。有 $w(T)=\mathbb E\max_i g_i\asymp\sqrt{\log n}$,$\operatorname{rad}(T)=1$。Theorem 9.7.2 给出 $\operatorname{conv}\{g_i\}$ 夹在半径 $\sqrt{\log n}\pm C\sqrt m$ 的 balls 中。若 $m\le c\log n$ 且 $c$ 小,误差项被吸收,得到半径 $\asymp\sqrt{\log n}$ 的近圆性。

Exercise 9.43Actual random projections
状态:依赖前章结论证明依赖前章定理或标准工具,本页给出调用链。

证明目标:把 Gaussian projection 版 Dvoretzky-Milman 改成真正的随机正交投影 $P:\mathbb R^n\to E$,其中 $E$ 在 Grassmannian 上均匀。

完整证明:令 $G$ 为 $m\times n$ standard Gaussian matrix。它的 row space 在 Grassmannian $G_{n,m}$ 上均匀分布。把 $G$ 写成极分解

$$G=R P,$$

其中 $P$ 是到随机 $m$ 维 row space 的正交投影再识别到 $\mathbb R^m$,$R=(GG^\top)^{1/2}$ 是 $\mathbb R^m$ 上的随机正定算子;$P$ 的方向与径向部分独立。Wishart 奇异值集中给出,当 $m\le c\varepsilon^2 n$ 时,以高概率

$$(1-\varepsilon)\sqrt n\,I_m\preceq R\preceq(1+\varepsilon)\sqrt n\,I_m.$$

也就是说,$G T=R(PT)$ 与 $\sqrt n\,PT$ 只差一个 $(1\pm\varepsilon)$ 的线性畸变。

对 Gaussian projection $G$ 应用 Theorem 9.7.2 或 Exercise 9.41 的高概率版本,若 $m\le c\varepsilon^2 w(T)^2/\operatorname{rad}(T)^2$,则

$$(1-\varepsilon)w(T)B_2^m\subset \operatorname{conv}(GT)\subset(1+\varepsilon)w(T)B_2^m.$$

用上面的奇异值比较把 $G T$ 除以 $\sqrt n$ 转换为 $PT$。由于 $w_s(T)=w(T)/\sqrt n$,并吸收两个 $(1\pm\varepsilon)$ 因子,得到

$$(1-C\varepsilon)w_s(T)B_2^m\subset \operatorname{conv}(PT)\subset(1+C\varepsilon)w_s(T)B_2^m.$$

重新命名 $\varepsilon$,得到题设的 true random projection 版本:

$$(1-\varepsilon)B\subset\operatorname{conv}(PT)\subset(1+\varepsilon)B,$$

其中 $B$ 是半径 $w_s(T)$ 的 Euclidean ball。

易混点

易混点 正确理解
Matrix deviation 是否只是 operator norm bound 不是。operator norm 是 $T=S^{n-1}$ 的特例。
$w(T)$ 与 $\gamma(T)$ 是否相同 对对称或含原点情形同阶;绝对值过程更自然使用 $\gamma(T)$。
$M^*$ 与 escape theorem 是否同一个结论 一个控制截面直径,一个控制是否相交。
Exact recovery 是否只靠 RIP 不只。Escape theorem / tangent cone 是更几何的证明。
Dvoretzky-Milman 是否说原集合是圆的 不是。它说随机低维像或截面在合适维度下近似圆。

公式卡片

场景 公式
Matrix deviation $\mathbb E\sup_{x\in T}|\|Ax\|_2-\sqrt m\|x\|_2|\le CK^2\gamma(T)$
Covariance $\mathbb E\|\Sigma_m-\Sigma\|\le CK^4(\sqrt{r/m}+r/m)\|\Sigma\|$
$M^*$ $\mathbb E\operatorname{diam}(T\cap\ker A)\le CK^2w(T)/\sqrt m$
Escape $m\gtrsim K^4w(T)^2\Rightarrow T\cap\ker A=\varnothing$
Constrained recovery $\mathbb E\|\hat x-x\|_2\lesssim K^2w(T)/\sqrt m$
Sparse width $w(S_{n,s})\asymp\sqrt{s\log(en/s)}$
General deviation $\mathbb E\sup_T|f(Ax)-\mathbb Ef(Ax)|\le Cb\gamma(T)$
Dvoretzky $r_-B_2^m\subset\operatorname{conv}(AT)\subset r_+B_2^m$

学习检查表

检查点 你应能完成的动作
Theorem 9.1.2 解释为什么 squared norm 差值可由 Bernstein 控制。
Matrix deviation 从 subgaussian increments 一步推出主定理。
Covariance estimation 把 sample covariance 写成 ellipsoid 上 quadratic deviation。
$M^*$ 用 $T-T$ 和 $\ker A$ 让 $\|A(x-y)\|$ 消失。
Escape 用 high-probability deviation 与 $\|x\|=1$ 制造矛盾。
Sparse recovery 从 $\ell^1$ optimality 推出 cone constraint。
RIP 把 sparse set 的 Gaussian width 代入 matrix deviation。
Dvoretzky 用 support function 判断 convex body 是否接近 Euclidean ball。

后续衔接

第 9 章完成了第 6-9 章的随机矩阵与随机过程工具链:Hanson-Wright/decoupling、Gaussian comparison、chaining、matrix deviation 逐层衔接。后续若继续制作第 10 章,可以把这里的 matrix deviation 视为进入更高级随机矩阵与几何泛函分析应用的基础模板。