HDP 读书笔记
设置
字号 标准
精校翻译 Ch.9 矩阵偏差

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

第 9 章集合上的随机矩阵偏差

一个 $m\times n$ 随机矩阵会怎样作用在一般集合 $T\subset\mathbb R^n$ 上?第 9.1 节证明本章的核心结果:矩阵偏差不等式。随后我们把它应用到一系列高维问题中,其中既有前面已经见过的问题,也有新的问题。

第 9.2 节会快速推出随机矩阵的双边界、集合随机投影的更尖锐界、低维数据的协方差估计,以及 Johnson-Lindenstrauss 引理的无限集合版本。

第 9.3 节包含两个优雅结果:随机切片如何缩小高维集合,也就是 $M^*$ 界;以及随机子空间如何完全避开某些集合,也就是逃逸定理。

第 9.4-9.5 节把这些想法用于数据科学中的基础问题:学习结构化的高维线性模型。

第 9.6 节把矩阵偏差不等式推广到一般范数;第 9.7 节再用它强化 Chevet 不等式,并推出 Dvoretzky-Milman 定理。后者说明:高维集合的随机低维投影看起来近似圆。

习题会继续展开矩阵偏差和随机过程偏差界,例如 Exercise 9.5;也会进入高维估计方法,包括稀疏回归中的 Lasso,见 Exercises 9.20-9.21;还包括 Garnaev-Gluskin 定理关于 cross-polytope 的随机切片,见 Exercise 9.28;关于其他 $\ell^p$ 球的 Exercise 9.14;关于一般范数 Johnson-Lindenstrauss 引理的 Exercises 9.37-9.39,等等。

Tips:第 9 章把前面工具合流到应用:矩阵偏差不等式是核心引擎,$M^*$ 界和逃逸定理是几何翻译器,稀疏/低秩恢复是数据科学应用,Dvoretzky-Milman 则展示这些工具在凸几何中的终点。

9.1 矩阵偏差不等式

Tips:矩阵偏差不等式的作用是把“固定向量上的范数集中”升级成“整个集合上的一致控制”。读证明时重点看两个层次:固定方向的次高斯增量,以及用第 8 章的随机过程工具控制上确界。

取一个 $m\times n$ 随机矩阵 $A$,其行独立、各向同性且次高斯。范数集中(Theorem 3.1.1)告诉我们:对任意固定向量 $x\in\mathbb R^n$,近似

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

以高概率成立。

现在问一个更大的问题:是否可以让 (9.1) 同时对许多向量 $x\in\mathbb R^n$ 成立?为了量化“许多”,任取一个有界集合 $T\subset\mathbb R^n$,并要求该近似同时对所有 $x\in T$ 成立。结果表明,最大误差大约由 $\gamma(T)$ 控制;$\gamma(T)$ 是 $T$ 的高斯复杂度,它与第 7.5.3 节的高斯宽度密切相关。

Theorem 9.1.1矩阵偏差不等式

令 $A$ 为 $m\times n$ 随机矩阵,其行 $A_i$ 独立、各向同性且次高斯。则对任意子集 $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)$ 是高斯复杂度,$K=\max_i\|A_i\|_{\psi_2}$。

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

证明计划是把

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

看成由 $x\in\mathbb R^n$ 索引的随机过程,然后从 Talagrand 比较不等式(Corollary 8.5.8)推出结论。为此只需要检查这个随机过程具有次高斯增量。

Theorem 9.1.2次高斯增量

令 $A$ 为 $m\times n$ 随机矩阵,其行 $A_i$ 独立、各向同性且次高斯。那么随机过程 (9.2) 具有次高斯增量:

$$ \|Z_x-Z_y\|_{\psi_2}\le CK^2\|x-y\|_2, \qquad x,y\in\mathbb R^n. \tag{9.3} $$

这里 $K=\max_i\|A_i\|_{\psi_2}$。

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

一旦 Theorem 9.1.2 成立,就把它代入 Talagrand 比较不等式(Corollary 8.5.8,更精确地说是 Exercise 8.37(a)),得到

$$ \mathbb E\sup_{x\in T}|Z_x| \le CK^2\gamma(T), $$

这正是 Theorem 9.1.1。因此剩下的任务只是证明 Theorem 9.1.2;这会容易一些,因为它只涉及固定的 $x,y$。

Proof of Theorem 9.1.2从固定向量增量到全局偏差

证明稍长,但思路是从简单情形逐步推广。

步骤 1:单位向量 $x$ 与零向量 $y$。 假设 $\|x\|_2=1$ 且 $y=0$。这时 (9.3) 变成

$$ \left\|\|Ax\|_2-\sqrt m\right\|_{\psi_2}\le CK^2. \tag{9.4} $$

随机向量 $Ax\in\mathbb R^m$ 的坐标为 $\langle A_i,x\rangle$,它们独立、次高斯,并且由各向同性性有 $\mathbb E\langle A_i,x\rangle^2=1$。于是 (9.4) 直接来自范数集中(Theorem 3.1.1)。

步骤 2:单位向量 $x,y$ 与平方过程。 现在假设 $\|x\|_2=\|y\|_2=1$。此时 (9.3) 等价于

$$ \left\|\|Ax\|_2-\|Ay\|_2\right\|_{\psi_2} \le CK^2\|x-y\|_2. \tag{9.5} $$

平方 $\ell^2$ 范数没有根号,更容易处理。直觉上,若 $\|Ax\|_2$ 和 $\|Ay\|_2$ 都大约为 $\sqrt m$,那么

$$ \begin{aligned} \|Ax\|_2^2-\|Ay\|_2^2 &=\bigl(\|Ax\|_2+\|Ay\|_2\bigr) \bigl(\|Ax\|_2-\|Ay\|_2\bigr)\\ &\approx \sqrt m\,\|x-y\|_2. \end{aligned} \tag{9.6} $$

展开矩阵乘法:

$$ \|Ax\|_2^2-\|Ay\|_2^2 =\sum_{i=1}^m \left(\langle A_i,x\rangle^2-\langle A_i,y\rangle^2\right) =\sum_{i=1}^m \langle A_i,x+y\rangle\langle A_i,x-y\rangle. $$

两边除以 $\|x-y\|_2$,记

$$ \Delta := \frac{\|Ax\|_2^2-\|Ay\|_2^2}{\|x-y\|_2} = \sum_{i=1}^m\langle A_i,u\rangle\langle A_i,v\rangle, \tag{9.7} $$

其中

$$ u:=x+y,\qquad v:=\frac{x-y}{\|x-y\|_2}. $$

目标是说明 $|\Delta|\lesssim\sqrt m$ 以高概率成立。式 (9.7) 是独立随机变量之和。每一项均值为零,因为

$$ \langle A_i,u\rangle\langle A_i,v\rangle = \frac{\langle A_i,x\rangle^2-\langle A_i,y\rangle^2}{\|x-y\|_2}, $$

而各向同性性给出 $\mathbb E[\langle A_i,x\rangle^2-\langle A_i,y\rangle^2]=1-1=0$。并且这些项是次指数:由 Lemma 2.8.6 和 $A_i$ 的次高斯假设,

$$ \begin{aligned} \|\langle A_i,u\rangle\langle A_i,v\rangle\|_{\psi_1} &\le \|\langle A_i,u\rangle\|_{\psi_2} \|\langle A_i,v\rangle\|_{\psi_2}\\ &\le K\|u\|_2\cdot K\|v\|_2 \le 2K^2. \end{aligned} $$

最后一步用了 $\|u\|_2\le\|x\|_2+\|y\|_2\le2$ 且 $\|v\|_2=1$。Bernstein 不等式(Theorem 2.9.1)于是给出:对任意 $0\le t\le\sqrt m$,

$$ \mathbb P\{|\Delta|\ge t\sqrt m\} \le 2\exp\left[ -c\min\left(\frac{t^2}{K^4},\frac{t\sqrt m}{K^2}\right) \right] \le 2\exp\left(-\frac{c_1t^2}{K^4}\right). \tag{9.8} $$

步骤 3:单位向量 $x,y$ 与原始过程。 现在去掉平方,证明单位向量情形的 (9.5)。根据次高斯范数的定义,(9.5) 可写成尾界

$$ p(s):= \mathbb P\left\{ \frac{\left|\|Ax\|_2-\|Ay\|_2\right|}{\|x-y\|_2} \ge s \right\} \le 4\exp\left(-\frac{cs^2}{K^4}\right), \qquad s\gt 0. \tag{9.9} $$

当 $s\le2\sqrt m$ 时,利用步骤 2。将定义 $p(s)$ 的不等式两边乘以 $\|Ax\|_2+\|Ay\|_2$,并用 $\Delta$ 的定义,得到

$$ p(s) \le \mathbb P\{|\Delta|\ge s\|Ax\|_2\}. $$

由 (9.4),$\|Ax\|_2$ 以高概率接近 $\sqrt m$。分成两种情形:常见事件 $\|Ax\|_2\ge\sqrt m/2$,以及罕见事件 $\|Ax\|_2\lt \sqrt m/2$。于是

$$ p(s)\le \mathbb P\left\{|\Delta|\ge\frac{s\sqrt m}{2}\right\} + \mathbb P\left\{\|Ax\|_2\lt \frac{\sqrt m}{2}\right\} =:p_1(s)+p_2(s). $$

步骤 2 控制 $p_1$,步骤 1 控制 $p_2$,合起来得到 (9.9)。

当 $s\gt 2\sqrt m$ 时,由三角不等式

$$ \left|\|Ax\|_2-\|Ay\|_2\right| \le \|A(x-y)\|_2. $$

令 $u=(x-y)/\|x-y\|_2$,则

$$ p(s) \le \mathbb P\{\|Au\|_2\ge s\} \le \mathbb P\{\|Au\|_2-\sqrt m\ge s/2\} \le 2\exp\left(-\frac{cs^2}{K^4}\right), $$

其中最后一步再次使用步骤 1。因此单位向量情形成立。

步骤 4:一般情形。 现在令 $x,y\in\mathbb R^n$ 任意。由缩放,不妨假设 $\|x\|_2=1$ 且 $\|y\|_2\ge1$。把 $y$ 投影到单位球面,记

$$ \bar y:=\frac{y}{\|y\|_2}. $$

三角不等式给出

$$ \|Z_x-Z_y\|_{\psi_2} \le \|Z_x-Z_{\bar y}\|_{\psi_2} + \|Z_{\bar y}-Z_y\|_{\psi_2}. $$

第一项由步骤 3 控制,因为 $x$ 和 $\bar y$ 都是单位向量:

$$ \|Z_x-Z_{\bar y}\|_{\psi_2} \le CK^2\|x-\bar y\|_2. $$

第二项利用 $\bar y$ 与 $y$ 共线以及齐次性:

$$ \|Z_{\bar y}-Z_y\|_{\psi_2} = \|\bar y-y\|_2\,\|Z_{\bar y}\|_{\psi_2}. $$

由于 $\bar y$ 是单位向量,步骤 1 给出 $\|Z_{\bar y}\|_{\psi_2}\le CK^2$。因此

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

这看起来不妙,因为我们希望右边被 $\|x-y\|_2$ 控制,而通常三角不等式方向相反。幸运的是,在当前几何关系中,三角不等式可以近似反向使用:

$$ \|x-\bar y\|_2+\|\bar y-y\|_2 \le \sqrt2\,\|x-y\|_2. $$

这正是 Exercise 9.1 要验证的内容。代入 (9.10),得到

$$ \|Z_x-Z_y\|_{\psi_2} \le \sqrt2\,CK^2\|x-y\|_2, $$

Theorem 9.1.2 得证。

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矩阵均值偏差

一个快速的 centering 技巧可以把 Theorem 9.1.1 变成围绕均值 $\mathbb E\|Ax\|_2$ 的 deviation inequality。请在 Exercise 9.2 中验证。

查看 Exercise 9.2
Remark 9.1.4矩阵偏差:高概率界

我们只把 Theorem 9.1.1 表述为期望界;但借助 Talagrand inequality 的高概率版本(见 Exercise 8.37(b)),它会自动升级为尾界。对任意 $u\ge0$,事件

$$ \sup_{x\in T} \left|\|Ax\|_2-\sqrt m\|x\|_2\right| \le CK^2\left[w(T)+u\,\operatorname{rad}(T)\right] \tag{9.11} $$

以至少 $1-2\exp(-u^2)$ 的概率成立。这里 $\operatorname{rad}(T)=\sup_{x\in T}\|x\|_2$ 是 $T$ 的半径。你可以检查为什么 (9.11) 推出期望界。

Remark 9.1.5矩阵平方偏差

如果关心二次过程 $\|Ax\|_2^2$ 的偏差,也可以容易地从 Theorem 9.1.1 推出:

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

请在 Exercise 9.3 中检查这个推导。

查看 Exercise 9.3

作为练习,可以把矩阵偏差不等式推广到经验过程(Exercise 9.5),并证明随机投影的一个版本(Exercise 9.6,这题有挑战性)。

9.2 随机矩阵、协方差估计与 Johnson-Lindenstrauss

矩阵偏差不等式有很多直接后果。下面几节会依次快速说明。

9.2.1 随机矩阵的奇异值

把矩阵偏差不等式应用于欧氏单位球面 $T=S^{n-1}$,就得到第 4.6 节的奇异值界。

快速检查如下:对球面有

$$ \operatorname{rad}(T)=1,\qquad w(T)\le\sqrt n. $$

因此矩阵偏差不等式 (9.11) 说明事件

$$ \sqrt m-CK^2(\sqrt n+u) \le \|Ax\|_2 \le \sqrt m+CK^2(\sqrt n+u), \qquad x\in S^{n-1} $$

以至少 $1-2\exp(-u^2)$ 的概率成立。利用 (4.14),这等价于

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

从而重新得到 Theorem 4.6.1。我们此前用另一种方法证明过它。

9.2.2 集合的随机投影

Proposition 9.2.1集合随机投影的大小

令 $T\subset\mathbb R^n$ 为有界集合,令 $A$ 为 $m\times n$ 矩阵,其行 $A_i$ 独立、各向同性且次高斯。则缩放矩阵

$$ P=\frac1{\sqrt n}A $$

作为次高斯投影满足

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

这里 $K=\max_i\|A_i\|_{\psi_2}$,$w_s(T)$ 是 $T$ 的球面宽度。

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

这比我们之前的界(Theorem 7.6.1 和 Exercise 7.25)略尖锐:在 $\sqrt{m/n}$ 前面没有额外常数。

Proof由半径界推出直径界

Theorem 9.1.1 与三角不等式给出

$$ \mathbb E\sup_{x\in T}\|Ax\|_2 \le \sqrt m\sup_{x\in T}\|x\|_2 + CK^2\gamma(T). $$

用半径写就是

$$ \mathbb E\operatorname{rad}(AT) \le \sqrt m\,\operatorname{rad}(T)+CK^2\gamma(T). $$

把这个界应用到差集 $T-T$,得到

$$ \mathbb E\operatorname{diam}(AT) \le \sqrt m\,\operatorname{diam}(T)+2CK^2w(T), $$

其中使用了 Lemma 7.5.11(a),从高斯复杂度转到高斯宽度。两边除以 $\sqrt n$ 即得结论。

现在可以尝试两件事:把 Proposition 9.2.1 的期望界升级为高概率界(Exercise 9.7);以及在 Exercise 9.8 中分析真正的随机投影,而不是这里的 “高斯” 投影。

9.2.3 低维分布的协方差估计

Theorem 9.2.2低维分布的协方差估计

若 $X$ 是次高斯随机向量,$\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 完整证明

证明回到第 4.7 节的协方差估计问题。我们希望从 $m$ 个独立样本估计 $n$ 维分布的协方差矩阵

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

一般情形中,$m=O(n\log n)$ 个样本足够(第 5.6 节);对次高斯分布,$m=O(n)$ 足够(第 4.7 节)。如果分布近似低维,例如集中在某个 $r$ 维子空间附近,Remark 5.6.3 说明 $m=O(r\log n)$ 足够。下面的定理说明:对次高斯分布,实际上 $m=O(r)$ 就够。

Proof把协方差估计化为 ellipsoid 上的二次偏差

像 Theorem 4.7.1 的证明一样,把分布放到各向同性位置:写

$$ X=\Sigma^{1/2}Z, \qquad X_i=\Sigma^{1/2}Z_i, $$

其中 $Z$ 和 $Z_i$ 是各向同性。于是

$$ \begin{aligned} \|\Sigma_m-\Sigma\| &= \|\Sigma^{1/2}R_m\Sigma^{1/2}\|, \qquad R_m=\frac1m\sum_{i=1}^mZ_iZ_i^{\mathsf T}-I_n\\ &= \max_{x\in S^{n-1}} \left|x^{\mathsf T}\Sigma^{1/2}R_m\Sigma^{1/2}x\right|\\ &= \max_{x\in T}|x^{\mathsf T}R_mx|, \qquad T:=\Sigma^{1/2}S^{n-1}\\ &= \frac1m \max_{x\in T} \left| \|Ax\|_2^2-m\|x\|_2^2 \right|, \end{aligned} $$

其中 $A$ 是以 $Z_i$ 为行的 $m\times n$ 矩阵。与 Theorem 4.7.1 中一样,$Z_i$ 各向同性,并满足 $\|Z_i\|_{\psi_2}\lesssim1$;为了简化记号,这里把对 $K$ 的依赖隐藏起来。

对 $A$ 使用矩阵偏差不等式的二次型版本(Exercise 9.3),得到

$$ \mathbb E\|\Sigma_m-\Sigma\| \lesssim \frac1m \left(\gamma(T)^2+\sqrt m\,\operatorname{rad}(T)\gamma(T)\right). $$

椭球 $T=\Sigma^{1/2}S^{n-1}$ 的半径和高斯复杂度满足

$$ \operatorname{rad}(T)=\|\Sigma\|^{1/2}, \qquad \gamma(T)\le(\operatorname{tr}\Sigma)^{1/2}. $$

因此

$$ \mathbb E\|\Sigma_m-\Sigma\| \lesssim \frac1m \left( \operatorname{tr}\Sigma + \sqrt{m\|\Sigma\|\operatorname{tr}\Sigma} \right). $$

代入 $\operatorname{tr}(\Sigma)=r\|\Sigma\|$ 并化简,即得定理。

Remark 9.2.3协方差估计:高概率保证

和之前一样(见 Remarks 4.7.3 与 5.6.5),Theorem 9.2.2 的期望界可以升级为高概率界。对任意 $u\ge0$,有

$$ \|\Sigma_m-\Sigma\| \le CK^4 \left( \sqrt{\frac{r+u}{m}}+\frac{r+u}{m} \right) \|\Sigma\| $$

以至少 $1-2e^{-u}$ 的概率成立。请在 Exercise 9.9 中证明这个结论。

查看学习笔记:Exercise 9.9

9.2.4 Johnson-Lindenstrauss 引理(无限集合)

矩阵偏差不等式很快可以重新推出第 5.3 节的 Johnson-Lindenstrauss 引理,并把它推广到一般甚至无限集合。

先看有限集合。固定一个 $N$ 点集 $\mathcal X\subset\mathbb R^n$,考虑归一化差集

$$ T= \left\{ \frac{x-y}{\|x-y\|_2}: x,y\in\mathcal X,\ x\ne y \right\}. $$

和 Example 7.5.9 一样,$T$ 的高斯复杂度满足

$$ \gamma(T)\le C\sqrt{\log N}. \tag{9.12} $$

Theorem 9.1.1 说明,以高概率,

$$ \sup_{x,y\in\mathcal X} \left| \frac{\|Ax-Ay\|_2}{\|x-y\|_2} -\sqrt m \right| \lesssim \sqrt{\log N}. $$

改写后,随机矩阵 $Q=A/\sqrt m$ 是 $\mathcal X$ 上的近似 isometry:

$$ (1-\varepsilon)\|x-y\|_2 \le \|Qx-Qy\|_2 \le (1+\varepsilon)\|x-y\|_2, \qquad x,y\in\mathcal X, $$

其中 $\varepsilon\asymp\sqrt{\log(N)/m}$。等价地,如果固定 $\varepsilon\gt 0$ 并取

$$ m\gtrsim \varepsilon^{-2}\log N, $$

则 $Q$ 以高概率成为 $\mathcal X$ 上的 $\varepsilon$-等距映射。这就恢复了经典 Johnson-Lindenstrauss 引理(Theorem 5.3.1)的一个版本。

上面的论证并不真正依赖 $\mathcal X$ 有限;关键量是高斯宽度。因此可推广为下面的 additive 版本。

Lemma 9.2.4加性 Johnson-Lindenstrauss 引理

令 $\mathcal X\subset\mathbb R^n$ 为有界集,令 $A$ 是 $m\times n$ 随机矩阵,其行 $A_i$ 独立、各向同性且次高斯。则以高概率(例如 $0.99$),缩放后的矩阵 $Q=A/\sqrt m$ 满足

$$ \left| \|Qx-Qy\|_2-\|x-y\|_2 \right| \le \delta \quad\text{for all }x,y\in\mathcal X, $$

其中

$$ \delta=\frac{CK^2w(\mathcal X)}{\sqrt m}, \qquad K=\max_i\|A_i\|_{\psi_2}. $$ 查看学习笔记:Lemma 9.2.4 完整证明
Proof加性 Johnson-Lindenstrauss 引理

对差集

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

应用矩阵偏差不等式(Theorem 9.1.1)。以高概率,

$$ \sup_{x,y\in\mathcal X} \left| \|Ax-Ay\|_2-\sqrt m\|x-y\|_2 \right| \le CK^2\gamma(\mathcal X-\mathcal X). $$

由 Lemma 7.5.11(a),$\gamma(\mathcal X-\mathcal X)=2w(\mathcal X)$。两边除以 $\sqrt m$ 即得结论。

和经典 Johnson-Lindenstrauss 引理对有限集给出的相对误差不同,Lemma 9.2.4 给的是绝对误差 $\delta$。这个差别看起来小,但一般不可避免;Exercise 9.11 会要求你说明原因。

Remark 9.2.5有效维数

为了更好地理解加性 Johnson-Lindenstrauss 引理,可以用数据的有效维数来重写它。令

$$ d(\mathcal X) \asymp \frac{w(\mathcal X)^2}{\operatorname{diam}(\mathcal X)^2}, $$

参见 Definition 7.5.12。若选择

$$ m\gtrsim \varepsilon^{-2}d(\mathcal X) $$

并暂时忽略对 $K$ 的依赖,则 Lemma 9.2.4 中的误差可以做到

$$ \delta=\varepsilon\,\operatorname{diam}(\mathcal X). $$

也就是说,$Q$ 会把所有距离保持到数据直径的一个小比例误差以内;换句话说,它把数据的维数降到了有效维数的量级。

9.3 随机截面:$M^*$ 界与逃逸定理

9.3.1 $M^*$ 界

Theorem 9.3.1$M^*$ 界

令 $E=\ker A$,其中 $A$ 的行独立、各向同性且次高斯。则任意有界 $T\subset\mathbb R^n$ 满足

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

这里有一个有些惊人的高维事实:若用一个余维为 $m$ 的随机子空间 $E$ 去切一个凸集 $T\subset\mathbb R^n$,切片 $T\cap E$ 往往很小,即使 $m\ll n$、$E$ 几乎仍是全维的。下面看它如何从矩阵偏差不等式推出。

建模随机子空间的方便方式是把它写成随机矩阵的核:

$$ E=\ker A. $$

总有 $\dim(E)\ge n-m$;如果 $A$ 有连续分布,则几乎必然 $\dim(E)=n-m$。例如,当 $A$ 是 i.i.d. $N(0,1)$ 高斯矩阵时,由旋转不变性,

$$ E=\ker(A)\sim\operatorname{Unif}(G_{n,n-m}). $$

Proof$M^*$ 界

对差集 $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) \le 2CK^2w(T). $$

若 $x,y\in T\cap E$,则 $z=x-y\in T-T$ 且 $Az=0$。因此

$$ \sqrt m\|x-y\|_2 = \left|\|Az\|_2-\sqrt m\|z\|_2\right|. $$

对 $x,y\in T\cap E$ 取上确界,再取期望,就得到

$$ \mathbb E\operatorname{diam}(T\cap E) \le \frac{CK^2w(T)}{\sqrt m}. $$
Example 9.3.2cross-polytope

把 $M^*$ 界应用于 cross-polytope $B_1^n$,也就是 $\ell^1$ 范数的单位球。由 (7.19),它的高斯宽度约为 $\sqrt{\log n}$,所以

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

例如,如果 $m=0.01n$,则

$$ \mathbb E\operatorname{diam}(B_1^n\cap E) \lesssim \sqrt{\frac{\log n}{n}}. \tag{9.13} $$

也就是说,一个随机的 $0.99n$ 维切片会把 cross-polytope 切得非常小。

直观上,回忆 Figure 7.4a 中对 cross-polytope 的示意:$B_1^n$ 的主体集中在半径 $1/\sqrt n$ 的内切欧氏球附近,而其余部分沿坐标轴伸出细长尖刺。随机子空间通常会避开这些尖刺,只穿过主体,所以切片直径大约是 $O(1/\sqrt n)$,至多多一个如 (9.13) 中的对数因子。这个直觉也适用于一般凸集。

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

为了获得更多直觉,可以用有效维数重写 $M^*$ 界。令

$$ d(T)\asymp \frac{w(T)^2}{\operatorname{diam}(T)^2} $$

参见 Definition 7.5.12。$M^*$ 界说明,只要

$$ m\gtrsim d(T), $$

随机切片就会显著缩小直径,例如

$$ \mathbb E\operatorname{diam}(T\cap E) \le 0.01\,\operatorname{diam}(T) $$

在常数适当时成立。由于 $\dim(E)=n-m$,这个条件等价于

$$ \dim(E)+c\,d(T)\le n. $$

这与线性代数直觉一致:如果 $T$ 是某个子空间 $F\subset\mathbb R^n$ 中以原点为中心的欧氏球,那么只有当 $\dim E+\dim F\le n$ 时,切片才可能真正缩小 $T$ 的直径。

练习 9.12、9.13、9.14 分别要求把 $M^*$ 界推广到仿射截面、推出高概率版本,并计算 $\ell^p$ 球的随机切片直径。

9.3.2 逃逸定理

$M^*$ 界说明随机切片通常很小。另一个相关问题是:随机子空间什么时候会完全避开一个集合?

Theorem 9.3.4逃逸定理

令 $T\subset S^{n-1}$ 为任意集合,令 $A$ 为 $m\times n$ 矩阵,其行 $A_i$ 独立、各向同性且次高斯。若

$$ m\ge CK^4w(T)^2, \tag{9.14} $$

则随机子空间 $E=\ker A$ 满足

$$ T\cap E=\varnothing $$

的概率至少为 $1-2\exp(-cm/K^4)$。这里 $K=\max_i\|A_i\|_{\psi_2}$。

查看学习笔记:Theorem 9.3.4 完整证明
Proof用高概率矩阵偏差排除交点

使用矩阵偏差不等式的高概率版本(见 Remark 9.1.4):以至少 $1-2\exp(-u^2)$ 的概率,

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

假设事件 (9.15) 发生。若 $T\cap E\ne\varnothing$,则对任意 $x\in T\cap E$ 都有 $Ax=0$,于是

$$ \sqrt m \le C_1K^2\bigl(w(T)+u\bigr). $$

$$ u=\frac{\sqrt m}{2C_1K^2}, $$

上式化简为

$$ \sqrt m\le 2C_1K^2w(T). $$

如果 (9.14) 中的绝对常数 $C$ 取得足够大,这与假设矛盾。因此,在这个 $u$ 的选择下,事件 (9.15) 蕴含 $T\cap E=\varnothing$。其概率至少为

$$ 1-2\exp\left(-\frac{m}{4C_1^2K^4}\right) \ge 1-2\exp\left(-\frac{cm}{K^4}\right), $$

这就证明了逃逸定理。

Escape theorem
Figure 9.3:逃逸定理描述随机子空间何时避开球面子集。

9.4 应用:高维线性模型

Tips:本节把几何误差界翻译成估计误差:可行集合 $T$ 表示先验结构,$w(T)$ 表示结构复杂度,Theorem 9.4.4 说明样本数足够时,随机观测不会把结构化信号混淆得太严重。

现在把工具用于一个经典的数据科学问题:在高维中学习线性模型。设有一个未知向量 $x\in\mathbb R^n$,我们想从 $m$ 个线性、可能带噪的观测中学习它:

$$ y_i=\langle A_i,x\rangle+w_i,\qquad i=1,\ldots,m. $$

这里 $A_i\in\mathbb R^n$ 已知,$w_i$ 是未知噪声。用矩阵写就是

$$ y=Ax+w, \tag{9.16} $$

其中 $A$ 是已知的 $m\times n$ 矩阵,$w\in\mathbb R^m$ 是未知噪声。目标是从 $y$ 和 $A$ 尽可能准确地恢复 $x$。

我们假设 $A$ 的行 $A_i$ 随机且独立。这在许多统计场景中是合理的,例如 i.i.d. observations;同时也非常适合使用高维概率工具。

High-dimensional linear model
Figure 9.4:高维线性模型 $y=Ax+w$。
Example 9.4.1音频采样

在信号处理中,$x$ 可以是数字化音频信号,而 $y$ 是在 $m$ 个随机时间点采样得到的结果。

音频采样恢复
Figure 9.5:音频采样中从少量随机采样恢复信号。
Example 9.4.2线性回归

统计学中的核心问题之一是线性回归:我们希望从 $m$ 个样本中学习 $n$ 个预测变量与响应变量之间的线性关系。模型写成

$$ Y=X\theta+w, $$

其中 $X$ 是 $m\times n$ 预测矩阵,$Y\in\mathbb R^m$ 是响应向量,$\theta\in\mathbb R^n$ 是要学习的参数向量,$w$ 是噪声。例如在遗传学中,可以用基因表达数据预测疾病:收集 $m$ 个病人的数据,$X_{ij}$ 表示第 $i$ 个病人第 $j$ 个基因的活跃程度,$Y_i$ 表示该病人是否患病。目标是学习参数向量 $\theta$,它描述每个基因如何影响疾病。

Remark 9.4.3高维情形

现代问题中,数据量经常少于参数量:

$$ m\ll n. $$

例如典型遗传学研究可能只有约 $100$ 个病人,却有约 $10,000$ 个基因。在这种高维设定中,即使是无噪声问题 $y=Ax$ 也无法直接求解,因为解太多,通常形成一个维数至少为 $n-m$ 的大子空间。

不过,如果我们有关于 $x$ 结构的先验信息,仍然可能恢复 $x$。把这种信息写成

$$ x\in T \tag{9.17} $$

其中 $T\subset\mathbb R^n$ 已知。例如,若 $x$ 稀疏,就取 $T$ 为稀疏向量集合;若 $x$ 是带限信号,就取 $T$ 为带限信号的集合。下面看这种结构如何提供帮助。

9.4.1 约束恢复

先考虑 noiseless case:

$$ y=Ax,\qquad x\in T. $$

如何求解这个高维约束线性问题?一个简单想法是:在 $T$ 中找任意一个与观测匹配的向量:

$$ \text{find }x': \quad y=Ax', \quad x'\in T. \tag{9.18} $$

如果 $T$ 是 convex 集合,那么这是一个 convex program,可以用许多数值算法求解。下面检查这个解有多准确。

Theorem 9.4.4约束恢复

假设 $A$ 的行 $A_i$ 独立、各向同性且次高斯。则程序 (9.18) 的任意解 $\hat x$ 满足

$$ \mathbb E\|\hat x-x\|_2 \le \frac{CK^2w(T)}{\sqrt m}, \tag{9.19} $$

其中 $K=\max_i\|A_i\|_{\psi_2}$。

查看学习笔记:Theorem 9.4.4 完整证明
Proof约束恢复来自仿射 $M^*$ 界

由于 $x,\hat x\in T$ 且 $Ax=A\hat x=y$,有

$$ x,\hat x\in T\cap E_x, \qquad E_x=x+\ker A. $$

于是仿射 $M^*$ 界(Exercise 9.12)给出

$$ \mathbb E\|\hat x-x\|_2 \le \mathbb E\operatorname{diam}(T\cap E_x) \le \frac{CK^2w(T)}{\sqrt m}. $$
Constrained high-dimensional linear problem
Figure 9.6:约束恢复问题的几何图像。
Remark 9.4.5有效维数

为了获得直觉,把 (9.19) 用有效维数

$$ d(T)\asymp \frac{w(T)^2}{\operatorname{diam}(T)^2} $$

重写。只要观测数满足

$$ m\gtrsim d(T), $$

就能得到非平凡误差界。由于 $d(T)$ 可能远小于环境维数 $n$,即使在 $m\ll n$ 的高维情形中,恢复也可能成功。

Remark 9.4.6凸松弛

如果 $T$ 非凸,可以把它替换为凸包 $\operatorname{conv}(T)$。这样 (9.18) 变成凸优化问题;而 Theorem 9.4.4 的恢复保证不变,因为 Proposition 7.5.2(c) 给出 $w(\operatorname{conv}(T))=w(T)$。

这个方法很灵活。可以在 Exercise 9.17 中分析均方误差,在 Exercise 9.18 中分析 optimization recovery,在 Exercise 9.19 中处理 noisy case (9.16)。

Remark 9.4.7无约束优化

像 (9.18) 那样强制要求 $y=Ax'$ 且 $x'\in T$ 可能过于僵硬:噪声或者错误的 $T$ 都可能导致无解。更稳妥的做法是放松约束,用惩罚项惩罚违反约束的程度。

给定带噪观测

$$ y=Ax+w, $$

一个常用恢复方式是求解无约束凸优化问题:

$$ \text{minimize } \|y-Ax'\|_2^2+\lambda\|x'\|_T \quad\text{over all }x'\in\mathbb R^n. \tag{9.20} $$

这里 $\|\cdot\|_T$ 是选择的范数,$\lambda\gt 0$ 控制拟合观测 $y$ 与保持解结构化之间的权衡。Exercise 9.20 会证明:若 $A$ 是随机矩阵且 $\lambda$ 选择合适,那么解 $\hat x$ 满足

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

其中 $T$ 是 $\|\cdot\|_T$ 的单位球。也就是说,如果 $x$ 结构好、$\|x\|_T$ 小,并且噪声 $w$ 小,那么从约 $m\asymp d(T)$ 个观测中就能准确恢复。

9.4.2 例:稀疏恢复

有时我们相信 $x$ 是稀疏的,即多数坐标为零或近似为零。例如在遗传学问题中,也许只有约 $10$ 个基因真正影响疾病;在其他应用中,$x$ 也可能在某个基底下稀疏:Example 9.4.1 中的音频信号在时域不稀疏,但它的 Fourier transform 可能集中在少量频率上。

用非零坐标数衡量稀疏性:

$$ \|x\|_0:=|\operatorname{supp}(x)| =|\{i:x_i\ne0\}|. \tag{9.21} $$

若 $\|x\|_0\le s$,称 $x$ 为 $s$-稀疏。这个 $\ell^0$ “范数” 不是真正的范数,但它可看作 $\ell^p$ 范数在 $p\to0$ 时的极限(见 Exercise 9.23)。

快速维数计数说明:若 $A$ 处于一般位置且观测数足够多,$m\ge2\|x\|_0$ 时可以从 $y=Ax$ 恢复 $x$(见 Exercise 9.22)。听起来很好,但若不知道 $x$ 的支撑集,直接搜索所有支撑集计算上不可行:需要检查的支撑集至少有 $\binom ns\ge2^s$ 个。随机观测可以帮助解决这个问题。

如果相信 $x$ 是 $s$-稀疏,自然想取

$$ T=\{x\in\mathbb R^n:\|x\|_0\le s\} $$

作为先验。但 $T$ 高度非凸,求解 (9.18) 仍然很难。一个简单修复是把 $\ell^0$ “范数” 换成最接近它且真正为范数的 $\ell^1$ 范数。由于所有满足 $\|x\|_2\le1$ 的 $s$-稀疏向量都有 $\|x\|_1\le\sqrt s$,取凸集合

$$ T:=\sqrt s\,B_1^n \tag{9.22} $$

作为先验是合理的。程序 (9.18) 变成

$$ \text{find }x': \quad y=Ax', \quad \|x'\|_1\le\sqrt s, \tag{9.23} $$

这是 convex 且可计算的。

Corollary 9.4.8稀疏恢复

假设 $A$ 的行 $A_i$ 独立、各向同性且次高斯。若未知 $s$-稀疏向量 $x\in\mathbb R^n$ 满足 $\|x\|_2\le1$,则程序 (9.23) 的任意解 $\hat x$ 满足

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

其中 $K=\max_i\|A_i\|_{\psi_2}$。

查看学习笔记:Corollary 9.4.8 完整证明
Proof由 $\ell^1$ 球的高斯宽度得到恢复误差

取 $T=\sqrt s\,B_1^n$。结论由 Theorem 9.4.4 和 $\ell^1$ 球的高斯宽度界 (7.19) 得到:

$$ w(T)=\sqrt s\,w(B_1^n)\le C\sqrt{s\log n}. $$
Remark 9.4.9观测数几乎随稀疏度线性增长

Corollary 9.4.8 说明,只要

$$ m\gtrsim s\log n, \tag{9.24} $$

误差就会很小。这意味着当 $s\ll n$ 时,可以用远少于完整维数 $n$ 的观测数有效恢复稀疏向量。Exercise 9.24 会把结论推广到近似稀疏向量。

Remark 9.4.10对数因子改进

单位 $s$-稀疏向量的集合 $S_{n,s}$ 可以更紧地 convexify。与其用 (9.22) 中的 $T=\sqrt s\,B_1^n$,可以用截断 $\ell^1$ 球

$$ T_{n,s}:=\sqrt s\,B_1^n\cap B_2^n. $$

这个 relaxation 很紧。Exercise 9.25 会证明

$$ \operatorname{conv}(S_{n,s}) \subset T_{n,s} \subset 2\operatorname{conv}(S_{n,s}). $$

这个更紧的集合把 (9.24) 中的对数因子改进为

$$ m\gtrsim s\log(en/s). $$

详见 Exercises 9.26-9.27。Exercise 9.28 给出 Garnaev-Gluskin 定理,它刻画 $\ell^1$ 球的随机切片;Exercise 9.21 则介绍 Lasso,它是 (9.20) 在 $\ell^1$ 范数下的特例。

9.4.3 例:低秩恢复

再看一个高维线性问题:从 $m$ 个线性观测恢复一个 $d\times d$ 矩阵 $X$:

$$ y_i=\langle A_i,X\rangle,\qquad i=1,\ldots,m, \tag{9.25} $$

其中 $A_i$ 是已知的独立随机矩阵,内积为 $\langle A,B\rangle=\operatorname{tr}(A^{\mathsf T}B)$。一般需要 $d^2$ 个观测,也就是每个元素一个观测。若想用更少观测,就需要 $X$ 有结构;常见结构是 low rank。

rank 是奇异值的 $\ell^0$ “范数”,高度非凸。与稀疏恢复中的做法类似,把它 convex relax 为奇异值的 $\ell^1$ 范数,也就是 nuclear 范数:

$$ \|X\|_*:=\sum_{i=1}^d s_i(X). $$

若矩阵 $X$ 的 rank 至多 $r$ 且 $\|X\|_F\le1$,则它的奇异值向量至多有 $r$ 个非零坐标且欧氏范数至多 $1$,因此 $\|X\|_*\le\sqrt r$。于是取 $T=\sqrt r\,B_*$,其中

$$ B_*:=\{X\in\mathbb R^{d\times d}:\|X\|_*\le1\} $$

是核范数单位球。恢复程序为

$$ \text{find }X': \quad y_i=\langle A_i,X'\rangle\ \forall i=1,\ldots,m, \quad \|X'\|_*\le\sqrt r. \tag{9.26} $$

Corollary 9.4.11低秩矩阵恢复

设 $A_i$ 是独立高斯随机矩阵,所有元素独立同分布为 $N(0,1)$。若未知 $d\times d$ 矩阵 $X$ 的 rank 至多为 $r$ 且 $\|X\|_F\le1$,则程序 (9.26) 的任意解 $\hat X$ 满足

$$ \mathbb E\|\hat X-X\|_F \le C\sqrt{\frac{rd}{m}}. $$ 查看学习笔记:Corollary 9.4.11 完整证明
Proof核范数球的高斯宽度

由核范数与算子范数的对偶性(Exercise 7.18(a)),对 $d\times d$ 高斯矩阵 $G$,

$$ w(B_*) = \mathbb E\sup_{\|X\|_*\le1}\langle G,X\rangle = \mathbb E\|G\| \le 2\sqrt d, $$

最后一步来自 Theorem 7.3.1。对 $T=\sqrt r\,B_*$ 使用 Theorem 9.4.4 即得结论。

Remark 9.4.12少量观测下恢复低秩矩阵

Corollary 9.4.11 在

$$ m\gtrsim rd $$

时给出小误差。若 $r\ll d$,这远少于元素数量 $d^2$。这和第 6.5 节的矩阵补全类似,后者用大约 $rd\log d$ 个随机观测恢复低秩矩阵。Exercise 9.29 会把低秩恢复推广到矩形矩阵和近似低秩矩阵。

9.5 应用:精确稀疏恢复

Tips:精确稀疏恢复有两条证明线:逃逸定理从几何上排除非零误差,RIP 从所有稀疏方向的近似等距性出发。两条线结论相近,但训练的工具不同。

在 noiseless case 中,还可以做得更强:从 $y=Ax$ 精确恢复稀疏向量 $x$,而且算法有效。我们看两条路线:

  1. 使用逃逸定理(Theorem 9.3.4)。
  2. 找到一个保证精确恢复的确定性条件,也就是限制等距性质,再证明随机矩阵以高概率满足它。

9.5.1 基于逃逸定理的精确恢复

先建立几何直觉。设要通过凸优化程序 (9.23) 从 $y=Ax$ 恢复未知 $s$-稀疏单位向量 $x$。解 $\hat x$ 位于先验集合 $T=\sqrt s\,B_1^n$ 与仿射子空间 $E_x=x+\ker A$ 的交中。

$T$ 是 cross-polytope,而 $x$ 位于它的某个 $(s-1)$ 维边上。若随机子空间 $E_x$ 在 $x$ 处与多胞体相切,那么 $T$ 与 $E_x$ 只在 $x$ 相交,从而解必然精确:

$$ \hat x=x. $$

要证明这个论证,需要说明随机子空间 $E_x$ 以高概率与 $\ell^1$ 球在 $x$ 处相切。放大看 $x$ 附近:这等价于切锥的球面部分与 $E_x$ 不相交,而逃逸定理正是用来保证这种不相交的工具。

现在形式化。我们希望从

$$ y=Ax $$

通过优化问题

$$ \text{minimize }\|x'\|_1 \quad\text{subject to }y=Ax' \tag{9.27} $$

恢复 $x$。

Theorem 9.5.1精确稀疏恢复

设 $A$ 是 $m\times n$ 随机矩阵,其行 $A_i$ 独立、各向同性且次高斯。若

$$ m\ge CK^4s\log n, $$

则以至少 $1-2\exp(-cm/K^4)$ 的概率,对任意 $s$-稀疏向量 $x\in\mathbb R^n$,程序 (9.27) 的解是精确的:

$$ \hat x=x. $$

这里 $K=\max_i\|A_i\|_{\psi_2}$。

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

证明要说明 recovery error 为零:

$$ h:=\hat x-x=0. $$

先证明一个弱结论:误差 $h$ 在 $x$ 的支撑集上更重。

Exact sparse recovery
Figure 9.7:精确稀疏恢复的切锥几何。
Lemma 9.5.2$x$ 的支撑集上误差更重

令 $S:=\operatorname{supp}(x)$,并令 $h_S$ 表示 $h$ 限制在 $S$ 上的向量,$h_{S^c}$ 类似。则

$$ \|h_{S^c}\|_1\le\|h_S\|_1. $$ Lemma 9.5.2
Proof利用 $\ell^1$ 最小性

因为 $\hat x$ 是 (9.27) 的 minimizer,

$$ \|\hat x\|_1\le\|x\|_1. \tag{9.28} $$

另一方面,由三角不等式且 $x_S=x$、$x_{S^c}=0$,有

$$ \|\hat x\|_1 = \|x+h\|_1 = \|x_S+h_S\|_1+\|x_{S^c}+h_{S^c}\|_1 \ge \|x\|_1-\|h_S\|_1+\|h_{S^c}\|_1. $$

与 (9.28) 合并并化简,即得结论。

Lemma 9.5.3误差近似稀疏

误差向量满足

$$ \|h\|_1\le2\sqrt s\,\|h\|_2. $$ Lemma 9.5.3
Proof由支撑集控制到近似稀疏

由 Lemma 9.5.2 和 Cauchy-Schwarz inequality,

$$ \|h\|_1 = \|h_S\|_1+\|h_{S^c}\|_1 \le 2\|h_S\|_1 \le 2\sqrt s\,\|h_S\|_2 \le 2\sqrt s\,\|h\|_2. $$
Proof of Theorem 9.5.1逃逸定理排除非零误差

假设 $h=\hat x-x\ne0$。Lemma 9.5.3 给出

$$ \frac{h}{\|h\|_2} \in T_s:= \left\{ z\in S^{n-1}:\|z\|_1\le2\sqrt s \right\}. $$

又因为 $Ah=A\hat x-Ax=y-y=0$,所以

$$ \frac{h}{\|h\|_2} \in T_s\cap\ker A. \tag{9.29} $$

逃逸定理(Theorem 9.3.4)说明,只要 $m\ge C_1K^4w(T_s)^2$,该交集就以高概率为空。另一方面,由 $T_s\subset2\sqrt s\,B_1^n$ 和 (7.19),

$$ w(T_s) \le 2\sqrt s\,w(B_1^n) \le C_2\sqrt{s\log n}. \tag{9.30} $$

因此若 $m\ge CK^4s\log n$,则 (9.29) 中的交集以高概率为空,与 $h\ne0$ 矛盾。所以 $h=0$,证明完成。

Remark 9.5.4改进对数因子

稍微改进 (9.30),可把 Theorem 9.5.1 中充分的观测数改进为

$$ m\ge CK^4s\log(en/s). $$

这来自 Exercise 9.26。Exercise 9.30 会让你几何化理解精确恢复证明,Exercise 9.31 处理带噪声情形,Exercise 9.32 则给出精确恢复的零空间性质。

9.5.2 受限等距性

下面寻找一个保证稀疏恢复的确定性条件,并证明随机矩阵以高概率满足它。这个条件称为限制等距性质。

Definition 9.5.5RIP

一个 $m\times n$ 矩阵 $A$ 满足参数为 $\alpha,\beta,s$ 的限制等距性质 (RIP),若对所有至多有 $s$ 个非零坐标的 $v\in\mathbb R^n$,都有

$$ \alpha\|v\|_2\le\|Av\|_2\le\beta\|v\|_2. $$

RIP 等价于说 $A$ 的所有 $m\times s$ 子矩阵 $A_I$ 的奇异值满足

$$ \alpha\le s_s(A_I)\le s_1(A_I)\le\beta. \tag{9.31} $$

若 $\alpha\approx\beta\approx1$,则所有这些子矩阵都是近似等距映射。

Theorem 9.5.6RIP 推出精确恢复

假设 $m\times n$ 矩阵 $A$ 满足参数 $\alpha,\beta,(1+\lambda)s$ 的 RIP,其中 $\lambda\gt (\beta/\alpha)^2$。那么每个 $s$-稀疏向量 $x\in\mathbb R^n$ 都可由 (9.27) 从 $y=Ax$ 精确恢复,即 $\hat x=x$。

查看学习笔记:Theorem 9.5.6 完整证明
ProofRIP 排除非零误差

和 Theorem 9.5.1 的证明一样,需要证明误差 $h=\hat x-x$ 为零。

步骤 1:分解支撑集。 令 $I_0=\operatorname{supp}(x)$。令 $I_1$ 为 $h_{I_0^c}$ 中绝对值最大的 $\lambda s$ 个元素的指标集合,$I_2$ 为接下来 $\lambda s$ 个,依此类推。记 $I_{01}=I_0\cup I_1$。由于 $Ah=A\hat x-Ax=0$,三角不等式给出

$$ 0=\|Ah\|_2 \ge \|A_{I_{01}}h_{I_{01}}\|_2 - \|A_{I_{01}^c}h_{I_{01}^c}\|_2. \tag{9.32} $$

步骤 2:使用 RIP。 因为 $|I_{01}|\le(1+\lambda)s$,RIP 给出

$$ \|A_{I_{01}}h_{I_{01}}\|_2 \ge \alpha\|h_{I_{01}}\|_2. $$

同时由三角不等式和 RIP,

$$ \|A_{I_{01}^c}h_{I_{01}^c}\|_2 \le \sum_{i\ge2}\|A_{I_i}h_{I_i}\|_2 \le \beta\sum_{i\ge2}\|h_{I_i}\|_2. $$

代入 (9.32),得到

$$ \beta\sum_{i\ge2}\|h_{I_i}\|_2 \ge \alpha\|h_{I_{01}}\|_2. \tag{9.33} $$

步骤 3:求和估计。 按照 $I_i$ 的定义,$h_{I_i}$ 的每个元素绝对值都不超过前一块平均绝对值,即对 $i\ge2$,

$$ \|h_{I_i}\|_2 \le \frac1{\sqrt{\lambda s}}\|h_{I_{i-1}}\|_1. $$

求和并使用 Lemma 9.5.2,得到

$$ \sum_{i\ge2}\|h_{I_i}\|_2 \le \frac1{\sqrt{\lambda s}}\|h_{I_0^c}\|_1 \le \frac1{\sqrt{\lambda s}}\|h_{I_0}\|_1 \le \frac1{\sqrt\lambda}\|h_{I_0}\|_2 \le \frac1{\sqrt\lambda}\|h_{I_{01}}\|_2. $$

代入 (9.33):

$$ \frac{\beta}{\sqrt\lambda}\|h_{I_{01}}\|_2 \ge \alpha\|h_{I_{01}}\|_2. $$

由于 $\beta/\sqrt\lambda\lt \alpha$,只能有 $h_{I_{01}}=0$。而 $I_{01}$ 包含 $h$ 的最大元素,因此 $h=0$。

虽然我们不知道如何构造具有好参数的 deterministic RIP 矩阵,但可以证明随机矩阵以高概率满足 RIP。

Theorem 9.5.7随机矩阵满足 RIP

设 $A$ 是 $m\times n$ 随机矩阵,其行 $A_i$ 独立、各向同性且次高斯。若

$$ m\ge CK^4s\log(en/s), $$

其中 $K=\max_i\|A_i\|_{\psi_2}$,则以至少 $1-2\exp(-cm/K^4)$ 的概率,$A$ 满足参数 $\alpha=0.9\sqrt m$、$\beta=1.1\sqrt m$ 和 $s$ 的 RIP。

查看学习笔记:Theorem 9.5.7 完整证明
Proof固定支撑集后并集界

需要对所有 $m\times s$ 子矩阵 $A_I$ 检查 (9.31)。先固定 $I$。由 Theorem 4.6.1,

$$ 0.9\sqrt m \le s_s(A_I) \le s_1(A_I) \le 1.1\sqrt m \tag{9.34} $$

以至少 $1-2\exp(-2cm/K^4)$ 的概率成立;这里是在 Theorem 4.6.1 中取 $t=\sqrt{2cm}/K^2$,并用 $m$ 的假设适当选择 absolute constants $c,C$。

现在对所有 $\binom ns$ 个 $s$ 元子集 $I\subset\{1,\ldots,n\}$ 做并集界。于是 (9.34) 同时对所有 $I$ 成立的概率至少为

$$ 1 - 2\exp\left(-\frac{2cm}{K^4}\right)\binom ns. $$

再由

$$ \binom ns\le\exp(s\log(en/s)) $$

以及 $m\ge CK^4s\log(en/s)$,上面的概率下界进一步变成

$$ 1-2\exp\left(-\frac{cm}{K^4}\right). $$

因此所有 $m\times s$ 子矩阵同时满足 (9.34),也就是 $A$ 满足所需 RIP。

Theorem 9.5.1 的第二种证明由 RIP 推出精确稀疏恢复

由 Theorem 9.5.7,随机矩阵 $A$ 以高概率对 $3s$ 阶满足 RIP,且参数可取

$$ \alpha=0.9\sqrt m,\qquad \beta=1.1\sqrt m. $$

在 Theorem 9.5.6 中取 $\lambda=2$,因为此时 $(1+\lambda)s=3s$,并且在常数选择下满足 $\lambda\gt (\beta/\alpha)^2$。于是 Theorem 9.5.6 保证 (9.27) 的解精确恢复 $x$。

所以 Theorem 9.5.1 成立;事实上,这一路线还给出 Remark 9.5.4 中的对数因子改进。

Exercise 9.33 会让你证明随机投影满足 RIP。

9.6 一般范数下的随机矩阵偏差

下面把矩阵偏差不等式推广到任意范数,而不仅是欧氏范数。事实上,甚至不要求函数非负;正齐次性和三角不等式型的次可加性已经足够。

Definition 9.6.1正齐次与次可加函数

线性空间 $V$ 上的实值函数 $f$ 称为正齐次,若对所有 $\alpha\ge0$ 和 $x\in V$,

$$ f(\alpha x)=\alpha f(x). $$

称为次可加,若对所有 $x,y\in V$,

$$ f(x+y)\le f(x)+f(y). $$
Example 9.6.2正齐次与次可加函数

以下函数都是正齐次且次可加:

(a) 任意范数;

(b) 任意实值 linear 函数,也叫 linear 泛函;

(c) 特别地,对固定 $y\in\mathbb R^m$,函数 $f(x)=\langle x,y\rangle$;

(d) 任意有界集合 $S\subset\mathbb R^m$ 的支撑函数

$$ f(x):=\sup_{y\in S}\langle x,y\rangle,\qquad x\in\mathbb R^m. \tag{9.35} $$

下面的 Theorem 9.1.1 版本适用于所有范数,甚至适用于正齐次、次可加函数;代价是这里只处理高斯矩阵。

Theorem 9.6.3一般矩阵偏差不等式

令 $A$ 为 $m\times n$ 随机矩阵,其元素独立同分布为 $N(0,1)$。设 $f:\mathbb R^m\to\mathbb R$ 是有界、正齐次且次可加的函数,并且存在 $b\in\mathbb R$ 使得

$$ f(x)\le b\|x\|_2,\qquad x\in\mathbb R^m. \tag{9.36} $$

则对任意 $T\subset\mathbb R^n$,

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

和第 9.1 节一样,只要证明随机过程

$$ Z_x:=f(Ax)-\mathbb Ef(Ax) \tag{9.37} $$

具有次高斯增量,就可以由 Talagrand 比较推出 Theorem 9.6.3。

Theorem 9.6.4一般范数下的次高斯增量

令 $A$ 为 $m\times n$ 高斯随机矩阵,元素独立同分布为 $N(0,1)$。设 $f:\mathbb R^m\to\mathbb R$ 正齐次、次可加且满足 (9.36)。则过程 (9.37) 满足

$$ \|Z_x-Z_y\|_{\psi_2} \le Cb\|x-y\|_2,\qquad x,y\in\mathbb R^n. \tag{9.38} $$ 查看学习笔记:Theorem 9.6.4 完整证明
Proof of Theorem 9.6.4用高斯集中制造增量控制

不妨设 $b=1$。先假设 $\|x\|_2=\|y\|_2=1$,此时要证

$$ \|f(Ax)-f(Ay)\|_{\psi_2} \le C\|x-y\|_2. \tag{9.39} $$

步骤 1:制造独立性。

$$ u:=\frac{x+y}{2},\qquad v:=\frac{x-y}{2}. \tag{9.40} $$

于是 $x=u+v$、$y=u-v$,从而

$$ Ax=Au+Av,\qquad Ay=Au-Av. $$

由于 $u$ 与 $v$ 正交,高斯随机向量 $Au$ 和 $Av$ 独立。

步骤 2:使用高斯集中。 条件在 $a:=Au$ 上,考察

$$ f(Ax)=f(a+Av). $$

由独立性,$a+Av$ 的条件分布可写成

$$ a+Av=a+\|v\|_2g, \qquad g\sim N(0,I_m). $$

函数 $z\mapsto f(a+\|v\|_2z)$ 的 Lipschitz 范数至多为 $\|v\|_2$。事实上,对任意 $t,s\in\mathbb R^m$,由次可加性、正齐次性和 (9.36),

$$ \begin{aligned} f(a+\|v\|_2t)-f(a+\|v\|_2s) &\le f(\|v\|_2(t-s))\\ &=\|v\|_2f(t-s)\\ &\le\|v\|_2\|t-s\|_2. \end{aligned} $$

高斯集中(Theorem 5.2.3)给出条件尾界

$$ \left\| f(a+Av)-\mathbb E_a f(a+Av) \right\|_{\psi_2(a)} \le C\|v\|_2. \tag{9.41} $$

步骤 3:去掉条件。 随机向量 $a-Av$ 与 $a+Av$ 同分布,因此也满足

$$ \left\| f(a-Av)-\mathbb E_a f(a-Av) \right\|_{\psi_2(a)} \le C\|v\|_2. \tag{9.42} $$

将 (9.42) 从 (9.41) 中相减,使用三角不等式,并注意两个条件期望相同,得到

$$ \|f(a+Av)-f(a-Av)\|_{\psi_2(a)} \le 2C\|v\|_2. $$

该条件界对每个固定 $a$ 都成立,因此无条件也成立。回到 $x,y$ 记号,并用 $\|v\|_2=\|x-y\|_2/2$,得到 (9.39)。一般 $x,y$ 的情形与 Theorem 9.1.2 证明中的步骤 4 相同。

Orthogonal u and v from x and y
Figure 9.8:由 $x,y$ 构造正交向量 $u=(x+y)/2$ 和 $v=(x-y)/2$。
Remark 9.6.5开放问题

Theorem 9.6.3 是否对一般次高斯矩阵 $A$ 成立,是一个开放问题。Exercises 9.35 和 9.36 分别给出各向异性版本和高概率版本。

9.7 双侧 Chevet 不等式与 Dvoretzky-Milman 定理

Tips:本节是“理论回流到几何”的部分:一般范数偏差给出双侧 Chevet,不只控制上界,也控制下界;Dvoretzky-Milman 再说明随机低维投影会把任意高维凸体圆化。

和本章原始的矩阵偏差不等式一样,更一般的 Theorem 9.6.3 有很多应用。比如,你现在可以在任意范数中得到 Johnson-Lindenstrauss 型引理,而不仅是欧氏范数;见 Exercises 9.37-9.39。

9.7.1 双侧 Chevet 不等式

general 矩阵偏差的另一个结果是 Chevet 不等式的更尖锐版本。我们在第 8.6 节已经见过 Chevet 不等式。

Theorem 9.7.1双侧 Chevet 不等式

令 $A$ 为 $m\times n$ 高斯随机矩阵,元素独立同分布为 $N(0,1)$。令 $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 完整证明

由三角不等式可见,Theorem 9.7.1 给出 Chevet 不等式(Theorem 8.6.1)的更强双侧版本。

Proof支撑函数代入一般偏差不等式

对 $S$ 的支撑函数使用 Theorem 9.6.3:

$$ f(x)=\sup_{y\in S}\langle x,y\rangle. $$

由 Cauchy-Schwarz inequality,

$$ f(x) \le \sup_{y\in S}\|x\|_2\|y\|_2 = \operatorname{rad}(S)\|x\|_2. \tag{9.43} $$

另一方面,$Ax$ 与 $g\|x\|_2$ 同分布,其中 $g\sim N(0,I_m)$。由正齐次性,

$$ \begin{aligned} \mathbb Ef(Ax) &= \|x\|_2\mathbb Ef(g)\\ &= \|x\|_2\mathbb E\sup_{y\in S}\langle g,y\rangle\\ &= \|x\|_2w(S). \tag{9.44} \end{aligned} $$

把 (9.43) 与 (9.44) 代入 Theorem 9.6.3,即得结论。

9.7.2 Dvoretzky-Milman 定理

现在证明一个令人惊讶的结果:把任意有界集合随机投影到低维空间后,它会以高概率看起来近似圆。

Theorem 9.7.2Dvoretzky-Milman 定理

令 $A$ 为 $m\times n$ 高斯随机矩阵,元素独立同分布为 $N(0,1)$,令 $T\subset\mathbb R^n$ 有界。则以至少 $0.99$ 的概率,

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

其中 $B_2^m$ 是 $\mathbb R^m$ 中的欧氏单位球,

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

左侧包含关系只在 $r_-\ge0$ 时有意义;右侧包含关系总成立。

查看学习笔记:Theorem 9.7.2 完整证明
Dvoretzky-Milman random projection
Figure 9.9:8 维 cube 与 $10^4$ 个高斯 points 投影到平面的近圆形效果。
Proof由双侧 Chevet 到支撑函数夹逼

先把 two-sided Chevet 不等式(Theorem 9.7.1)写成如下形式:

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

这是把 Theorem 9.7.1 应用于 $A^{\mathsf T}$ 并交换 $T,S$ 得到的。

令 $S=S^{m-1}$。由于 $\gamma(S)\le\sqrt m$,Markov inequality 给出:以至少 $0.99$ 的概率,对所有 $y\in S^{m-1}$,

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

由 $r_\pm$ 的定义,这意味着

$$ r_- \le \sup_{x\in T}\langle Ax,y\rangle \le r_+, \qquad y\in S^{m-1}. $$

把 $\sup_{x\in T}\langle Ax,y\rangle$ 写成 $\sup_{z\in AT}\langle z,y\rangle$,并利用齐次性,得到对所有 $y\in\mathbb R^m$,

$$ r_-\|y\|_2 \le \sup_{z\in AT}\langle z,y\rangle \le r_+\|y\|_2. $$

按支撑函数的对偶性,这正等价于 (9.45)。Exercise 9.40 会让你写出这个对偶性细节。

Remark 9.7.3有效维数

设 $T$ 有界、convex 且包含原点,且

$$ m\le cd(T), $$

其中

$$ d(T)\asymp \frac{w(T)^2}{\operatorname{rad}(T)^2}. $$

若常数 $c$ 足够小,则 $C\sqrt m\,\operatorname{rad}(T)\le0.01w(T)$,因此 Theorem 9.7.2 给出

$$ 0.99B\subset AT\subset1.01B, $$

其中 $B=w(T)B_2^m$。也就是说,把有界凸集 $T$ 随机投影到维数约为 $d(T)$ 的子空间后,它看起来几乎像圆球。

Example 9.7.4立方体的近圆投影

考虑 cube $T=[-1,1]^n$。由 (7.18),$w(T)=\sqrt{2/\pi}\,n$,而 $\operatorname{diam}(T)=2\sqrt n$,所以有效维数满足 $d(T)\asymp n$。因此若 $m\le cn$,则以高概率

$$ 0.99B \subset A[-1,1]^n \subset 1.01B, $$

其中 $B$ 是半径为 $\sqrt{2/\pi}\,n$ 的欧氏球。也就是说,把 $n$ 维立方体投影到 $m=cn$ 维子空间后,它看起来几乎像圆球。Figure 9.9 展示了这一现象。

Remark 9.7.5随机投影总结

第 7.6 节和第 9.2.2 节说明:集合 $T$ 到 $\mathbb R^n$ 中 $m$ 维随机子空间的随机投影 $P$ 会经历相变。

在高维情形,也就是 $m\gtrsim d(T)$ 时,投影会按 $\sqrt{m/n}$ 的量级缩小直径:

$$ \operatorname{diam}(PT) \asymp \sqrt{\frac mn}\operatorname{diam}(T). $$

此外,additive Johnson-Lindenstrauss Lemma 9.2.4 说明,在这个 regime 中,随机投影近似保持 $T$ 的几何,即所有点间距离大约按同一比例缩小。

在 low-dimensional regime,也就是 $m\lesssim d(T)$ 时,缩小停止:

$$ \operatorname{diam}(PT) \asymp w_s(T) \asymp \frac{w(T)}{\sqrt n}, $$

无论 $m$ 多小都如此。Dvoretzky-Milman 定理解释了原因:此时 $PT$ 近似为一个半径约为 $w_s(T)$ 的圆球,而圆球在进一步投影下不会继续按原集合直径缩小。Exercises 9.41-9.43 会分别处理 Dvoretzky-Milman 的高概率版本、高斯点云的低维近圆性,以及真正随机投影版本。

9.8 Notes

矩阵偏差不等式(Theorem 9.1.1)及其证明来自 [213],不过很多相关结果更早就已经存在。对高斯矩阵 $A$ 以及单位球面上的集合 $T$,它可以从高斯比较不等式推出:上界来自 Sudakov-Fernique(Theorem 7.2.8),下界来自 Gordon(Theorem 7.2.9)。Schechtman [295] 证明了高斯 $A$ 与一般范数情形下的矩阵偏差不等式;我们在第 9.6 节中介绍了这一版本。对于次高斯 $A$,相关版本见 [187, 236, 103];可参见 [213, Section 3] 中的比较。关于二次型过程,一个特别干净的结果见 [79, Theorem 3.2.1]。其他推广还包括稀疏 $A$ 的版本 [54]、$\ell^p$ 范数的版本 [301],以及独立列的版本 [277]。

Theorem 9.1.1 中关于 $K$ 的二次依赖后来在 [176] 中改进为最优的 $K\sqrt{\log K}$。这会自动改进所有从 Theorem 9.1.1 推出的结果里对次高斯范数的依赖,例如 Theorems 3.1.1、4.6.1、5.3.1、9.3.1、9.3.4、9.4.4,Proposition 9.2.1,Corollary 9.4.8,等等。

关于集合随机投影的界(Proposition 9.2.1)的一个版本可追溯到 V. Milman [245];见 [21, Proposition 5.7.1]。

低维分布协方差估计的 Theorem 9.2.2 归功于 V. Koltchinskii 与 K. Lounici [190];他们使用了另一种同样基于主控测度定理的方法。R. van Handel [331] 对高斯分布给出了另一个证明,使用解耦、条件化和 Slepian 引理。对高斯矩阵,Theorem 9.2.2 的界是紧的;现在已有许多推广,参考文献见第 4 章末尾。

加性 Johnson-Lindenstrauss 引理(Lemma 9.2.4)的一个版本来自 [213]。

$M^*$ 界(Theorem 9.3.1)在几何泛函分析中已有较长历史。早期版本来自 V. Milman [243, 244];Pajor 和 Tomczak [269] 证明了一个对 $m$ 依赖正确的版本,Gordon 后来给出了具有精确常数的更尖锐形式 [143]。更多内容,包括证明与变体,可见 [21, Sections 7.3-7.4, 9.3]、[143, 233, 341]。本书 Theorem 9.3.1 的版本来自 [213]。

逃逸定理(Theorem 9.3.4)也称为 “网格逃逸”。它最早由 Y. Gordon [143] 对高斯矩阵证明,并利用他的比较不等式(Theorem 7.2.9)在 (9.14) 中得到尖锐常数。对球面凸集,已有匹配的下界 [308, 18];在这种情形中,甚至可以用积分几何的工具精确计算击中概率 [18]。Oymak 和 Tropp [267] 说明了如何把该结果推广到高斯情形之外。本书 Theorem 9.3.4 的版本来自 [213]。

第 9.4-9.5 节中的应用来自高维统计学与信号处理,特别是压缩感知。教程 [341] 对这两类问题给出了统一处理,本章遵循了其中的思路。书籍 [67, 157, 344, 137] 讨论统计方面,综述 [91] 与书籍 [127] 侧重信号处理方面。

第 9.4.1 节中基于 $M^*$ 界的恢复方法源于 [341],其中包含 Theorem 9.4.4 与 Corollary 9.4.8 的多种版本。

综述 [93] 对第 9.4.3 节中讨论的低秩矩阵恢复问题作了全面综述。我们的表述基于 [341, Section 10]。

第 9.5 节中的精确稀疏恢复发现于压缩感知领域;见 [91] 以及书籍 [127]。通过逃逸定理的方法(第 9.5.1 节)最早发现于 [290];这里我们大致遵循 [341, Section 9]。另见 [80, 307],尤其是 [322],其中有逃逸定理在稀疏恢复中的应用。关于精确恢复所需观测数的尖锐界最早由 [112] 证明,随后在 [111, 108, 109, 110] 中推广。更一般的集合 $T$ 和矩阵 $A$ 的相变见 [18, 266, 267]。

基于 RIP 的稀疏恢复方法(第 9.5.2 节)由 E. Candes 和 T. Tao [73] 开创;全面介绍见 [127, Chapter 6]。Theorem 9.5.6 的一个版本出现在他们的工作中,而我们的证明基于 Y. Plan 的一个论证,类似于 [70]。随机矩阵满足 RIP(Theorem 9.5.7)是压缩感知的基石之一;见 [127, 340]。

关于二次型过程的偏差(Exercise 9.5),见 [79, Theorem 3.2.1],那是一个非常相近的结果。

Exercises 9.19-9.21 讨论稀疏线性回归中的一个常用工具,即统计文献中的 Lasso(least absolute shrinkage and selection operator)。它由 R. Tibshirani [319] 开创。书籍 [157, 67, 344] 对带稀疏约束的统计问题作了全面介绍;这些书也讨论 Lasso 及其许多变体。

基于零空间性质的稀疏恢复(Exercise 9.32)可追溯到 Cohen、Dahmen 和 DeVore [88],见 [281, 127, 344]。

Garnaev-Gluskin 界(Exercise 9.28)最早由 [132] 证明,另见 [222] 与 [127, Chapter 10]。

一般矩阵偏差不等式(Theorem 9.6.3)及其证明归功于 G. Schechtman [295]。

Chevet 不等式的原始版本由 S. Chevet [84] 证明,其中的常数因子由 Y. Gordon [140] 改进;另见 [21, Section 9.4]、[210, Theorem 3.20] 与 [321, 7]。我们在 Theorem 9.7.1 中陈述的 Chevet 不等式版本可以从 Y. Gordon [140, 142] 的工作中重构出来;见 [210, Corollary 3.21]。

Dvoretzky-Milman 定理在泛函分析中有很长历史。A. Dvoretzky [117, 118] 证明了 A. Grothendieck 的一个猜想:任意 $n$ 维赋范空间都有一个 $m$ 维近似欧氏子空间,其中 $m=m(n)$ 随 $n$ 增大而趋于无穷。V. Milman 给出了该定理的概率证明 [242],并开创了对最佳依赖 $m(n)$ 的研究。Theorem 9.7.2 归功于 V. Milman [242];它是最优的 [247],见 [21, Theorem 5.3.3]。教程 [23] 对 Dvoretzky-Milman 定理作了轻量介绍。关于 Dvoretzky-Milman 定理及其许多影响的完整阐述,可见 [21, Chapter 5 and Section 9.2]、[210, Section 9.1] 及其中参考文献。Dvoretzky-Milman 定理还有一种“分布式”版本,问题是 $n$ 维分布的随机 $m$ 维边缘分布是否近似正态;这一路线由 B. Klartag [186] 在对数凹分布中开创(凸体的中心极限定理),并由 E. Meckes [230] 在离散分布中开创。

Remark 9.7.5 中提到的相变现象由 V. Milman [245] 提出;见 [21, Proposition 5.7.1]。

Exercises

Tips:第 9 章习题建议按应用线索读:9.1-9.6 加深矩阵偏差本身,9.12-9.17 训练 $M^*$/逃逸定理,9.19-9.33 进入稀疏和低秩恢复,9.37-9.43 连接一般范数 JL 与 Dvoretzky-Milman。
Exercise 9.1反向三角不等式

验证 Theorem 9.1.2 证明中用到的几何观察。设 $x,y\in\mathbb R^n$ 满足 $1=\|x\|_2\le\|y\|_2$,并令 $\bar y=y/\|y\|_2$。证明

$$ \|x-y\|_2 \le \|x-\bar y\|_2+\|\bar y-y\|_2 \le \sqrt2\,\|x-y\|_2. $$ 查看学习笔记:Exercise 9.1 证明
Exercise 9.2矩阵偏差不等式:均值偏差

从 Theorem 9.1.1 推出矩阵偏差的中心化版本:

$$ \mathbb E\sup_{x\in T} \left| \|Ax\|_2-\mathbb E\|Ax\|_2 \right| \le CK^2\gamma(T). $$ 查看学习笔记:Exercise 9.2 证明
Exercise 9.3二次型矩阵偏差

证明 Remark 9.1.5 中的结论:在 Theorem 9.1.1 的条件下,

$$ \mathbb E\sup_{x\in T} \left| \|Ax\|_2^2-m\|x\|_2^2 \right| \le CK^4\gamma(T)^2 + CK^2\sqrt m\,\operatorname{rad}(T)\gamma(T). $$ 查看学习笔记:Exercise 9.3 证明
Exercise 9.4各向异性矩阵偏差

Theorem 9.1.1 假设随机矩阵的行是各向同性。现在去掉这个假设。令 $B$ 为 $m\times n$ 随机矩阵,其独立行 $B_i$ 满足

$$ \mathbb E B_iB_i^\mathsf T=\Sigma,\qquad \|\langle B_i,x\rangle\|_{\psi_2} \le K\|\langle B_i,x\rangle\|_{L^2} \quad\text{for any }x\in\mathbb R^n. $$

证明对任意 $T\subset\mathbb R^n$,有

$$ \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.4 证明
Exercise 9.5二次型经验过程

Exercise 8.36 曾对次高斯过程证明经验均值偏离样本均值的界。现在通过扩展 Theorem 9.1.1,为 $L^2$ 均值做同类估计。

设 $\mathcal F$ 是定义在某个定义域 $\Omega$ 上的一族实值函数,并且是 star-shaped:若 $f\in\mathcal F$,则任意 $a\in[0,1]$ 都有 $af\in\mathcal F$。令 $X,X_1,\ldots,X_m$ 是 $\Omega$ 中的 i.i.d. 随机点。考虑 $\mathcal F$ 上的 $L^2$ metric

$$ d(f,g):=\bigl(\mathbb E(f(X)-g(X))^2\bigr)^{1/2}, $$

并假设

$$ \|f(X)-g(X)\|_{\psi_2}\le Kd(f,g) \quad\text{for all }f,g\in\mathcal F. $$

证明

$$ \mathbb E\sup_{f\in\mathcal F} \left| \left(\frac1m\sum_{i=1}^m f(X_i)^2\right)^{1/2} - \left(\mathbb Ef(X)^2\right)^{1/2} \right| \le \frac{CK^2\gamma_2(\mathcal F,d)}{\sqrt m}, $$

右侧的 $\gamma_2$ 是第 8.5 节中的泛函。说明当 $T\subset S^{n-1}$ 时,这个结果如何包含 Theorem 9.1.1。

查看学习笔记:Exercise 9.5 证明
Exercise 9.6随机投影偏差

证明 Theorem 9.1.1 对随机投影的版本。令 $P$ 是 $\mathbb R^n$ 到一个 $m$ 维随机子空间上的正交投影,该子空间在 Grassmannian $G_{n,m}$ 上均匀分布。证明对任意 $T\subset\mathbb R^n$,

$$ \mathbb E\sup_{x\in T} \left| \|Px\|_2-\sqrt{\frac mn}\|x\|_2 \right| \le \frac{C\gamma(T)}{\sqrt n}. $$ 查看学习笔记:Exercise 9.6 证明
Exercise 9.7投影大小:高概率界

证明 Proposition 9.2.1 的高概率版本。令 $T\subset\mathbb R^n$ 为有界集,令 $A$ 是 $m\times n$ 随机矩阵,其行 $A_i$ 独立、各向同性且次高斯,并考虑次高斯投影 $P=A/\sqrt n$。证明对 $\varepsilon\gt 0$,

$$ \operatorname{diam}(PT) \le (1+\varepsilon)\sqrt{\frac mn}\operatorname{diam}(T) + CK^2w_s(T) $$

以概率至少 $1-\exp(-c\varepsilon^2m/K^4)$ 成立。这里 $K=\max_i\|A_i\|_{\psi_2}$,$w_s(T)$ 是 $T$ 的球面宽度。

查看学习笔记:Exercise 9.7 证明
Exercise 9.8投影到随机子空间

为第 7.6 节中的原始模型证明 Proposition 9.2.1 的版本:令 $P$ 是到随机 $m$ 维子空间 $E\sim\operatorname{Unif}(G_{n,m})$ 上的正交投影。

查看学习笔记:Exercise 9.8 证明
Exercise 9.9高概率协方差估计

验证 Remark 9.2.3 中提到的协方差估计高概率保证。

查看学习笔记:Exercise 9.9 证明
Exercise 9.10从矩阵偏差到 Johnson-Lindenstrauss

使用矩阵偏差不等式给出 Exercise 5.14 的另一种解法。请量化成功概率以及对次高斯范数的依赖。

查看学习笔记:Exercise 9.10 证明
Exercise 9.11加性 Johnson-Lindenstrauss 引理

说明 Lemma 9.2.4 结论中的误差一般必须是绝对误差,而不能改成相对误差。

查看学习笔记:Exercise 9.11 证明
Exercise 9.12随机仿射截面的 $M^*$ 界

我们已经对经过原点的随机截面证明了 $M^*$ 界(Theorem 9.3.1)。把它推广到所有仿射截面:

$$ \mathbb E\max_{z\in\mathbb R^n} \operatorname{diam}\bigl(T\cap(z+\ker A)\bigr) \le \frac{CK^2w(T)}{\sqrt m}. $$ 查看学习笔记:Exercise 9.12 证明
Exercise 9.13高概率 $M^*$ 界

证明 $M^*$ 界(Theorem 9.3.1)的高概率版本。

查看学习笔记:Exercise 9.13 证明
Exercise 9.14切割 $\ell^p$ 球

令 $B_p^n$ 是 $\mathbb R^n$ 中的 unit $\ell^p$ 球。令 $E$ 是随机 $k$ 维子空间,其中 $1\le k\le0.99n$。

(a) 证明对所有 $p\in(1,\infty)$,

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

其中 $\asymp_p$ 隐含的正数常数只允许依赖于 $p$。

(b) 验证:当 $p\le2$ 时,这个量等价于 $B_p^n$ 的内切欧氏球半径;当 $p\ge2$ 时,它等价于外接欧氏球半径。你能给出直观解释吗?

查看学习笔记:Exercise 9.14 证明
Exercise 9.15逃逸定理的紧性

通过令 $T$ 为 $\mathbb R^n$ 中某个子空间里的单位球面,说明逃逸定理(Theorem 9.3.4)对所有 $m\le n$ 在一般情形下都是最优的。

查看学习笔记:Exercise 9.15 证明
Exercise 9.16在球面上贴标签

这是逃逸定理的另一个版本。令 $T\subset S^{n-1}$ 为任意集合,并令 $\mathcal X\subset S^{n-1}$ 为一个 $N$ 点集。假设

$$ \sigma_{n-1}(T)\lt \frac1N, $$

其中 $\sigma_{n-1}$ 是球面上的 normalized 表面积。证明存在一个 rotation $U\in O(n)$,使得

$$ UT\cap\mathcal X=\varnothing. $$ 查看学习笔记:Exercise 9.16 证明
Exercise 9.17约束恢复:均方误差

把 Theorem 9.4.4 的误差界推广到更强的均方误差:

$$ \mathbb E\|\hat x-x\|_2^2. $$ 查看学习笔记:Exercise 9.17 证明
Exercise 9.18通过优化恢复

令 $T$ 是 $\mathbb R^n$ 中某个范数 $\|\cdot\|_T$ 的单位球。证明 Theorem 9.4.4 的结论也适用于如下优化程序:

$$ \text{minimize }\|x'\|_T \quad\text{subject to}\quad y=Ax'. $$ 查看学习笔记:Exercise 9.18 证明
Exercise 9.19约束优化

把约束恢复结果(Theorem 9.4.4)推广到带噪声模型 (9.16):

$$ y=Ax+w,\qquad x\in T, $$

其中 $w$ 是未知噪声向量,甚至可以依赖于 $A$。为了恢复 $x$,在 $T$ 中找一个尽可能拟合观测 $y$ 的向量:

$$ \text{minimize }\|y-Ax'\|_2 \quad\text{subject to}\quad x'\in T. $$

证明解 $\hat x$ 满足

$$ \mathbb E\|\hat x-x\|_2 \lesssim \frac{K^2w(T)+\|w\|_2}{\sqrt m}. $$ 查看学习笔记:Exercise 9.19 证明
Exercise 9.20无约束优化

把 Remark 9.4.7 更严谨地写出来。令 $x\in\mathbb R^n$ 为任意向量,令 $A$ 是 $m\times n$ 随机矩阵,其行 $A_i$ 独立、各向同性且次高斯。考虑带噪线性模型

$$ y=Ax+w, $$

其中 $w\in\mathbb R^m$ 是任意噪声向量,可以依赖于 $A$。从 $y$ 和 $A$ 恢复 $x$ 时,考虑无约束凸优化问题

$$ \text{minimize }\|y-Ax'\|_2^2+\lambda\|x'\|_T \quad\text{over all }x'\in\mathbb R^n, $$

其中 $\|\cdot\|_T$ 是任意范数,$\lambda\gt 0$ 是调节参数。证明如果选择

$$ \lambda\asymp\frac{\|w\|_2^2}{\|x\|_T}, $$

则解 $\hat x$ 满足

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

其中 $T$ 是范数 $\|\cdot\|_T$ 的单位球,$K=\max_i\|A_i\|_{\psi_2}$。

查看学习笔记:Exercise 9.20 证明
Exercise 9.21Lasso

把 Exercise 9.20 专门化到 $\ell^1$ 范数,分析稀疏回归中常用的 Lasso。设 $x\in\mathbb R^n$ 是 $s$-稀疏,且 $A$ 是 $m\times n$ 随机矩阵,其行 $A_i$ 独立、各向同性且次高斯。考虑带噪线性模型

$$ y=Ax+w. $$

从 $y$ 和 $A$ 恢复 $x$ 时,考虑无约束凸优化问题

$$ \text{minimize }\|y-Ax'\|_2^2+\lambda\|x'\|_1 \quad\text{over all }x'\in\mathbb R^n. $$

证明如果选择 $\lambda\asymp\|w\|_2^2/\|x\|_1$,则解 $\hat x$ 满足

$$ \mathbb E\|\hat x-x\|_2 \lesssim \frac{K^2\sqrt{s\log n}+\|w\|_2}{\sqrt m}, $$

其中 $K=\max_i\|A_i\|_{\psi_2}$。

查看学习笔记:Exercise 9.21 证明
Exercise 9.22稀疏恢复问题适定

令 $A$ 是一个处于一般位置的 $m\times n$ 矩阵。你可以自行选择一个方便的一般位置定义。

(a) 证明如果 $m\ge2\|x\|_0$,那么方程 $y=Ax$ 的解如果存在就是唯一的。

(b) 在这种情况下,如果你已知 $x$ 的支撑集,也就是哪些元素非零,那么如何高效地算法求出 $x$?

查看学习笔记:Exercise 9.22 证明
Exercise 9.23$0\le p\lt 1$ 时的 $\ell^p$ “范数”

在 (9.21) 中,我们把 $\|x\|_0$ 定义为向量 $x$ 的非零坐标个数。

(a) 验证 $\|\cdot\|_0$ 不是 $\mathbb R^n$ 上的范数。

(b) 验证当 $0\lt p\lt 1$ 时,$\|\cdot\|_p$ 不是 $\mathbb R^n$ 上的范数。Figure 9.10 说明此时单位球不是 convex。

(c) 证明对每个 $x\in\mathbb R^n$,

$$ \|x\|_0 = \lim_{p\to0+}\|x\|_p^p. $$
lp balls for p<1
Figure 9.10:$\ell^p$ 单位球在 $0\lt p\lt 1$ 时非凸,因此 $\|\cdot\|_p$ 不是范数。
查看学习笔记:Exercise 9.23 证明
Exercise 9.24近似稀疏恢复

把稀疏恢复(Corollary 9.4.8)推广到近似稀疏信号。

(a) 证明任意 $s$-稀疏信号 $x$ 都可以从观测 $y=Ax$ 中通过求解

$$ \text{minimize }\|x'\|_1 \quad\text{subject to}\quad y=Ax' $$

来恢复,其解 $\hat x$ 满足

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

(b) 说明近似稀疏信号也有类似结果。请陈述并证明一个这样的保证。

查看学习笔记:Exercise 9.24 证明
Exercise 9.25稀疏向量集合的凸化

考虑 unit $s$-稀疏向量的集合

$$ S_{n,s}:= \{x\in\mathbb R^n:\|x\|_0\le s,\ \|x\|_2\le1\}, \tag{9.46} $$

以及 truncated $\ell^1$ 球

$$ T_{n,s}:=\sqrt s B_1^n\cap B_2^n = \{x\in\mathbb R^n:\|x\|_1\le\sqrt s,\ \|x\|_2\le1\}. \tag{9.47} $$

证明

$$ \operatorname{conv}(S_{n,s}) \subset T_{n,s} \subset 2\operatorname{conv}(S_{n,s}). $$

为了证明第二个包含关系,固定 $x\in T_{n,s}$,并把 $x$ 的支撑集分解成互不相交的子集 $I_1,I_2,\ldots$,其中 $I_1$ 对应 $x$ 中绝对值最大的 $s$ 个系数,$I_2$ 对应接下来 $s$ 个系数,依此类推。证明

$$ \sum_{i\ge1}\|x_{I_i}\|_2\le2, $$

其中 $x_I$ 表示 $x$ 限制在坐标集合 $I$ 上的向量。

查看学习笔记:Exercise 9.25 证明
Exercise 9.26稀疏恢复中的对数改进

使用 Exercise 9.25 证明

$$ w(T_{n,s}) \le 2w(S_{n,s}) \le C\sqrt{s\log(en/s)}. $$

把稀疏恢复保证(Corollary 9.4.8)中的对数因子改进为

$$ \mathbb E\|\hat x-x\|_2 \le CK^2\sqrt{\frac{s\log(en/s)}{m}}, $$

从而说明 $m\gtrsim s\log(en/s)$ 个观测足够。

查看学习笔记:Exercise 9.26 证明
Exercise 9.27稀疏向量的高斯宽度

证明 unit $s$-稀疏向量集合 (9.46) 及其近似凸包 (9.47) 具有如下高斯宽度:

$$ w(T_{n,s}) \asymp w(S_{n,s}) \asymp \sqrt{s\log(en/s)}, $$

其中 $\asymp$ 隐含 absolute positive constants。

查看学习笔记:Exercise 9.27 证明
Exercise 9.28Garnaev-Gluskin 定理

改进 Example 9.3.2 中随机切片 $\ell^1$ 球的对数因子:

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

特别地,(9.13) 中的对数因子可以被去掉。

查看学习笔记:Exercise 9.28 证明
Exercise 9.29低秩矩阵恢复的扩展

为第 9.4.3 节中的低秩恢复推出若干版本。

(a) 证明 Corollary 9.4.11 的结论也适用于如下优化程序:

$$ \text{minimize }\|X'\|_* \quad\text{subject to}\quad y_i=\langle A_i,X'\rangle,\quad i=1,\ldots,m. $$

(b) 把矩阵恢复结果推广到近似低秩矩阵。

(c) 把矩阵恢复结果推广到 $d_1\times d_2$ 矩阵。

查看学习笔记:Exercise 9.29 证明
Exercise 9.30精确稀疏恢复的几何

给出 Theorem 9.5.1 证明的几何解释,参见 Figure 9.7b。该证明对切锥 $T(x)$ 及其球面部分 $S(x)$ 说明了什么?

查看学习笔记:Exercise 9.30 证明
Exercise 9.31带噪观测

把精确稀疏恢复结果(Theorem 9.5.1)推广到带噪观测 $y=Ax+w$。请相应修改恢复程序 (9.27)。

查看学习笔记:Exercise 9.31 证明
Exercise 9.32零空间性质

下面是保证精确恢复的一个有用条件。称一个 $m\times n$ 矩阵 $A$ 满足阶 $s$ 的零空间性质,如果对任意非零 $h\in\ker(A)$ 和任意 $s$ 元子集 $S\subset\{1,\ldots,n\}$,都有

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

(a) 证明 $A$ 满足零空间性质当且仅当每个 $s$-稀疏向量 $x\in\mathbb R^n$ 都是 (9.27) 在 $y=Ax$ 时的唯一解。

(b) 证明当 $m\gtrsim s\log(en/s)$ 时,随机矩阵 $A$ 以高概率满足零空间性质。请像 Theorem 9.5.7 那样把该结论精确表述出来。

查看学习笔记:Exercise 9.32 证明
Exercise 9.33随机投影满足 RIP

令 $P$ 是 $\mathbb R^n$ 到一个随机 $m$ 维子空间上的正交投影,该子空间在 Grassmannian 上均匀分布。

(a) 证明 $P$ 满足 RIP,形式与 Theorem 9.5.7 类似,只差一个归一化。

(b) 推出来自随机投影的 Theorem 9.5.1 精确恢复版本。

查看学习笔记:Exercise 9.33 证明
Exercise 9.34次可加性

令 $f:V\to\mathbb R$ 是向量空间 $V$ 上的次可加函数。证明

$$ f(x)-f(y)\le f(x-y) \quad\text{for all }x,y\in V. $$ 查看学习笔记:Exercise 9.34 证明
Exercise 9.35各向异性分布的一般矩阵偏差

把 Theorem 9.6.3 推广到 $m\times n$ 矩阵 $A$ 的行是独立 $N(0,\Sigma)$ 随机向量的情形。证明

$$ \mathbb E\sup_{x\in T} \left| f(Ax)-\mathbb Ef(Ax) \right| \le Cb\,\gamma(\Sigma^{1/2}T). $$ 查看学习笔记:Exercise 9.35 证明
Exercise 9.36一般矩阵偏差的高概率界

证明 Theorem 9.6.3 的高概率版本。

查看学习笔记:Exercise 9.36 证明
Exercise 9.37一般范数下的 Johnson-Lindenstrauss 引理

使用一般矩阵偏差不等式(Theorem 9.6.3),为 $\mathbb R^m$ 上任意范数推出 Johnson-Lindenstrauss 引理的版本,而不仅限于欧氏范数。

查看学习笔记:Exercise 9.37 证明
Exercise 9.38$\ell^1$ 范数下的 Johnson-Lindenstrauss 引理

把 Exercise 9.37 专门化到 $\ell^1$ 范数。取 $\mathcal X$ 为 $\mathbb R^n$ 中 $N$ 个点组成的集合,令 $A$ 为 $m\times n$ 高斯矩阵,其元素独立且服从 $N(0,1)$,并取任意 $\varepsilon\in(0,1)$。假设

$$ m\ge C(\varepsilon)\log N. $$

证明以高概率,矩阵

$$ Q=\sqrt{\frac\pi2}\,\frac1m A $$

满足

$$ (1-\varepsilon)\|x-y\|_2 \le \|Qx-Qy\|_1 \le (1+\varepsilon)\|x-y\|_2 \quad\text{for all }x,y\in\mathcal X. $$ 查看学习笔记:Exercise 9.38 证明
Exercise 9.39Johnson-Lindenstrauss 嵌入到 $\ell^\infty$

把 Exercise 9.37 专门化到 $\ell^\infty$ 范数。取 $\mathcal X$ 为 $\mathbb R^n$ 中 $N$ 个点组成的集合,令 $A$ 为 $m\times n$ 高斯矩阵,其元素独立且服从 $N(0,1)$,并取任意 $\varepsilon\in(0,1)$。假设

$$ m\ge N^{C(\varepsilon)}. $$

证明以高概率,矩阵

$$ Q=C(\log m)^{-1/2}A $$

满足

$$ (1-\varepsilon)\|x-y\|_2 \le \|Qx-Qy\|_\infty \le (1+\varepsilon)\|x-y\|_2 \quad\text{for all }x,y\in\mathcal X, $$

其中绝对常数 $C$ 可适当选择。注意在此情形下 $m\ge N$,所以 $Q$ 不是降维映射,而是嵌入。

查看学习笔记:Exercise 9.39 证明
Exercise 9.40对偶性

证明:对闭且有界的集合 $V\subset\mathbb R^m$,它的支撑函数接近欧氏范数当且仅当其凸包接近欧氏球。具体地,对 $r_-,r_+\ge0$,证明

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

当且仅当

$$ r_-\|y\|_2 \le \sup_{x\in V}\langle x,y\rangle \le r_+\|y\|_2 \quad\text{for all }y\in\mathbb R^m. $$ 查看学习笔记:Exercise 9.40 证明
Exercise 9.41高概率 Dvoretzky-Milman 定理

陈述并证明 Dvoretzky-Milman 定理的高概率版本。

查看学习笔记:Exercise 9.41 证明
Exercise 9.42高斯点云近似圆

考虑 i.i.d. 随机向量 $g_1,\ldots,g_n\sim N(0,I_m)$。假设

$$ m\le c\log n. $$

证明这些点的凸包以高概率近似为半径 $\asymp\sqrt{\log n}$ 的欧氏球,参见 Figure 9.9。

查看学习笔记:Exercise 9.42 证明
Exercise 9.43真实随机投影

我们把 Dvoretzky-Milman 定理表述为 “高斯 projections” 的版本。请证明它对真正的随机投影也成立:令 $P$ 是 $\mathbb R^n$ 到一个随机 $m$ 维子空间上的正交投影。在相同假设下,结论应为

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

其中 $B$ 是半径为 $w_s(T)$ 的欧氏球,$w_s(T)$ 是 $T$ 的球面宽度(回忆 Definition 7.5.4)。

查看学习笔记:Exercise 9.43 证明
学习笔记 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 视为进入更高级随机矩阵与几何泛函分析应用的基础模板。