最速降线问题——欧拉-拉格朗日方程的推导

MedivhMai

已于 2022-08-11 17:24:38 修改

阅读量4.8k

点赞数 7

文章标签：算法

于 2022-08-11 17:21:48 首次发布

本文链接：https://blog.csdn.net/m0_63039263/article/details/126277590

版权

由机械能守恒：

$mgy=\frac{1}{2}mv^2$

得 $v$ ：

$v=\sqrt{2gy}$

由勾股定理得：

$ds = \sqrt{dx^{2}+dy^{2}} = \sqrt{1+\frac{dy^2}{dx^2}}dx$

对弧长积分得：

$T = \int_{a}^{b}\frac{ds}{v} = \int_{x_{a}}^{x_{b}} \frac{ \sqrt{1+\frac{dy^2}{dx^2}}}{\sqrt{2gy}}dx= \int_{x_{a}}^{x_{b}}\sqrt{\frac{1+(y')^2}{2gy}}dx$

函数y是任意一条两端固定的曲线函数，将T表示成函数的函数：

$T(y)=\int_{x_{a}}^{x_{b}}f(y,y' ,x)dx$

假设y是一个极大值, 当y发生微小的变化时, y的变化率趋近于0，这种微小的变化记作 $\delta{T}$ .

$\delta{T} = T(y + \delta{y}) - T(y) = \int_{x_{a}}^{x_{b}}[f(x,y + \delta y,y' + \delta y' ) - f(x, y, y')]dx$

根据多元积分链式法则得：

$f(x,y + \delta y,y' + \delta y' ) = f(x,y,y') + \frac{\partial f}{\partial y}\delta y + \frac{\partial f}{\partial y'}\delta y'$

带入上式得：

$\delta T = \int_{x_{a}}^{x_{b}}(\frac{\partial f}{\partial y}\delta y + \frac{\partial f}{\partial y'}\delta y')dx$

如果将 $\delta y$ , $\delta y'$ 理解为与原函数的差函数, 那么变化量和 $dx$ 无关.

由于两端固定, 从数学上可证明微分、变分的交换性质：

$d(\delta y) = \delta(dy)$

进而得：

$\delta y' = \frac{d}{dx} \delta y$

带入上式得：

$\delta T = \int_{x_{a}}^{x_{b}}(\frac{\partial f}{\partial y}\delta y + \frac{\partial f}{\partial y'}(\delta y)')dx$

积分第二项分部积分得：

$\int_{x_{b}}^{x_{a}} \frac{\partial f}{\partial y'}(\delta y)'dx = \frac{\partial f}{\partial y'}\delta y \bigg|_{x_{a}}^{x_{b}} - \int_{x_{b}}^{x_{a}} ( \frac{d}{dx} \frac{\partial f}{\partial y'}\delta y) dx$