分析力学基本原理介绍7.6：用变分原理推导哈密顿方程

$\bullet$ 变分原理/哈密顿原理既可以推导出拉格朗日方程，也可用来推导哈密顿正则方程。为了推出后者，我们先要对原始的变分原理做出一些调整。

原始的哈密顿原理指出，系统偏离正确路径的虚位移导致的作用量的虚变分必须为零：

$\delta I \equiv \delta\int_{t_1}^{t_2}\mathscr{L} \;dt = 0\\$

为了推出哈密顿方程，变量 $q$ 和 $p$ 都必须为独立变量。所以作用量的积分路径需从由 $n$ 个广义坐标张成的位形空间转移至由 $2n$ 个广义坐标以及正则动量张成的相空间。借助勒让德变换，我们得到：

$\boxed{\delta I = \delta\int _{t_1}^{t_2}(\sum_ip_i\dot{q_i} - \mathscr{H}(\mathbf{q},\mathbf{p},t))dt = 0}\\$

这是相空间的哈密顿原理，因为是由原始的哈密顿原理得来，所以也被称作修改的哈密顿原理(modified Hamilton's principle)。该理论也常常出现在与变换理论(transformation theory)相关的内容中。

$\bullet$ 修改的哈密顿原理与高维变分原理其实并无太大差别。在 $2n$ 维相空间中的变分原理具有形式：

$\delta I = \delta\int_{t_1}^{t_2}f(\mathbf{q},\mathbf{\dot{q}},\mathbf{p},\mathbf{\dot{p}},t)\;dt = 0\\$

其中 $f(\mathbf{q},\mathbf{\dot{q}},\mathbf{p},\mathbf{\dot{p}},t)= \sum_ip_i\dot{q_i} - \mathscr{H}(\mathbf{q},\mathbf{p},t)$

于是

$\mathbf{\delta} I = \sum_j \int_{t_1}^{t_2}\left(\frac{\partial f}{\partial q_j}\delta q_j + \frac{\partial f}{\partial \dot{q_j}}\delta\dot{q_j}\right)\;dt + \sum_j \int_{t_1}^{t_2}\left(\frac{\partial f}{\partial p_i}\delta p_i + \frac{\partial f}{\partial \dot{p_i}}\delta\dot{p_i}\right)\;dt\\$

$\implies \left.\frac{\partial I}{\partial \alpha}\right|_{\alpha = 0} = \sum_j \int_{t_1}^{t_2}\left(\frac{\partial f}{\partial q_j}\frac{\partial q_i}{\partial \alpha} + \frac{\partial f}{\partial \dot{q_j}}\frac{\partial\dot{q_i}}{\partial\alpha}\right)\;dt + \sum_j \int_{t_1}^{t_2}\frac{\partial}{\partial p_j}(\sum_ip_i\dot{q_i} - \mathscr{H}(\mathbf{q},\mathbf{p},t))\frac{\partial p_j}{\partial\alpha}\;dt$

$\implies \sum_j\int_{t_1}^{t_2}\left[\frac{\partial f}{\partial q_i} - \frac{d}{dt}\left(\frac{\partial f}{\partial \dot{q_i}}\right)\right]\frac{\partial q_i}{\partial\alpha}\;dt + \sum_j \int_{t_1}^{t_2}\left(\dot{q_j} - \frac{\partial \mathscr{H}}{\partial p_j}\right)\frac{\partial p_j}{\partial\alpha}\;dt = 0$

根据变分法基本引理，得到 $2n$ 个欧拉-拉格朗日方程：

$\frac{\partial f}{\partial q_i} - \frac{d}{dt}\left(\frac{\partial f}{\partial \dot{q_i}}\right) = 0;\quad \dot{q_j} - \frac{\partial \mathscr{H}}{\partial p_j}=0\\$

将 $f$ 代入第一个等式

$\frac{\partial f}{\partial \dot{q_j}} = \sum_i p_i \frac{\partial \dot{q_i}}{\partial \dot{q_j}} = p_i \delta_ j^i = p_j\\$

$\frac{\partial f}{\partial q_j} = -\frac{\partial\mathscr{H}}{\partial q_j}\\$

$\implies \boxed{\dot{p_j} + \frac{\partial\mathscr{H}}{\partial q_j} = 0}\\$

结合第二个等式，

$\boxed{\dot{q_j} - \frac{\partial \mathscr{H}}{\partial p_j}=0}\\$

我们得到了哈密顿正则方程。

$\bullet$ 就上述步骤而言，推导看上去似乎过于简单了。很自然地，你可能会问：会不会有什么隐含的物理条件？修改的哈密顿原理是否真的等价于原始哈密顿原理？修改的哈密顿原理只保证了在相空间中能够导出正确的正则方程，并非用于位形空间。而且，即便对于原始的哈密顿原理，也没有再出现任何额外条件了。所以，这些疑问很大程度上都是不相关的：只要勒让德变换能保证拉格朗日函数和哈密顿函数具有相同的物理意义，使用修改后的哈密顿原理构造哈密顿方程就是有意义的。

$\bullet$ 有关路径作用量边界的问题反倒更值得注意。对于位形空间里的路径而言，推导出拉格朗日方程只需要路径的两个端点保持固定，即 $\delta q_{i}(t_1) = \delta q_{i}(t_2) = 0$ 。但对于修改后的哈密顿原理，在相空间里的路径虚变分为零，我们不仅需要 $\delta q_{i} = 0$ ，还需要 $\delta p_{i} = 0$ 。条件的差异是否会影响推导的有效性？

答案是不会，在修改的哈密顿原理下，作用量的被积函数 $f$ 其实并不存在对变量 $\dot{p_i}$ 的依赖关系， $\frac{\partial f}{\partial \dot{p_i}} = 0$ ，所以在推导第二组哈密顿正则方程时，我们根本不需要使用条件 $\delta p_{i} = 0$ 或 $\delta \dot{p_i} = 0$ 。

$\bullet$ 当然，边界条件的引入也存在优势，它可以将哈密顿方程的成立条件推广至更一般的情况：

当相空间中路径的端点固定时，被积函数即使加上或减去一个连续的二次可微函数 $F(\mathbf{q},\mathbf{p},t)$ 关于时间的全导项也不会影响最后的哈密顿方程的形式。

上述结论很容易证明。

设函数 $F(\mathbf{q},\mathbf{p},t) = \sum_iq_i p_i$ ，使用修改的哈密顿原理：

$\begin{align*} \delta I &= \delta\int_{t_1}^{t_2}\left[\sum_i p_i\dot{q_i} - \mathscr{H}(\mathbf{q},\mathbf{p},t) - \frac{d}{dt}(\sum_i q_i p_i)\right]\;dt\\&= \delta\int_{t_1}^{t_2}\left[-\mathscr{H}(\mathbf{q},\mathbf{p},t)-\sum_iq_i\dot{p_i}\right]\;dt\\\end{align*}\\$

我们将此时的被积函数用 $g$ 表示： $g(\mathbf{q},\mathbf{\dot{q}},\mathbf{p},\mathbf{\dot{p}},t) = -\mathscr{H}(\mathbf{q},\mathbf{p},t)- \sum_iq_i\dot{p_i}$ ，于是

$\delta \int_{t_1}^{t_2}g\;dt = \sum_j \int_{t_1}^{t_2}\left(\frac{\partial g}{\partial q_j}\delta q_j + \frac{\partial g}{\partial \dot{q_j}}\delta\dot{q_j}\right)\;dt + \sum_j \int_{t_1}^{t_2}\left(\frac{\partial g}{\partial p_i}\delta p_i + \frac{\partial g}{\partial \dot{p_i}}\delta\dot{p_i}\right)\;dt \\$

$\implies \left.\frac{\partial I}{\partial\alpha}\right|_{\alpha = 0} = \sum_j\int_{t_1}^{t_2}\left(\frac{\partial g}{\partial q_j}\frac{\partial q_j}{\partial \alpha} + \frac{\partial g}{\partial \dot{q_j}}\frac{\partial^2 q_j}{\partial t\partial\alpha}\right)\;dt + \sum_j\int_{t_1}^{t_2}\left(\frac{\partial g}{\partial p_j}\frac{\partial p_j}{\partial\alpha}+\frac{\partial g}{\partial \dot{p_j}}\frac{\partial^2p_j}{\partial t\partial \alpha}\right)\;dt$

$\implies \sum_j\int_{t_1}^{t_2}\left[\frac{\partial g}{\partial q_j} - \frac{d}{dt}\left(\frac{\partial g}{\partial \dot{q_j}}\right)\right]\frac{\partial q_j}{\partial \alpha}\;dt + \sum_j\int_{t_1}^{t_2}\left[\frac{\partial g}{\partial p_j} - \frac{d}{dt}\left(\frac{\partial g}{\partial \dot{p_j}}\right)\right]\frac{\partial p_j}{\partial \alpha}\;dt$

$\implies \begin{cases}\frac{\partial g}{\partial q_j} - \frac{d}{dt}\left(\frac{\partial g}{\partial \dot{q_j}}\right) = -\frac{\partial \mathscr{H}}{\partial q_j} - \sum_i \dot{p_i}\delta_j^i = \boxed{-\frac{\partial\mathscr{H}}{\partial q_j} - \dot{p_j} = 0}\\ \frac{\partial g}{\partial p_j} - \frac{d}{dt}\left(\frac{\partial g}{\partial \dot{p_j}}\right) = -\frac{\partial\mathscr{H}}{\partial p_j} - \frac{d}{dt}\left(-\sum_i q_i\delta_j^i\right) = \boxed{-\frac{\partial\mathscr{H}}{\partial p_j} + \dot{q_j} = 0}\end{cases} \\$

$\bullet$ 这时的函数 $g$ 既不是拉格朗日函数，也不属于任何可以通过位形空间的拉格朗日函数点变换得到函数。可见，边界条件的引入去除了拉格朗日函数广义坐标、速度与哈密顿函数正则变量之间的关系。借助边界条件，即使不从拉格朗日函数出发，我们也可以得到哈密顿方程。

最后编辑于：2020.02.26 12:37:12

人面猴
序言：七十年代末，一起剥皮案震惊了整个滨河市，随后出现的几起案子，更是在滨河造成了极大的恐慌，老刑警刘岩，带你破解...
沈念sama阅读 219,753评论 6赞 508
死咒
序言：滨河连续发生了三起死亡事件，死亡现场离奇诡异，居然都是意外死亡，警方通过查阅死者的电脑和手机，发现死者居然都...
沈念sama阅读 93,668评论 3赞 396
救了他两次的神仙让他今天三更去死
文/潘晓璐我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
开封第一讲书人阅读 166,090评论 0赞 356
道士缉凶录：失踪的卖姜人
文/不坏的土叔我叫张陵，是天一观的道长。经常有香客问我，道长，这世上最难降的妖魔是什么？我笑而不...
开封第一讲书人阅读 59,010评论 1赞 295
港岛之恋（遗憾婚礼）
正文为了忘掉前任，我火速办了婚礼，结果婚礼上，老公的妹妹穿的比我还像新娘。我一直安慰自己，他们只是感情好，可当我...
茶点故事阅读 68,054评论 6赞 395
恶毒庶女顶嫁案：这布局不是一般人想出来的
文/花漫我一把揭开白布。她就那样静静地躺着，像睡着了一般。火红的嫁衣衬着肌肤如雪。梳的纹丝不乱的头发上，一...
开封第一讲书人阅读 51,806评论 1赞 308
城市分裂传说
那天，我揣着相机与录音，去河边找鬼。笑死，一个胖子当着我的面吹牛，可吹牛的内容都是我干的。我是一名探鬼主播，决...
沈念sama阅读 40,484评论 3赞 420
双鸳鸯连环套：你想象不到人心有多黑
文/苍兰香墨我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
开封第一讲书人阅读 39,380评论 0赞 276
万荣杀人案实录
序言：老挝万荣一对情侣失踪，失踪者是张志新（化名）和其女友刘颖，没想到半个月后，有当地人在树林里发现了一具尸体，经...
沈念sama阅读 45,873评论 1赞 319
护林员之死
正文独居荒郊野岭守林人离奇死亡，尸身上长有42处带血的脓包…… 初始之章·张勋以下内容为张勋视角年9月15日...
茶点故事阅读 38,021评论 3赞 338
白月光启示录
正文我和宋清朗相恋三年，在试婚纱的时候发现自己被绿了。大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
茶点故事阅读 40,158评论 1赞 352
活死人
序言：一个原本活蹦乱跳的男人离奇死亡，死状恐怖，灵堂内的尸体忽然破棺而出，到底是诈尸还是另有隐情，我是刑警宁泽，带...
沈念sama阅读 35,838评论 5赞 346
日本核电站爆炸内幕
正文年R本政府宣布，位于F岛的核电站，受9级特大地震影响，放射性物质发生泄漏。R本人自食恶果不足惜，却给世界环境...
茶点故事阅读 41,499评论 3赞 331
男人毒药：我在死后第九天来索命
文/蒙蒙一、第九天我趴在偏房一处隐蔽的房顶上张望。院中可真热闹，春花似锦、人声如沸。这庄子的主人今日做“春日...
开封第一讲书人阅读 32,044评论 0赞 22
一桩弑父案，背后竟有这般阴谋
文/苍兰香墨我抬头看了看天上的太阳。三九已至，却和暖如春，着一层夹袄步出监牢的瞬间，已是汗流浃背。一阵脚步声响...
开封第一讲书人阅读 33,159评论 1赞 272
情欲美人皮
我被黑心中介骗来泰国打工，没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留，地道东北人。一个月前我还...
沈念sama阅读 48,449评论 3赞 374
代替公主和亲
正文我出身青楼，却偏偏与公主长得像，于是被迫代替她去往敌国和亲。传闻我的和亲对象是个残疾皇子，可洞房花烛夜当晚...
茶点故事阅读 45,136评论 2赞 356

分析力学基本原理介绍7.6：用变分原理推导哈密顿方程

推荐阅读更多精彩内容