图像生成发展起源：从VAE、VQ-VAE、扩散模型DDPM、DETR到ViT、Swin transformer 转自 七月_vqvae和transformer

前言
2018年我写过一篇博客，叫：《一文读懂目标检测：R-CNN、Fast R-CNN、Faster R-CNN、YOLO、SSD》，该文相当于梳理了2019年之前CV领域的典型视觉模型，比如

2014 R-CNN
2015 Fast R-CNN、Faster R-CNN
2016 YOLO、SSD
2017 Mask R-CNN、YOLOv2
2018 YOLOv3
随着2019 CenterNet的发布，特别是2020发布的DETR(End-to-End Object Detection with Transformers)之后，自此CV迎来了生成式下的多模态时代(本文介绍其中有下划线的部分，其他部分下一篇介绍)

1月    3月    4月    5月    6月    8月    9月    10月    11月
2020                DETR    DDPM            
DDIM

VisionTransf.. 

2021    
CLIP

DALL·E

SwinTransf..                            
MAE

SwinTransf..V2

2022    BLIP        DALL·E 2            
StableDiffusion 

BEiT-3

Midjourney V3

2023    BLIP2    
VisualChatGPT 

GPT4

Midjourney V5

SAM(Segment Anything Model)                DALLE3        
但看这些模型接二连三的横空出世，都不用说最后爆火的GPT4，便可知不少CV同学被卷的不行

说到GPT4，便不得不提ChatGPT，实在是太火了，改变了很多行业，使得国内外绝大部分公司的产品、服务都值得用LLM全部升级一遍(比如微软的365 Copilot、阿里所有产品、金山WPS等等)

而GPT4相比GPT3.5或GPT3最本质的改进就是增加了多模态的能力，使得ChatGPT很快就能支持图片的输入形式，从而达到图生文和文生图的效果，而AI绘画随着去年stable diffusion和Midjourney的推出，使得文生图火爆异常，各种游戏的角色设计、网上店铺的商品/页面设计都用上了AI绘画这样的工具，更有不少朋友利用AI绘画取得了不少的创收，省时省力还能赚钱，真香

但面对这么香的技术，其背后的一系列原理到底是什么呢，本文特从头开始，不只是简单的讲一下扩散模型的原理，而是在反复研读相关论文之后，准备把20年起相关的CV多模态模型全部梳理一遍，从VE、VAE、DDPM到ViT/Swin transformer、CLIP/BLIP，再到stable diffusion/Midjourney、GPT4，当然，实际写的时候，会分成两篇甚至多篇文章，比如

第一篇，即本文《AI绘画能力的起源：从VAE、扩散模型DDPM、DETR到ViT/MAE/Swin transformer》
第二篇，即下篇《CV多模态和AIGC的原理解析：从CLIP、BLIP到Stable Diffusion、Midjourney》
就当2020年之后的CV视觉发展史了，且过程中会尽可能写透彻每一个模型的原理，举两个例子

网上介绍VAE的文章都太数学化(更怕那种表面正确其实关键的公式是错的误导人)，如果更边推导边分析背后的理论意义(怎么来的 出发点是什么 为什么要这么做 这么做的意义是什么)，则会更好理解，这就跟变介绍原理边coding实现 会更好理解、理解更深 一个道理
如果完全展开DDPM推导的所有细节，假定需要100步的话，本文正在朝展开80步而努力，截止5月份之前，绝大部分的中文资料只展开了60步(正在因为只展开了60%，让很多初学者卡到中途)，所以你害怕的不是公式，你只是怕公式的展开不够细致，毕竟对每一个人而言，公式展开越细致 越不怕
(如果本文有任何一个公式展开的不够细致、不够一目了然，请随时指出，一定及时二次展开)
第一部分 编码器VE与变分自编码器VAE
1.1 AE：编码器(数据X
�
压缩为低维表示z
�
)-解码器(低维表示恢复为原始数据X^
�
^
)架构
自编码器(Autoencoder，简称AE)是一种无监督学习的神经网络，用于学习输入数据的压缩表示。具体而言，可以将其分为两个部分：编码器和解码器

编码器：编码器是一个神经网络，负责将输入数据X
�
（如图像、文本等）压缩为一个低维表示z
�
，且表示为z=g(X)
�
=
�
(
�
)

解码器：解码器是另一个神经网络，负责将编码器生成的低维表示恢复为原始数据X^
�
^
，且表示为X^=f(z)
�
^
=
�
(
�
)

从而最终完成这么一个过程：X→z→X^
�
→
�
→
�
^
，而其训练目标即是最小化输入数据X
�
与解码器重建数据X^
�
^
之间的差异，所以自编码器常用的一个损失函数为l=||X−X^||2
�
=
|
|
�
−
�
^
|
|
2

这个自编码的意义在于

模型训练结束后，我们就可以认为编码z
�
囊括了输入数据X
�
的大部分信息，也因此我们可以直接利用z
�
表达原始数据，从而达到数据降维的目的
解码器只需要输入某些低维向量z
�
，就能够输出高维的图片数据X^
�
^
，那我们能否把解码器模型直接当做生成模型，在低维空间中随机生成某些向量z
�
，再喂给解码器f(z)
�
(
�
)
来生成图片呢？
对于第二点，理论上可以这么做，但问题在于

绝大多数随机生成的z,f(z)
�
,
�
(
�
)
只会生成一些没有意义的噪声，之所以如此，原因在于没有显性的对z
�
的分布p(z)
�
(
�
)
进行建模，我们并不知道哪些z
�
能够生成有用的图片
而且我们用来训练f(z)
�
(
�
)
的数据是有限的，f
�
可能只会对极有限的z
�
有响应。而整个低维空间又是一个比较大的空间，如果只在这个空间上随机采样的话，我们自然不能指望总能恰好采样到能够生成有用的图片的z
�
有问题自然便得探索对应的解决方案，而VAE(自变分编码器，Variational Autoencoders)则是在AE的基础上，显性的对z
�
的分布p(z)
�
(
�
)
进行建模(比如符合某种常见的概率分布)，使得自编码器成为一个合格的生成模型

1.2 Variational AutoEncoder (VAE)：学习一个分布
1.2.1 VAE：标数据的分布X^
�
^
和目标分布X
�
尽量接近
VAE和GAN一样，都是从隐变量Z
�
生成目标数据， 具体如下图所示(本1.2节的部分图来自苏剑林)：

先用某种分布随机生成一组隐变量Z={Z1,Z2,⋯,Zk}
�
=
{
�
1
,
�
2
,
⋯
,
�
�
}
(假设隐变量服从正态分布)
然后这个Z
�
隐变量经过一个生成器生成一组数据X^={X1^,X2^,⋯,Xk^}
�
^
=
{
�
1
^
,
�
2
^
,
⋯
,
�
�
^
}
而VAE和GAN都希望这组生成数据的分布X^
�
^
和目标分布X
�
尽量接近，看似美好，但有两个问题

一方面，“尽量接近”并没有一个确定的关于X^
�
^
 和 X
�
的相似度的评判标准，比如KL散度便不行，原因在于KL散度是针对两个已知的概率分布求相似度的，而 X^
�
^
和 X
�
 的概率分布目前都是未知(只有一批采样数据 没有分布表达式)
二方面，经过采样出来的每一个Zk
�
�
，不一定对应着每一个原来的Xk
�
�
，故最后没法直接最小化D2(Xk,Xk^)
�
2
(
�
�
,
�
�
^
)
实际是怎么做的呢？如苏剑林所说，与自动编码器由编码器与解码器两部分构成相似，VAE利用两个神经网络建立两个概率密度分布模型：

一个用于原始输入数据X={X1,...,Xk}
�
=
{
�
1
,
.
.
.
,
�
�
}
的变分推断，生成隐变量Z
�
的变分概率分布p(Z|Xk)
�
(
�
|
�
�
)
，称为推断网络
而VAE的核心就是，我们不仅假设 p(Z)
�
(
�
)
是正态分布，而且假设每个p(Z|Xk)
�
(
�
|
�
�
)
 也是正态分布。什么意思呢？即针对每个采样点Xk
�
�
获得一个专属于它和 Z
�
 的一个正态分布p(Z|Xk)
�
(
�
|
�
�
)

换言之，有 k
�
 个 X
�
 sample，就有k
�
个正态分布 p(Z|Xk)
�
(
�
|
�
�
)
，毕竟没有任何两个采样点是完全一致的，而后面要训练一个生成器 Xk^=f(Z)
�
�
^
=
�
(
�
)
，希望能够把从分布 p(Z|Xk)
�
(
�
|
�
�
)
采样出来的一个 Zk
�
�
 还原为 Xk
�
�

→
→
  而如果从 p(Z)
�
(
�
)
中采样一个 Zk
�
�
，没法知道这个 Zk
�
�
 对应于真实的 Xk
�
�
 呢？
→
→
  现在 p(Z|Xk)
�
(
�
|
�
�
)
 专属于Xk
�
�
，我们有理由说从这个分布采样出来的 Zk
�
�
 可以还原到对应的Xk
�
�
 中去

而如何确定这 k
�
 个正态分布呢，众所周知，确定一个正太分布只需确定其均值u
�
和方差σ2
�
2
 即可，故可通过已知的Xk
�
�
 和 假设的Z
�
 去确定均值和方差，具体可以构建两个神经网络μk=f1(Xk)
�
�
=
�
1
(
�
�
)
，logσ2k=f2(Xk)
�
�
�
�
�
2
=
�
2
(
�
�
)
去计算「值得一提的是，选择拟合logσ2k
�
�
�
�
�
2
而不是直接拟合σ2k
�
�
2
，是因为σ2k
�
�
2
总是非负的，需要加激活函数处理，而拟合logσ2k
�
�
�
�
�
2
不需要加激活函数，因为它可正可负」
另一个根据生成的隐变量Z
�
的变分概率分布p(Z)
�
(
�
)
，还原生成原始数据的近似概率分布p(X^|Z)
�
(
�
^
|
�
)
，称为生成网络 
因为已经学到了这 k
�
 个正态分布，那可以直接从专属分布p(Z|Xk)
�
(
�
|
�
�
)
中采样一个Zk
�
�
出来，然后经过一个生成器得到Xk^=f(Zk)
�
�
^
=
�
(
�
�
)
，那接下来只需要最小化方差 D2(Xk,Xk^)
�
2
(
�
�
,
�
�
^
)
 就行

仔细理解的时候有没有发现一个问题？为什么在文章最开头，我们强调了没法直接比较 X
�
 与 X^
�
^
 的分布，而在这里，我们认为可以直接比较这俩？注意，这里的 Zk
�
�
 是专属于或针对于Xk
�
�
的隐变量，那么和 Xk^
�
�
^
本身就有对应关系，因此右边的蓝色方框内的“生成器”，是一一对应的生成。

另外，大家可以看到，均值和方差的计算本质上都是encoder。也就是说，VAE其实利用了两个encoder去分别学习均值和方差

1.2.2 VAE的Variational到底是个啥：需要方差持续存在从而带来噪声
这里还有一个非常重要的问题，如苏剑林所说，由于我们通过最小化D2(Xk,Xk^)
�
2
(
�
�
,
�
�
^
)
来训练右边的生成器，最终模型会逐渐使得 Xk
�
�
 和Xk^
�
�
^
趋于一致。但是注意，因为 Zk
�
�
是重新随机采样过的，而不是直接通过均值和方差encoder学出来的，这个生成器的输入 Z
�
是有噪声的

仔细思考一下，这个噪声的大小其实就用方差来度量。为了使得分布的学习尽量接近，我们希望噪声越小越好，所以我们会尽量使得方差趋于 0
但是方差不能为 0，因为我们还想要给模型一些训练难度。如果方差为 0，模型永远只需要学习高斯分布的均值，这样就丢失了随机性，VAE就变成AE了……这就是为什么VAE要在AE前面加一个Variational：我们希望方差能够持续存在，从而带来噪声
那如何解决这个问题呢？其实保证有方差就行，但是VAE给出了一个优雅的答案：不仅需要保证有方差，还要让所有 p(Z|X)
�
(
�
|
�
)
趋于标准正态分布N(0,1)
�
(
0
,
1
)
，根据定义可知

这个式子的关键意义在于告诉我吗：如果所有p(Z|X)
�
(
�
|
�
)
都趋于N(0,1)
�
(
0
,
1
)
，那么我们可以保证p(Z)
�
(
�
)
也趋于N(0,1)
�
(
0
,
1
)
，从而实现先验的假设，这样就形成了一个闭环！那怎么让所有p(Z|X)
�
(
�
|
�
)
趋于N(0,1)
�
(
0
,
1
)
呢？还是老套路：加loss

到此为止，我们可以把VAE进一步画成：

现在我们来回顾一下VAE到底做了啥。VAE在AE的基础上

一方面，对均值的encoder添加高斯噪声(正态分布的随机采样)，使得decoder(即生成器)有噪声鲁棒性
二方面，为了防止噪声消失，将所有p(Z|X)
�
(
�
|
�
)
趋近于标准正态分布，将encoder的均值尽量降为 0，而将方差尽量保持住
这样一来，当decoder训练的不好的时候，整个体系就可以降低噪声；当decoder逐渐拟合的时候，就会增加噪声

1.3 VAE的改进：VQ-VAE/VQ-VAE2
1.3.1 什么是VQ-VAE
VQ即Vector Quantised，它编码出的向量是离散的，也就是把VAE做量化，所以VQ-VAE最后得到的编码向量的每个元素都是一个整数

现实生活中，很多信息(声音、图片)都是连续的，你的大部分任务都是一个回归任务。但是等你真正将其表示出来或真正解决这些任务的时候，我们都将其离散化了。图像变成了像素，语音也抽样过了，大部分工作的很好的也都是分类模型(回归任务转换成分类任务)
如果还是之前VAE的模式，就不好把模型做大，分布也不好学
故最终，取而代之的不是去直接预测分布z
�
，而是用一个codebook代替。codebook可以理解为聚类的中心，大小一般是K*D（K=8192，Dim=512/768），也就是有8192个长为D的向量

1.3.2 VQ-VAE的算法流程
x
�
输入编码器得到高宽分别为(h,w)
(
ℎ
,
�
)
的特征图f
�
然后计算特征图里的向量和codebook里的向量(聚类中心)的相似性
接着把和特征图最接近的聚类中心向量的编号(比如1-8192)存到矩阵z
�
里面
训练完成之后，不再需要编码特征f
�
，而是取出矩阵z
�
中的编号对应的codebook里面的向量，生成一个新的特征图fq
�
�
 (经过量化后的特征，即quantised feature)
最后和之前一样，使用fq
�
�
解码重构原图
此时这个量化特征就非常可控了，因为它们永远都是从codebook里面来的，而非随机生成，这样优化起来相对容易，如下图所示

左图：VQ-VAE的模型结构
其中，红色线的梯度▽zL
▽
�
�
迫使encoder在下一次forword时改变其输出(参数更新)
由于编码器的输出和解码器的输入共享D维空间，梯度包含了编码器如何改变参数以降低损失的有效信息
右图：embedding space可视化。编码器输出z(x)
�
(
�
)
会mapped到最相近（nearest）的点e2
�
2
VQ-VAE也可以用来做CV领域的自监督学习，比如BEIT就是把DALL·E训练好的codebook拿来用。将图片经过上面同样的过程quantise成的特征图作为ground truth，自监督模型来训练一个网络

后续还有VL-BEIT (vision language BEIT)的工作，也是类似的思路，只不过是用一个Transformer编码器来做多模态的任务

1.3.2 VQ-VAE2
// 待更

第二部分 扩散模型DDPM：先前向加噪后反向去噪从而建立噪声估计模型
在写本部分之前，我反复看了网上很多阐述DDPM的文章，实话说，一开始看到那种一上来就一堆公式的，起初基本看不下去，虽然后来慢慢的都看得下去了，但如果对于一个初次接触DDPM的初学者来说，一上来就一堆公式确实容易把人绕晕，但如果没有公式，则又没法透彻理解背后的算法步骤

两相权衡，本文将侧重算法每一步的剖析，而公式更多为解释算法原理而服务，说白了，侧重原理 其次公式，毕竟原理透彻了，公式也就自然而然的能推出来了

言归正传，咱们先来了解下扩散模型的极简发展史

2.1 扩散模型发展史：DDPM、improved DDPM、Latent Diffusion Model到DALL·E/DALL·E2
2.1.1 从扩散模型概念的提出到DDPM
1  2015年，Sohl-Dickstein提出「扩散模型」的概念
2015年，斯坦福大学的一博士后Sohl-Dickstein通过此篇论文《Deep Unsupervised Learning using Nonequilibrium Thermodynamics》提出扩散模型的概念
简单来讲，扩散模型的灵感来自非平衡热力学，通过定义了一个扩散步骤的马尔可夫链，以缓慢地将「符合高斯分布的随机噪声」添加到数据中，然后反转扩散过程以从噪声中构建所需的数据样本

2  2019年，斯坦福的宋飏等人估计数据分布梯度的生成模型
随后，2019年，斯坦福一在读博士宋飏和其导师通过此文《Generative Modeling by Estimating Gradients of the Data Distribution》提出了一种新方法来构建生成模型：即不需要估计数据的概率分布(数据概率的分布类似高维曲面)，相反，它估计的是分布的梯度(分布的梯度可以看成是高维曲面的斜率)

「顺带说一下，后来宋飏等人推出了扩散模型的改进，即DALLE 3的解码器之一致性模型Consistency Models(详见：AI绘画神器DALLE 3的解码器：一步生成的扩散模型之Consistency Models)」

3  2020年，UC Berkeley的Jonathan Ho等人正式提出：DDPM
再之后，2020年6月，UC Berkeley的Jonathan Ho等人意识到宋飏的工作可以改进 Sohl-Dickstein的扩散模型，很快，便通过论文《Denoising Diffusion Probabilistic Models》正式提出对于普通扩散模型的改进版：DDPM(全称即论文名称：Denoising Diffusion Probabilistic Models)

DDPM主要有两个贡献

一方面，从预测转换图像改进为预测噪声 (即如DiT论文所说，reformulating diffusion models to predict noise instead of pixel，可惜强调这点的文章太少了，可它是DDPM的关键，更是DDPM的本质）
作者认为，每次直接从xt
�
�
预测xt−1
�
�
−
1
，这种图像到图像的转化不太好优化，所以直接去预测从xt
�
�
到xt−1
�
�
−
1
这一步所添加的噪声ε
�
，这样就简化了问题：毕竟噪声一旦被预测出来，xt
�
�
减去噪声即得xt−1
�
�
−
1

这种操作就有点类似ResNet的残差结构。每次新增一些层，模型不是直接从x
�
 去预测y
�
，而是让新增的层去预测(y−x
�
−
�
)。这样新增层不用全部重新学习，而是学习原来已经学习到的x
�
 和真实值y
�
 之间的残差就行(residual)

DDPM采用了一个U-Net 结构的Autoencoder来对t
�
 时刻的高斯噪声z
�
进行预测，训练目标即希望预测的噪声fθ(xt,t)
�
�
(
�
�
,
�
)
和真实的噪声z
�
一致，所以目标函数为fθ(xt,t)
�
�
(
�
�
,
�
)
和z
�
的L1
�
1
 Loss：
p(xt−1∣xt)=∥z−fθ(xt,t)∥
�
(
�
�
−
1
∣
�
�
)
=
‖
�
−
�
�
(
�
�
,
�
)
‖

→
→
  这里的标签z
�
是正向扩散过程中，我们每一步真实添加的噪声(所以噪声是已知的，预测噪声时，其可以拿来当作Ground truth)
→
→
  这里的fε
�
�
就对应了U-Net 模型结构
→
→
  至于t
�
 就是U-Net 的输入xt
�
�
之外的另一个输入time embedding「类似transformer里的正弦位置编码，主要用于告诉 U-Net模型，现在到了反向过程的第几步，相当于引导U-Net生成」
最终，通过这个简单的L1
�
1
损失函数，模型就可以训练起来了
二方面，DDPM只预测正态分布的均值
虽然正态分布由均值和方差决定，但作者在这里发现，其实模型不需要学方差，只需要学习均值就行。逆向过程中高斯分布的方差项直接使用一个常数，模型的效果就已经很好。所以就再一次降低了模型的优化难度
为方便大家更好的理解本文，特地解释下什么叫U-Net网络(对应论文为：U-Net: Convolutional Networks for Biomedical Image Segmentation)

在目前绝大部分的图像或视频扩散方法中，主导的骨干网络一般是由一系列卷积和自注意力层构成的 U-Net 架构

它总共有23个卷积层的网络，由一个收缩路径(左侧)和一个扩张路径(右侧)组成

收缩路径遵循卷积网络的典型架构。 它由两个3x3卷积(无填充卷积)的重复结构组成，每个卷积后面跟着ReLU和一个2x2最大池化操作，步长为2，用于下采样，在每个下采样步骤中，我们将特征通道的数量加倍
扩张路径中的每个步骤都包括特征图的上采样，然后是一个2x2卷积(“上卷积”)，将特征通道数量减半，与收缩路径中相应裁剪的特征图进行连接(a concatenation with the correspondingly cropped feature map from the contracting path)，以及两个3x3卷积，每个卷积后面跟着一个ReLU。由于每次卷积都会丢失边界像素(the loss of border pixels in every convolution)，所以裁剪是必要的
在最后一层，使用1x1卷积将每个64个分量的特征向量映射到所需的类别数
人们之所以偏好 U-Net，是因为 Transformer 中全注意力机制的内存需求会随输入序列长度而二次方增长，而在处理视频这样的高维信号时，这样的增长模式会让计算成本变得非常高

当然，在此文《视频生成的原理解析：从Gen2、Emu Video到PixelDance、SVD、Pika 1.0、W.A.L.T》的第六部分你会看到，已有最新的研究把Transformer用做扩散模型的骨干网络

DDPM也有些类似VAE，也可以将其当做一个encoder-decoder的结构，但是有几点区别：

扩散过程是编码器一步步的走到xt
�
�
，而且是一个固定的过程；而VAE的编码器是可以学习的；
DDPM的每一步输出输出都是同样维度大小的，但对一般的自编码器（AE/VAE等），往往中间的bottleneck特征会比输入小很多
扩散模型有步数step的概念（time step、time embedding），模型要经过很多步才能生成图片，且在所有step中，U-Net都是共享参数的
2.1.2 improved DDPM、Diffusion Model Beat GANs、Latent Diffusion Model到DALL·E、DALL·E2
4  2020年底，OpenAI提出improved DDPM
DDPM使得扩散模型可以在真实数据集上work得很好之后，一下子吸引了很多人的兴趣。因为DDPM在数学上简洁美观，无论正向还是逆向，都是高斯分布，可以做很多推理证明，而且还有很多不错的性质，2020年底左右，OpenAI推出了 improved DDPM

improved DDPM相比DDPM做了几点改动：

DDPM的逆向过程中，高斯分布的方差项直接使用一个常数而不用学习
improved DDPM作者就觉得如果对方差也进行学习的话，效果应该会更好，改了之后果然取样和生成效果都好了很多
DDPM添加噪声时采用的线性的variance schedule改为余弦schedule，效果更好(类似学习率从线性改为余弦)
简单尝试了scale大模型之后，生成效果更好
5  2021年5月 Diffusion Model Beat GANs：使用classifier guidance的方法，引导模型进行采样和生成
上面第三点对OpenAI来说，无疑是个好消息。所以improved DDPM的二作和三作马上着手研究，发布了《Diffusion Models Beat GANs on Image Synthesis》这篇论文，比之前的improved DDPM又做了一些改进：

使用大模型：加大加宽网络、使用更多的自注意力头attention head，加大自注意力scale (single-scale attention改为multi-scale attention)
提出了新的归一化方式——Adaptive Group Normalization，相当于根据步数进行自适应的归一化，这个方法是对group归一化的一个改进：AdaGN(h,y=[ys,yb])=ysGroupNorm(h)+yb
AdaGN
(
ℎ
,
�
=
[
�
�
,
�
�
]
)
=
�
�
GroupNorm
(
ℎ
)
+
�
�

上面公式中的h
ℎ
是残差块激活函数的输出，y
�
是一个线性层对时步和后面用到的类别信息的嵌入。组归一化是对输入的通道方向进行分组归一化的归一化方法，可以理解为局部LayerNorm

使用classifier guidance的方法，引导模型进行采样和生成
这样不仅使生成的图片更逼真，而且加速了反向采样过程。论文中，只需要25次采样，就可以从噪声生成图片

所谓classifier guided diffusion
1 即在反向过程训练U-Net的同时，也训练一个简单的图片分类器。这个分类器是在ImageNet上训练的，只不过图片加了很多噪声 (毕竟扩散模型的输入始终是加了很多噪声的，跟真实的ImageNet图片是很不一样的，是从头训练的)
2 当采样xt
�
�
之后，直接扔给分类器，就可以看到图片分类是否正确，这时候就可以算一个交叉熵目标函数，对应的就得到了一个梯度。之后使用分类器对xt
�
�
的梯度信息指导扩散模型的采样和生成
3 这个梯度暗含了当前图片是否包含物体，以及这个物体是否真实的信息。通过这种梯度的引导，就可以帮助U-Net将图片生成的更加真实，要包含各种细节纹理，而不是意思到了就行，要和真实物体匹配上
 
当然，除了最简单最原始的classifier guidance之外，还有很多其它的引导方式
→
→
  CLIP guidance：将简单的分类器换成CLIP之后，文本和图像就联系起来了。此时不光可以利用这个梯度引导模型采用和生成，而且可以利用文本指导其采样和生成
→
→
  image侧引导：除了利用图像重建进行像素级别的引导，还可以做图像特征和风格层面的引导，只需要一个gram matrix就行
→
→
  text 侧：可以用训练好的NLP大模型做引导
以上所有引导方式，都是下面目标函数里的y
�
，即模型的输入不光是xt
�
�
和time embedding，还有condition，加了condition之后，可以让模型的生成又快又好p(xt−1∣xt)=∥z−fθ(xt,t,y)∥
�
(
�
�
−
1
∣
�
�
)
=
‖
�
−
�
�
(
�
�
,
�
,
�
)
‖

且值得一提的是，额外引入一个网络来指导，推理的时候比较复杂 (扩散模型需要反复迭代，每次迭代都需要额外算一个分数)，所以引出了后续2022年7月的一个工作：classifier free guidance

6  2021年12月 潜在扩散空间Latent Diffusion Model
2021年年底，此篇论文《High-Resolution Image Synthesis with Latent Diffusion Models》提出了潜在扩散模型，也是后续奠定stable diffusion的核心论文(关于SD详见本系列的另一篇文章《AI绘画原理解析：从CLIP到DALLE1/2、DALLE 3、Stable Diffusion》)

为何要弄这么个隐空间或潜在空间呢？

原因很简单，为了使扩散模型在有限的计算资源上训练，并且保留它们的质量和灵活性，故首先训练了一个强大的预训练自编码器，这个自编码器所学习到的是一个潜在的空间，这个潜在的空间要比像素空间要小的多(可以简单粗暴的理解为就是一个被压缩或被降维的空间)，把扩散模型在这个潜在的空间去训练，大大的降低了对算力的要求，这也是Stable Diffusion比原装Diffusion速度快的原因

7  2022年7月 Classifier-Free Diffusion Guidance
所谓classifier free guidance的方式(对应论文为《Jonathan Ho and Tim Salimans. Classifier-free diffusion guidance》)，只是改变了模型输入的内容，除了 conditional输入外(随机高斯噪声输入加引导信息)，还有 unconditional 的采样输入，两种输入都会被送到同一个 diffusion model，从而让其能够具有无条件和有条件生成的能力

得到有条件输出fθ(xt,t,y)
�
�
(
�
�
,
�
,
�
)
和无条件输出fθ(xt,t,ϕ)
�
�
(
�
�
,
�
,
�
)
后，就可以用前者监督后者，来引导扩散模型进行训练了
最后反向扩散做生成时，我们用无条件的生成，也能达到类似有条件生成的效果，这样一来就摆脱了分类器的限制，所以叫classifier free guidance
比如在训练时使用图像-文本对，这时可以使用文本做指导信号，也就是训练时使用文本作为y
�
生成图像。然后y
�
把去掉，替换为一个空集ϕ
�
（空的序列），生成另外的输出
总之，扩散模型本来训练就很贵了，classifier free guidance这种方式在训练时需要生成两个输出，所以训练更贵了。但是这个方法确实效果好，所以在GLIDE 、DALL·E2和Imagen里都用了，而且都提到这是一个很重要的技巧，用了这么多技巧之后，GLIDE终于是一个很好的文生图模型了，只用了35亿参数，生成效果和分数比120亿参数的DALL·E还要好

8  2021-2022年 DALL·E/DALL·E2：条件引导生成
2021年，OpenAI一看GLIDE这个方向靠谱，就马上跟进，不再考虑DALL·E的VQ-VAE路线了，而是将GLIDE改为层级式生成（56→256→1024）并加入prior网络等等，于是

Jonathan Ho和他在谷歌研究中心的同事Tim Salimans，与其他地方的团队合作，展示了如何结合大型语言模型的信息与图像生成扩散模型，即用文本（比如“金鱼在海滩上喝可口可乐”）指导扩散过程，从而生成图像，最终得到了DALL·E2 (其具体的训练细节在下一篇文章《AI绘画与多模态原理解析：从CLIP到DALLE 3、Stable Diffusion、MDJ》会重点讲解 ）
总之，DALL·E2 这样的从文本到图像模型成功的背后原因，就是这种“引导扩散”(guided diffusion)过程

回到DDPM，每一个噪声都是在前一时刻增加噪声而来的，从最开始的x0
�
0
时刻开始，最终得到xT
�
�
时刻的纯噪声图像。不过问题来是为什么要加噪声？

Diffusion的最终目标是去噪以生成图片，而为了推导出逆向的去噪方法，必须了解增加噪声的原理。同时，添加噪声的过程其实就是不断构建标签的过程。如果在前一时刻可以预测出来后一时刻的噪声，便能很方便地实现还原操作 (就和人走路一样，不管你从哪来，哪怕走过万水千山，最后都可按原路返回至原出发点)
说白了 当你学会了怎么加噪(前向扩散)，就一定能知道怎么去噪(逆向生成)，毕竟知道怎么来 也必知道怎么回

且在噪声的添加过程中，每一步都要保持尽量相同的噪声扩散幅度。比如，在给上图加噪的过程中，前期的分布非常均匀，添加一些噪声便可以将原始分布改变，但到后期，需要添加更多的噪声，方可保证噪声扩散幅度相同(这就像往水中加糖，为了使糖的甜味增长相同，后期需要加更多的糖)
所以DDPM为了从随机噪声中直接生成图片，首先需要训练一个噪声估计模型，然后将输入的随机噪声还原成图片，相当于就两个关键，一个是训练过程，一个是推理过程

训练过程：随机生成噪声ϵ
�
，经过T
�
步将噪声扩散到输入原始图片x0
�
0
中，破坏后的图片xT
�
�
，学习破坏图片的预估噪声ϵθ(xt,t)
�
�
(
�
�
,
�
)
，并用L2 loss约束预估噪声ϵθ(xt,t)
�
�
(
�
�
,
�
)
与原始输入噪声的距离
推理过程：即输入噪声，经过预估噪声模型还原成图片
2.2 DDPM的两个过程：从前向过程到逆向过程
2.2.1 前向过程(加噪)：通过高斯噪音随机加噪ϵ
�
 ——给图片打马赛克
前向过程(forward process)也称为扩散过程(diffusion process)，简单理解就是对原始图片x0
�
0
通过逐步添加「方差为βt
�
�
的高斯噪声」变成xT
�
�
，从而达到破坏图片的目的，如下图 

在从xt−1
�
�
−
1
到xt
�
�
的过程中，其对应的分布q(xt|xt−1)
�
(
�
�
|
�
�
−
1
)
是一个正太分布，且其均值是ut=1−βt−−−−−√xt−1
�
�
=
1
−
�
�
�
�
−
1
，方差为βt
�
�
，则有 

q(xt|xt−1)=N(xt;ut=1−βt−−−−−√xt−1,βtI)
�
(
�
�
|
�
�
−
1
)
=
�
(
�
�
;
�
�
=
1
−
�
�
�
�
−
1
,
�
�
�
)

对于这个公式，解释下3点

正态分布的概率密度函数具有以下形式：

f(x)=1σ2π√e−(x−μ)22σ2
�
(
�
)
=
1
�
2
�
�
−
(
�
−
�
)
2
2
�
2

顺带帮你再回顾下：正态分布有两个参数：均值μ
�
和方差σ2
�
2
。其中，μ
�
是分布的均值，决定了分布的中心位置，σ
�
是标准差，决定了分布的宽度

上面的方差之所以表示为βtI
�
�
�
，原因在于我们一般处于多维情况下，而I
�
是单位矩阵，表明每个维度有相同的标准偏差βt
�
�

且βt∈(0,1)
�
�
∈
(
0
,
1
)
是事先给定的常量，代表从xt−1
�
�
−
1
到xt
�
�
这一步的方差，且正因为βt
�
�
设置的比较小，所以使得q(xt|xt−1)
�
(
�
�
|
�
�
−
1
)
的均值在xt−1
�
�
−
1
附近
换言之，xt
�
�
相当于就是在xt−1
�
�
−
1
的基础上加了一些噪声，而且是渐进式逐步增加/扩散的，当然 从加噪大小的角度上讲，前期加噪较弱，后期加噪加强，所以在DDPM的论文中，作者取β1,⋯βT
�
1
,
⋯
�
�
为从0.0001到0.02的线性递增序列
此外，值得一提的是，因为是马尔可夫链，所以其联合分布便是：
q(x1:T∣x0)=q(x1,x2,…,xT∣x0)=∏Tt=1q(xt∣xt−1)
�
(
�
1
:
�
∣
�
0
)
=
�
(
�
1
,
�
2
,
…
,
�
�
∣
�
0
)
=
∏
�
=
1
�
�
(
�
�
∣
�
�
−
1
)

如下图所示

接下来，如果我们定义 αt=1−βt
�
�
=
1
−
�
�
,   且{αt}Tt=1
{
�
�
}
�
=
1
�
『被称为Noise schedule，通常是一些列很小的值』，以及 ϵt−1∼N(0,1)
�
�
−
1
∼
�
(
0
,
1
)
是高斯噪声，便可以得到xt
�
�
的采样值

xt=αt−−√xt−1+1−αt−−−−−√ϵt−1
�
�
=
�
�
�
�
−
1
+
1
−
�
�
�
�
−
1

把上述公式迭代变换下，可以直接得出 x0
�
0
 到 xt
�
�
 的公式，如下：

xt=α¯t−−√x0+1−α¯t−−−−−√ϵ0
�
�
=
�
¯
�
�
0
+
1
−
�
¯
�
�
0

其中 αt¯=∏ti=1αi
�
�
¯
=
∏
�
=
1
�
�
�
，ϵ∼N(0,1)
�
∼
�
(
0
,
1
)
 也是一个高斯噪声

换言之，所以 xt
�
�
 在 x0
�
0
条件下的分布就是均值为 αt¯−−√x0
�
�
¯
�
0
 , 方差为 1−αt¯
1
−
�
�
¯
的正态分布 (下式的意义在于，只需要给出x0
�
0
，便可以计算出任意时刻t
�
的xt
�
�
)

q(xt∣x0)=N(xt;α¯t−−√x0,(1−α¯t)I)
�
(
�
�
∣
�
0
)
=
�
(
�
�
;
�
¯
�
�
0
,
(
1
−
�
¯
�
)
�
)

考虑到可能会有读者对这个x0
�
0
 到 xt
�
�
 的一步到位感到困惑，而一般的同类文章不会展开的特别细，故本文细致展开下(能拆10步则10步，确保阅读无障碍)

首先通过xt=αt−−√xt−1+1−αt−−−−−√ϵt−1
�
�
=
�
�
�
�
−
1
+
1
−
�
�
�
�
−
1
可知，xt−1=αt−1−−−−√xt−2+1−αt−1−−−−−−−√ϵt−2
�
�
−
1
=
�
�
−
1
�
�
−
2
+
1
−
�
�
−
1
�
�
−
2
，把这个代入到xt
�
�
的表达式后，再展开即可得
xt=αt−−√xt−1+1−αt−−−−−√ϵ∗t−1=αt−−√(αt−1−−−−√xt−2+1−αt−1−−−−−−−√ϵ∗t−2)+1−αt−−−−−√ϵ∗t−1=αtαt−1−−−−−√xt−2+αt−αtαt−1−−−−−−−−−√ϵ∗t−2+1−αt−−−−−√ϵ∗t−1
�
�
=
�
�
�
�
−
1
+
1
−
�
�
�
�
−
1
∗
=
�
�
(
�
�
−
1
�
�
−
2
+
1
−
�
�
−
1
�
�
−
2
∗
)
+
1
−
�
�
�
�
−
1
∗
=
�
�
�
�
−
1
�
�
−
2
+
�
�
−
�
�
�
�
−
1
�
�
−
2
∗
+
1
−
�
�
�
�
−
1
∗
考虑到「两个独立正态分布的随机变量之和是正态的，其均值是两个均值之和，其方差是两个方差之和（即标准差的平方是标准差的平方），比如两个方差不同的高斯分布N(0,σ21I)
�
(
0
,
�
1
2
�
)
和N(0,σ22I)
�
(
0
,
�
2
2
�
)
相加等于一个新的高斯分布N(0,(σ21+σ22)I)
�
(
0
,
(
�
1
2
+
�
2
2
)
�
)
」，然后再通过重参数技巧可得
xt=αtαt−1−−−−−√xt−2+​αt−αtαt−1−−−−−−−−−√2+1−αt−−−−−√2−−−−−−−−−−−−−−−−−−−−−√ϵt−2=αtαt−1−−−−−√xt−2+αt−αtαt−1+1−αt−−−−−−−−−−−−−−−−√ϵt−2=αtαt−1−−−−−√xt−2+1−αtαt−1−−−−−−−−−√ϵt−2=…
�
�
=
�
�
�
�
−
1
�
�
−
2
+
​
�
�
−
�
�
�
�
−
1
2
+
1
−
�
�
2
�
�
−
2
=
�
�
�
�
−
1
�
�
−
2
+
�
�
−
�
�
�
�
−
1
+
1
−
�
�
�
�
−
2
=
�
�
�
�
−
1
�
�
−
2
+
1
−
�
�
�
�
−
1
�
�
−
2
=
…

对此，本文参考文献中的这篇《Understanding Diffusion Models: A Unified Perspective》也解释了这几个步骤
最后定义一个累积混合系数， αt¯=∏ti=1αi
�
�
¯
=
∏
�
=
1
�
�
�
，即αt¯−−√=αt−−√αt−1−−−−√⋯α1−−√
�
�
¯
=
�
�
�
�
−
1
⋯
�
1
，可得
xt=∏ti=1αi−−−−−−√x0+1−∏ti=1αiϵ0−−−−−−−−−−−√=α¯t−−√x0+1−α¯t−−−−−√ϵ0∼N(xt;α¯t−−√x0,(1−α¯t)I)
�
�
=
∏
�
=
1
�
�
�
�
0
+
1
−
∏
�
=
1
�
�
�
�
0
=
�
¯
�
�
0
+
1
−
�
¯
�
�
0
∼
�
(
�
�
;
�
¯
�
�
0
,
(
1
−
�
¯
�
)
�
)
2.2.2 逆向过程(去噪)：求解真实后验分布q(xt−1|xt)
�
(
�
�
−
1
|
�
�
)
—— 复原被加噪的图片使之清晰化
2.2.2.1 没法直接通过x0=1α¯t√(xt−1−α¯t−−−−−√ϵ)
�
0
=
1
�
¯
�
(
�
�
−
1
−
�
¯
�
�
)
去做逆向
逆向过程就是通过估测噪声，多次迭代逐渐将被破坏的 xt
�
�
恢复成 x0
�
0
 ，如下图

更具体而言，正向扩散和逆扩散过程都是马尔可夫，唯一的区别就是正向扩散里每一个条件概率的高斯分布的均值和方差都是已经确定的（依赖于 βt
�
�
 和 x0
�
0
），而逆扩散过程里面的均值和方差需要通过网络学出来，怎么个学法呢？

有人可能要说，直接把上一节得到的xt=αt¯−−√x0+1−αt¯−−−−−√ϵ0
�
�
=
�
�
¯
�
0
+
1
−
�
�
¯
�
0
移个项不就行了(先把带x0
�
0
的项移到等式左边，然后所有项各自除以α¯−−√
�
¯
，最后把等式右边的1/α¯−−√
1
/
�
¯
提取到括号外边即可)？
x0=1α¯t√(xt−1−α¯t−−−−−√ϵ0)
�
0
=
1
�
¯
�
(
�
�
−
1
−
�
¯
�
�
0
)

但问题在于xt=αt¯−−√x0+1−αt¯−−−−−√ϵ0
�
�
=
�
�
¯
�
0
+
1
−
�
�
¯
�
0
中的ϵ
�
是个随机变量，意味着 xt
�
�
 也是个随机变量，其具体取值由 ϵ
�
实际取值决定「相当于现在我们有一个具体的 xt
�
�
，它对应着 ϵ
�
 的某个取值，但是什么值我们并不知道」，所以我们只能以前向过程的 ϵ
�
 取值为标签，训练一个模型去估计它，即：
xθ(xt,t):=1α¯t√(xt−1−α¯t−−−−−√ϵθ(xt,t))
�
�
(
�
�
,
�
)
:=
1
�
¯
�
(
�
�
−
1
−
�
¯
�
�
�
(
�
�
,
�
)
)

其中
→
→
  ϵθ(xt,t)
�
�
(
�
�
,
�
)
就是所谓的模型，用来近似真实的(即前向过程采样出来的)ϵ
�

→
→
  相应地，xθ(xt,t)
�
�
(
�
�
,
�
)
就是 x0
�
0
 的近似，或者，你也可以无视ϵ
�
，直接把 xθ(xt,t)
�
�
(
�
�
,
�
)
 视为模型

为了训练它，最直接的想法就是用 L2 损失∥ϵ−ϵθ(xt,t)∥2
‖
�
−
�
�
(
�
�
,
�
)
‖
2
或者∥x0−xθ(xt,t)∥2
‖
�
0
−
�
�
(
�
�
,
�
)
‖
2

理论上没问题，但是实际效果很差，为什么呢？如果直接用 xθ(xt,t)
�
�
(
�
�
,
�
)
，那么中间的x2,x3,…,xT−1
�
2
,
�
3
,
…
,
�
�
−
1
都没用了，整个 DDPM 就退化成了 VAE 的结构，但是
→
→
  VAE 的生成模型和后验都是自己学习出来的，二者双向奔赴共同优化去寻找最优解
→
→
  而 DDPM 的后验是人为指定的（即x0=1/α¯t−−√(xt−1−α¯t−−−−−√ϵ0)
�
0
=
1
/
�
¯
�
(
�
�
−
1
−
�
¯
�
�
0
)
），并且由于 α¯t→0
�
¯
�
→
0
，q(xt|x0)
�
(
�
�
|
�
0
)
基本上就是一个标准正态分布，磨灭掉了几乎所有的输入信息，全靠生成模型这一边去恢复，难度未免过大..

所以，实际应用中，我们是一点一点来的，比如先生成xt−1
�
�
−
1
、然后xt−2
�
�
−
2
……由于每一步的变化都比较小，保留了上一步足够的信息，生成模型的负担就轻了很多
如果我们能够逆转前向过程并从真实分布 q(xt−1|xt)
�
(
�
�
−
1
|
�
�
)
采样，就可以从高斯噪声 xt∼N(0,I)
�
�
∼
�
(
0
,
�
)
还原出原图分布 x0∼q(x)
�
0
∼
�
(
�
)
 

因为我们可以证明如果前向q(xt|xt−1)
�
(
�
�
|
�
�
−
1
)
满足高斯分布且 βt
�
�
 足够小，其逆向q(xt−1|xt)
�
(
�
�
−
1
|
�
�
)
仍然是一个高斯分布
那样，我们便可以使用「参数为 θ 的U-Net+attention 结构pθ
�
�
」去预测这样的一个逆向的分布(类似VAE)：
pθ(X0:T)=p(xt)∏Tt=1pθ(xt−1∣xt)
�
�
(
�
0
:
�
)
=
�
(
�
�
)
∏
�
=
1
�
�
�
(
�
�
−
1
∣
�
�
)

pθ(xt−1∣xt)=N(xt−1;μθ(xt,t),Σθ(xt,t))
�
�
(
�
�
−
1
∣
�
�
)
=
�
(
�
�
−
1
;
�
�
(
�
�
,
�
)
,
Σ
�
(
�
�
,
�
)
)

不过在DDPM的论文中，作者把条件概率 pθ(xt−1|xt)
�
�
(
�
�
−
1
|
�
�
)
 的方差直接取了βt
�
�
，而不是上面说的需要网络去估计的 Σθ(xt,t)
Σ
�
(
�
�
,
t
)
，所以说实际上只有均值需要网络去估计

然现在的问题是，我们无法直接去推断 q(xt−1|xt)
�
(
�
�
−
1
|
�
�
)
，即q(xt−1|xt)
�
(
�
�
−
1
|
�
�
)
 is unknown

所以，接下来的问题 自然而然 就转换成了我们希望求解q(xt−1∣xt)
�
(
�
�
−
1
∣
�
�
)
，因为我们已知前向过程q(xt∣xt−1)
�
(
�
�
∣
�
�
−
1
)
，所以自然想到使用贝叶斯公式：

q(xt−1∣xt)=q(xt∣xt−1)q(xt−1)q(xt)
�
(
�
�
−
1
∣
�
�
)
=
�
(
�
�
∣
�
�
−
1
)
�
(
�
�
−
1
)
�
(
�
�
)

可惜 q(xt)
�
(
�
�
)
和 q(xt−1)
�
(
�
�
−
1
)
是未知的，事情到这里似乎走入了僵局，好在我们发现q(xt|x0)
�
(
�
�
|
�
0
)
和 q(xt−1|x0)
�
(
�
�
−
1
|
�
0
)
是已知的，这样一变换，下述等式右边的三项就都可知了
q(xt−1∣xt)=q(xt∣xt−1)q(xt−1|x0)q(xt|x0)
�
(
�
�
−
1
∣
�
�
)
=
�
(
�
�
∣
�
�
−
1
)
�
(
�
�
−
1
|
�
0
)
�
(
�
�
|
�
0
)

相当于如果给上式加上 x0
�
0
 为条件，则立马柳暗花明，而一旦知道了 x0,q(xt−1|xt,x0)
�
0
,
�
(
�
�
−
1
|
�
�
,
�
0
)
，便可以直接写出q(xt−1∣xt,x0)=N(xt−1;μ~(xt,x0),β~tI)
�
(
�
�
−
1
∣
�
�
,
�
0
)
=
�
(
�
�
−
1
;
�
~
(
�
�
,
�
0
)
,
�
~
�
�
)

2.2.2.2 因q(xt−1|xt)
�
(
�
�
−
1
|
�
�
)
无法直接求解，故加上x0
�
0
：问题转换成了求解q(xt−1∣xt,x0)
�
(
�
�
−
1
∣
�
�
,
�
0
)
接下来，我们便好好推导下

解释下上面7.1~7.5这5个步骤的推导

7.1依据的是
P(A|B)=P(AB)P(B)
�
(
�
|
�
)
=
�
(
�
�
)
�
(
�
)

7.2中，分母部分依据的是
P(AB)=P(A)P(B|A)
�
(
�
�
)
=
�
(
�
)
�
(
�
|
�
)

分子部分依据的是
P(ABC)=P(A)P(B|A)P(C|AB)
�
(
�
�
�
)
=
�
(
�
)
�
(
�
|
�
)
�
(
�
|
�
�
)

注，此处的A B与上面7.1的A B非同一个具体的指向，只是公式层面的原有表达
7.3依据的是分子分母同时除以q(x0)
�
(
�
0
)
至于7.3到7.4
=q(xt∣xt−1,x0)q(xt−1∣x0)q(xt∣x0)7.3∝exp(−12((xt−αt√xt−1)2βt+(xt−1−α¯t−1√x0)21−a¯t−1−(xt−α¯t√x0)21−a¯t))7.4
=
�
(
�
�
∣
�
�
−
1
,
�
0
)
�
(
�
�
−
1
∣
�
0
)
�
(
�
�
∣
�
0
)
7.3
∝
exp
⁡
(
−
1
2
(
(
�
�
−
�
�
�
�
−
1
)
2
�
�
+
(
�
�
−
1
−
�
¯
�
−
1
�
0
)
2
1
−
�
¯
�
−
1
−
(
�
�
−
�
¯
�
�
0
)
2
1
−
�
¯
�
)
)
7.4

依据的是
ρ(x)=12πσ√e−12(x−μσ)2
�
(
x
)
=
1
2
�
�
e
−
1
2
(
x
−
�
�
)
2

且由前向扩散过程的特性『别忘了2.2.1节中，有 q(xt∣x0)=N(xt;α¯t−−√x0,(1−α¯t)I)
�
(
�
�
∣
�
0
)
=
�
(
�
�
;
�
¯
�
�
0
,
(
1
−
�
¯
�
)
�
)
』，可知
q(xt∣xt−1,x0)=q(xt∣xt−1)=N(xt;1−βt−−−−−√xt−1,βtI)
�
(
�
�
∣
�
�
−
1
,
�
0
)
=
�
(
�
�
∣
�
�
−
1
)
=
�
(
�
�
;
1
−
�
�
�
�
−
1
,
�
�
�
)

q(xt−1∣x0)=N(xt−1;α¯t−1−−−−√x0,(1−α¯t−1)I)
�
(
�
�
−
1
∣
�
0
)
=
�
(
�
�
−
1
;
�
¯
�
−
1
�
0
,
(
1
−
�
¯
�
−
1
)
�
)

q(xt∣x0)=N(xt;α¯t−−√x0,(1−α¯t)I)
�
(
�
�
∣
�
0
)
=
�
(
�
�
;
�
¯
�
�
0
,
(
1
−
�
¯
�
)
�
)
最后，再解释下怎么从7.4到的7.5
=q(xt∣xt−1,x0)q(xt−1∣x0)q(xt∣x0)7.3∝exp(−12((xt−αt√xt−1)2βt+(xt−1−α¯t−1√x0)21−a¯t−1−(xt−α¯t√x0)21−a¯t))=exp⎛⎝⎜⎜⎜⎜⎜−12((αtβt+11−α¯t−1)x2t−1xt−1−(2αt√βtxt+2a¯t−1√1−α¯t−1x0)xt−1xt−1+C(xt,x0))⎞⎠⎟⎟⎟⎟⎟⋅7.5
=
�
(
�
�
∣
�
�
−
1
,
�
0
)
�
(
�
�
−
1
∣
�
0
)
�
(
�
�
∣
�
0
)
7.3
∝
exp
⁡
(
−
1
2
(
(
�
�
−
�
�
�
�
−
1
)
2
�
�
+
(
�
�
−
1
−
�
¯
�
−
1
�
0
)
2
1
−
�
¯
�
−
1
−
(
�
�
−
�
¯
�
�
0
)
2
1
−
�
¯
�
)
)
=
exp
⁡
(
−
1
2
(
(
�
�
�
�
+
1
1
−
�
¯
�
−
1
)
�
�
−
1
2
⏟
�
�
−
1
−
(
2
�
�
�
�
�
�
+
2
�
¯
�
−
1
1
−
�
¯
�
−
1
�
0
)
�
�
−
1
⏟
�
�
−
1
+
�
(
�
�
,
�
0
)
⏟
)
)
⋅
7.5

先举一个最简单的例子，比如对于 12(Ax2+Bx+C)
1
2
(
�
�
2
+
�
�
+
�
)
，稍加转化下即是12A(x+B2A)2+C
1
2
�
(
�
+
�
2
�
)
2
+
�
，而这个
→
→
  A
�
则对应于7.5中的αtβt+11−α¯t−1
�
�
�
�
+
1
1
−
�
¯
�
−
1

→
→
  B
�
则对应于7.5中的−(2αt√βtxt+2a¯t−1√1−α¯t−1x0)
−
(
2
�
�
�
�
�
�
+
2
�
¯
�
−
1
1
−
�
¯
�
−
1
�
0
)

且其均值为−B2A
−
�
2
�
，方差为1A
1
�
，从而有
μt(xt,x0)=−B2A
�
�
(
�
�
,
�
0
)
=
−
�
2
�

=(αt√βtxt+α¯t−1√1−α¯t−1x0)/(αtβt+11−α¯t−1)=(αt√βtxt+α¯t−1√1−α¯t−1x0)1−α¯t−11−α¯t⋅βt=αt√(1−α¯t−1)1−α¯txt+α¯t−1√βt1−α¯tx0
=
(
�
�
�
�
�
�
+
�
¯
�
−
1
1
−
�
¯
�
−
1
�
0
)
/
(
�
�
�
�
+
1
1
−
�
¯
�
−
1
)
=
(
�
�
�
�
�
�
+
�
¯
�
−
1
1
−
�
¯
�
−
1
�
0
)
1
−
�
¯
�
−
1
1
−
�
¯
�
⋅
�
�
=
�
�
(
1
−
�
¯
�
−
1
)
1
−
�
¯
�
�
�
+
�
¯
�
−
1
�
�
1
−
�
¯
�
�
0

βt~=1A=1/(αtβt+11−α¯t−1)=1/(αt−α¯t+βtβt(1−α¯t−1))=1−α¯t−11−α¯t⋅βt
�
�
~
=
1
�
=
1
/
(
�
�
�
�
+
1
1
−
�
¯
�
−
1
)
=
1
/
(
�
�
−
�
¯
�
+
�
�
�
�
(
1
−
�
¯
�
−
1
)
)
=
1
−
�
¯
�
−
1
1
−
�
¯
�
⋅
�
�
2.2.2.3 q(xt−1|xt,x0)
�
(
�
�
−
1
|
�
�
,
�
0
)
中方差和均值的推导
好，接下来关键来了

根据xt=α¯t−−√x0+1−α¯t−−−−−√ϵ0
�
�
=
�
¯
�
�
0
+
1
−
�
¯
�
�
0
，可知x0=1α¯t√(xt−1−α¯t−−−−−√ϵ0)
�
0
=
1
�
¯
�
(
�
�
−
1
−
�
¯
�
�
0
)
，代入上面μt(xt,x0)
�
�
(
�
�
,
�
0
)
的表达式 可得

μ(xt,x0)=1αt√(xt−1−αt1−α¯t√ϵ0)
�
(
�
�
,
�
0
)
=
1
�
�
(
�
�
−
1
−
�
�
1
−
�
¯
�
�
0
)

大部分文章对上面这个的推导都是一步到位的，但本文为细致起见，故还是一步步来推导下

首先直接把x0=1α¯t√(xt−1−α¯t−−−−−√ϵ0)
�
0
=
1
�
¯
�
(
�
�
−
1
−
�
¯
�
�
0
)
和βt=1−αt
�
�
=
1
−
�
�
代入进去，可得
μ(xt,x0)=αt−−√(1−α¯t−1)xt+α¯t−1−−−−√(1−αt)x01−α¯t=αt−−√(1−α¯t−1)xt+α¯t−1−−−−√(1−αt)xt−1−α¯t√ϵ0α¯t√1−α¯t
�
(
�
�
,
�
0
)
=
�
�
(
1
−
�
¯
�
−
1
)
�
�
+
�
¯
�
−
1
(
1
−
�
�
)
�
0
1
−
�
¯
�
=
�
�
(
1
−
�
¯
�
−
1
)
�
�
+
�
¯
�
−
1
(
1
−
�
�
)
�
�
−
1
−
�
¯
�
�
0
�
¯
�
1
−
�
¯
�

接下来，我们可以进一步观察到 分子中的后半部分有α¯t−1−−−−√(1−αt)xt−1−α¯t√ϵ0α¯t√
�
¯
�
−
1
(
1
−
�
�
)
�
�
−
1
−
�
¯
�
�
0
�
¯
�
这一项，怎么进一步化简呢？
接下来非常关键(截止23年5月份之前，暂时没看到有其他中英文资料解释了这个细节)

好在之前有定义：αt¯=∏ti=1αi
�
�
¯
=
∏
�
=
1
�
�
�
，即αt¯−−√=αt−−√αt−1−−−−√⋯α1−−√
�
�
¯
=
�
�
�
�
−
1
⋯
�
1
，从而有αt¯−−√=αt−−√αt−1¯−−−−√
�
�
¯
=
�
�
�
�
−
1
¯
所以我们可以针对这一项α¯t−1−−−−√(1−αt)xt−1−α¯t√ϵ0α¯t√
�
¯
�
−
1
(
1
−
�
�
)
�
�
−
1
−
�
¯
�
�
0
�
¯
�
的分子分母同时除以αt−1¯−−−−√
�
�
−
1
¯
，得到

μt(xt,x0)=αt√(1−α¯t−1)xt+(1−αt)xt−1−α¯t√ϵ0αt√1−α¯t
�
�
(
�
�
,
�
0
)
=
�
�
(
1
−
�
¯
�
−
1
)
�
�
+
(
1
−
�
�
)
�
�
−
1
−
�
¯
�
�
0
�
�
1
−
�
¯
�
之后的推导就比较简单了

以下分别对上面的三行公式做解释说明：
→
→
  接着把上阶段2得到的式子的分子拆成三项，且三项中最后两项的分子分母同时乘以αt−−√
�
�

→
→
  然后再把上一步骤中分子三项中的前两项通过提取出xt
�
�
从而实现合并
→
→
  前两项合并之后，再对前两项中第一项的分子分母同时乘以αt−−√
�
�
，然后对第三项的分子分母同时除以1−αt−−−−−√
1
−
�
�
，即可得1−αt1−α¯t√αt√ϵ0
1
−
�
�
1
−
�
¯
�
�
�
�
0
，原因很简单，因为：1−αt¯=(1−αt¯−−−−−√)2
1
−
�
�
¯
=
(
1
−
�
�
¯
)
2

 接下来，针对上面阶段3得到的式子的前两项再做合并，合并中用到了一个细节，即αt×αt−1¯=αt¯
�
�
×
�
�
−
1
¯
=
�
�
¯
，原因也同样很简单，根据上面阶段2出现的这个式子αt¯−−√=αt−−√αt−1¯−−−−√
�
�
¯
=
�
�
�
�
−
1
¯
而来，再之后就更eazy 便不再赘述了

总之，从最终得到的结果可以看出，在给定 x0
�
0
 的条件下，后验条件高斯分布的均值只和超参数αt
�
�
、xt
�
�
、ϵ0
�
0
有关，即

μ(xt,x0)=1αt√(xt−1−αt1−α¯t√ϵ0)
�
(
�
�
,
�
0
)
=
1
�
�
(
�
�
−
1
−
�
�
1
−
�
¯
�
�
0
)

方差只与超参数α
�
有关，即

βt~=1A=1/(αtβt+11−α¯t−1)=1/(αt−α¯t+βtβt(1−α¯t−1))=1−α¯t−11−α¯t⋅βt
�
�
~
=
1
�
=
1
/
(
�
�
�
�
+
1
1
−
�
¯
�
−
1
)
=
1
/
(
�
�
−
�
¯
�
+
�
�
�
�
(
1
−
�
¯
�
−
1
)
)
=
1
−
�
¯
�
−
1
1
−
�
¯
�
⋅
�
�

从而通过以上的方差和均值，我们就得到了q(xt−1|xt,x0)
�
(
�
�
−
1
|
�
�
,
�
0
)
的解析形式

2.3 DDPM如何训练：最小化「噪声估计模型ϵθ(xt,t)
�
�
(
�
�
,
�
)
估计的噪声」与「真实噪声」之间的差距
继续下文之前，先总结一下

生成模型的本质是根据给定的样本(训练数据)生成新样本

具体而言，给定一批训练数据X
�
，假设其服从某种复杂的真实分布p(x)
�
(
�
)
，则给定的训练数据可视为从该真实分布中采样的观测样本x
�
如果能从这些观测样本x
�
中估计出训练数据的真实分布，相当于就可以从该分布(估计出的接近真实分布的分布)中不断的采样出新的样本了，故说白了，生产模型的目标就是估计训练数据的真实分布，并假定其真实分布为q(x)
�
(
�
)
从而问题自然而然就变成了尽可能缩小估计的分布q(x)
�
(
�
)
与真实分布p(x)
�
(
�
)
之间的差距
接下来介绍这个模型要怎么优化，即网络该怎么训练：去估计分布q(xt−1∣xt)
�
(
�
�
−
1
∣
�
�
)
的条件概率pθ(xt−1|xt)
�
�
(
�
�
−
1
|
�
�
)
的均值uθ(xt,t)
�
�
(
�
�
,
�
)
 和方差Σθ(xt,t)
Σ
�
(
�
�
,
t
)

与之前介绍的VAE相比，扩散模型的隐变量是和原始数据是同维度的，而且encoder（即扩散过程）是固定的

2.3.1 确立目标函数
既然扩散模型是隐变量模型，那么我们可以基于变分推断来得到variational lower bound（VLB，又称ELBO）作为最大化优化目标，当然实际训练时一般对VLB取负，即我们要最小化目标分布的负对数似然：

−logpθ(x0)≤−logpθ(x0)+DKL(q(x1:T∣x0)∥pθ(x1:T∣x0))=−logpθ(x0)+Ex1:T∼q(x1:T∣x0)[logq(x1:T∣x0)pθ(x0:T)/pθ(x0)]=−logpθ(x0)+Eq[logq(x1:T∣x0)pθ(x0:T)+logpθ(x0)]=Eq[logq(x1:T∣x0)pθ(x0:T)]
−
log
⁡
�
�
(
�
0
)
≤
−
log
⁡
�
�
(
�
0
)
+
�
K
L
(
�
(
�
1
:
�
∣
�
0
)
‖
�
�
(
�
1
:
�
∣
�
0
)
)
=
−
log
⁡
�
�
(
�
0
)
+
�
�
1
:
�
∼
�
(
�
1
:
�
∣
�
0
)
[
log
⁡
�
(
�
1
:
�
∣
�
0
)
�
�
(
�
0
:
�
)
/
�
�
(
�
0
)
]
=
−
log
⁡
�
�
(
�
0
)
+
�
�
[
log
⁡
�
(
�
1
:
�
∣
�
0
)
�
�
(
�
0
:
�
)
+
log
⁡
�
�
(
�
0
)
]
=
�
�
[
log
⁡
�
(
�
1
:
�
∣
�
0
)
�
�
(
�
0
:
�
)
]

考虑到本文的定位起见，逐一解释下上面推导的每一行

第一行：由 KL 散度的非负性质（KL 散度始终大于等于零），我们得到如下不等式：−logpθ(x0)≤−logpθ(x0)+DKL(q(x1:T∣x0)∥pθ(x1:T∣x0))
−
log
⁡
�
�
(
�
0
)
≤
−
log
⁡
�
�
(
�
0
)
+
�
K
L
(
�
(
�
1
:
�
∣
�
0
)
‖
�
�
(
�
1
:
�
∣
�
0
)
)

第二行：将 KL 散度的定义代入上式可得

其中 E
�
 表示期望，即对分布 q(x1:T|x0)
�
(
�
1
:
�
|
�
0
)
 中的所有可能值求期望

第三行：对上式进行简化，将−logpθ(x0)
−
�
�
�
�
�
(
�
0
)
项移到期望内部

其中 Eq
�
�
表示对分布 q(x1:T|x0)
�
(
�
1
:
�
|
�
0
)
 中的所有可能值求期望

第四行：−logpθ(x0)
−
�
�
�
�
�
(
�
0
)
 和 +logpθ(x0)
+
�
�
�
�
�
(
�
0
)
 相互抵消可得

令

 Let LVLB=Eq(x0:T)[logq(x1:T∣x0)pθ(x0:T)]≥−Eq(x0)logpθ(x0)
 Let 
�
V
L
B
=
�
�
(
�
0
:
�
)
[
log
⁡
�
(
�
1
:
�
∣
�
0
)
�
�
(
�
0
:
�
)
]
≥
−
�
�
(
�
0
)
log
⁡
�
�
(
�
0
)

所以 LVLB
�
�
�
�
就是我们的上界，我们要最小化它，接着进行变形

老规矩，上面整个推导总计九行，下面逐行解释下上面推导的每一行(纵使其他所有文章都不解释，本文也要给你解释的明明白白)

第一行，直接给出了 LVLB
�
V
L
B
 的定义，即计算概率分布 q
�
 和pθ
�
�
之间的对数比值的期望(注意，这是咱们的目标)

LVLB=Eq(x0:T)[logq(x1:T∣x0)pθ(x0:T)]
�
V
L
B
=
�
�
(
�
0
:
�
)
[
log
⁡
�
(
�
1
:
�
∣
�
0
)
�
�
(
�
0
:
�
)
]

考虑到七月SD课的一学员对第一步有疑问，再多解释一下
在第1步中，公式的意思是对于所有可能的的x0,x1,…,xT
�
0
,
�
1
,
…
,
�
�
路径，我们要计算中括号内部表达式的期望值

这里是一个联合分布，表示所有xt(t=0,1,...,T)
�
�
(
�
=
0
,
1
,
.
.
.
,
�
)
一起的分布，而期望是在这个分布下计算的

这意味着我们在所有这个分布支持的路径上平均这个中括号里的量
所以，如果​x0,x1,…,xT
�
0
,
�
1
,
…
,
�
�
都服从分布q，那么我们确实是在对中括号里的式子求整体期望

第二行，将条件概率 q(x1:T∣x0)
�
(
�
1
:
�
∣
�
0
)
 和联合概率 pθ(x0:T)
�
�
(
�
0
:
�
)
 展开为一系列条件概率的乘积
q(x1:T∣x0)=∏Tt=1q(xt∣xt−1)
�
(
�
1
:
�
∣
�
0
)
=
∏
�
=
1
�
�
(
�
�
∣
�
�
−
1
)

考虑到pθ(x0:T)
�
�
(
�
0
:
�
)
实际上就是pθ(x1:T∣x0)
�
�
(
�
1
:
�
∣
�
0
)
，所以有
pθ(x0:T)=pθ(x0)∏Tt=1pθ(xt|x0:t−1)=pθ(xT)∏Tt=1pθ(xt−1∣xt)
�
�
(
�
0
:
�
)
=
�
�
(
�
0
)
∏
�
=
1
�
�
�
(
�
�
|
�
0
:
�
−
1
)
=
�
�
(
�
�
)
∏
�
=
1
�
�
�
(
�
�
−
1
∣
�
�
)

然后把上述结果分别分别代入q(x1:T∣x0)
�
(
�
1
:
�
∣
�
0
)
和 pθ(x0:T)
�
�
(
�
0
:
�
)
，即可得到第二行的结果

Eq[log∏Tt=1q(xt∣xt−1)pθ(xT)∏Tt=1pθ(xt−1∣xt)]
�
�
[
log
⁡
∏
�
=
1
�
�
(
�
�
∣
�
�
−
1
)
�
�
(
�
�
)
∏
�
=
1
�
�
�
(
�
�
−
1
∣
�
�
)
]

第三行，将乘积转换为求和，并将pθ(xT)
�
�
(
�
�
)
项移到前面

Eq[−logpθ(xT)+∑Tt=1logq(xt∣xt−1)pθ(xt−1∣xt)]
�
�
[
−
log
⁡
�
�
(
�
�
)
+
∑
�
=
1
�
log
⁡
�
(
�
�
∣
�
�
−
1
)
�
�
(
�
�
−
1
∣
�
�
)
]

第四行，调整求和的范围，使其从2开始，从而达到将 t=1
�
=
1
 的项分离出来的目的

Eq[−logpθ(xT)+∑Tt=2logq(xt∣xt−1)pθ(xt−1∣xt)+logq(x1∣x0)pθ(x0∣x1)]
�
�
[
−
log
⁡
�
�
(
�
�
)
+
∑
�
=
2
�
log
⁡
�
(
�
�
∣
�
�
−
1
)
�
�
(
�
�
−
1
∣
�
�
)
+
log
⁡
�
(
�
1
∣
�
0
)
�
�
(
�
0
∣
�
1
)
]

第五行，将 t
�
 项的对数比值分解为两个对数比值的和，其中一个涉及xt−1
�
�
−
1
和 xt
�
�
，另一个涉及xt
�
�
和 x0
�
0
，相当于补了个x0
�
0

Eq[−logpθ(xT)+∑Tt=2log(q(xt−1∣xt,x0)pθ(xt−1∣xt)⋅q(xt∣x0)q(xt−1∣x0))+logq(x1∣x0)pθ(x0∣x1)]
�
�
[
−
log
⁡
�
�
(
�
�
)
+
∑
�
=
2
�
log
⁡
(
�
(
�
�
−
1
∣
�
�
,
�
0
)
�
�
(
�
�
−
1
∣
�
�
)
⋅
�
(
�
�
∣
�
0
)
�
(
�
�
−
1
∣
�
0
)
)
+
log
⁡
�
(
�
1
∣
�
0
)
�
�
(
�
0
∣
�
1
)
]

这里得着重解释下
把第四行的第二项的分子和分母都乘以q(xt−1|x0)
�
(
�
�
−
1
|
�
0
)
，即得
logq(xt∣xt−1)pθ(xt−1∣xt)=log(q(xt−1∣xt,x0)pθ(xt−1∣xt)⋅q(xt∣x0)q(xt−1∣x0))
log
⁡
�
(
�
�
∣
�
�
−
1
)
�
�
(
�
�
−
1
∣
�
�
)
=
log
⁡
(
�
(
�
�
−
1
∣
�
�
,
�
0
)
�
�
(
�
�
−
1
∣
�
�
)
⋅
�
(
�
�
∣
�
0
)
�
(
�
�
−
1
∣
�
0
)
)

这里面的关键是，即同时乘以q(xt−1|x0)
�
(
�
�
−
1
|
�
0
)
后，怎么就得到上式了呢，分母部分一目了然，直接乘上的q(xt−1|x0)
�
(
�
�
−
1
|
�
0
)
，但分子部分呢，明明应该是q(xt|xt−1)q(xt−1|x0)
�
(
�
�
|
�
�
−
1
)
�
(
�
�
−
1
|
�
0
)
，则就变成了这个呢：q(xt−1|xt,x0)q(xt|x0)
�
(
�
�
−
1
|
�
�
,
�
0
)
�
(
�
�
|
�
0
)
？好问题! 原因在于这两个式子是等价的，即(定义为等式1)

q(xt|xt−1)q(xt−1|x0)=q(xt−1|xt,x0)q(xt|x0)
�
(
�
�
|
�
�
−
1
)
�
(
�
�
−
1
|
�
0
)
=
�
(
�
�
−
1
|
�
�
,
�
0
)
�
(
�
�
|
�
0
)

为何等价呢，或者说上面这个等式1是怎么来的？其实也简单，因有

p(A, B, C) = p(A|B, C) p(B, C) = p(A|B, C) p(B|C) p(C)

p(A, B, C) = p(B|A, C) p(A, C) = p(B|A, C) p(A|C) p(C)

故有『下面五个等式先后依据：马尔科夫假设倒推、条件概率定义、分母中联合概率定义、分子中联合概率定义、分子分母同时约掉q(x0)
�
(
�
0
)
』

q(xt∣xt−1)=q(xt∣xt−1,x0)=q(xt,xt−1,x0)q(xt−1,x0)=q(xt−1∣xt,x0)q(xt∣x0)q(x0)q(xt−1,x0)=q(xt−1∣xt,x0)q(xt∣x0)q(x0)q(xt−1∣x0)q(x0)=q(xt−1∣xt,x0)q(xt∣x0)q(xt−1∣x0)
�
(
�
�
∣
�
�
−
1
)
=
�
(
�
�
∣
�
�
−
1
,
�
0
)
=
�
(
�
�
,
�
�
−
1
,
�
0
)
�
(
�
�
−
1
,
�
0
)
=
�
(
�
�
−
1
∣
�
�
,
�
0
)
�
(
�
�
∣
�
0
)
�
(
�
0
)
�
(
�
�
−
1
,
�
0
)
=
�
(
�
�
−
1
∣
�
�
,
�
0
)
�
(
�
�
∣
�
0
)
�
(
�
0
)
�
(
�
�
−
1
∣
�
0
)
�
(
�
0
)
=
�
(
�
�
−
1
∣
�
�
,
�
0
)
�
(
�
�
∣
�
0
)
�
(
�
�
−
1
∣
�
0
)

第六行，将第五行的中间项一分为二，即拆分为两个求和项

Eq[−logpθ(xT)+∑Tt=2logq(xt−1∣xt,x0)pθ(xt−1∣xt)+∑Tt=2logq(xt∣x0)q(xt−1∣x0)+logq(x1∣x0)pθ(x0∣x1)]
�
�
[
−
log
⁡
�
�
(
�
�
)
+
∑
�
=
2
�
log
⁡
�
(
�
�
−
1
∣
�
�
,
�
0
)
�
�
(
�
�
−
1
∣
�
�
)
+
∑
�
=
2
�
log
⁡
�
(
�
�
∣
�
0
)
�
(
�
�
−
1
∣
�
0
)
+
log
⁡
�
(
�
1
∣
�
0
)
�
�
(
�
0
∣
�
1
)
]

第七行，将第五行中间部分得到的两个求和项的第二个求和项的最后一项t=T
�
=
�
分离出来，说白了，将第二个求和项的范围调整为从1到T−1
�
−
1
，啥意思呢
首先，第五行中间部分的两个求和项可以表示为

∑Tt=2logq(xt−1∣xt,x0)pθ(xt−1∣xt)+∑Tt=2logq(xt∣x0)q(xt−1∣x0)=∑Tt=2[logq(xt−1∣xt,x0)pθ(xt−1∣xt)+logq(xt∣x0)q(xt−1∣x0)]
∑
�
=
2
�
log
⁡
�
(
�
�
−
1
∣
�
�
,
�
0
)
�
�
(
�
�
−
1
∣
�
�
)
+
∑
�
=
2
�
log
⁡
�
(
�
�
∣
�
0
)
�
(
�
�
−
1
∣
�
0
)
=
∑
�
=
2
�
[
log
⁡
�
(
�
�
−
1
∣
�
�
,
�
0
)
�
�
(
�
�
−
1
∣
�
�
)
+
log
⁡
�
(
�
�
∣
�
0
)
�
(
�
�
−
1
∣
�
0
)
]

接下来，关键的一步在于，上面中括号里的第二个求和项在求和过程中相邻两项会相互抵消 「依据：logaM/N=logaM−logaN
�
�
�
�
�
/
�
=
�
�
�
�
�
−
�
�
�
�
�
」
具体地，当t=k
�
=
�
时的q(xk∣x0)
�
(
�
�
∣
�
0
)
会和当t=k+1
�
=
�
+
1
时的q(xk∣x0)
�
(
�
�
∣
�
0
)
相互抵消，这样的抵消会发生在每一对相邻的项上，从2到T−1
�
−
1
，最后，只剩下t=T
�
=
�
和t=2
�
=
2
时的两项，即logq(xT∣x0)q(x1∣x0)
log
⁡
�
(
�
�
∣
�
0
)
�
(
�
1
∣
�
0
)

从而得到最终整个第7行所示的结果，如下

Eq[−logpθ(xT)+∑Tt=2logq(xt−1∣xt,x0)pθ(xt−1∣xt)+logq(xT∣x0)q(x1∣x0)+logq(x1∣x0)pθ(x0∣x1)]
�
�
[
−
log
⁡
�
�
(
�
�
)
+
∑
�
=
2
�
log
⁡
�
(
�
�
−
1
∣
�
�
,
�
0
)
�
�
(
�
�
−
1
∣
�
�
)
+
log
⁡
�
(
�
�
∣
�
0
)
�
(
�
1
∣
�
0
)
+
log
⁡
�
(
�
1
∣
�
0
)
�
�
(
�
0
∣
�
1
)
]

第八行，上一行第7行总共4个带log的项，把最后两个log项拆开成4个式子，抵消两个，还分别剩一个logq(xT|x0)
�
�
�
�
(
�
�
|
�
0
)
、一个−logpθ(x0|x1)
−
�
�
�
�
�
(
�
0
|
�
1
)
，然后logq(xT|x0)
�
�
�
�
(
�
�
|
�
0
)
与最初4项中的第1项−logpθ(xT)
−
�
�
�
�
�
(
�
�
)
合并，即可得到整个第八行的结果

Eq[logq(xT∣x0)pθ(xT)+∑Tt=2logq(xt−1∣xt,x0)pθ(xt−1∣xt)−logpθ(x0∣x1)]
�
�
[
log
⁡
�
(
�
�
∣
�
0
)
�
�
(
�
�
)
+
∑
�
=
2
�
log
⁡
�
(
�
�
−
1
∣
�
�
,
�
0
)
�
�
(
�
�
−
1
∣
�
�
)
−
log
⁡
�
�
(
�
0
∣
�
1
)
]

第九行，将最后一项中的负号移到对数里面，并将整个表达式重写为一系列 KL 散度项的和，这些项分别为 LT
�
�
、Lt−1
�
�
−
1
 和 L0
�
0

Eq[DKL(q(xT∣x0)∥pθ(xT))LT+∑Tt=2DKL(q(xt−1∣xt,x0)∥pθ(xt−1∣xt))Lt−1−logpθ(x0∣x1)L0]
�
�
[
�
K
L
(
�
(
�
�
∣
�
0
)
‖
�
�
(
�
�
)
)
⏟
�
�
+
∑
�
=
2
�
�
K
L
(
�
(
�
�
−
1
∣
�
�
,
�
0
)
‖
�
�
(
�
�
−
1
∣
�
�
)
)
⏟
�
�
−
1
−
log
⁡
�
�
(
�
0
∣
�
1
)
⏟
�
0
]

还是考虑到七月SD课的一学员对第一步有疑问，再多解释一下
对于第9步，KL散度是一种衡量两个概率分布p和q之间差异的方法，这里的公式是在计算q和p之间的KL散度的期望值
尽管KL散度自己就是一个期望值的形式，但这里的Eq
�
�
指的是在分布q下对KL散度本身求期望
这意味着你不仅计算了q和p之间的差异，而且你要考虑所有q分布可能产生的不同序列x0,x1,…,xT
�
0
,
�
1
,
…
,
�
�
，对这些序列的KL散度进行平均

最后得到的表达式表示了最初我们想求解的LVLB
�
V
L
B
最终是一系列 KL 散度项之和，我们可以利用这个结果进行参数优化，使得两个概率分布之间的差异最小

2.3.2 拉近估计分布pθ(xt−1∣xt)
�
�
(
�
�
−
1
∣
�
�
)
 和真实后验分布q(xt−1|xt,x0)
�
(
�
�
−
1
|
�
�
,
�
0
)
 
对于上面公式最后第九行得到的结果

首先，LT
�
�
是和优化无关的(由于前向过程 q
�
 没有可学习参数，而 xT
�
�
 则是纯高斯噪声，因此 LT
�
�
 可以当做常量忽略)，所以不用管，只用看右边的Lt−1
�
�
−
1
然后，Lt−1
�
�
−
1
 是KL散度，则可以看做拉近估计分布pθ(xt−1∣xt)
�
�
(
�
�
−
1
∣
�
�
)
 和真实后验分布q(xt−1|xt,x0)
�
(
�
�
−
1
|
�
�
,
�
0
)
 这两个分布之间的距离： 
对于真实后验分布 q(xt−1|xt,x0)
�
(
�
�
−
1
|
�
�
,
�
0
)
 ，我们已经在上一节2.2.2节推导出其解析形式，这是一个高斯分布，其均值和方差为
μ~t=1αt√(xt−1−αt1−α¯t√ϵ0)
�
~
�
=
1
�
�
(
�
�
−
1
−
�
�
1
−
�
¯
�
�
0
)

β~t=1−α¯t−11−α¯t⋅βt
�
~
�
=
1
−
�
¯
�
−
1
1
−
�
¯
�
⋅
�
�

对于估计分布 pθ(xt−1∣xt)
�
�
(
�
�
−
1
∣
�
�
)
，其是我们网络期望拟合的目标分布，也是一个高斯分布，均值用网络估计，方差被设置为了一个和 βt
�
�
 有关的常数
pθ(xt−1∣xt)=N(xt−1;μθ(xt,t),Σθ(xt,t))
�
�
(
�
�
−
1
∣
�
�
)
=
�
(
�
�
−
1
;
�
�
(
�
�
,
�
)
,
�
�
(
�
�
,
�
)
)
考虑到，如果有两个分布 p,q 都是高斯分布，则他们的KL散度为
KL(p,q)=logσ2σ1+σ21+(μ1−μ2)22σ22−12
�
�
(
�
,
�
)
=
log
⁡
�
2
�
1
+
�
1
2
+
(
�
1
−
�
2
)
2
2
�
2
2
−
1
2

然后因为这两个分布的方差全是常数，和优化无关『说白了，去掉方差σ
�
相关的项，只留下均值相关的 (μ1−μ2)2
(
�
1
−
�
2
)
2
』，所以其实优化目标就是两个分布均值的二范数，从而可得
“ 怎么来的？我再细致解释下，对于这个公式而言
Lt=Eq[∥μ~t(xt,x0)−μθ(xt,t)∥2]
�
�
=
�
�
[
‖
�
~
�
(
�
�
,
�
0
)
−
�
�
(
�
�
,
�
)
‖
2
]

这里的 Eq
�
�
是在分布 q
�
 下的期望，当我们μt~
�
�
~
的表达式代入后，得到：
Lt=Eq[∥∥∥1αt√(xt−βt1−α¯t√ϵ)−μθ(xt,t)∥∥∥2]
�
�
=
�
�
[
‖
1
�
�
(
�
�
−
�
�
1
−
�
¯
�
�
)
−
�
�
(
�
�
,
�
)
‖
2
]

在这个式子中，ϵ
�
 是一个服从标准正态分布的随机变量，而xt
�
�
则取决于x0
�
0
和ϵ
�
。因此，这个期望 Eq
�
�
实际上是在 x0
�
0
和ϵ
�
 的联合分布下的期望(在 x0
�
0
和ϵ
�
 的所有可能值上取平均)，于是我们得到：
Lt=Ex0,ϵ[∥∥∥1αt√(xt(x0,ϵ)−βt1−α¯t√ϵ)−μθ(xt(x0,ϵ),t)∥∥∥2]
�
�
=
�
�
0
,
�
[
‖
1
�
�
(
�
�
(
�
0
,
�
)
−
�
�
1
−
�
¯
�
�
)
−
�
�
(
�
�
(
�
0
,
�
)
,
�
)
‖
2
]

Ex0,ϵ
�
�
0
,
�
 代表就是在 x0
�
0
和ϵ
�
 的联合分布下的期望，ϵ
�
 依然是从标准正态分布N(0,1)
�
(
0
,
1
)
中采样的噪声 ”

这个时候我们可以直接整个网络出来直接学习 uθ(xt,t)
�
�
(
�
�
,
�
)
 ，然后uθ(xt,t)
�
�
(
�
�
,
�
)
再去预测

因为xt
�
�
是uθ
�
�
 的输入，其它的量都是常数，所以其中的未知量其实只有ϵ
�
，所以我们干脆把需要学习的uθ(xt,t)
�
�
(
�
�
,
�
)
定义成：
μθ(xt,t)=1αt√(xt−βt1−α¯t√ϵθ(xt,t))
�
�
(
�
�
,
�
)
=
1
�
�
(
�
�
−
�
�
1
−
�
¯
�
�
�
(
�
�
,
�
)
)

也就是说，不用网络预测ut~(xt,x0)
�
�
~
(
�
�
,
�
0
)
，而是用网络ϵθ(xt,t)
�
�
(
�
�
,
�
)
先预测噪声ϵ
�
 (注意，这是个关键步骤)，然后把预测出来的噪声带入到定义好的表达式中去计算出预测的均值即可
所以，最终把μθ(xt,t)
�
�
(
�
�
,
�
)
这个公式，代入到步骤3得到的公式中，可得：

经过这样一番推导之后就是个 L2 loss，网络的输入是一张和噪声线性组合的图片，然后要通过ϵθ(α¯t−−√x0+1−α¯t−−−−−√ϵ,t)
�
�
(
�
¯
�
�
0
+
1
−
�
¯
�
�
,
�
)
估计出来这个噪声
由上可知，DDPM的关键是训练 ϵθ(xt,t)
�
�
(
�
�
,
�
)
模型，使其预测的 ϵ^
�
^
与真实用于破坏的ϵ
�
 相近，用L2距离刻画相近程度就好，总之，我们的Loss就是如下公式『相当于训练时，网络输入为 xt
�
�
(由 x0
�
0
 和噪声 ϵ
�
 线性组合而成) 和时刻t
�
 ，输出要尽可能的拟合输入的噪声 ϵ
�
 (通过L2 loss约束)』

下图可以总结噪声估计模型的训练过程 (依然是经典的那一套：对比预测噪声predicted noise与真实噪声true noise ε
�
之间的差距 建loss 反向传播，训练好之后，好预测噪声，毕竟模糊的图片减掉噪声 不就得到清晰的图片了么)

而整个训练过程可如下图描述

DDPM论文中对应的伪代码为

2.4 如何通过训练好的DDPM生成图片
通过上文2.2节的最后，我们得知

从最终得到的结果可以看出，在给定 x0
�
0
 的条件下

后验条件高斯分布的均值只和超参数αt
�
�
、xt
�
�
、ϵ0
�
0
有关，即

μ(xt,x0)=1αt√(xt−1−αt1−α¯t√ϵ0)
�
(
�
�
,
�
0
)
=
1
�
�
(
�
�
−
1
−
�
�
1
−
�
¯
�
�
0
)

方差只与超参数α
�
有关，即

βt~=1A=1/(αtβt+11−α¯t−1)=1/(αt−α¯t+βtβt(1−α¯t−1))=1−α¯t−11−α¯t⋅βt
�
�
~
=
1
�
=
1
/
(
�
�
�
�
+
1
1
−
�
¯
�
−
1
)
=
1
/
(
�
�
−
�
¯
�
+
�
�
�
�
(
1
−
�
¯
�
−
1
)
)
=
1
−
�
¯
�
−
1
1
−
�
¯
�
⋅
�
�

从而通过以上的方差和均值，我们就得到了q(xt−1|xt,x0)
�
(
�
�
−
1
|
�
�
,
�
0
)
的解析形式

通过2.2节的最后，我们估计到了噪声估测模型 ϵθ(xt,t)
�
�
(
�
�
,
�
)
 ，接下来要生成模型就很简单了。从N(0,1)中随机生成一个噪声作为XT
�
�
，然后再用该模型逐步从估测噪声，并用去噪公式逐渐恢复到 x0
�
0
 即可，见如下伪代码

用通俗的语言来说，如何去噪生成清晰图片呢？如上述第4行代码所述，比较模糊的图片xt
�
�
 减掉噪声估计器预测出来的噪声 (只是实际操作时，减掉的是与一个常数1−αt1−α¯t√
1
−
�
�
1
−
�
¯
�
相乘之后的噪声)，得到的结果乘以1/at−−√
1
/
�
�
之后，最后再加上一个带σt
�
�
的z
�

换言之，推理时，我们从各项独立的高斯分布xt
�
�
 开始，一共 T
�
 步，每一步其实都是用了一次reparameterization trick

每一步具体来说，我们有了 xt
�
�
 ， 想要得到 xt−1
�
�
−
1
，因为我们之前逆扩散过程建模有：  

pθ(xt−1∣xt)=N(xt−1;μθ(xt,t),Σ2θ(xt,t))
�
�
(
�
�
−
1
∣
�
�
)
=
�
(
�
�
−
1
;
�
�
(
�
�
,
�
)
,
Σ
�
2
(
�
�
,
�
)
)

=N(xt−1;1αt√(xt−βt1−α¯t√ϵθ(xt,t)),βt)
=
�
(
�
�
−
1
;
1
�
�
(
�
�
−
�
�
1
−
�
¯
�
�
�
(
�
�
,
�
)
)
,
�
�
)

所以由reparameterization trick我们有：

每一轮都这样迭代，最终就得到了生成的图片，如下图所示

第三部分 DETR：首次通过结合CNN+Transformer端对端解决object detection
注，本小节的内容主要参考科技猛兽此文的相关部分

一般目标检测的任务是预测一系列的Bounding Box的坐标以及Label，而大多数检测器的具体做法是

要么基于proposal，比如RCNN系列的工作，类似Faster R-CNN、Mask R-CNN
要么基于anchor，比如YOLO
把问题构建成为一个分类和回归问题来间接地完成这个任务，但最后都会生成很多个预测框(确定框的坐标及框内是什么物体)，从而不可避免的出现很多冗余的框，而要去除这些冗余的框，则都需要做一个NMS(non-maximum suppersion，非极大值抑制)的后处理(使得最后调参不易、部署不易)，所以如果要是有一个端对端的模型，不需要做NMS之类的后处理 也不需要太多先验知识则该有多好

而通过论文《End-to-End Object Detection with Transformers》提出的DETR则满足了大家这个期待，其取代了现在的模型需要手工设计的工作，效果不错且可扩展性强(在DETR上加个专用的分割头便可以做全景分割)，其解决的方法是把检测问题看做是一个集合预测的问题(即set prediction problem，说白了，各种预测框本质就是一个集合)，其基本流程如下图所示

CNN抽特征且拉直
全局建模，给到transformer-encoder去进一步学习全局信息
通过借助Transformer中的的self-attention机制，可以显式地对一个序列中的所有elements两两之间的interactions进行建模或交互，如此就知道了图片中哪块是哪个物体，从而对于同一个物体只需出一个预测框即可
接着通过不带掩码机制的transformer-decoder生成很多预测框
注意是并行预测(即并行出框，而不是像原始transformer预测下一个token时一个一个往外蹦)
相当于一次性生成 N
�
 个box prediction，其中 N
�
是一个事先设定的远远大于image中object个数的一个整数(比如100)
预测框和真实框做二分图匹配
最后通过bipartite matching loss的方法，基于预测的100个boxex和ground truth boxes的二分图做匹配，计算loss的大小，从而使得预测的box的位置和类别更接近于ground truth
当然，这第4步更多是做模型训练的时候用，如果模型训练好了去做推理时，该第4步就不需要了，可以直接在预测的100个框中设定个阈值，比如置信度大于0.7的预测框保留下来，其他视为背景物体而 舍弃
3.1 DETR整体结构：backbone，encoder，decoder和FFN
3.1.1 DETR结构之前两部分：backbone与encoder
更细致的讲，DETR整体结构可以分为四个部分：backbone，encoder，decoder和FFN，如下图所示

对于前两部分而言

一开始的backbone面对的是 ximg∈B×3×H0×W0
�
�
�
�
∈
�
×
3
×
�
0
×
�
0
维的图像，首先把它转换为f∈RB×C×H×W
�
∈
�
�
×
�
×
�
×
�
维的feature map（一般来说，通道数C=2048
�
=
2048
或256，H=H032
�
=
�
0
32
，W=W032
�
=
�
0
32
）

然后由于encoder的输入是f∈RB×C×H×W
�
∈
�
�
×
�
×
�
×
�
维的feature map，故正式输入encoder之前还需要依次进行以下过程(图源：科技猛兽)：

通道数压缩(其实就是降维操作)：用 1×1 convolution处理将通道数channels数量从 C
�
 压缩到d
�
，即得到z0∈RB×d×H×W
�
0
∈
�
�
×
�
×
�
×
�
维的新feature map
转化为序列化数据：将空间的维度(高H
�
和宽W
�
)压缩为一个维度HW
�
�
，即把上一步得到的z0∈RB×d×H×W
�
0
∈
�
�
×
�
×
�
×
�
(d=256
�
=
256
)维的feature map通过reshape成(HW,B,256
�
�
,
�
,
256
)维的feature map
这步相当于把编码矩阵的维度是 (B,256,H,W)
(
�
,
256
,
�
,
�
)
 ，序列化成维度为 (HW,B,256)
(
�
�
,
�
,
256
)
维的张量
位置编码：在上一步得到了z0∈RB×d×H×W
�
0
∈
�
�
×
�
×
�
×
�
维的feature map之后，再对 z0∈RB×d×H×W
�
0
∈
�
�
×
�
×
�
×
�
维的feature map做positional encoding，最后也做下reshape：高H
�
和宽W
�
压缩为一个维度HW
�
�
，使得其与上面input embedding维度是一致的
3.1.2 详解DETR的位置编码
对于上节第三步的位置编码，再好好解释说明下
首先，通过此文《从零实现Transformer的简易版与强大版：从300多行到3000多行》的1.1.2节可知，原始transformer中的Positional Encoding的表达式为：

其中，d
�
 就是这个 d×HW
�
×
�
�
维的feature map的第一维，pos∈[1,HW]
�
�
�
∈
[
1
,
�
�
]
表示token在sequence中的位置，sequence的长度是HW
�
�
，例如第一个token 的 pos=0
�
�
�
=
0
，第二个token的pos=1
�
�
�
=
1

i
�
，或者准确意义上是2i
2
�
和2i+1
2
�
+
1
表示了Positional Encoding的维度，i
�
 的取值范围是：[0,...,dmodel/2]
[
0
,
.
.
.
,
�
�
�
�
�
�
/
2
]
，所以当 pos
�
�
�
 为1时，对应的Positional Encoding可以写成(注意到d=256
�
=
256
)：

PE1=[sin(1100000256),cos(1100000256),sin(1100002256),cos(1100002256),sin(1100004256),cos(1100004256),...,sin(110000255256),cos(110000255256)]
�
�
1
=
[
�
�
�
(
1
10000
0
256
)
,
�
�
�
(
1
10000
0
256
)
,
�
�
�
(
1
10000
2
256
)
,
�
�
�
(
1
10000
2
256
)
,
�
�
�
(
1
10000
4
256
)
,
�
�
�
(
1
10000
4
256
)
,
.
.
.
,
�
�
�
(
1
10000
255
256
)
,
�
�
�
(
1
10000
255
256
)
]

其次，DETR与原版transformer中的位置编码有两点不同

第一点不同的是，原版Transformer只考虑 x
�
 方向的位置编码，但是DETR考虑了 xy
�
�
 方向的位置编码，因为图像特征是2-D特征。采用的依然是 sin cos 模式，但是需要考虑 xy
�
�
 两个方向。不是类似vision transoformer做法简单的将其拉伸为 d×HW
�
×
�
�
，然后从 [1,HW]
[
1
,
�
�
]
进行长度为256的位置编码，而是考虑了 xy
�
�
 方向同时编码，每个方向各编码128维向量，这种编码方式更符合图像特点

Positional Encoding的输出张量是：(B,d,H,W)
(
�
,
�
,
�
,
�
)
,d=256
�
=
256
 ，其中 d
�
 代表位置编码的长度， H,W
�
,
�
代表张量的位置。意思是说，这个特征图上的任意一个点 (H1,W1)
(
�
1
,
�
1
)
有个位置编码，这个编码的长度是256，其中，前128维代表 H1
�
1
 的位置编码，后128维代表W1
�
1
 的位置编码
a)PE(px,2i)=sin(posx/100002i/128) b)PE(p0,2i+1)=cos(posx/100002i/128)c)PE(posy,2i)=sin(posy/100002i/128)d)PE(posy,2i+1)=cos(posy/100002i/128)
�
)
�
�
(
�
�
,
2
�
)
=
sin
⁡
(
pos
�
⁡
/
10000
2
�
/
128
)
 
�
)
�
�
(
�
0
,
2
�
+
1
)
=
cos
⁡
(
pos
�
⁡
/
10000
2
�
/
128
)
�
)
�
�
(
�
�
�
�
,
2
�
)
=
sin
⁡
(
pos
�
⁡
/
10000
2
�
/
128
)
�
)
�
�
(
�
�
�
�
,
2
�
+
1
)
=
cos
⁡
(
pos
�
⁡
/
10000
2
�
/
128
)

假设你想计算任意一个位置(posx,posy)
(
�
�
�
�
,
�
�
�
�
)
 ，posx∈[1,HW]
�
�
�
�
∈
[
1
,
�
�
]
，posy∈[1,HW]
�
�
�
�
∈
[
1
,
�
�
]
的Positional Encoding
→
→
 先把posx
�
�
�
�
代入上面4个公式中的a
�
 式和 b
�
 式可以计算得到128维的向量，它代表 posx
�
�
�
�
 的位置编码
→
→
 再把 posy
�
�
�
�
 代入上面4个公式中的 c
�
 式和d
�
 式可以计算得到128维的向量，它代表 posy
�
�
�
�
 的位置编码
→
→
 把这2个128维的向量拼接起来，就得到了一个256维的向量，它代表(posx,posy)
(
�
�
�
�
,
�
�
�
�
)
的位置编码
从而计算所有位置的编码，就得到了(256,H,W)
(
256
,
�
,
�
)
的张量，代表这个batch的位置编码
第二点不同的是，原版Transformer只在Encoder之前使用了Positional Encoding，而且是在输入上进行Positional Encoding，再把输入经过transformation matrix变为Query，Key和Value这几个张量
但是DETR在Encoder的每一个Multi-head Self-attention之前都使用了Positional Encoding，且只对Query和Key使用了Positional Encoding，即：只把维度为(HW,B,256)
(
�
�
,
�
,
256
)
维的位置编码与维度为(HW,B,256)
(
�
�
,
�
,
256
)
维的Query和Key相加，而不与Value相加

下图为DETR的Transformer的详细结构，读者可以对比下原版Transformer的结构 可以发现，除了Positional Encoding设置的不一样外，Encoder其他的结构是一致的。每个Encoder Layer包含一个multi-head self-attention 的module和一个前馈网络Feed Forward Network

 所以，了解了DETR的位置编码之后，你应该明白了其实input embedding和位置编码维度其实是一样的，从而也就可以直接相加，使得Encoder最终输出的是(HW,b,256)
(
�
�
,
�
,
256
)
维的编码矩阵Embedding，按照原版Transformer的做法，把这个东西给Decoder

3.1.3 DETR结构的后两部分：decoder和FFN
通过上节最后的对比图，可知DETR的Decoder和原版Transformer的decoder也是不太一样的

对于原版Transformer，其decoder的最后一个框(上节最后对比图的左图右上角所示)：output probability，代表我们一次只产生一个单词的softmax，根据这个softmax得到这个单词的预测结果，即：predicts the output sequence one element at a time
不同的是，DETR的Transformer Decoder是一次性处理全部的object queries(上节最后对比图的右图右上角所示)，即一次性输出全部的predictions(而不像原始的Transformer是auto-regressive的，从左到右一个词一个词地输出)，即：decodes the N objects in parallel at each decoder layer
至于DETR的Decoder主要有两个输入：

第一个输入是Transformer Encoder输出的Embedding与 position encoding(在下图右侧第二个multi-head self-attention处)相加之后 给到K
�

其中的Embedding就是上文提到的(HW,B,256)
(
�
�
,
�
,
256
)
的编码矩阵
第二个输入是Object queries
所谓Object queries是一个维度为(100,b,256)
(
100
,
�
,
256
)
维的张量，数值类型是nn.Embedding(意味着是可以学习的)。Object queries矩阵内部通过学习建模了100个物体之间的全局关系，例如房间里面的桌子旁边(A类)一般是放椅子(B类)，而不会是放一头大象(C类)，那么在推理时候就可以利用该全局注意力更好的进行解码预测输出

关于上图右侧第1个multi-head self-attention的Q K V
如上图所示，它的Q K是这么来的：Decoder原本的输入一开始初始化成维度为(100,b,256)
(
100
,
�
,
256
)
维的全部元素都为0的张量，然后和Object queries加在一起之后充当第1个multi-head self-attention的Query和Key，至于Value则是Decoder原本的输入，也就是全0的张量

关于上图右侧第2个multi-head self-attention的Q K V
它的Key和Value来自Encoder的输出张量，维度为(hw,b,256)
(
ℎ
�
,
�
,
256
)
，其中Key值还进行位置编码(正如上面第一个输入所述)
至于其Query值一部分来自第1个Add and Norm的输出，维度为(100,b,256)
(
100
,
�
,
256
)
的张量，另一部分来自Object queries，充当可学习的位置编码
所以，第2个multi-head self-attention的Key和Value的维度为(hw,b,256)
(
ℎ
�
,
�
,
256
)
，而Query的维度为(100,b,256)
(
100
,
�
,
256
)
每个Decoder的输出维度为(1,b,100,256)
(
1
,
�
,
100
,
256
)
，送入后面的前馈网络

到这里你会发现：Object queries充当的其实是位置编码的作用，只不过它是可以学习的位置编码，所以，我们对Encoder和Decoder的每个self-attention的Query和Key的位置编码做个归纳，如下图所示，Value没有位置编码

3.1.4 损失函数部分解读
得到了Decoder的输出以后，如前文所述，应该是输出维度为(b,100,256)
(
�
,
100
,
256
)
的张量。接下来要送入2个前馈网络FFN得到class和Bounding Box(如再度引用的下图所示)，它们会得到 N=100
�
=
100
 个预测目标 包含类别和Bounding Box(当然这个100肯定是大于图中的目标总数的，如果不够100，则采用背景填充

所以，DETR输出张量的维度为(b,100,class+1)
(
�
,
100
,
�
�
�
�
�
+
1
)
，和(b,100,4)
(
�
,
100
,
4
)
『对应COCO数据集来说：class+1=92
�
�
�
�
�
+
1
=
92
 ， 4 指的是每个预测目标归一化的(cx,cy,w,h)
(
�
�
,
�
�
,
�
,
ℎ
)
，归一化就是除以图片宽高进行归一化』，对于这两个维度

前者代表分类分支，指100个预测框的类型
后者代表回归分支，指100个预测框的Bounding Box(仅仅计算有物体位置，背景集合忽略)
但是读者可能会有疑问：预测框和真值是怎么一一对应的？换句话说：你怎么知道第47个预测框对应图片里的狗，第88个预测框对应图片里的车？..

这就需要用到经典的双边匹配算法了，也就是常说的匈牙利算法，该算法广泛应用于最优分配问题
一幅图片，我们把第 i
�
 个物体的真值表达为 yi=(ci,bi)
�
�
=
(
�
�
,
�
�
)
，其中，ci
�
�
表示它的 class ，bi
�
�
 表示它的 Bounding Box
然后定义y^={y^i}Ni=1
�
^
=
{
�
^
�
}
�
=
1
�
为网络输出的N
�
 个预测值
对于第i
�
 个真值GT
�
�
 ，σ(i)
�
(
�
)
为匈牙利算法得到的与真值GTi
�
�
�
 对应的预测值prediction的索引，举个例子，比如 i=3
�
=
3
，σ(i)=18
�
(
�
)
=
18
 ，意思就是：与第3个真值对应的预测值是第18个
那如何根据匈牙利算法找到与每个真值对应的预测值到底是哪个呢？
σ^=argminσ∈ΣN∑NiLmatch (yi,y^σ(i))
�
^
=
arg
⁡
min
�
∈
Σ
�
∑
�
�
�
match 
(
�
�
,
�
^
�
(
�
)
)

对于某一个真值yi
�
�
 ，假设已经找到这个真值对应的预测值y^σ(i)
�
^
�
(
�
)
，这里的ΣN
Σ
�
是所有可能的排列(代表从真值索引到预测值索引的所有的映射)，然后用Lmatch
�
�
�
�
�
ℎ
最小化yi
�
�
和y^σ(i)
�
^
�
(
�
)
的距离
这个Lmatch
�
�
�
�
�
ℎ
具体是：
−1{ci≠∅}p^σ(i)(ci)+1{ci≠∅}Lbox(bi,b^σ(i))
−
1
{
�
�
≠
∅
}
�
^
�
(
�
)
(
�
�
)
+
1
{
�
�
≠
∅
}
�
�
�
�
(
�
�
,
�
^
�
(
�
)
)

意思是：假设当前从真值索引到预测值索引的所有映射为σ
�
，对于图片中的每个真值 i
�

→
→
  先找到对应的预测值σ(i)
�
(
�
)
，再看看分类网络的结果p^σ(i)(ci)
�
^
�
(
�
)
(
�
�
)
，取反作为 Lmatch
�
�
�
�
�
ℎ
的第1部分
→
→
  再计算回归网络的结果b^σ(i)
�
^
�
(
�
)
与真值的 Bounding Box 的差异，即 Lbox(bi,b^σ(i))
�
�
�
�
(
�
�
,
�
^
�
(
�
)
)
，作为Lmatch
�
�
�
�
�
ℎ
的第2部分
所以，可以使得Lmatch
�
�
�
�
�
ℎ
最小的排列 σ^
�
^
 就是我们要找的排列，即：对于图片中的每个真值i
�
来讲，σ^(i)
�
^
(
�
)
就是这个真值所对应的预测值的索引
上述这个匈牙利算法的过程与Anchor或Proposal有异曲同工的地方，只是此时我们找的是一对一匹配
接下来就是使用上一步得到的排列σ^
�
^
，计算匈牙利损失
LHungarian (y,y^)=∑Ni=1[−logp^σ^(i)(ci)+1{ci≠∅}Lbox(bi,b^σ^(i))]
�
Hungarian 
(
�
,
�
^
)
=
∑
�
=
1
�
[
−
log
⁡
�
^
�
^
(
�
)
(
�
�
)
+
1
{
�
�
≠
∅
}
�
�
�
�
(
�
�
,
�
^
�
^
(
�
)
)
]

式中的Lbox
�
�
�
�
具体为：

最常用的L1
�
1
 loss对于大小 Bounding Box 会有不同的标度，即使它们的相对误差是相似的。为了缓解这个问题，作者使用了L1
�
1
 loss和广义IoU损耗Liou
�
�
�
�
的线性组合，它是比例不变的
Hungarian意思就是匈牙利，也就是前面的Lmatch
�
�
�
�
�
ℎ
，上述意思是需要计算M
�
个 GTBounding Box 和N
�
个输预测出集合两两之间的广义距离，距离越近表示越可能是最优匹配关系，也就是两者最密切，广义距离的计算考虑了分类分支和回归分支
第四部分 从ViT到MAE
4.1 Vision Transformer：用标准的Transformer直接干CV任务
继DETR、DDPM之后，Google于2020年10月发布的这篇论文《An Image is Worth 16x16 Words: Transformers for Image Recognition at Scale》提出的ViT彻底引燃了多模态的火热，更是直接挑战了此前CNN在视觉领域长达近10年的绝对统治地位

这个工作是怎么一步步出来的呢？自从Google在2017年发布的transformer在NLP领域大杀四方的时候，就一直不断有人想把如此强大且充满魔力的transformer用到CV领域中，但前路曲折啊

一开始面对的问题就是，当把transformer中对NLP的各个token之间两两互相做相似度计算的self-attention引入到图片各个像素点之间两两做self-attention时，你会发现计算复杂度瞬间爆炸(transformer 的计算复杂度是序列长度 n 的 平方 O(n2)
�
(
�
2
)
)
原因很简单，一句话才多少个token(顶多几百而已)，但一张图片呢？比如一张像素比较低的224×224
224
×
224
分辨率的图片，就已经达到了 50176 个像素点，再考虑到RGB三个维度，直接就是15万起步，你品..
兴高采烈之下理想遇挫，才发现没那么简单，咋办呢，那就降低序列长度呗。比如要么像CVPR Wang et al. 2018的工作把网络中间的特征图当做transformer的输入，毕竟ResNet 50 最后一个 stage, res4 的 feature map也就14×14=196
14
×
14
=
196
的大小，要么借鉴CNN的卷积机制用一个局部窗口去抓图片的特征，从而降低图片的复杂度，但这一系列工作虽然逻辑上通畅，但因为硬件上无法加速等，导致模型没法太大
总之，之前的工作要么把CNN和self-attention结合起来，要么把self-attention取代CNN，但都没取得很好的扩展效果，看来得再次冲击transformer for CV (为何如此执着？还不是因为理想实在过于美丽，要不然人到中年 每每听到beyond的“原谅我这一生不羁放纵爱自由”便如此共鸣强烈)，但还是得回头开头的老大难问题：如何处理图片的复杂度

于此，Vision Transformer(ViT)来了：不再一个一个像素点的处理，而是把整个图片切分成一个个图片块(比如分割为九宫格)，这些图片块作为transformer的输入
下面，我们来仔细探究下ViT到底是怎么做的

4.1.1 ViT的架构：Embedding层 + Transformer Encoder + MLP Head
简单而言，Vision Transformer(ViT)的模型框架由三个模块组成：

Embedding层（线性投射层Linear Projection of Flattened Patches）
以ViT_base_patch16为例，一张224 x 224的图片先分割成 16 x 16 的 patch ，很显然会因此而存在 (224×224/16×16)2=196
(
224
×
224
/
16
×
16
)
2
=
196
 个 patch(这个patch数如果泛化到一般情况就是图片长宽除以patch的长宽，即N=HW/P2
�
=
�
�
/
�
2
)
且图片的长宽由原来的224  x 224 变成：14  x 14(因为224/16 = 14)

你可能还没意识到这个操作的价值，这相当于把图片需要处理的像素单元从5万多直接降到了14x14 = 196个像素块，如果一个像素块当做一个token，那针对196个像素块/token去做self-attention不就轻松多了么(顺带提一句，其实在ViT之前，已经有人做了类似的工作，比如ICLR 2020的一篇paper便是针对CIFAR-10中32×32
32
×
32
的图片抽 2×2
2
×
2
的像素块)
16*16    16*16    16*16    16*16    16*16    16*16    16*16    16*16    16*16    16*16    16*16    16*16    16*16    16*16
16*16                                                    
16*16                                                    
16*16                                                    
...                                                    
而对于图片而言，还得考虑RGB channel这个因素，故每个 patch 的 维度便是 [16, 16, 3]，但标准Transformer的输入是一个一维的token序列，所以需要把这个三维的维度通过线性映射linear projection成一维的维度，从而使得每个 patch 最终的输出维度是：16×16×3=768
16
×
16
×
3
=
768

这样图片就由原来的 [224, 224, 3] 的大小变成了 [14, 14, 768] 的大小，相当于之前图片横竖都是224个像素点的，现在横竖只是14个像素点了，而每个像素点的维度(相当于每个token的序列长度)为768

之后经过 Flatten 就得到 [196,768]
[
196
,
768
]
，接着再经过一个维度为768×768
768
×
768
的Linear projection (本质上就是一个全连接层，用大写 E
�
表示，这个768的维度即embedding_dim可以变，简写为D
�
，比如原始的transformer设置的维度为512)，故最终的维度还是为196×768
196
×
768

 
至于在代码实现中，可通过一个卷积层来实现，卷积核大小为16，步长为16，输入维度是3，通过对整个图片进行卷积操作：[224, 224, 3] -> [14, 14, 768]，然后把H以及W两个维度展平即可[14, 14, 768] -> [196, 768]

Conv2d(in_c, embed_dim, kernel_size=patch_size, stride=patch_size)
接下来，为了做最后的分类，故在所有tokens的前面加一个可以通过学习得到的 [class] token作为这些patchs的全局输出，相当于BERT中的分类字符CLS (这里的加是concat拼接)，得益于self-attention机制，所有token两两之间都会做交互，故这个[class] token也会有与其他所有token交互的信息
且为了保持维度一致，[class] token的维度为 [1, 768] ，通过Concat操作，[196, 768]  与 [1, 768] 拼接得到 [197, 768]

由于self-attention本身没有考虑输入的位置信息，无法对序列建模。而图片切成的patches也是有顺序的，打乱之后就不是原来的图片了，故随后和transformer一样，就是对于这些 token 添加位置信息，也就是 position embedding，VIT的做法是在每个token前面加上位置向量position embedding(这里的加是直接向量相加即sum，不是concat)，这里和 transformer 一致，都是可训练的参数，因为要加到所有 token 上，所以维度也是 [197, 768]
Transformer Encoder

维度为[197, 768]的embedded patches进来后，先做一次Layer Norm
然后再做Multi-head Attention，通过乘以三个不同的Q K V矩阵得到三个不同的Q K V向量，且ViT_base_patch16设计的是12个头，故每个头的维度为：[197, 768/12] = [197, 64]，最后把12个头拼接起来，会再次得到[197, 768]的维度
接着再做Norm
最后是MLP，维度上先放大4倍到[197, 3072]，之后又缩小回去恢复到[197, 768]的维度
MLP Head（最终用于分类的层结构）
MLP里面，是用tanh作为一个非线性的激活函数，去根据[class] token做分类的预测
再通过小绿豆根据ViT的源码画的这个图总结一下

 此外，上述这三个阶段的过程最终可以用如下4个公式表达

对于公式(1)，xNp
�
�
�
表示图片patch，总共N
�
个patch，xNpE
�
�
�
�
则表示patch embedding，xclass
�
�
�
�
�
�
表示class embedding，因为需要用它做最后的输出/分类
对于公式(2)，z′l
�
�
′
则表示多头注意力的结果(zl−1
�
�
−
1
先Norm再multi-head attention，得到的结果再与zl−1
�
�
−
1
做残差连接)
对于公式(3)，zl
�
�
表示最终整个transformer decoder的输出(z′l
�
�
′
先做Norm再做MLP，得到的结果再与z′l
�
�
′
做残差连接)

4.1.2 ViT与CNN在先验知识上的对比
值得一提的是，由于ViT不像CNN那样对图像有比较多的先验知识，即没有用太多的归纳偏置

具体来说，CNN的局部性locality (以滑动窗口的形式一点一点在图片上进行卷积，故会假设图片上相邻的区域会有近似的特征)，和平移等变性translation equivariance「无论先做平移还是先做卷积，最后的 结果都是一样的，类似f(g(x))=g(f(x))
�
(
�
(
�
)
)
=
�
(
�
(
�
)
)
，毕竟卷积核就像一个模板一样，输入一致的情况下，不论图片移动到哪里，最后的输出都是一样的」贯穿整个CNN模型的始终
而对于ViT而言，也就在最后的MLP用到了局部且平移等变性，以及针对每个图片的patch加了位置编码，除这两点之外，ViT没有再专门针对CV问题做任何额外的处理，说白了，就是干：直接拿transformer干CV
所以在中小型的数据集上训练的结果不如CNN也是可以理解的。既如此，transformer的全局建模能力比较强，而CNN又不用那么多的训练数据，那是否可以把这个模型的优势给结合起来呢？比如做一个混合网络，前头是CNN，后头是transformer呢，答案是：也是可以的! 但这是不是就类似上文介绍过的DETR呢？读者可以继续深入思考下。

4.2 MAE

 // 待更

第五部分 Swin Transformer
5.1 Swin Transformer：多尺度的ViT
swin transformer作为多尺度的ViT更加适合处理视觉问题，且证明transformer不但能在ViT所证明的分类任务上取得很好的效果，在检测、分割上也能取得很好的效果，而在结构上，swin transformer主要做了以下两点改进

获取图像多尺寸的特征
对于ViT而言，经过12层每一层的transformer都是16✖️16的patch块(相当于16倍下采样率)，虽然通过transformer全局的自注意力操作可以达到全局的建模能力，但它对多尺寸特征的把握则相对弱些，而这个多尺寸的特征有多重要呢，比如对于目标检测而言，用的比较广的一个方法叫FPN，这个方法用的一个分层式的CNN，而每一个卷积层因为不同的感受野则会获取到不同尺寸的特征
而swin transformer就是为解决ViT只有单一尺寸的特征而来的
降低序列长度是图像处理中一个很关键的问题，虽然ViT把整张图片打成了16✖️16的patch，但但图像比较大的时候，计算的复杂度还是比较高
而Swin Transformer使用窗口Window的形式将16✖️16的特征图划分成了多个不相交的区域(比如16个4✖️4的，4个8✖️8的)，并且只在每个小窗口(4✖️4或8✖️8)内进行多头注意力计算，大大减少计算量
之所以能这样操作的依据在于借鉴了CNN中locality先验知识(CNN是以滑动窗口的形式一点一点地在图片上进行卷积，原因在于图片上相邻的区域会有相邻的特征)，即同一个物体的不同部位在一个小窗口的范围内是临近着的，从而在小窗口内算自注意力够用
且swin transformer使用patch merging，可以把相邻的四个小的patch合成一个大的patch(即patch merging)，提高了感受野，这样就能获取多尺度的特征(类似CNN中的池化效果)，这些特征通过FPN结构就可以做检测，通过UNet结构就可以做分割了
5.1.1 如何让4个互相独立不重叠的窗口彼此交互做自注意力操作：经过一系列shift操作
接下来，我们看下swin transformer中移动窗口的设计，图中灰色的小patch就是4✖️4的大小，然后最左侧的4个红色小窗口中均默认有7✖️7=49个小patch(当然，示意图中只展示了16个小patch)，如果做接下来几个操作

把最左侧的整个大窗口layer 1向右下角整体移动两个小patch，并把移动之后的大窗口的最右侧的宽为2个小patch、高为6个小patch的部分平移到大窗口之外的左侧
且同时把移动之后的大窗口最底部宽为8个小patch、高为2个小patch部分的一分为二：
→
→
  一方面 其左侧部分(宽为6个小patch 高为2个小patch)整体平移到大窗口之外的最上方
→
→
  二方面 最后遗留下来的右下角的2个小patch移动到大窗口之外的最左上角

则成为图中右侧所示的大窗口layer 1+1，从而使得之前互相独立不重叠的4个小窗口在经过这一系列shift操作之后，彼此之间可以进行互动做自注意力的计算了

貌似还是有点抽象是不？没事，我画个图 就一目了然了

如下所示，在右侧加粗的4个新的小窗口内部，每个小窗口都有其他小窗口的信息了(每个小窗口都由之前的单一颜色的patch组成，变成了由4种不同的颜色patch组成，相当于具备了全局的注意力，够直观吧?!)

5.1.2 Swin Transformer模型总览图
以下是整个swin transformer模型的总览图

 从左至右走一遍整个过程则是

stage 1
patch partition
对于一张原始图片224✖️224✖️3，打成4✖️4的patch(则每个patch的维度是4✖️4✖️3 = 48，其中3是图片的RGB通道)，从而会存在2242/42=56×56
224
2
/
4
2
=
56
×
56
个patch，相当于整个图片由[224,224,3]的维度变成了[56,56,48]的维度
linear embedding
为了变成transformer能接受的值，[56,56,48]的维度变成[56,56,96]，最前面的两个维度拉直之后，则维度变成了[3136,96]，很明显在ViT里在这一步对应的维度是[196,768]，故这个3136的维度太大了，咋办呢？
swin transformer block：基于7✖️7个小patch的小窗口计算自注意力
好在swin transformer引入了基于窗口的自注意力机制，而每个窗口默认只有七七四十九个patch，所以序列长度就只有49了，也就解决了计算复杂度的问题
stage 2
patch mergeing(很像lower level任务中很常用的一个上采样方式：pixel shuffle)
patch mergeing的作用在于把临近的小patch合并成一个大patch，比如针对下图中维度为[H,W,C]的张量
由于下采样两倍，所以选点的时候，是每隔一个点选一个，最终整个大张量变成了4个小张量(每隔张量大小为[H/2,W/2])
之后把这4个张量的在C的维度上拼接起来，拼接之后的张量的维度则就变成了[H/2,W/2,4C]
接着，在C这个维度上通过一个1✖️1的卷积操作把张量的维度降了下来(这个1✖️1的卷积操作类似于linear的作用)，最终变成了[H/2,W/2,2C]

所以对应到模型总览图上则是：[56,56,96]的维度变成了[28,28,192]的维度

stage 3
维度上从[28,28,192]变成[14,14,384]
stage 4
维度上从[14,14,384]变成[7,7,768] 
整个前向传播过程走完了之后，可能有读者问，swin transformer如何做分类呢？它为了和CNN保持一致，没有像ViT在输入序列上加一个用于最后分类的CLS token，而是在得到最后的特征图之后，用了一个golbal average polling(即全局池化)的操作，直接把[7,7,768]中的7✖️1取平均并拉直成1，使得最终的维度变成[1,768]

//待更新..

参考文献与推荐阅读
变分自编码器（一）：原来是这么一回事
VAE原始论文
关于VAE的几篇文章：一文理解变分自编码器（VAE）、机器学习方法—优雅的模型（一）：变分自编码器（VAE）、
苏剑林关于扩散模型的几篇文章：（一）：DDPM = 拆楼 + 建楼、（二）：DDPM = 自回归式VAE
怎么理解今年 CV 比较火的扩散模型（DDPM）？
知乎上关于扩散模型的几篇文章：全网最简单的扩散模型DDPM教程、Diffusion扩散模型大白话讲解、扩散生成模型: 唯美联姻物理概念与机器学习 

Understanding Diffusion Models: A Unified Perspective(写于2022年8月，此文写的非常细致，另，这是其PDF版本)
扩散模型是如何工作的：从0开始的数学原理
What are Diffusion Models?，写于2021年7月
Introduction to Diffusion Models for Machine Learning
关于扩散模型的几篇论文
CVPR 2022 Tutorial: Denoising Diffusion-based Generative Modeling: Foundations and Applications
Diffusion Models Beat GANs on Image Synthesis
关于DDPM的几篇文章(如果相关文章与本文有冲突，建议以本文为准，因为有些文章有笔误或错误)：DDPM概率扩散模型（原理+代码)、Denoising Diffusion Probabilistic Models (DDPM)、从VAE到DDPM、扩散模型原理解析
Stable Diffusion一周年：扩散模型编年简史
DALL·E2（生成模型串讲，从GANs、VE/VAE/VQ-VAE/DALL·E到扩散模型DDPM/ADM）
AI生成艺术的底层原理：非平衡物理的扩散模型
Vision Transformer 超详细解读 (原理分析+代码解读) (一)
DETR论文的解读之一
End-to-End Object Detection with Transformers，DETR原始论文
Vision Transformer 超详细解读 (原理分析+代码解读) (二)
AN IMAGE IS WORTH 16X16 WORDS: TRANSFORMERS FOR IMAGE RECOGNITION AT SCALE
ViT原始论文，挑战CNN的在CV领域的统治地位
ViT论文逐段精读，这是针对该视频解读的笔记之一(神洛华)、针对该视频解读的笔记之二(MT_Joy)
Vision Transformer 论文 + 详解，Vision Transformer详解by 小绿豆
Masked Autoencoders Are Scalable Vision Learners
MAE原始论文
MAE 论文逐段精读，此文则为对该视频解读所做的笔记
Swin Transformer: Hierarchical Vision Transformer using Shifted Windows
Swin transformer原始论文
Swin Transformer论文精读，此文是针对该视频解读所做的笔记(神洛华)、Swin-Transformer网络结构详解(by小绿豆)
图解Swin Transformer、如何看待微软亚洲研究院的Swin Transformer？
李宏毅关于生成式AI模型diffusion model/stable diffusion概念讲解