[论文解读] Continuous versus discrete structures II -- Discrete Hamiltonian systems and Helmholtz conditions

Jacky Cresson, Frédéric Pierret|arXiv (Cornell University)|Jan 13, 2015

Numerical methods for differential equations参考文献 13被引用 4

一句话总结

本文在离散嵌入框架内建立了离散哈密顿系统，为连续哈密顿力学提供了自洽的离散类比。它推导出离散系统成为哈密顿系统的显式赫尔姆霍兹条件，证明仅有特定离散系统——如满足特定步长约束的线性或修正谐振子——才具有哈密顿结构，从而解决了离散哈密顿设定下的逆问题。

ABSTRACT

We define discrete Hamiltonian systems in the framework of discrete embeddings. An explicit comparison with previous attempts is given. We then solve the discrete Helmholtz's inverse problem for the discrete calculus of variation in the Hamiltonian setting. Several applications are discussed.

研究动机与目标

• 通过离散嵌入定义一致的离散哈密顿系统框架，确保与连续哈密顿力学的一致性。
• 解决哈密顿设定下的离散赫尔姆霍兹逆问题，识别出离散系统源自变分原理的条件。
• 将所提框架与先前方法（如West、Opris）进行比较，突出其在公式化与结构一致性方面的差异。
• 提供离散系统（如线性、阻尼或修正谐振子）具备哈密顿结构的显式条件。

提出的方法

• 使用离散嵌入形式化定义时间网格 T = hZ ∩ [a,b] 上的离散导数 Δ 与 ∇，以及离散辛标量积。
• 应用离散分部积分与离散反导数，推导离散哈密顿原理与变分结构。
• 通过要求离散微分算子关于离散辛形式自伴，推导出离散赫尔姆霍兹条件（CH1 与 CH2）。
• 引入有限差分算子的离散伴随算子，并利用其刻画哈密顿结构所必需的自伴性条件。
• 将该框架应用于具体系统：牛顿方程、线性阻尼与修正谐振子，通过 ρ(P) = P - h∇P 的代换分析一致性。
• 利用向量场的雅可比矩阵必须对称（CH1）且系统必须保持辛结构（CH2）的条件，验证哈密顿实现性。

实验结果

研究问题

• RQ1 在何种条件下，离散二阶差分系统可表示为离散哈密顿系统？
• RQ2 离散嵌入框架如何确保离散与连续哈密顿力学之间的一致性？
• RQ3 离散系统成为哈密顿系统的必要且充分赫尔姆霍兹型条件是什么？它们与连续情形有何不同？
• RQ4 为何某些系统（如线性阻尼或修正谐振子）无法具备离散哈密顿结构？
• RQ5 离散赫尔姆霍兹条件能否应用于具有非线性或非对称项的系统？会引出何种约束？

主要发现

• 由牛顿方程经离散嵌入导出的系统 (68)，当且仅当势能 U 为线性函数，即 U''(Q) = 0 时，具备哈密顿结构。
• 对于线性阻尼系统 (70)，仅当阻尼系数 γ = 0 时，离散哈密顿结构才存在，与连续情形一致。
• 对于系统 (73)，哈密顿结构仅在 γ = h/m 时存在，表明时间步长 h 必须与系统参数匹配，以保证哈密顿一致性。
• 修正谐振子 (74) 不满足赫尔姆霍兹条件 (CH1)，因其向量场的雅可比矩阵不对称，从而证明其不具备离散哈密顿结构。
• 离散赫尔姆霍兹条件 (CH1 与 CH2) 对于离散系统在哈密顿设定下源自变分原理而言，既是必要也是充分条件。
• 该框架可推广至时间尺度微积分，并可能通过 Bismut 的随机变分法扩展至随机哈密顿系统。