[論文レビュー] A class of Hamilton-Jacobi equations with constraint: uniqueness and constructive approach

Sepideh Mirrahimi, Jean‐Michel Roquejoffre|arXiv (Cornell University)|May 22, 2015

Evolution and Genetic Dynamics参考文献 23被引用数 25

ひとこと要約

本稿は、解の上限が常にゼロに保たれる制約を伴う時間に依存するハミルトニアン・ヤコビ方程式のクラスについて、粘性解の一意性および構成的存在を確立する。動的計画法を用い、最大値に関するODEを導出することで、著者らは解が古典的であり、一意に定まることを証明し、微小な拡散を伴う集団動態における選択・変異モデルの解析を前進させる。

ABSTRACT

We discuss a class of time-dependent Hamilton-Jacobi equations, where an unknown function of time is intended to keep the maximum of the solution to the constant value 0. Our main result is that the full problem has a unique viscosity solution, which is in fact classical. The motivation is a selection-mutation model which, in the limit of small diffusion, exhibits concentration on the zero level set of the solution of the Hamilton-Jacobi equation. Uniqueness is obtained by noticing that, as a consequence of the dynamic programming principle, the solution of the Hamilton-Jacobi equation is classical. It is then possible to write an ODE for the maximum of the solution, and treat the full problem as a nonstandard Cauchy problem.

研究の動機と目的

• 解のesssupが常にゼロに保たれる最大値制約を伴うハミルトニアン・ヤコビ方程式の解の一意性に関する長年の未解決問題を解決すること。
• 先行研究で用いられた粘性近似に依存しない、制約付き系の構成的存在証明を提供すること。
• データに凹性仮定をおくと、解が粘性解ではなく古典的解であることを確立すること。
• 選択・変異の複雑な集団動的モデルにおけるハミルトニアン・ヤコビアプローチの発展を支援すること。
• 微小変異極限における、対応する選択・変異モデルの収束結果を強化し、漸近展開を可能にすること。

提案手法

• 著者らは、初期データおよび関数 $ R(x,I) $ に凹性仮定をおくと、動的計画法を用いてハミルトニアン・ヤコビ方程式の解が古典的であることを示す。
• 空間上の解の最大値に関するODEを導出し、全問題を $ u(t,x) $ とレギュレータ $ I(t) $ を含む非標準的な初期値問題として扱う。
• 制約 $ \max_x u(t,x) = 0 $ は、この条件を維持するために調整されるラグランジュ乗数 $ I(t) $ としてモデル化される。
• 粘性近似を用いないで、凹性から生じる正則性に依存する直接的な方程式の解法により、構成的アプローチが開発される。
• 一意性の証明は、2つの解の二階微分を比較し、特性曲線を用いて差を評価することで行われる。
• 2次の場合には、オイラー＝ラグランジュ方程式と行列指数関数を用いて明示解が導出され、理論的結果を確認する。

実験結果

リサーチクエスチョン

• RQ1 関数 $ R $ および $ u_0 $ に凹性仮定をおくと、最大値制約を伴うハミルトニアン・ヤコビ方程式は、解が一意に存在するか？
• RQ2 先行研究で用いられた粘性近似に依存せずに、構成的存在証明を提供できるか？
• RQ3 凹性仮定のもとで、動的計画法はどのように解の古典的正則性を導くか？
• RQ4 特に2次の場合に、解および制約関数 $ I(t) $ の長時間挙動はいかなるものか？
• RQ5 解の正則性および一意性を用いて、対応する選択・変異モデルの収束結果を強化できるか？

主な発見

• 関数 $ R $ および $ u_0 $ に凹性仮定をおくと、動的計画法のおかげで、制約付きハミルトニアン・ヤコビ方程式の解 $ u(t,x) $ は粘性解ではなく古典的解であることが保証される。
• 解のペア $ (u, I) $ の一意性は、$ u $ の最大値が明確に定義されたODEに従って進化することを示すことにより証明され、任意の2つの解は一致する。
• 制約 $ \max_x u(t,x) = 0 $ は、$ I(t) $ がラグランジュ乗数として機能することで維持され、$ I(t) $ は動的法則によって一意に定まる。
• 2次の場合、ヘッセ行列 $ D^2u(t,x) $ は時間に一様にゼロから離れており、有界である。解は定常状態に漸近的に収束する。
• 大きな $ t $ に対して、最適軌道 $ \overline{x}(t) $ は $ A_1^{-1}b $ に収束し、$ I(t) $ は $ I_0 + \frac{1}{2}A_1^{-1}b \cdot b $ に収束する。これにより、長期的挙動が確認される。
• 本手法により、粘性解の漸近的展開が可能となり、微小変異モデルにおける表現型分布の近似が可能になる。