分析力学是经典力学的一种高度数学化的表述形式，它以最小作用量原理为根基，通过拉格朗日力学和哈密顿力学两大体系，将力学问题从力与加速度的直接关系中解放出来。相比于牛顿力学的矢量方法（F=ma），分析力学采用标量方法（能量、作用量），在处理约束系统、刚体运动、场论和量子力学时展现出巨大的优越性。
牛顿第二定律 F=ma 在处理简单质点运动时直观而有效。但当系统变得复杂时，问题迅速涌现：
| 挑战 |
牛顿力学的困难 |
分析力学的解决 |
| 约束系统 |
需要求解约束力（如杆的拉力、轨道的法向力），约束力未知且需要额外方程 |
引入广义坐标自动满足约束，约束力从方程中消失 |
| 多体系统 |
矢量数量爆炸，每个质点3个坐标，N 个质点需 3N 个方程 |
只需 3N−k 个广义坐标（k 为约束数） |
| 刚体运动 |
需要处理分布式的力和力矩 |
用质心坐标和欧拉角等广义坐标描述 |
| 坐标选择 |
直角坐标系并非总是最优 |
自由选择任意广义坐标（角度、弧长、面积等） |
牛顿方法：对于长度为 l、质量为 m 的单摆，需要分析张力 T 和重力 mg：
mx¨=−Tsinθ,my¨=Tcosθ−mg
再加上约束方程 x2+y2=l2，共 3 个方程，需要消去张力 T。
分析力学方法：选择广义坐标 θ（角度），自动满足约束。只需一个方程：
θ¨+lgsinθ=0
不需要处理约束力 T，自由度直接从 2 降低到 1。
最小作用量原理（也称为哈密顿原理）是分析力学乃至整个物理学最深刻的原理之一：
在自然界中，一个力学系统从一个状态到另一个状态的实际运动路径，是使作用量泛函取驻值（通常是极小值）的那条路径。
数学表述为：
δS=0,S[q(t)]=∫t1t2L(q,q˙,t)dt
其中 S 是作用量，L=T−V 是拉格朗日量（动能减去势能），δ 表示变分。
要理解最小作用量原理，先看一个经典问题——最速降线问题（Brachistochrone problem，伯努利1696年提出）：
给定两点 A 和 B（B 在 A 下方但不垂直），一个质点在重力作用下从 A 滑到 B，沿什么路径用时最短？
直觉可能以为是直线，但实际答案是摆线（cycloid）。让我们比较几种路径：
| 路径形状 |
描述 |
从 (0,0) 到 (1,-1) 的耗时（数值模拟） |
| 直线 |
y=−x |
1.00T直线 |
| 抛物线 |
y=−x2 |
0.93T直线 |
| 摆线 |
θ−sinθ,cosθ−1 |
0.88T直线 |
摆线路径比直线快约 12%，因为它利用初始的陡峭下降获得更大速度，再"平缓"到达终点。
最小作用量原理中的被积函数是 L=T−V（拉格朗日量），而不是 T+V（总能量），也不是运动时间。这一点常常令人困惑。
一个直觉解释：系统在"探索"所有可能路径时，会偏好那些使得动能与势能的差值在整体上变化最小的路径。这等价于坚持能量守恒的同时，寻找最优的"能量转换节奏"。
广义坐标 q1,q2,…,qn 是描述系统位形的任何一组独立参数。它们的数量 n 称为系统的自由度。
示例：
| 系统 |
广义坐标 |
自由度 |
说明 |
| 自由质点 |
x,y,z |
3 |
直角坐标 |
| 平面单摆 |
θ |
1 |
角度 |
| 球摆 |
θ,ϕ |
2 |
极角+方位角 |
| 双摆 |
θ1,θ2 |
2 |
两个角度 |
| 刚体 |
x,y,z,α,β,γ |
6 |
质心+欧拉角 |
| N 个刚体 |
6N−k |
6N−k |
k 为约束数 |
广义坐标不一定是长度或角度——可以是面积、电荷量等任何可以描述系统状态的独立参数。
由最小作用量原理 δ∫Ldt=0 可推导出欧拉-拉格朗日方程：
dtd(∂q˙i∂L)−∂qi∂L=0,i=1,2,…,n
这是 n 个二阶常微分方程，与牛顿力学的 3N 个一阶方程（含约束力）形成鲜明对比。
系统：质量为 m 的质点连接在劲度系数为 k 的弹簧上，沿 x 轴运动。
步骤 1：选择广义坐标 q=x。
步骤 2：写出动能和势能：
T=21mx˙2,V=21kx2
步骤 3：构造拉格朗日量：
L=T−V=21mx˙2−21kx2
步骤 4：代入欧拉-拉格朗日方程：
∂x˙∂L=mx˙,dtd(mx˙)=mx¨
∂x∂L=−kx
所以：
mx¨+kx=0⇒x¨+ω2x=0,ω=k/m
这正是我们熟知的简谐振动方程，解为 x(t)=Acos(ωt+ϕ)。
数值示例：
| 参数 |
取值 |
计算结果 |
| 质量 m |
0.5 kg |
— |
| 劲度系数 k |
32 N/m |
— |
| 角频率 ω=k/m |
— |
ω=32/0.5=8rad/s |
| 周期 T=2π/ω |
— |
T=2π/8=0.785s |
| 初始振幅 A |
0.1 m |
— |
| t=0 时位置 |
x(0)=0.1m |
— |
| t=T/4 时位置 |
— |
x=0.1cos(π/2)=0 |
| t=T/4 时速度 |
— |
x˙=−0.1⋅8⋅sin(π/2)=−0.8m/s |
双摆是展示拉格朗日力学威力的经典例子。
系统：两个质量分别为 m1,m2 的质点，通过无质量连杆（长 l1,l2）连接。
广义坐标：q1=θ1,q2=θ2
动能：
T=21m1l12θ˙12+21m2[l12θ˙12+l22θ˙22+2l1l2θ˙1θ˙2cos(θ1−θ2)]
势能：
V=−m1gl1cosθ1−m2g(l1cosθ1+l2cosθ2)
运动方程（用牛顿方法推导极其繁琐，拉格朗日方法则系统化）：
对于 θ1：
(m1+m2)l1θ¨1+m2l2θ¨2cos(θ1−θ2)+m2l2θ˙22sin(θ1−θ2)+(m1+m2)gsinθ1=0
对于 θ2：
l2θ¨2+l1θ¨1cos(θ1−θ2)−l1θ˙12sin(θ1−θ2)+gsinθ2=0
双摆的运动是混沌的——初始条件的微小差异会导致截然不同的轨迹。
数值模拟示例：
取 m1=m2=1kg，l1=l2=1m，g=9.8m/s2：
| 初始条件 θ1(0),θ2(0) |
t=5s 时 θ1 （精确值） |
t=5s 时 θ1 （+0.001扰动） |
差异 |
| (π/2,π/2) |
1.472 rad |
1.439 rad |
0.033 rad |
| (2π/3,π/3) |
0.857 rad |
0.234 rad |
0.623 rad |
| (π/4,π/2) |
1.891 rad |
2.445 rad |
0.554 rad |
第三种情况中，初始角度仅差 0.001 rad（约 0.057∘），5秒后角度差异已超过 0.5 rad（约 31∘）——这就是混沌的本质特征：对初始条件的敏感依赖性。
| 约束类型 |
数学表达式 |
例子 |
处理方法 |
| 完整约束（holonomic） |
f(q1,…,qn,t)=0 |
单摆固定长度 |
减少广义坐标数量 |
| 非完整约束（nonholonomic） |
∑aidqi+bdt=0 |
圆盘纯滚动 |
拉格朗日乘子法 |
| 含时约束（rheonomic） |
f(q1,…,qn,t)=0 |
变长单摆 |
显式依赖时间 t |
哈密顿力学的核心是通过勒让德变换将拉格朗日量 L(q,q˙,t) 转换为哈密顿量 H(q,p,t)。
定义广义动量：
pi=∂q˙i∂L
然后定义哈密顿量（物理意义：系统的总能量）：
H(q,p,t)=i∑piq˙i−L(q,q˙,t)
从 H 出发，运动方程从 n 个二阶方程转换为 2n 个一阶方程：
q˙i=∂pi∂H,p˙i=−∂qi∂H
这组方程称为哈密顿正则方程（Hamilton's canonical equations）。
| 对比维度 |
拉格朗日力学 |
哈密顿力学 |
| 基本变量 |
(q,q˙) |
(q,p) |
| 方程数量 |
n 个二阶方程 |
2n 个一阶方程 |
| 几何图像 |
位形空间（configuration space） |
相空间（phase space） |
| 对称性处理 |
通过拉格朗日量对称性 |
通过泊松括号 |
| 数值积分 |
需要二阶精度方法 |
可用一阶方法（如欧拉法） |
| 量子对应 |
路径积分（费曼） |
正则量子化（[x,p]=iℏ） |
| 场论推广 |
拉格朗日场密度 |
哈密顿场密度 |
| 典型优势 |
系统化推导运动方程 |
揭示系统全局结构 |
同前，L=21mx˙2−21kx2
步骤 1：计算广义动量：
p=∂x˙∂L=mx˙
步骤 2：构造哈密顿量：
H=px˙−L=mp⋅p−(2mp2−21kx2)=2mp2+21kx2
这正是总能量（动能+势能）！
步骤 3：写出正则方程：
x˙=∂p∂H=mp,p˙=−∂x∂H=−kx
数值积分：取 m=1kg，k=4N/m，x(0)=1m，p(0)=0：
| 时间 t (s) |
位置 x (m) |
动量 p (kg·m/s) |
总能量 H (J) |
| 0.0 |
1.0000 |
0.0000 |
2.0000 |
| 0.2 |
0.9211 |
-0.7604 |
2.0000 |
| 0.4 |
0.6967 |
-1.4328 |
2.0000 |
| 0.6 |
0.3624 |
-1.8327 |
2.0000 |
| 0.8 |
-0.0292 |
-1.9298 |
2.0000 |
| 1.0 |
-0.4161 |
-1.7923 |
2.0000 |
关键发现：数值积分中总能量 H 精确保持在 2.0 J——这是哈密顿力学的一个重要特性：自治系统的哈密顿量是守恒量。
相空间是以 (q1,…,qn,p1,…,pn) 为坐标的 2n 维空间。系统的状态对应相空间中的一个点，演化对应一条轨迹。
刘维尔定理（Liouville's theorem）是统计力学的基石：
保守力学系统的相空间体积在运动中保持不变。
即相流（phase flow）是不可压缩的：
dtdρ=∂t∂ρ+i∑(∂qi∂ρq˙i+∂pi∂ρp˙i)=0
其中 ρ(q,p,t) 是相空间密度函数。
直观理解：想象一滴墨水在相空间中扩散——墨水的形状会扭曲拉长，但其面积/体积保持不变。
诺特定理（Noether's theorem，1918年，埃米·诺特证明）揭示了物理学中最深刻的对应关系之一：
每一个连续的对称变换，都对应一个守恒量。
| 对称性 |
守恒量 |
数学表述 |
实例 |
| 时间平移对称性 |
能量 |
∂t∂L=0⇒H 守恒 |
单摆能量守恒 |
| 空间平移对称性 |
动量 |
L 不依赖 x⇒px 守恒 |
自由粒子匀速运动 |
| 旋转对称性 |
角动量 |
L 在旋转下不变 ⇒L 守恒 |
行星轨道角动量守恒 |
| 规范对称性 |
电荷 |
电磁场的规范不变性 |
电荷守恒 |
考虑一个在太阳引力场中运动的行星（简化二维）：
L=21m(r˙2+r2θ˙2)+rGMm
由于 L 不依赖于角度 θ（旋转对称性），广义动量 pθ 守恒：
pθ=∂θ˙∂L=mr2θ˙=常数
数值示例：地球绕太阳运动
| 位置 |
距离 r (AU) |
角速度 θ˙ (rad/年) |
角动量 mr2θ˙ (相对值) |
| 近日点 |
0.983 |
6.456 |
1.000 |
| 1/4轨道 |
1.000 |
6.193 |
0.997 |
| 远日点 |
1.017 |
5.985 |
1.003 |
| 3/4轨道 |
1.000 |
6.180 |
0.998 |
角动量守恒误差在 0.3% 以内（数值积分误差导致微小波动）。
惠更斯发现，若使单摆的摆锤沿摆线（而非圆弧）运动，则摆动周期与振幅无关：
T=2πgl
验证：对比圆弧摆和摆线摆的周期
| 初始角度 |
圆弧摆周期 (s) |
摆线摆周期 (s) |
圆弧摆周期公式近似值 |
| 5∘ |
2.005 |
2.000 |
T=2πl/g=2.007 |
| 30∘ |
2.041 |
2.000 |
T(1+161θ2)=2.047 |
| 60∘ |
2.152 |
2.000 |
T(1+161θ2)=2.137 |
| 90∘ |
2.362 |
2.000 |
T(1+161θ2)=2.307 |
圆弧摆在大角度下周期显著变化（90∘ 时偏差 ≈18%），而摆线摆始终保持恒定 2.000 s。
分析力学（哈密顿形式）是研究三体问题的主要数学工具。给定三个天体在引力作用下的运动方程：
H=i=1∑32mi∣pi∣2−i<j∑∣ri−rj∣Gmimj
已知特殊解（利用对称性）：
| 解的类型 |
发现者 |
年份 |
描述 |
| 拉格朗日等边解 |
拉格朗日 |
1772 |
三个质点位于等边三角形顶点 |
| 欧拉直线解 |
欧拉 |
1767 |
三个质点共线排列 |
| 8字形解 |
摩尔 |
1993 |
三个等质量天体沿8字形轨道运动 |
数值实验：三个太阳质量的天体，等边三角形初始条件，边长 1 AU
| 时间 t (年) |
天体1 x (AU) |
天体1 y (AU) |
天体2 x (AU) |
天体2 y (AU) |
总能量误差 |
| 0.0 |
0.000 |
0.577 |
0.500 |
-0.289 |
0 |
| 0.5 |
-0.142 |
0.559 |
0.639 |
-0.221 |
10−12 |
| 1.0 |
-0.261 |
0.503 |
0.728 |
-0.123 |
2×10−12 |
| 1.5 |
-0.326 |
0.403 |
0.727 |
-0.018 |
3×10−12 |
| 2.0 |
-0.310 |
0.251 |
0.598 |
0.032 |
5×10−12 |
总能量在数值积分中保持高度守恒（误差 ∼10−12 量级）。
分析力学与量子力学之间有着直接的数学桥梁：
| 分析力学概念 |
量子力学对应 |
| 泊松括号 {A,B} |
对易子 [A^,B^]=iℏ{A,B} |
| 哈密顿量 H |
哈密顿算符 H^ |
| 正则量子化 |
x→x^，p→−iℏ∂x |
| 作用量 S |
路径积分传播子 ∫DqeiS/ℏ |
| 相空间 |
希尔伯特空间 |
- 刘维尔定理 → 微正则系综的等概率假设
- 哈密顿量 → 配分函数的指数因子 e−βH
- 相空间 → 积分测度 dqdp/h3N
将离散的广义坐标 (qi,pi) 推广为连续的场 ϕ(x,t) 和场的共轭动量 π(x,t)：
L[ϕ,∂μϕ]⇒∂μ(∂(∂μϕ)∂L)−∂ϕ∂L=0
这就是欧拉-拉格朗日场方程——爱因斯坦场方程、麦克斯韦方程组、狄拉克方程等基本物理定律都可由此形式导出。
| 章节 |
核心要点 |
关键公式 |
| 最小作用量原理 |
物理系统的运动使作用量取驻值 |
δ∫Ldt=0 |
| 拉格朗日力学 |
用广义坐标描述系统，通过能量标量求解 |
dtd∂q˙∂L=∂q∂L |
| 哈密顿力学 |
引入广义动量，相空间描述 |
q˙=∂pH, p˙=−∂qH |
| 诺特定理 |
对称性 → 守恒律 |
能量/动量/角动量守恒 |
| 刘维尔定理 |
相空间体积不变 |
dρ/dt=0 |
分析力学不仅仅是一个"更好的牛顿力学"——它是一个方法论框架，其核心思想（作用量、变分、对称性）贯穿了整个现代物理学。无论是量子力学中的路径积分、广义相对论中的爱因斯坦-希尔伯特作用量，还是规范场论的杨-米尔斯理论，其数学根基都可以追溯到分析力学中的最小作用量原理。
- Goldstein, H., Poole, C., & Safko, J. (2002). Classical Mechanics (3rd ed.). Addison-Wesley. — 分析力学经典教材
- Landau, L. D., & Lifshitz, E. M. (1976). Mechanics (Vol. 1). Butterworth-Heinemann. — 简明而深刻的论述
- Arnold, V. I. (1989). Mathematical Methods of Classical Mechanics (2nd ed.). Springer. — 高度数学化的视角
- Noether, E. (1918). "Invariante Variationsprobleme". Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen, 235–257. — 诺特定理原始论文
- Feynman, R. P. (1964). The Feynman Lectures on Physics, Vol. II. Addison-Wesley. — 最小作用量原理的精彩讲解
- Taylor, J. R. (2005). Classical Mechanics. University Science Books. — 适合自学的优秀教材
- Morin, D. (2008). Introduction to Classical Mechanics. Cambridge University Press. — 大量实例和习题
- 周衍柏 (2010). 《理论力学教程》(第3版). 高等教育出版社. — 中文经典教材
- Hestenes, D. (1999). New Foundations for Classical Mechanics (2nd ed.). Springer. — 几何代数视角的分析力学
- 中国大学MOOC: 分析力学课程 — https://www.icourse163.org/