Apunte Edo Cualitativas
Apunte Edo Cualitativas
Apunte Edo Cualitativas
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.
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
Ejemplo x = x
2
1
1
2
2.
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
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
= 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
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
3.
= f (t, x, y)
= g(t, x, y)
x
y
dY
=
dt
x
y
F (t, Y ) =
f (t, x, y)
g(t, x, y)
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
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.
Teorema 3. Sea
dY
= F (t, Y )
dt
Y (t0 ) = Y0
f (t, x, y)
g(t, x, y)
= 2R 1,2RF
= F + 0,9RF
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
10
k
y
m
k
m
= 1:
dy
dt
dv
dt
y(0) =
y (0) =
(1)
(y(0), v(0)) =
(1, 0)
(2)
4.
Sistemas desacoplados
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.
2x + 3y
4y
x(t) = k1 e k2 e4t
2
y(t) = k2 e4t
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
(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
2t
e e4t
2 e4t
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
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
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
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.
No se entiende
k1 , k2 R
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
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
k1 , k2 R
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
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
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
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
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
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
8.
r
r cos r sen
r sen + r cos
=
=
x cos + y sen
1
sen )
r (y cos x
24
=
=
ar
b
con b > 0
Su solucin es:
r(t)
(t)
=
=
r(0) eat
(0) + bt
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
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.
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
27