viernes, 13 de octubre de 2023

(887) - Trigonometría elíptica de Jacobi

¿Alguien se ha preguntado el lector por qué cuando se ve el péndulo simple en bachillerato, siempre se aproxima, nunca dando su solución exacta, aun siendo un problema ideal? Esto a veces ocurre hasta en los cursos inferiores de universidad. El problema no es que no tenga una solución analítica cerrada, sino que la introducción matemática previa para poder comprenderlo, es demasiado a veces.

En física hay ciertas ecuaciones diferenciales cuyas soluciones hacen necesario emplear una familia de funciones: las funciones elípticas de Jacobi. Algunos de estos casos son el péndulo simple (como ya hemos comentado), el oscilador de Duffing, o la solución de la I Ley de Kepler en Relatividad general.

Esta familia de funciones puede ser muy laboriosa y engorrosa de trabajar y lidiar con ellas, ya que son funciones univariables que se definen en función de un parámetro a través de integrales. Algunas se definen simplemente como la inversa de otra, añadiendo otro grado de dificultad en algunos puntos.

Ecuación diferencial y solución para el ángulo de un péndulo simple: $$ \frac{\mathrm{d}^2}{\mathrm{d}t^2} \theta(t) + {\omega_0}^2 \sin\!\big(\theta(t)\big) = 0 \implies \theta(t) = 2\operatorname{am}\!\left( \frac{\sqrt{2+c_1\,}}{2}(\omega_0t+c_2) \Big| \frac{4}{2+c_1} \right) $$ Las funciones elípticas de Jacobi se definen como un conjunto de funciones integrales paramétricas: Veamos algunos ejemplos de cómo varían según este parámetro: $$ \lim_{m\to 0} \operatorname{sn}(u|m) = \sin(u) \qquad \lim_{m\to 1} \operatorname{sn}(u|m) = \operatorname{tgh}(u) $$ $$ \lim_{m\to 0} \operatorname{cn}(u|m) = \cos(u) \qquad \lim_{m\to 1} \operatorname{cn}(u|m) = \frac{1}{\cosh(u)} $$ $$ \lim_{m\to 0} \operatorname{dn}(u|m) = 1 \qquad \lim_{m\to 1} \operatorname{dn}(u|m) = \frac{1}{\cosh(u)} $$ $$ \lim_{m\to 0} \operatorname{am}(u|m) = u \qquad \lim_{m\to 1} \operatorname{am}(u|m) = \operatorname{gd}(u) $$ Donde $\operatorname{gd}(u)$ es la función gundermaniana, definida como $$ \operatorname{gd}(u) \overset{\mathrm{def}}{=} \int_0^u \frac{1}{\cosh(t)}\;\mathrm{d}t = \arctan\!\big(\sinh(u)\big) $$


Autor: Đɑvɪẟ Ƒernández-De la Cruʒ.