8.3: Bifurcaciones bidimensionales
- Page ID
- 116980
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)una bifurcación Hopf. La bifurcación Hopf viene en dos tipos: bifurcación Hopf supercrítica y bifurcación Hopf subcrítica. Para la bifurcación supercrítica de Hopf, a medida que\(\mu\) aumenta ligeramente por encima de cero, la oscilación resultante alrededor del punto fijo ahora inestable se estabiliza rápidamente a pequeña amplitud. Esta órbita estable se denomina ciclo límite. Para la bifurcación subcrítica de Hopf, a medida que\(\mu\) aumenta ligeramente por encima de cero, el ciclo límite salta inmediatamente a gran amplitud.
Bifurcación supercrítica de Hopf
Un ejemplo simple de una bifurcación supercrítica de Hopf se puede dar en coordenadas polares:
\[\begin{aligned}\overset{.}{r}&=\mu r-r^3 , \\ \overset{.}{\theta}&=\omega +br^2,\end{aligned}\]dónde\(x = r \cos\theta\) y\(y = r \sin\theta\). El parámetro\(\mu\) controla la estabilidad del punto fijo en el origen, el parámetro\(\omega\) es la frecuencia de oscilación cerca del origen, y el parámetro\(b\) determina la dependencia de la frecuencia de oscilación a oscilaciones de mayor amplitud. Aunque incluimos\(b\) por generalidad, nuestro análisis cualitativo de estas ecuaciones será independiente de\(b\).
La ecuación para el radio es de la forma de bifurcación de horca supercrítica. Los puntos fijos son\(r_* = 0\) y\(r_* =\sqrt{\mu}\) (nótese que\(r > 0\)), y el primer punto fijo es estable para\(\mu < 0\) y el segundo es estable para\(\mu > 0\). La transición de los valores propios del jacobiano de la parte real negativa a la parte real positiva se puede ver si transformamos estas ecuaciones en coordenadas cartesianas. Tenemos usando\(r^2 = x^2 + y^2\),\[\begin{aligned}\overset{.}{x}&=\overset{.}{r}\cos\theta -\overset{.}{\theta}r\sin\theta \\ &=(\mu r-r^3)\cos\theta -(\omega +br^2)r\sin\theta \\ &=\mu x-(x^2+y^2)x-\omega y-b(x^2+y^2)y \\ &=\mu x-\omega y-(x^2+y^2)(x+by); \\ \overset{.}{y}&=\overset{.}{r}\sin\theta +\overset{.}{\theta}r\cos\theta \\ &=(\mu r-r^3)\sin\theta +(\omega +br^2)r\cos\theta \\ &=\mu y-(x^2+y^2)y+\omega x+b(x^2+y^2)x \\ &=\omega x+\mu y-(x^2+y^2)(y-bx).\end{aligned}\]
La estabilidad del origen está determinada por la matriz jacobiana evaluada en el origen. Los términos no lineales en la ecuación desaparecen y la matriz jacobiana en el origen viene dada por\[J=\left(\begin{array}{cc}\mu &-\omega \\ \omega &\mu\end{array}\right).\nonumber\]
Los valores propios son las soluciones de\((\mu −\lambda)^2 + \omega^2 = 0\), o\(\lambda = \mu ± i\omega\). A medida que\(\mu\) aumenta de valores negativos a positivos, las oscilaciones amortiguadas exponencialmente cambian en oscilaciones de crecimiento exponencial. Los términos no lineales en las ecuaciones estabilizan las oscilaciones crecientes en un ciclo límite.
Bifurcación subcrítica de Hopf
El ejemplo análogo de una bifurcación subcrítica Hopf viene dado por\[\begin{aligned}\overset{.}{r}&=\mu r+r^3-r^5, \\ \overset{.}{\theta}&=\omega +br^2.\end{aligned}\]
Aquí, la ecuación para el radio es de la forma de la bifurcación de horca subcrítica. A medida que\(\mu\) aumenta de valores negativos a positivos, el origen se vuelve inestable y las oscilaciones de crecimiento exponencial aumentan hasta que el radio alcanza un punto fijo estable lejos del origen. En la práctica, puede ser difícil determinar analíticamente si una bifurcación Hopf es supercrítica o subcrítica a partir de las ecuaciones de movimiento. La solución computacional, sin embargo, puede distinguir rápidamente entre los dos tipos.