Schrödinger Bridge:最优传输、熵正则与生成模型的统一数学框架
一、从最优传输到 Schrödinger Bridge
1.1 最优传输问题回顾
Monge 最优传输问题(1781):
给定两个概率分布
其中:
是代价函数(cost function) 表示 将 映射到 (push-forward) - 典型选择:
物理直觉:
这相当于将一堆沙土(
Brenier 定理(当
如果
其中
1.2 Schrödinger 问题的引入
需区分三类不同的 Schrödinger Bridge:
- 经典统计物理 SB(1931 原始):平衡态假设,固定稳态分布,路径熵最小化
- 动力学 SB:以路径运动动能
为优化目标 - 生成建模非平衡 SB:当前 AI 主流,用于时序/插值分布演化,本文聚焦此版本
原始 Schrödinger 问题(1931):
设
其中
动力学版本 SB(动作泛函形式):
重要澄清:式(1.3)与式(1.4)仅在单位扩散系数、平衡稳态假设下弱等价,并非无条件等价的原始定义。式(1.4)是动力学版本的 SB,而非 1931 年薛定谔原始定义。
物理直觉:
在所有连接给定端点分布的随机过程中,寻找”最像布朗运动”的路径。代价函数
1.3 Schrödinger Bridge 与最优传输的联系
定理(Schrödinger Bridge 的 OT 等价性):
当路径无噪声(即
证明思路:
无噪声时,路径是确定性的:
动作泛函:
由积分型柯西-施瓦茨不等式:
因此:
等号成立当且仅当
因此,无噪声 SB 的最优路径是直线,这正是 Rectified Flow 的设计。
二、Schrödinger Bridge 的数学理论
2.1 路径空间与 Radon-Nikodym 导数
路径空间定义:
设
对于维纳过程
** Schrödinger Bridge 是
其中
2.2 Fokker-Planck 方程约束
定理(SB 的必要性条件):
设
其中
对应的 Fokker-Planck 方程(仅适用于单位扩散系数
边际分布
其中
2.3 互熵最小化视角
定义(路径 KL 散度):
设
SB 的变分形式(最小化路径 KL 而非互熵):
通过拉格朗日乘子法求解:
引入边际约束的拉格朗日乘子
通过变分法求解,得到最优路径的分解形式(式 2.2)。
三、Simulation-Free Schrödinger Bridges(SFSB)
3.1 核心问题:为什么需要 Simulation-Free?
原始 SB 的计算瓶颈:
直接求解 Schrödinger Bridge 需要:
- 在路径空间上优化分布
- 通过 Fokker-Planck 方程反向计算 score
- 这需要大量的 SDE 模拟(simulation),算力成本极高
“Simulation-Free” 的真正核心内涵:
SFSB 舍弃全局路径遍历,仅利用两端静态分布和时序插值完成训练,彻底摆脱传统 SB 依赖 SDE 轨迹模拟的算力瓶颈。不通过显式模拟 SDE 路径来训练,而是通过端点驱动的条件分布进行训练。
3.2 SFSB 的数学框架
关键观察:
在独立采样的假设下(
定理(SFSB 的损失函数):
设
其中:
3.3 与 Score Matching 的等价性
SFSB 等价去噪分数匹配需满足强约束:
当以下全部三个前置条件同时满足时,SFSB 损失(式 3.1)等价于条件去噪 score matching:
- 端点独立采样:
, 相互独立 - 高斯插值核:
- 固定时序调度:
预先确定(非学习参数)
满足约束时的推导:
条件分布
Score 函数:
因此训练目标变为预测
脱离约束条件时,SFSB 与 Score Matching 仅为形式相似,数学期望层面并不等价。
3.4 与 Flow Matching 的联系
当
SFSB 退化为条件速度场回归:
这正是 Flow Matching 的条件流匹配损失!
3.5 统一视角下的模型谱系
Schrödinger Bridge(统一框架)
├── 无噪声路径(σ=0)
│ └── Flow Matching / Rectified Flow
└── 有噪声路径(σ≠0)
├── Score Matching(条件版本)
└── DDPM(特殊噪声调度)
四、Entropy-Regularized Schrödinger Bridge
4.1 熵正则化的引入
标准 SB 的问题:
当
熵正则化 SB:
第二项是路径熵的正则化,鼓励路径更加”均匀”分布。
4.2 Sinkhorn 算法的联系
通过路径空间分解:
设
其中
与 Schrödinger Bridge 的联系:
熵正则化 OT 的最优传输映射等价于某个熵正则化 SB 的终点映射。
4.3 训练中的应用
Practical Entropic SB 训练:
- 初始化:通过 Sinkhorn 计算初始配对
- 路径采样:沿配对
采样 对 - 路径优化:在采样路径上训练 score/velocity 模型
- 迭代改善:重复上述步骤直到收敛
五、Schrödinger Bridge 与 Diffusions 的深层联系
5.1 互为对偶的关系
定理(SB-Diffusion 对偶性):
设
其中
直观理解:
- Diffusion 模型:从数据分布出发,通过加噪过程趋近于噪声分布
- Schrödinger Bridge:从噪声分布出发,通过”最像布朗运动”的路径趋近于数据分布
两者是同一个数学对象的不同方向描述。
5.2 路径测度的表示
克喇莫尔表示(仅适用于高斯端点分布):
设
即在给定端点分布的条件下,布朗运动的条件期望。
注意:此显式表示仅严格成立于纯高斯分布场景,对于复杂多峰数据无通用显式表达式,需通过数值方法近似。
5.3 与 Score SDE 的等价性
定理(SB = Score SDE):
对于任意的 Schrödinger Bridge
构造:
给定
其中
定义 Score SDE:
可以验证,这个 SDE 的边际分布与
六、训练过程与算法
6.1 SFSB 的训练算法
Algorithm: Simulation-Free Schrödinger Bridge Training
输入:
- 数据分布
- 噪声分布
- 调度函数
- 网络
(score network)或 (velocity network)
for iteration
- 采样批次
,其中 , (独立) - 采样时间步
- 采样噪声
- 计算插值点:
- 计算目标 score:
- 计算损失:
- 反向传播并更新
输出:训练好的网络
6.2 带 OT 配对的改进版本
Algorithm: OT-Regularized SFSB
初始化:
- 通过 Sinkhorn 算法计算
和 之间的最优配对 - 令
,其中 是 对应的映射
主循环(与上面相同,但配对固定)
收敛性:
当配对是最优传输映射时,直线插值路径就是 SB 的最优路径。
6.3 多步 Rectification(Re-flow 的 SB 解释)
Re-flow 是 SB 最优路径的工程数值近似解法:
Re-flow(重边缘化整流匹配)并非等价于 SB,而是 SB 最优路径的工程数值近似实现。二者关系:
- SB:理论最优路径定义(路径熵最小化)
- Re-flow:通过迭代数值近似求解 SB 的工程落地方法
具体而言:
- 初始化:假设独立配对
- 迭代:通过当前配对训练
,生成新轨迹,用轨迹终点更新配对 - 收敛:当配对趋近于 SB 最优传输映射时,路径接近 SB 最优
七、潜在问题与解决方案
7.1 路径爆炸问题
问题:当 OT 配对不是双射时,路径可能交叉或”爆炸”。
数学表现:
- 学习到的映射
不是光滑的双射 - 某些区域密度过高或过低
- 轨迹曲率在某些点变得非常大
解决方案:
-
Unbalanced OT:允许边际分布不完全匹配
-
正则化:在损失中添加曲率惩罚
7.2 训练不稳定性
问题:当
原因分析:
时,条件分布退化 在边界处可能奇异- score 的幅度
在边界处爆炸
解决方案:
- 平滑调度:使用
代替纯 - 截断梯度:对边界区域的梯度进行裁剪
- 双时间步策略:分别训练
和 的模型
7.3 维度灾难
问题:当维度
表现:
- Sinkhorn 迭代收敛慢:
- 轨迹方差随
线性增长 - 模型需要拟合
维 score 向量
解决方案:
- 子空间建模:在低维子空间中进行 OT 和插值
- 分块处理:将高维数据分块独立处理
- 对比学习:使用对比损失辅助学习
八、与现有模型的关系总结
8.1 模型对应关系图
生成模型(概率传输)
│
├── 最优传输(OT)
│ ├── Monge 问题
│ └── Brenier 定理(凸势函数)
│
└── Schrödinger Bridge
├── 标准 SB(路径熵最小化)
├── 熵正则化 SB(+OT 正则)
└── Simulation-Free SB
├── 当 σ=0:Flow Matching / Rectified Flow
├── 当 σ≠0:Score Matching / DDPM
└── 中间情况:条件插值
8.2 核心公式对照
| 概念 | 数学表达 | 物理意义 |
|---|---|---|
| Monge OT | 确定性映射,最小化传输代价 | |
| SB | 找到最像布朗运动的路径 | |
| 熵正则化 SB | OT 与最大熵的平衡 | |
| SFSB | 端点驱动的 score 训练 |
8.3 采样效率对比
| 方法 | 采样步数 | 理论依据 |
|---|---|---|
| DDPM | 50-1000 | 多步去噪链 |
| Score SDE | 10-1000 | SDE 求解 |
| Flow Matching | 1-10 | 路径直线化 |
| SB(精确) | 1(理论上) | 最优传输映射 |
| SFSB | 1-10 | 近似 SB |
九、数学推导速查
9.1 Schrödinger Bridge 核心方程
路径 KL 散度(而非互熵):
SB 变分形式:
Fokker-Planck 约束(仅适用于单位扩散系数
9.2 SFSB 损失函数
条件 score matching:
条件 velocity matching(当
9.3 OT 配对公式
Sinkhorn 势函数更新:
最优配对:
十、总结
Schrödinger Bridge 的核心价值:
-
统一理论框架:将最优传输、扩散模型、score matching、F low matching 统一在一个变分框架下
-
最优路径解释:SB 提供了”什么是最好的从噪声到数据的路径”的理论答案
-
Simulation-Free 训练:SFSB 使得在不模拟 SDE 路径的情况下训练成为可能
与生成模型实践的联系:
- SD3、Flux 等大模型都可以在 SB 框架下理解
- Re-flow 是 SB 的迭代近似求解
- 最优传输配对是 SB 与工业实践的关键桥梁
核心洞察:
Schrödinger Bridge 回答了”最优生成路径应该是什么样的”这个问题。 它告诉我们,在所有可能的随机过程中,最优的那条是”最像布朗运动”且满足端点约束的路径。这个结论将最优传输的确定性与扩散模型的随机性统一了起来。
延伸阅读:
- Schrödinger, “Über die Umkehrung der Naturgesetze” (1931)
- Föllmer, “Snell Envelope and Schrödinger Bridge” (1988)
- Léonard, “From Schrödinger to McCann’s Optimal Transport” (2014)
- Wang et al., “Schrödinger Bridge for Generative Modeling” (2024)
- De Bortoli et al., “Simulation-Free Schrödinger Bridges via Score and Flow Matching” (2024)