[논문 리뷰] An ergodic problem for Mean Field Games: qualitative properties and numerical simulations

Simone Cacace, Fabio Camilli외 2명|arXiv (Cornell University)|2018. 01. 26.

Advanced Mathematical Modeling in Engineering참고 문헌 4인용 수 6

한 줄 요약

이 논문은 소음이 작은 조건에서 평균장 게임(Mean Field Games, MFG) 시스템의 에르고딕 극한을 조사하며, 거듭제곱형 및 로그형 비선형성에 대해 효과적 해밀토니안 $\bar{H}$와 드리프트 $\bar{b}$의 정성적 성질을 분석한다. 논문은 로그형의 경우 ($\nabla_P\bar{H} = \bar{b}$) MFG 구조가 유지되나, 거듭제곱형의 경우 ($\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$의 수치 오차 수렴 차수를 확인한다.

제안 방법

• 거듭제곱형 $V(y,m) = v(y) + m^q$와 로그형 $V(y,m) = v(y) + \log m$의 잠재력에 대해, $v$가 1주기이고 리프시츠 조건을 만족하는 조건에서 토러스 $\mathbb{T}^n$ 상에서 에르고딕 MFG 셀 문제를 해결한다.
• 에르고딕 시스템 (3)의 해 $(u,m)$을 이용하여 $\bar{H}(P,\alpha)$와 $\bar{b}(P,\alpha)$를 계산하며, $\bar{b}(P,\alpha) = \int_{\mathbb{T}^n} (\nabla u + P) m\,dy$를 사용한다.
• 선형화된 셀 문제 (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차 유한차분 스킴을 적용하여 PDE를 이산화하며, $\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)$는 매개변수 공간 전반에서 어떻게 행동하는가? 그리고 로그형의 경우에 이 값은 0이 되는가?

주요 결과

• 거듭제곱형 잠재력 $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$에서 잔차가 0이 되며($\mathcal{R}(P,\alpha) = 0$), $V_m = 1/m$ 이고 $\tilde{m}$의 평균이 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차 수렴을 보이며 이는 이론적 항등식의 타당성을 검증한다.