Apunte Edo Cualitativas

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 27

1.

Comportamiento cualitativo de las EDO

La mayora de las ED no puede resolverse explcitamente. Interesa sin embargo describir cmo se comportan las soluciones, por ejemplo
1. al aumentar t crecen sin cota las soluciones?
2. tienden las soluciones a algn valor (0 por ejemplo)
3. oscilan entre ciertos valores?
Definicin 1. Sea D R R, F : D R. Una expresin de la forma
dx
= F (t, x)
dt
o x = F (t, x) es una EDO de primer orden.
Observacin 1. Si F (t, x) = F (x), i.e., F depende slo de la variable x, entonces
x = F (x),
y en este caso la ecuacin se dice autonma.
Teorema 1 (existencia de soluciones). Sea F : R una funcin continua
en el dominio . Para todo punto (x0 , y0 ) existe una funcin : I R
definida en algn intervalo I, solucin del problema de valor inicial llamado de
Cauchy
dx
= F (t, x), tal que (t0 ) = x0
dt
La figura 1 ilustra el teorema con dos soluciones 1 , 2
x

x = 2 (t)

x = 1 (t)
t

Figura 1

F
son funciones continuas en el dox
minio . Para todo punto (t0 , x0 ) existe una nica funcin : I R
definida en algn intervalo I de R, solucin nica del problema de valor inicial
de Cauchy.

Teorema 2. Sea F : R tal que F,

Ntese que si en (2), F es un polinomio en la variable x, entonces la ED


admite solucin nica (t, x) R2 . En lo que sigue, siempre supondremos existencia y unicidad de los problemas de Cauchy.
Una solucin x(t) de x = F (t, x) se representa geomtricamente en el plano
tx. Como dado cualquier (t0 , x0 ) D existe una solucin que pasa por (t0 , x0 ),
ello quiere decir que las soluciones de la ED se representan por una familia de
curvas solucin en D y que ! curva solucin que pasa por un punto determinado.
Ejemplos
1. x = x
}t
| {z
f (x,t)

dx
Cualitativamente
= 0 x(t) al- Analticamente
dt
3
canza sus puntos crticos en xt = 0 Mtodo 1 x x = t ED lineal
 Z

R
R
dt
dt
2
x(t) = e
te
dt + C
2

= et t et + et +C
1
1
= t + 1 + C et
0
2
-1
-2
-2

Mtodo 2 u = x t

2
-1

du
dx
=
1
dt
dt
du
du

+1=u
=0
dt
u1
ln(u 1) = t + k

u(t) 1 = C et

x(t) t = C et +1
x
dx
dt
2. x = , t 6= 0
=
t
x
t
ln |x| = ln |t| + C
1 dx
x
|x(t)| = k
= 0 = 0 x = 0 (solucin trivial)
t dt
t
Pero
x > 0 t > 0 x < 0
x > 0 t < 0 x > 0
x < 0 t > 0 x > 0
x < 0 t < 0 x < 0
2

t
3. x = , x 6= 0 x dx = t dt x2 + t2 = C
x

4. Autnoma: x =

1 2
(x 1) x = 0 x(t) = 1 o x = 1 + et/2
2

2
1

1
2

Vemos entonces que no es necesario tener ni la solucin explcita ni un dibujo


exacto para poder describir cualitativamente el comportamiento de la solucin
dx
de una ED. Ms an, si
= F (x) es una ED autnoma, las soluciones de ella
dt
tienen la siguiente propiedad:
Si : R R es una solucin, entonces cualquier traslacin de la grfica de
en la direccin del eje t es la grfica de otra solucin de la ecuacin.
En efecto, sea : R R con (t) = (t c), c R constante
c) F ((t c)) F ((t))
= (t
tambin es una solucin

Ejemplo x = x

2
1

1
2

Observacin 2. Notemos que el comportamiento est determinado por f (x)


1. Si c R : f (c) = 0, entonces x(t) = c es una solucin
4

2. Si f (x) 6= 0, entonces x(t) es creciente o decreciente, dependiendo del


sgnf (x)
Esta informacin puede representarse en una recta real que representa a x.

2.

Anlisis analtico vs cualitativo


dy
dy
= 4y(1 y)
= dt
dt
4y(1 y)

Podramos encontrar la solucin tras la integracin, pero las frmulas obtenidas


no son fciles de interpretar. Por ello, los mtodos geomtricos y cualitativa permiten obtener mucha informacin con poco trabajo. Esto es importante sobre
todo para ED difciles.

Considerar
2
dy
= ey /10 sen2 y
Z  dt
Z
1
2

ey /10 sen2 y
dy = dt

Pero cualitativamente
f (y) = ey

/10

sen2 y > 0

excepto para y = n que nos dan las soluciones de equilibrio.

2
1
9 8 7 6 5 4 3 2 1
1

2
La ecuacin del ejemplo tiene slo tres soluciones esencialmente diferentes:
crecientes, decreciente y constante.
5

Este resumen de comportamientos cualitativos se puede ilustrar geomtricamente definiendo una lnea de fases

Definicin 2. Sea x = F (x), y sea x0 R : F (x0 ) = 0. Entonces diremos que


x0 es una singularidad o un punto de equilibrio de la ED
Notar que x(t) = x0 es una solucin (cte.) de la ED, pues
x(t)

= 0 y F (x0 ) = 0
Observacin 3. La lnea de fase se puede obtener directamente de x(t)

= F (x),
simplemente mirando cuando F (x) 0.
Ejemplo Supongamos que Gr F viene dado por

Entonces, la lnea de fase

en donde la direccin de la flecha indica si la funcin crece () o decrece ()


Observacin 4. Dependiendo del comportamiento local de las soluciones de la
ecuacin de una singularidad aislada en la Lnea de Fases, se distinguen los
siguientes tipos de singularidades
1. Repulsor o fuente

2. Atractor o sumidero
3. Atractor-repulsor
4. Repulsor-atractor
Ejemplo. Describir cualitativamente los tipos de soluciones que admite la
4
EDO x = (x 1)3 ecos(x 1) .
Solucin. Es inmediato que x = 1 es la nica singularidad (aislada) de la
4
4
EDO. Como ecos(x 1) > 0, x R, la grfica de F (x) = (x 1) ecos(x 1) , la
lnea de fases y los tipos de soluciones de la EDO, son como ilustra la figura,
respectivamente.
La ecuacin del ejercicio anterior, es muy fcil de integrar bajo el punto de
vista del clculo. En consecuencia es difcil conocer explcitamente las ecuaciones
de las soluciones que indica la figura. Sin embargo, la lnea de fases permite dar
una descripcin cualitativa del tipo de soluciones y punto de equilibrio.
Definicin 3. Dos ecuaciones diferenciales autnomas, se dicen cualitativamente equivalentes si, y slo si, tienen la misma lnea de fases, en el sentido del
mismo nmero de puntos singulares, de la misma naturaleza y distribuidos en
el mismo orden.
Ejemplos
4

1. Las ED x = x, x = (x 1) ecos(x 1) , son equivalentes, pues solo tiene un


nico punto atractor en sus respectivas lneas de fases.
2. x = (x+2)(x+1) x = 12 (x2 1). Pero x = (x+2)(x+1) 6 x = 12 (x2 1)
(pues en este ltimo caso el atractor y el repulsor se encuentran en distinto
orden).

3.

Sistemas de ED de primer orden (en el plano)


Consideremos el sistema
x
y

= f (t, x, y)
= g(t, x, y)

que podemos expresar matricialmente en la forma


dY
= F (t, Y ),
dt
donde
Y =

 
x
y

dY
=
dt

 
x
y

F (t, Y ) =



f (t, x, y)
g(t, x, y)

Una condicin inicial para este sistema viene dada por


x(t0 ) = x0
y(t0 ) = y0

 
x0
Y (t0 ) =
= Y0
y0

Una solucin de este sistema es una matriz Y (t) que satisface el sistema.

Ejemplo. Comprobar que Y (t) = e4t 3 e2t , 3 e4t +3 e2t es una
solucin de


dY
x + y
=
3x 5y
dt

Solucin. Derivar y reemplazar.


x(t) = e4t 3 e2t

y(t) = 3 e4t +3 e2t

Definicin 4.

1. Un sistema de la forma

dY
= F (Y ) se llama autnomo
dt

dY
= F (Y ) si F (Y0 ) = 0. La funcin
dt
Y (t) = Y0 t es una solucin de equilibrio.

2. Y0 es un punto de equilibrio para

3. El plano de fase (equivalente a lnea de fase pero en dos dimensiones)


determina el comportamiento cualitativo de x(t) e y(t), si t crece
Ejemplos
1.
2.


dY
f (x, y)
Sea Y (t) = (x(t), y(t)) solucin de
=
. La curva C : x =
g(x, y)
dt
x(t), y = y(t) tiene la propiedad que en el punto (x, y) C, el vector tangente
a la curva en (x, y) es (f (x, y), g(x, y)).
Esto sugiere dibujar el campo de vectores del sistema, donde cada punto
P = (x, y) del plano se dibuja el vector (f(x,y),g(x,y)) con punto inicial P .
Problema. Los vectores varan en magnitud, esto dificulta el dibujo, por lo
tanto, se utiliza el campo de direcciones, se dibujan vectores de igual magnitud,
(f (x, y), g(x, y))
por ejemplo,
si (f (x, y), g(x, y)) 6= (0, 0)
k(f (x, y), g(x, y))k

Teorema 3. Sea

dY
= F (t, Y )
dt
Y (t0 ) = Y0

un problema de valor inicial


F (t, Y ) =



f (t, x, y)
g(t, x, y)

Si F es de clase C 1 , es decir, si f y g son de clase C 1 , entonces > 0 y una


funcin Y (t) definida para t ]t0 , t0 + [ tal que Y (t) satisface el problema
de valor inicial. Adems esta solucin es nica.
dY
= F (Y ), con F de clase
dt
1
C , entonces, por el teorema, existen nicas dos funciones distintas y1 (t) e y2 (t)
no pueden intersectarse.

Observacin 5. Si el sistema es autnomo, es decir,

Supongamos que y1 (t) e y2 (t) se intersectan en el punto Y0 , es decir Y1 (t1 ) =


Y0 = Y2 (t2 ). Como el campo vectorial F (Y ) no cambia con el tiempo, obtenemos
las mismas curvas solucin de dos soluciones diferentes que empiezan en el mismo
punto Y0 en tiempos diferentes Y1 (t1 +t) = Y2 (t2 +t) t. No ocurre lo mismo
para sistemas no autnomos.

Ejemplos de sistemas de ED de primer orden


Ejemplo. (Sistema depredador-presa) Sea R(t) la poblacin de presas presentes en el instante t. F (t) la poblacin de depredadores presentes en el instante
t. Se supone R(t), F (t) > 0
dR
dt
dF
dt

= 2R 1,2RF
= F + 0,9RF

El trmino 2R representa el crecimiento exponencial de las presas en ausencia


de depredadores y 1,2RF corresponde al efecto negativo sobre la presa por la
interaccin de ambas especies.
F corresponde a la hiptesis de que los depredadores mueren si no hay presas
que comer y 0,9RF es el efecto positivo sobre los depredadores por la misma
interaccin.
Solucin de equilibrio:
(2 1,2F )R =
(1 + 0,9R)F =

0
0

10
5
(R(t) = 0, F (t) =) (R(t) =
, F (t) = ) el sistema est en
9
3
equilibrio perfecto.
Computador. Campo de direcciones y algunas soluciones.
Ejemplo. Otro modelo
dR
dt
dF
dt



R
= 2R 1
1,2RF
2
= F + 0,9RF

Capacidad de soporte de R es 2.
Ejemplo. Movimiento de una masa unida a un resorte

y=0
posicin de reposo

Figura 2

Una masa unida a un resorte se desliza sobre una mesa sin friccin.
Problema. Calcular su movimiento horizontal cuando el resorte se estira (o
comprime) y luego se suelta.
Suponer que lo nico que acta sobre la masa es la fuerza del resorte
Ley de Hooke: La fuerza restauradora ejercida por un resorte es linealmente proporcional al desplazamiento del resorte desde su posicin de reposo y
est dirigida hacia esa misma posicin.
F = ky

d y
= ky
dt2

k = cte. del resorte


d2 y
k
+
y=0
2
dt
m
oscilador armnico

10

(ecuacin diferencial de 2.o orden)

Problema de valor inicial: se necesitan dos condiciones iniciales y(0) = y0


y (0) = v0 (pues, estirarlo y liberarlo 6= estirarlo y despus empujarle con una
cierta velocidad inicial). Sea v(t) =velocidad de la masa en el tiempo t. La
ecuacin de 2.o orden se puede escribir como un sistema:
dy
dt
dv
dt

k
y
m

Analizar el problema de valor inicial, con


d2 y
dt2

k
m

= 1:
dy
dt
dv
dt

y(0) =

y (0) =

(1)

(y(0), v(0)) =

(1, 0)

(2)

Solucin de (1): y(t) = cos t.


Solucin de (2): (y(t), v(t)) = (cos t, sen t).
En el plano yv la solucin es una circunferencia con centro en el origen y
radio 1, es decir, la solucin es peridica, la masa oscila para siempre de ida ya
vuelta a travs de su posicin de reposo.

4.

Sistemas desacoplados

Definicin 5. Un sistema de ED est desacoplado si la razn de cambio de una


ms de las variables dependientes, depende slo de su propio valor.
Ejemplo.
dx
dt
dy
dt

2x

11


Solucin general (x(t), y(t)) = k1 e2t , k2 et . Si se da la condicin inicial
y(0) = (1, 1) se tiene k1 = 1 k2 = 1 y(t) = e2t , et es la solucin del
problema de valor inicial y(t) es una parametrizacin de x = y 2 con y > 0.
x, y

y(t)
x(t)

5.

Sistema parcialmente desacoplado


dx
dt
dy
dt

2x + 3y

4y

Solucin para y: y(t) = k2 e4t . Sustituyendo en la primera ecuacin:


dx
= 2x + 3k2 e4t
dt

ecuacin diferencial lineal de primer orden


Z

2t
2t
4t
x(t) = e
e
3k2 e
dt + K
1
x(t) = k1 e2t k2 e4t
2
1
2t

x(t) = k1 e k2 e4t
2
y(t) = k2 e4t

Si ponemos la condicin inicial x(0) = 0, y(0) = 1, se tiene que


k1

1
2

k2
k2

=
=

0
1

k1 =

1
k2 = 1
2

1 2t 1 4t
e e
2
2
y(t) = e4t

x(t) =

12


Si la condicin inicial es (x(0), y(0)) = 21 , 1 k1 = 0 k2 = 1.
x(t) =

1 4t
e
2

y(t) = e4t
y(t)
= 2 t la solucin se encuentra
x(t)
sobre una recta que pasa por el origen y que tiende a (0, 0) cuanto t y
crece cuando t .
En este ltimo caso observemos que

y(t)

x(t)

6.

Sistemas Lineales
dx
dt
dy
dt

ax + by

cx + dy

sistema lineal (en el plano) con coeficientes constantes

(autnoma)

6.1.

Notacin matricial


a
Sea A =
c


 
b
x
(matriz de los coeficientes) e Y =
, el sistema queda
d
y

dY
= AY .
dt
Anlogamente para un sistema n-dimensional.
Observacin 6. Si det A 6= 0, el nico punto de equilibrio para el sistema es el
origen. En este caso se dice que el sistema es no degenerado.

6.2.
Sea

Principio de linealidad
dY
= AY
dt

1. Si Y (t) es una solucin kY (t) es una solucin k R.


2. Si Y1 (t)Y2 (t) son dos soluciones Y1 (t)+ Y2 (t) es tambin una solucin.
13

Consecuencia Y1 (t) Y2 (t) son soluciones k1 Y1 (t) + k2 Y2 (t) es solucin


k1 , k2 R


dY
2 3
Por ejemplo,
=
Y.
0 4
dt
Ya sabemos que
 2t 
 4t 
e
e
Y1 (t) =
Y2 (t) =
0
2 e4t
entonces
Y1 (t) + Y2 (t) =

 2t

e e4t
2 e4t

es tambin una solucin.

y1 (t) + y2 (t)
y2 (t)

y1 (t)
Figura 3

dY
= AY (A
dt
M (2 2)). Si Y1 (0) e Y2 (0) son l.i., entonces para cualquier condicin inicial
Y (0) = (x0 , y0 ), k1 , k2 R tales que k1 Y1 (t) + k2 Y2 (t) es la solucin del problema de valor inicial
 
dY
x0
= AY, Y (0) =
y0
dt

Teorema 4. Sean Y1 (t) e Y2 (t) dos soluciones del sistema lineal

Demostracin. {Y1 (0), Y2 (0)} l.i. {Y1 (0), Y2 (0)} es base de R2 !k1 , k2 R
tales que
 
x0
k1 Y1 (0) + k2 Y2 (0) =
y0
 
x0
Y (t) = k1 Y1 (t) + k2 Y2 (t) es solucin del sistema que satisface Y (0) =
.
y0
Observacin 7.
1. Si Y1 (t) e Y2 (t) son dos soluciones del sistema tales que
Y1 (0) e Y2 (0) son l.i. Y1 (t) e Y2 (t) son l.i.

14

2. Teorema el conjunto de las soluciones del sistema es un espacio vectorial


de dimensin dos, donde B = {Y1 (t), Y2 (t)} es una base. En este caso, se
dice que k1 Y1 (t) + k2 Y2 (t) con k1 , k2 R se llama la solucin general del
sistema df racdY dt = AY .
Ejemplos
 2t 
 4t 
dY
e
e
=
1. Se sabe que Y1 (t) =
e Y2 (t) =
son soluciones de
4t
0
2
e
dt


 
2 3
3
. Determinar la solucin Y (t) tal que Y (0) =
.
0 4
4
 
 
1
1
Solucin. Y1 (0) =
e Y2 =
{Y1 (0), Y2 (0)} l.i. {Y1 (t), Y2 (t)}
0
2
base del espacio solucin.
 
   
1
1
3
k1
+ k2
=
0
2
4
k1 k2 = 3
k1 = 1 k2 = 2
2k2 = 4
 2t

e
2 e4t
Y (t) =
es la solucin pedida
4 e4t

d2 Y
2. Oscilador armnico:
= y. Es fcil ver que y1 (t) = cos t e y2 (t) =
dt2
sen t son soluciones de la ecuacin diferencial. El sistema asociado es:
dy
=
v


dt
dY
0 1

=
Y
dv
1 0
dt
= y
dt

Y1 (t) =

 


 

y1 (t)
cos t
y2 (t)
sen t
=
Y2 (t) =
=
v1 (t)
sen t
v2 (t)
cos t

son soluciones del sistema.


 
 
1
0
Y1 (0) =
Y2 (0) =
son l.i.
0
1




cos t
sen t
Y (t) = k1
+ k2
solucin general del sistema
sen t
cos t
y(t) = k1 cos t + k2 sen t

k1 , k2 R es solucin general de la ecuacin diferencial.

15

Observacin8. A veces
en lneas
 los sistemas tienen
 2tsoluciones

 rectas.
Por ejemdY
2 3
e
e4t
plo:
=
Y tiene Y1 (t) =
Y2 (t) =
solucin
0 4
0
2 e4t
dt
en lnea recta:

y2 (t)

y1 (t)
Figura 4
 
0
Y1 (t) est sobre el eje x: lm y1 (t) lm Y1 (t) =
. Y2 (t) =
0
t
t
 
 
1
1
e4t
. Y2 (t) est sobre la recta de direccin
(en el II cuadrante)
2
2
 
0
lm Y2 (t) =
.
0
t
Pero no todos los sistemas tiene solucin en lnea recta, el ejemplo 2 del
oscilador armnico no las tiene.
Cmo determinar si las tiene?
6.2.1.

Soluciones en lnea recta


dY
= AY . Si V = (x, y) est sobre una solucin en lnea recta, entonces
Sea
dt
el campo vectorial en (x, y) debe apuntar en la misma direccin o en la direccin
opuesta que el vector de (0, 0) a (x, y). Por lo tanto debemos buscar V no nulos
tales que AV = V . Es decir, debemos determinar los valores propios de A.
Supongamos que ???A y V es vector propio asociado a . Sea Y (t) =
et V .
dY
= et V = Y (t)
dt
Por otra parte
AY (t) = A et V = et AV = et V = Y (t)
Y (t) = et V

es solucin del sistema.


16

No se entiende

Teorema 5. Sea valor propio de A M (2, R) con V un vector propio asocidY


ado, entonces el sistema lineal
= AY tiene la solucin en lnea recta:
dt
Y (t) = et V.
Si adems A tiene dos valores propios distintos 1 , 2 con V1 , V2 vectores propios
asociados, respectivamente, se tiene que Y1 (t) = e1 t V1 Y2 (t) = e2 t V2 son
dos soluciones en lnea recta.
Adems Y1 (t) e Y2 (t) son l.i.
Y (t) = k1 e1 t V1 + k2 e2 t V2

k1 , k2 R

es solucin general del sistema.


Ejemplos


dY
2 1
1. Sea
= BY =
Y , pB () = 2 5 + 4 = ( 4)( 1), B
2 3
dt
tiene dos valores propios 1 = 4 2 = 1.

   
2 1
x
0
V1 vector propio asociado a 1 :
=
2 1
y
0

V2

2x + y = 0 y = 2x
 
1
V1 =
2

   
1 1
x
0
vector propio asociado a 2 :
=
x+y =0
2 2
y
0
 
1
V2 =
1
 
 
1
4t 1
t
k1 , k2 R
Y (t) = k1 e
+ k2 e
2
1

es la solucin general del sistema.

17

2.

d2 y
dy
+7
+ 10y = 0. El sistema correspondiente es:
dt2
dt

 
 
dY
v
0
1
y
=
=
= CY
7v 10y
10 7
v
dt

pC () = 2 + 7 + 10 = ( + 2)( + 5) 1 = 2 2 = 5

   
2
1
x
0
V1 :
=
10 5
y
0
 
1
2x + y = 0 y = 2x V1 =
2

   
5
1
x
0
V2 :
=
10 2
y
0
 
1
y = 5x V2 =
5
 
 
1
1
Y (t) = k1 e2t
+ k2 e5t
5
2

y(t) = k1 e

7.

2t

k1 , k2 R
+k2 e

5t

solucin general del sistema

k1 , k2 R

solucin general de la ecuacin diferencial.

Retrato de fase para sistemas lineales con valores propios reales

dY
= AY, A M (2, R) y det A 6= 0
dt
la nica singularidad es el origen.
Si A tiene dos valores propios 1 , 2 R {0}, pueden ocurrir los casos
siguientes:
1. 1 6= 2 .
2. 1 = 2 = .
Caso 1:

1 6= 2

1 < 0 < 2 En este caso se dice que el origen es un punto silla. Hay dos
lneas que corresponden a la solucin en lnea recta. A lo largo de una lnea las
soluciones tienden a (0, 0) cuanto t + y a lo largo de la otra, se alejan de
(0, 0) cuanto t +. Las otras soluciones llegan del y se van al infinito.
Sean Vi , i = 1, 2 vectores propio asociados a i . Soluciones en lnea recta
Y1 (t) = k1 e1 t V1 Y2 (t) = k2 e2 t V2 . Solucin general
Yg (t) = k1 e 1 tV1 + k2 e2 t V2
k1 , k2 R
lm Y1 (t) = (0, 0) lm Y2 (t) = (0, 0)
t+

18

V1

V2

Figura 5

Ejemplos
dY
1.
= AY donde A =
dt


2
0


0
3

1 = 2, = 3,
Y (t) = k1 e

2t

V1 =

 
 
1
0
, V2 =
0
1



k1 e2t
V1 + k2 e V2 =
k2 e3t
3t

Sea Y (t) tal que k1 , k2 6= 0. Sea |t| grande:


a) si t > 0, Y (t) k2 e3t V2 Y (t) tiene al eje y como asntota.
19

b) si t < 0, Y (t) k1 e2t V1 Y (t) tiene al eje x como asntota.


 2t 
e
Por ejemplo, sea Y (t) tal que Y (0) = (1, 1), entonces Y (t) =
=
e3t


x(t)
y(t)
lm x(t) = 0

t+

lm y(t) = +

t+

lm x(t) = +

lm y(t) = 0

y(t)

x(t)

dY
2.
= BY , donde B =
dt


8
6


11
9

pB () = 2 + 6 = ( + 3)( 2) 1 = 3 2 = 2
 
 
1
11
V1 =
V2 =
1
6
Solucin en lnea recta Y1 (t) = e3t V1 Y2 (t) = e2t V2 . Solucin
general Y (t) = k1 e3t V1 + k2 e2t V2 , k1 , k2 R. Si k1 , k2 6= 0, cuando
t +, Y (t) tiene como asntota la recta de direccin V2 . Cuando t
, Y (t)
2 < 1 < 0
Y (t) = k1 e1 t V1 + k2 e2 t V2 ,

k1 , k2 R

Cuando t +, todas las soluciones tienden a (0, 0) nica singularidad del


sistema, el (0, 0) es un pozo o resumidero.
Sea Y (t) tal
 que
 k1 , k2 6=
 0Cmo se acerca Y (t) al origen?
a
c
Sean V1 =
y V2 =
, con ad bc 6= 0
b
d
x(t) = ak1 e1 t +ck2 e2 t
y(t) = bk1 e1 t +dk2 e2 t
20

V1

V2

Figura 6

Sea a 6= 0
dy
dy
1 bk1 e1 t +2 dk2 e2 t
1 bk1 + 2 dk2 e(2 1 )t
= dt =
=

t
1
2
dx
dx
1 ak1 e +2 ck2 e
1 ak1 + 2 ck2 e(2 1 )t
dt
dy
b
lm
= .
t+ dx
a
Es decir, Y (t), cuando t +, tiende al origen con pendiente que tiende a ab ,
es decir, en direccin tangencial a la lnea de vectores propios correspondientes
al valor propio 1 .
0 < 1 < 2
Y (t) = k1 e1 t V1 + k2 e2 t V2

k1 , k2 R

Si (k1 , k2 ) 6= (0, 0), Y (t) se aleja del origen cuanto t +, es decir, el


origen es una fuente.
Todas las soluciones con k1 k2 6= 0, se acercan al origen cuanto t + en
forma tangencial a la lnea de vectores propios correspondientes al valor propio
1 .


 
dY
1 3
1
Por ejemplo
= AY donde A =
, 1 = 1, 2 = 2. V1 =
, V2 =
0
2
0
dt
 
3
1
21

Caso 2: 1 = 2 = 6= 0
Debemos distinguir dos subcasos:
1) Si V = R2 (V subespacio propio asociado a )
2) Si dim V = 1
1)
Y (t) = et (k1 V1 + k2 V2 )
k1 , k2 R
 
a
= et
a, b R
b

(a) modo estelar atractor o


sumidero estelar si < 0

(b) modo estelar repulsor o


fuente estelar si > 0

dY
Forma general de la solucin Sea
= AY, A M (2, R), donde A
dt
tiene un valor propio de multiplicidad dos, pero dim V = 1. Si V1 es vector
propio:
Y1 (t) = et V1 es una solucin en lnea recta
Conjetura
Y2 (t) = et (tV1 + V2 )
22

es solucin, donde V2 es un vector por determinar

Si Y2 (t) es solucin

dY2
= AY2
dt

t. Pero

dY2
= et (tV1 V2 ) + et V1 = et tV1 + et V1 + et V2
dt
AY2 = et (tAV1 + AV2 ) = et (tV1 + AV2 )
Debemos tener
et (tV1 + V1 + V2 ) = et (tV1 + AV2 )
V1 + V2 = AV2
(A I)V2 = V1
sistema lineal en V2 que tiene infinitas soluciones
Y (t) = k1 e1 t V1 + k2 et (tV1 + V2 ) es solucin general.

Ejemplo.
dy
dt
dv
dt

=
=

dy
d2 y
+2 2
+ 2y = 0
2
dt
dt
v

dY

dt
2y 2 2 v


0
1

Y
2 2

 
y
con Y =
v

2 + 2 2 + 2 = 0 = 2

   

2
1
x
0

=
2x + y = 0 y = 2x
V1 :
y
0
2 2





1
1
2t

V1 =
Y1 (t) = e
.
2
2
Determinacin de V2 :

2
2

  


1
1
x

=
2x+y = 1
y
2
2
1
x = (1 y)
2

23

Por ejemplo:
 
0
V2 =
1
Solucin general:
 

  

1
1
0

+ k2 e 2t t
+
1
2
2



k1 +tk2
k1 , k2 R
Y (t) = e 2t
2k1 2tk2 + k2

Y1 (t) = k1 e

2t

Definicin 6. Un punto de equilibrio (x0 , y0 ) es estable, si las soluciones que


parten de puntos situados en una vecindad de (x0 , y0 ) tienden a (x0 , y0 ) cuando
t +, en caso contrario son inestables.
Por ejemplo, los sumideros son estables; las sillas y fuentes son inestables.

8.

Sistemas en que pA () tiene races complejas

Consecuencia El sistema no tiene soluciones en lnea recta. Por ejemplo:




dY
a b
=
Y,
b>0
b a
dt
Usando coordenadas polares: x = r cos , y = r sen (r > 0) se tiene
x =
y =

r
r cos r sen

r sen + r cos

=
=

x cos + y sen
1
sen )
r (y cos x

r = cos (ar cos br sen ) + sen (br cos + ar sen ) = ar


1
= (cos (br cos + ar sen ) sen (ar cos br sen )) = b
r

24

El sistema en coordenadas polares es:


r

=
=

ar
b

con b > 0

Su solucin es:
r(t)
(t)

=
=

r(0) eat
(0) + bt

Si a = 0 r(t) = r(0) t las rbitas son circunferencias con centro en el


origen todas las soluciones son peridicas de perodo 2/b y frecuencia b/2.
En este caso se dice que el origen es un centro.

a b
Nota 1. Si A no es de la forma b
con b 6= 0 pero pA () tiene races
a
complejas con parte real nula, todas las soluciones son peridicas. Por ejemplo,
elipses.
Si a < 0, entonces lm r(t) = 0 y las rbitas son espirales que se enrollan
t+

en el origen. El origen es un foco atractor o sumidero en espiral.


Si a > 0, entonces lm r(t) = 0 y las rbitas son espirales que se alejan del
t

origen. El origen es un foco repulsor o fuente en espiral.


Resolucin de un sistema de este tipo (sin usar coordenadas polares)
dY
Sea
= AY tal que pA () tiene races complejas = a + bi, b 6= 0. Se
dt
considera A M (2, C) y se resuelve el sistema en C2 . Sea Y (t) una solucin
compleja del sistema, se tiene
Teorema 6. Supongamos que Y (t) = YRe (t) + iYIm (t), donde YRe (t) e YIm (t)
son las funciones reales, parte real y parte imaginaria de Y (t). Entonces YRe (t)
e YIm (t) son soluciones del sistema original.
Demostracin. Y (t) = YRe (t) + iYIm (t) es solucin de

dY
= AY en C2
dt

dY
dYRe
dYIm
=
+i
= AY (t) = AYRe (t) + iAYIm (t)
dt
dt
dt
dYRe (t)
dYIm (t)
= AYRe (t)
= AYIm (t)

dt
dt

Ejemplos


dY
2 3
1.
= AY con A =
, = 2 3i (origen es un pozo en
3 2
dt
espiral). Si V es vector propio complejo de 2 + 3i, entonces V es vector
25

propio complejo de 2 3i. Calculemos V :



   
3i 3
x
0
3ix 3y = 0
=

3
3i
y
0
3x 3iy = 0
 
i
x = iy V =
1
 
(2+3i)t i
Y (t) = e
es una solucin del sistema
1
 
i
2t
Y (t) = e (cos 3t + i sen 3t)
1


2t sen 3t + i cos 3t
=e
cos 3t + i sen 3t




2t sen 3t
2t cos 3t
YRe (t) = e
,
YIm (t) = e
cos 3t
sen 3t
YRe (t) YIm (t) son l.i. Solucin general del sistema
 



sen 3t
cos 3t
Y (t) = e2t k1
+ k2
k1 , k2 R
cos 3t
sen 3t



dY
0 1
2.
=
Y , = i 2 (origen es un centro).
2 0
dt





1
1
V =
Y (t) = (cos 2 t + sen 2 t)
i 2
i 2





cos 2t
sen
2t

YRe (t) =
, YIm (t) =
2 sen 2 t
2 cos 2 t
Solucin general del sistema:


 

2 t
x(t)
k1 cos 2 t + k2 sen

Y (t) =
=
y(t)
k1 2 sen 2 t + k2 2 cos 2 t
x(t) e y(t) son las ecuaciones paramtricas de una elipse (si k1 k2 6= 0).
En efecto:
2x2 + y 2 = 2k12 + 2k22 = k32
x2
y2
+
=1
k32
2k3
Para la direccin, notar que el sistema
dx
dy
=y
= 2x
dt
dt
Dibujamos el campo de direcciones
~v (1, 1) = (1, 2), ~v(0, 1) = (1, 0), ~v (1, 1) = (1, 2)
26

9.

Casos con valor propio cero

dY
Supongamos que
= AY , donde A tiene 2 valores propios 1 = 02 6= 0.
dt
Sean V1 V2 vectores propios asociados a 1 2 , respectivamente.
Solucin.
Y (t) = k1 e1 t V1 + k2 e2 t V2
Y (t) = k1 V1 + k2 e2 t V2

k1 , k2 R

Si k2 = 0, Y (t) = k1 V1 todos los puntos de la recta que pasa por el origen


con direccin V1 son puntos de equilibrio.
Si 2 < 0, Y (t) = k1 V1 + k2 e2 t V2 tiende a k1 V1 cuando t + a lo largo
de una recta k V2 .
Si 2 > 0, Y (t) se aleja de la lnea de puntos de equilibrio cuando t crece.


dY
3 1
Ejemplo.
= AY con A =
pA () = 2 + 4 = 0
3 1
dt
1 = 0, 2 = 4.

   
 
3 1
x
0
1
V1 :
=
3x + y = 0 y = 3x V1 =
3 1
y
0
3
 

   
1 1
x
0
1
V2 :
=
x = y V2 =
3 3
y
0
1
 
 
1
1
4t
Solucin. Y (t) = k1
+ k2 e
La recta : y = 3x es la lnea de
3
1
puntos de equilibrio. Cualquier otra solucin tiende a un punto de siguiendo
una semirecta k a y.

27

También podría gustarte