Dyson 方程的最优控制问题是如何求解的？

Dyson 方程描述平均场谱演化

Dyson 方程描述大随机矩阵谱的平均场演化。当矩阵规模 N 趋于无穷时，Dyson Brownian 运动特征值的经验测度收敛到该方程的解；Dyson Brownian 运动是一个带排斥对数库仑相互作用的随机粒子系统。按论文记号，经验测度为 m_{N,t} = (1/N) 乘以特征值处 Dirac 质量之和，其 N 趋于无穷的极限满足带 Hilbert 变换项 H[m] 的 Dyson 方程。

这里研究的控制问题在该演化中加入一个可测速度 alpha。给定把时间与空间映射到实轴的控制 alpha，受控测度按一个连续性方程演化，其通量为 m 乘以 (H[m] + alpha)。需要最小化的收益由二次型运行成本、F[m_s] 运行成本与终端成本 U_T(m_T) 组成，其中二次型运行成本是 alpha 平方的一半对测度积分。优化在满足受控方程、测度路径在 2-Wasserstein 空间中连续且能量有限的所有对 (alpha, m) 上进行。

正则与奇异终端数据需分别处理

论文区分两类终端条件。第一类中，终端成本 U_T 是有限二阶矩概率测度空间上的连续函数，作者把这一情形与软大偏差联系起来。第二类中，U_T 是单个目标测度 nu 的指示函数，在 nu 处取零、在其他处取无穷；这一情形与转移概率的大偏差渐近相联系。

这一区分很重要，因为目标问题固定终端测度，要求该终端测度具有有限的自由熵。作者指出，假设终端测度具有有限熵，结果证明这是解存在的必要条件。因此两类情形不可互换：一类允许终端成本有分布，另一类把端点钉在特定测度上，并对该测度施加额外的可容许条件。

可容许控制被构造且最优控制存在

一个核心贡献是关于一般控制 alpha 下受控 Dyson 方程解存在性的一组新结果。作者指出，一般控制下受控方程的适定性此前基本是开放问题，因为只研究过光滑向量场的情形。他们构造了可容许控制，并证明最优控制存在。

支撑这一点的紧性结果是命题 3.1：对于可容许集 A(mu, t) 中的对 (alpha, m)，存在仅依赖于能量、时间跨度与初始测度二阶矩的常数 C，使得 m_s 的二阶矩保持有界，自由熵 E(m_s) 被 C 乘以 (1 + (s - t)^(-1/2)) 所界。同一命题给出一个能量恒等式，把 m 的三次方在时间区间上的积分用两倍控制能量加上与熵下降有关的项来界定。换言之，有界的控制能量换来了测度路径的紧性与正则性。

Hilbert 变换项使 HJB 方程奇异

在没有 Hilbert 变换且 F 恒为零时，该最小化问题就是 Benamou 和 Brenier 在质量输运中研究的问题，其 Wasserstein 空间上的 HJB 方程已在早期工作中得到分析。作者指出，当 H 出现时，问题更复杂也更奇异，证明 HJB 方程的适定性以及值函数确实求解该方程，需要在 Wasserstein 空间上的 Hamilton-Jacobi 方程方面发展新思想，也需要带漂移的 Dyson 方程的新结果。

论文导出的 HJB 方程在 Wasserstein 空间与时间区间的乘积上按时间反向。它涉及值函数的水平导数（即 Wasserstein 导数）D_mu U、其平方范数对测度的积分，以及把 D_mu U 与 Hilbert 变换 H[mu] 配对的项。终端条件为 U(., T) = U_T。值函数 U(mu, t) 定义为收益在可容许控制上的下确界，论文证明它是该方程的唯一粘性解。

粘性解在测度空间中定义

由于 HJB 方程定义在无穷维测度空间上，论文采用测度空间中的粘性解理论。作者把提出合适的粘性解概念并建立其适定性——尤其是在 U_T 为目标测度指示函数的奇异终端情形——描述为论文的两项主要贡献。

分析使用具有有限二阶矩概率测度的 2-Wasserstein 空间、该空间上函数的水平导数，以及一类 C^{1,1} 型正则性。一个关键技术工具是用正则化映射逼近 Wasserstein 距离平方与最优输运成本，这使作者能够处理 Wasserstein 距离平方未必处处可微的问题。自由熵 E(m) = -(1/2) 对 log|x - y| dm(x) dm(y) 的二重积分也起核心作用；它关于测度的弱收敛是下半连续的，并沿 Dyson 流递减。

值函数是唯一的粘性解

对于局部一致连续的 F 与 U_T，定理 6.2 指出 U 是 HJB 方程满足 U(mu, T) = U_T(mu)（对所有 mu）的唯一连续粘性解。对于固定终端测度 nu 的目标问题，定理 6.3 在附加条件下给出刻画：候选 V 必须局部一致连续，被 C^{-1} W_2^2(., nu) - C 从下方界住，满足 V(mu, nu, t) 至少为 -C(T - t)，当 t 趋于 T 时沿对角线 V(nu, nu, t) 趋于零，并且当 t 趋于 T 时在远离 nu 处局部一致地趋于无穷。在这些条件下，V 就是目标问题的值函数。

目标情形的证明使用比较原理与动态规划关系。作者还通过惩罚论证构造了一族 U_epsilon，其终端数据为 (1/(2 epsilon)) W_2^2(., nu)；命题 6.6 表明当 epsilon 趋于零时，该族收敛到目标问题的值函数。奇异终端条件正是以正则终端条件的极限方式达到的。

结果覆盖二次成本与有限二阶矩

论文把注意力限制在具有有限二阶矩的概率测度空间中的初始数据，并表明对于有限成本的控制，动力学仍留在该空间中。作者指出，他们关于 Dyson 方程控制及所得 HJB 方程的结果可推广到比二次型能量与分离型运行成本更一般的情形，即对 alpha 一致凸的成本 L(alpha, m)，并把细节留给感兴趣的读者。该推广只是被陈述，在现有摘录中并未证明。

动机是研究 Dyson Brownian 运动的大偏差。作者表示，在后续论文中他们将给出大随机矩阵谱大偏差的完整直接证明，并研究 Dyson Brownian 运动转移概率的渐近行为。因此本文是控制论与 HJB 基础，大偏差结论本身被推迟到那篇配套工作中，并未在本文建立。

奇异终端情形是最尖锐的检验

本文中最具决策相关性的区分是正则终端情形与奇异终端情形。在正则情形中，只要 F 与 U_T 局部一致连续，连续粘性解的唯一性即可得到，证明被描述为早期结果的推论。在奇异情形中，U_T 是目标测度 nu 的指示函数，证明更复杂，因为它依赖比较原理，并且终端测度必须具有有限的自由熵。

因此，对于判断本文实际贡献的读者，优先关注点是奇异终端问题：正是在这里，新的粘性解概念及其适定性最为需要，nu 的有限熵条件成为真正的可容许阈值，而终端数据为 (1/(2 epsilon)) W_2^2(., nu) 的惩罚族 U_epsilon 提供了从正则数据到奇异数据的桥梁。正则情形是基线；奇异情形才是改变 HJB 方程必须如何提出与求解的那一类。

本文引用的文献