变分法理解3——例题求解

1,082 阅读2分钟
原文链接: zhuanlan.zhihu.com

本文是变分法理解系列的第3篇文章

在第2篇文章中,推导出了欧拉-拉格朗日方程,本文用它来解决第1篇文章中的两个简单的泛函问题

两点之间的最短路径

问题的具体描述见第1篇文章,在第1篇文章中求的是从坐标原点(0,0)到点(a,b)的连接曲线是 y = y(x),不失一般性,这里将其改成求点 (x_0,y_0) 到点 (x_1,y_1) 的最短连线。

两点之间最短路径问题就是如何找出曲线y(x),使得曲线的总弧长,也就是泛函

J[y(x)] = \int_{x_0}^{x_1} (1 + y'^2 )^{1/2}dx

最小

因为 F = (1+y'^2)^{1/2} ,

所以

\frac{\partial F}{\partial y}=0,\qquad \frac{\partial F}{\partial y'} = \frac{y'}{(1+y'^2)^{1/2}}

代入欧拉-拉格朗日方程为:

\frac{\partial F}{\partial y}-\frac{d}{dx}\frac{\partial F}{\partial y'}=0

因此:

\frac{\partial F}{\partial y'}= C = \frac{y'}{(1+y'^2)^{1/2}}

其中C是常数,因此:

y'=\frac{C}{(1-C^2)^{1/2}}

可见曲线y(x)的斜率 y' 是一个常数,设 y(x)=ax+b ,由已知条件 y ( x_0 ) = y_0 , y (x_1) = y_1 ,可得两点间最短连线是:

y(x)= \frac{y_1 - y_0}{ x_1 - x_0 } (x-x_0) + y_0

最速降线问题

问题的具体描述见第1篇文章,重物由 O 点运动到 A 点所需时间 t 是 泛函

t=J[y(x)]=\int_0^a \sqrt{\frac{(1 + y'^2)}{2gy}}dx

并且满足条件:

y ( 0 ) = 0 , y ( a ) = b

最速降线的问题就是在所有连续函数 y(x) 中,求出一个函数 y 使时间 t 取最小值

首先证明一个定理

若F只依赖于 y 和 y',F(y,y') 具有二阶偏导数,y(x)也具有二阶偏导数,那么:

F(y,y') - y' \frac{\partial F}{\partial y'} = C

证明:

\begin{aligned}\frac{d}{dx}[F(y,y') - y' \frac{\partial F}{\partial y'}] &= \frac{\partial F}{\partial y} y' + \frac{\partial F}{\partial y'}y''- y''\frac{\partial F}{\partial y'}-y'\frac{d}{dx}(\frac{\partial F}{\partial y'}) \ &= y' \frac{d}{dx}(\frac{\partial F}{\partial y'})-y'\frac{d}{dx}(\frac{\partial F}{\partial y'}) \ &= 0 \end{aligned}

因此

\displaystyle C = F(y,y') - y' \frac{\partial F}{\partial y'}

定理得证

\displaystyle F(y,y')=\sqrt{\frac{(1 + y'^2)}{2gy}} 代入定理,有:

\begin{aligned} C=\sqrt{\frac{(1 + y'^2)}{2gy}}-y'\frac{y'}{\sqrt{2gy}\sqrt{(1 + y'^2)}} = \frac{1}{\sqrt{2gy}\sqrt{(1 + y'^2)}} \end{aligned}

\displaystyle \frac{1}{2gC^2}=r ,因此:

y(1+y'^2)=r

y'=\cot \theta ,则:

y=\frac{r}{1+y'^2}=r \sin^2 \theta =\frac{r}{2}(1-\cos 2\theta)

可以求出:

dx=\frac{dy}{y'}=\frac{r \sin 2\theta d\theta}{\cot \theta}=\frac{r \sin \theta \cos \theta d\theta}{\cot \theta}=r(1-\cos 2 \theta) d \theta

积分可得:

x=\frac{r}{2} (2 \theta -\sin 2 \theta) + C1

根据曲线过原点 ( 0 , 0 ) 及 ( a , b ) 可求出 C1 = 0 及 r ,这样,因此令 \displaystyle t=2\theta,d=\frac{r}{2} ,最速降线的解为:

x=d(t-\sin t)

y=d(1-\cos t)

由此可知,最速降线是旋转轮线的一段,旋转轮线又叫摆线,是由一个圆沿一条直线运动时,圆边界上一定点所形成的轨迹。