Este es el enfoque loco ... Vamos a encontrar el ecuación ¡de la curva de la Flor!
La curva es la sobre de segmentos de línea determinados por los dos puntos del planeta, que podemos parametrizar como $$P = p \,( \cos 8 t, \sin 8 t ) \qquad Q = q\, ( \cos 13 t, \sin 13 t )$$ donde $p$ y $q$ son los radios de las órbitas. La línea $PQ$ es entonces $$f(t) := x \left( p \sin 8 t - q \sin 13 t \right) + y ( - p \cos 8 t + q \cos 13 t ) + p q \sin 5 t \tag{1}$$ para que $$f^\prime(t) := x \left( 8 p \cos 8 t - 13 q \cos 13 t \right) + y \left( 8 p \sin 8 t - 13 q \sin 13 t \right) + 5 p q \cos 5 t \tag{2}$$ Para derivar la ecuación de la envolvente, "todo lo que tenemos que hacer" es eliminar $t$ de $f(t)$ y $f^\prime(t)$ . Con muchas envolventes clásicas (cardioides, astroides, parábolas, etc.), el paso de eliminación es bastante sencillo. Aquí, con $t$ atrapados dentro de los distintos senos y cosenos, las cosas se complican un poco. Ayuda un poco reescribir las ecuaciones en términos de la exponencial compleja, $\omega := \exp(it)$ ; ayudará a definir aún más $z := x + i y$ y $\bar{z} := x - i y$ . Entonces, hacemos las sustituciones $$\cos kt = \frac12(\omega^k+\omega^{-k}) \qquad \sin kt = \frac1{2i}(\omega^k-\omega^{-k}) \qquad x = \frac12(z+\bar{z}) \qquad y = \frac1{2i}(z-\bar{z})$$
Entonces $(1)$ y $(2)$ se convierten, después de despejar los denominadores, $$\omega^{26} q \bar{z} - \omega^{21} p \bar{z} - \omega^{18} p q + \omega^8 p q + \omega^5 p z - q z \tag{3}$$ $$13 \omega^{26} q \bar{z} - 8 \omega^{21} p \bar{z} - 5 \omega^{18} p q - 5 \omega^8 p q - 8 \omega^5 p z + 13 q z \tag{4}$$ y ahora "todo lo que tenemos que hacer" es eliminar $\omega$ . De nuevo, esto es conceptualmente es sencillo: nos limitamos a invocar el método de resultantes que Mathematica implementa gentilmente como su Resultant[]
función. Sin embargo, la complejidad computacional de la búsqueda de la resultante está relacionada con la suma de los grados de los polinomios ---aquí 52--- lo que provoca Resultant[]
se atasque en mi ordenador. Felizmente, pude conseguir una implementación más rápida al tener Mathematica construir el Matriz Sylvester para los polinomios, y tomar su determinante (que equivale a la resultante). El resultado(ant) es :
$$\begin{align} &2814749767106560000000000 p^{36} q^{42} \\ + &\cdots \;\text{($760$ terms)}\; \cdots \\ + &18258084432456195379316013815625 p^{26} q^{16} z^{13} \bar{z}^{23} \end{align}\tag{5}$$
Ahora, en lugar de restaurar $z$ y $\bar{z}$ volver a $x$ y $y$ a la forma polar compleja más adecuada: $$z = r \exp(i\theta) \qquad \bar{z} = r\exp(-i\theta)$$ para conseguir $$\begin{align} &15879378388503914086400000 \;\left(e^{15i\theta}+e^{-15i\theta} \right)\; p^{39} q^{24} r^{15} \\ +&\phantom{1}3 \left(e^{10i\theta}+e^{-10i\theta} \right)\;p^{26} q^{16} r^{10} (3289829448089600000000000 p^{12} q^{14} + \cdots) \\ -&12 \left(e^{5i\theta}+e^{-5i\theta} \right)\; p^{13} q^8 r^5 (1179648000000000000000 p^{26} q^{26} + \cdots) \\ -&\phantom{1}2\;(1407374883553280000000000 p^{36} q^{42} + \cdots) \\ \end{align} \tag{6}$$ Convenientemente, las exponenciales complejas se emparejan para formar cosenos (con un factor de $2$ ), por lo que la ecuación polar de la Flor de Venus tiene la forma
$$\begin{align} 0 \;=\;\; &31758756777007828172800000 \;\cos 15\theta\; p^{39} q^{24} r^{15} \\ +&\phantom{1}6 \cos 10\theta \;p^{26} q^{16} r^{10} (3289829448089600000000000 p^{12} q^{14} + \cdots) \\ -&24 \cos 5\theta \; p^{13} q^8 r^5 (1179648000000000000000 p^{26} q^{26} + \cdots) \\ -&\phantom{1}2\;(1407374883553280000000000 p^{36} q^{42} + \cdots) \\ \end{align}\tag{$\star$}$$
Como los argumentos de los cosenos son todos múltiplos de $5\theta$ la curva debe tener $5$ -simetría rotacional doble. Fácil, fácil. $\square$
Por supuesto, esto ad hoc enfoque no proporciona ninguna idea de por qué el $5$ -aparece la doble simetría, lo que ocurre con otras "diferencias de proporción", etc. En realidad sólo quería comprobar por mí mismo lo disparatada que es la ecuación de la curva de Flor. Sospecho que debe haber una forma de anticiparse a la simetría de $(1)$ y $(2)$ (o $(3)$ y $(4)$ ) sin tener que calcular realmente la resultante; si es así, entonces (presumiblemente) esa estrategia debería aplicarse a las diferencias de proporción arbitrarias.
0 votos
Probablemente no sea el mejor lugar para preguntar, pero ¿podrías compartir cómo hiciste eso en geogebra? tengo curiosidad
1 votos
Por supuesto :) He introducido dos números complejos de la forma $r cis \theta$ . Luego creé las circunferencias para las trayectorias de los dos números complejos para que sus trayectorias fueran visibles en la animación. Luego creé un segmento entre los dos números, activé el trazado, creé algunos deslizadores, activé la animación y ¡voilá!