Las leyes de Kepler fueron enunciadas por Johannes Kepler para describir matemáticamente el movimiento de los planetas en sus órbitas alrededor del Sol.^[1]​

Primera ley (1609)

Todos los planetas se desplazan alrededor del Sol describiendo órbitas elípticas. El Sol se encuentra en uno de los focos de la elipse.

Segunda ley (1609)

El radio vector que une un planeta y el Sol recorre áreas iguales en tiempos iguales.

La ley de las áreas es equivalente a la constancia del momento angular, es decir, cuando el planeta está más alejado del Sol (afelio) su velocidad es menor que cuando está más cercano al Sol (perihelio).

El afelio y el perihelio son los dos únicos puntos de la órbita en los que el radio vector y la velocidad son perpendiculares. Por ello solo en esos 2 puntos el módulo del momento angular ${\displaystyle L}$ se puede calcular directamente como el producto de la masa del planeta por su velocidad y su distancia al centro del Sol.

${\displaystyle L=m\cdot r_{a}\cdot v_{a}=m\cdot r_{p}\cdot v_{p}\,}$

En cualquier otro punto de la órbita distinto del Afelio o del Perihelio el cálculo del momento angular es más complicado, pues como la velocidad no es perpendicular al radio vector, hay que utilizar el producto vectorial.

${\displaystyle \mathbf {L} =m\cdot \mathbf {r} \times \mathbf {v} \,}$

Tercera ley (1619)

Para cualquier planeta, el cuadrado de su período orbital es directamente proporcional al cubo de la longitud del semieje mayor de su órbita elíptica.

${\displaystyle {\frac {T^{2}}{a^{3}}}=C={\text{constante}}}$

Donde, T es el período orbital (tiempo que tarda en dar una vuelta alrededor del Sol), a  la distancia media del planeta con el Sol y C la constante de proporcionalidad.

Estas leyes se aplican a otros cuerpos astronómicos que se encuentran en mutua influencia gravitatoria, como el sistema formado por la Tierra y el Sol.

Antes de que se redactaran las leyes de Kepler hubo otros científicos como Claudio Ptolomeo, Nicolás Copérnico y Tycho Brahe cuyas principales contribuciones al avance de la ciencia estuvieron en haber conseguido medidas muy precisas de las posiciones de los planetas y de las estrellas. Kepler, que fue discípulo de Tycho Brahe, aprovechó todas estas mediciones para poder formular su tercera ley.

Kepler logró describir el movimiento de los planetas. Utilizó los conocimientos matemáticos de su época para encontrar relaciones entre los datos de las observaciones astronómicas obtenidas por Tycho Brahe y con ellos logró componer un modelo heliocéntrico del universo. Comenzó trabajando con el modelo tradicional del cosmos, planteando trayectorias excéntricas y movimientos en epiciclos, pero encontró que los datos de las observaciones lo situaban fuera del esquema que había establecido Copérnico, lo que lo llevó a concluir que los planetas no describían una órbita circular alrededor del Sol. Ensayó otras formas para las órbitas y encontró que los planetas describen órbitas elípticas, las cuales tienen al Sol en uno de sus focos. Analizando los datos de Brahe, Kepler también descubrió que la velocidad de los planetas no es constante,^[2]​ sino que el radio vector que une al Sol (situado en uno de los focos de la trayectoria elíptica) con un planeta determinado, describe áreas iguales en tiempos iguales. En consecuencia, la velocidad de los planetas es mayor cuando están próximos al Sol (perihelio) que cuando se mueven por las zonas más alejadas (afelio). Esto da origen a las tres Leyes de Kepler sobre el movimiento planetario.

Las leyes de Kepler representan una descripción cinemática del sistema solar.

• Primera Ley de Kepler: Todos los planetas se mueven alrededor del Sol siguiendo órbitas elípticas. El Sol está en uno de los focos de la elipse.

• Segunda Ley de Kepler: Los planetas se mueven con velocidad areolar constante. Es decir, el vector posición r de cada planeta con respecto al Sol barre áreas iguales en tiempos iguales.

Se puede demostrar que el momento angular es constante lo que nos lleva a las siguientes conclusiones:

Las órbitas son planas y estables.

Se recorren siempre en el mismo sentido.

La fuerza que mueve los planetas es central.

• Tercera Ley de Kepler: Se cumple que para todos los planetas, la razón entre el periodo de revolución al cuadrado y el semieje mayor de la elipse al cubo se mantiene constante. Esto es:

${\displaystyle {\frac {T^{2}}{a^{3}}}=C}$

El estudio de Newton de las leyes de Kepler condujo a su formulación de la ley de la gravitación universal.

La formulación matemática de Newton de la tercera ley de Kepler para órbitas circulares es:

La fuerza gravitacional crea la aceleración centrípeta necesaria para el movimiento circular de radio a:

${\displaystyle {\frac {GMm}{a^{2}}}=m\omega ^{2}a}$

recordando la expresión que relaciona la velocidad angular y el período de revolución:

${\displaystyle \omega ={\frac {2\pi }{T}}}$

de donde se deduce que el cuadrado del tiempo de una órbita completa o periodo es:

${\displaystyle T^{2}={\frac {4\pi ^{2}}{GM}}a^{3}}$,

y despejando:

${\displaystyle {\frac {T^{2}}{a^{3}}}={\frac {4\pi ^{2}}{GM}}=C}$,

donde ${\displaystyle C}$ es la constante de Kepler, T  es el periodo orbital, a  el semieje mayor de la órbita, M es la masa del cuerpo central y G  una constante denominada Constante de gravitación universal cuyo valor marca la intensidad de la interacción gravitatoria y el sistema de unidades a utilizar para las otras variables de esta expresión. Esta expresión es válida tanto para órbitas circulares como elípticas.

En realidad, esta última expresión es solo una aproximación de la expresión más general que se deduce con todo rigor de las Leyes de Newton y que es:

${\displaystyle {\frac {T^{2}}{a^{3}}}\ (M+m)={\frac {4\pi ^{2}}{G}}}$

Donde ${\displaystyle M}$ es la masa del cuerpo central, ${\displaystyle m}$ la del astro que gira en torno a él y ${\displaystyle a}$ sería el semieje mayor con respecto al centro de masas del sistema. Como en el sistema solar la masa del Sol es muy superior a la de cualquier planeta, ${\displaystyle m\ll M}$, la expresión simplificada se obtiene de la más general haciendo ${\displaystyle M+m\simeq M}$

La demostración de la primera y segunda ley de Kepler se basa en las leyes de Newton y la ley de gravitación universal.

En primer lugar, se fija un sistema de referencia de coordenadas polares:

${\textstyle {\overrightarrow {x}}(t)=r(t)\cdot (cos(\theta (t)),sen(\theta (t)))}$, ${\displaystyle {\overrightarrow {u_{r}}}(t)=(cos(\theta (t)),sen(\theta (t)))}$, ${\displaystyle {\overrightarrow {u_{\theta }}}(t)=(-sen(\theta (t)),cos(\theta (t)))}$,

donde ${\displaystyle {\overrightarrow {x}}(t)}$ denota la posición del cuerpo con masa ${\displaystyle m}$ en el instante ${\displaystyle t}$; el cuerpo con masa ${\displaystyle M}$ está quieto y en el origen; y ${\displaystyle {\overrightarrow {u_{r}}}}$ y ${\displaystyle {\overrightarrow {u_{\theta }}}}$ son vectores unitarios en las direcciones radial y circunferencial, respectivamente; y ${\displaystyle \theta (t)}$ es el ángulo que forma ${\displaystyle r(t)}$ con el eje polar (eje de referencia desde el que se mide el ángulo polar).

${\displaystyle {\overrightarrow {u_{r}}}}$ y ${\displaystyle {\overrightarrow {u_{\theta }}}}$ satisfacen las siguientes propiedades:

${\displaystyle {d{\overrightarrow {u_{r}}} \over d\theta }={\overrightarrow {u_{\theta }}}}$; ${\displaystyle \quad {d{\overrightarrow {u_{r}}} \over dt}={\overrightarrow {u_{r}}}'={\overrightarrow {u_{\theta }}}\theta '}$; ${\displaystyle \quad {\overrightarrow {u_{r}}}\bot {\overrightarrow {u_{\theta }}}}$.

La fuerza ${\displaystyle {\overrightarrow {F}}}$sobre el cuerpo de masa ${\displaystyle m}$ se descompone en: ${\displaystyle {\overrightarrow {F}}=F_{r}{\overrightarrow {u_{r}}}+F_{\theta }{\overrightarrow {u_{\theta }}}}$. Además, como ${\displaystyle {\overrightarrow {F}}}$ es una fuerza central, ${\displaystyle F_{\theta }\equiv 0}$.

Por lo tanto, aplicando la segunda ley de Newton,

${\displaystyle F_{r}{\overrightarrow {u_{r}}}={\overrightarrow {F}}=m{\overrightarrow {x}}''(t)}$. ${\displaystyle [1]}$

La velocidad del planeta es la derivada de la posición:

${\displaystyle x'(t)=r'(t){\overrightarrow {u_{r}}}+r(t)\theta '(t){\overrightarrow {u_{\theta }}}}$, ${\displaystyle [2]}$

y su aceleración es la derivada de la velocidad:

${\textstyle {\begin{aligned}x''(t)&=r''(t){\overrightarrow {u_{r}}}+r'(t)\theta '(t){\overrightarrow {u_{\theta }}}+r'(t)\theta '(t){\overrightarrow {u_{\theta }}}+r(t)\theta ''(t){\overrightarrow {u_{\theta }}}-r(t)\theta '(t)\theta '(t){\overrightarrow {u_{r}}}\\&=(r''(t)-r(t)(\theta '(t))^{2}){\overrightarrow {u_{r}}}+(2r'(t)\theta '(t)+r(t)\theta ''(t)){\overrightarrow {u_{\theta }}}.\\\end{aligned}}}$ ${\displaystyle [3]}$

Usando ${\displaystyle [1]}$ y ${\displaystyle [3]}$:

${\displaystyle {\begin{cases}F_{r}=[r''(t)-r(t)(\theta '(t))^{2}]m\\0=2r'(t)\theta '(t)+r(t)\theta ''(t)\end{cases}}}$ ${\displaystyle {\begin{matrix}[4]\\{[5]}\end{matrix}}}$

Multiplicando por ${\displaystyle r}$ a ambos lados de ${\displaystyle [5]}$:

${\displaystyle 0=2r(t)r'(t)\theta '(t)+r^{2}(t)\theta ''(t)=(r^{2}(t))'\theta '(t)+r^{2}(t)\theta ''(t)=(r^{2}(t)\theta '(t))'}$.

Así que ${\displaystyle r^{2}(t)\theta '(t)=c}$ (constante). ${\displaystyle [6]}$

Por otra parte, sea ${\displaystyle A_{\theta _{1},\theta _{2}}}$ el área del sector barrido entre los ángulos ${\displaystyle \theta _{1}}$ y ${\displaystyle \theta _{2}}$:

${\displaystyle A_{\theta _{1},\theta _{2}}={1 \over 2}\textstyle \int \limits _{\theta _{1}}^{\theta _{2}}\displaystyle r^{2}(\theta )d\theta }$.

Tomemos ${\displaystyle \theta _{1}=0}$ por simplicidad y denotemos ${\displaystyle A_{\theta _{1},\theta _{2}}=A}$.

Por el teorema fundamental del cálculo, ${\displaystyle {dA \over d\theta }={r^{2}(\theta ) \over 2}}$. Como ${\displaystyle \theta }$ es función de ${\displaystyle t}$, por ${\displaystyle [6]}$:

${\displaystyle A'(t)={dA(\theta (t)) \over d\theta }\theta '(t)={r^{2}(\theta (t)) \over 2}\theta '(t)=c}$ (constante).

Por lo tanto, ${\displaystyle A(t)=ct+k}$, para alguna constante ${\displaystyle k}$.

Se obtiene: ${\displaystyle A(t)=ct}$.

Aplicando esto a dos intervalos de tiempo de igual longitud, ${\displaystyle [t_{1},t_{2}]}$ y ${\displaystyle [t_{3},t_{4}]}$:

${\displaystyle A(t_{2}-t_{1})=(t_{2}-t_{1})c=(t_{4}-t_{3})c=A(t_{4}-t_{3})}$. ■

Imponiendo que ${\displaystyle F_{r}}$ cumpla la ley universal de gravitación:

${\displaystyle F_{r}={-{GmM \over r^{2}(t)}}}$. ${\displaystyle [7]}$

Aplicando la segunda ley de Newton y la ley de gravitación universal:

${\displaystyle {\overrightarrow {x''}}(t)=-{\Bigl (}G{M \over r^{2}(t)}{\Bigr )}{\overrightarrow {u_{r}}}}$. ${\displaystyle [8]}$

Igualando ${\displaystyle [4]}$ y ${\displaystyle [7]}$:

${\displaystyle -{GM \over r^{2}(t)}=r''(t)-r(t)(\theta '(t))^{2}}$ . ${\displaystyle [9]}$

Despejando de la ecuación ${\displaystyle [6]}$ se obtiene: ${\displaystyle \theta '(t)={c \over r^{2}(t)}}$, para c constante. ${\displaystyle [10]}$

Se puede reescribir la ecuación ${\displaystyle [9]}$ , usando ${\displaystyle [10]}$, como:

${\displaystyle -{GM \over r^{2}(t)}=r''(t)-{c^{2} \over r^{3}(t)}}$. ${\displaystyle [11]}$

Haciendo el cambio ${\displaystyle z(t)={1 \over r(t)}}$ y derivando dos veces, obtenemos lo siguiente:

${\displaystyle r'(t)=-{1 \over z^{2}(t)}z'(t)=-{1 \over z^{2}(t)}{dz(t) \over d\theta (t)}\theta '(t)=-c{dz(t) \over d\theta (t)}}$,

${\displaystyle r''(t)=-c{d \over dt}{\biggl (}{dz(t) \over d\theta (t)}{\Biggr )}=-c{d^{2}z(t) \over d\theta ^{2}(t)}\theta '(t)}$.

Usando ${\displaystyle [10]}$:

${\displaystyle r''(t)=-c^{2}z^{2}(t){d^{2}z(t) \over d\theta ^{2}(t)}}$.

Sustituyendo ${\displaystyle r''(t)}$ en el lado derecho de ${\displaystyle [11]}$:

${\displaystyle r''(t)-{c^{2} \over r^{3}(t)}=-c^{2}z^{2}(t){d^{2}z(t) \over d\theta ^{2}(t)}-{c^{2} \over r^{3}(t)}=-c^{2}z^{2}(t){d^{2}z(t) \over d\theta ^{2}(t)}-c^{2}z^{3}(t)}$,

y en el lado izquierdo aplicando de nuevo el cambio ${\displaystyle z(t)={1 \over r(t)}}$:

${\displaystyle -{GM \over r^{2}(t)}=-GMz^{2}(t)}$.

De esta forma, la ecuación ${\displaystyle [9]}$ se puede escribir como:

${\displaystyle z+{d^{2}z \over d\theta ^{2}}={GM \over c^{2}}}$.

Esta ecuación diferencial tiene como únicas soluciones:

${\displaystyle z(\theta )=Acos(\theta -\alpha )+{GM \over c^{2}}}$, donde ${\displaystyle A}$ y ${\displaystyle \alpha }$ son constantes.

Eligiendo el eje polar de manera que ${\displaystyle \alpha =0}$:

${\displaystyle r(\theta )={1 \over Acos(\theta )+{GM \over c^{2}}}={{c^{2} \over GM} \over Bcos(\theta )+1}}$, donde ${\displaystyle B={Ac^{2} \over GM}}$.

Haciendo los cambios ${\displaystyle e=B}$ y ${\displaystyle p={1 \over A}}$, se obtiene la ecuación de una cónica con foco en el origen:

${\displaystyle r(\theta )={ep \over ecos(\theta )+1}}$, donde ${\displaystyle e}$ es la excentricidad y ${\displaystyle p}$ es la distancia del foco a la directriz.

Según el valor de ${\displaystyle e}$, esta cónica puede ser una elipse, una hipérbola o una parábola.

En este caso, si la trayectoria de los cuerpos celestes está acotada, el único caso posible es que sea una elipse, ${\displaystyle 0<e<1}$. ■

Johannes Kepler descubrió sus leyes gracias a un considerable trabajo de análisis de las observaciones astronómicas realizadas por Tycho Brahe, mucho más precisas que las ya conocidas; se basó en particular en las posiciones de Marte, cuyo movimiento estudió a partir de 1600. Estaba convencido de que el Sol era de algún modo el "verdadero" centro del sistema solar (para los planetas exteriores, como Marte, Copérnico utilizaba un punto ficticio cercano al Sol como centro de un círculo sobre el que giraba a velocidad uniforme el centro de un pequeño epiciclo que llevaba el planeta). Guiado por esta creencia y tras mucho divagar, acabó descubriendo que el movimiento de los planetas es elíptico, con el sol situado en un foco de la elipse. Sus resultados y el modo en que llegó a ellos se recogen en su obra principal, la Astronomia nova, que apareció en 1609, pero que en realidad se terminó a finales de 1605.^[3]​.

Sus propias leyes permitieron afinar la investigación astronómica y revelar las irregularidades de los movimientos de los cuerpos conocidos mediante una asombrosa progresión de los análisis.

El ejemplo más espectacular fue el de las irregularidades de Urano, que condujeron al descubrimiento de Neptuno por John Couch Adams (1819 - 1892) y Urbain Le Verrier (1811 - 1877), mediante cálculo: descubrimiento confirmado por la observación de Johann Gottfried Galle (1812 - 1910) en 1846.

• SIMMONS, G. F., & KRANTZ, S. G. ECUACIONES DIFERENCIALES: TEORÍA, TÉCNICA Y PRÁCTICA/GEORGE F. SIMMONS Y STEVEN G. KRANTZ (No. QA371. S46 2007.).
• Perez, I. I. C. (2015). Leyes de Kepler.

• Movimiento circular
• Gravedad
• Órbita de Kepler
• Problema de Kepler
• Ecuación de Kepler
• Vector de Runge-Lenz
• Momento angular relativo específico
• Ecuación del tiempo
• Problema de los dos cuerpos
• Revolución de Copérnico

• Crowell, Benjamin, Conservation Laws, http://www.lightandmatter.com/area1book2.html Archivado el 1 de junio de 2020 en Wayback Machine., un libro electrónico que ofrece una prueba de la primera ley sin el uso del cálculo. (section 5.2, p. 112)
• David McNamara and Gianfranco Vidali, Kepler's Second Law - Java Interactive Tutorial, http://www.phy.syr.edu/courses/java/mc_html/kepler.html, an interactive Java applet that aids in the understanding of Kepler's Second Law.
• Audio - Cain/Gay (2010) Astronomy Cast Johannes Kepler and His Laws of Planetary Motion
• University of Tennessee's Dept. Physics & Astronomy: Astronomy 161 page on Johannes Kepler: The Laws of Planetary Motion [1]
• Equant compared to Kepler: interactive model [2] Archivado el 26 de diciembre de 2008 en Wayback Machine.
• Kepler's Third Law:interactive model [3] Archivado el 26 de diciembre de 2008 en Wayback Machine.
• Solar System Simulator (Interactive Applet)
• Kepler and His Laws, educational web pages by David P. Stern
• [4], Descripción detallada de la segunda ley de Kepler.
• [5], PUEMAC. Demostración de las leyes de Kepler a partir de las de Newton.
• [6], Demostración de las leyes de Kepler.