8.2: Bifurcaciones unidimensionales

Una bifurcación ocurre en una ecuación diferencial no lineal cuando un pequeño cambio en un parámetro da como resultado un cambio cualitativo en la solución a largo plazo. Ejemplos de bifurcaciones son cuando se crean o destruyen puntos fijos, o cambian su estabilidad.

Ahora consideramos cuatro bifurcaciones clásicas de ecuaciones diferenciales no lineales unidimensionales: bifurcación de nodo de silla de montar, bifurcación transcrítica, bifurcación de horca supercrítica y bifurcación de horca subcrítica. La ecuación diferencial correspondiente se escribirá como\[\overset{.}{x}=f_r(x),\nonumber\] donde el subíndice\(r\) representa un parámetro que da como resultado una bifurcación cuando se varía a través de cero. Las ecuaciones diferenciales más simples que exhiben estas bifurcaciones se denominan las formas normales, y corresponden a un análisis local (es decir, expansión de la serie Taylor) de ecuaciones diferenciales más generales alrededor del punto fijo, junto con una posible reescalación de\(x\).

Bifurcación de nodo de silla de montar

Ver tutorial en YouTube

Figura\(\PageIndex{1}\): Bifurcación de Saddlenode. (a)\(\overset{.}{x}\) versus\(x\); (b) diagrama de bifurcación.

La bifurcación del nodo de silla da como resultado la creación o destrucción de puntos fijos. La forma normal para una bifurcación de nodo de silla de montar viene dada por\[\overset{.}{x}=r+x^2.\nonumber\]

Los puntos fijos son\(x_* = ±\sqrt{-r}\). Claramente, existen dos puntos fijos reales cuando\(r < 0\) y no existen puntos fijos reales cuando\(r > 0\). La estabilidad de los puntos fijos cuando\(r < 0\) están determinados por la derivada de\(f(x) = r + x^2\), dada por\(f' (x) = 2x\). Por lo tanto, el punto fijo negativo es estable y el punto fijo positivo es inestable.

Gráficamente, podemos ilustrar esta bifurcación de dos maneras. Primero, en la Fig. \(\PageIndex{1}\)(a), trazamos\(\overset{.}{x}\) versus\(x\) para los tres valores de parámetros correspondientes a\(r < 0\),\(r = 0\) y\(r > 0\). Los valores en los que\(\overset{.}{x} = 0\) corresponden a los puntos fijos, y se dibujan flechas que indican cómo\(x(t)\) evoluciona la solución (a la derecha si\(\overset{.}{x}> 0\) y a la izquierda si\(\overset{.}{x}< 0\)). El punto fijo estable se indica mediante un círculo relleno y el punto fijo inestable por un círculo abierto. Tenga en cuenta que cuando\(r = 0\), las soluciones convergen al origen desde la izquierda, pero divergen del origen de la derecha. Segundo, en la Fig. \(\PageIndex{1}\)(b), trazamos un diagrama de bifurcación que ilustra el punto fijo\(x_*\) versus el parámetro de bifurcación\(r\). El punto fijo estable se denota con una línea continua y el punto fijo inestable por una línea discontinua. Tenga en cuenta que los dos puntos fijos chocan y aniquilan en\(r = 0\), y no hay puntos fijos para\(r > 0\).

Bifurcación transcrítica

Ver tutorial en YouTube

Figura\(\PageIndex{2}\): Bifurcación transcrítica. (a)\(\overset{.}{x}\) versus\(x\); (b) diagrama de bifurcación.

Una bifurcación transcrítica ocurre cuando hay un intercambio de estabilidades entre dos puntos fijos. La forma normal para una bifurcación transcrítica viene dada por\[\overset{.}{x}=rx-x^2.\nonumber\]

Los puntos fijos son\(x_* = 0\) y\(x_* = r\). El derivado del lado derecho es\(f' (x) = r − 2x\), para que\(f' (0) = r\) y\(f ' (r) = −r\). Por lo tanto\(r < 0\), for,\(x_* = 0\) es estable y\(x_* = r\) es inestable, mientras que for\(r > 0\),\(x_* = r\) es estable y\(x_* = 0\) es inestable. Los dos puntos fijos intercambian así la estabilidad a medida que\(r\) pasa por cero. La bifurcación transcrítica se ilustra en la Fig. \(\PageIndex{2}\).

Bifurcación supercrítica

Ver tutorial en YouTube

Figura\(\PageIndex{3}\): Bifurcación supercrítica de tridente. (a)\(\overset{.}{x}\) versus\(x\); (b) diagrama de bifurcación.

Las bifurcaciones de horca ocurren en modelos físicos donde aparecen puntos fijos y desaparecen en pares debido a alguna simetría intrínseca del problema. Las bifurcaciones de horca pueden venir en uno de dos tipos. En la bifurcación supercrítica, se crean un par de puntos fijos estables en el punto de bifurcación (o crítico) y existen después de (súper) la bifurcación. En la bifurcación subcrítica, se crea un par de puntos fijos inestables en el punto de bifurcación y existen antes de (sub) la bifurcación.

La forma normal para la bifurcación supercrítica de horca viene dada por\[\overset{.}{x}=rx-x^3.\nonumber\]

Tenga en cuenta que el término lineal da como resultado un crecimiento exponencial cuando\(r > 0\) y el término no lineal estabiliza este crecimiento. Los puntos fijos son\(x_* = 0\) y\(x_* = ±\sqrt{r}\), estos últimos puntos fijos existentes sólo cuando\(r > 0\). El derivado de\(f\) es\(f' (x) = r − 3x^2\) para que\(f' (0) = r\) y\(f' (±\sqrt{r}) = −2r\). Por lo tanto, el punto fijo\(x_* = 0\) es estable para\(r < 0\) e inestable\(r > 0\) mientras los puntos fijos\(x = ±\sqrt{r}\) existen y son estables para\(r > 0\). Observe que el punto fijo\(x_* = 0\) se vuelve inestable a medida que\(r\) cruza cero y\(x_* = ±\sqrt{r}\) nacen dos nuevos puntos fijos estables. La bifurcación supercrítica de horca se ilustra en la Fig. \(\PageIndex{3}\).

Bifurcación supercrítica

Ver tutorial en YouTube

Figura\(\PageIndex{4}\): Bifurcación subcrítica de horca.

En el caso subcrítico, el término cúbico es desestabilizador. La forma normal (a la orden\(x^3\)) es\[\overset{.}{x}=rx+x^3.\nonumber\]

Los puntos fijos son\(x_* = 0\) y\(x_* = ±\sqrt{-r}\), estos últimos puntos fijos existentes sólo cuando\(r ≤ 0\). El derivado del lado derecho es\(f' (x) = r +3x^2\) así que\(f' (0) = r\) y\(f' (±\sqrt{-r}) = −2r\). Por lo tanto, el punto fijo\(x_* = 0\) es estable para\(r < 0\) e inestable\(r > 0\) mientras los puntos fijos\(x = ±\sqrt{-r}\) existen y son inestables para\(r < 0\). No hay puntos fijos estables cuando\(r > 0\).

La ausencia de puntos fijos estables para\(r > 0\) indica que el descuido de términos de orden\(x\) superior en que\(x^3\) en la forma normal puede ser injustificado. Manteniendo la simetría intrínseca de las ecuaciones (solo potencias impares de\(x\)) podemos agregar un término estabilizador no lineal proporcional a\(x^5\). La forma normal extendida (por orden\(x^5\)) es\[\overset{.}{x}=rx+x^3-x^5,\nonumber\] y es algo más difícil de analizar. Los puntos fijos son soluciones de\[x(r+x^2-x^4)=0.\nonumber\]

El punto fijo\(x_* = 0\) existe para todos\(r\), y se pueden encontrar cuatro puntos fijos adicionales a partir de las soluciones de la ecuación cuadrática en\(x^2\):

\[x_*=\pm\sqrt{\frac{1}{2}(1\pm\sqrt{1+4r})}.\nonumber\]

Estos puntos fijos existen sólo si\(x_*\) es real. Claramente, para que la raíz cuadrada interior sea real,\(r ≥ −1/4\). También observa que\(1 −\sqrt{1+4r}\) se vuelve negativo para\(r > 0\). Por lo tanto, tenemos tres intervalos\(r\) a considerar, y estas regiones y sus puntos fijos son\[\begin{array}{rlll} r<-\frac{1}{4}: & x_*=0 &\text{(one fixed point);} \\ -\frac{1}{4}<r<0: & x_*=0, &x_*=\pm\sqrt{\frac{1}{2}(1\pm\sqrt{1+4r})} &\text{(five fixed points);} \\ r>0: &x_*=0, &x_*=\pm\sqrt{\frac{1}{2}(1+\sqrt{1+4r})} &\text{(three fixed points).}\end{array}\nonumber\]

La estabilidad se determina a partir de\(f' (x) = r + 3x^2 − 5x^4\). Ahora,\(f' (0) = r\) así\(x_* = 0\) es estable para\(r < 0\) e inestable para\(r > 0\). El cálculo para las otras cuatro raíces puede simplificarse señalando que\(x_*\) satisface\(r + x_*^2 − x_*^4 = 0\), o\(x_*^4 = r + x_*^2\). Por lo tanto,\[\begin{aligned}f'(x_*)&=r+3x_*^2-5x_*^4 \\ &=r+3x_*^2-5(r+x_*^2) \\ &=-4r-2x_*^2 \\ &=-2(2r+x_*^2).\end{aligned}\]

Con\(x_*^2=\frac{1}{2}(1\pm\sqrt{1+4r})\), tenemos\[\begin{aligned} f'(x_*)&=-2\left(2r+\frac{1}{2}(1\pm\sqrt{1+4r})\right) \\ &=-\left( (1+4r)\pm\sqrt{1+4r}\right) \\ &=-\sqrt{1+4r}\left(\sqrt{1+4r}\pm 1\right).\end{aligned}\]

Claramente, la raíz plus siempre es estable desde entonces\(f' (x_*) < 0\). La raíz menos existe solo para\(−\frac{1}{4}< r < 0\) y es inestable desde entonces\(f' (x_*) > 0\). Resumimos la estabilidad de los diversos puntos fijos:

\[\begin{array}{rll} r<-\frac{1}{4}:&x_*=0&\text{(stable);} \\ -\frac{1}{4}<r<0:&x_*=0,&\text{(stable)} \\ &x_*=\pm\sqrt{\frac{1}{2}(1+\sqrt{1+4r})}&\text{(stable);} \\ &x_*=\pm\sqrt{\frac{1}{2}(1-\sqrt{1+4r})}&\text{(unstable);} \\ r>0:&x_*=0&\text{(unstable)} \\ &x_*=\pm\sqrt{\frac{1}{2}(1+\sqrt{1+4r})}&\text{(stable).}\end{array}\nonumber\]

El diagrama de bifurcación se muestra en la Fig. \(\PageIndex{4}\), y consiste en una bifurcación subcrítica de horca en\(r = 0\) y dos bifurcaciones de nodo de silla de montar en\(r = −1/4\). Podemos imaginar lo que sucede a\(x\) medida que\(r\) aumenta a partir de valores negativos, suponiendo que haya algún pequeño ruido en el sistema para\(x = x(t)\) que diverja de puntos fijos inestables. Para\(r < −1/4\), el valor de equilibrio de\(x\) es\(x_* = 0\). A\(r\) medida que aumenta en el rango\(−1/4 < r < 0\),\(x\) se mantendrá en\(x_* = 0\). Sin embargo, una catástrofe ocurre tan pronto como\(r > 0\). El punto\(x_* = 0\) fijo se vuelve inestable y la solución saltará hacia arriba (o hacia abajo) al único punto fijo estable restante. Tal comportamiento se llama bifurcación de salto. Una catástrofe similar puede ocurrir como\(r\) disminuciones a partir de valores positivos. En este caso, el salto ocurre tan pronto como sea\(r < −1/4\).

Dado que el valor de equilibrio estable de\(x\) depende de si estamos aumentando o disminuyendo\(r\), decimos que el sistema exhibe histéresis. La existencia de una bifurcación de horca subcrítica puede ser muy peligrosa en aplicaciones de ingeniería, ya que un pequeño cambio en los parámetros de un problema puede resultar en un gran cambio en el estado de equilibrio. Físicamente, esto puede corresponder a un colapso de una estructura, o al fallo de un componente.