LLM WIKI · 课程精读

LEARNING UNIT · 18

扩散过程的高维几何

用旋转插值、球壳集中和非均匀角度理解 DDPM 的加噪与采样。

已整理章节
4 节
单元来源
3 条视频
总时长
11:02
状态
已发布
学习位置
18 / 20
01

主题讲解 · 03:52

DDPM 前向加噪:从单位圆系数到椭圆轨迹

学习目标

  • 写出 DDPM 任意时间步的闭式前向加噪公式。
  • 解释“九十度旋转”首先描述的是哪一对系数。
  • 从二维张成子空间理解高维向量的前向轨迹。
  • 区分单位圆系数、椭圆轨迹、几何极角与刚体旋转。
  • 说清近似正交、等范数和终点为纯噪声各自需要的条件。

前置与衔接

DDPM 训练时, 先从数据分布采样干净样本 x0x_0, 再采样标准高斯噪声:

ϵN(0,I).\epsilon\sim\mathcal N(0,I).

随后随机选择一个时间步 tt, 构造带噪样本 xtx_t, 让网络根据 xtx_ttt 预测所加噪声。

图 1

DDPM 训练对每个样本抽取一个时间步和一份噪声,再让网络预测这份噪声。

原视频 · 00:20 ↗

常见的简化噪声预测目标为:

Lmathrmsimple=Ex0,ϵ,t[ϵϵθ(xt,t)22].\mathcal L_{mathrm{simple}} = \mathbb E_{x_0,\epsilon,t} \left[ \left\| \epsilon-\epsilon_\theta(x_t,t) \right\|_2^2 \right].

这里有一个容易混淆的点:

训练并不需要先生成 x1x_1、 再从 x1x_1 生成 x2x_2, 一路模拟到 xtx_t

只要抽取一份 epsilonepsilon, 就能直接从闭式边缘分布采样 xtx_t

本课要研究的, 正是这个闭式采样公式背后的几何图像。

核心讲解

1. 任意时间步的闭式前向公式

令每一步的信号保留率为:

αt=1βt,\alpha_t=1-\beta_t,

累计信号保留率为:

αˉt=s=1tαs.\bar\alpha_t = \prod_{s=1}^{t}\alpha_s.

DDPM 的前向边缘分布可以写成:

q(xtx0)=N(αˉtx0,(1αˉt)I).q(x_t\mid x_0) = \mathcal N \left( \sqrt{\bar\alpha_t}x_0, (1-\bar\alpha_t)I \right).

因此可重参数化为:

xt=αˉtx0+1αˉtϵ.x_t = \sqrt{\bar\alpha_t}x_0 + \sqrt{1-\bar\alpha_t}\epsilon.
图 2

闭式前向边缘分布把任意时间步样本写成干净样本与标准高斯噪声的线性组合。

原视频 · 01:00 ↗

这条公式同时说明:

  • αˉt\sqrt{\bar\alpha_t} 控制保留多少干净信号。
  • 1αˉt\sqrt{1-\bar\alpha_t} 控制注入多少噪声。
  • 对固定的 x0x_0epsilonepsilon,所有 xtx_t 都位于二者张成的子空间。

x0x_0epsilonepsilon 线性无关, 这个子空间最多是二维平面:

span{x0,ϵ}.\operatorname{span}\{x_0,\epsilon\}.

2. 为什么系数落在单位圆上

定义系数:

at=αˉt,bt=1αˉt.a_t=\sqrt{\bar\alpha_t}, \qquad b_t=\sqrt{1-\bar\alpha_t}.

它们满足精确恒等式:

at2+bt2=αˉt+(1αˉt)=1.a_t^2+b_t^2 = \bar\alpha_t+(1-\bar\alpha_t) =1.
图 3

两项系数的平方和为 1,因此系数对沿第一象限单位圆移动。

原视频 · 01:20 ↗

因为两个系数都非负, 可引入一个系数角 θt\theta_t

at=cosθt,bt=sinθt,θt[0,π/2].a_t=\cos\theta_t, \qquad b_t=\sin\theta_t, \qquad \theta_t\in[0,\pi/2].

于是:

xt=cosθtx0+sinθtϵ.x_t = \cos\theta_t\,x_0 + \sin\theta_t\,\epsilon.

αˉt\bar\alpha_t 从接近 11 降到接近 00 时, 系数点从 (1,0)(1,0) 移向 (0,1)(0,1)

这就是“九十度旋转插值”最严格的含义:

它尚未证明高维向量 xtx_t 本身做了保持长度的刚体旋转。

3. 高维近似正交的用途与边界

为了把张成平面画成直角坐标系, 视频使用了高维近似正交直觉。

对固定的非零 x0Rdx_0\in\mathbb R^d, 若:

ϵN(0,Id),\epsilon\sim\mathcal N(0,I_d),

则内积:

x0ϵN(0,x022).x_0^\top\epsilon \sim \mathcal N(0,\|x_0\|_2^2).

归一化余弦为:

x0ϵx02ϵ2.\frac{x_0^\top\epsilon} {\|x_0\|_2\|\epsilon\|_2}.

由于 ϵ2\|\epsilon\|_2 通常是 d\sqrt d 量级, 该余弦的典型幅度是 O(d1/2)O(d^{-1/2})

所以维数很大时, 随机噪声方向通常接近与固定样本垂直。

图 4

固定样本与各向同性高斯噪声在高维中通常近似正交,但并非逐次采样都严格垂直。

原视频 · 02:00 ↗

但需要保留三个边界:

  1. 这是概率意义下的近似,不是每份噪声都严格正交。
  2. 正交不等于两者范数相等。
  3. 即使不正交,轨迹仍位于同一个二维张成子空间。

4. 从系数单位圆到向量椭圆

固定一对 x0x_0epsilonepsilon, 考虑线性映射:

L(a,b)=ax0+bϵ.L(a,b)=ax_0+b\epsilon.

单位圆上的系数点 (a,b)(a,b), 经 LL 映射到数据空间中的 xtx_t

图 5

在干净样本与噪声张成的平面内,二者可作为两条方向来表示前向轨迹。

原视频 · 02:40 ↗

线性映射通常把圆变为椭圆。

若进一步采用正交近似, 并取单位方向:

u=x0x02,v=ϵϵ2,uv0,u=\frac{x_0}{\|x_0\|_2}, \qquad v=\frac{\epsilon}{\|\epsilon\|_2}, \qquad u^\top v\approx0,

xtx_t(u,v)(u,v) 坐标中的分量约为:

Xt=αˉtx02,X_t = \sqrt{\bar\alpha_t}\|x_0\|_2,
Yt=1αˉtϵ2.Y_t = \sqrt{1-\bar\alpha_t}\|\epsilon\|_2.

消去 αˉt\bar\alpha_t 得:

Xt2x022+Yt2ϵ22=1.\frac{X_t^2}{\|x_0\|_2^2} + \frac{Y_t^2}{\|\epsilon\|_2^2} =1.
图 6

单位圆系数经过由两条向量定义的线性映射后形成椭圆轨迹;端点分别接近干净样本和纯噪声。

原视频 · 03:20 ↗

椭圆两条半轴分别由:

x02,ϵ2\|x_0\|_2, \qquad \|\epsilon\|_2

决定。

只有当二者近似正交且范数相等时, 这个椭圆才退化为圆, 此时才更接近通常意义上的等半径旋转。

5. 系数角不一定是几何极角

系数角定义为:

tanθt=1αˉtαˉt.\tan\theta_t = \sqrt{\frac{1-\bar\alpha_t}{\bar\alpha_t}}.

在正交坐标近似下, 实际向量 xtx_t 在张成平面中的几何极角 ϕt\phi_t 满足:

tanϕt=YtXt=1αˉtαˉtϵ2x02.\tan\phi_t = \frac{Y_t}{X_t} = \sqrt{\frac{1-\bar\alpha_t}{\bar\alpha_t}} \frac{\|\epsilon\|_2}{\|x_0\|_2}.

因此只有在:

x02=ϵ2\|x_0\|_2=\|\epsilon\|_2

时, ϕt\phi_t 才与 θt\theta_t 一致。

x0x_0ϵ\epsilon 不正交, 还需使用它们的 Gram 矩阵描述几何, 不能直接把两条向量当成直角坐标轴。

6. 两个端点只是“接近”还是“等于”

理想起点若满足:

αˉ0=1,\bar\alpha_0=1,

则:

x0=1x0+0ϵ.x_0=1\cdot x_0+0\cdot\epsilon.

有限终点 TT 通常只要求:

αˉT0.\bar\alpha_T\approx0.

于是:

xTϵ.x_T\approx\epsilon.

除非调度明确令 αˉT=0\bar\alpha_T=0, 否则终点仍保留极少量信号, 不能把近似写成严格等号。

跟练与练习

跟练:从系数恒等式推到椭圆

已知:

xt=atx0+btϵ,at2+bt2=1.x_t=a_tx_0+b_t\epsilon, \qquad a_t^2+b_t^2=1.

假设 x0ϵx_0\perp\epsilon, 并使用单位方向:

u=x0x02,v=ϵϵ2.u=\frac{x_0}{\|x_0\|_2}, \qquad v=\frac{\epsilon}{\|\epsilon\|_2}.

完成以下三步:

  1. 写出 xtx_tu,vu,v 上的坐标。
  2. 从坐标解出 at,bta_t,b_t
  3. 代回 at2+bt2=1a_t^2+b_t^2=1

最终应得到:

Xt2x022+Yt2ϵ22=1.\frac{X_t^2}{\|x_0\|_2^2} + \frac{Y_t^2}{\|\epsilon\|_2^2} =1.
编者练习

x0ϵx_0\perp\epsilon, 但: ϵ2=2x02,\|\epsilon\|_2=2\|x_0\|_2, 当系数角 θt=45\theta_t=45^\circ 时, 实际几何极角 ϕt\phi_t 是多少?

查看参考答案

θt=45\theta_t=45^\circ 得:
tanθt=1.\tan\theta_t=1.
几何极角满足:
tanϕt=tanθtϵ2x02=2.\tan\phi_t = \tan\theta_t \frac{\|\epsilon\|_2}{\|x_0\|_2} =2.
因此:
ϕt=arctan263.4.\phi_t=\arctan2\approx63.4^\circ.
这说明系数走到单位圆的 4545^\circ
并不意味着数据空间中的向量也转到 4545^\circ

常见误区

误区 1:训练必须逐步执行所有加噪转移

训练可用闭式边缘分布直接构造任意 xtx_t

逐步马尔可夫链定义了过程, 但不要求每个训练样本都实际模拟前 tt 步。

误区 2:xtx_t 本身一定在做刚体旋转

严格落在单位圆上的是系数对:

(αˉt,1αˉt).(\sqrt{\bar\alpha_t},\sqrt{1-\bar\alpha_t}).

数据空间轨迹是该单位圆经过线性映射后的像, 通常为椭圆。

误区 3:高维随机向量严格正交

高维集中只说明归一化内积通常接近零。

它不是代数恒等式, 更不表示任意两条高维向量都正交。

误区 4:正交就等于圆周运动

正交只保证可以使用直角坐标轴。

若两条方向的范数不同, 轨迹仍是椭圆而不是圆。

误区 5:有限 TT 时一定得到严格纯噪声

常用调度令 αˉT\bar\alpha_T 很小, 所以 xTx_T 近似标准高斯噪声。

是否严格等于噪声, 取决于终点调度是否真的令 αˉT=0\bar\alpha_T=0

本课小结

  • DDPM 可直接采样 xt=αˉtx0+1αˉtϵx_t=\sqrt{\bar\alpha_t}x_0+\sqrt{1-\bar\alpha_t}\epsilon
  • 两个非负系数的平方和为一,因此位于第一象限单位圆。
  • 对固定 x0,ϵx_0,\epsilon,所有前向样本位于二者张成的二维子空间。
  • 高维近似正交让二维坐标图更直观,但不是严格前提。
  • 系数单位圆经线性映射通常成为椭圆,半轴尺度由两条向量的范数决定。
  • 系数角、实际几何极角和刚体旋转只有在额外条件下才会重合。
  • “九十度旋转插值”应视为有用的系数与平面几何图像,而不是无条件的字面旋转。
02

主题讲解 · 03:30

DDPM 噪声调度:为何均匀时间不等于均匀角度

学习目标

  • 区分 βt\beta_tαt\alpha_tαˉt\bar\alpha_t
  • 从累计乘积解释线性 beta 为何不会产生线性角度。
  • 区分系数角与数据平面中的几何极角。
  • 推导一个让系数角均匀变化的简化余弦构造。
  • 识别简化余弦构造和 Improved DDPM cosine schedule 的差别。

前置与衔接

上一课把 DDPM 前向边缘采样写成:

xt=αˉtx0+1αˉtϵ,ϵN(0,I).x_t=\sqrt{\bar\alpha_t}x_0+\sqrt{1-\bar\alpha_t}\epsilon, \qquad \epsilon\sim\mathcal N(0,I).

两个系数满足:

αˉt+(1αˉt)=1.\bar\alpha_t+(1-\bar\alpha_t)=1.

因此可以定义系数角 θt\theta_t

cosθt=αˉt,sinθt=1αˉt.\cos\theta_t=\sqrt{\bar\alpha_t}, \qquad \sin\theta_t=\sqrt{1-\bar\alpha_t}.

本课的问题不是轨迹是否位于椭圆上, 而是:

为了画出二维几何图, 视频继续采用两条高维向量近似正交的直觉。

图 1

前向轨迹位于干净样本与噪声张成的二维平面,高维近似正交只用于简化坐标解释。

原视频 · 01:00 ↗

需要先保留边界:

  • 轨迹位于 span{x0,ϵ}\operatorname{span}\{x_0,\epsilon\} 是精确的。
  • x0x_0ϵ\epsilon 当作直角坐标轴是高维近似。
  • 把系数角等同于实际几何角还需要范数相等。

核心讲解

1. 椭圆上的位置由累计系数决定

若使用正交单位方向:

u=x0x02,v=ϵϵ2,u=\frac{x_0}{\|x_0\|_2}, \qquad v=\frac{\epsilon}{\|\epsilon\|_2},

xtx_t 的二维坐标约为:

Xt=αˉtx02,X_t=\sqrt{\bar\alpha_t}\|x_0\|_2,
Yt=1αˉtϵ2.Y_t=\sqrt{1-\bar\alpha_t}\|\epsilon\|_2.

消去时间步可得椭圆:

Xt2x022+Yt2ϵ22=1.\frac{X_t^2}{\|x_0\|_2^2} + \frac{Y_t^2}{\|\epsilon\|_2^2} =1.
图 2

时间步增加时,样本沿椭圆从干净端点向噪声端点推进,但推进角度未必匀速。

原视频 · 01:20 ↗

随着 tt 增加, 通常有:

αˉt,\bar\alpha_t\downarrow,

所以 XtX_t 减小、YtY_t 增大。

这只确定了行进方向, 没有说明每一步跨过相同角度或相同弧长。

2. 三个 alpha 符号不能混用

DDPM 中常见三个相关量:

βt=第 t 步加入的噪声方差比例,\beta_t = \text{第 }t\text{ 步加入的噪声方差比例},
αt=1βt,\alpha_t=1-\beta_t,
αˉt=s=1tαs=s=1t(1βs).\bar\alpha_t = \prod_{s=1}^{t}\alpha_s = \prod_{s=1}^{t}(1-\beta_s).
图 3

累计信号保留率是各步 alpha 的乘积,线性 beta 会经过乘积变成非线性的累计曲线。

原视频 · 02:00 ↗

即使 βt\beta_ttt 是线性的, αˉt\bar\alpha_t 也是许多线性项的乘积, 不会对 tt 保持线性。

取对数更容易看清:

logαˉt=s=1tlog(1βs).\log\bar\alpha_t=\sum_{s=1}^{t}\log(1-\beta_s).

βs\beta_s 很小时:

log(1βs)βs,\log(1-\beta_s)\approx-\beta_s,

所以:

αˉtexp(s=1tβs).\bar\alpha_t\approx\exp\left(-\sum_{s=1}^{t}\beta_s\right).

βs\beta_s 线性增长, 其累加大致是关于 tt 的二次函数, 再经过指数映射后更不可能成为线性曲线。

3. 传统线性 beta 调度

早期 DDPM 常使用线性 beta:

βt=βmin+t1T1(βmaxβmin).\beta_t =\beta_{min}+\frac{t-1}{T-1}(\beta_{max}-\beta_{min}).

经典像素空间设置常取:

βmin=104,βmax=0.02.\beta_{min}=10^{-4}, \qquad \beta_{max}=0.02.

字幕把起点识别成了“十的四次方”, 但画面和 DDPM 上下文均显示应为:

104.10^{-4}.

这个数值是具体离散步数和数据表示下的经典设置, 不是所有扩散模型都必须照搬的常数。

4. 系数角为何不均匀

由定义:

tanθt=sinθtcosθt=1αˉtαˉt.\tan\theta_t =\frac{\sin\theta_t}{\cos\theta_t} =\sqrt{\frac{1-\bar\alpha_t}{\bar\alpha_t}}.

因此:

θt=arctan1αˉtαˉt.\theta_t =\arctan\sqrt{\frac{1-\bar\alpha_t}{\bar\alpha_t}}.

ttθt\theta_t 至少经过:

  1. βs\beta_s 的逐步调度。
  2. s(1βs)\prod_s(1-\beta_s) 的累计乘积。
  3. 平方根比值。
  4. 反正切映射。

即使第一步是线性的, 后面三步也会让角度变成非线性函数。

图 4

传统线性 beta 调度不会让累计信号系数或几何角度随时间线性变化。

原视频 · 02:40 ↗

这就是均匀抽取时间步并不代表均匀抽取系数角的原因。

5. 几何极角还有一个范数比

视频在椭圆图上把角度记为 θ\theta

为了避免符号混淆, 这里把实际二维坐标的极角记为 ϕt\phi_t

在正交近似下:

tanϕt=YtXt=1αˉtαˉtϵ2x02.\tan\phi_t = \frac{Y_t}{X_t} = \sqrt{ \frac{1-\bar\alpha_t}{\bar\alpha_t} } \frac{\|\epsilon\|_2}{\|x_0\|_2}.
图 5

在正交坐标近似下,几何角的正切还包含噪声与干净样本范数之比。

原视频 · 02:20 ↗

所以:

ϕt=θt\phi_t=\theta_t

还需要:

x02=ϵ2.\|x_0\|_2=\|\epsilon\|_2.

若二者范数不同, 系数角均匀也不保证几何极角均匀。

若二者不正交, 连上述直角坐标正切公式也需要修正。

6. 如何构造均匀的系数角

若只希望系数角从 00π/2\pi/2 均匀变化, 可以直接指定:

θt=tTπ2.\theta_t = \frac{t}{T}\frac{\pi}{2}.

再令:

αˉt=cosθt,\sqrt{\bar\alpha_t}=\cos\theta_t,
1αˉt=sinθt.\sqrt{1-\bar\alpha_t}=\sin\theta_t.

于是:

αˉt=cos2(tTπ2).\bar\alpha_t = \cos^2 \left( \frac{t}{T}\frac{\pi}{2} \right).
图 6

令系数角与时间成正比可得到余弦和正弦系数,但这只是简化构造。

原视频 · 03:00 ↗

对应的逐步 beta 不能直接拿余弦值代替, 而应由累计量反推:

αt=αˉtαˉt1,\alpha_t = \frac{\bar\alpha_t}{\bar\alpha_{t-1}},
βt=1αˉtαˉt1.\beta_t = 1- \frac{\bar\alpha_t}{\bar\alpha_{t-1}}.

7. 与 Improved DDPM cosine schedule 的边界

视频把上述简单构造连接到 Improved DDPM 的 cosine schedule。

这个联系抓住了余弦累计信号率的核心直觉, 但论文实现并不是无偏移的:

αˉt=f(t)f(0),\bar\alpha_t = \frac{f(t)}{f(0)},

其中:

f(t)=cos2(t/T+s1+sπ2).f(t) = \cos^2 \left( \frac{t/T+s}{1+s} \frac{\pi}{2} \right).

常用:

s=0.008.s=0.008.

随后仍需通过累计量之比得到 βt\beta_t, 并在实现中对过大的 beta 做截断。

因此应分成两层理解:视频的 cos2(πt/2T)\cos^2(\pi t/2T) 是教学用的均匀系数角构造;Improved DDPM 使用带偏移、归一化和截断的实际 schedule。二者概念相连,但公式不能逐项画等号。

跟练与练习

跟练:从累计 alpha 反推 beta

给定一条期望的累计曲线 αˉt\bar\alpha_t, 逐步调度可按下面顺序恢复:

  1. 先检查 αˉ0=1\bar\alpha_0=1
  2. 计算 αt=αˉt/αˉt1\alpha_t=\bar\alpha_t/\bar\alpha_{t-1}
  3. 计算 βt=1αt\beta_t=1-\alpha_t
  4. 检查 0<βt<10<\beta_t<1

例如简化余弦累计量:

αˉt=cos2(πt2T)\bar\alpha_t = \cos^2 \left( \frac{\pi t}{2T} \right)

对应:

βt=1cos2(πt/2T)cos2(π(t1)/2T).\beta_t = 1- \frac{ \cos^2(\pi t/2T) }{ \cos^2(\pi(t-1)/2T) }.
编者练习

假设 x0ϵx_0\perp\epsilon, 且: ϵ2=2x02.\|\epsilon\|_2=2\|x_0\|_2. 若采用均匀系数角, 在 t=T/2t=T/2 时, 系数角与几何极角各是多少?

查看参考答案

均匀系数角给出:
θT/2=12π2=π4=45.\theta_{T/2} = \frac{1}{2}\frac{\pi}{2} = \frac{\pi}{4} =45^\circ.
而几何极角满足:
tanϕT/2=tanθT/2ϵ2x02=2.\tan\phi_{T/2} = \tan\theta_{T/2} \frac{\|\epsilon\|_2}{\|x_0\|_2} =2.
因此:
ϕT/2=arctan263.4.\phi_{T/2} = \arctan2 \approx63.4^\circ.
均匀的是系数角,
不是一般椭圆上的几何极角或弧长。

常见误区

误区 1:βt\beta_t 就是 αˉt\bar\alpha_t

βt\beta_t 是单步噪声比例, αt=1βt\alpha_t=1-\beta_t 是单步信号比例, αˉt\bar\alpha_t 是多步信号比例的累计乘积。

三者不能互换。

误区 2:线性 beta 意味着线性加噪效果

真正进入闭式前向公式的是 αˉt\bar\alpha_t

它由连乘产生, 通常会呈现明显非线性。

误区 3:均匀时间意味着均匀角度

从时间到系数角需要经过累计乘积和反三角函数。

均匀时间步只表示索引间隔相同。

误区 4:系数角就是椭圆几何极角

几何极角还包含两条轴的尺度比。

只有正交且等范数时, 二者才重合。

误区 5:余弦调度就是一个固定公式

教学用的无偏移 cos2\cos^2 构造便于理解。

Improved DDPM 的实际 schedule 还包含偏移、归一化和 beta 截断。

实现时必须核对所用代码库的具体定义。

本课小结

  • DDPM 的角度进度由累计量 αˉt\bar\alpha_t 决定,而非由单步 βt\beta_t 直接决定。
  • 线性 beta 经连乘、平方根和反正切后,系数角不会随时间均匀变化。
  • 数据平面中的几何极角还受 ϵ2/x02\|\epsilon\|_2/\|x_0\|_2 影响。
  • 直接令 θt=πt/(2T)\theta_t=\pi t/(2T) 可构造均匀系数角,并得到简化的 cos2\cos^2 累计曲线。
  • 逐步 beta 应从相邻累计 alpha 的比值反推。
  • Improved DDPM cosine schedule 与简化构造直觉一致,但实际公式带偏移、归一化和截断。
  • “旋转是否均匀”必须先说明讨论的是系数角、几何极角还是椭圆弧长。
03

主题讲解 · 03:40

高维高斯薄壳:噪声为何集中在根号 d 附近

学习目标

  • 从标准多元高斯的坐标分布推导噪声范数尺度。
  • 区分 Eϵ22\mathbb E\|\epsilon\|_2^2(Eϵ2)2(\mathbb E\|\epsilon\|_2)^2
  • 写出高斯半径的精确均值和高维渐近式。
  • 用绝对波动与相对波动解释薄壳集中。
  • 说明“在球壳上选点”是近似图像而非固定半径采样。
  • ϵ2d\|\epsilon\|_2\approx\sqrt d 接回 DDPM 椭圆几何。

前置与衔接

扩散模型经常采样:

ϵN(0,Id).\epsilon\sim\mathcal N(0,I_d).

其均值向量为:

E[ϵ]=0,\mathbb E[\epsilon]=0,

协方差矩阵为:

Cov(ϵ)=Id.\operatorname{Cov}(\epsilon)=I_d.
图 1

标准多元高斯的均值为零、协方差为单位阵,各坐标方差为一且彼此不相关。

原视频 · 00:20 ↗

单位阵意味着:

Var(ϵi)=1,\operatorname{Var}(\epsilon_i)=1,

并且当 iji\ne j 时:

Cov(ϵi,ϵj)=0.\operatorname{Cov}(\epsilon_i,\epsilon_j)=0.

因为这里是联合高斯分布, 零协方差进一步推出各坐标相互独立:

ϵiiidN(0,1).\epsilon_i\overset{\mathrm{iid}}\sim\mathcal N(0,1).

本课关注的随机变量是高斯向量到原点的半径:

R=ϵ2.R=\|\epsilon\|_2.

核心问题是:

核心讲解

1. 平方半径服从卡方分布

欧氏范数平方为:

R2=ϵ22=i=1dϵi2.R^2=\|\epsilon\|_2^2=\sum_{i=1}^{d}\epsilon_i^2.
图 2

高斯向量的平方半径是各坐标平方之和,服从自由度为 d 的卡方分布。

原视频 · 01:20 ↗

每个 ϵi2\epsilon_i^2 服从自由度为 11 的卡方分布, 独立求和得到:

R2χd2.R^2\sim\chi_d^2.

卡方分布给出两个精确结果:

E[R2]=d,\mathbb E[R^2]=d,
Var(R2)=2d.\operatorname{Var}(R^2)=2d.

也可以只用二阶矩推导第一个等式。

对任一坐标:

E[ϵi2]=Var(ϵi)+(Eϵi)2=1.\mathbb E[\epsilon_i^2] =\operatorname{Var}(\epsilon_i)+(\mathbb E\epsilon_i)^2=1.

由期望的线性性:

E[R2]=i=1dE[ϵi2]=d.\mathbb E[R^2] =\sum_{i=1}^{d}\mathbb E[\epsilon_i^2]=d.
图 3

每个标准正态坐标的二阶矩为一,所以平方半径的期望恰好等于维数 d。

原视频 · 01:40 ↗

2. 平方半径期望不等于半径期望的平方

从:

E[R2]=d\mathbb E[R^2]=d

不能直接推出:

E[R]=d.\mathbb E[R]=\sqrt d.

平方根是凹函数, Jensen 不等式给出:

E[R]=E[R2]E[R2]=d.\mathbb E[R] =\mathbb E[\sqrt{R^2}] \le\sqrt{\mathbb E[R^2]}=\sqrt d.

对有限 dd, 严格均值实际上略小于 d\sqrt d

由于 RR 服从自由度为 dd 的 chi 分布, 其精确均值为:

E[R]=2Γ((d+1)/2)Γ(d/2).\mathbb E[R] =\sqrt2\,\frac{\Gamma((d+1)/2)}{\Gamma(d/2)}.

高维渐近展开为:

E[R]=d(114d+O(d2)).\mathbb E[R] =\sqrt d\left(1-\frac{1}{4d}+O(d^{-2})\right).
图 4

高维时半径均值接近根号 d,但平方根与期望不能直接交换。

原视频 · 02:00 ↗

因此视频中的 E[R]d\mathbb E[R]\approx\sqrt d 是高维近似,而不是从 E[R2]=d\mathbb E[R^2]=d 直接开平方得到的恒等式。

3. 为什么半径方差趋近常数

利用:

Var(R)=E[R2](E[R])2,\operatorname{Var}(R) =\mathbb E[R^2]-(\mathbb E[R])^2,

可得精确表达:

Var(R)=d2[Γ((d+1)/2)Γ(d/2)]2.\operatorname{Var}(R) =d-2\left[\frac{\Gamma((d+1)/2)}{\Gamma(d/2)}\right]^2.

代入高维渐近式:

Var(R)12.\operatorname{Var}(R) \to \frac12.

因此标准差约为:

Std(R)12.\operatorname{Std}(R) \approx\frac{1}{\sqrt2}.
图 5

半径的绝对波动约为常数量级,而平均半径随根号 d 增长,因此相对厚度趋近于零。

原视频 · 02:20 ↗

注意它并没有趋近于零。

高维集中来自相对波动:

Std(R)E[R]12d0.\frac{\operatorname{Std}(R)}{\mathbb E[R]} \approx \frac{1}{\sqrt{2d}} \to0.

即半径中心约为 d\sqrt d, 绝对厚度仍是常数量级, 但相对于不断增大的半径, 这一厚度越来越薄。

4. “薄壳”到底是什么意思

可把大部分高斯概率质量想象成集中在:

RdR\approx\sqrt d

附近的一层球壳中。

这个图像有三条重要边界。

第一, 每次采样的半径仍是随机的, 不是严格固定为 d\sqrt d

第二, 高斯样本不是在某个单一球面上均匀采样。

更准确地说, 各向同性高斯的方向:

U=ϵϵ2U=\frac{\epsilon}{\|\epsilon\|_2}

在单位球面上均匀, 且方向 UU 与半径 RR 相互独立; 半径则服从 chi 分布。

第三, 多元高斯密度函数在原点最大, 但概率质量不只由点密度决定。

半径为 rr 的球面面积随:

rd1r^{d-1}

增长, 径向密度与:

rd1er2/2r^{d-1}e^{-r^2/2}

成正比。

其众数位于:

rmode=d1,d1.r_{\mathrm{mode}}=\sqrt{d-1}, \qquad d\ge1.

这解释了为什么原点密度最高, 大部分质量却远离原点并集中在 d\sqrt d 尺度。

5. 薄壳与高维近似正交

固定非零向量 x0x_0, 标准高斯噪声满足:

x0ϵN(0,x022).x_0^\top\epsilon \sim \mathcal N(0,\|x_0\|_2^2).

归一化余弦为:

cosψ=x0ϵx02R.\cos\psi = \frac{x_0^\top\epsilon} {\|x_0\|_2R}.

分子是 x02\|x_0\|_2 量级, 分母中的 RRd\sqrt d 量级, 所以:

cosψ=Op(d1/2).\cos\psi=O_p(d^{-1/2}).

薄壳提供了范数尺度, 各向同性提供了随机方向; 二者共同形成高维近似正交的直觉。

但“近似正交”仍是概率性陈述, 不能替代一次具体采样的真实内积。

6. 接回 DDPM 的椭圆轨迹

DDPM 前向样本为:

xt=αˉtx0+1αˉtϵ.x_t = \sqrt{\bar\alpha_t}x_0 + \sqrt{1-\bar\alpha_t}\epsilon.

x0ϵx_0\perp\epsilon 的近似下, 取二者单位方向为坐标轴, 可写出椭圆:

Xt2x022+Yt2ϵ22=1.\frac{X_t^2}{\|x_0\|_2^2} + \frac{Y_t^2}{\|\epsilon\|_2^2} =1.

利用高维薄壳近似:

ϵ22d,\|\epsilon\|_2^2\approx d,

就得到:

Xt2x022+Yt2d1.\frac{X_t^2}{\|x_0\|_2^2} + \frac{Y_t^2}{d} \approx1.
图 6

用高维噪声范数平方近似 d,可把 DDPM 前向轨迹写成噪声轴尺度为根号 d 的椭圆。

原视频 · 03:20 ↗

这里的近似来自两层:

  1. x0x_0ϵ\epsilon 近似正交。
  2. ϵ22\|\epsilon\|_2^2 集中在 dd 附近。

因此该椭圆是一种高维典型几何, 不是对每一份样本和噪声都成立的精确坐标恒等式。

跟练与练习

跟练:区分三个半径统计量

R=ϵ2R=\|\epsilon\|_2, 分别填写:

E[R2]=,\mathbb E[R^2]=\underline{\qquad},
E[R]=,\mathbb E[R]=\underline{\qquad},
Var(R)=.\operatorname{Var}(R)=\underline{\qquad}.

正确层次是:

  • 第一个是精确值 dd
  • 第二个是 gamma 函数比,高维近似 d\sqrt d
  • 第三个高维趋近于 1/21/2

不要把三个量都用“约等于根号 dd”代替。

编者练习

假设 d=4096d=4096, 使用高维近似: E[R]d,Std(R)12.\mathbb E[R]\approx\sqrt d, \qquad \operatorname{Std}(R)\approx\frac{1}{\sqrt2}. 估算平均半径、标准差和相对标准差。

查看参考答案

平均半径约为:
4096=64.\sqrt{4096}=64.
标准差约为:
120.707.\frac{1}{\sqrt2}\approx0.707.
相对标准差约为:
0.707640.0110=1.10%.\frac{0.707}{64} \approx0.0110 =1.10\%.
半径仍有约 0.70.7 的绝对波动,
但相对于半径 6464 已非常薄。

常见误区

误区 1:E[R2]=d\mathbb E[R^2]=d 所以 E[R]=d\mathbb E[R]=\sqrt d

非线性函数一般不能与期望交换。

有限维下应使用 chi 分布的精确均值; d\sqrt d 是高维近似。

误区 2:方差约为 1/21/2,所以没有集中

薄壳强调的是相对波动趋近于零。

绝对标准差保持常数量级, 而平均半径按 d\sqrt d 增长。

误区 3:高斯就是在半径 d\sqrt d 的球面均匀采样

方向确实在球面上均匀, 但半径仍是服从 chi 分布的随机变量。

高斯质量集中在一层有厚度的球壳, 不是单一球面。

误区 4:原点密度最大,所以样本都在原点附近

单点密度没有计入半径增大带来的球面面积。

径向概率质量由密度与体积因子共同决定, 其峰值在 d1\sqrt{d-1} 附近。

误区 5:DDPM 椭圆中的噪声半轴严格等于 d\sqrt d

对具体噪声样本, 半轴是实际范数 ϵ2\|\epsilon\|_2

d\sqrt d 替换它是高维集中近似。

本课小结

  • 标准高斯噪声的平方半径服从 χd2\chi_d^2,且 Eϵ22=d\mathbb E\|\epsilon\|_2^2=d
  • 半径均值有精确 gamma 函数表达式,高维时才近似 d\sqrt d
  • 半径方差趋近 1/21/2,相对标准差约为 1/2d1/\sqrt{2d}
  • 薄壳表示概率质量在相对意义上集中,不表示半径严格固定。
  • 各向同性高斯的方向在球面上均匀,且与随机半径独立。
  • 原点密度最大与概率质量集中在 d\sqrt d 附近并不矛盾,体积因子连接了两者。
  • 在 DDPM 椭圆中用 dd 替换 ϵ22\|\epsilon\|_2^2 是高维典型近似,而非逐样本恒等式。
04

单元综合

从单位圆系数到高维薄壳:DDPM 前向加噪的几何账本

单元能力目标

完成本单元后,应能在不混淆三种“角度”和两种“半径”的前提下解释 DDPM 前向过程。

具体需要做到:

  • 写出从 x0x_0 直接采样 xtx_t 的闭式公式;
  • 说明系数平方和为 1 与样本范数不变不是同一结论;
  • 判断固定 x0,ϵx_0,\epsilon 的轨迹位于哪个子空间;
  • 区分系数角、数据平面几何极角与椭圆弧长;
  • αˉt\bar\alpha_t 推导角度和单步 βt\beta_t
  • 说明线性 βt\beta_t 为何不产生均匀角速度;
  • χd2\chi^2_dχd\chi_d 分布解释高维高斯薄壳;
  • 区分原点密度最大与概率质量远离原点。

概念连接

1. DDPM 前向过程的闭式采样

定义单步保留系数

αt=1βt,\alpha_t=1-\beta_t,

累计保留系数

αˉt=s=1tαs.\bar\alpha_t = \prod_{s=1}^{t}\alpha_s.

则任意时刻可直接从干净样本 x0x_0 采样:

xt=αˉtx0+1αˉtϵ,ϵN(0,Id).x_t = \sqrt{\bar\alpha_t}x_0 + \sqrt{1-\bar\alpha_t}\epsilon, \qquad \epsilon\sim\mathcal N(0,I_d).

它把数据分量与独立高斯噪声按时间相关系数组合。

2. 两个系数落在单位圆上

at=αˉt,bt=1αˉt.a_t=\sqrt{\bar\alpha_t}, \qquad b_t=\sqrt{1-\bar\alpha_t}.

at2+bt2=1,a_t^2+b_t^2=1,

at,bt0a_t,b_t\ge0

所以系数对 (at,bt)(a_t,b_t) 位于第一象限单位圆。

可以参数化为

at=cosθt,bt=sinθt,a_t=\cos\theta_t, \qquad b_t=\sin\theta_t,

其中

θt=arccosαˉt=arctan1αˉtαˉt.\theta_t = \arccos\sqrt{\bar\alpha_t} = \arctan \sqrt{ \frac{1-\bar\alpha_t} {\bar\alpha_t} }.

3. 系数单位圆不保证样本在球面上

样本范数为

xt2=at2x02+bt2ϵ2+2atbtx0Tϵ.\|x_t\|^2 = a_t^2\|x_0\|^2 +b_t^2\|\epsilon\|^2 +2a_tb_t x_0^T\epsilon.

只有在额外条件下,例如

x0Tϵ=0x_0^T\epsilon=0

x0=ϵ,\|x_0\|=\|\epsilon\|,

范数才与角度无关。

因此“系数在单位圆”是精确结论,“样本绕圆旋转”通常只是条件化直觉。

4. 固定两端时,轨迹位于二维子空间

对固定的 x0x_0ϵ\epsilon

xt=atx0+btϵ.x_t=a_tx_0+b_t\epsilon.

所有 xtx_t 都位于

span{x0,ϵ}\operatorname{span}\{x_0,\epsilon\}

中。

即使数据维度为百万,这条固定噪声实例的轨迹最多也在二维平面内。

若两向量线性相关,轨迹还会退化为一维。

5. 系数圆经线性映射通常变成椭圆

从系数空间到数据空间的映射为

[ab]ax0+bϵ.\begin{bmatrix} a\\b \end{bmatrix} \mapsto ax_0+b\epsilon.

单位圆在一般线性映射下变成椭圆。

x0x_0ϵ\epsilon 正交,椭圆半轴尺度与

x0,ϵ\|x_0\|, \qquad \|\epsilon\|

相关。

若不正交,椭圆还会旋转或剪切。

所以 DDPM 的“旋转插值”更准确地说是系数圆在两向量张成平面中的线性像。

6. 三种进度必须分别命名

讨论“是否均匀旋转”前,要说明是哪一种量:

  1. 系数角 θt\theta_t:由 (at,bt)(a_t,b_t) 在单位圆上的角度决定;
  2. 几何极角xtx_t 在数据平面所选坐标中的方向角;
  3. 椭圆弧长:沿实际轨迹走过的几何距离。

三者只有在额外对称条件下才可能一致。

7. 线性 beta 先经过累计连乘

常见线性 schedule 让

βt\beta_t

tt 线性变化。

但系数角由累计量

αˉt=s=1t(1βs)\bar\alpha_t = \prod_{s=1}^{t}(1-\beta_s)

决定。

βt\beta_tθt\theta_t 依次经过:

  1. αt=1βt\alpha_t=1-\beta_t
  2. 连乘得到 αˉt\bar\alpha_t
  3. 开平方得到 at,bta_t,b_t
  4. 反三角函数得到 θt\theta_t

所以线性 βt\beta_t 一般不会产生线性 θt\theta_t

8. 均匀系数角可反推累计 alpha

若教学上直接规定

θt=π2tT,\theta_t = \frac{\pi}{2} \frac{t}{T},

αˉt=cos2θt=cos2(πt2T).\bar\alpha_t = \cos^2\theta_t = \cos^2 \left( \frac{\pi t}{2T} \right).

单步系数由相邻累计量反推:

αt=αˉtαˉt1,\alpha_t = \frac{\bar\alpha_t} {\bar\alpha_{t-1}},
βt=1αˉtαˉt1.\beta_t = 1- \frac{\bar\alpha_t} {\bar\alpha_{t-1}}.

这是一种均匀系数角的简化构造。

9. 不能把简化余弦构造直接当作实际 schedule

Improved DDPM 的 cosine schedule 与余弦累计量直觉相关,但实际公式通常包含:

  • 时间偏移与起点归一化;
  • 数值截断和离散步长边界。

课程中的 cos2\cos^2 构造用于解释角度,不应不经核对替代具体实现。

视频画面确认线性 β\beta 起点示例为 10410^{-4},不是 10410^4

10. 几何极角还依赖向量范数

x0x_0ϵ\epsilon 正交,并以二者单位方向为坐标轴,则实际坐标为

(atx0,btϵ).(a_t\|x_0\|, b_t\|\epsilon\|).

几何极角 ϕt\phi_t 满足

tanϕt=btϵatx0.\tan\phi_t = \frac{b_t\|\epsilon\|} {a_t\|x_0\|}.

只有两向量等范数时,ϕt=θt\phi_t=\theta_t

因此即使系数角均匀,数据平面几何角也未必均匀。

11. 标准高斯半径的平方服从卡方分布

ϵN(0,Id).\epsilon \sim \mathcal N(0,I_d).

每个分量独立服从标准正态。

平方半径为

R2=ϵ2=i=1dϵi2χd2.R^2 = \|\epsilon\|^2 = \sum_{i=1}^{d}\epsilon_i^2 \sim \chi_d^2.

所以

E[R2]=d,Var(R2)=2d.\mathbb E[R^2]=d, \qquad \operatorname{Var}(R^2)=2d.

12. 半径本身服从 chi 分布

R=ϵχd.R=\|\epsilon\| \sim\chi_d.

精确均值为

E[R]=2Γ((d+1)/2)Γ(d/2).\mathbb E[R] = \sqrt2 \frac{ \Gamma((d+1)/2) }{ \Gamma(d/2) }.

高维渐近为

E[R]d1/2d.\mathbb E[R] \approx \sqrt{d-1/2} \approx \sqrt d.

所以“半径为 d\sqrt d”是高维典型近似,不是逐样本恒等式。

13. 薄壳是相对集中

高维时

Var(R)12,\operatorname{Var}(R) \to\frac12,

半径的绝对标准差保持 O(1)O(1),而均值为 O(d)O(\sqrt d)

相对标准差约为

12d.\frac1{\sqrt{2d}}.

所以概率质量集中在半径约 d\sqrt d、相对厚度很小的球壳内。

薄壳不表示所有样本严格落在同一球面。

14. 原点密度最大与薄壳不矛盾

标准高斯点密度

p(x)ex2/2p(x) \propto e^{-\|x\|^2/2}

在原点最大。

但半径为 rr 的球壳体积因子随

rd1r^{d-1}

增长。

径向概率密度近似为

pR(r)rd1er2/2.p_R(r) \propto r^{d-1}e^{-r^2/2}.

密度下降与球壳面积增长共同作用,使概率质量峰值远离原点,位于约 d\sqrt d

15. 高斯方向与半径分离

标准高斯各向同性。

可写为 ϵ=Ru\epsilon=Ru,其中 u=ϵ/ϵu=\epsilon/\|\epsilon\| 在单位球面上均匀,且与随机半径 RR 独立。

这解释了为什么高维噪声既没有偏好方向,又在半径上形成薄壳。

16. 薄壳怎样修正 DDPM 椭圆直觉

对高维标准高斯,常用

ϵ2d\|\epsilon\|^2\approx d

描述典型样本。

把这个近似代入 DDPM 轨迹范数,可得到更稳定的典型几何图像。

但它仍是概率集中结论:

  • 不同噪声样本半径不同;
  • x0x_0 范数与数据分布有关;
  • x0Tϵx_0^T\epsilon 不是逐样本严格为零;
  • 轨迹不是固定半径球面刚体旋转。

对比与决策

1. 看到“旋转”时先问三个问题

  1. 指的是系数对 (at,bt)(a_t,b_t) 吗?
  2. x0x_0ϵ\epsilon 是否正交、等范数?
  3. 讨论的是角度、极角还是弧长?

只有回答清楚,才能判断“九十度”“匀速”或“圆轨迹”是否准确。

2. 设计 schedule 时从目标累计量出发

若目标是控制信噪比或系数角,应先设计 αˉt\bar\alpha_t,再反推单步 βt\beta_t

直接让 βt\beta_t 线性,只保证单步噪声方差按线性表变化,不保证累计几何进度线性。

3. 使用高维近似前说明量级

  • 精确:R2χd2R^2\sim\chi_d^2E[R2]=d\mathbb E[R^2]=d
  • 渐近:RdR\approx\sqrt d,相对波动约 1/2d1/\sqrt{2d}
  • 错误:每个样本都满足 R=dR=\sqrt d

综合训练

编者练习

x0x_0ϵ\epsilon 正交,x0=2\|x_0\|=2ϵ=6\|\epsilon\|=6,系数角 θ=π/4\theta=\pi/4。求 xtx_t 在两条单位基方向上的坐标、范数平方与几何极角,并判断几何极角是否等于 π/4\pi/4

查看参考答案

a=b=1/2a=b=1/\sqrt2。坐标为 (2/2,6/2)=(2,32)(2/\sqrt2,6/\sqrt2)=(\sqrt2,3\sqrt2),范数平方为 2+18=202+18=20。几何极角满足 tanϕ=(32)/(2)=3\tan\phi=(3\sqrt2)/(\sqrt2)=3,所以 ϕ=arctan3π/4\phi=\arctan3\neq\pi/4。系数角相等不代表数据平面几何角相等,因为两条基向量范数不同。

编者练习 2

希望 TT 步的系数角严格满足 θt=πt/(2T)\theta_t=\pi t/(2T)。写出 αˉt\bar\alpha_tβt\beta_t

查看参考答案

αˉt=cos2(πt/(2T))\bar\alpha_t=\cos^2(\pi t/(2T))。由 αt=αˉt/αˉt1\alpha_t=\bar\alpha_t/\bar\alpha_{t-1},得到 βt=1cos2(πt/(2T))/cos2(π(t1)/(2T))\beta_t=1-\cos^2(\pi t/(2T))/\cos^2(\pi(t-1)/(2T))。实际实现需处理末端数值、偏移与 clipping;这只是无偏移简化构造。

编者练习 3

d=10000d=10000 时,用高维近似估计标准高斯半径均值、半径标准差和相对标准差。

查看参考答案

均值约为 d=100\sqrt d=100,半径方差约为 1/21/2,标准差约为 1/20.7071/\sqrt2\approx0.707。相对标准差约为 0.707/1000.007070.707/100\approx0.00707,即约 0.7%。这说明半径相对集中,但不严格固定。

进入下一单元前

  • 已能写出 DDPM 前向闭式采样与累计 αˉt\bar\alpha_t
  • 已能区分系数单位圆、数据平面椭圆与刚体旋转。
  • 已能区分系数角、几何极角与弧长进度。
  • 已能从目标 αˉt\bar\alpha_t 反推 βt\beta_t
  • 已能用 χd2\chi_d^2χd\chi_d 解释高维薄壳。
  • 已能说明原点密度最大与 d\sqrt d 附近质量集中为何不冲突。
  • 若仍会把线性 β\beta 当作匀角度,回看 P82 的累计链路。
  • 若仍会把噪声看成固定半径球面采样,回看 P83 的精确与渐近结论。
  • 若仍会把系数角当作实际几何极角,回看 P81 的范数与内积条件。