Ponente
Descripción
Este trabajo propone un modelo matemático dinámico de la respuesta sexual humana en parejas, fundamentado en ecuaciones diferenciales ordinarias (EDOs). Se parte de la necesidad de superar los enfoques descriptivos clásicos mediante una representación cuantitativa que capture las interacciones continuas, simultáneas y recíprocas entre dos individuos durante un encuentro sexual. Los modelos secuenciales de Masters y Johnson (Masters, W.H. & Johnson, V.E., Human sexual response, Little, Brown and Company, 1966) y Kaplan (Kaplan, H.S., The new sex therapy, Brunner/Mazel, 1979) segmentan la respuesta en fases lineales (excitación, meseta, orgasmo y resolución), pero no representan la naturaleza bidireccional y retroalimentada de la interacción real, donde los cambios en un individuo afectan dinámicamente al otro.
El estado del arte revela que los modelos descriptivos han sentado las bases de la sexología clínica, pero no abordan explícitamente las dinámicas simultáneas ni los bucles de retroalimentación presentes en la interacción real. Recientemente, Blyuss y Kyrychko (Blyuss, K.B. & Kyrychko, Y.N., Chaos 33, 043136, 2023) desarrollaron un modelo basado en un sistema excitable acoplado inspirado en FitzHugh-Nagumo, que representa la excitación fisiológica u y psicológica v de un individuo mediante ecuaciones no lineales con función tipo N invertida. El sistema unipersonal de Blyuss y Kyrychko (2023) se describe mediante u˙=f(u)−v+Eu y v˙=ϕ[au−bv+(Ev−E0v)] , donde f(u)=u−u3 . Este avance motivó la presente extensión hacia un sistema bipersonal.
La propuesta introduce variables de excitación fisiológica y psicológica para cada individuo, (u1,v1) y (u2,v2) , e incorpora términos de acoplamiento que modelan la atracción mutua y la sincronización entre ambos. El sistema resultante consta de cuatro EDOs no lineales acopladas:
\begin{equation} \dot{u}_1 = f(u_1) - v_1 + E_1^u + \alpha(u_2 - u_1) \end{equation}
\begin{equation} \dot{u}_2 = f(u_2) - v_2 + E_2^u + \alpha(u_1 - u_2) \end{equation}
\begin{equation} \dot{v}_1 = \phi [a u_1 - b_1 v_1 + E_1^v] \end{equation}
\begin{equation} \dot{v}_2 = \phi [a u_2 - b_2 v_2 + E_2^v] \end{equation}
Los parámetros ajustan la intensidad de la interacción, la estimulación externa física y mental, y la retroalimentación entre los miembros de la pareja. El coeficiente de acoplamiento α representa específicamente la atracción entre las excitaciones fisiológicas de ambos individuos.
Los objetivos principales son: desarrollar el modelo dinámico bipersonal integrando componentes fisiológicos y psicológicos; simular diversos escenarios mediante la variación de parámetros clave como atracción mutua, estimulación externa y condiciones iniciales; y evaluar si el sistema reproduce comportamientos esperados como picos de excitación tipo orgasmo y fenómenos de sincronización o desincronización de la pareja.
Los resultados de 1000 simulaciones aleatorias y cinco escenarios específicos demuestran que los orgasmos dobles sincronizados son poco frecuentes bajo estimulación baja o asimétrica. Mientras ambos individuos reciban alta estimulación externa y el coeficiente α es elevado, la dinámica cambia significativamente: la excitación se acumula en paralelo y emergen picos simultáneos, favoreciendo la sincronización. Esto sugiere que la coincidencia de entradas altas y un acoplamiento fuerte constituye una condición necesaria para el orgasmo doble sincronizado.
La modelación permite integrar conocimientos fisiológicos y psicológicos dispersos y realizar experimentos en escenarios difíciles de ensayar en la realidad.