[论文解读] 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

一句话总结

本文建立了具有最大值约束的时变Hamilton-Jacobi方程的粘性解的唯一性与构造性存在性，该约束确保解的最大值保持为零。通过利用动态规划原理并推导解最大值的常微分方程（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.

研究动机与目标

• 解决长期悬而未决的开放问题：在解的上确界被强制保持为零的约束下，Hamilton-Jacobi方程的唯一性问题。
• 为约束系统提供一种构造性存在性证明，避免先前工作中依赖粘性逼近的方法。
• 在数据满足凹性假设的条件下，证明解为经典解，而不仅仅是粘性解。
• 支持Hamilton-Jacobi方法在复杂选择与突变种群动力学模型中的发展。
• 在小突变极限下，实现更强的收敛结果与渐近展开，以改进底层选择-突变模型的分析。

提出的方法

• 作者利用动态规划原理，在初始数据与函数 $ R(x,I) $ 满足凹性假设的条件下，证明Hamilton-Jacobi方程的解为经典解。
• 他们推导了解在空间上最大值的常微分方程，将完整问题视为一个涉及 $ u(t,x) $ 与调节器 $ I(t) $ 的非标准初值问题。
• 通过将 $ I(t) $ 建模为拉格朗日乘子，使其随时间调整以维持约束 $ \max_x u(t,x) = 0 $。
• 提出一种构造性方法，直接求解方程而不使用粘性逼近，其正则性依赖于凹性。
• 唯一性证明基于对两个解的二阶导数进行比较，并利用特征曲线控制其差值。
• 在二次型情形下，通过Euler-Lagrange方程与矩阵指数显式求解，验证了理论结果。

实验结果

研究问题

• RQ1 在 $ R $ 与 $ u_0 $ 满足凹性假设的条件下，具有最大值约束的Hamilton-Jacobi方程是否具有唯一解？
• RQ2 能否在不使用粘性逼近（如先前工作所采用的）的前提下，提供一种构造性存在性证明？
• RQ3 在凹性假设下，动态规划原理如何导致解的经典正则性？
• RQ4 在长时间极限下，解与约束函数 $ I(t) $ 的行为如何，特别是在二次型情形下？
• RQ5 解的正则性与唯一性是否可用于加强底层选择-突变模型的收敛性结果？

主要发现

• 在 $ R $ 与 $ u_0 $ 满足凹性假设的条件下，由于动态规划原理的作用，约束Hamilton-Jacobi方程的解 $ u(t,x) $ 为经典解，而不仅仅是粘性解。
• 通过证明解的最大值按一个明确定义的常微分方程演化，唯一性得到证明：任意两个解必须完全重合。
• 约束 $ \max_x u(t,x) = 0 $ 通过 $ I(t) $ 实现，$ I(t) $ 作为拉格朗日乘子，其演化由动力学唯一确定。
• 在二次型情形下，Hessian矩阵 $ 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 $，确认了长期行为。
• 该方法允许对粘性解进行渐近展开，从而可近似小突变模型中的表型分布。