分析力学基本原理介绍4.2:变分法的简单应用

通过前篇内容我们对变分原理有了一定了解,并利用变分法得到了能使形如

的泛函取得驻值的条件:

接下来,我打算展示一些关于这个方程的简单应用。

需要注意的是,由于习惯,我们通常将函数视为因变量,为自变量,所以我为了保持一致,将上面积分中的考虑成了单向的,但这个习惯与大部分问题中的变量和充当的角色刚好相反(通常才是单向的,并非),导致有的问题可以使用原来的表示习惯得到正确答案,有的却不能。所以在下面的例子中,我会变换着使用二者。


1.平面内两点间的最短距离

对,这个方程可被用来证明平面内两点之间的最短距离——一个小学生都知道的结论——直线。

平面内连接了任意两点,的曲线长度可以用积分表示为:

若曲线表达式采用,则

在笛卡尔平面内,长度微分元

代入积分中:

我们发现,该积分正好具有之前所讨论的泛函的形式。所以对应地,

于是

代入之前得到的方程我们可以发现:

其中是一个常数。

解得:

等式右边均是常数,不妨令其整体为,则

这正是直线方程的表达式。

其中常数,可通过代入端点值,求解。

顺带一提,利用上面的方法,我们不仅证明了该直线方程能使泛函取得极值,而且是极小值。


2.最小回转曲面:

我们把空间中连接了两定点的任意曲线环绕转轴(比如轴)弧度形成的曲面叫做回转曲面(surface of revolution)或者旋转曲面。

同样地,使用开头给出的方程我们也可以得到空间中两点间一条使得回转曲面面积最小的曲线。

考虑转轴为轴,曲线经过了点和。

使用极坐标,若曲线形成的回转曲面上任意一点距轴的距离为,转动微小角度后形成的环带面积为:,其中是之前的长度微分元(这个例子比较特殊,使用在最后也可以得到同样的结果,不信你自己试)。

则:

对比之前的泛函可得:

于是

代入方程得到:

,其中是一个常数。

简单的代数操作后可以得到微分方程:

或者,拆分变量后:

,其中是第二个常数。

这是一个分母含根号的基本积分,结果是可以用初等函数表示的。大部分同类型积分都可以考虑配完全平方或者使用三角、双曲换元的方法来计算。

将常数从根号中提取出来,

利用关系,可使用双曲余弦函数换元:,

得到:

或者写成

这是悬链线(catenary)方程,常数,可通过代入端点坐标解得。

顺带一提,任何满足之前方程的曲线都是按和的倍数变化,独立变量也暗示了经过同样的两定点可以存在一条甚至多条悬链线,在特殊情况下也许一条也会不存在。更为复杂的情形就需要再考虑函数的连续性以及一二阶导函数的存在情况,毕竟这些是限制满足之间方程的解的条件。


3.最速落径

最速落径(The brachistochrone),又被称为捷线。

我把它放在最后就是想强调一下它的地位——就是那种“the last but not the least”的地位。嗯,它很有名,曾经被一位叫做“不努力”的数学家分析过,属于变分法的开创性问题,基本也属于那种每个学过高物或者高数的人章口就莱的陈词滥调。

至于为什么叫“不努力”,呸,至于为什么叫“The brachistochrone”,其实该词源于古希腊语“”,意思是“最短的时间”,为什么是“最短”的“时间”?不应该是“最速”的“落径”吗?希腊语的命名其实是很准确的——至少在我看来——因为它忠于题设:

寻找一条连接两个定点的曲线,使得一个只受重力影响的微粒,沿该曲线从较高点抵达较低点所用的时间最短

与前面两题做法相似,我们首先要找到那个对应的函数。

既然要使时间最短,我们的优化对象肯定就是一个表示时间的泛函了:

微粒的速度为:

微分元可以表示成:

所以

同样是之前的长度微分元,但这次我将使用:

至于微粒的速度,一个如此简单的系统,用初中的知识足以解决。

使用牛顿定律或者能量守恒定律随便搞两下就可以得到:

,其中是重力加速度。

代入积分,

对比本篇最上方的积分,

于是

代入方程,

我们需要计算积分:

对于这类积分,我们考虑使用三角函数的换元。

令,则

于是可以得到:

为了让结果更简洁,可选择将替换为;将替换为:

,其中和均为常数。

同时,

最后我们得到的是一组含的参数方程:

曲线是最速落径,是属于摆线(cycloid)的一种可以被形象地理解成是由一块粘在行进的自行车车胎上的泡泡糖形成的运动轨迹。

更多关于摆线的分析可自行查阅。

你可能感兴趣的:(分析力学基本原理介绍4.2:变分法的简单应用)