3. 拉格朗日方程：广义坐标与自由度、拉格朗日函数、拉格朗日方程推导、保守系统与非保守系统

各位好，今天我们聊一个动力学建模里非常核心的工具——拉格朗日方程。说实话，我刚入行那会儿，面对复杂的多体系统，牛顿力学那一套写受力分析写得我头大。后来接触到拉格朗日方法，感觉像是打开了新世界的大门。它不关心系统内部那些乱七八糟的约束力，只盯着能量看，清爽多了。

3.1 广义坐标与自由度

先说说自由度。你想想看，一个自由质点在三维空间里乱飞，需要几个坐标来描述？三个，对吧？x、y、z。但如果这个质点被限制在一个球面上运动呢？那就只需要两个坐标，比如经纬度。这个“最少需要的独立坐标数目”，就是自由度。

那广义坐标是什么？说白了，就是用来描述系统位形的任意一组独立参数。不一定非得是笛卡尔坐标，可以是角度、距离，甚至是你自己定义的某个量。我习惯用 q₁, q₂, ..., qₙ 来表示，n 就是系统的自由度。

举个例子。一个单摆，摆长 l，摆角 θ。自由度是 1。广义坐标我选 θ，而不是 x 和 y。为什么？因为约束条件 x² + y² = l² 自动被广义坐标满足了，不用再费心去处理它。

3.2 拉格朗日函数

拉格朗日函数 L，定义很简单：

L = T - V

其中 T 是系统的动能，V 是势能。就这么简单。但你别小看它，整个系统的动力学信息都浓缩在这个标量函数里了。

为什么是 T - V 而不是 T + V？嗯，这个问题当年我也困惑过。其实从变分原理出发，自然就会得到这个形式。你暂时可以把它当作一个定义来接受，后面推导方程时你会看到它的妙处。

写 L 的时候，要注意几点：

• T 是广义速度 q̇ᵢ 的函数，也可能显含广义坐标 qᵢ（比如旋转动能与位置有关）
• V 通常是广义坐标 qᵢ 的函数，不显含速度（保守力情况下）
• L 也可能显含时间 t，如果系统有随时间变化的约束或外力

3.3 拉格朗日方程推导

好，重头戏来了。拉格朗日方程长这样：

d/dt (∂L/∂q̇ᵢ) - ∂L/∂qᵢ = Qᵢ

其中 Qᵢ 是广义力。如果系统是保守的，Qᵢ = 0，方程简化为：

d/dt (∂L/∂q̇ᵢ) - ∂L/∂qᵢ = 0

这个方程怎么来的？从哈密顿原理（也叫最小作用量原理）推导而来。简单说，系统真实的运动路径，使得作用量 S = ∫L dt 取极值。通过变分法，就能得到上面的方程。

我不打算在这里做严格的数学推导，那会占用太多时间。但我想强调一下这个方程的结构：

• ∂L/∂q̇ᵢ 是广义动量 pᵢ
• d/dt (∂L/∂q̇ᵢ) 是广义动量的变化率
• ∂L/∂qᵢ 是广义力（保守部分）

所以拉格朗日方程本质上就是广义的牛顿第二定律：动量变化率 = 力。

3.4 保守系统与非保守系统

保守系统，就是机械能守恒的系统。说白了，没有摩擦，没有阻尼，没有非保守外力。这时候拉格朗日方程右边是零，形式最简洁。

非保守系统呢？有耗散力（比如阻尼）、有外力（比如电机驱动力）。这时候就需要在方程右边加上广义力 Qᵢ。

广义力怎么算？虚功法：

δW = Σ Qᵢ δqᵢ

其中 δW 是所有非保守力在虚位移上做的虚功。把虚功表达成广义坐标变分的形式，系数就是广义力。

举个例子，一个带阻尼的单摆：

• 动能 T = ½ m l² θ̇²
• 势能 V = m g l (1 - cos θ)
• 阻尼力 F_d = -c l θ̇（与速度相反）

拉格朗日函数 L = T - V。广义力 Q 来自阻尼力做的虚功：δW = -c l² θ̇ δθ，所以 Q = -c l² θ̇。

代入拉格朗日方程：

d/dt (m l² θ̇) - (-m g l sin θ) = -c l² θ̇
→ m l² θ̈ + m g l sin θ = -c l² θ̇
→ m l² θ̈ + c l² θ̇ + m g l sin θ = 0

你看，阻尼项自然就出现在方程里了。

3.5 知识体系总览

下面这张图，是我自己总结的拉格朗日方程知识框架，帮你理清思路：

3.6 实用建议

最后，给各位几个实战中的小建议：

1. 先判断自由度：拿到一个系统，第一件事就是数自由度。数错了，后面全白搭。
2. 选好广义坐标：选得好，方程简洁；选得差，自找麻烦。多花点时间在这步上，值得。
3. 写 L 要仔细：动能和势能表达式写错了，方程肯定不对。我建议写完 L 后，检查一下量纲是否一致。
4. 非保守力单独处理：先把保守系统的拉格朗日方程写出来，再通过虚功加广义力。分步走，不容易出错。
5. 验证结果：对于简单系统（比如单摆、弹簧振子），拉格朗日方程应该退化成你熟悉的牛顿方程。如果对不上，说明哪里出问题了。

---