11.1: Puntos Fijos y Estabilidad
- Page ID
- 119189
\( \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}}} \)
11.1.1. Una dimensión
Considere la ecuación diferencial unidimensional para\(x=x(t)\) dada por
\[\dot{x}=f(x) . \nonumber \]
Decimos que\(x_{*}\) es un punto fijo, o punto de equilibrio, de la Ecuación\ ref {11.1} si\(f\left(x_{*}\right)=0\). En el punto fijo,\(\dot{x}=0\). Se utiliza la terminología punto fijo ya que la solución a la Ecuación\ ref {11.1} con condición inicial\(x(0)=x_{*}\) es\(x(t)=x_{*}\) para todos los tiempos\(t\).
Un punto fijo, sin embargo, puede ser estable o inestable. Se dice que un punto fijo es estable si una pequeña perturbación de la solución desde el punto fijo decae en el tiempo; se dice que es inestable si una pequeña perturbación crece en el tiempo. Podemos determinar la estabilidad mediante un análisis lineal. Let\(x=x_{*}+\epsilon(t)\), donde\(\epsilon\) representa una pequeña perturbación de la solución desde el punto fijo\(x_{*}\). Porque\(x_{*}\) es una constante,\(\dot{x}=\dot{\epsilon}\); y porque\(x_{*}\) es un punto fijo,\(f\left(x_{*}\right)=0\). Serie Taylor expandiéndose sobre\(\epsilon=0\), tenemos
\[\begin{aligned} \dot{\epsilon} &=f\left(x_{*}+\epsilon\right) \\ &=f\left(x_{*}\right)+\epsilon f^{\prime}\left(x_{*}\right)+\ldots \\ &=\epsilon f^{\prime}\left(x_{*}\right)+\ldots \end{aligned} \nonumber \]
Los términos omitidos en la expansión de la serie Taylor son proporcionales a\(\epsilon^{2}\), y pueden hacerse despreciables en un corto intervalo de tiempo con respecto al término guardado, proporcionales a\(\epsilon\), tomando\(\epsilon(0)\) suficientemente pequeños. Por lo tanto, al menos en tiempos cortos, la ecuación diferencial a considerar\(\dot{\epsilon}=f^{\prime}\left(x_{*}\right) \epsilon\), es lineal y tiene por ahora la solución familiar
\[\epsilon(t)=\epsilon(0) e^{f^{\prime}\left(x_{*}\right) t} \nonumber \]
La perturbación de la solución de punto fijo decae\(x(t)=x_{*}\) así exponencialmente si\(f^{\prime}\left(x_{*}\right)<0\), y decimos que el punto fijo es estable. Si\(f^{\prime}\left(x_{*}\right)>0\), la perturbación crece exponencialmente y decimos que el punto fijo es inestable. Si\(f^{\prime}\left(x_{*}\right)=0\), decimos que el punto fijo es marginalmente estable, y se debe considerar el siguiente término de orden superior en la expansión de la serie Taylor. Ejemplo: Encontrar todos los puntos fijos de la ecuación logística\(\dot{x}=x(1-x)\) y determinar su estabilidad.
Hay dos puntos fijos en los que\(\dot{x}=0\), dado por\(x_{*}=0\) y\(x_{*}=1\). La estabilidad de estos puntos de equilibrio se puede determinar considerando la derivada de\(f(x)=x(1-x)\). Tenemos\(f^{\prime}(x)=1-2 x\). Por lo\(f^{\prime}(0)=1>0\) tanto, entonces ese\(x_{*}=0\) es un punto fijo inestable, y\(f^{\prime}(1)=-1<0\) así ese\(x_{*}=1\) es un punto fijo estable. En efecto, previamente hemos encontrado que todas las soluciones se acercan asintóticamente al punto fijo estable.
11.1.2. Dos dimensiones
La idea de puntos fijos y estabilidad puede extenderse a sistemas de odas de orden superior. Aquí, consideramos un sistema bidimensional y necesitaremos hacer uso de la expansión bidimensional de la serie Taylor de una función\(F(x, y)\) sobre el origen. En general, la serie Taylor de\(F(x, y)\) está dada por
\[F(x, y)=F+x \frac{\partial F}{\partial x}+y \frac{\partial F}{\partial y}+\frac{1}{2}\left(x^{2} \frac{\partial^{2} F}{\partial x^{2}}+2 x y \frac{\partial^{2} F}{\partial x \partial y}+y^{2} \frac{\partial^{2} F}{\partial y^{2}}\right)+\ldots, \nonumber \]
donde la función\(F\) y todas sus derivadas parciales en el lado derecho se evalúan en el origen. Tenga en cuenta que la serie Taylor está construida de manera que todas las derivadas parciales del lado izquierdo coincidan con las del lado derecho en el origen.
La matriz de dos por dos en la Ecuación\ ref {11.3} se llama la matriz jacobiana en el punto fijo. Un análisis de valores propios de la matriz jacobiana normalmente dará dos valores propios\(\lambda_{1}\) y\(\lambda_{2}\). Estos valores propios pueden ser reales y distintos, pares conjugados complejos o repetidos. El punto fijo es estable (todas las perturbaciones decaen exponencialmente) si ambos valores propios tienen partes reales negativas. El punto fijo es inestable (algunas perturbaciones crecen exponencialmente) si al menos uno de los valores propios tiene una parte real positiva. Los puntos fijos se pueden clasificar además como nodos estables o inestables, puntos de sillín inestables, puntos espirales estables o inestables, o nodos inestables o inestables.
Ejemplo: Encontrar todos los puntos fijos del sistema no lineal\(\dot{x}=x(3-x-2 y)\),\(\dot{y}=y(2-x-y)\), y determinar su estabilidad.
Los puntos fijos se determinan resolviendo
\[f(x, y)=x(3-x-2 y)=0, \quad g(x, y)=y(2-x-y)=0 \nonumber \]
Evidentemente,\((x, y)=(0,0)\) es un punto fijo. Por un lado, aunque sólo sea\(x=0\), entonces la ecuación\(g(x, y)=0\) rinde\(y=2\). Por otro lado, aunque sólo sea\(y=0\), entonces la ecuación\(f(x, y)=0\) rinde\(x=3\). Si ambos\(x\) y\(y\) son distintos de cero, entonces debemos resolver el sistema lineal
\[x+2 y=3, \quad x+y=2 \text {, } \nonumber \]
y la solución se encuentra fácilmente para ser\((x, y)=(1,1)\). De ahí que hayamos determinado los cuatro puntos fijos\(\left(x_{*}, y_{*}\right)=(0,0),(0,2),(3,0),(1,1)\). La matriz jacobiana viene dada por
\[\left(\begin{array}{ll} \frac{\partial f}{\partial x} & \frac{\partial f}{\partial y} \\ \frac{\partial g}{\partial x} & \frac{\partial g}{\partial y} \end{array}\right)=\left(\begin{array}{cc} 3-2 x-2 y & -2 x \\ -y & 2-x-2 y \end{array}\right) . \nonumber \]
La estabilidad de los puntos fijos puede considerarse a su vez. Con\(J_{*}\) la matriz jacobiana evaluada en el punto fijo, tenemos
\[\left(x_{*}, y_{*}\right)=(0,0): \quad J_{*}=\left(\begin{array}{ll} 3 & 0 \\ 0 & 2 \end{array}\right) \nonumber \]
Los valores propios de\(J_{*}\) son\(\lambda=3,2\) para que el punto fijo\((0,0)\) sea un nodo inestable. Siguiente,
\[\left(x_{*}, y_{*}\right)=(0,2): \quad \mathrm{J}_{*}=\left(\begin{array}{rr} -1 & 0 \\ -2 & -2 \end{array}\right) . \nonumber \]
Los valores propios de\(\mathrm{J}_{*}\) son\(\lambda=-1,-2\) para que el punto fijo\((0,2)\) sea un nodo estable. Siguiente,
\[\left(x_{*}, y_{*}\right)=(3,0): \quad J_{*}=\left(\begin{array}{rr} -3 & -6 \\ 0 & -1 \end{array}\right) \text {. } \nonumber \]
Los valores propios de\(J_{*}\) son\(\lambda=-3,-1\) para que el punto fijo\((3,0)\) sea también un nodo estable. Por último,
\[\left(x_{*}, y_{*}\right)=(1,1): \quad J_{*}=\left(\begin{array}{ll} -1 & -2 \\ -1 & -1 \end{array}\right) \text {. } \nonumber \]
La ecuación característica de\(\mathrm{J}_{*}\) viene dada por\((-1-\lambda)^{2}-2=0\), así que eso\(\lambda=-1 \pm \sqrt{2}\). Dado que un valor propio es negativo y el otro positivo, el punto fijo\((1,1)\) es un punto de sillín inestable. A partir de nuestro análisis de los puntos fijos, se puede esperar que todas las soluciones sean asíntotas a uno de los puntos fijos estables\((0,2)\) o\((3,0)\), dependiendo de las condiciones iniciales.
Es de interés esbozar el diagrama de espacio de fase para este sistema no lineal. Los vectores propios asociados con el punto de sillín inestable\((1,1)\) determinan las direcciones del flujo hacia y desde este punto fijo. El vector propio asociado con el valor propio positivo se\(\lambda_{1}=-1+\sqrt{2}\) puede determinar a partir de la primera ecuación de\(\left(\mathrm{J}_{*}-\lambda_{1} \mathrm{I}\right) \mathrm{V}_{1}=0\), o
\[-\sqrt{2} v_{11}-2 v_{12}=0 \nonumber \]
para que\(v_{12}=-(\sqrt{2} / 2) v_{11}\). El vector propio asociado con el valor propio negativo\(\lambda_{1}=-1-\sqrt{2}\) satisface\(v_{22}=(\sqrt{2} / 2) v_{21}\). Los vectores propios dan la pendiente de las líneas con origen en el punto fijo para trayectorias entrantes (valor propio negativo) y salientes (valor propio positivo). Las trayectorias salientes tienen pendiente negativa\(-\sqrt{2} / 2\) y las trayectorias entrantes tienen pendiente positiva\(\sqrt{2} / 2\). Un boceto aproximado del diagrama de espacio de fase se puede hacer a mano (como se demuestra en clase). Aquí, una gráfica generada por computadora obtenida a partir de la solución numérica de las odas acopladas no lineales se presenta en la Fig. 11.1. La curva que comienza desde el origen y en el infinito, y que termina en el punto inestable de la silla de montar se llama separatriz. Esta curva separa el espacio de fase en dos regiones: condiciones iniciales para las cuales la solución asíntota al punto fijo\((0,2)\), y condiciones iniciales para las cuales la solución asíntota al punto fijo\((3,0)\).