跳到正文
格致开物MATHWIKI

哈密顿向量场

哈密顿向量场(Hamiltonian vector field)由一个哈密顿函数和辛结构共同确定,用来给出相空间中的运动方向。在标准位置—动量坐标 (q,p) 中,它产生哈密顿方程 q˙=∂H∂p,p˙=−∂H∂q. 自治系统沿轨迹保持哈密顿函数的值;在力学模型中,适当的 H 常表示总能量。相空间中的方向场让能量等值线与实际运动连接起来。

从弹簧振子得到相轨迹

取质量为一、弹性系数为四的无阻尼一维振子。位置为 q,动量为 p=q˙,哈密顿函数为 H(q,p)=p22+2q2. 因此 q˙=p、p˙=−4q。初值 q(0)=1,p(0)=0 的解为 q(t)=cos⁡2t,p(t)=−2sin⁡2t. 将其代回,得到 H=2。相轨迹满足 q2+p24=1. 在点 (1,0),向量场为 (0,−4),所以运动沿相平面中的椭圆顺时针进行。椭圆上的点描述同一个系统在不同状态下的位置与动量,不是物体在二维物理空间中的路径。

振子在q,p相平面的能量椭圆,箭头顺时针,初点一零处速度竖直向下
初始能量为二,动量最大绝对值为二。相轨迹同时满足微分方程与能量恒等式。

若初值改为 q(0)=0,p(0)=1,则 q(t)=12sin⁡2t,p(t)=cos⁡2t, 能量变成 1/2,轨迹为较小的椭圆。这一变式改变振幅与能量,线性振子的周期仍为 π。

辛形式与符号约定

在 ℝ2n 中,取标准辛形式 ω=∑j=1ndqj∧dpj. 对光滑函数 H,本文用 ω(XH,Y)=dH(Y)对所有切向量 Y 定义 XH。辛形式的非退化性保证这个向量唯一。若写 XH=(a,b)、Y=(u,v),则 ω(XH,Y)=a⋅v−b⋅u,dH(Y)=Hq⋅u+Hp⋅v. 比较任意 u,v 的系数,得到 a=Hp、b=−Hq,即开头的哈密顿方程。不同教材若把辛形式或定义式取相反符号,应同时调整对应约定。[1]

在一般辛流形上,给定闭的非退化二形式 ω 后,同一定义仍成立。这使坐标中的方程具有内在的几何意义;本文后续计算均在标准正则坐标中进行。

自治能量守恒

若 H=H(q,p) 不显含时间,沿解求导: dHdt=Hq⋅q˙+Hp⋅p˙=Hq⋅Hp−Hp⋅Hq=0. 等价地,dH(XH)=ω(XH,XH)=0,因为辛形式是反对称的。向量场因此切于能量等值面。

若 H 显含时间,则总导数还包含 ∂H/∂t。例如 H(q,p,t)=tp 给出 q˙=t,p˙=0,所以 H(t)=tp(0) 一般随时间改变。自治性正是前述能量守恒计算使用的条件。

正则坐标中的体积保持

设 H∈C2。向量场的散度为 div⁡XH=∑j=1n(Hpjqj−Hqjpj)=0. 混合偏导相等使每项抵消。若局部流 Φt 存在,变分方程与行列式求导给出 ddtdet⁡DΦt=(div⁡XH)∘Φtdet⁡DΦt=0. 初始流为恒等,故 det⁡DΦt=1:相空间中的标准体积保持。这是刘维尔体积定理的正则坐标形式。它与能量守恒描述不同性质;具有显式时间依赖的哈密顿场,在每个时刻仍有相同的零散度计算。

零散度单独还不足以在高维判定某场为哈密顿场。例如在四维坐标 (q1,q2,p1,p2) 中,向量场 (q1,−q2,0,0) 的散度为零。若它来自某个 H,前两分量迫使 H=q1p1−q2p2+h(q1,q2). 其动量方向分量却会包含 −p1−hq1,无法对所有 p1 恒为零,产生矛盾。哈密顿结构比体积保持包含更多条件。

单摆的振荡与旋转

在质量、重力加速度、摆长均取一的无量纲单摆模型中, H(q,p)=p22+1−cos⁡q. 方程为 q˙=p、p˙=−sin⁡q。固定能量 E 后, p=±2(E−1+cos⁡q). 当 0<E<2 时,动量在转向点为零,角度在有限范围内往返。例如 E=1 时,转向角为 ±π/2。沿能量二的轨迹趋近倒立平衡点时,速度趋于零,形成分界轨道;能量大于二时,动量不再变号,运动对应持续旋转。

单摆能量一的闭合振荡轨道、能量二的分界轨道及能量三的旋转支路
图把角度展开在负圆周率至圆周率之间。真实角变量按二倍圆周率周期识别,旋转轨道应在圆柱相空间中理解。

阻尼项会改变结构。例如加上 −ηp、η>0 后,能量变化成为 H˙=−ηp2,原来沿等值线运行的相轨迹会向低能量区域移动。标准无阻尼哈密顿模型的能量守恒与体积保持,不能直接用在这个耗散系统上。

参考文献

  1. ↑ Daniel Liberzon,Calculus of Variations and Optimal Control Theory,§7.1.3“Symplectic geometry and Hamiltonian flows”,式(7.6)–(7.8)。