在分析力学中，初积分（First Integral）是一个贯穿理论结构始终的核心概念。它不仅是求解运动方程的技术工具，更是连接对称性、守恒律、相空间几何与可积性理论的枢纽。以下从定义出发，逐步揭示初积分在分析力学中扮演的多重角色。
一、初积分的定义
对于哈密尔顿正则方程
q ˙ k = ∂ H ∂ p k , p ˙ k = − ∂ H ∂ q k \dot{q}_k = \frac{\partial H}{\partial p_k}, \quad \dot{p}_k = -\frac{\partial H}{\partial q_k} q˙k=∂pk∂H,p˙k=−∂qk∂H
若函数 f ( q , p , t ) f(q, p, t) f(q,p,t) 沿任意真实运动轨迹满足
d f d t = ∂ f ∂ t + ∑ k ( ∂ f ∂ q k q ˙ k + ∂ f ∂ p k p ˙ k ) = 0 \frac{df}{dt} = \frac{\partial f}{\partial t} + \sum_k\left(\frac{\partial f}{\partial q_k}\dot{q}_k + \frac{\partial f}{\partial p_k}\dot{p}_k\right) = 0 dtdf=∂t∂f+k∑(∂qk∂fq˙k+∂pk∂fp˙k)=0
则称 f = 常数 f = \text{常数} f=常数 为系统的一个初积分。它意味着系统演化被约束在相空间的一张不变曲面上。
最常见的例子包括：当 ∂ H / ∂ t = 0 \partial H/\partial t = 0 ∂H/∂t=0 时， H = E H = E H=E 为能量积分；当 q k q_k qk 为循环坐标（ ∂ H / ∂ q k = 0 \partial H/\partial q_k = 0 ∂H/∂qk=0）时， p k = c k p_k = c_k pk=ck 为广义动量守恒。
我们来严格重新给出初积分的定义，并基于此定义逐条论证 ∂ H / ∂ t = 0 \partial H/\partial t = 0 ∂H/∂t=0 和 ∂ H / ∂ q k = 0 \partial H/\partial q_k = 0 ∂H/∂qk=0 如何与之匹配。
1.1. 初积分的严格定义
考虑由 2 n 2n 2n 个一阶常微分方程组成的哈密顿系统：
q ˙ k = ∂ H ∂ p k , p ˙ k = − ∂ H ∂ q k , k = 1 , 2 , ... , n \dot{q}_k = \frac{\partial H}{\partial p_k}, \quad \dot{p}_k = -\frac{\partial H}{\partial q_k}, \quad k = 1, 2, \dots, n q˙k=∂pk∂H,p˙k=−∂qk∂H,k=1,2,...,n
定义（初积分 / 首次积分 / First Integral）：
设函数 Φ ( q 1 , ... , q n ; p 1 , ... , p n ; t ) \Phi(q_1, \dots, q_n; p_1, \dots, p_n; t) Φ(q1,...,qn;p1,...,pn;t) 在相空间的某个开区域 D ⊂ R 2 n + 1 D \subset \mathbb{R}^{2n+1} D⊂R2n+1 上连续可微。若将系统的任意解 ( q ( t ) , p ( t ) ) (q(t), p(t)) (q(t),p(t)) 代入 Φ \Phi Φ 后，所得函数 Φ ( q ( t ) , p ( t ) , t ) \Phi(q(t), p(t), t) Φ(q(t),p(t),t) 恒等于某一常数（该常数可依赖于具体解，即由初始条件确定），即：
d d t Φ ( q ( t ) , p ( t ) , t ) ≡ 0 \frac{d}{dt}\Phi(q(t), p(t), t) \equiv 0 dtdΦ(q(t),p(t),t)≡0则称 Φ \Phi Φ 为该哈密顿系统的一个初积分。
等价表述：
d Φ d t ∣ 沿正则方程 = ∂ Φ ∂ t + ∑ k = 1 n ( ∂ Φ ∂ q k ∂ H ∂ p k − ∂ Φ ∂ p k ∂ H ∂ q k ) = 0 \frac{d\Phi}{dt}\bigg|{\text{沿正则方程}} = \frac{\partial \Phi}{\partial t} + \sum{k=1}^n \left(\frac{\partial \Phi}{\partial q_k}\frac{\partial H}{\partial p_k} - \frac{\partial \Phi}{\partial p_k}\frac{\partial H}{\partial q_k}\right) = 0 dtdΦ 沿正则方程=∂t∂Φ+k=1∑n(∂qk∂Φ∂pk∂H−∂pk∂Φ∂qk∂H)=0
即：
∂ Φ ∂ t + { Φ , H } = 0 \frac{\partial \Phi}{\partial t} + \{\Phi, H\} = 0 ∂t∂Φ+{Φ,H}=0
其中 { Φ , H } \{\Phi, H\} {Φ,H} 为泊松括号。
1.2. ∂ H / ∂ t = 0 \partial H/\partial t = 0 ∂H/∂t=0 与初积分定义的匹配
1.2.1 候选函数
取 Φ = H ( q 1 , ... , q n ; p 1 , ... , p n ; t ) \Phi = H(q_1, \dots, q_n; p_1, \dots, p_n; t) Φ=H(q1,...,qn;p1,...,pn;t)。
1.2.2 验证初积分条件
计算 Φ = H \Phi = H Φ=H 的全导数：
d H d t = ∂ H ∂ t + ∑ k = 1 n ( ∂ H ∂ q k q ˙ k + ∂ H ∂ p k p ˙ k ) \frac{dH}{dt} = \frac{\partial H}{\partial t} + \sum_{k=1}^n \left(\frac{\partial H}{\partial q_k}\dot{q}_k + \frac{\partial H}{\partial p_k}\dot{p}_k\right) dtdH=∂t∂H+k=1∑n(∂qk∂Hq˙k+∂pk∂Hp˙k)
代入正则方程 q ˙ k = ∂ H ∂ p k \dot{q}_k = \frac{\partial H}{\partial p_k} q˙k=∂pk∂H， p ˙ k = − ∂ H ∂ q k \dot{p}_k = -\frac{\partial H}{\partial q_k} p˙k=−∂qk∂H：
d H d t = ∂ H ∂ t + ∑ k = 1 n ( ∂ H ∂ q k ∂ H ∂ p k + ∂ H ∂ p k ( − ∂ H ∂ q k ) ) = ∂ H ∂ t + 0 = ∂ H ∂ t \frac{dH}{dt} = \frac{\partial H}{\partial t} + \sum_{k=1}^n \left(\frac{\partial H}{\partial q_k}\frac{\partial H}{\partial p_k} + \frac{\partial H}{\partial p_k}\left(-\frac{\partial H}{\partial q_k}\right)\right) = \frac{\partial H}{\partial t} + 0 = \frac{\partial H}{\partial t} dtdH=∂t∂H+k=1∑n(∂qk∂H∂pk∂H+∂pk∂H(−∂qk∂H))=∂t∂H+0=∂t∂H
1.2.3 应用条件 ∂ H / ∂ t = 0 \partial H/\partial t = 0 ∂H/∂t=0
若：
∂ H ∂ t = 0 \frac{\partial H}{\partial t} = 0 ∂t∂H=0
则：
d H d t = 0 \frac{dH}{dt} = 0 dtdH=0
1.2.4 匹配结论
| 初积分定义要求 | H H H 的满足情况 |
|---|---|
| Φ \Phi Φ 是 ( q , p , t ) (q, p, t) (q,p,t) 的连续可微函数 | H ( q , p , t ) H(q, p, t) H(q,p,t) 满足 |
| 沿任意解 轨道 d Φ d t = 0 \frac{d\Phi}{dt} = 0 dtdΦ=0 | d H d t = ∂ H ∂ t = 0 \frac{dH}{dt} = \frac{\partial H}{\partial t} = 0 dtdH=∂t∂H=0，对任意解成立 |
| Φ \Phi Φ 沿解保持常数 | H ( q ( t ) , p ( t ) ) = H ( q 0 , p 0 ) = E H(q(t), p(t)) = H(q_0, p_0) = E H(q(t),p(t))=H(q0,p0)=E |
严格结论： 当 ∂ H / ∂ t = 0 \partial H/\partial t = 0 ∂H/∂t=0 时， H H H 满足初积分的全部定义条件，因此 H H H 是一个初积分，称为能量积分 或广义能量积分。
1.3. ∂ H / ∂ q k = 0 \partial H/\partial q_k = 0 ∂H/∂qk=0 与初积分定义的匹配
1.3.1 候选函数
取 Φ = p k \Phi = p_k Φ=pk（第 k k k 个广义动量）。
1.3.2 验证初积分条件
计算 Φ = p k \Phi = p_k Φ=pk 的全导数：
d p k d t = ∂ p k ∂ t + ∑ j = 1 n ( ∂ p k ∂ q j q ˙ j + ∂ p k ∂ p j p ˙ j ) \frac{dp_k}{dt} = \frac{\partial p_k}{\partial t} + \sum_{j=1}^n \left(\frac{\partial p_k}{\partial q_j}\dot{q}_j + \frac{\partial p_k}{\partial p_j}\dot{p}_j\right) dtdpk=∂t∂pk+j=1∑n(∂qj∂pkq˙j+∂pj∂pkp˙j)
由于 p k p_k pk 与 q j q_j qj、 t t t 无关，且 ∂ p k ∂ p j = δ k j \frac{\partial p_k}{\partial p_j} = \delta_{kj} ∂pj∂pk=δkj：
d p k d t = 0 + ∑ j = 1 n ( 0 ⋅ q ˙ j + δ k j ⋅ p ˙ j ) = p ˙ k \frac{dp_k}{dt} = 0 + \sum_{j=1}^n \left(0 \cdot \dot{q}j + \delta{kj} \cdot \dot{p}_j\right) = \dot{p}_k dtdpk=0+j=1∑n(0⋅q˙j+δkj⋅p˙j)=p˙k
由正则方程：
p ˙ k = − ∂ H ∂ q k \dot{p}_k = -\frac{\partial H}{\partial q_k} p˙k=−∂qk∂H
因此：
d p k d t = − ∂ H ∂ q k \frac{dp_k}{dt} = -\frac{\partial H}{\partial q_k} dtdpk=−∂qk∂H
1.3.3 应用条件 ∂ H / ∂ q k = 0 \partial H/\partial q_k = 0 ∂H/∂qk=0
若：
∂ H ∂ q k = 0 \frac{\partial H}{\partial q_k} = 0 ∂qk∂H=0
则：
d p k d t = 0 \frac{dp_k}{dt} = 0 dtdpk=0
1.3.4 匹配结论
| 初积分定义要求 | p k p_k pk 的满足情况 |
|---|---|
| Φ \Phi Φ 是 ( q , p , t ) (q, p, t) (q,p,t) 的连续可微函数 | p k p_k pk 作为相空间坐标，满足 |
| 沿任意解 轨道 d Φ d t = 0 \frac{d\Phi}{dt} = 0 dtdΦ=0 | d p k d t = − ∂ H ∂ q k = 0 \frac{dp_k}{dt} = -\frac{\partial H}{\partial q_k} = 0 dtdpk=−∂qk∂H=0，对任意解成立 |
| Φ \Phi Φ 沿解保持常数 | p k ( t ) = p k ( 0 ) = p k 0 p_k(t) = p_k(0) = p_{k0} pk(t)=pk(0)=pk0 |
严格结论： 当 ∂ H / ∂ q k = 0 \partial H/\partial q_k = 0 ∂H/∂qk=0 时， p k p_k pk 满足初积分的全部定义条件，因此 p k p_k pk 是一个初积分，称为循环积分 或广义动量积分。
1.4. 统一视角：泊松括号判据的验证
1.4.1 对 H H H 的验证
∂ H ∂ t + { H , H } = ∂ H ∂ t + 0 = ∂ H ∂ t \frac{\partial H}{\partial t} + \{H, H\} = \frac{\partial H}{\partial t} + 0 = \frac{\partial H}{\partial t} ∂t∂H+{H,H}=∂t∂H+0=∂t∂H
当 ∂ H ∂ t = 0 \frac{\partial H}{\partial t} = 0 ∂t∂H=0 时，上式为零， H H H 满足初积分条件。
1.4.2 对 p k p_k pk 的验证
∂ p k ∂ t + { p k , H } = 0 + ∑ j = 1 n ( ∂ p k ∂ q j ∂ H ∂ p j − ∂ p k ∂ p j ∂ H ∂ q j ) \frac{\partial p_k}{\partial t} + \{p_k, H\} = 0 + \sum_{j=1}^n \left(\frac{\partial p_k}{\partial q_j}\frac{\partial H}{\partial p_j} - \frac{\partial p_k}{\partial p_j}\frac{\partial H}{\partial q_j}\right) ∂t∂pk+{pk,H}=0+j=1∑n(∂qj∂pk∂pj∂H−∂pj∂pk∂qj∂H)
= ∑ j = 1 n ( 0 ⋅ ∂ H ∂ p j − δ k j ⋅ ∂ H ∂ q j ) = − ∂ H ∂ q k = \sum_{j=1}^n \left(0 \cdot \frac{\partial H}{\partial p_j} - \delta_{kj} \cdot \frac{\partial H}{\partial q_j}\right) = -\frac{\partial H}{\partial q_k} =j=1∑n(0⋅∂pj∂H−δkj⋅∂qj∂H)=−∂qk∂H
当 ∂ H ∂ q k = 0 \frac{\partial H}{\partial q_k} = 0 ∂qk∂H=0 时，上式为零， p k p_k pk 满足初积分条件。
1.5. 总结：匹配关系的完整图景
┌─────────────────────────────────────────────────────────────┐
│ 初积分定义 │
│ Φ(q,p,t) 沿任意解满足 dΦ/dt = ∂Φ/∂t + {Φ,H} = 0 │
├─────────────────────────────────────────────────────────────┤
│ │
│ 条件：∂H/∂t = 0 │
│ ↓ │
│ 候选：Φ = H │
│ ↓ │
│ 验证：dH/dt = ∂H/∂t + {H,H} = 0 + 0 = 0 │
│ ↓ │
│ 结论：H = E 是初积分（能量积分） │
│ │
├─────────────────────────────────────────────────────────────┤
│ │
│ 条件：∂H/∂q_k = 0 │
│ ↓ │
│ 候选：Φ = p_k │
│ ↓ │
│ 验证：dp_k/dt = ∂p_k/∂t + {p_k,H} = 0 - ∂H/∂q_k = 0 │
│ ↓ │
│ 结论：p_k = const 是初积分（循环积分/广义动量积分） │
│ │
└─────────────────────────────────────────────────────────────┘核心要点：
- ∂ H / ∂ t = 0 \partial H/\partial t = 0 ∂H/∂t=0 使 H H H 本身成为初积分，因为 H H H 对 t t t 的显式依赖消失后，其全导数沿任意解恒为零。
- ∂ H / ∂ q k = 0 \partial H/\partial q_k = 0 ∂H/∂qk=0 使 p k p_k pk 成为初积分，因为正则方程直接给出 p ˙ k = 0 \dot{p}_k = 0 p˙k=0，即 p k p_k pk 沿任意解守恒。
- 两者都是充分条件：满足条件即可构造出满足定义的初积分函数。
二、作为求解工具：降维与可积性
初积分的最直接作用是降维 。每找到一个独立初积分，相空间维度降低一维。 2 n 2n 2n 维相空间被逐层"压缩"，最终缩为系统真实的运动轨迹。
若系统拥有 n n n 个相互对合 （泊松括号两两为零）的独立初积分，则称该系统刘维尔可积 。此时可通过正则变换引入作用量-角变量 ，将问题化为求积形式，从而显式解出运动。能量积分 H = E H = E H=E 是最常见也最基本的初积分，它将 2 n 2n 2n 维问题降到 ( 2 n − 1 ) (2n-1) (2n−1) 维的能量面上。
然而，必须强调的是：初积分不等于 q k ( t ) q_k(t) qk(t) 或 p k ( t ) p_k(t) pk(t) 的显式解。它是运动常数，是相空间中的约束条件。它的作用不是直接给出轨迹函数，而是揭示轨道被限制在哪些曲面上。
三、超越求解：初积分的深层物理角色
将初积分的作用限定为"求解方程的技术手段"，是极大的低估。初积分的意义至少体现在以下五个层面：
（1）对称性与守恒律的化身
根据Noether 定理 ，每一个连续对称性对应一个初积分。时间平移对称 → \rightarrow → 能量守恒；空间平移对称 → \rightarrow → 动量守恒；转动对称 → \rightarrow → 角动量守恒。初积分是物理定律不变性的数学体现，它们的存在揭示了系统内在的结构------这是物理本质，而非计算技巧。
（2）可积性与混沌的判据
初积分的数量直接决定系统的定性行为。找到足够多的对合初积分，系统规则运动（轨迹约束在环面上）；反之，系统可能进入混沌 ------对初值极端敏感。初积分在此不是"帮你算出 q ( t ) q(t) q(t)"，而是告诉你这个系统本质上能不能被精确理解。
（3）相空间拓扑的刻画
每个初积分 f ( q , p ) = c f(q,p) = c f(q,p)=c 定义了相空间中的一张不变流形。系统演化始终被限制在这些曲面的交线上。即使永远不计算具体解，仅从初积分的存在性和取值，就能判断轨道的定性行为------是有界还是散射，是闭合还是开放。
（4）数值校验的标准
在数值模拟中，初积分是检验算法的"金标准"。辛算法之所以被推崇，正是因为它能长时间保持能量等初积分。若能量在数值演化中严重漂移，则算法或步长存在问题。
（5）通向量子化的桥梁
经典初积分（特别是作用量变量 J J J）是半经典量子化的基础。Bohr-Sommerfeld 量子化条件 ∮ p d q = n h \oint p\,dq = nh ∮pdq=nh 完全建立在经典初积分之上。没有经典守恒量，量子化无从谈起。
四、实例一：开普勒问题
哈密尔顿量
H = p r 2 2 m + p θ 2 2 m r 2 − k r H = \frac{p_r^2}{2m} + \frac{p_\theta^2}{2mr^2} - \frac{k}{r} H=2mpr2+2mr2pθ2−rk
初积分 1：能量 H = E H = E H=E
- E < 0 ⇒ E < 0 \Rightarrow E<0⇒ 有界运动（椭圆）
- E = 0 ⇒ E = 0 \Rightarrow E=0⇒ 抛物线（恰好逃逸）
- E > 0 ⇒ E > 0 \Rightarrow E>0⇒ 双曲线（自由散射）
仅凭初积分的取值，未解任何微分方程，就已判定轨道的拓扑类型。
初积分 2：角动量 L = p θ = 常数 L = p_\theta = \text{常数} L=pθ=常数
- 运动被限制在垂直于 L ⃗ \vec{L} L 的平面内
- 结合 H = E H = E H=E 构造有效势能 V eff ( r ) = L 2 2 m r 2 − k r V_{\text{eff}}(r) = \frac{L^2}{2mr^2} - \frac{k}{r} Veff(r)=2mr2L2−rk
- 从 V eff V_{\text{eff}} Veff 的图像直接读出：圆轨道半径、振动频率、椭圆轨道的近日点与远日点
初积分 3：Laplace-Runge-Lenz 矢量
A ⃗ = p ⃗ × L ⃗ − m k r ⃗ r \vec{A} = \vec{p} \times \vec{L} - mk\frac{\vec{r}}{r} A =p ×L −mkrr
这是 1 / r 1/r 1/r 势特有的"隐藏"对称性（SO(4) 动力学对称性）的产物。其方向指向椭圆长轴，大小给出偏心率。
三个初积分相互对合，系统完全可积。Bertrand 定理进一步表明：正是初积分的"超额存在"，逆向约束了势函数必须是 1 / r 1/r 1/r 形式，才能保证所有有界轨道闭合。
五、实例二：一维谐振子
哈密尔顿量
H = p 2 2 m + 1 2 m ω 2 q 2 H = \frac{p^2}{2m} + \frac{1}{2}m\omega^2 q^2 H=2mp2+21mω2q2
初积分 H = E H = E H=E 直接在相空间 ( q , p ) (q, p) (q,p) 中给出一族同心椭圆 ：
p 2 2 m E + q 2 2 E / ( m ω 2 ) = 1 \frac{p^2}{2mE} + \frac{q^2}{2E/(m\omega^2)} = 1 2mEp2+2E/(mω2)q2=1
无需解微分方程，相空间几何图像已完全清晰。
引入正则变换到作用量-角变量 ：
J = 1 2 π ∮ p d q = E ω , θ = ω t + θ 0 J = \frac{1}{2\pi}\oint p\,dq = \frac{E}{\omega}, \quad \theta = \omega t + \theta_0 J=2π1∮pdq=ωE,θ=ωt+θ0
哈密尔顿量变为 H = ω J H = \omega J H=ωJ，运动方程彻底平凡化： J ˙ = 0 \dot{J} = 0 J˙=0， θ ˙ = ω \dot{\theta} = \omega θ˙=ω。一个常数，一个匀速转动。
进一步，经典初积分 J J J 直接通向量子化：
J = n ℏ ⇒ E n = n ℏ ω J = n\hbar \quad \Rightarrow \quad E_n = n\hbar\omega J=nℏ⇒En=nℏω
与严格量子力学结果仅差零点能。此处初积分不是求解工具，而是连接经典与量子的桥梁。
六、总结
| 角色层面 | 核心功能 |
|---|---|
| 技术层面 | 降维、辅助求解、判定刘维尔可积性 |
| 物理层面 | 体现对称性与守恒律（Noether 定理） |
| 几何层面 | 刻画相空间的不变流形与拓扑结构 |
| 动力学层面 | 区分规则运动与混沌 |
| 实践层面 | 数值算法的守恒性校验 |
| 理论层面 | 半经典量子化的基石 |
初积分是分析力学中最富有多重身份的概念之一。它既是相空间中的运动常数 ，也是对称性的物理化身 ；既是求解方程的阶梯 ，也是理解系统本质行为的窗口 。掌握初积分，常常比显式解出 q ( t ) q(t) q(t) 更有价值------因为显式解可能是一堆复杂的特殊函数，而初积分直接告诉我们：轨道长什么样，它稳不稳定，能不能闭合，以及如何通向量子世界。