[论文解读] Computational methods for nonlocal mean field games with applications

Siting Liu, Matt Jacobs+3 more|arXiv (Cornell University)|Apr 25, 2020

Stochastic processes and financial applications参考文献 32被引用 8

一句话总结

本文提出了一种基于核函数的计算框架，用于非局部平均场博弈（MFG）系统，采用特征空间展开与鞍点重构方法，通过原始-对偶算法（如PDHG）实现高效求解。该方法实现线性计算复杂度，并能准确求解非局部相互作用问题，在时变与稳态MFG系统上均得到验证，收敛至已知基准解。

ABSTRACT

We introduce a novel framework to model and solve mean-field game systems with nonlocal interactions. Our approach relies on kernel-based representations of mean-field interactions and feature-space expansions in the spirit of kernel methods in machine learning. We demonstrate the flexibility of our approach by modeling various interaction scenarios between agents. Additionally, our method yields a computationally efficient saddle-point reformulation of the original problem that is amenable to state-of-the-art convex optimization methods such as the primal-dual hybrid gradient method (PDHG). We also discuss potential applications of our methods to multi-agent trajectory planning problems.

研究动机与目标

• 开发计算高效且灵活的数值方法，用于具有非局部相互作用的非局部平均场博弈（MFG）系统。
• 通过基于核函数的特征展开，解耦相互作用，克服传统方法在非局部MFG中产生的超线性计算成本。
• 通过将问题重构为鞍点形式，使原始-对偶凸优化方法（如PDHG）能够应用于非局部MFG系统。
• 通过可调制的核基函数设计，灵活建模多样化相互作用情景（如异质性或子域特定相互作用）。
• 通过长时间演化逼近稳态MFG系统，利用Turnpike性质。

提出的方法

• 使用核函数 $ K(x,y) $ 和 $ S(x,y) $ 表征非局部相互作用，并通过特征基展开降低维度。
• 基于原始-对偶框架，将MFG系统重构为鞍点问题，从而可高效使用PDHG求解。
• 通过特征空间展开解耦系统，使PDHG算法中的邻近步骤可并行化，实现线性计算成本。
• 通过显式邻近更新核展开系数，保持线性计算复杂度。
• 采用时间缩放与长时间逼近方法，通过在大 $ T $ 下嵌入时间依赖框架，求解稳态MFG系统。
• 在周期性区域中，使用三角多项式（如阶数至 $ n=2 $）逼近积分算子，实现谱精度。

实验结果

研究问题

• RQ1 能否通过基于核函数的特征展开，高效求解具有通用相互作用核的非局部平均场博弈？
• RQ2 所提出的鞍点重构方法是否允许在非局部MFG中使用PDHG实现线性计算成本？
• RQ3 该方法通过长时间动力学近似稳态MFG系统解的精度如何？
• RQ4 能否通过核基函数设计灵活建模异质性或子域特定的代理相互作用？
• RQ5 当 $ T \to \infty $ 时，时间依赖近似解向稳态MFG解的收敛行为如何？

主要发现

• 由于特征空间展开带来的解耦邻近步骤，该方法在网格点数量上实现线性计算复杂度。
• 时间依赖系统（6.4）在 $ T=10 $ 时的解弱收敛至稳态MFG解 $ \bar{\rho} $，满足 $ \|V[\rho_T(\cdot,t)] - V[\hat{\rho}]\|_{L^\infty} \leq C/T^{1/2} $。
• 该方法准确重构了解，与文献[5]中的结果相当，且由于采用一阶形式，扩散效应略小。
• 值函数 $ \phi_T(x,t) - \int \phi_T(y,t)dy $ 近似于稳态解 $ \bar{\phi} $，证实与理论预期一致。
• 通过调节核基函数，该方法成功建模复杂相互作用模式，实现异质性相互作用区域的灵活模拟。
• 对于 $ d=2 $，$ \mu=200 $，且 $ Q(x_1,x_2) = -\sin(2\pi x_2) - \sin(2\pi x_1) - \cos(4\pi x_1) $ 的情形，数值结果在 $ t=0.4 $ 时显示稳定且精确的解。