Diagramas de Fase

Sentido de Giro de las Espirales para el Analisis
Cualitativo de SNL. 2da Version

Francisco J. Bravo S.
1 de septiembre de 2011
1. Introduccion
El poder graficar correctamente el sentido de giro de los puntos crticos
en los SNL cuando estos resultan ser espirales debe ser una de las dudas mas
comunes que surgen entre los alumnos. Una de las soluciones mas faciles es
hacer calzar el sentido de giro con el de los otros puntos crticos conocidos
en el diagrama de fase (esto dado que las soluciones del sistema lineal son
continuas). En esta gua me propongo explicar el fundamento y metodologa
necesaria para poder determinar el sentido y forma de los puntos crticos
espirales.
Recordemos que estos puntos crticos se hacen para sistemas 2 2, y que
en dichos puntos la matriz del sistema linealizado posee 2 valores propios
complejos conjugados a ib.
La primera parte es la justificacion del procedimiento hecho en la seccion
Algoritmo. Si no esta interesado en la teora, pude perfectamente pasar a
esta ultima seccion
2. Desarrollo teorico
Teorema 1 Sea A una matriz real de 2 2 con valor propio complejo =
a + ib (b 6= 0) y vector propio asociado v = u1 + iu2 C2 . Entonces

1 a b
A = P CP ; con P = [u1 ; u2 ] C =
b a
Demostracion: Si = a + ib, a, b R es un autovalor de A (matriz cuadrada
real) y v = u1 + iu2 (u1 , u2 Rm ) es un autovector de A asociado a ,
entonces v = u1 iu2 es autovector de A asociado a = a ib y, por tanto,
tenemos las igualdades
Av = v = (a + bi)(u1 + iu2 )) Au1 + iAu2 = (au1 bu2 ) + i(bu1 + au2 )

Av = v = (a bi)(u1 iu2 )) Au1 iAu2 = (au1 bu2 ) i(bu1 + au2 )
1
y as, identificando las partes real e imaginaria en cualquiera de las dos
igualdades anteriores tenemos

Au1 = au1 bu2
Au2 = bu1 + au2
Expresando estas igualdades de forma matricial obtenemos

a b
A[u1 , u2 ] = [u1 , u2 ]
b a
Este teorema nos dice que la matriz A del sistema es semejante a la

matriz C
Para el sistema X 0 = AX, con A cumpliendo las condiciones del teorema,
podemos hacer
X 0 = AX
X 0 = P CP 1 X
1 0 1
| {zX} = C P
P | {zX}
Y0 Y
Y 0 = CY
Luego, basta resolver el sistema Y 0 = CY , y volver a la variable original

X = P Y . De este modo, la trayectoria de la espiral que obtengamos al re-
solver el sistema para Y , sera la trayectoria que dibujaremos, pero sobre el
eje formado por los vectores < u1 , u2 >.
Sabemos que la solucion del sistema

0 a b
Y = Y
b a
es
y01

at cos(bt) sin(bt)
Y =e
sin(bt) cos(bt) y02
2
Para este sistema sabemos cuales son sus trayectorias para los valores de
ayb
Figura 1: A la izquierda a < 0, b > 0. A la derecha a > 0, b > 0
3. Algoritmo
Luego, el algoritmo para dibujar diagramas para un sistema X 0 = AX
es el siguiente:
1. Calcular los valores propios = a ib
2. Escoger el valor propio con b > 0
3. Calcular el vector propio v = u1 + iu2 asociado al valor propio elegido
4. Dado el valor de a, el diagrama es uno de los dos de la Fig (1)
5. Dibujar el diagrama final, graficando los vectores u1 y u2 en el plano,

con el dibujo de la espiral
Ejemplo:
0 1
Dibujemos el diagrama del sistema X = X
1 0
Los valores propios del sistema
son =i. Elegimos
el valor propio
1 1 0
= i, cuyo vector propio es = +i
i 0 1
| {z } | {z }
u1 u2
3
Luego, como elegimos el valor propio con parte imaginaria positiva, las
trayectorias giran en sentido horario, en el diagrama con ejes u1 u2 ( Fig (2))
Figura 2: Trayectoria del sistema en el espacio de base < u1 , u2 >
Ahora trasladamos este diagrama.

El vector u1 descansa sobre el eje x.

0
El vector u2 tiene coordenadas , luego al graficarlo en el diagrama
1
final, la figura se invierte, obteniendo las trayectorias de nuestro sistema
inicial
Figura 3: Trayectorias del sistema X=AX
Notar que el sentido de giro cambio a antihorario, debido a la rotacion

del eje que tuvimos que realizar.
4
Con este algoritmo, note que solo es necesario conocer 2 diagramas para
las espirales, los del sistema Y 0 = CY , con b > 0
Referencias
[1] Departamento de Matematica Aplicada II. Universidad de Sevi-
lla. Tema 11.- Autovalores y Autovectores http://personal.us.es/
ealgaba/algebra%20=04-05/TEMA11.pdf
[2] A. Osses. Apunte del curso: Ecuaciones diferenciales ordinarias

Diagramas de Fase

Cargado por

Copyright:

Formatos disponibles

Diagramas de Fase

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Diagramas de Fase

Cargado por

Copyright:

Formatos disponibles

Sentido de Giro de las Espirales para el Analisis

Cualitativo de SNL. 2da Version

Av = v = (a + bi)(u1 + iu2 )) Au1 + iAu2 = (au1 bu2 ) + i(bu1 + au2 )

Este teorema nos dice que la matriz A del sistema es semejante a la

Luego, basta resolver el sistema Y 0 = CY , y volver a la variable original

Sabemos que la solucion del sistema

Figura 1: A la izquierda a < 0, b > 0. A la derecha a > 0, b > 0

1. Calcular los valores propios = a ib

2. Escoger el valor propio con b > 0

3. Calcular el vector propio v = u1 + iu2 asociado al valor propio elegido

4. Dado el valor de a, el diagrama es uno de los dos de la Fig (1)

5. Dibujar el diagrama final, graficando los vectores u1 y u2 en el plano,

Figura 2: Trayectoria del sistema en el espacio de base < u1 , u2 >

Ahora trasladamos este diagrama.

Figura 3: Trayectorias del sistema X=AX

Notar que el sentido de giro cambio a antihorario, debido a la rotacion

[2] A. Osses. Apunte del curso: Ecuaciones diferenciales ordinarias

También podría gustarte