第 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。