[論文レビュー] An ergodic problem for Mean Field Games: qualitative properties and numerical simulations

Simone Cacace, Fabio Camilli+2 more|arXiv (Cornell University)|Jan 26, 2018

Advanced Mathematical Modeling in Engineering参考文献 4被引用数 6

ひとこと要約

本稿は、小ノイズ下での平均場ゲーム（MFG）系の漸近的極限を調査し、べき乗型および対数型非線形性における有効ハミルトニアン $\bar{H}$ およびドリフト $\bar{b}$ の定性的性質を分析する。理論的に、対数型の場合にはMFG構造が保持される（$\nabla_P\bar{H} = \bar{b}$）が、べき乗型の場合には失われる（$\nabla_P\bar{H} \neq \bar{b}$）ことが示され、有限差分スキームおよび乖離度 $\mathcal{R}$ の収束解析により数値的に検証されている。

ABSTRACT

This paper is devoted to some qualitative descriptions and some numerical results for ergodic Mean Field Games systems which arise, e.g., in the homogenization with a small noise limit. We shall consider either power type potentials or logarithmic type ones. In both cases, we shall establish some qualitative properties of the effective Hamiltonian $\bar H$ and of the effective drift $\bar b$. In particular we shall provide two cases where the effective system keeps/looses the Mean Field Games structure, namely where $ abla_P \bar H(P,α)$ coincides or not with $\bar b(P, α)$. On the other hand, we shall provide some numerical tests validating the aforementioned qualitative properties of $\bar H$ and $\bar b$. In particular, we provide a numerical estimate of the discrepancy $ abla_P \bar H(P,α)-\bar b(P, α)$.

研究の動機と目的

• 平均場ゲームの小ノイズ極限における有効ハミルトニアン $\bar{H}$ およびドリフト $\bar{b}$ の定性的な振る舞いを研究すること。
• 漸近的有効系がMFG構造を保持する条件、すなわち $\nabla_P\bar{H} = \bar{b}$ が成り立つかどうかを特定すること。
• 周期的領域における有限差分スキームを用いて、特に乖離度 $\nabla_P\bar{H} - \bar{b}$ の数値的妥当性を検証すること。
• 異なる非線形性について、$|P| \to \infty$ における $\bar{H}$ および $\bar{b}$ の漸近的振る舞いを分析すること。
• グリッドの細分化に伴う、関係式 $\partial_P\bar{H} = \bar{b} - \alpha \int V_m \tilde{m} m\,dy$ における数値誤差の収束次数を確認すること。

提案手法

• 1周期的かつリプシッツ連続な関数 $v$ を用いて、$V(y,m) = v(y) + m^q$（べき乗型）および $V(y,m) = v(y) + \log m$（対数型）とする、トーラス $\mathbb{T}^n$ 上の漸近的MFGセル問題を解く。
• 漸近的系（3）の解 $(u,m)$ を用いて、$\bar{b}(P,\alpha) = \int_{\mathbb{T}^n} (\nabla u + P) m\,dy$ を用いて $\bar{H}(P,\alpha)$ および $\bar{b}(P,\alpha)$ を計算する。
• 線形化セル問題（6）の $\tilde{m}_i$ を数値的に解き、乖離度 $\mathcal{R}(P,\alpha) = \left| \alpha \int V_m \tilde{m} m\,dy \right|$ を計算する。これは $\nabla_P\bar{H}$ と $\bar{b}$ の差を測る指標である。
• 一様グリッド上に2次精度の有限差分スキームを適用し、$\mathbb{T}^1$ 上に $N$ 個のノードを配置し、$P \in [-10,10]$、$\alpha \in [1,20]$ の $51 \times 51$ パラメータグリッドを用いる。
• $\mathcal{E}_\infty = \left\| \partial_P\bar{H} - \bar{b} + \alpha \int V_m \tilde{m} m\,dy \right\|_\infty$ の誤差を推定し、グリッドの細分化に伴う収束率を分析する。
• 恒等式における微分項の妥当性を検証するため、直接的な残差計算と $\|m_{(P,\alpha)}\|_{L^{q+1}}^{q+1}$ の有限差分近似の両方を用いる。

実験結果

リサーチクエスチョン

• RQ1 ポテンシャル $V(y,m)$ がどのような条件下で、有効系が平均場ゲーム構造（$\nabla_P\bar{H} = \bar{b}$）を保持するか。
• RQ2 べき乗型非線形性において、$|P| \to \infty$ のとき、有効ハミルトニアン $\bar{H}(P,\alpha)$ はどのように漸近的に振る舞うか。
• RQ3 恒等式 $\partial_P\bar{H} = \bar{b} - \alpha \int V_m \tilde{m} m\,dy$ における誤差 $\mathcal{E}_\infty$ の数値的収束次数は何か。
• RQ4 対数型ポテンシャル $V(y,m) = v(y) + \log m$ は、$\bar{H}(P,\alpha)$ に分離可能な構造をもたらすか。また、$\bar{b}$ は $\alpha$ に依存しないか。
• RQ5 パラメータ空間全体にわたり、残差 $\mathcal{R}(P,\alpha)$ はどのように振る舞い、対数型の場合にゼロに収束するか。

主な発見

• べき乗型ポテンシャル $V(y,m) = v(y) + m^q$ の場合、MFG構造は失われる：数値シミュレーションにより非ゼロの残差 $\mathcal{R}(P,\alpha)$ が確認され、$\nabla_P\bar{H}(P,\alpha) \neq \bar{b}(P,\alpha)$ であることが示された。
• 対数型の場合 $V(y,m) = v(y) + \log m$ では、$V_m = 1/m$ かつ $\tilde{m}$ が平均ゼロであるため、残差がゼロ（$\mathcal{R}(P,\alpha) = 0$）となり、$\nabla_P\bar{H} = \bar{b}$ が確認された。
• $q=1$ および $q=2$ の両方において、数値誤差 $\mathcal{E}_\infty$ はグリッドの細分化に伴い1次収束を示し、有限差分スキームの1次収束性が裏付けられた。
• $|P| \to \infty$ のとき、$\bar{H}(P,\alpha)/P^2 \to 1/2$ かつ $|\bar{b}(P,\alpha) - P|/P \to 0$ となることが確認され、命題2.1-(iv)が予測する漸近的挙動と整合的であった。
• 対数型の場合、$\bar{H}(P,\alpha)$ は $P$ に関して凸であり、$\bar{b}(P,\alpha)$ は $\alpha$ に依存しないことが確認され、理論的分離変数構造が裏付けられた。
• 対数型の場合、$\mathcal{E}_\infty$ の $L^\infty$ ノルムは $\alpha$ にかかわらず一定であり、1次収束を示した。これにより、理論的恒等式が妥当であることが検証された。