Hamilton 系统 | Codeboy's Blog

形式陈述

固定整数 $d\ge1$ ，在相空间 $\mathbb R^{2d}$ 中写状态 $z=(q,p)$ ，其中 $q,p\in\mathbb R^d$ 分别表示广义位置与动量。给定连续可微的 Hamiltonian

H:\mathbb R^{2d}\to\mathbb R,

canonical Hamilton 方程为

\dot q=\nabla_pH(q,p), \qquad \dot p=-\nabla_qH(q,p).

令标准辛矩阵

J= \begin{pmatrix} 0&I_d\\ -I_d&0 \end{pmatrix}, \qquad J^{\mathsf T}=-J,

则方程可紧写为常微分方程

\dot z=J\nabla H(z).

若 $H\in C^2$ ，向量场局部 Lipschitz，Picard–Lindelöf 定理便给出每个初值的局部唯一解。只要解存在，由链式法则可见沿轨道的 Hamiltonian 保持不变：

\frac{d}{dt}H(z(t)) =\nabla H(z(t))^{\mathsf T}\dot z(t) =\nabla H(z(t))^{\mathsf T}J\nabla H(z(t))=0,

因为任意实向量 $v$ 对反对称矩阵都满足 $v^{\mathsf T}Jv=0$ 。若 $H$ 显含时间， $H=H(t,q,p)$ ，则沿解有 $dH/dt=\partial H/\partial t$ ，一般不再守恒。

在 $H\in C^2$ 的条件下，局部流 $\Phi_t$ 对初值连续可微，并保持辛形式：

(D\Phi_t(z))^{\mathsf T}J^{-1}D\Phi_t(z)=J^{-1}.

这个恒等式可直接核验。记 $Y(t)=D\Phi_t(z)$ ，变分方程为

\dot Y=J S(t)Y,\qquad S(t)=\nabla^2H(\Phi_t(z)),\qquad Y(0)=I.

Hessian $S(t)$ 对称，故

\frac d{dt}(Y^{\mathsf T}J^{-1}Y) =Y^{\mathsf T}\bigl(SJ^{\mathsf T}J^{-1}+J^{-1}JS\bigr)Y=0.

初值给出所需恒等式。取行列式得到 $(\det Y)^2=1$ ，再由 $\det Y(0)=1$ 与连续性得到 $\det Y=1$ 。局部流在其存在域上以 $\Phi_{-t}$ 为逆，因此相空间体积保持。辛性还保留任意两条切向量的辛配对，比只约束总体积更强。

直觉

梯度 $\nabla H$ 指向 Hamiltonian 增长最快的方向，乘以 $J$ 相当于在每个 $(q_i,p_i)$ 平面旋转四分之一圈。于是速度始终与能量梯度正交，轨道沿着等能面运动而不横穿它。这个几何结构把“守恒”写进了向量场的生成方式，而不是事后观察到的数值巧合。

位置变化由能量对动量的敏感度决定，动量变化则由能量对位置的负梯度决定。对机械系统 $H=T(p)+V(q)$ ，第一式把动量转成速度，第二式就是由势能产生的力。

例子与边界

一维质量—弹簧系统的 Hamiltonian 为

H(q,p)=\frac{p^2}{2m}+\frac{kq^2}{2}, \qquad m,k>0.

Hamilton 方程给出

\dot q=\frac pm, \qquad \dot p=-kq.

消去 $p$ 得 $m\ddot q+kq=0$ 。若初值为 $q(0)=a,p(0)=0$ ，则

q(t)=a\cos(\omega t), \qquad p(t)=-ma\omega\sin(\omega t), \qquad \omega=\sqrt{k/m}.

代回可见 $H(t)=ka^2/2$ ；相轨道是椭圆，而非向平衡点螺旋收缩。

自由粒子取 $H=p^2/(2m)$ 、 $m>0$ ，得到 $p(t)=p_0$ 、 $q(t)=q_0+p_0t/m$ 。能量守恒不意味着状态静止，它只把轨道限制在固定能级上。

边界与失败情形

任意守恒 ODE 不自动是给定坐标下的 canonical Hamilton 系统。三维刚体或带约束系统可能具有 Poisson 结构，需要状态依赖的反对称张量；耗散系统 $\dot q=p,\ \dot p=-q-\gamma p$ （ $\gamma>0$ ）则使机械能满足

\frac{d}{dt}\frac{q^2+p^2}{2}=-\gamma p^2\le0,

其向量场散度为 $-\gamma<0$ ，使相面积收缩，因而不是这组坐标下的标准 Hamilton 流。

辛性与精确能量守恒不是同一性质。辛数值方法精确保持相应的离散辛形式，却未必逐步满足 $H(z_{n+1})=H(z_n)$ ；长期能量误差仍须结合步长、轨道范围和正则性另作估计。反过来，一个专门投影回能量面的算法可能精确保能，却破坏相空间的辛结构。

能量守恒也不保证解全局存在。例如 $H(q,p)=q^2p$ 给出 $\dot q=q^2,\dot p=-2qp$ ，初值 $(1,0)$ 的解为 $(q(t),p(t))=(1/(1-t),0)$ ，能量恒为零却在 $t=1$ 爆破。若 Hamiltonian 的等能面不紧，轨道可能逃向无穷远；要从守恒量推出全局有界，还需能量的强制性等附加条件。平衡点由 $\nabla H=0$ 给出，但它可以是稳定中心、鞍点或更退化类型，不能仅凭“能量不变”判定稳定。

若 $H$ 不够光滑，局部唯一性和可微流都可能失败。能量沿解的链式推导至少需要可用的可微性，辛流结论通常要求更高正则性。

推论与应用

Hamilton 结构连接经典力学、几何光学、分子动力学与偏微分方程的半离散系统。Poisson括号

\{F,G\}=\nabla F^{\mathsf T}J\nabla G

把任意观测量沿流的变化写为 $dF/dt=\{F,H\}$ ；因 $\{H,H\}=0$ ，能量守恒成为反对称性的直接结果。

Hamilton–Jacobi方程把一族相轨道进一步压成标量函数的梯度图像。作用量的标签导数等于动量作用在位置变化上，位置投影可逆时便恢复PDE；谐振子的一族初值在π/2时位置同时坍缩，动量却仍分开，说明相流唯一与经典势函数能否继续是两项不同检查。

在一般辛流形上，辛场与Hamilton场通过 $\iota_X\omega$ 的闭性与恰当性区分。环面平移保持辛形式，却因周期积分非零而没有全局Hamilton函数；因此不能把局部的canonical方程直接当成全局函数已经存在的证据。

结构保持积分器以离散映射逼近 $\Phi_h$ ，重点保留辛几何而不只追求单步截断误差。对长时间轨道，几何性质往往比每一步的能量精确相等更能阻止非物理漂移。

数值积分可以直接要求一步映射保持辛结构，也可以从离散作用量导出这种结构。两条路线都提示：原能量、辛形式与对称性动量是不同对象，离散时需要分别核验。

参考资料

V. I. Arnold, Mathematical Methods of Classical Mechanics, 2nd ed., Springer, 1989，Ch. 8–9。
Jerrold E. Marsden and Tudor S. Ratiu, Introduction to Mechanics and Symmetry, 2nd ed., Springer, 1999，Ch. 2。
Ernst Hairer, Christian Lubich, and Gerhard Wanner, Geometric Numerical Integration, 2nd ed., Springer, 2006，Ch. VI。
Ernst Hairer, Geometric Numerical Integration, Lecture 1: Hamiltonian systems, TU München, 2010，§4, Theorem 2 与变分方程证明： $C^2$ Hamiltonian 的流在其存在域上保持辛形式。

形式陈述 ​

边界与失败情形 ​

反向链接

形式陈述

边界与失败情形