第一次接触 Hamilton 力学时，常见的感受是：公式看起来很整齐，但不知道它到底比牛顿方程多说了什么。再往后看到 Poisson bracket，又容易觉得这只是一个“定义出来的代数运算”，和动力系统本身没有直接关系。其实这两件事是连在一起的。

Hamilton 形式做的事情可以概括成一句话：

不再直接写“受力等于质量乘加速度”，而是把动力系统写成一种由能量驱动的几何演化。

而 Poisson Bracket 则是在这套结构里回答另一个问题：

给定一个物理量，它随时间怎么变？

这篇文章只讲最基本的有限维情形。目标有三个：

• 解释 Hamilton 方程到底在描述什么；
• 说明 Poisson Bracket 为什么会自然出现；
• 顺手说明它为什么和守恒律、PDE 结构有关。

从相空间开始

在 Hamilton 力学里，一个系统的状态通常不用位置和速度写，而是用位置和动量写成

$$(q,p).$$

这里

• $q$ 表示位置；
• $p$ 表示动量。

所有可能的 $(q,p)$ 组成的空间，叫做相空间。对一个 $n$ 自由度系统来说，$q\in\mathbb{R}^n$，$p\in\mathbb{R}^n$，所以相空间是 $2n$ 维的。

这一步看起来只是换变量，但含义其实很重要。牛顿方程通常写成二阶方程，而 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(q,p)$，系统的演化规律就能直接从它读出来。

Hamilton 方程

Hamilton 力学把系统的时间演化写成

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

这就是 Hamilton 方程。

它的结构很值得注意：

• 位置的变化率由 $H$ 对动量的导数决定；
• 动量的变化率由 $H$ 对位置的导数决定，并且带一个负号。

这说明动力学不再是“分别给 $q$ 和 $p$ 写方程”，而是由一个统一的标量函数 $H$ 同时生成。

用最经典的例子看一遍

若

$$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).$$

第一式就是

$$v=\frac{p}{m},$$

第二式则是

$$F=-\nabla V.$$

再结合 $p=m\dot q$，就能恢复熟悉的牛顿方程。所以 Hamilton 形式并不是另一套物理学，而是对同一个动力系统的另一种写法。它真正的价值在于：这种写法把系统背后的几何结构暴露了出来。

Poisson Bracket 为什么会出现

到这里自然会冒出一个问题。

如果系统状态由 Hamilton 方程控制，那么一个一般的物理量

$$F(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。

因此 Hamilton 系统中任意物理量的演化都可以统一写成

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

这句话非常重要。它说明：

Hamiltonian $H$ 不是只负责“能量解释”，它还是生成时间演化的那个函数。

Hamilton 方程可以从 Poisson Bracket 反推回来

上面的公式并不是抽象装饰，它真的把原来的动力学统一了。

取

$$F=q_i,$$

则

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

所以

$$\dot q_i=\{q_i,H\}=\frac{\partial H}{\partial p_i}.$$

再取

$$F=p_i,$$

则

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

所以

$$\dot p_i=\{p_i,H\}=-\frac{\partial H}{\partial q_i}.$$

也就是说，Hamilton 方程本身就是

$$\dot F=\{F,H\}$$

在坐标函数 $q_i,p_i$ 上的特例。

Poisson Bracket 的几个基本性质

Poisson Bracket 之所以重要，不只是因为它能把链式法则压缩成一个符号，还因为它满足一组和动力系统密切相关的结构性质。

双线性

Poisson Bracket 对两个输入都线性：

$$\{aF+bG,H\}=a\{F,H\}+b\{G,H\},$$

以及

$$\{F,aG+bK\}=a\{F,G\}+b\{F,K\}.$$

这保证了它确实像一个“可计算的代数运算”。

反对称性

$$\boxed{ \{F,G\}=-\{G,F\} }$$

这条性质看起来简单，但它立刻带出一个核心结论。取 $F=G=H$，就有

$$\{H,H\}=0,$$

于是

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

也就是说，在没有显含时间的 Hamilton 系统里，Hamiltonian 自动守恒。能量守恒不是额外加上的，而是这套结构内生的结果。

Leibniz 规则

Poisson Bracket 对乘积满足类似导数的规则：

$$\{FG,H\}=F\{G,H\}+G\{F,H\}.$$

这说明它和“时间导数”之间的行为非常一致。

Jacobi 恒等式

$$\boxed{ \{F,\{G,H\}\} + \{G,\{H,F\}\} + \{H,\{F,G\}\} =0 }$$

这是 Poisson 结构最深的一条性质。它保证这种括号不是随便写出来的双线性算子，而是和几何上的相空间结构真正兼容。

从矩阵形式看，结构会更清楚

把相空间坐标记成

$$z=(q,p),$$

再引入标准辛矩阵

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

那么 Poisson Bracket 可以写成

$$\boxed{ \{F,G\}=\nabla F^T J \nabla G }$$

相应地，Hamilton 方程就变成

$$\dot z=J\nabla H.$$

写成这个形式以后，Hamilton 系统的几何味道就更明显了：演化方向由梯度 $\nabla H$ 决定，但真正推动系统流动的是那个反对称矩阵 $J$。这就是标准辛几何结构在坐标里的体现。

为什么说 Casimir 不依赖具体 Hamiltonian

有些量的守恒并不是因为某个特定 Hamiltonian 的形式，而是因为 Poisson 结构本身。

如果一个量 $C$ 满足

$$\boxed{ \{C,G\}=0 \qquad \forall G }$$

那么 $C$ 称为一个 Casimir。

它的含义是：无论系统选哪个 Hamiltonian，只要 Poisson 结构不变，这个量都会守恒。于是 Casimir 反映的就不是“某个能量函数的特殊性质”，而是几何结构本身的退化方向或不变量。

这类量在流体和等离子体问题里尤其常见。

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

前面讲的是有限维系统，也就是状态由有限个坐标 $(q,p)$ 描述的情形。到了 PDE 里，状态变量不再是有限维向量，而是函数，例如

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

这时普通偏导数会换成泛函导数：

$$\frac{\partial}{\partial q} \quad\longrightarrow\quad \frac{\delta}{\delta u}.$$

但核心思路没有变：仍然是

• 先写出一个 Hamiltonian；
• 再通过某种 Poisson 结构生成演化。

换句话说，PDE 里的 Hamilton 结构不是“完全不同的理论”，而是把有限维的思路推广到无穷维函数空间。

以浅水方程为例

浅水方程里，Hamiltonian 通常对应总能量：

$$\boxed{ H = \int \left( \frac12 h|u|^2 + \frac12 gh^2 \right)\,dA }$$

其中

• 第一项是动能；
• 第二项是重力势能。

如果继续往 Hamilton 结构的方向写下去，就需要给出一个适合浅水系统的 Poisson Bracket。此时括号不再是有限维里那种对 $(q,p)$ 的简单偏导，而会涉及速度、厚度、位涡等变量之间的耦合。

在离散浅水系统中，常见的形式会写成类似

$$\{F,H\} = \left\langle \frac{\delta F}{\delta u_h}, -q_h \left( \frac{\delta H}{\delta u_h} \right)^\perp \right\rangle + \left\langle \nabla\cdot\frac{\delta F}{\delta u_h}, \frac{\delta H}{\delta D_h} \right\rangle - \left\langle \frac{\delta F}{\delta D_h}, \nabla\cdot\frac{\delta H}{\delta u_h} \right\rangle.$$

这类表达式看起来复杂，但本质仍然是同一件事：

Hamiltonian 负责告诉你系统储存了什么能量，Poisson 结构负责告诉你这些能量怎样在各个变量之间交换，并生成整个动力学。

也正因为如此，Hamilton 结构常常和下面这些性质绑在一起出现：

• 能量守恒；
• 位涡结构；
• 涡量或 enstrophy 之类的不变量；
• 数值格式中的结构保持性。

Almost-Poisson bracket 是什么

在很多 PDE 或离散系统里，人们会先写出一个“看起来像 Poisson Bracket”的东西，但还没有证明它满足 Jacobi 恒等式。这时更准确的叫法是 almost-Poisson bracket。

也就是说，它可能已经满足：

• 双线性；
• 反对称性；
• Leibniz 规则。

但如果 Jacobi 恒等式还没证出来，就不能直接把它叫成真正的 Poisson Bracket。这个区分在结构保持离散里很重要，因为很多守恒性质和几何解释最终都依赖 Jacobi 恒等式是否成立。

小结

把全文压缩成最核心的几句话，大概就是：

• Hamiltonian 描述系统的能量；
• Hamilton 方程把动力学写成相空间中的一阶流；
• Poisson Bracket 把“任意物理量如何演化”统一写成 $\dot F=\{F,H\}$；
• 反对称性带来能量守恒，Jacobi 恒等式保证这套结构真正是几何性的；
• 到 PDE 里，变量从有限维坐标变成函数，但“能量 + Poisson 结构生成演化”的主线不变。

如果只记一句话，可以记成：

Hamiltonian 决定系统储存什么能量，Poisson Bracket 决定系统怎样演化。

版权所有

版权归属：Guisong Wu

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

下一页什么是地球流体力学