La PizarraMatemática.

Máster y doctorado

Bifurcaciones y caos en sistemas no lineales

Una bifurcación aparece al variar un parámetro y cambiar el retrato cualitativo del sistema. Puede nacer una pareja de equilibrios o un ciclo. Empezaremos con x′=μ−x², donde todo se puede calcular; después distinguiremos esos cambios locales de una afirmación de caos.

1. Mover un equilibrio no siempre cambia su tipo de comportamiento

Con el parámetro fijo, estabilidad pregunta qué hacen las órbitas cercanas a un objeto. Una bifurcación compara sistemas vecinos de una familia. Cambiar solo la posición de un equilibrio puede conservar su tipo; crear o perder una pareja de equilibrios sí cambia el retrato cualitativo.

En x′=f(x,μ), si fₓ≠0 en un equilibrio, el teorema implícito lo continúa como una rama suave y el signo de fₓ persiste para cambios pequeños. Una bifurcación local de equilibrios necesita pérdida de hiperbólicidad en este marco. La pérdida por sí sola no garantiza que haya bifurcación.

El equilibrio se mueve sin perder atracción

x′=μ−x.

  1. Equilibrio

    0=μ−x da x*=μ.
  2. Derivada

    fₓ=−1 para todos los parámetros.
  3. Tipo

    La solución se acerca exponencialmente a μ; variar μ traslada el equilibrio pero no cambia este tipo local.

Resultado

No todo movimiento de una rama es una bifurcación cualitativa.

Lo comprobamos

El ejemplo tiene una derivada no nula en cada equilibrio.

Una derivada cero no obliga a un cambio de familia

x′=−x³ para todos los valores de μ.

  1. Equilibrio

    x=0 tiene derivada cero.
  2. Parámetro

    El campo ni siquiera depende de μ.

  3. Comparar

    Los sistemas vecinos son iguales; no aparece un cambio al atravesar un parámetro.

Resultado

Falta de hiperbólicidad y bifurcación no son afirmaciones equivalentes.

Lo comprobamos

Hay que revisar los términos y la dependencia real del parámetro.

2. Silla-nodo crea una pareja con estabilidad opuesta

En x′=μ−x², los equilibrios cumplen x²=μ.

Con μ>0 existen ±√μ; con μ<0 no hay ninguno; en μ=0 se unen en cero. La derivada fₓ=−2x es negativa en la rama positiva y positiva en la negativa.

La pareja de equilibrios nace en el parámetro crítico
x′=μ−x²μxestableinestableμ<0: ninguno; μ>0: dos.En μ=0: un equilibrio semiestable.
Diagrama de ramas para x prima igual a mu menos x al cuadrado. El eje horizontal es mu y el vertical es el estado x. La rama superior x igual a raíz de mu es atractora y se dibuja continua en verde; la inferior menos raíz es repulsora y se dibuja discontinua en rojo. Las ramas solo existen para mu no negativo y se unen en el origen. El dibujo es un diagrama de equilibrios al variar el parámetro, no trayectorias de una sola solución.

La rama positiva atrae localmente y la negativa repele. Una silla-nodo escalar no degenerada requiere f=0, fₓ=0, fₓₓ≠0 y fμ≠0 en el punto crítico. En el plano, añadir y′=−y produce una silla en una rama y un nodo atractivo en la otra.

Clasificar con un parámetro concreto

μ=9.

  1. Resolver

    x²=9 da x=3 y x=−3.
  2. Derivadas

    fₓ(3)=−6 y fₓ(−3)=6.
  3. Clasificar

    Tres es atractivo y menos tres repulsor localmente.

Resultado

Los signos se calculan después de hallar los equilibrios.

Lo comprobamos

El signo de la coordenada no es una definición de estabilidad: aquí decide la derivada del campo.

Comprobar el punto de bifurcación

f(x,μ)=μ−x² en (0,0).

  1. Equilibrio degenerado

    f=0 y fₓ=−2x=0.
  2. Segundo orden

    fₓₓ=−2≠0.
  3. Parámetro

    fμ=1≠0; se cumplen las condiciones de no degeneración de esta forma normal.

Resultado

En μ=0 confluyen las dos ramas.

Lo comprobamos

El dibujo de ramas usa parámetro horizontal y estado vertical; no es una trayectoria temporal.

3. En el parámetro crítico atrae por un lado y repele por el otro

Con μ=0 queda x′=−x².

El campo es negativo a ambos lados del origen: una solución positiva baja hacia cero, pero una negativa se hace más negativa. El equilibrio es semiestable, no estable bilateralmente en el sentido de Lyapunov.

La solución de dato x₀ es x(t)=x₀/(1+x₀t) mientras el denominador no se anule.

Un dato negativo puede explotar hacia menos infinito en tiempo positivo. Un campo polinómico suave garantiza existencia local, no existencia global de toda órbita. El retrato de una bifurcación local tampoco resuelve todo el comportamiento lejano.

Acercarse desde la derecha

x₀=2 en el parámetro crítico.

  1. Fórmula

    x(t)=2/(1+2t).
  2. Signo

    El denominador es positivo para todo t≥0.
  3. Límite

    x(t) decrece a cero; su decaimiento es algebraico.

Resultado

El origen atrae desde la derecha.

Lo comprobamos

No se obtiene una tasa exponencial de una derivada lineal que aquí vale cero.

Alejarse desde la izquierda

x₀=−1/2.

  1. Fórmula

    x(t)=−1/(2−t).
  2. Tiempo

    Está definida en este sentido hasta t=2.
  3. Límite

    Cuando t se acerca a dos por debajo, x(t) tiende a menos infinito.

Resultado

No todas las soluciones tienen existencia global positiva.

Lo comprobamos

Incluso datos negativos arbitrariamente cercanos a cero terminan saliendo de una vecindad fija.

4. Transcrítica y horquilla organizan ramas de maneras diferentes

En x′=μx−x², las ramas de equilibrio son cero y μ.

Sus derivadas son μ y −μ respectivamente: intercambian estabilidad al cruzar μ=0. La rama cero se conserva para todos los parámetros, una estructura que una perturbación general puede destruir.

En x′=μx−x³ hay simetría x↦−x.

Cero atrae si μ<0 y repele si μ>0; para μ>0 nacen ±√μ con derivada −2μ<0. Es la horquilla supercrítica. En μx+x³ las ramas no nulas aparecen para μ<0 y son repulsoras: el signo cúbico cambia el diagrama.

Ver el intercambio transcrítico

f=μx−x².

  1. Con μ=1

    Cero tiene derivada uno, repulsora; x=1 tiene derivada menos uno, atractora.

  2. Con μ=−1

    Cero tiene derivada menos uno, atractora; x=−1 tiene derivada uno, repulsora.

  3. Ramas

    Las curvas x=0 y x=μ se cruzan y cambian sus papeles.

Resultado

El cambio se verifica con las derivadas de cada rama.

Lo comprobamos

No se ha creado una pareja ausente a un lado como en silla-nodo.

Cambiar el signo cúbico cambia qué ramas atraen

Comparar μx−x³ y μx+x³.

  1. Supercrítica

    Con μ=4, ±2 son equilibrios y su derivada es 4−3·4=−8, atractora.
  2. Subcrítica

    Con μ=−4, ±2 son equilibrios del otro modelo y su derivada es −4+3·4=8, repulsora.

  3. Origen

    La rama cero tiene derivada μ en ambos modelos.

Resultado

No basta con dibujar tres ramas sin indicar el lado del parámetro y su estabilidad.

Lo comprobamos

La simetría de estos modelos es parte de la estructura que organiza la horquilla.

5. Hopf puede crear un ciclo cuando una pareja compleja cruza el eje imaginario

En sistemas con pocas direcciones críticas y otras hiperbólicas, una variedad central local permite reducir el estudio bajo un teorema de regularidad aplicable. No consiste en borrar términos sin justificación y no tiene necesariamente la unicidad del conjunto estable hiperbólico.

Hopf crea un ciclo y las órbitas cercanas pueden espiralar hacia él
z′=(1+2i)z−|z|²zr′=r−r³; θ′=2El círculo r=1 atrae radialmente.Periodo: π; el origen es repulsor.
Soluciones exactas en coordenadas polares del modelo complejo con mu igual a uno y omega igual a dos. La curva azul parte de radio un cuarto y ángulo menos pi sobre dos; la roja de radio tres medios y ángulo cero. Sus radios siguen uno sobre raíz de uno más r inicial elevado a menos dos menos uno por e elevado a menos dos t, y los ángulos aumentan dos t. Ambas se aproximan al círculo verde de radio uno, con periodo pi. Las curvas se dibujan en el mismo plano y escala; el origen permanece como equilibrio inestable.

En Hopf, un par complejo cruza transversalmente el eje imaginario con frecuencia no nula. Con las condiciones de no degeneración, incluida una contribución cúbica de amplitud no nula, aparece una rama de ciclos.

El modelo z′=(μ+iω)z−|z|²z, ω≠0, se escribe en polares como r′=μr−r³ y θ′=ω.

Para μ>0, el círculo r=√μ es una órbita periódica atractora radialmente.

Para μ<0 atrae el origen.

Un par imaginario sin cruce transversal no basta para esta conclusión.

Calcular radio, periodo y estabilidad radial

μ=1 y ω=2.

  1. Radio

    r′=r(1−r²) se anula en r=1; para 0<r<1 es positivo y para r>1 es negativo.
  2. Periodo

    θ′=2 da una vuelta de 2π en tiempo π.
  3. Atracción

    La derivada radial en r=1 es 1−3=−2; el ciclo atrae radialmente.
Resultado

El origen sigue siendo equilibrio, pero es repulsor y el nuevo ciclo es atractivo.

Lo comprobamos

En el gráfico las espirales tienen radios y velocidades angulares calculados con este modelo.

Un par imaginario sin cruce no garantiza un ciclo

z′=(−μ²+2i)z−|z|²z.

  1. Parte real

    En cero es −μ² y toca cero sin cambiar de signo.

  2. Radio

    r′=−μ²r−r³<0 para todo r>0.
  3. Consecuencia

    No hay radio positivo de equilibrio y no aparece un ciclo de este tipo al cruzar μ=0.

Resultado

Falta la transversalidad del cruce de Hopf.

Lo comprobamos

Mirar únicamente los autovalores imaginarios en el parámetro crítico habría sido insuficiente.

6. Un mapa discreto se clasifica por módulo y un ciclo por su retorno completo

En xₙ₊₁=Fμ(xₙ), la perturbación en un punto fijo se multiplica por Fμ′.

Atrae si el módulo es menor que uno. Puede ser negativo y acercarse alternando lados. En una ecuación diferencial escalar x′=f(x), la condición lineal atractora es en cambio fₓ<0.

Un multiplicador que cruza +1 o −1 puede producir distintas bifurcaciones bajo condiciones adicionales. El cruce de −1 puede dar duplicación de periodo. Para un ciclo de k puntos se usa (Fᵏ)′, producto de las derivadas en todo el ciclo. Los mapas de retorno transversales permiten estudiar así ciclos de flujos.

El mismo número positivo tiene dos lecturas temporales

Comparar x′=x/2 y xₙ₊₁=xₙ/2.

  1. Flujo

    La solución es eᵗᐟ²x₀ y se aleja del origen si x₀≠0.
  2. Mapa

    La solución es 2⁻ⁿx₀ y se acerca.

  3. Criterio

    En el primero importa la parte real del generador; en el segundo el módulo del multiplicador.

Resultado

Antes de clasificar hay que fijar si se itera o se integra en tiempo.

Lo comprobamos

Una derivada negativa de módulo mayor que uno puede ser repulsora en un mapa.

Un ciclo puede atraer aunque uno de sus pasos expanda

Dos puntos de un ciclo tienen derivadas −2 y −1/5.

  1. Retorno

    La derivada de F² es el producto (−2)(−1/5)=2/5.
  2. Módulo

    2/5<1.
  3. Clasificar

    El retorno de dos pasos es atractivo localmente, bajo el contexto diferenciable del ciclo.

Resultado

Se analiza la vuelta completa.

Lo comprobamos

Mirar solo el primer factor de módulo dos habría dado una conclusión equivocada sobre el ciclo.

7. El mapa logístico permite calcular la primera duplicación

Lᵣ(x)=rx(1−x) tiene puntos fijos cero y, si r≠0, 1−1/r.

En la rama no nula la derivada es 2−r. Dentro del rango físico [0,1], atrae para 1<r<3.

Al llegar a r=3 el multiplicador es −1.
Después de la primera duplicación pueden atraer dos valores alternos
L(x)=(16/5)x(1−x)01a≈0,5130 → b≈0,7995 → a(L²)′(a)=4/25: ciclo atractivo.
Gráfica exacta del mapa logístico con parámetro dieciséis quintos, junto a la diagonal identidad. El rectángulo de telaraña azul conecta los dos puntos del ciclo: a aproximadamente cero coma quinientos trece milésimos y b aproximadamente cero coma setecientos noventa y nueve milésimos. Se sube de a en la diagonal a L de a igual b, se pasa a la diagonal y se repite para volver a a. Son puntos distintos y L al cuadrado los fija. El multiplicador conjunto es cuatro veinticincoavos, de módulo menor que uno. No se trata de una nube simulada ni de una afirmación de caos para todo parámetro.

Para el nuevo ciclo, quitar las raíces fijas de Lᵣ²(x)−x deja r²x²−r(r+1)x+r+1=0.

Sus raíces son x_\pm=\frac{r+1\pm\sqrt{(r-3)(r+1)}}{2r}.

Cerca de tres por arriba forman un ciclo de dos puntos; su multiplicador conjunto es −r²+2r+4, atractivo para 3<r<1+√6.

Un punto fijo atractivo antes del umbral

r=5/2.

  1. Punto

    x*=1−2/5=3/5.
  2. Derivada

    2−5/2=−1/2.
  3. Módulo

    Vale 1/2<1; las aproximaciones pueden alternar lados mientras se acercan.

Resultado

La rama no nula es atractora.

Lo comprobamos

La derivada negativa no es por sí sola prueba de repulsión en tiempo discreto.

Un ciclo de periodo dos atractivo después del umbral

r=16/5.

  1. Puntos

    La fórmula da aproximadamente a=0,5130 y b=0,7995, distintos y con L(a)=b, L(b)=a.
  2. Retorno

    −r²+2r+4=−256/25+160/25+100/25=4/25.
  3. Estabilidad

    |4/25|<1: el ciclo atrae localmente.
Resultado

Ahora alternan dos valores, no un único punto fijo atractivo.

Lo comprobamos

Este intervalo de periodo dos no demuestra caos en todos los parámetros posteriores.

8. Una definición de caos necesita un conjunto y unas propiedades precisas

No hay una sola definición universal. Una formulación topológica habitual combina transitividad, puntos periódicos densos y sensibilidad en un conjunto invariante. Otra mide complejidad mediante entropía positiva. Un exponente de Lyapunov positivo aislado no prueba todas estas propiedades ni recurrencia.

Una ecuación escalar autónoma con unicidad no admite ciclos no constantes: entre equilibrios, el campo tiene signo constante y las soluciones son monótonas. En el plano, Poincaré–Bendixson restringe un conjunto omega-límite compacto sin equilibrios a una órbita periódica bajo sus hipótesis. Mapas, forzamiento temporal y dimensiones mayores requieren otro análisis.

Una recta de flujo no puede producir un ciclo no constante

x′=μ−x² entre dos equilibrios, con μ>0.

  1. Signo

    Entre −√μ y √μ, μ−x² es positivo.

  2. Movimiento

    Toda solución allí crece mientras permanezca en ese intervalo.

  3. Retorno

    Una función estrictamente creciente no puede volver a un valor anterior y cerrar una órbita no constante.

Resultado

Una silla-nodo escalar no es una demostración de caos.

Lo comprobamos

La conclusión usa autonomía y unicidad; no se traslada sin más a mapas discretos.

Oscilar periódicamente no basta para declarar caos

x′=y, y′=−x.

  1. Cantidad

    La derivada de x²+y² es 2xy−2xy=0.
  2. Órbitas

    Un dato no nulo recorre una circunferencia periódica.

  3. Interpretar

    Tener oscilación y muchas órbitas no verifica una definición de caos.

Resultado

La palabra caos requiere propiedades adicionales declaradas.

Lo comprobamos

Este ejemplo también muestra por qué no se decide el comportamiento solo por un dibujo que gira.

9. Una herradura necesita expansión, contracción y retorno

La herradura de Smale estira una dirección, contrae otra y devuelve dos franjas a un rectángulo. Los puntos que permanecen dentro para todos los tiempos forman un conjunto invariante tipo Cantor. En la herradura completa de dos ramas, la dinámica es conjugada al desplazamiento bilateral de dos símbolos.

Una intersección homoclínica transversal de una silla periódica de un difeomorfismo produce una herradura para una iteración adecuada, por Smale–Birkhoff. El retorno global y la transversalidad importan. La existencia de una silla local sin retorno no contiene por sí sola esa conclusión.

La codificación cuenta muchas historias

Herradura completa de dos ramas.

  1. Tres tiempos

    Hay 2³=8 palabras binarias de longitud tres.
  2. Periodos

    La palabra 01 repetida da una órbita de periodo dos del desplazamiento.

  3. Crecimiento

    Las 2ⁿ palabras dan tasa log2 de complejidad por iteración de ese modelo.

Resultado

La codificación describe comportamiento recurrente en un conjunto invariante.

Lo comprobamos

Si la herradura aparece para una potencia del mapa, se ajusta la unidad de tiempo al comparar entropías.

Una silla sin retorno tiene muy poca dinámica acotada en ambos sentidos

A(x,y)=(2x,y/2) y un rectángulo acotado alrededor de cero.

  1. Hacia delante

    Permanecer siempre dentro exige x=0.
  2. Hacia atrás

    A⁻ⁿ expande la coordenada y; permanecer dentro exige y=0.
  3. Intersección

    El único punto que queda dentro para todos los tiempos enteros es el origen.

Resultado

Contracción y expansión sin un mecanismo de retorno no forman esta herradura.

Lo comprobamos

El modelo lineal local no sustituye las hipótesis globales del teorema homoclínico.

10. Una conjugación exacta puede justificar un caso complejo que la simulación solo sugiere

El mapa tienda Q(y)=1−|2y−1| tiene dos ramas que cubren [0,1]. h(y)=sin²(πy/2) es un homeomorfismo creciente del intervalo y satisface L₄∘h=h∘Q. La identidad 4sin²a cos²a=sin²(2a), junto a la simetría de las dos ramas, verifica esa conjugación.

L₄ hereda las propiedades topológicas correspondientes del mapa tienda. La codificación binaria requiere atender a puntos de frontera; no afirmamos que todo el intervalo sea homeomorfo a un espacio de Cantor de sucesiones. Una gráfica numérica de bifurcación puede mostrar transitorios largos y ventanas periódicas; no es una prueba de todos los parámetros.

Además, una implementación determinista en un conjunto finito de estados termina repitiendo alguno y después se hace periódica. Eso describe la precisión finita del cálculo y no demuestra que la órbita real correspondiente tenga ese periodo.

Comprobar la conjugación con un punto sencillo

y=1/3.

  1. Cambiar coordenada

    h(1/3)=sin²(π/6)=1/4.
  2. Logística

    L₄(1/4)=4·(1/4)·(3/4)=3/4.
  3. Tienda

    Q(1/3)=2/3 y h(2/3)=sin²(π/3)=3/4.
Resultado

Los dos caminos dan la misma imagen.

Lo comprobamos

La identidad trigonométrica y la reflexión de la segunda rama prueban la igualdad para todos los puntos, no solo este.

Un periodo calculado puede ser un efecto de precisión finita

Una implementación usa N estados posibles en una región invariante.

  1. Contar

    En N+1 visitas algún estado computado se repite.

  2. Determinismo

    Desde el mismo estado la implementación sigue la misma sucesión: aparece un ciclo.

  3. Alcance

    No se ha probado un ciclo idéntico de la órbita definida con números reales exactos.

Resultado

El cálculo ayuda a explorar; la conjugación y los teoremas aportan otra clase de evidencia.

Lo comprobamos

Descartar transitorios y revisar precisión mejora la exploración, pero no elimina esta diferencia lógica.

Ahora lo resolvemos de principio a fin.

EJEMPLO RESUELTO 1

Parámetro positivo

Para μ=4, halla y clasifica equilibrios.

  1. Modelo

    Usamos x′=μ−x² y μ=4.
  2. Equilibrios

    Se impone 0=4−x², por tanto x²=4.
  3. Raíces

    Los dos valores son dos y menos dos.

  4. Derivada

    fₓ=−2x: en dos vale −4 y en menos dos vale cuatro.

  5. Clasificar

    La derivada negativa hace atractivo a dos; la positiva hace repulsor a menos dos.

  6. Comprobar

    Las flechas del campo apuntan hacia dos localmente y se alejan de menos dos.

Resultado

2 estable y −2 inestable.

Lo comprobamos

Las flechas del campo apuntan hacia 2 localmente.

EJEMPLO RESUELTO 2

Parámetro negativo

Para μ=−1, ¿hay equilibrios reales?

  1. Modelo

    x′=μ−x², ahora μ=−1.
  2. Equilibrio

    Pedir x′=0 da −1−x²=0.
  3. Reordenar

    Eso exigiría x²=−1.
  4. Reales

    El cuadrado de un real es no negativo; no hay solución real.

  5. Campo

    −1−x² es estrictamente negativo para todo x.

  6. Conclusión

    No hay equilibrios reales; la pareja de ramas solo existe en el lado μ>0.

Resultado

No hay equilibrios reales.

Lo comprobamos

El campo es estrictamente negativo en toda R.

Los símbolos que hemos usado

x′=f(x,μ)
Familia de ecuaciones diferenciales autónomas con estado x y parámetro real μ.
x*(μ)
Rama de equilibrios, definida por f(x*(μ),μ)=0.
fₓ y fμ
Derivadas parciales del campo respecto del estado y del parámetro.
Hiperbólicidad
Ausencia de autovalores con parte real cero para un equilibrio de un flujo; para un punto fijo de un mapa se excluye módulo uno.
Forma normal
Modelo local simplificado que conserva los términos decisivos para un tipo de bifurcación.
Multiplicador
Autovalor de la derivada de un mapa o del mapa de retorno de una órbita periódica.
Lᵣ(x)=rx(1−x)
Mapa logístico discreto; para0≤r≤4 preserva el intervalo[0,1].
Σ₂ y σ
Espacio de sucesiones de dos símbolos y aplicación que desplaza cada sucesión una posición.
Homoclínica transversal
Intersección transversal, fuera del punto periódico base, entre sus variedades estable e inestable.

Ahora te toca a ti.

Tienes cinco ejercicios. Inténtalos; las pistas te ayudan si te atascas.

Practicar este tema →

Para ampliar