Featured image of post Почему кратчайший путь не самый быстрый: брахистохрона и рождение вариационного исчисления

Почему кратчайший путь не самый быстрый: брахистохрона и рождение вариационного исчисления

Почему бусинка без трения скатывается по циклоиде быстрее, чем по прямой? От вызова Иоганна Бернулли 1696 года и ночного решения Ньютона до оптической аналогии Ферма, уравнения Эйлера — Лагранжа и таутохроны Гюйгенса.

1. Исторический пролог и физическая интуиция: почему кривая быстрее прямой?

Представьте две точки $A(0,0)$ и $B(x_1, y_1)$, закрепленные на разной высоте в вертикальной плоскости. Какую форму должен иметь гладкий желоб без трения между $A$ и $B$, чтобы стальная бусинка, отпущенная без начальной скорости из точки $A$ в однородном поле тяжести $g$, достигла точки $B$ за кратчайшее возможное время? Геометрическая интуиция подсказывает: «По прямой, ведь прямая — кратчайшее расстояние между двумя точками!». Однако физика и математика доказывают обратное: по круто ныряющей кривой бусинка приходит к финишу намного раньше, а если точка $B$ достаточно удалена по горизонтали, самая быстрая кривая даже опускается ниже уровня точки $B$, чтобы затем подняться к финишу!

Эта классическая задача называется задачей о брахистохроне (от греч. brakhistos — «кратчайший» и khronos — «время»). Почему кратчайший путь не самый быстрый? Потому что в начальной точке $A$ скорость равна нулю ($v_0 = 0$). На пологой прямой наклонной плоскости касательное ускорение $a = g\sin\alpha$ невелико, и бусинка тратит огромное время на начальном участке, двигаясь с ничтожной скоростью. Напротив, траектория, отвесно падающая вниз на старте, мгновенно переводит потенциальную энергию $mgy$ в высокую скорость $v = \sqrt{2gy}$. В 1638 году Галилей в «Беседах и математических доказательствах» показал, что дуга окружности быстрее прямой хорды, но ошибочно счел окружность абсолютным оптимумом. В июне 1696 года Иоганн Бернулли бросил вызов математикам Европы; вечером 29 января 1697 года, вернувшись со службы на Королевском монетном дворе в Лондоне, Исаак Ньютон решил задачу за несколько часов перед сном и опубликовал решение анонимно, заставив Бернулли воскликнуть: «Tanquam ex ungue leonem (По когтю узнаю льва)!».


2. Функционал времени спуска и оптическая аналогия Иоганна Бернулли (закон Снеллиуса)

Поместив начало координат в точку $A(0,0)$ с горизонтальной осью $x$ и направленной вертикально вниз осью $y$, из закона сохранения механической энергии $\frac{1}{2}mv^2 = mgy$ находим скорость $v(y) = \sqrt{2gy}$. С учетом элемента дуги $ds = \sqrt{1 + (y'(x))^2}\,dx$ полное время спуска $T[y]$ представляет собой функционал (отображение, ставящее в соответствие целой кривой $y(x)$ вещественное число):

$$ T[y] = \int_{A}^{B} \frac{ds}{v(y)} = \frac{1}{\sqrt{2g}} \int_{0}^{x_1} \sqrt{\frac{1 + (y'(x))^2}{y(x)}}\,dx, \qquad \frac{\sin\alpha(y)}{\sqrt{2gy}} = \text{const} $$

3. Рождение вариационного исчисления: уравнение Эйлера — Лагранжа, тождество Бельтрами и циклоида

Иоганн Бернулли в 1696 году нашел гениальное физическое решение с помощью принципа минимального времени Ферма: частица, движущаяся со скоростью $v(y) = \sqrt{2gy}$, следует по той же траектории, что и луч света в слоистой оптической среде с показателем преломления $n(y) = c/v(y) \propto y^{-1/2}$! Согласно закону преломления Снеллиуса, угол $\alpha$ с вертикалью удовлетворяет условию $\frac{\sin\alpha}{\sqrt{2gy}} = \text{const}$. Подстановка $\sin\alpha = \frac{dx}{ds} = \frac{1}{\sqrt{1+(y')^2}}$ сразу дает дифференциальное уравнение $y(1 + (y')^2) = 2R$.

Для систематической минимизации общих функционалов $J[y] = \int L(x,y,y')\,dx$ Леонард Эйлер (1744) и Жозеф Луи Лагранж (1755) создали вариационное исчисление. Приравнивание к нулю первой вариации $\delta J[y;\eta] = 0$ для произвольной вариации $\eta(x)$ с закрепленными концами после интегрирования по частям приводит к уравнению Эйлера — Лагранжа. Поскольку $L(y,y') = \sqrt{\frac{1+(y')^2}{y}}$ не зависит явно от $x$, тождество Бельтрами $L - y'\frac{\partial L}{\partial y'} = C$ сразу дает $y(1 + (y')^2) = 2R$. Параметрическая подстановка $y = 2R\sin^2(\theta/2) = R(1 - \cos\theta)$ доказывает, что единственной брахистохроной является циклоида:

$$ \frac{\partial L}{\partial y} - \frac{d}{dx}\!\left(\frac{\partial L}{\partial y'}\right) = 0 \implies L - y'\frac{\partial L}{\partial y'} = C \implies y\Big(1 + (y')^2\Big) = 2R $$
  flowchart TD
    A["1696: Brachistochrone Problem (Minimize T[y] = ∫ √(1+y'^2)/√(2gy) dx)"] --> B["Snell's Law Optical Analogy (sin α / √y = const)"]
    A --> C["First Variation δT[y] = 0 → Euler-Lagrange Equation & Beltrami Identity"]
    B --> D["ODE: y (1 + y'^2) = 2R"]
    C --> D
    D --> E["Parametric Cycloid: x = R(θ - sin θ), y = R(1 - cos θ)"]
    E --> F["Huygens Tautochrone Property (T = π√(R/g) Independent of Starting Height)"]
    C --> G["Hamilton's Principle of Least Action δS = 0 (Modern Physics & Optimal Control)"]

4. Численное сравнение, таутохрона Гюйгенса и наследие принципа наименьшего действия

Между точками $A(0,0)$ и $B(\pi, 2)$ дуга циклоиды ($S = 4{,}000\,\text{м}$) на $7{,}4\%$ длиннее отрезка прямой ($3{,}724\,\text{м}$), однако бусинка проходит циклоиду на $18{,}5\%$ быстрее! Более того, еще в 1659 году Христиан Гюйгенс доказал, что эта же самая циклоида является и таутохроной (изохроной): с какой бы начальной высоты на циклоиде ни отпустили бусинку, время спуска до нижней точки всегда строго одинаково и равно $T = \pi\sqrt{R/g}$! Из вариационного исчисления брахистохроны выросли принцип наименьшего действия Гамильтона ($\delta S = 0$), теорема Нётер, уравнения геодезических общей теории относительности Эйнштейна, континуальный интеграл Фейнмана и современная теория оптимального управления.

Кривая-кандидатУравнение $y(x)$Начальный наклон $y'(0)$Длина дуги $S$ ($A(0,0)\to B(\pi,2)$)Нормированное время ($\sqrt{g}=1$)Проигрыш циклоиде по времени
1. Прямая (Кратчайшее расстояние)$y = \frac{2}{\pi}x$$0{,}6366$ ($32{,}5^\circ$)$3{,}724\,\text{м}$ (Минимум)$3{,}7242$Медленнее на $+18{,}5\%$
2. Парабола$y = 2\sqrt{x/\pi}$$\infty$ (вертикально)$3{,}960\,\text{м}$$3{,}2185$Медленнее на $+2{,}4\%$
3. Дуга окружности (Галилей)Окружность, касательная к оси $y$$\infty$ (вертикально)$4{,}062\,\text{м}$$3{,}1829$Медленнее на $+1{,}3\%$
4. Циклоида (Брахистохрона)$x=\theta-\sin\theta,\; y=1-\cos\theta$$\infty$ (вертикально)$4{,}000\,\text{м}$$\pi \approx 3{,}1416$ (Абсолютный минимум)Теоретический оптимум ($0{,}0\%$)
$$ \boxed{\; x(\theta) = R(\theta - \sin\theta), \qquad y(\theta) = R(1 - \cos\theta), \qquad T_{\text{tautochrone}}(\theta_0) = \pi\sqrt{\frac{R}{g}} \;} $$
comments powered by Disqus