Featured image of post Warum der kürzeste Weg nicht der schnellste ist: Die Brachistochrone und die Geburt der Variationsrechnung

Warum der kürzeste Weg nicht der schnellste ist: Die Brachistochrone und die Geburt der Variationsrechnung

Warum erreicht eine reibungsfreie Perle auf einer Zykloide das Ziel schneller als auf einer Geraden? Von Johann Bernoullis Herausforderung 1696 und Newtons Lösung über Fermats optische Analogie bis zur Herleitung der Euler-Lagrange-Gleichung.

1. Historisches Vorspiel und physikalische Intuition: Warum schlägt die Kurve die Gerade?

Stellen wir uns zwei Punkte $A(0,0)$ und $B(x_1, y_1)$ auf unterschiedlicher Höhe in einer vertikalen Ebene vor. Welche Form muss eine reibungsfreie Schiene zwischen $A$ und $B$ haben, damit eine in $A$ aus der Ruhe losgelassene Kugel unter dem Einfluss der Schwerkraft $g$ in der kürzestmöglichen Zeit das Ziel $B$ erreicht? Unsere geometrische Intuition antwortet sofort: „Eine Gerade, denn die Gerade ist die kürzeste Verbindung zwischen zwei Punkten!“ Doch Physik und Mathematik beweisen das Gegenteil: Auf einer steil abfallenden Kurve gewinnt die Kugel mit großem Vorsprung – und liegt $B$ horizontal weit genug entfernt, taucht die schnellste Kurve sogar unter das Niveau von $B$ hinab, um am Ende wieder nach oben ins Ziel zu rollen!

Dieses Problem heißt das Brachistochronen-Problem (von griechisch brakhistos „kürzester“ und khronos „Zeit“). Warum ist der kürzeste Weg nicht der schnellste? Weil die Anfangsgeschwindigkeit bei $A$ gleich null ist ($v_0 = 0$). Auf einer flachen schiefen Ebene ist die Tangentialbeschleunigung $a = g\sin\alpha$ klein, sodass die Kugel zu Beginn quälend viel Zeit bei minimaler Geschwindigkeit verliert. Stürzt die Bahn hingegen am Start vertikal nach unten, wandelt sie potentielle Energie $mgy$ sofort in hohe Geschwindigkeit $v = \sqrt{2gy}$ um. Galileo Galilei erkannte 1638 in seinen Discorsi, dass ein Kreisbogen schneller ist als die gerade Sehne, hielt jedoch fälschlich den Kreis für das absolute Optimum. Im Juni 1696 stellte Johann Bernoulli das Problem als öffentliche Herausforderung an die Mathematiker Europas; am Abend des 29. Januar 1697 löste Isaac Newton nach einem langen Arbeitstag in der Royal Mint die Aufgabe binnen weniger Stunden vor dem Schlafengehen und veröffentlichte sie anonym – woraufhin Bernoulli ausrief: „Tanquam ex ungue leonem (Am Klauenabdruck erkenne ich den Löwen)!“


2. Das Zeit-Funktional und Johann Bernoullis optischer Geniestreich (Snellius-Gesetz)

Mit dem Ursprung in $A(0,0)$, horizontaler $x$-Achse und vertikal nach unten gerichteter $y$-Achse liefert der Energieerhaltungssatz $\frac{1}{2}mv^2 = mgy$ die allein von der Tiefe abhängige Geschwindigkeit $v(y) = \sqrt{2gy}$. Mit dem Bogenelement $ds = \sqrt{1 + (y'(x))^2}\,dx$ wird die gesamte Laufzeit $T[y]$ zu einem Funktional (einer Funktion, die einer ganzen Kurve $y(x)$ eine reelle Zahl zuordnet):

$$ 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. Die Geburt der Variationsrechnung: Euler-Lagrange-Gleichung, Beltrami-Identität und Zykloide

Johann Bernoulli fand 1696 eine brillante physikalische Abkürzung über das Fermatsche Prinzip der kürzesten Zeit: Eine Kugel mit der Geschwindigkeit $v(y) = \sqrt{2gy}$ folgt exakt derselben Bahn wie ein Lichtstrahl in einem optischen Medium mit tiefenabhängigem Brechungsindex $n(y) = c/v(y) \propto y^{-1/2}$! Nach dem Snelliusschen Brechungsgesetz gilt für den Winkel $\alpha$ zur Vertikalen $\frac{\sin\alpha}{\sqrt{2gy}} = \text{const}$. Mit $\sin\alpha = \frac{dx}{ds} = \frac{1}{\sqrt{1+(y')^2}}$ folgt sofort die Differentialgleichung $y(1 + (y')^2) = 2R$.

Um allgemeine Funktionale $J[y] = \int L(x,y,y')\,dx$ systematisch zu minimieren, entwickelten Leonhard Euler (1744) und Joseph-Louis Lagrange (1755) die Variationsrechnung. Setzt man die erste Variation $\delta J[y;\eta] = \left.\frac{d}{d\varepsilon}J[y+\varepsilon\eta]\right|_{\varepsilon=0} = 0$ für beliebige Testfunktionen $\eta(x)$ mit $\eta(0)=\eta(x_1)=0$, liefert partielle Integration die Euler-Lagrange-Gleichung. Da $L(y,y') = \sqrt{\frac{1+(y')^2}{y}}$ nicht explizit von $x$ abhängt, reduziert die Beltrami-Identität $L - y'\frac{\partial L}{\partial y'} = C$ das Problem direkt auf $y(1 + (y')^2) = 2R$. Die Substitution $y = 2R\sin^2(\theta/2) = R(1 - \cos\theta)$ liefert als eindeutige Lösung die Zykloide (Rollkurve):

$$ \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. Numerischer Vergleich, Huygens’ Tautochrone und das Prinzip der kleinsten Wirkung

Für den Weg von $A(0,0)$ nach $B(\pi, 2)$ ist die Zykloide ($S = 4{,}000\,\text{m}$) zwar um $7{,}4\%$ länger als die Gerade ($3{,}724\,\text{m}$), erreicht das Ziel aber um $18{,}5\%$ schneller! Überdies bewies Christiaan Huygens bereits 1659, dass genau dieselbe Zykloide auch die Tautochrone (Isochrone) ist: Egal aus welcher Starthöhe $y_0$ eine Kugel auf einer Zykloidenbahn losgelassen wird, sie benötigt bis zum tiefsten Punkt immer exakt dieselbe Zeit $T = \pi\sqrt{R/g}$! Aus der Variationsrechnung der Brachistochrone erwuchsen das Hamiltonsche Prinzip der kleinsten Wirkung ($\delta S = 0$), das Noether-Theorem, Einsteins Geodätengleichung der Allgemeinen Relativitätstheorie, Feynmans Pfadintegral der Quantenmechanik und die moderne optimale Steuerungstheorie.

KandidatenkurveGleichung $y(x)$Startsteigung $y'(0)$Bogenlänge $S$ ($A(0,0)\to B(\pi,2)$)Normierte Laufzeit ($\sqrt{g}=1$)Zeitverlust gegenüber Zykloide
1. Gerade (Kürzeste Strecke)$y = \frac{2}{\pi}x$$0{,}6366$ ($32{,}5^\circ$)$3{,}724\,\text{m}$ (Minimum)$3{,}7242$$+18{,}5\%$ langsamer
2. Parabel$y = 2\sqrt{x/\pi}$$\infty$ (vertikal)$3{,}960\,\text{m}$$3{,}2185$$+2{,}4\%$ langsamer
3. Kreisbogen (Galileis Tipp)Kreis tangential zur $y$-Achse$\infty$ (vertikal)$4{,}062\,\text{m}$$3{,}1829$$+1{,}3\%$ langsamer
4. Zykloide (Brachistochrone)$x=\theta-\sin\theta,\; y=1-\cos\theta$$\infty$ (vertikal)$4{,}000\,\text{m}$$\pi \approx 3{,}1416$ (Absolutes Minimum)Theoretisches Optimum ($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