[论文解读] Mean field games equations with quadratic Hamiltonian: a specific approach

Olivier Guéant|arXiv (Cornell University)|Jun 16, 2011

Stochastic processes and financial applications参考文献 6被引用 70

一句话总结

本文提出了一种新颖的构造性单调格式，用于求解具有二次哈密顿量的平均场博弈（MFG）方程，通过变量变换将原系统转化为一对耦合的热方程。该方法在变换变量 φ 和 ψ 上采用迭代更新，通过单调性保证收敛性，并利用隐式有限差分格式结合牛顿求解器实现有效数值求解，提供了一种与以往基于优化或直接 u/m 格式不同的新方法。

ABSTRACT

Mean field games models describing the limit of a large class of stochastic differential games, as the number of players goes to $+\\infty$, have been introduced by J.-M. Lasry and P.-L. Lions. We use a change of variables to transform the mean field games (MFG) equations into a system of simpler coupled partial differential equations, in the case of a quadratic Hamiltonian. This system is then used to exhibit a monotonic scheme to build solutions of the MFG equations. Effective numerical methods based on this constructive scheme are presented and numerical experiments are carried out.

研究动机与目标

• 开发一种用于求解具有二次哈密顿量的平均场博弈方程的构造性且单调的数值格式，此类方程在经济与博弈论建模中十分常见。
• 通过引入一种新颖的变换，将 MFG 系统解耦为两个带有共享源项的热方程类方程，以克服现有数值方法的局限性。
• 通过变换变量 φ 和 ψ 中的单调序列，确保求解过程的收敛性与稳定性。
• 提供一种计算高效且鲁棒的数值算法，避免在迭代过程中出现质量守恒违反，优于以往的方案。
• 通过数值实验验证方法，展示其收敛速率、计算复杂度以及在不同参数（如 σ）下的行为表现。

提出的方法

• 该方法应用变量变换：φ = exp(u/σ²) 和 ψ = m exp(-u/σ²)，将原始 MFG 系统转化为两个耦合的抛物型方程，其形式类似于带有源项的热方程。
• 该格式构建两组序列：φ^{n+1/2} 和 ψ^{n+1}，通过隐式有限差分离散化迭代求解变换后的方程。
• 在每一步中，φ^{n+1/2} 沿时间向前求解，终端条件为 φ(T) = exp(u_T/σ²)，而 ψ^{n+1} 沿时间向后求解，初始条件为 ψ(0) = m₀ / φ^{n+1/2}(0)。
• 迭代过程确保了两组序列对真实 φ 和 ψ 的单调收敛，其依据是在 f 的标准假设下系统的单调性特性。
• 离散格式为全隐式，每一步均使用牛顿法求解，从而在数值实现中保证鲁棒性与稳定性。
• 该算法在中间步骤不强制质量守恒，从而实现更快的收敛速度与更优的数值行为。

实验结果

研究问题

• RQ1 能否通过变量变换构造一种用于具有二次哈密顿量的平均场博弈的单调迭代格式，以简化系统？
• RQ2 将系统变换为 φ 和 ψ 方程后，相较于直接求解 u 和 m，其在数值稳定性与收敛性方面有何改善？
• RQ3 所提格式在时间与空间离散化下的收敛速率与计算复杂度如何？
• RQ4 为何该格式在迭代过程中无需强制质量守恒仍能保持单调性，这种特性对解的精度有何影响？
• RQ5 参数 σ 如何影响算法的收敛速度与计算成本？

主要发现

• 所提格式在时间和空间上均实现一阶收敛，误差范数随 Δt 和 Δx 线性减小，表明收敛性成立。
• 计算时间与 O((Δx Δt)^{-1}) 成比例，当 Δx 与 Δt 成比例时，表明具有二次复杂度，与该格式的理论预期一致。
• 近似 m^n+1 序列的总质量在早期迭代中减少，但随着解的收敛而趋于稳定，反映出迭代过程的非守恒特性。
• σ 值越小，计算成本显著增加，原因在于 φ 和 ψ 对 σ² 具有指数敏感性，导致需要更多迭代次数及每步更多牛顿迭代。
• 在实际应用中，该格式收敛迅速，通常在标准参数下 5–6 次迭代内即可满足停止准则（||φ^{n+1/2}ψ^{n+1} - φ^{n-1/2}ψ^n||_∞ < 10^{-7}）。
• 数值实验表明，解的动力学行为符合预期：参与者最初集中在稳态均衡附近，但当终端激励发生变化时，其分布趋于在时间终点附近扩散。