分析力学基本原理介绍7.7:最小作用量原理

\bullet现在我们讨论另一个与哈密顿体系有关的变分原理。它涉及另一种变分——\Delta-变分。用于推导拉格朗日方程的哈密顿原理所用的是\delta-变分。对于\delta-变分,我们会固定路径两个端点,即,\delta q_i(t_1) = \delta q_i(t_2) = 0

相比之下,\Delta-变分在使用条件上会宽松一些。通常情况下,对于\Delta-变分,错误路径两个端点的时间间隔跟正确路径相比会有所不同,而且端点所对应的坐标也可能改变。

\bullet对于\Delta-变分,我们仍然可以使用与\delta-变分相同的参化方程:

q_i(t,\alpha) = q_i(t,0) + \alpha \eta_i(t)

它代表一个路径曲线家族。

\bullet作用量的\Delta-变分被定义为:

\Delta \int_{t_1}^{t_2}\mathscr{L}\;dt \equiv \int_{t_1 + \Delta t_1}^{t_2 + \Delta t_2}\mathscr{L}(\alpha)\;dt - \int_{t_1}^{t_2}\mathscr{L}(0)\;dt

其中\mathscr{L}(\alpha)代表剩下无数条错误路径;\mathscr{L}(0)代表正确路径。

\bullet作用量的\Delta-变分主要由两部分造成:

1)积分边界的变分(路径端点的变分)

2)被积函数的变分(偏离正确路径)

定义中的第一项近似到一阶小量就等于被积函数乘以边界的变分;第二项则与之前的\delta-变分类似:

\Delta \int_{t_1}^{t_2}\mathscr{L}\;dt = \mathscr{L}(t_2)\Delta t_2 - \mathscr{L}(t_1)\Delta t_1 + \delta\int_{t_1}^{t_2}\mathscr{L}\;dt

\bullet让我们先考虑第二项。变化的路径仍可以使用跟之前相同的参化方程,不同的是,现在路径端点处q_i\delta-变分并不等于零,所以我们要保留由分部积分法得到的端点项:

\begin{align*}\mathbf{\delta} \int_{t_1}^{t_2}\mathscr{L}\;dt &= \sum_i\int_{t_1}^{t_2}\frac{\partial\mathscr{L}}{\partial q_i}\delta q_i + \frac{\partial\mathscr{L}}{\partial \dot{q_i}}\delta\dot{q_i}\\&= \sum_i \int_{t_1}^{t_2}\frac{\partial \mathscr{L}}{\partial q_i}\delta q_idt + \sum_i\int_{t_1}^{t_2}\frac{\partial \mathscr{L}}{\partial \dot{q_i}}\delta\dot{q_i}dt\\&= \sum_i\int_{t_1}^{t_2}\frac{\partial \mathscr{L}}{\partial q_i}\delta q_idt + \sum_i \left.\frac{\partial\mathscr{L}}{\partial\dot{q_i}}\delta q_i\right|_{t_1}^{t_2} - \sum_i\int_{t_1}^{t_2}\frac{d}{dt}\left(\frac{\partial\mathscr{L}}{\partial \dot{q_i}}\right)\delta q_idt\\&= \sum_i\int_{t_1}^{t_2}\left[\frac{\partial \mathscr{L}}{\partial q_i} - \frac{d}{dt}\left(\frac{\partial\mathscr{L}}{\partial \dot{q_i}}\right)\right]\delta q_idt + \sum_i \left.\frac{\partial \mathscr{L}}{\partial \dot{q_i}}\delta q_i\right|_{t_1}^{t_2}\\&= \sum_i \left.\frac{\partial \mathscr{L}}{\partial \dot{q_i}}\delta q_i\right|_{t_1}^{t_2}\end{align*}

根据拉格朗日方程,带方括号的第一项消失。

代入作用量变分中,得

\begin{align*}\Delta \int_{t_1}^{t_2}\mathscr{L}\;dt &= \mathscr{L}(t_2)\Delta t_2 - \mathscr{L}(t_1)\Delta t_1 + \delta\int_{t_1}^{t_2}\mathscr{L}\;dt\\&= \left.(\mathscr{L}\Delta t + \sum_ip_i\delta q_i)\right|_{t_1}^{t_2}\end{align*}

\bullet接下来我们需要将q_i在端点t_1t_2\delta-变分用q_i在正确路径以及错误路径端点处的\Delta-变分\Delta q_i来表示。我们需要找出q_i的两种变分之间的关系。

路径家族的参化方程:

q_i(t_1,\alpha) = q_i(t_1,0) + \alpha\eta_i(t_1)

现在考虑端点t_1。端点变化\Delta t_1导致路径的变分为:

\begin{align*}\Delta q_i(t_1) &= q_i(t_1 + \Delta t_1,\alpha) - q_i(t_1,0)\\&= q_i(t_1 + \Delta t_1,0) + \alpha\eta_i(t_1 + \Delta t_1) - q_i(t_1,0)\end{align*}

\Delta t_1\alpha均足够小时(<\!< 1),

可对q_i(t_1 + \Delta t_1,0) 以及\eta_i(t_1 + \Delta t_1) 使用最优线性展开:

q_i(t_1 + \Delta t_1,0) = q_i(t_1,0) + \frac{\partial q_i}{\partial t}\Delta t_1 = q_i(t_1,0) + \dot{q_i}\Delta t_1

\eta_i(t_1 + \Delta t_1) = \eta_i(t_1) + \frac{\partial \eta_i}{\partial t}\Delta t_1 = \eta_i(t_1) (已将二阶导舍去)

于是

\Delta q_i(t_1) = \dot{q_i}\Delta t_1 + \alpha\eta_i(t_1)

根据\delta-变分的定义,\delta q_i = \frac{\partial q_i}{\partial \alpha}\alpha\alpha <\!< 1),有\delta q_i = \eta_i(t) \alpha

替换后我们就得到了两种变分的关系式:

\Delta q_i(t_1) = \dot{q_i}\Delta t_1 + \delta q_i(t_1)

同理,对于端点t_2,有

\Delta q_i(t_2) = \dot{q_i}\Delta t_2 + \delta q_i(t_2)

所以

\boxed{\begin{cases}\delta q_i(t_1) = \Delta q_i(t_1) - \dot{q_i}\Delta t_1\\\delta q_i(t_2) = \Delta q_i(t_2) - \dot{q_i}\Delta t_2\end{cases}}

再代入作用量变分:

\begin{align*}\Delta \int_{t_1}^{t_2}\mathscr{L}\;dt &= \left[\mathscr{L}\Delta t + \sum_i p_i(\Delta q_i - \dot{q_i}\Delta t)\right]_{t_1}^{t_2}\\&= \left.(\sum_i p_i\Delta q_i - \mathscr{H}\Delta t)\right|_{t_1}^{t_2}\end{align*}

如果没有额外的条件,至此就不能再继续了。

\bullet为了得到最小作用量原理,系统需要满足下列三个条件:

1)系统的拉格朗日函数\mathscr{L}以及哈密顿函数\mathscr{H}不能显含时间,这也意味着\mathscr{H}是一个守恒量。

2)变分需保证不管是正确还是错误路径的哈密顿函数\mathscr{H}都是守恒的。

3)所有错误路径端点坐标的变分需为零(时间的变分\Delta t不为零!)。

第三条是间接为第二条服务的。哈密顿函数\mathscr{H} = p_i\dot{q_i} - \mathscr{L},若想要保证它始终是一个守恒量,q_i(t)就不能改变。由于不同的路径会有不同的速度\dot{q_i}(t),为保证\mathscr{H}守恒,时间必须要变化。

\bullet满足上述三个条件的变分原理则简化为:

\begin{align*}\Delta \int_{t_1}^{t_2}\mathscr{L}\;dt &= -\mathscr{H}(\Delta t_2 - \Delta t_1)\\\Delta\int_{t_1}^{t_2}(\sum_i p_i\dot{q_i} - \mathscr{H})\;dt &= -\mathscr{H}(\Delta t_2 - \Delta t_1)\end{align*}

\implies \boxed{\Delta\int_{t_1}^{t_2}\sum_i p_i\dot{q_i}\;dt = 0}

这个简洁的表达式就是最小作用量原理(principle of least action)

积分\int_{t_1}^{t_2}\sum_i p_i\dot{q_i}\;dt = 0在上世纪的相关书籍中被称为作用量或作用量积分。现在人们口头常说的作用量通常指由哈密顿原理得来的积分\int_{t_1}^{t_2}\mathscr{L}\;dt。所以为了区分,前者有时也被称为简缩作用量(abbreviated action)

\bullet最小作用量原理形式多样。我们知道,在非相对论力学中,如果广义坐标不显含时间,动能T的所有齐次展开项中只有速度的二次项不为零:

T = \frac{1}{2}M_{jk}(q)\dot{q}_j\dot{q}_k

如果势函数不显含速度,拉格朗日函数的齐次展开项中速度的零次项\mathscr{L}_0 = -V,于是:

\mathscr{L} = \mathscr{L}_2 + \mathscr{L}_1 + \mathscr{L}_0 = T - V

正则动量

p_i = \frac{\partial\mathscr{L}}{\partial \dot{q_i}} = \frac{\partial(T - V)}{\partial \dot{q_i}} = \frac{\partial T}{\partial \dot{q_i}}

利用欧拉齐次函数定理

\sum_i p_i\dot{q_i} = \sum_i \frac{\partial T}{\partial \dot{q_i}}\dot{q_i} = 2T

代入最小作用量原理中,可以得到简单形式:

\boxed{\Delta\int_{t_1}^{t_2}T\;dt = 0}

而如果这时系统所受合外力也为零(哈密顿函数不显含时间),总能量守恒:

\frac{d\mathscr{H}}{dt} = \frac{d(T + V)}{dt} = \frac{dT}{dt} = 0

动能也同样不含有时间。

于是最小作用量原理又可被进一步简化为:

\Delta\int_{t_1}^{t_2}T\;dt = T\Delta\int_{t_1}^{t_2}\;dt = T\Delta(t_2 - t_1)

\implies \boxed{\Delta (t_2 - t_1) = 0}

这个惊人结果告诉我们,在位形空间的所有路径中,系统会选择一条使得通过两端点用时最短(严格来讲是取得极值)的路径运动。该形式的最小作用量原理又被称为费马几何光学原理(Fermat's principle in geometrical optics)或最短时间原理(principle of shortest time),它是光线传播需要遵守的法则:一条介于两点之间的光束总是会选择一条用时最短的路径传播。

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

推荐阅读更多精彩内容