【原】著名的数学难题——最速落径问题，求解这个问题莱布尼茨用了6个多月，伯努利用了两个星期，牛顿只用了一个晚上

1696年6月，著名数学家约翰·伯努利（ Johann Bernoulli ）)在德国一份科学期刊《博学报》（Acta Eruditorum）上发表了以下问题：

在铅直平面上两点A，B之间要连一条曲线，使得不受摩擦的质点在重力的作用下沿这条曲线由A运动到B所需要的时间最少？

• 图1：从A到B，三种可能的最优路径

下图显示了约翰·伯努利和1696年6月用拉丁文在《博学报》上对该问题的表述。

• 图2

这一数学难题被称为捷线（最速落径）。尽管约翰·伯努利自己已经知道如何解决这个问题，但他还是挑战了欧洲的其他数学家，并给他们6个月的时间来解决这个问题。然而，在那之后，没有人给出任何答案。就连历史上最伟大的知识分子之一戈特弗里德·莱布尼茨（Gottfried Leibniz）也要求延迟最后期限。1697年1月29日下午，艾萨克·牛顿在他的邮件中（一封来自伯努利的信）发现了这个问题。然后，他在夜间解出了这个难题，并以匿名方式寄回了答案。

下面是牛顿手写的答案。这个故事让我们对牛顿的天赋有了一些了解，因为约翰·伯努利花了两周的时间才解出它。

• 图3：牛顿手写的捷线问题解决方案。

牛顿手写解的翻译是：

从给定点A出发，画一条平行于水平面的无界直线APCZ，在这条直线上描述任意摆线AQP，在Q点上与直线AB相交（并在必要时延伸），然后另一个摆线ADC的底和高[as AC: AP]应分别为前一个的底和高AB到AQ。这条最近的摆线将穿过B点，成为一条曲线，在这条曲线上，一个重物在自身重量的作用下，最迅速地从A点到达B点。

要了解牛顿对上述解的详细过程，请私信我。

最速落径曲线

最速落径曲线是一条位于二维平面上的曲线，有一个初始点A和一个终点B，仅受重力作用的一个质点从A点到B点时间最短的路径。

求曲线的问题有以下假设：

• 曲线上没有摩擦

• 质点开始时是静止的

• 引力场是常数是g

• 图4：A和B之间的一条可能路径Γ

现代解

假设解是函数y=y(x)，为了方便起见，我们选择初始点A =(0,0)。最后一个点定义为B = (a, b)。由于质点最初处于静止状态，由能量守恒将得到：

• 式1：能量守恒

然后我们把dt写成：

• 式2：用x和y表示的无穷小区间dt。

质点从A =(0,0)到B = (a, b)的总时间则为：

• 式3：质点从(0,0)到(a,b)的总时间间隔T。

数学对象T依赖于函数y(x)，因此它被称为泛函（函数的函数）。泛函只依赖于（一个或多个）变量，而不依赖于完整的函数。

我们要解决的问题是找出函数y(x)使总时间t最小。为此，我们需要学习一个叫做变分法的数学。

变分法

考虑一个函数ψ(x)，ψ满足以下条件，即ψ(x_0)=y_0和ψ(x_1)=y_1。考虑第二个非常接近第一个的函数，把它写成:

• 式4：第二个函数，非常接近第一个，u(x)所满足的条件。

请注意，关于ψ(x)的条件必须满足上述关于u(x)的条件。

• 图5：函数ψ(x)和另一个函数。

现在考虑以下函数：

• 式5：一个函数，其被积函数L显式地依赖于x, y和y'。

注意，通过改变L(x, y, y ')，我们得到了定积分S[y(x)]的不同值。现在我们考虑随ψ(x)变化而变化的L：

• 式6：ψ(x)和ψ ' (x)变化时L的变化。

对两边积分，对第二项进行分部积分，利用u(x)所满足的条件，得到积分S的变化量如下：

• 式7：积分S经过一个小的变化后的变化。

如果S是最小值，δS=0。由于u(x)是任意的，必须有：

• 式8：δS=0的必要条件。

当y(x)等于使L为极值的函数ψ(x)时，括号内的表达式消失。简化符号，我们得到了著名的欧拉-拉格朗日方程：

• 方程9：欧拉-拉格朗日方程。

我们用它求出式3中最短的时间，其中：

• 式10：式3的被积函数，代入欧拉-拉格朗日方程式。

经过代数的几步，我们得到以下微分方程及其相应的解：

• 式11：表示摆线的参数方程。

式中k为某常数（依赖于边界条件），变量的变化如下：

• 式12：用来推导方程式11的变量的改变。

• 图6：滚动圆周长上的一点产生摆线

这些参数方程描述了一个摆线，它是使T最小化的曲线，如下图所示。

• 图7：这个图显示了最速落径的曲线是摆线。

牛顿的解

1699年，数学家、自然哲学家、天文学家、发明家、宗教活动家法蒂奥（ Fatio）发表了一篇论文“关于最速落径曲线的双重几何研究”，其中包含了另一种解决捷线问题的方法。

• 图8：法蒂奥的肖像

大卫·格雷戈里要求牛顿简化法蒂奥的解。这一节，我将描述牛顿的解。

• 图9：牛顿和他发送给大卫·格雷戈里的解。

在图10中定义了相关的量。我们首先写出：

• 图10：在法蒂奥-牛顿的解中使用的构造

现在，我们从基本运动学得知，下落的质点在x处的速度为：

• 式13：下落质点的速度与高度x有关。

质点沿着ENG移动所需的时间正比于：

• 式14：质点从E到G的时间。

现在定义：

• 式15：R²和S²的定义。

使总的时间t相对于q最小，经过一些简单的代数运算，我们得到了一个摆线的微分方程，由式11给出。