¿Todos los puntos de equilibrio de un mapeo discreto aparecen en el diagrama de bifurcación?

La pregunta en el título quizás esté vagamente planteada, así que incluiré el ejemplo concreto que me está molestando.

Supongamos que tenemos un mapeo dado por

norte t + 1 = norte t Exp ( r ( 1 norte t PAG norte t / ( α 2 + norte t 2 ) ) ) ,
dónde α y r son algunas constantes fijas. Si trazamos este mapeo no lineal, obtenemos algo como esto (trazado para r = 0.35 , α = 0.1 y, de arriba abajo, PAG 1 = 0.165 , PAG 2 = 0.235 , PAG 3 = 0.35 ). La línea azul corresponde a norte t + 1 = norte t

Gráfico 1

A partir de aquí se puede ver que hay 3 regímenes en los que puede estar el sistema, dependiendo del valor del parámetro PAG :

  • Para valores bajos de PAG , el sistema tiene 2 puntos fijos (mostrado en 0 y en cca. 0.8)
  • Para valores medios de PAG , el sistema tiene 4 puntos fijos (mostrado en 0, A , B y C )
  • Para valores altos de PAG , el sistema tiene 2 puntos fijos nuevamente (mostrado en 0 y en cca. 0.03).

Me tomé la libertad de dibujar diagramas de telaraña para todos ellos (excepto norte = 0 , que parece ser un punto fijo inestable para los valores dados), pero mostraré solo los importantes para mi pregunta:

  1. Caso PAG = 0.235 , punto de partida norte 0 < norte B , 2048 iteraciones muestran que es un punto de órbita C .

Gráfico 2

  1. Mismo caso, esta vez el punto de partida es norte 0 > norte B por lo que el sistema tiende hacia el punto A , pero mucho más rápido.

Gráfico 3

Esto parecía extraño al principio, así que decidí trazar el diagrama de bifurcación y el exponente de Lyapunov de los valores de PAG :

Diagrama de bifurcación (los puntos son bastante densos después de PAG = 0.4 , por lo que no es tan confiable para la comparación con el exponente de Lyapunov después de ese punto):

Gráfico 4

Exponente de Lyapunov:

Gráfico 5

La mayoría de las conclusiones que hice antes de dibujar esos 2 últimos diagramas se confirmaron, pero queda un problema: en el medio PAG régimen ( PAG = 0.235 ), el diagrama de bifurcación muestra que el sistema tiene solo un estado estable (el más alto, punto A ), pero se niega a reconocer al otro ( C ), a pesar de que el diagrama de telaraña muestra que probablemente sea estable. No solo eso, el exponente de Lyapunov está alrededor de -2 para ese punto, y estoy bastante seguro de que el gráfico es preciso (está trazado para los primeros 1000 términos en la suma de Lyapunov, con un paso discreto de 1/10000).

Entonces, mi problema principal es: ¿El diagrama de bifurcación muestra todos los puntos de equilibrio a los que se puede acceder desde una multitud de puntos de partida para un valor dado del parámetro de bifurcación? En este caso, si el sistema aterrizó en A o C dependía únicamente de la condición inicial, a saber norte 0 < norte B o norte 0 > norte B .

Para los regímenes de período doble y caótico obtuve los diagramas de telaraña correctos (el sistema salta entre 2 puntos distintos y no se asienta en ningún punto, respectivamente), entre otras predicciones correctas, por lo que estoy convencido de que el error no está en mi codificación. Cualquier idea en este asunto sería muy apreciada.

EDITAR: siéntase libre de migrar esto a math.se si no lo encuentra apropiado. La pregunta surge de un problema de física, pero entiendo si se considera fuera de tema

EDIT#2: Como se sugirió, aquí está la versión ampliada del diagrama de bifurcación en el intervalo de PAG = 0.11 a PAG = 0.5

bifurcación agrandada

¿Qué tipo de software o algoritmo utiliza para dibujar el diagrama de bifurcación? Esto parece ser un problema numérico/de interpretación. También sería interesante ver el diagrama de bifurcación solo desde PAG = 0.25 porque la parte caótica de Low P se roba el espectáculo.
Sí, es un problema de interpretación, aunque es posible que no entienda lo suficientemente bien el tema. Dibujé todos los gráficos (incluido el diagrama de bifurcación) en Sage, un paquete de programación simbólica de código abierto similar a Python. Si es necesario, puedo proporcionar el código, pero no creo que sea necesario. También agregaré el diagrama de bifurcación para un rango más pequeño a la pregunta, gracias por señalarlo.

Respuestas (1)

Vale, no conozco el algoritmo de Sage pero voy a ofrecer una conjetura de lo que está pasando. Tienes que verificar la conjetura mediante más investigaciones numéricas. Supongo que el algoritmo de Sage funciona de manera óptima para las bifurcaciones de un solo equilibrio y puede encontrarse con problemas como los que vemos aquí cuando se trata de equilibrios (también conocidos como puntos fijos) asociados con más de un punto.

Lo que ves en tu diagrama de bifurcación para PAG > 0.26 es en realidad una oscilación de 2 puntos alrededor del equilibrio inestable en C (un algoritmo de bifurcación básico no mostrará puntos fijos inestables). La razón por la cual esta oscilación parece converger en C es probablemente porque la oscilación siempre está muy cerca de C .

En PAG 0.26 la curva del mapa norte t + 1 ( norte t ) toca el norte t + 1 = norte t línea, los equilibrios A y B surgir como una "bifurcación tangente" para PAG < 0.26 y esto confunde el algoritmo de Sage para cambiar de la oscilación alrededor C a A . Esta confusión se muestra como la línea vertical en PAG 0.26 en tu diagrama. (Alternativamente, podría haber una breve ventana de caos). Una vez más, el diagrama de bifurcación no puede mostrar el equilibrio inestable B que se mostraría como una "rama inferior" en la bifurcación tangente como se ve para PAG < 0.26 . Para verificar esta conjetura, puedes comprobar que en PAG 0.26 emergen tanto A como B.

Entonces, el problema probablemente sea que el algoritmo de Sage está tratando de ahorrar tiempo de cálculo al buscar solo en la vecindad de los puntos anteriores en el diagrama de bifurcación, pero al hacerlo podría omitir partes del diagrama que se deben por completo a un equilibrio diferente. Para solucionar este problema, simplemente debe configurar Sage de una manera diferente o usar un software diferente. (No es difícil escribir un programa ineficiente pero casi a prueba de fallas que dibuje el diagrama simplemente disparando un montón de condiciones iniciales aleatorias densamente colocadas y simplemente mostrando dónde se establecen después de una gran cantidad de iteraciones para cada PAG .)

Buena respuesta. Sus primeros dos párrafos son algunas de las conclusiones que ya hice, pero no las incluí en la pregunta, ya que quería reducir el tema. Aún así, es tranquilizador ver pensamientos similares de otra persona. En cuanto a la respuesta: actualmente no tengo tiempo, pero a primera hora de la mañana estoy revisando el algoritmo de trazado. No estoy realmente convencido de que el problema esté ahí (mi intuición me dice que es un problema conceptual, no un error de trazado que se origina en la llegada discontinua de un nuevo equilibrio), pero vale la pena intentarlo. Por ahora, +1 de mi parte.
Después de unos años, volví a esta pregunta y concluí que su respuesta era correcta: después de experimentar un poco, resultó que el punto C es en realidad un punto fijo inestable (con una órbita estable a su alrededor, pero que no se muestra en el gráfico debido al algoritmo simple).