【深度学习】扩散模型(Diffusion Model)详解

文章目录

  • 【深度学习】扩散模型(Diffusion Model)详解
    • 1. 介绍
    • 2. 具体方法
      • 2.1 扩散过程
      • 2.2 逆扩散过程
      • 2.3 损失函数
    • 3. 总结
    • 4. 参考

1. 介绍

扩散模型有两个过程:

  • 扩散过程:如上图所示,扩散过程为从右到左 X0→ XT X_0 \rightarrow X_TX0XT 的过程,表示对图片逐渐加噪,且 X t + 1 X_{t+1}Xt+1是在 Xt X_{t}Xt上加躁得到的,其只受 Xt X_{t}Xt的影响。因此扩散过程是一个马尔科夫过程

    • X 0X_0 X0表示从真实数据集中采样得到的一张图片,对 X 0X_0 X0添加 TT T 次噪声,图片逐渐变得模糊。当 TT T 足够大时, X TX_T XT为标准正态分布。在训练过程中,每次添加的噪声是已知的,即 q( X t∣ X t−1 )q(X_t|X_{t-1}) q(XtXt1) 是已知的。根据马尔科夫过程的性质,我们可以递归得到 q( X t∣ X 0)q(X_t|X_0) q(XtX0),即 q( X t∣ X 0)q(X_t|X_0) q(XtX0) 是已知的。

    其中,扩散过程最主要的是 q( X t∣ X 0)q(X_t|X_0) q(XtX0)q( X t∣ X t−1 )q(X_t|X_{t-1}) q(XtXt1)的推导。

  • 逆扩散过程:如上图所示,逆扩散过程为从左到右 XT→ X0 X_T \rightarrow X_0XTX0 的过程,表示从噪声中逐渐复原出图片。如果我们能够在给定 Xt X_tXt 条件下知道 X t − 1 X_{t-1}Xt1 的分布,即如果我们可以知道 q ( X t − 1∣ Xt)q(X_{t-1}|X_t)q(Xt1Xt),那我们就能够从任意一张噪声图片中经过一次次的采样得到一张图片而达成图片生成的目的。

    • 显然我们很难知道 q( X t−1 ∣ X t)q(X_{t-1}|X_t) q(Xt1Xt),因此我们才会用 p Θ( X t−1 ∣ X t)p_{Θ}(X_{t-1}|X_t) pΘ(Xt1Xt) 来近似 q( X t−1 ∣ X t)q(X_{t-1}|X_t) q(Xt1Xt),而 p Θ( X t−1 ∣ X t)p_{Θ}(X_{t-1}|X_t) pΘ(Xt1Xt) 就是我们要训练的网络,在原文中就是个U-Net。而很妙的是,虽然我们不知道 q( X t−1 ∣ X t)q(X_{t-1}|X_t) q(Xt1Xt),但是 q( X t−1 ∣ X t X 0)q(X_{t-1}|X_tX_0) q(Xt1XtX0) 却是可以用 q( X t∣ X 0)q(X_t|X_0) q(XtX0)q( X t∣ X t−1 )q(X_t|X_{t-1}) q(XtXt1) 表示的,即 q( X t−1 ∣ X t X 0)q(X_{t-1}|X_tX_0) q(Xt1XtX0) 是可知的,因此我们可以用 q( X t−1 ∣ X t X 0)q(X_{t-1}|X_tX_0) q(Xt1XtX0) 来指导 p Θ( X t−1 ∣ X t)p_{Θ}(X_{t-1}|X_t) pΘ(Xt1Xt) 进行训练。

    其中,逆扩散过程最主要的是 q( X t−1 ∣ X t X 0)q(X_{t-1}|X_tX_0) q(Xt1XtX0)的推导。

2. 具体方法

在上面的介绍中,我们已经明确了要训练 pΘ( X t − 1∣ Xt)p_{Θ}(X_{t-1}|X_t)pΘ(Xt1Xt),但是目标函数如何确定?

有两个很直接的想法:

  • 负对数的最大似然概率,即 −log p Θ( X 0)-logp_{Θ}(X_0) logpΘ(X0)
  • 真实分布与预测分布的交叉熵,即 − E q( X 0) log p Θ( X 0)-E_{q(X_0)}logp_{Θ}(X_0) Eq(X0)logpΘ(X0)

但是这两种方法都很难去求解(求积分)和优化。因此扩散模型参考了VAE,不去优化这两个东西,而是优化他们的变分上界(variational lower bound),定义 L V L B L_{VLB}LVLB,如下:

L V L B L_{VLB}LVLB 减小就代表着 − l o g pΘ( X0)-logp_{Θ}(X_0)logpΘ(X0) − E q ( X0)l o g pΘ( X0)-E_{q(X_0)}logp_{Θ}(X_0)Eq(X0)logpΘ(X0) 的上界减小。且经过推导,$L_{VLB}又可写成如下形式:

由上式不难发现, Lt L_{t}Lt就是逆扩散过程中 q ( Xt∣ X t + 1X0)q(X_{t}|X_{t+1}X_0)q(XtXt+1X0) pΘ( Xt∣ X t + 1)p_{Θ}(X_{t}|X_{t+1})pΘ(XtXt+1) K LKLKL 散度,即用 q ( Xt∣ X t + 1X0)q(X_{t}|X_{t+1}X_0)q(XtXt+1X0) 来指导 pΘ( Xt∣ X t + 1)p_{Θ}(X_{t}|X_{t+1})pΘ(XtXt+1) 进行训练。这部分主要就是(1)式和(2)式的推导,细节部分见下文的损失函数。

2.1 扩散过程

如上图所示,扩散过程为从右到左 X0→ XT X_0 \rightarrow X_TX0XT 的过程,表示对图片逐渐加噪。

  • X t+1 X_{t+1} Xt+1是在 X tX_{t} Xt上加躁得到的,其只受 X tX_{t} Xt的影响。因此扩散过程是一个马尔科夫过程

下面,我们对扩散过程进行推导:

由于每一步扩散的步长受变量 { βt∈ ( 0 , 1 ) } t = 1T \{β_{t} \in (0,1)\}_{t=1}^{T}{βt(0,1)}t=1T 的影响。 q ( Xt∣ X t − 1)q(X_{t}|X_{t-1})q(XtXt1) 可写为如下形式,即给定 X t − 1 X_{t-1}Xt1 的条件下, Xt X_{t}Xt服从均值为 1 − βt X t − 1 \sqrt{1-β_{t}}X_{t-1}1βt Xt1,方差为 βtIβ_{t}IβtI的正态分布:

用重参数化技巧表示 Xt X_{t}Xt,令 αt= 1 − βt α_{t}=1-β_{t}αt=1βt ,令 Zt∼ N ( 0 , I ) , t ≥ 0Z_{t} \sim N(0,I), t \ge 0ZtN(0,I),t0,即:

其中,

可以得到,令 α ˉt=∏ i=1t α i \bar{α}_{t}= {\textstyle \prod_{i=1}^{t}α_{i}}αˉt=i=1tαi

设随机变量 Z ˉ t − 1 \bar{Z}_{t-1}Zˉt1 为:

Z ˉ t − 1 \bar{Z}_{t-1}Zˉt1 的期望和方差如下:

因此,

至此,我们推出了q( X t∣ X t−1 )q(X_{t}|X_{t-1}) q(XtXt1)q( X t∣ X 0)q(X_{t}|X_{0}) q(XtX0),完成了扩散过程。

2.2 逆扩散过程

如上图所示,逆扩散过程为从左到右 XT→ X0 X_T \rightarrow X_0XTX0 的过程,表示从噪声中逐渐复原出图片。

  • 如果我们能够在给定 X tX_t Xt 条件下知道 X t−1 X_{t-1} Xt1 的分布,即如果我们可以知道 q( X t−1 ∣ X t)q(X_{t-1}|X_t) q(Xt1Xt),那我们就能够从任意一张噪声图片中经过一次次的采样得到一张图片而达成图片生成的目的。
  • 显然我们很难知道 q( X t−1 ∣ X t)q(X_{t-1}|X_t) q(Xt1Xt),因此我们才会用 p Θ( X t−1 ∣ X t)p_{Θ}(X_{t-1}|X_t) pΘ(Xt1Xt) 来近似 q( X t−1 ∣ X t)q(X_{t-1}|X_t) q(Xt1Xt),而 p Θ( X t−1 ∣ X t)p_{Θ}(X_{t-1}|X_t) pΘ(Xt1Xt) 就是我们要训练的网络。

虽然我们不知道 q ( X t − 1∣ Xt)q(X_{t-1}|X_t)q(Xt1Xt),但是 q ( X t − 1∣ XtX0)q(X_{t-1}|X_tX_0)q(Xt1XtX0) 却是可以用 q ( Xt∣ X0)q(X_t|X_0)q(XtX0) q ( Xt∣ X t − 1)q(X_t|X_{t-1})q(XtXt1) 表示的,即 q ( X t − 1∣ XtX0)q(X_{t-1}|X_tX_0)q(Xt1XtX0) 是可知的。

  • 因此我们可以用 q( X t−1 ∣ X t X 0)q(X_{t-1}|X_tX_0) q(Xt1XtX0) 来指导 p Θ( X t−1 ∣ X t)p_{Θ}(X_{t-1}|X_t) pΘ(Xt1Xt) 进行训练。

下面我们对逆扩散过程进行推导:

先对 q ( X t − 1∣ XtX0)q(X_{t-1}|X_tX_0)q(Xt1XtX0)进行推导:

现在,我们已经把 q ( X t − 1∣ XtX0)q(X_{t-1}|X_tX_0)q(Xt1XtX0) q ( Xt∣ X0)q(X_t|X_0)q(XtX0) q ( Xt∣ X t − 1)q(X_t|X_{t-1})q(XtXt1) 进行表示,下面对 q ( X t − 1∣ XtX0)q(X_{t-1}|X_tX_0)q(Xt1XtX0) 的表达式进行推导:


至此,我们得到了 q( X t−1 ∣ X t X 0)q(X_{t-1}|X_tX_0) q(Xt1XtX0)的分布表达式,完成了逆扩散过程。

2.3 损失函数

我们已经明确了要训练 pΘ( X t − 1∣ Xt)p_{Θ}(X_{t-1}|X_t)pΘ(Xt1Xt),那要怎样进行训练?有两个很直接的想法:

  • 一个是负对数的最大似然概率,即 −log p Θ( X 0)-logp_{Θ}(X_0) logpΘ(X0)
  • 另一个是真实分布与预测分布的交叉熵,即 − E q( X 0) log p Θ( X 0)-E_{q(X_0)}logp_{Θ}(X_0) Eq(X0)logpΘ(X0)

然而,类似于VAE,由于我们很难对噪声空间进行积分,因此直接优化 − l o g pΘ -logp_{Θ}logpΘ E q ( X0)l o g pΘ( X0)E_{q(X_0)}logp_{Θ}(X_0)Eq(X0)logpΘ(X0)都是很困难的。

因此我们不直接优化它们,而是优化它们的变分上界 L V L B L_{VLB}LVLB,其定义如下:

下面证明 L V L B L_{VLB}LVLB − l o g pΘ( X0)-logp_{Θ}(X_0)logpΘ(X0) − E q ( X0)l o g pΘ( X0)-E_{q(X_0)}logp_{Θ}(X_0)Eq(X0)logpΘ(X0) 的上界,即证明 L V L B≥ − l o g pΘ( X0)L_{VLB} \ge -logp_{Θ}(X_0)LVLBlogpΘ(X0) L V L B≥ − E q ( X0)l o g pΘ( X0)L_{VLB} \ge -E_{q(X_0)}logp_{Θ}(X_0)LVLBEq(X0)logpΘ(X0)


至此,证明了 L V L B L_{VLB}LVLB − l o g pΘ( X0)-logp_{Θ}(X_0)logpΘ(X0) − E q ( X0)l o g pΘ( X0)-E_{q(X_0)}logp_{Θ}(X_0)Eq(X0)logpΘ(X0)的上界。进而,对 L V L B L_{VLB}LVLB进行化简,得到:

Lt L_{t}Lt 即可看出,对 pΘ( Xt∣ X t + 1)p_{Θ}(X_{t}|X_{t+1})pΘ(XtXt+1) 的监督就是最小化 pΘ( Xt∣ X t + 1)p_{Θ}(X_{t}|X_{t+1})pΘ(XtXt+1) q ( Xt∣ X t + 1X0)q(X_t|X_{t+1}X_0)q(XtXt+1X0) 的KL散度。

3. 总结

总结来说,扩散模型的目的是希望学习出一个 pΘ( X t − 1∣ Xt)p_{Θ}(X_{t-1}|X_t)pΘ(Xt1Xt),即能够从噪声图恢复出原图。
为了达到这一个目的,我们使用 q ( X t − 1∣ XtX0)q(X_{t-1}|X_tX_0)q(Xt1XtX0) 来监督 pΘ( X t − 1∣ Xt)p_{Θ}(X_{t-1}|X_t)pΘ(Xt1Xt) 进行训练,而 q ( X t − 1∣ XtX0)q(X_{t-1}|X_tX_0)q(Xt1XtX0) 是可以用 q ( Xt∣ X0)q(X_t|X_0)q(XtX0) q ( Xt∣ X t − 1)q(X_t|X_{t-1})q(XtXt1) 进行表示的,即 q ( X t − 1∣ XtX0)q(X_{t-1}|X_tX_0)q(Xt1XtX0)是已知的。

4. 参考

【1】https://blog.csdn.net/Little_White_9/article/details/124435560
【2】https://lilianweng.github.io/posts/2021-07-11-diffusion-models/
【3】https://arxiv.org/abs/2105.05233
【4】https://arxiv.org/abs/1503.03585
【5】https://arxiv.org/abs/2006.11239