变分法证明两点之间线段最短

最新推荐文章于 2024-02-18 14:06:44 发布

丨Modred丨

最新推荐文章于 2024-02-18 14:06:44 发布

阅读量1.5w

点赞数 5

分类专栏：算法

算法专栏收录该内容

24 篇文章 1 订阅

订阅专栏

传送门https://zhuanlan.zhihu.com/yueaptx

变分法简介Part 1.（Calculus of Variations）

Dr.Stein

计算力学

关注他

283 人赞了该文章

泛函数 (Functionals)

简而言之，泛函数是函数的函数，即它的输入是函数，输出是实数。而这个输出值取决于一个或多个函数（输入）在一整个路径上的积分而非像一般函数一样取决于离散的变量。这样说可能还是比较抽象，不过坚持看到下文的Example 1就可以更好理解了。

通常在变分法中，泛函数是一个积分，记做。

$I(y)=\int_{x1}^{x2}Fdx$

其中我们想要通过选择被积函数来最大化或最小化泛函数的值。同时我们称为拉格朗日函数(Lagrange function)。可以是函数 y(x) 和 y(x) 各阶导数的函数(以下 y(x) 均简写成)。为了说明方便，我们先姑且设是和的函数，所以我们可以进一步将泛函数写成：

$I(y)=\int_{x1}^{x2}F(y,y';x)dx$

积分里面我用分号;将和前面的隔开代表和是的函数。一般和的函数关系是已知的，所以想要最大或最小化泛函数，实际上是通过选择适当的函数 y(x) 。

为了透彻理解这个概念，我们可以来看一个简单的例子。

Example 1.

一个最简单直观的例子是求两个固定点之间的最短路径。当然大家都知道两点之间直线最短，这里可以用变分法做出解释。

如上图所示路径是一任意路径，我们取区中一小段微元，可以容易计算微元断的长度为：

ds^2=dx^2+dy^2=[1+(y')^2]dx^2 ,即：

$ds=\sqrt{1+(y')^2}dx$

积分得到总的路径长度为：

$L=\int_{x1}^{x2}ds=\int_{x1}^{x2}\sqrt{1+(y')^2}dx$

这个例子中，是泛函数， $\sqrt{1+(y')^2}$ 是拉格朗日函数，我们想要找一个函数 y(x) 使得泛函数的值最小。这次Part 1.的任务就是为解决这个问题做准备。Part 2.中我们会用变分法证明这个 y(x) 确实是直线的方程。

2. 泛函数的极值

这里重申下，泛函数在区间 [x_1,x_2] 上的值取决于积分路径的选择，即取决于函数 y(x) 的选择。我们有理由假设存在一个这样的 y(x) ，可以使得泛函数取到极值。而在这个 y(x) 附近的任意路径我们记做 $\tilde{y}(x)$ 。另外，我们假设 y(x) 两阶可微。通过引入一个微小量 $\epsilon\ll 1$ 和一个任意可微函数 $\eta(x)$ ，我们可以用 y(x) 表示 $\tilde{y}(x)$ :

$\tilde{y}(x)=y(x)+\epsilon\eta(x)$

这样做的好处是对于一个给定的 $\eta(x)$ ，我们可以通过改变 $\epsilon$ 的值来得到无穷多的路径，同时对于任何 $\eta(x)$ ，当 $\epsilon=0$ 的时候， $\tilde{y}(x)$ 和 y(x) 重合。

图像直观表示如下图：

由于在边界条件的限制， $\eta(x_1)=\eta(x_2)=0$ 。这样就能保证 $\tilde{y}(x)$ 可以通过两个固定端点。

这时我们可以说， y(x) 所对应的泛函数的值是泛函数 $\tilde{I}=\int_{x_1}^{x_2}F(\tilde{y},\tilde{y}';x)dx$ 的极值。我们可以进一步用 $\epsilon$ 表示 $\tilde{I}$ ：

$\tilde{I}=\int_{x_1}^{x_2}F(\tilde{y},\tilde{y}';x)dx=\int_{x1}^{x_2}F(y+\epsilon\eta,y'+\epsilon\eta';x)dx$

虽然 y(x) 未知，但是根据之前的合理假设， y(x) 是一个存在的确定函数。所以根据上式，如果给定一个特定的 $\eta(x)$ ， $\tilde{I}$ 的变化只取决于 $\epsilon$ 的变化。所以我们现在可以把 $\tilde{I}$ 看做是 $\epsilon$ 的函数。用泰勒展开公式将 $\tilde{I}$ 在 $\epsilon=0$ 处展开得到：

$\tilde{I}(\epsilon)=\tilde{I}|_{\epsilon=0}+(\frac{d\tilde{I}}{d\epsilon})\Big|_{\epsilon=0}\cdot\epsilon+ (\frac{d^2\tilde{I}}{d\epsilon^2})\Big|_{\epsilon=0}\cdot\frac{\epsilon^2}{2!} +\cdot \cdot \cdot =\tilde{I}_0+\tilde{I}_1\epsilon+\tilde{I}_2\epsilon^2+\cdot \cdot \cdot$

很明显，当 $\epsilon=0$ 时， $\tilde{I}|_{\epsilon=0}=I$ ,带入上式可得到：

$\tilde{I}-I=\tilde{I}_1\epsilon+\tilde{I}_2\epsilon^2+\cdot \cdot \cdot$

这里我们记 $\delta I=\tilde{I}_1\epsilon=\frac{d\tilde{I}}{d\epsilon}\Big|_{\epsilon=0}\cdot\epsilon$ ，并称之为一阶变分。同理二阶变分为 $\delta I^2=\tilde{I}_2\epsilon^2$ 。

（这里插一句变分和微分的区别。变分在上图的直观解释是 $\tilde{y}$ 和在竖直方向上的距离，称之为 $\delta y$ ，所以这个差是在同一个x上计算的。而微分则是由于x的微小变动引起的y的变动。）

然后我们可以类比求函数极值时的做法。求函数极值时，我们会令函数的一阶导数为零。这里同样，为了求泛函数 $\tilde{I}$ 的极值，我们令一阶变分 $\delta I=0$ 。现在我们计算化简 $\delta I$ :

$\delta I=(\int_{x_1}^{x_2}\frac{d\tilde{F}}{d\epsilon}\Big|_{\epsilon=0} dx)\cdot\epsilon$
$\frac{d\tilde{F}}{d\epsilon}\Big|_{\epsilon=0} \cdot\epsilon=(\frac{\partial\tilde{F}}{\partial\tilde{y}}\cdot\frac{d\tilde{y}}{d\epsilon}+ \frac{\partial\tilde{F}}{\partial\tilde{y'}}\cdot\frac{d\tilde{y'}}{d\epsilon})\Big|_{\epsilon=0}\cdot\epsilon$

因为 $\tilde{y}(x)=y(x)+\epsilon\eta(x)$ , 不难得到： $\frac{d\tilde{y}}{d\epsilon}=\eta$ , $\frac{d\tilde{y'}}{d\epsilon}=\eta'$ ,另外我们有 $\delta y=\epsilon \eta,\delta y'=\epsilon\eta'$

又因为当 $\epsilon\rightarrow 0$ 时， $\tilde{F}\rightarrow 0, \tilde{y}\rightarrow y,\tilde{y}'\rightarrow y'$ ,将这些式子带入原式可以得到：

$\delta I=\int_{x_1}^{x_2}(\frac{\partial F}{\partial y}\delta y+\frac{\partial F}{\partial y'}\delta y')dx$

终于到最后一步啦，分部积分一下得到：

$\delta I=\int_{x_1}^{x_2}(\frac{\partial F}{\partial y}-\frac{d}{dx}(\frac{\partial F}{\partial y'} ))\delta ydx+\frac{\partial F}{\partial y'}\delta y\Big|_{x_1}^{x_2}$

另 $\delta I=0$ 就可以解得最小化泛函数的y啦。我们注意到 $\delta I$ 有两个部分。对于第一个积分部分，由于 $\delta y$ 是任意的，所以要想使这个部分等于零，需要保证 $^{[1]}$ ：

$\frac{\partial F}{\partial y}-\frac{d}{dx}(\frac{\partial F}{\partial y'} )=0$ ( $x_1\leq x\leq x_2$ )

这就是传说中的欧拉-拉格朗日方程(E-L equation)。

而第二部分等于零则是边界条件。

在Part 2.，我们会以用这次介绍的内容和上述方程解决两点之间直线最短的问题为开头，继续介绍变分法。

---------------------------------------------------------------------------------------------------------------------------

注[1]:

假设 x_1 和 x_2 是给定的常数， $\phi(x)$ 是一个特定的在 $x_1\leq x\leq x_2$ 上连续的函数，那么如果对于任意连续可微的函数 $\eta(x)$ 都成立 $\int_{x_1}^{x_2}\phi(x)\eta(x)dx=0$ ，则 $\phi(x)=0$ ( $x_1\leq x\leq x_2$ )。