Hamilton 力学与 Poisson bracket

约 2790 字大约 9 分钟

2026-05-15

如果你没有系统学过 Hamilton 力学，第一次看到它时很容易有两个疑问：

这不就是把 Newton 力学换了个写法吗？
后面的 Poisson bracket 看起来像硬定义出来的符号，它到底有什么用？

先说结论：Hamilton 力学确实是在重写 Newton 力学，但它不是随便改写，而是把“系统怎么动”这件事和“系统的能量结构”直接绑在了一起。 而 Poisson bracket 的作用，就是把“任意一个物理量随时间怎么变”统一地写出来。

这篇只讲最基本的有限维情形。你可以把它当成一篇“第一次上手”的直观介绍。

从熟悉的 Newton 力学出发

在 Newton 力学里，我们通常想的是：

粒子在哪里；
它速度多大；
它受了什么力；
然后由 $m\ddot q = F$ 算出它接下来怎么运动。

这是很自然的想法，但它有一个特点：方程通常是二阶的，因为里面直接出现了加速度 $\ddot q$ 。

Hamilton 力学换了一个看问题的角度。它不再把重点放在“加速度由什么力决定”，而是改成：

我只要知道系统当前的状态，以及系统的总能量长什么样，就能写出它怎么演化。

这里“状态”不再写成“位置 + 速度”，而是写成“位置 + 动量”。

什么是相空间

Hamilton 力学里，一个系统的状态通常记成

(q,p).

这里

$q$ 是位置；
$p$ 是动量。

所有可能的 $(q,p)$ 组成的空间，叫做相空间。

如果系统有 $n$ 个自由度，那么

q\in\mathbb{R}^n,\qquad p\in\mathbb{R}^n,

所以相空间是 $2n$ 维的。

你可以先别把“相空间”想得太玄。它本质上就是一个“记录系统状态的坐标空间”。
Newton 力学里你会说“给我初始位置和初始速度”；Hamilton 力学里则说“给我初始位置和初始动量”。

Hamiltonian 是什么

Hamiltonian 通常记作

H(q,p).

在最常见的情形里，它就是系统的总能量：

H(q,p)=T(p)+V(q),

其中

$T$ 是动能；
$V$ 是势能。

最经典的例子是单个粒子：

H(q,p)=\frac{|p|^2}{2m}+V(q).

因为动量满足

p=m\dot q,

所以

\frac{|p|^2}{2m}

刚好就是动能。

这里最重要的不是“ $H$ 等于总能量”这个标签，而是下面这句话：

一旦你写出了 $H(q,p)$ ，系统怎么动就能从 $H$ 里直接读出来。

Hamilton 方程

Hamilton 力学把时间演化写成

\dot q=\frac{\partial H}{\partial p}, \qquad \dot p=-\frac{\partial H}{\partial q}.

这就是 Hamilton 方程。

第一次看会觉得有点突然，所以不妨把它翻译成白话：

位置怎么变，由能量对动量的变化率决定；
动量怎么变，由能量对位置的变化率决定，并且带一个负号。

这和 Newton 写法最大的区别是：现在整个动力学由一个函数 $H$ 统一生成。

它和 Newton 方程其实是同一回事

还是看最熟悉的例子：

H(q,p)=\frac{p^2}{2m}+V(q).

代入 Hamilton 方程，得到

\dot q=\frac{\partial H}{\partial p}=\frac{p}{m},

以及

\dot p=-\frac{\partial H}{\partial q}=-\nabla V(q).

第一条就是

p=m\dot q,

第二条就是“力等于势能的负梯度”：

F=-\nabla V(q).

再把两句合起来，你就会回到熟悉的 Newton 方程。

所以这里没有第二套物理学。
Hamilton 力学和 Newton 力学描述的是同一个系统，只是组织方程的方式不同。

为什么这种改写值得学

如果只是为了算一个小球怎么动，Newton 写法已经够用了。
Hamilton 写法真正厉害的地方在于，它把很多“隐藏在方程背后的结构”变得非常清楚，比如：

能量守恒为什么这么自然；
哪些量会一起演化；
为什么很多系统能写成统一的几何形式；
为什么这种结构能推广到 PDE 和结构保持数值格式。

所以它不只是“改个记号”，而是把动力系统整理成了一个更有组织的框架。

Poisson bracket

讲到这里，自然会问一个问题：

系统状态是 $(q,p)$ ，那如果我关心的不是 $q$ 或 $p$ 本身，而是某个一般的物理量

F(q,p),

它随时间怎么变？

做法其实很朴素：直接用链式法则。

\frac{dF}{dt} = \sum_{i=1}^n \left( \frac{\partial F}{\partial q_i}\dot q_i + \frac{\partial F}{\partial p_i}\dot p_i \right).

再把 Hamilton 方程代进去：

\frac{dF}{dt} = \sum_{i=1}^n \left( \frac{\partial F}{\partial q_i}\frac{\partial H}{\partial p_i} - \frac{\partial F}{\partial p_i}\frac{\partial H}{\partial q_i} \right).

这个式子出现得太频繁了，于是人们干脆给它起个名字：

\boxed{ \{F,G\} = \sum_{i=1}^n \left( \frac{\partial F}{\partial q_i}\frac{\partial G}{\partial p_i} - \frac{\partial F}{\partial p_i}\frac{\partial G}{\partial q_i} \right) }

这就是 Poisson bracket。

于是上面的时间演化立刻变成

\boxed{ \frac{dF}{dt}=\{F,H\} }

Poisson bracket 不是凭空发明出来的，而是链式法则加上 Hamilton 方程的自然结果。

这个公式为什么重要

公式

\frac{dF}{dt}=\{F,H\}

可以理解成:

只要你知道 Hamiltonian $H$ ，那系统里任何物理量 $F$ 的时间变化都能统一算出来。

也就是说， $H$ 不只是“能量函数”，它还是“生成时间演化的函数”。

如果你取

F=q_i,

那么

\{q_i,H\}=\frac{\partial H}{\partial p_i},

所以

\dot q_i=\{q_i,H\}.

如果取

F=p_i,

那么

\{p_i,H\}=-\frac{\partial H}{\partial q_i},

所以

\dot p_i=\{p_i,H\}.

这说明 Hamilton 方程本身，其实就是

\dot F=\{F,H\}

在最基本坐标函数上的特例。

Poisson bracket 的性质

1. 反对称性

\{F,G\}=-\{G,F\}.

把 $F=G=H$ 代进去，就得到

\{H,H\}=0.

于是

\dot H=\{H,H\}=0.

这说明在没有显含时间的情形下，Hamiltonian 自己就是守恒的。
也就是说，能量守恒在这套框架里几乎是“自带的”。

2. Jacobi 恒等式

\{F,\{G,H\}\} + \{G,\{H,F\}\} + \{H,\{F,G\}\} =0.

这条式子第一次看通常不会有直观感觉。你可以先把它当成一句“结构一致性条件”。
它保证这个 bracket 不是随便凑出来的记号，而是真的和 Hamilton 系统背后的几何结构兼容。

矩阵写法

如果把相空间坐标合在一起记成

z=(q,p),

再引入矩阵

J= \begin{pmatrix} 0&I\\ -I&0 \end{pmatrix},

那么 Hamilton 方程可以压缩成

\dot z=J\nabla H,

Poisson bracket 则可以写成

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

这组写法的好处是：你能一眼看出，系统演化不是随意的，而是被一个固定的反对称结构 $J$ 控制着。

为什么 PDE 里也会出现 Hamilton 结构

到了 PDE，状态变量不再是有限个数，而是函数，比如

u(x),\quad h(x).

这时普通偏导数会换成泛函导数，但主线其实没变，仍然是：

先找系统的总能量；
再用某种 Poisson 结构把演化写出来。

所以 PDE 里的 Hamilton 结构，不是完全另一门新理论，而是把有限维的想法推广到了“函数作为变量”的情形。

以浅水方程做个例子

对浅水方程来说，Hamiltonian 往往就是总能量：

H = \int \left( \frac12 h|u|^2 + \frac12 gh^2 \right)\,dA.

这里

$\frac12 h|u|^2$ 是动能；
$\frac12 gh^2$ 是重力势能。

浅水系统的 Poisson bracket 长什么样

如果把状态变量取成速度 $u(x)$ 和水深 $h(x)$ ，那么浅水系统常见的一种 Poisson bracket 写成

\boxed{ \{F,G\} = -\int q\, \frac{\delta F}{\delta u}\cdot \left(\frac{\delta G}{\delta u}\right)^\perp \,dA + \int \left( \nabla\cdot\frac{\delta F}{\delta u}\, \frac{\delta G}{\delta h} - \nabla\cdot\frac{\delta G}{\delta u}\, \frac{\delta F}{\delta h} \right)\,dA }

这里

$F,G$ 是依赖于 $u,h$ 的泛函；
$\delta F/\delta u,\delta F/\delta h$ 是泛函导数；
$q=(\zeta+f)/h$ 是浅水位涡；
$(\cdot)^\perp$ 表示向量逆时针旋转 $90^\circ$ 。

这个公式第一次看通常会有点重，但它其实可以粗略拆成两部分去理解：

第一项和位涡 $q$ 有关，负责旋转结构；
后面两项把速度 $u$ 和层厚 $h$ 耦合起来。

所以它不是随便拼起来的复杂表达式，而是在编码浅水系统最核心的两件事：旋转和 层厚输运。

它怎样生成浅水方程

Hamilton 结构的用法和有限维时完全一样，仍然是

\partial_t F = \{F,H\}.

对浅水系统，Hamiltonian 取总能量

H = \int \left( \frac12 h|u|^2 + \frac12 gh^2 \right)\,dA.

它对两个状态变量的泛函导数分别是

\frac{\delta H}{\delta u}=hu, \qquad \frac{\delta H}{\delta h}=\frac12|u|^2+gh.

把它们代回 bracket 里，就会得到浅水方程的 vector-invariant form：

\partial_t u + q(hu)^\perp + \nabla\left(\frac12|u|^2+gh\right)=0,

\partial_t h + \nabla\cdot(hu)=0.

再用

q=\frac{\zeta+f}{h}

化简第一式，就得到更熟悉的写法

\partial_t u + (\zeta+f)u^\perp + \nabla\left(\frac12|u|^2+gh\right)=0.

也就是说，这个 bracket 的作用正是：

把旋转效应组织成 $(\zeta+f)u^\perp$ ；
把质量守恒组织成 $\partial_t h+\nabla\cdot(hu)=0$ ；
让整个系统由同一个能量泛函 $H$ 统一生成。

为什么它不是正则 bracket

这里的 bracket 和有限维粒子系统里

\{F,G\} = \sum_i \left( \frac{\partial F}{\partial q_i}\frac{\partial G}{\partial p_i} - \frac{\partial F}{\partial p_i}\frac{\partial G}{\partial q_i} \right)

长得很不一样。原因是浅水系统的变量 $(u,h)$ 不是一组标准正则坐标。

这类 bracket 通常叫做 noncanonical Poisson bracket。
“noncanonical” 的意思不是它不合法，而是说它不是最简单的正则形式；它仍然有 Poisson 结构，只是结构矩阵本身依赖于状态变量，比如这里就显式出现了 $q$ 和 $h$ 。

也正因为它是非正则的，很多流体里特有的不变量才会自然出现，例如 Casimir、位涡积分和 enstrophy 一类量。

这件事直观上在说什么

Hamiltonian 说清楚“系统把能量存在哪里”，Poisson 结构说清楚“这些变量怎样交换和传递这些能量”。

对浅水系统来说，这句话可以再具体一点：

Hamiltonian 告诉你动能和重力势能怎么存储；
Poisson bracket 告诉你速度场、层厚和位涡怎样耦合演化。

这也是为什么在流体、浅水方程、结构保持离散这些问题里，Hamilton 结构会反复出现。

小结

Hamilton 力学不是另一套物理学，而是 Newton 力学的一种更有结构的写法；
状态改成用 $(q,p)$ 描述，也就是位置和动量；
Hamiltonian $H$ 通常是总能量；
Hamilton 方程告诉你，系统演化可以直接从 $H$ 读出来；
Poisson bracket 则把“任意物理量怎么随时间变化”统一写成 $\dot F=\{F,H\}$ 。

Hamiltonian 负责描述能量，Poisson bracket 负责描述演化。

版权所有

版权归属：Guisong Wu

许可证：署名 4.0 国际 (CC-BY-4.0)