Kepler encontró la sustitución cuando intentaba resolver la ecuación $$\ell=mr^2\frac{d\nu}{dt}=\text{constante}$$ donde $\ell$ es el momento angular orbital, $m$ es la masa del cuerpo en órbita, la anomalía verdadera $\nu$ es el ángulo en la órbita pasado el periapsis, $t$ es el tiempo, y $r$ es la distancia al atractor. Esta es la segunda ley de Kepler, la ley de las áreas equivalente a la conservación del momento angular. Entonces, la primera ley de Kepler, la ley de la trayectoria, es $$r=\frac{a(1-e^2)}{1+e\cos\nu}$$ donde $a$ y $e$ son el semieje mayor y la excentricidad de la elipse. Para obtener $\nu(t)$, necesitas resolver la integral $$\int\frac{d\nu}{(1+e\cos\nu)^2}$$. Así que para relacionar el área barrida por un segmento de línea que une el cuerpo en órbita al atractor, Kepler dibujó un pequeño dibujo [![Figura 1][1]][1]
El atractor está en el foco de la elipse en $O$ que es el origen de coordenadas, el punto del periapsis está en $P$, el centro de la elipse está en $C$, el cuerpo en órbita está en $Q$, habiendo recorrido el área azul desde el periapsis y ahora en una anomalía verdadera de $\nu$. Para realizar la integral dada anteriormente, Kepler amplió el dibujo por un factor de $1/\sqrt{1-e^2}$ en la dirección $y$ para convertir la elipse en un círculo.
[![Figura 2][2]][2] [1]: https://i.sstatic.net/TlBAu.png [2]: https://i.sstatic.net/MVCCe.png
El cuerpo en órbita se ha movido hacia arriba a $Q^{\prime}$ a la altura $$y=\frac{a\sqrt{1-e^2}\sin\nu}{1+e\cos\nu}$$ Pero aún $$x=\frac{a(1-e^2)\cos\nu}{1+e\cos\nu}$$. Ahora pudo obtener el área de la región azul porque el sector $CPQ^{\prime}$ del círculo centrado en $C$, en $-ae$ en el eje $x$ y radio $a$ tiene un área $$\frac12a^2E$$ donde $E$ es la anomalía excéntrica y el triángulo $COQ^{\prime}$ tiene un área $$\frac12ae\cdot\frac{a\sqrt{1-e^2}\sin\nu}{1+e\cos\nu}=\frac12a^2e\sin E$$, por lo que el área del sector azul $OPQ^{\prime}$ es $$\frac12a^2(E-e\sin E)$$ y luego podemos volver y encontrar el área del sector $OPQ$ de la elipse original como $$\frac12a^2\sqrt{1-e^2}(E-e\sin E)$$. Así que si hacer una integral con un factor de $\frac1{1+e\cos\nu}$ a través de la anomalía excéntrica fue lo suficientemente bueno para Kepler, seguramente es lo suficientemente bueno para nosotros.
En el entero original, $$\int\frac{dx}{a+b\cos x}=\frac1a\int\frac{dx}{1+\frac ba\cos x}=\frac1a\int\frac{d\nu}{1+\left|\frac ba\right|\cos\nu}$$ donde $\nu=x+\pi$ (si $ab<0$) o $\nu=x$ (si $ab>0$). Ahora vemos que $e=\left|\frac ba\right|$ y podemos usar la anomalía excéntrica, $$\sin E=\frac{\sqrt{1-e^2}\sin\nu}{1+e\cos\nu}$$ $$\cos E=\frac{\cos\nu+e}{1+e\cos\nu}$$ $$d E=\frac{\sqrt{1-e^2}}{1+e\cos\nu}d\nu$$ y la integral es $$\begin{align}\int\frac{dx}{a+b\cos x}&=\frac1a\int\frac{d\nu}{1+e\cos\nu}=\frac12\frac1{\sqrt{1-e^2}}\int dE\\ &=\frac1a\frac1{\sqrt{1-e^2}}E+C=\frac{\text{sgn}\,a}{\sqrt{a^2-b^2}}\sin^{-1}\left(\frac{\sqrt{a^2-b^2}\sin\nu}{|a|+|b|\cos\nu}\right)+C\\&=\frac{1}{\sqrt{a^2-b^2}}\sin^{-1}\left(\frac{\sqrt{a^2-b^2}\sin x}{a+b\cos x}\right)+C\end{align}$$ Por supuesto, es otra historia si $\left|\frac ba\right|\ge1$, donde obtenemos una órbita no ligada, pero esa es una historia para otra hora de dormir.
2 votos
El nombre "sustitución de Weierstrass" es desafortunado, ya que Weierstrass no tuvo nada que ver con ello (a pesar del libro de cálculo de Stewart). $\qquad$