HJ方程（Hamilton-Jacobi Equation）是整个分析力学体系中最深刻的方程。它不是"又一种求解方法"，而是从经典力学内部自然生长出来的、揭示力学与几何光学、量子力学、辛几何之间深层联系的枢纽。
一、数学形式
对于Hamiltonian为 H(q,p,t)H(q, p, t)H(q,p,t) 的系统，HJ方程为：
H(q,∂S∂q,t)+∂S∂t=0\boxed{H\left(q, \frac{\partial S}{\partial q}, t\right) + \frac{\partial S}{\partial t} = 0}H(q,∂q∂S,t)+∂t∂S=0
即：
H(q1,...,qn,∂S∂q1,...,∂S∂qn,t)+∂S∂t=0H\left(q_1, \ldots, q_n, \frac{\partial S}{\partial q_1}, \ldots, \frac{\partial S}{\partial q_n}, t\right) + \frac{\partial S}{\partial t} = 0H(q1,...,qn,∂q1∂S,...,∂qn∂S,t)+∂t∂S=0
其中 S(q,t)S(q, t)S(q,t) 称为Hamilton主函数 （Hamilton's Principal Function），是一个关于 nnn 个广义坐标和时间的标量函数。
在 simplest 的自治系统情形（HHH 不显含 ttt），若进一步分离变量：
H(q,∂S∂q)=E=constH\left(q, \frac{\partial S}{\partial q}\right) = E = \text{const}H(q,∂q∂S)=E=const
则：
S(q,t)=W(q)−EtS(q, t) = W(q) - EtS(q,t)=W(q)−Et
W(q)W(q)W(q) 称为Hamilton特征函数 ，满足定态HJ方程 ：
H(q,∂W∂q)=EH\left(q, \frac{\partial W}{\partial q}\right) = EH(q,∂q∂W)=E
二、产生脉络：从 Lagrange → Hamilton → HJ
HJ方程不是凭空出现的。它在分析力学内部有一条清晰的"进化链"：
第一步：Lagrange力学------变分原理
从Hamilton原理出发：
δ∫t1t2L(q,q˙,t) dt=0\delta \int_{t_1}^{t_2} L(q, \dot{q}, t)\, dt = 0δ∫t1t2L(q,q˙,t)dt=0
得到Euler-Lagrange方程。这是二阶的ODE系统。
第二步：Legendre变换------进入Hamilton力学
通过Legendre变换：
pi=∂L∂q˙i,H(q,p,t)=∑piq˙i−Lp_i = \frac{\partial L}{\partial \dot{q}_i}, \quad H(q, p, t) = \sum p_i\dot{q}_i - Lpi=∂q˙i∂L,H(q,p,t)=∑piq˙i−L
进入 2n2n2n 维相空间，得到一阶 的正则方程：
q˙i=∂H∂pi,p˙i=−∂H∂qi\dot{q}_i = \frac{\partial H}{\partial p_i}, \quad \dot{p}_i = -\frac{\partial H}{\partial q_i}q˙i=∂pi∂H,p˙i=−∂qi∂H
力学从"切丛上的二阶方程"进化到"余切丛上的一阶方程组"。
第三步：正则变换------寻找循环坐标
Hamilton力学的核心技巧是正则变换 ：找新坐标 (Q,P)(Q, P)(Q,P) 使得 K=HnewK = H_{\text{new}}K=Hnew 尽可能简单。
理想情况：所有 QiQ_iQi 都是循环坐标，即 KKK 只依赖于 PPP：
K=K(P),P˙i=−∂K∂Qi=0,Q˙i=∂K∂Pi=constK = K(P), \quad \dot{P}_i = -\frac{\partial K}{\partial Q_i} = 0, \quad \dot{Q}_i = \frac{\partial K}{\partial P_i} = \text{const}K=K(P),P˙i=−∂Qi∂K=0,Q˙i=∂Pi∂K=const
这意味着 Pi=αi=constP_i = \alpha_i = \text{const}Pi=αi=const，Qi=ωit+βiQ_i = \omega_i t + \beta_iQi=ωit+βi，问题完全解决。
第四步：生成函数------正则变换的"母函数"
正则变换可以通过生成函数 （Generating Function）来构造。对于第二类生成函数 F2=S(q,P,t)F_2 = S(q, P, t)F2=S(q,P,t)，变换关系为：
pi=∂S∂qi,Qi=∂S∂Pi,K=H+∂S∂tp_i = \frac{\partial S}{\partial q_i}, \quad Q_i = \frac{\partial S}{\partial P_i}, \quad K = H + \frac{\partial S}{\partial t}pi=∂qi∂S,Qi=∂Pi∂S,K=H+∂t∂S
第五步：关键一跃------令 K=0K = 0K=0
如果能让新Hamiltonian恒为零：
K=0K = 0K=0
则所有新变量都是常数：
P˙i=0,Q˙i=0\dot{P}_i = 0, \quad \dot{Q}_i = 0P˙i=0,Q˙i=0
即 Pi=αiP_i = \alpha_iPi=αi，Qi=βiQ_i = \beta_iQi=βi，问题完全解耦。
由 K=H+∂S∂t=0K = H + \frac{\partial S}{\partial t} = 0K=H+∂t∂S=0 且 pi=∂S∂qip_i = \frac{\partial S}{\partial q_i}pi=∂qi∂S，立即得到：
H(q,∂S∂q,t)+∂S∂t=0\boxed{H\left(q, \frac{\partial S}{\partial q}, t\right) + \frac{\partial S}{\partial t} = 0}H(q,∂q∂S,t)+∂t∂S=0
这就是HJ方程。
核心思想 ：HJ方程的求解等价于寻找一个特殊的正则变换------把原来的Hamiltonian变为零的"终极简化"。一旦找到这个变换（即找到 SSS），整个力学问题就化为了平凡的常数运动。
三、深层含义
1. 力学 = 波传播的几何光学
HJ方程与几何光学 中的Eikonal方程（程函方程）形式完全相同：
(∂W∂x)2+(∂W∂y)2+(∂W∂z)2=n2(x,y,z)\left(\frac{\partial W}{\partial x}\right)^2 + \left(\frac{\partial W}{\partial y}\right)^2 + \left(\frac{\partial W}{\partial z}\right)^2 = n^2(x, y, z)(∂x∂W)2+(∂y∂W)2+(∂z∂W)2=n2(x,y,z)
其中 nnn 是折射率，WWW 是光程。
Hamilton本人在1830年代就发现了这个对应：
| 力学 | 几何光学 |
|---|---|
| Hamilton主函数 SSS | 光程函数（Eikonal） |
| 动量 ∇S=p\nabla S = p∇S=p | 光线方向 |
| Hamiltonian HHH | 折射率分布 |
| 等 SSS 面（波前） | 等相位面 |
| 质点轨迹 | 光线轨迹 |
这催生了Hamilton光学。后来de Broglie和Schrödinger正是沿着这条思路，把"波"的概念从光学推广到物质粒子，建立了量子力学。
2. 作用量作为"波前"
S(q,t)S(q, t)S(q,t) 的等值面 S=constS = \text{const}S=const 可以看作在构型空间中传播的波前。
动量 pi=∂S∂qip_i = \frac{\partial S}{\partial q_i}pi=∂qi∂S 表示波前的法线方向------即质点运动的方向。
整个力学演化可以看作一个波前在构型空间中的传播过程。
3. 完全可积性的判据
如果一个系统可以通过分离变量法求解HJ方程，它就是完全可积的。
Liouville-Arnold定理 说：如果系统有 nnn 个独立的、泊松对易的运动积分，则相空间被不变环面所纤维化，运动是准周期的------这等价于HJ方程可分离变量。
因此，HJ方程的可解性 = 系统的可积性。
4. 经典→量子的"通道"
这是HJ方程最深远的影响。
Schrödinger的推导思路（1926年）：
设波函数 Ψ=A(q,t)eiS(q,t)/ℏ\Psi = A(q, t)e^{iS(q, t)/\hbar}Ψ=A(q,t)eiS(q,t)/ℏ，代入Schrödinger方程，令 ℏ→0\hbar \to 0ℏ→0，领头阶就是HJ方程。
Schro¨dinger方程→ℏ→0HJ方程\text{Schrödinger方程} \xrightarrow{\hbar \to 0} \text{HJ方程}Schro¨dinger方程ℏ→0 HJ方程
HJ方程的 SSS 是量子力学相位 S/ℏS/\hbarS/ℏ 的经典极限。
| 量子力学 | 经典力学 |
|---|---|
| Schrödinger方程 | HJ方程 |
| 波函数 Ψ\PsiΨ | 作用量 SSS（相位） |
| $ | \Psi |
| 路径积分 ∫eiS/ℏDq\int e^{iS/\hbar}\mathcal{D}q∫eiS/ℏDq | 经典路径 δS=0\delta S = 0δS=0（驻相近似） |
四、求解HJ方程的方法
1. 分离变量法
如果HJ方程可以写成：
S=∑iSi(qi)−EtS = \sum_i S_i(q_i) - EtS=i∑Si(qi)−Et
且每个 SiS_iSi 只依赖于一个坐标，则方程可分离。例如中心力场问题：
- 径向坐标 rrr 和角坐标 θ,φ\theta, \varphiθ,φ 可分离
- S=Sr(r)+Sθ(θ)+Sφ(φ)−EtS = S_r(r) + S_\theta(\theta) + S_\varphi(\varphi) - EtS=Sr(r)+Sθ(θ)+Sφ(φ)−Et
2. Hamilton特征函数法（自治系统）
对于 HHH 不显含 ttt：
S(q,α,t)=W(q,α)−EtS(q, \alpha, t) = W(q, \alpha) - EtS(q,α,t)=W(q,α)−Et
其中 αi\alpha_iαi 是 nnn 个积分常数（新的"动量"）。
3. Jacobi定理------从 SSS 提取运动
一旦求出 S(q,α,t)S(q, \alpha, t)S(q,α,t)，运动的解由以下隐式方程给出：
βi=∂S∂αi,pi=∂S∂qi\beta_i = \frac{\partial S}{\partial \alpha_i}, \quad p_i = \frac{\partial S}{\partial q_i}βi=∂αi∂S,pi=∂qi∂S
其中 αi,βi\alpha_i, \beta_iαi,βi 是 2n2n2n 个由初始条件决定的常数。
五、在分析力学体系中的位置
Hamilton原理
δ∫L dt = 0
↓
Euler-Lagrange方程
(n个二阶ODE)
↓
Legendre变换
↓
Hamilton正则方程
(2n个一阶ODE)
↓
正则变换理论
(生成函数F₁, F₂, F₃, F₄)
↓
"终极正则变换"：令 K = 0
↓
Hamilton-Jacobi方程
(一个一阶PDE)
↓
┌─────────────────────────┐
↓ ↓
分离变量法 几何光学类比
(可积系统) (Eikonal方程)
↓ ↓
作用量-角变量 Schrödinger方程
(完全可积) (量子力学)
↓ ↓
KAM定理 Feynman路径积分六、一句话总结
Hamilton-Jacobi方程是从经典力学内部自然演化出的"终极方程"------它通过寻找一个把Hamiltonian变为零的正则变换，将 2n2n2n 个ODE的求解问题转化为一个PDE的求解问题。表面上看这是"升维"，但由于PDE的结构（特别是可分离变量性），问题反而可能变得可解。更深远的意义在于，HJ方程揭示了经典力学本质上是构型空间中的波传播问题，这正是从几何光学到物质波、从经典力学到量子力学的"黄金通道"。