Generative Modeling by Estimating Gradients of the Data Distribution

本文介绍宋飏提出的带噪声扰动的基于得分的生成模型。首先介绍基本的基于得分的生成模型的训练方法（得分匹配）和采样方法（朗之万动力学）。然后基于流形假设，分析了现有方法面临的低密度区域得分估计不准、标准朗之万动力学采样对于多 pattern 分布混合慢的问题。最后提出解决方法，使用多种等级的噪声对数据分布进行扰动，从而使得在低密度区域有更多的训练信号；并提出退火的朗之万动力学采样方法，基于一个噪声条件模型，在不同等级噪声扰动分布上进行采样，最终采样出数据样本。

一、背景知识：Score-based generative modeling

生成模型就是使用数据集来训练一个模型从 $p_{\text{data}}(\mathbf{x})$ 中生成新的样本。假设有从未知的数据分布 $p_{\text{data}}(\mathbf{x})$ 中 iid 采样的数据集 { x i ∈ R D } i = 1 N \{\mathbf{x}_i\in\mathbb{R}^D\}_{i=1}^N {xi​∈RD}i=1N​。我们定义分布 $p(\mathbf{x})$ 的分数 score 为其概率密度对数的梯度，即 ${\nabla }_{\mathbf{x}}\log p(\mathbf{x})$。得分网络 ${\mathbf{s}}_{\theta }:{\mathbb{R}}^D\to {\mathbb{R}}^D$ 是一个参数化网络，用于估计 $p_{\text{data}}(\mathbf{x})$ 的分数。score-based 生成模型有两个关键部分：得分匹配（score matching），和朗之万动力学（Langevin dynamics），分别对应训练过程和采样生成过程。score matching 用来训练一个神经网络估计数据分布的得分，而朗之万动力学则能够根据分布的得分来进行采样。以下分别进行介绍。

1 score matching for score estimation

score matching 最初是用于基于未知分布的样本点来训练非归一化统计模型（non-normalized statistical model）。之前已经开始有研究将其用于得分估计，本文也是这样做。使用 score matching，我们无需知道数据分布 $p_{\text{data}}(\mathbf{x})$ 本身，而是可以直接训练一个得分网络 ${\mathbf{s}}_{\theta }(\mathbf{x})$ 来估计数据分布的得分 ${\nabla }_{\mathbf{x}}\log p_{\text{data}}(\mathbf{x})$。

score matching 的目标函数：
$$\frac{1}{2}{\mathbb{E}}_{p_{\text{data}}(\mathbf{x})}[||{\mathbf{s}}_{\theta }(\mathbf{x})-{\nabla }_{\mathbf{x}}\log p_{\text{data}}(\mathbf{x})|{|}_2^2]$$
可以等价（至多差一个常数）为：
$${\mathbb{E}}_{p_{\text{data}}(\mathbf{x})}[\text{tr}({\nabla }_{\mathbf{x}}{\mathbf{s}}_{\theta }(\mathbf{x}))+\frac{1}{2}||{\mathbf{s}}_{\theta }(\mathbf{x})|{|}_2^2]$$
其中 ${\nabla }_{\mathbf{x}}{\mathbf{s}}_{\theta }(\mathbf{x})$ 表示 ${\mathbf{s}}_{\theta }(\mathbf{x})$ 的雅可比矩阵。已有研究工作显示，在一定条件下，方程（3）的最小化解（记作 ${\mathbf{s}}_{{\theta }^{\ast }}(\mathbf{x})$）也能满足 ${\mathbf{s}}_{{\theta }^{\ast }}(\mathbf{x})={\nabla }_{\mathbf{x}}\log p_{\text{data}}(\mathbf{x})$。实际中，可以使用数据样本快速估计上式对 $p_{\text{data}}(\mathbf{x})$ 的期望，但由于迹 $\text{tr}({\nabla }_{\mathbf{x}}{\mathbf{s}}_{\theta }(\mathbf{x}))$ 的计算复杂度太高，因此得分匹配并不适合直接用深层网络和高维数据训练。下面我们将讨论两种常见的大规模得分匹配的计算方法。

denoising score matching

denoising score matching 是 score matching 的一种变体，它完全绕过了对迹的计算。denoising score matching 中，首先使用预定义的噪声分布 $q_{\sigma }(\tilde{\mathbf{x}}|\mathbf{x})$ 对数据点 $\mathbf{x}$ 进行扰动，然后使用 score matching 估计扰动后的分布，即 $q_{\sigma }(\mathbf{x})\stackrel{△}{=}\int q_{\sigma }(\tilde{\mathbf{x}}|\mathbf{x})p_{\text{data}}(\mathbf{x})d\mathbf{x}$ ，的分数。这个目标函数为：
$$\frac{1}{2}{\mathbb{E}}_{q_{\sigma }(\tilde{\mathbf{x}}|\mathbf{x})p_{\text{data}}(\mathbf{x})}[||{\mathbf{s}}_{\theta }(\tilde{\mathbf{x}})-{\nabla }_{\mathbf{x}}\log q_{\sigma }(\tilde{\mathbf{x}}|\mathbf{x})|{|}_2^2)]$$
已有研究证明，最小化上式的最优解 ${\mathbf{s}}_{{\theta }^{\ast }}(\mathbf{x})$ 也能满足 ${\mathbf{s}}_{{\theta }^{\ast }}(\mathbf{x})={\nabla }_{\mathbf{x}}\log q_{\sigma }(\mathbf{x})$。但是注意，只有在噪声很小的时候，扰动后的分布才约等于数据分布。在生成时，当 $q_{\sigma }(\mathbf{x})\approx p_{\text{data}}(\mathbf{x})$ 才能使用噪声扰动分布的得分近似当做真实数据分布的得分。

siliced score matching

sliced score matching 使用一个随机映射来对 score matching 中迹 $\text{tr}({\nabla }_{\mathbf{x}}{\mathbf{s}}_{\theta }(\mathbf{x}))$ 的计算进行近似。目标函数是：
$${\mathbb{E}}_{p_{\mathbf{v}}}{\mathbb{E}}_{p_{\text{data}}}[{\mathbf{v}}^T{\nabla }_{\mathbf{x}}{\mathbf{s}}_{\theta }(\mathbf{x})\mathbf{v}+\frac{1}{2}||{\mathbf{s}}_{\theta }(\mathbf{x})|{|}_2^2]$$
其中 $p_{\mathbf{v}}$ 是一个简单的随机向量分布（如多元高斯分布）。${\mathbf{v}}^T{\nabla }_{\mathbf{x}}{\mathbf{s}}_{\theta }(\mathbf{x})\mathbf{v}$ 这一项可以由自动微分框架的前向模式进行高效地计算。不同于 denoising score matching 估计的是扰动后分布的分数，sliced score matching 估计的是未扰动的原始数据的分布，但这会带来大概四倍的额外计算量。

2 sampling with Langevin dynamics

朗之万动力学可以在仅知道分布的得分函数 ${\nabla }_{\mathbf{x}}\log p(\mathbf{x})$ 的情况下，采样出分布 $p(\mathbf{x})$ 的样本点。具体来说，给定固定的步长 $ϵ>0$ 和一个先验分布中采样的初始值 $\tilde{\mathbf{x}}\sim \pi (x)$ ，朗之万迭代进行如下计算：
$${\tilde{\mathbf{x}}}_t={\tilde{\mathbf{x}}}_{t-1}+\frac{ϵ}{2}{\nabla }_{\mathbf{x}}\log p({\tilde{\mathbf{x}}}_{t-1})+\sqrt{ϵ}{\mathbf{z}}_t$$
其中 ${\mathbf{z}}_t\sim \mathcal{N}(0,\mathbf{I})$。当 $ϵ\to 0,T\to \infty$ 时，${\tilde{\mathbf{x}}}_T$ 的分布等于 $p(\mathbf{x})$，此时（在一定条件下） ${\tilde{\mathbf{x}}}_T$ 就相当于是从分布 $p(\mathbf{x})$ 中采样出的一个样本点了。严格来说，在中间步，即 $ϵ>0,T<\infty$ 时，需要使用 Metropolis-Hastings 更新来修正上式的误差。但在实际中，这个误差经常是忽略的，本文也假设在 $ϵ$ 足够小，$T$ 足够大的条件下，该误差可以忽略。

在上述采样公式中，可以看到，朗之万动力学只需要知道得分函数 ${\nabla }_{\mathbf{x}}\log p(\mathbf{x})$ 就可以进行采样生成了。因此，要采样得到数据分布 $p_{\text{data}}(\mathbf{x})$ 的样本，我们只需有一个得分网络能估计数据分布的得分 ${\mathbf{s}}_{\theta }(\mathbf{x})\approx {\nabla }_{\mathbf{x}}\log p_{\text{data}}(\mathbf{x})$，然后就可以根据朗之万动力学，使用 ${\mathbf{s}}_{\theta }(\mathbf{x})$ 进行采样就行了。这就是基于得分的生成模型（score-based generative modeling）的核心思路。

二、面临挑战：challenges of score-based model generative modeling

以上介绍了 score-based generative modeling 的基本思想与做法，但实际中，这种生成模型也存在着几个问题，导致实际中应用还不广泛。本节将介绍并分析 score-based generative modeling 所面临的难点。

1 the manifold hypothesis

流形假设认为，现实世界中的数据大都集中在嵌入在高维空间中的低维流形上（换句话说，尽管数据可能存在于一个高维空间，但数据点实际上只在一个低维的子空间（流形）上）。该假设是流形学习的基石，在许多数据集上都能得到验证。在流形假设下，score-based generative models 面临着两个问题。首先，得分 ${\nabla }_{\mathbf{x}}\log p_{\text{data}}(\mathbf{x})$ 是高维空间分布的数据梯度，而当数据点 $\mathbf{x}$​ 只分布在低维流形上时，该梯度是未定义的；第二，score matching 的目标函数只有在数据分布充满整个空间中时，得分估计器才是准确的，当数据点只分布在低维流形上时，是不准确的。

数据分布在低维流形上，这对分数估计的影响可以在下图 sliced score matching 在 CIFAR10 训练过程的 ssm loss 曲线中看出。左侧是直接在 CIFAR10 原图上进行训练，可以看到，loss 曲线明显是震荡的，非常不稳定，难以收敛；而右侧是在原图上叠加了很小幅度 $\mathcal{N}(0,0.0001)$ 的高斯噪声，这个噪声幅度人眼几乎看不出来，但是却显著稳定了训练进行，loss 曲线可以收敛。

2 low data density regions

数据点在低密度区域的分布太少，会对 score-based generative models 中 score matching 训练出的噪声估计器的精确度和朗之万动力学的 MCMC 采样过程噪声带来困难。以下分别探讨数据低密度区域对 score matching 训练和朗之万动力学采样造成的影响。

2.1 inaccurate score estimation with score matching

在低数据密度区域，由于数据稀缺，score mathcing 没有足够的参考依据来准确地估计得分。score matching 的目标是最小化得分估计的平均误差，即 $\frac{1}{2}{\mathbb{E}}_{p_{\text{data}}}[||{\mathbf{s}}_{\theta }(\mathbf{x})-{\nabla }_{\mathbf{x}}\log p_{\text{data}}(\mathbf{x})|{|}_2^2]$。在实际中，关于数据分布的期望是通过对数据分布的独立同分布采样 { x i } i = 1 N ∼ i.i.d p data ( x ) \{\mathbf{x}_i\}_{i=1}^N\stackrel{\text{i.i.d}}{\sim}p_\text{data}(\mathbf{x}) {xi​}i=1N​∼i.i.dpdata​(x)。对任意区域 $\mathcal{R}\subset {\mathbb{R}}^D$，有 $p_{\text{data}}(\mathcal{R})\approx 0$。在大多数情况下，有 { x i } i = 1 N ∩ R = ∅ \{\mathbf{x}_{i}\}_{i=1}^N\cap \mathcal{R}=\empty {xi​}i=1N​∩R=∅，对于该区域中 score matching 没有足够的样本来准确地估计得分。

下图通过一个简易的实验验证了这一点，该实验中，作者用 sliced score matching 来拟合一个简单的数据分布 $p_{\text{data}}=\frac{1}{5}\mathcal{N}((-5,5),\mathbf{I})+\frac{4}{5}\mathcal{N}((5,5),\mathbf{I})$，下图左侧是数据分布真实的得分 ${\nabla }_{\mathbf{x}}\log p_{\text{data}}(\mathbf{x})$，右侧是模型估计的得分 ${\mathbf{s}}_{\theta }$。可以看到（看箭头），只有在数据密集的区域，得分估计才是准确的。

2.2 slow mixing of Langevin dynamics

当数据分布的两个 pattern 被低密度区域分隔开时，朗之万动力学对两个 pattern 的相对权重的预测很慢，从而很难使得采样结果收敛到真实分布。

考虑一个混合分布 $p_{\text{data}}(\mathbf{x})=\pi p_1(\mathbf{x})+(1-\pi )p_2(\mathbf{x})$，其中 $p_1(\mathbf{x})$ 和 $p_2(\mathbf{x})$ 是分布区域不相交的归一化分布，$\pi \in (0,1)$。在 $p_1(\mathbf{x})$ 的分布区域内，有 $\nabla_x \log p_{\text{data}}(x) = \nabla_x(\log \pi + \log p_1(x)) = \nabla_x \log p_1(x) $，同理在 $p_2(\mathbf{x})$ 的分布区域内，有 ${\nabla }_x\log p_{\text{data}}(\mathbf{x})={\nabla }_x(\log (1-\pi )+\log p_2(\mathbf{x}))={\nabla }_x\log p_2(\mathbf{x})$。也就是说，在这两种情况下，得分 ${\nabla }_x\log p_{\text{data}}(\mathbf{x})$ 都与 $\pi$ 无关。朗之万动力学使用 ${\nabla }_x\log p_{\text{data}}(\mathbf{x})$ 来从 $p_{\text{data}}(\mathbf{x})$ 中采样，所得到的样本将不依赖于权重系数 $\pi$。在实际中，虽然不同的 pattern 可能存在一定的相交分布区域，这部分相交区域的数据密度很低，因此结论是类似的，在这种情况下，理论上朗之万动力学还是可以生成正确的样本，但可能需要非常小的步长、非常大的步数。

为了验证这一分析，作者使用了与上节相同的简单数据分布进行了朗之万动力学采样实验，结果如下图。图中 (a)、(b) 分别是对数据分布直接采样、根据数据分布的真实 score 进行朗之万动力学采样（注意这里的实验设了一个简单的数据分布，也就是说数据分布的确切解析式我们是已知的），显然图 (b) 朗之万动力学得到的样本之间两种 pattern 的相对密度是不正确的，与上述分析相符。

三、提出方法：Noise Conditional Score Networks: learning and inference

如何解决 score-based generative models 的上述问题呢？

作者发现，使用随机高斯噪声对数据分布进行扰动，会使得数据分布更适合于 score-based generative modeling。首先，由于高斯分布遍布整个空间，因此使用其对数据分布进行扰动，分布就不会仅局限在低维流形上了，这使得数据分布的得分在整个空间上都是良定义的，消除了流形假设带来的难点。第二，大的高斯噪声可以填满原本未扰动数据分布的低密度区域，从而 score mathcing 的训练过程就更多的监督信号，模型估计的得分就更准。第三，使用多种噪声等级，我们可以获得一系列噪声扰动分布，这一系列噪声扰动分布逐渐逼近真实数据分布。通过使用模拟退火（simulated annealing）和退火重要性采样（annealed importance sampling）的思想，我们可以利用这些中间分布来提高多模态分布上朗之万动力学采样的速率。

基于以上思考，作者提出了对 score-based generative modeling 的改进措施：1）使用多种不同等级的噪声对数据分布进行扰动；2）训练单个（条件）神经网络来同时估计所有噪声等级所对应分布的得分。训练结束后，使用朗之万动力学来进行采样生成时，我们从最大噪声所对应分布的得分开始，逐步退火降低噪声等级，这有助于权衡大小噪声水平的益处，平稳地将大噪声水平的益处逐渐转移到低噪声水平，直到扰动后的数据几乎与原始数据无法区分。下面，我们将更详细地阐述我们方法的细节，包括得分网络的架构、训练目标和退火的朗之万动力学采样。

1 网络：noise conditional score networks

记 { σ i } i = 1 L \{\sigma_i\}_{i=1}^L {σi​}i=1L​ 是一个正的等比序列，满足：$\frac{{\sigma }_1}{{\sigma }_2}=\cdots =\frac{{\sigma }_{L-1}}{{\sigma }_L}>1$，它是我们的噪声等级序列。记 $q_{\sigma }(\mathbf{x})\stackrel{△}{=}\int p_{\text{data}}(\mathbf{t})\mathcal{N}(\mathbf{x}|\mathbf{t},{\sigma }^2\mathbf{I})d\mathbf{t}$ 表示经噪声扰动的数据分布。在我们选定的噪声等级序列中，要满足 1) ${\sigma }_1$ 足够大，能够消除流形假设中提到的低密度区域训练信号不足问题；2) ${\sigma }_L$ 足够小，在最后减少对真实数据分布的影响。我们要训练一个条件得分网络，来按条件估计不同噪声扰动等级下分布的得分，即 ∀ σ ∈ { σ i } i = 1 L , s θ ( x , σ ) ≈ ∇ x log ⁡ q σ ( x ) \forall\sigma\in\{\sigma_i\}_{i=1}^L,\ \ \mathbf{s_\theta}(\mathbf{x},\sigma)\approx\nabla_\mathbf{x}\log q_\sigma(\mathbf{x}) ∀σ∈{σi​}i=1L​,  sθ​(x,σ)≈∇x​logqσ​(x)。注意我们的模型输入输出尺寸是相同的，即当 $\mathbf{x}\in {\mathbb{R}}^D$ 时，${\mathbf{s}}_{\theta }(\mathbf{x},\sigma )\in {\mathbb{R}}^D$。作者将网络 ${\mathbf{s}}_{\theta }(\mathbf{x},\sigma )$ 命名为噪声条件得分网络（Noise Conditional Score Networks，NCSN）。

具体网络结构方面，由于 NCSN 的输入输出尺寸相同，作者借鉴密集预测类任务（如图像分割），选择了 UNet 作为主体网络结构，加入了空洞卷积、实例归一化等技术，并使用一种改进的条件实例归一化（conditional instance normalization）来将条件 $\sigma$ 注入到模型中。

2 训练：learning NCSNs via score matching

作者强调，之前介绍的 sliced score matching 和 denoising score matching 都可以用来训练 NCSN。本文中考虑到 dsm 训练更快且更适合我们估计噪声扰动数据分布的得分的任务，故而采用 dsm 进行训练。本文选用噪声分布为
$$q_{\sigma }(\tilde{\mathbf{x}}|\mathbf{x})=\mathcal{N}(\tilde{\mathbf{x}}|\mathbf{x},{\sigma }^2\mathbf{I})$$
从而有：
KaTeX parse error: Got function '\tilde' with no arguments as subscript at position 9: \nabla_\̲t̲i̲l̲d̲e̲{\mathbf{x}}\lo…
对于给定的 $\sigma$，denosing score matching 的目标函数为：
$$l(\theta ;\sigma )\stackrel{△}{=}\frac{1}{2}{\mathbb{E}}_{p_{\text{data}}(\mathbf{x})}{\mathbb{E}}_{\tilde{\mathbf{x}}\sim \mathcal{N}(\mathbf{x},{\sigma }^2\mathbf{I})}[||{\mathbf{s}}_{\theta }(\tilde{\mathbf{x}},\sigma )+\frac{\tilde{\mathbf{x}}-\mathbf{x}}{{\sigma }^2}|{|}_2^2]$$
对于所有的 σ ∈ { σ i } i = 1 L \sigma\in\{\sigma_i\}_{i=1}^L σ∈{σi​}i=1L​ 结合起来，统一的目标函数为：
L ( θ ; { σ i } i = 1 L ) = △ 1 L ∑ i = 1 L λ ( σ i ) l ( θ ; σ i ) \mathcal{L}(\theta;\{\sigma_i\}_{i=1}^L)\stackrel{\triangle}{=}\frac{1}{L}\sum_{i=1}^L\lambda(\sigma_i)l(\theta;\sigma_i) L(θ;{σi​}i=1L​)=△L1​i=1∑L​λ(σi​)l(θ;σi​)
其中 $\lambda ({\sigma }_i)>0$ 是一个关于 ${\sigma }_i$ 的加权函数。假设 ${\mathbf{s}}_{\theta }(\mathbf{x},\sigma )$ 训练得足够好，对所有的 $i\in [1,L]$ 都能最小化上述目标函数，那么模型就能准确地估计噪声扰动数据分布的得分，即 ${\mathbf{s}}_{{\theta }^{\ast }}(\mathbf{x},\sigma )={\nabla }_{\mathbf{x}}\log q_{{\sigma }_i}(\mathbf{x})$。

对于 $\lambda (\cdot )$ 的选择，可以有很多种。理想情况下，我们希望对所有的 ${\sigma }_i$，损失 $\lambda ({\sigma }_i)l(\theta ;{\sigma }_i)$ 的值大致在相同的数量级上。实验中，作者发现当网络训练得比较好时，大致有 $||{\mathbf{s}}_{\theta }(\mathbf{x},\sigma )|{|}_2\propto 1/\sigma$。作者据此选择 $\lambda (\sigma )={\sigma }^2$。这是因为在此选择下，有 $\lambda (\sigma )l(\theta ;\sigma )={\sigma }^{2}l((\theta ;\sigma ))=\frac{1}{2}\mathbb{E}[||\sigma {\mathbf{s}}_{\theta }(\mathbf{x};\sigma )+\frac{\tilde{\mathbf{x}}-\mathbf{x}}{\sigma }|{|}_2^2]$ ，而由于 $\frac{\tilde{\mathbf{x}}-\mathbf{x}}{\sigma }\sim \mathcal{N}(0,\mathbf{I})$、$||\sigma {\mathbf{s}}_{\theta }(\mathbf{x};\sigma )|{|}_2\propto 1$，因此可得此时 $\lambda (\sigma )l(\theta ;\sigma )$ 值的数量级与 $\sigma$ 无关。

3 采样：NCSN inference via annealed Langevin dynamics

在 NCSN ${\mathbf{s}}_{\theta }(\mathbf{x},\sigma )$ 训练完成后，我们就能够对不同噪声扰动数据分布 $q_{{\sigma }_i}(\mathbf{x})$ 的得分 ${\nabla }_{\mathbf{x}}{q}_{{\sigma }_1}(\mathbf{x})$ 进行较为准确的估计。这样，就能执行朗之万动力学采样，从各个分布 $q_{{\sigma }_i}(\mathbf{x})$ 中采样出样本。但是，如前所述，流形假设指出了朗之万动力学采样对多 pattern 数据分布采样时混合较慢的问题，并且我们的 NCSN 需要估计的是多个噪声扰动分布 $q_{{\sigma }_i}(\mathbf{x})$ 的得分，而不是数据分布 $p_{\text{data}}(\mathbf{x})$ 的得分。因此，本文提出了一种新的采用方法——退火的朗之万动力学（Annealed Langevin Dynamics），来使用 NCSN 进行采样生成。

具体算法如下图所示。我们从一个固定的先验分布（如均匀分布）中采样一个初始样本。首先对第一级的噪声扰动数据分布 $q_{{\sigma }_1}(\mathbf{x})$ 执行 $T$ 步朗之万动力学采样，步长为 ${\alpha }_1$。然后将对 $q_{{\sigma }_1}(\mathbf{x})$ 的最终采样结果作为初始样本，再对 $q_{{\sigma }_2}(\mathbf{x})$ 进行朗之万动力学采样，步长为 ${\alpha }_2$，注意 ${\alpha }_2<{\alpha }_1$。迭代进行这种方式，每次使用上一级分布 $q_{{\sigma }_{i-1}}(\mathbf{x})$ 朗之万动力学采样迭代的最后一个样本，作为下一级分布 $q_{{\sigma }_i}(\mathbf{x})$ 采样的起始样本，并逐渐调小步长 ${\alpha }_i=ϵ\cdot {\sigma }_i^2/{\sigma }_L^2$。最终，我们对最后一级噪声扰动分布 $q_{{\sigma }_L}(\mathbf{x})$ 执行朗之万动力学采样，得到最终的生成结果。由于 ${\sigma }_L\approx 0$ 足够小，因此 $q_{{\sigma }_L}(\mathbf{x})$ 已经非常接近真实的数据分布 $p_{\text{data}}(\mathbf{x})$。

由于噪声扰动分布 { q σ i } i = 1 L \{q_{\sigma_i}\}_{i=1}^L {qσi​​}i=1L​ 都是由高斯噪声进行扰动，因此这些分布的存在区域都是充满整个空间的，从而他们的得分都是良定义的，避免了前文所述的流形假设指出的分布得分未定义的问题。当 ${\sigma }_1$ 足够大时，分布 $q_{{\sigma }_1}(\mathbf{x})$ 的低密度区域就很小，并且不同 pattern 之间也没有那么分离。从而，得分估计就更准，并且朗之万动力学采样对不同 pattern 的混合也更快。因此，朗之万动力学从 $q_{{\sigma }_1}(\mathbf{x})$ 采样出的结果就是比较好的样本。这些样本很大概率是来自 $q_{{\sigma }_1}(\mathbf{x})$ 的高密度区域，而由于 $q_{{\sigma }_2}(\mathbf{x})$ 与 $q_{{\sigma }_1}(\mathbf{x})$ 非常接近，从而这些样本同样很大概率落在 $q_{{\sigma }_2}(\mathbf{x})$ 的高密度区域。而得分估计和朗之万动力学采样都是在分布的高密度区域工作得更好，也就是说第一级噪声扰动分布 $q_{{\sigma }_1}(\mathbf{x})$ 中采样出来的最终样本为第二级噪声扰动分布 $q_{{\sigma }_2}(\mathbf{x})$ 的朗之万动力学采样提供了很好的初始化。以此类推，每一级分布 $q_{{\sigma }_{i-1}}(\mathbf{x})$ 的样本都为下一级分布 $q_{{\sigma }_i}(\mathbf{x})$ 的采样提供了很好的初始化。最终我们能从分布 $q_{{\sigma }_L}(\mathbf{x})$ 中采样到高质量的结果。

步长 ${\alpha }_i$ 的设置也有很多可能。本文的选择是 ${\alpha }_i\propto {\sigma }_i^2$。这样设计的动机是想要固定朗之万动力学采样中信噪比 (Signal-to-Noise Ratio, SNR) $\frac{{\alpha }_i{\mathbf{s}}_{\theta }(\mathbf{x},{\sigma }_{\mathbf{i}})}{2\sqrt{{\alpha }_i}\mathbf{z}}$ 的尺度。注意有 $\mathbb{E}[||\frac{{\alpha }_i{\mathbf{s}}_{\theta }(\mathbf{x},{\sigma }_{\mathbf{i}})}{2\sqrt{{\alpha }_i}\mathbf{z}}|{|}_2^2]\approx \mathbb{E}[\frac{{\alpha }_i||{\mathbf{s}}_{\theta }(\mathbf{x},{\sigma }_i)|{|}_2^2}{4}]\propto \frac{1}{4}\mathbb{E}[||{\sigma }_i{\mathbf{s}}_{\theta }(\mathbf{x},{\sigma }_i)|{|}_2^2]$。我们之前提到，当 NCSN 训练得比较好的时候，有 $||{\mathbf{s}}_{\theta }(\mathbf{x},\sigma )|{|}_2\propto 1/\sigma$，此时有 $\mathbb{E}[||{\sigma }_i{\mathbf{s}}_{\theta }(\mathbf{x},{\sigma }_i)|{|}_2^2\propto 1$，从而信噪比 $\mathbb{E}[||\frac{{\alpha }_i{\mathbf{s}}_{\theta }(\mathbf{x},{\sigma }_{\mathbf{i}})}{2\sqrt{{\alpha }_i}\mathbf{z}}|{|}_2^2]\propto \frac{1}{4}\mathbb{E}[||{\sigma }_i{\mathbf{s}}_{\theta }(\mathbf{x},{\sigma }_i)|{|}_2^2]\propto \frac{1}{4}$，与 ${\sigma }_i$ 无关。

为了展示退火朗之万动力学的有效性，本文在前面介绍标准朗之万动力学混合慢的简单实验中，对比了本文所提方法。回顾一下，实验的目标是使用得分从一个简单的两个高斯混合的数据分布中进行采样。退火的朗之万动力学的采样结果如图 © 所示，明显可以看到，采样出的数据能够反映出两个高斯分布的比重差异，也就是说能够更真实地从数据分布中采样，而标准的朗之万动力学采样（图 (b）显然没有做到这一点。

结语

这篇 NCSN 的方法提出在 DDPM 之前，并且对 DDPM 的设计也有启发。更重要的是，笔者认为，基于得分的生成模型似乎更加直觉，更具有可解释性。从理解的角度来说，NSCN 更加顺畅。但是 DDPM 的效果似乎更好？因此一般认为 DDPM 是真正将 Diffusion 这一类方法做 work 的工作。宋飏博士在之后的 score based SDE 中也已经将 NCSN 和 DDPM 统一到了 SDE 的框架下。