Introducci´ on on a las ecuaciones diferenciales ordinarias
Noem´ No em´ ı Wolans Wolanski ki
α / β
γ / δ Diagrama de fases para dos poblaciones simbi´ oticas oticas
´ Indice General Preliminares
5
Cap´ Cap´ıtulo 1. Introducci´ on on 1. Generalidades. 2. Desc Descri ripci pci´ on o´n de algunos m´etodos etodos de resoluci´ on de ecuaciones de 1er. orden. Ejercicios
7 7 10 12
Cap´ Cap´ıtulo 2. Existencia y unicidad de soluci´ on Ejercicios
15 22
Cap´ Cap´ıtulo 3. Sistemas Sistemas lineales de 1er. orden y ecuaciones ecuaciones lineales de orden n 1. Generalidades y sistemas homog´eneos 2. Sistemas no homog´eneos
23 23 29
Cap´ Cap´ıtulo 4. Resoluci´ on de sistemas lineales con co eficientes constantes Ejercicios
33 45
Cap´ Cap´ıtulo 5. Ecuaciones lineales de orden n con coeficientes constantes Ejercicios
47 52
Cap´ Cap´ıtulo 6. Comportamiento asint´ asint´ otico de las soluciones 1. Diagramas de fases 2. Diagramas de fases de sistemas lineales a coefi oeficientes constantes 3. Linear Lineariza izaci´ ci´ on 4. Sistemas Conservativos Ejercicios
55 56 60 68 74 79
Agradecimientos
81
Bibliograf´ıa
83
3
Preliminares El objetivo de estas notas es dar una introducci´on on al tema tema de Ecuacio Ecuaciones nes Difere Diferenci nciales ales Ordina Ordinaria riass (en adelant adelantee ODE) ODE) a nivel nivel elemen elemental. tal. Las notas est´ an dirigidas a estudiantes de an la materia An´ alisis alisis II – Matem´ atica 3 de la Facultad de Ciencias Exactas y Naturales de la atica Unive Universi rsidad dad de Buenos Buenos Aires. Aires. Al dise˜ nar estas notas debemos tener en cuenta que en esta nar materia el tema de ODE se dicta en no m´as a s de 5 sema semana nas. s. Es por esta esta raz´ razon o´n que ciertos temas se dejan para desarrollar en los trabajos pr´ acticos. acticos. Entre Entre esos temas est´ an an los m´etodos eto dos de resoluci´on on de ecuaciones de primer orden y muchos ejemplos de aplicaciones que est´an como ejercicio para los alumnos. En estas notas discutiremos algunos problemas en los cuales aparecen ODE, daremos la demostraci´ on del Teorema de Existencia y Unicidad local de soluci´ on on y analizaremos el dominio on de definici´on on de las mismas. mismas. A fin de dar claridad claridad al texto, texto, daremos daremos las demost demostrac racione ioness bajo condiciones simples. Se dar´an an los m´etodos eto dos de resoluci r esoluci´ o´n de ecuaciones y sistemas lineales a coeficientes constantes on (tanto (ta nto homog´ homo g´eneos eneo s com comoo no homog´ hom og´eneos). eneo s). Por otro lado, se discutir´ a la noci´ on de diagrama de fases y su relaci´on on on con la posibilidad de predicci´on on del comportamiento de las soluciones sin conocer una f´ ormula ormula an´ alitica alitica de las mismas. mismas. Se ver´ vera´ c´ omo son los diagramas de fases de sistemas lineales a coeficientes constantes omo de dimensi´on on 2 y tambi´ tambi´en en para sistemas no lineales lineales conserv conservativ ativos. os. Se discutir´ discutir´ a la noci´ on o n de estabilidad lineal y se utilizar´ a para determinar la estabilidad de equilibrios de sistemas no lineales de dimensi´ on on 2.
5
CAP´ıTULO ıT ULO 1
Introducci´ on on 1. Generalid Generalidades. ades.
Sea V ( V (t,x,y,z) t,x,y,z ) un campo de velocidades velocidades correspondien correspondiente te a un fluido (por ejemplo). ejemplo). En el curso ya vimos que una part´ıcula ıcula que se mueve en el fluido sigue una trayectoria trayectoria σ (t) tal que su vector velocidad, σ (t), verifica σ (t) = V ( V (t, σ(t)) para todo tiempo t. Esto es un sistema de ecuaciones diferenciales de 1er. orden. A saber, si σ (t) = (x(t), y(t), z (t)) se debe tener para todo t,
(1.1)
x = V 1 (t,x,y,z) t,x,y,z ), y = V 2 (t,x,y,z) t,x,y,z ), z = V 3 (t,x,y,z) t,x,y,z ).
Claramente, Claramente, para determinar determinar la posici´ on on de d e una part´ıcula ıcula en un u n instante ins tante t debemos conocer tambi´ tamb i´en en su posici po sici´ on o´n en alg´ un un instante t0 ya que en un instante dado habr´a part´ pa rt´ıculas ıcul as en diferent dif erentes es puntos y por lo tanto sus trayectorias no son iguales. De modo que lo que nos plantearemos ser´a encontrar una soluci´ on o n de (1.1) sujeta a que 3 σ (t0 ) = X 0 donde t0 R y X 0 R son dados.
∈
∈
Por ejemplo, en una variable podr po dr´´ıamos intentar resolver el problema
x = x, x(0) = 1. 1.
x x (t) d Tenemos = 1, pero = log x(t). Por lo tanto, lo que queremos es que x x(t) dt d log x(t) = 1 dt
para to do t.
De aqu´ aqu´ı que se deba tener log x(t) = t + c para alguna constante c. Como para para t = 0 tenemos x(0) (0) = 1. Debe Debe ser ser log 1 = c. Esto Esto nos dice dice que que c = 0 y por lo tanto log x(t) = t o lo que es equivalente x(t) = et . Por otro lado, si tenemos
x = x, x(0) = a > 0, 7
´ 1. INTRODUCCION
8
la misma cuenta nos da log a = c. Por lo tanto, log x = t + log a x = et+log a = aet . Vemo emoss que a distin distintos tos datos datos inicial iniciales es le corres correspond ponden en distin distintas tas solucion soluciones es y adem´ adem´ as, si son distintas distintas en t = 0 son distintas para todo t. Verem eremos os m´ mas a´s adelante que este hecho es una propiedad general de las soluciones de ODE que dice que dos trayectorias de part´ part´ıculas diferentes no se cortan. Veamos otro ejemplo de sistema de ecuaciones. Supongamos Supongamos que tenemos una part´ part´ıcula de masa unitaria sujeta sujeta a un campo de fuerzas fuerzas F = (F 1 , F 2 , F 3 ). Entonces, como la fuerza es la masa por la aceleraci´on, on, si σ (t) es la trayectoria de la part´ part´ıcula, ıcula , se verifica σ (t) = F ( F (t, σ(t)) para to do t. Es decir,
x = F 1 (t,x,y,z) t,x,y,z ), y = F 2 (t,x,y,z) t,x,y,z ), z = F 3 (t,x,y,z) t,x,y,z ).
Ahora bien, si llamamos x0 = x , x1 = x , y0 = y , y1 = y , z0 = z , z1 = z . Enton Entonces ces,, obtenemos el siguiente sistema de primer orden:
x0 = x1 , x1 = F 1 (t, x0 , y0 , z0 ), y0 = y1 , y1 = F 2 (t, x0 , y0 , z0 ), z0 = z1 , z1 = F 3 (t, x0 , y0 , z0 ).
Este mismo enfoque permite tratar el caso en que el campo de fuerzas depende de la velocidad Esto es as´ as´ı cuando, cuando, por ejemplo, ejemplo, hay alg´ un un tipo de fricci´on on (como la resistencia del aire). aire). Esta Esta fuerza fuerza de fricci´ fricci´ on on es proporc proporcion ional al a la velocida velocidad d y con sentido sentido opuest opuesto. o. De modo que en general la fuerza ser´a de la forma F = F ( F (t,x,y,z,x , y , z ) y tendremos (reordenando las ecuaciones) x0 = x1 , (x , y , z ).
y0 = y1 , z0 = z1 , x1 = F 1 (t, x0 , y0 , z0 , x1 , y1 , z1 ), y1 = F 2 (t, x0 , y0 , z0 , x1 , y1 , z1 ), z1 = F 3 (t, x0 , y0 , z0 , x1 , y1 , z1 ).
Es decir, un sistema de ecuaciones de la forma ϕ = G(t, ϕ),
1. GENERALIDADES.
9
donde ϕ ahora ahora no es la traye trayecto ctoria ria de una part´ part´ıcula ıcula en el espacio espacio,, sino sino en lo que se llama llama el “Espacio de Fases” donde una fase es un par (σ, ( σ, σ ) donde σ = posici´ on on y σ = velocidad. En el espacio de fases ϕ es una trayectoria del campo G. De modo mo do que cualquier teor´ teor´ıa y cualquier informaci´on on que podamos recoger para sistemas de 1er orden, nos dar´ a informaci´ on on para sistemas de 2do. orden (mediante la reducci´ on descripta arriba). Pero ahora, si queremos on determinar la trayectoria σ de la part´ part´ıcula a partir de la trayectoria ϕ en el espacio de fases, necesitamos necesitamos datos iniciales iniciales para ϕ y ´estos est os son so n σ (t0 ) , σ (t0 ). Es decir, hay que dar la posici´on on y velocidad velocidad en un mismo tiempo t0 para obtener la trayectoria trayectoria de una part´ part´ıcula sujeta a un campo de fuerzas. fuerza s. Esto es bastante intuitivo intuiti vo desde el punto de vista vi sta f´ısico ısico dado que una part´ıcula ıcula sujeta s ujeta a un campo de fuerzas que empieza, digamos en el instante t = 0, en un cierto lugar, podr´ podr´ıa tener trayectorias distintas si originalmente su velocidad apunta en direcciones distintas. En general, si tengo una ecuaci´ on on de orden n: x(n) = f ( f (t,x,x , x ,
(1.2)
· · · , x(n−1)),
podemos reducirla a un sistema de n ecuaciones con n inc´ ognitas de la siguiente forma: Llamamos ognitas − ( n x0 = x , x1 = x , x2 = x , x3 = x , , xn−1 = x 1) . Median Mediante te este este proceso proceso (1.2) resulta resulta equivalente a
···
(1.3)
x0 = x1 , x1 = x2 , x2 = x3 , x3 = x4 , .. . xn −2 = xn−1 , xn −1 = f ( f (t, x0 , x1 , x2 ,
· · · , xn−1).
Luego, en este caso, vemos que tenemos que dar condiciones iniciales x(t0 ), x (t0 ), x (t0 ), x (t0 ), . . . , x(n−1) (t0 ), para determinar la trayectoria x(t). Un caso particula particularr de sistemas sistemas de ecuaci ecuacione oness de 1er. orden orden que result resultan an ser de especial especial n importancia importancia son los sistemas lineales, es decir aquellos sistemas sistemas X = V ( V (t, X ), ), X R en donde V es una funci´on on lineal de X para cada t y continua con respecto a t. Estos sistemas tienen la forma x1 = a11 x1 + a12 x2 + + a1n xn ,
∈
(1.4)
···
··· x2 = a21 x1 + a22 x2 + · · · + a2n xn ,
.. .
xn = an1 x1 + an2 x2 +
· · · + annxn.
Aqu´ı (x1 , x2 , , xn ) = X y la matriz (a (aij ) es la matriz asociada a la funci´ on on lineal V . V . Los aij son, en general, funciones continuas de la variable t. Cuando los coeficientes aij no dependen de t, decimos que se trata de un sistema lineal de 1er. orden con coeficientes constantes.
´ 1. INTRODUCCION
10
Uno de los motivos que da especial importancia al estudio de los sistemas lineales, es que los mismos pueden dar informaci´on on relevante sobre el comportamiento de las soluciones de sistemas m´ as generales (sistemas no lineales). Veremos esto con m´as as as detalle en el ´ultimo ult imo cap´ıtulo. ıtu lo.
2. Descripci´ Descripci´ on on de algunos m´ etodos etodos de resoluci´on de ecuaciones ecuaciones de 1er. orden. orden.
Para el caso particular de 1 ecuaci´on on de 1er. orden, existen varios m´ etodos etodos para hallar las soluciones. En estas notas s´olo olo mostraremos el m´ as sencillo de estos, que ya hemos usado, y es as el llamado m´etodo etodo de separaci´ separac i´ on de variables. variables. ¿En qu´e consiste? consis te? Supongamos que la ecuaci´ on on tiene la forma x = f ( f (x)g (t),
entonces, debe tenerse (si f ( f (x) = 0)
Sea F ( F (s) =
Sea ahora G(t) =
x (t) = g (t). f ( f (x(t))
ds 1 , es decir F (s) = . Entonces f ( f (s) f ( f (s)
d x (t) F ( F (x(t)) = F (x(t))x ))x (t) = . dt f ( f (x(t)) g (t) dt. dt. Entonces, d d F ( F (x(t)) = G(t). dt dt
De aqu aq u´ı que F ( F (x(t)) = G(t) + c (i.e. F ( F (x) = G(t) + c) y si podemos po demos despejar despe jar de d e aqu a qu´´ı x tendremos la soluci´ on on general x(t) dependiendo de una constante c a determinar por los datos iniciales. En la pr´actica, actica, la forma de escribir este razonamiento es como sigue: dx x = , dt por lo tanto, dx = f ( f (x)g (t). dt Llevamos todo lo que depende de x a la izquierda izquierda y lo que depende de t a la derecha operando con los diferenciales como si fueran n´ umeros. Entonces se obtiene umeros. dx = g (t) dt. f ( f (x) Integrando a ambos lados (olvidando que dependen de distinta variable, ya que son los diferenciales los que nos dicen respecto resp ecto de qu´e variable variable integramos)
dx = f ( f (x)
g (t) dt,
´ DE ALGUNOS M ETODOS ´ ´ DE 2. DESCRIPCION DE RESOLUCI ON DE ECUA ECUACI CION ONES ES DE DE 1ER. 1ER. ORD ORDEN EN..
11
o lo que es lo mismo F ( F (x) = G(t) + c. 1.1. Hallemos la soluci´ on on general de x = x2 . Aplicando el m´etodo etodo de separaci´ on on de variables, tenemos Ejemplo
dx = x2 dt
⇒
dx = dt x2
Si, por ejemplo, x(0) = 1 se tiene
⇒ − x1 = t + c ⇒
x=
− t +1 c .
−1 = c. Por lo tanto la soluci´on on es x=
1
1
− t.
Observemos que la soluci´on on hallada no est´a definida para todos los valores de t. El intervalo (que contiene al 0) en donde se encuentra definida la soluci´on on es ( , 1) y no puede extenderse 2 m´ as as all´a de d e ah a h´ı. En principio, princip io, de la ecuaci´ ecuacion o´n diferencial x = x , no hab´ıa ıa ning´ nin g´ un un elemento que 2 nos hiciera pensar que algo as´ as´ı podr´ıa ıa suceder, dado que la funci´ on on V ( V (x) = x es “tan buena” como uno quiere.
−∞
Este Este ejempl ejemploo nos dice dice que en el desarr desarrollo ollo de la teor teor´ıa genera generall no podemos podemos esperar esperar un resultado de existencia que nos diga que si V ( V (x) es regular entonces vaya a existir una soluci´on on definida para todo tiempo t. El resultado de existencia ser´a local . Ejemplo
1.2. Hallar, Halla r, por p or el m´etodo etod o de separaci´ separa ci´ on on de variables, ariables, las soluciones soluciones del problema problema x =
√x,
x(0) = 0. 0.
Aplicando el m´etodo etodo de separaci´ on de variables (suponiendo que x(t) = 0 para t > 0 para on poder dividir por x) obtenemos
√
√dxx = dt ⇒
√
2 x = t + c.
Como x(0) = 0 se sigue que c = 0 y por lo tanto 1 x(t) = t2 . 4 Pero observemos que, como x (0) = 0 podemos prolongar x a t < 0 como la funci´ on on id´ id´enti en ti-camente nula. Es decir, tenemos una soluci´on on dada por x(t) =
1 2 4t
0
si t > 0, si t 0.
≤
Por otro lado, si resolvemos resolvemos por p or el mismo m´etodo etodo este problema problema con dato inicial inicial 0 dado en un τ > 0 (es decir, pidiendo que x(τ ) τ ) = 0) y resolviendo para t > τ obtenemos 1 x¯(t) = (t 4
− τ ) τ )2
para t
≥ τ
´ 1. INTRODUCCION
12
y como antes antes podemos p odemos extender extender esta soluci´ on on a t < τ como la funci´ on on id´ i d´enticamente enticamente 0. Es decir, de cir, tenemos 1 (t τ ) τ )2 si t > τ, 4 x ¯(t) = 0 si t τ.
−
≤
Observemos que ambas funciones son soluci´on on de la ecuaci´on on x =
√x y satisfacen x(0) = 0.
Vemos con este ejemplo que no siempre existe una ´unica unica soluci´ on al problema de valores on iniciales. iniciales. Para Para obtener obtener unicidad, unicidad, que desde el punto punto de vista de las aplicaciones aplicaciones f´ f´ısicas es una propiedad importante, habr´ a que pedir hip´otesis otesis adicionales sobre la funci´ on on f ( f (t, x) del segundo miembro de la ecuaci´ on on que garanticen garanticen la unicidad. Observ Observar ar que f ( f (x) = x no es regular en x = 0. En el pr´oximo oximo cap´ cap´ıtulo estableceremos condiciones que aseguran unicidad de soluci´ on. on.
√
Ejercicios
(1) Para cada una de las ecuaciones diferenciales que siguen, encontrar la soluci´on on general general y, en los casos que se indica, la soluci´on on particular que satisfaga la condici´ on on dada: (a)
x
(c) x
1+x = 1+t
− 2tx = t, (e) x − x1/3 = 0,
(b) x(1) = 0
x
(d) x
1 + x2 = , 1 + t2
x(0) = 1
− tan t = cos t
x(0) = 0
En todos los casos dar el intervalo maximal de existencia de las soluciones y decir si son unicas. u ´ nicas. (2) Si y = y (t) denota el n´ umero de habitantes de una poblaci´ umero on on en funci´on on del tiempo, se denomina tasa de crecimiento de la poblaci´ on on a la funci´on on definida como el cociente y /y. /y. (a) Caracterizar (encontrar la ecuaci´ on) de las poblaciones con tasa de crecimiento on) constante. (b) Dibujar el gr´afico afico de y (t) para poblaciones con tasa de crecimiento constante, positiva y negativa. (c) ¿Cu´ ales son las poblaciones con tasa de crecimiento nula? ales (d) Una poblaci´ on on tiene tasa de crecimiento constante. El 1 de enero de 1992 ten´ıa ıa 1000 individuos, y cuatro meses despu´ es es ten´ ten´ıa 1020. Estimar el n´ umero umero de individuos que tendr´a el 1 de enero del a˜no no 2010, usando los resultados anteriores. (e) Caracterizar las poblaciones cuya tasa de crecimiento es una funci´on on lineal de t (at + b). (f) Caracterizar las poblaciones cuya tasa de crecimiento es igual a r cy, cy , donde r y c son constantes constantes positivas. positivas. Este es el llamado crecimiento crecimiento log´ log´ıstico, en tanto tanto que el correspondiente a tasas constantes es llamado crecimiento exponencial (por razones razones obvias ¿no?). ¿no?). Para poblaciones poblaciones peque˜ nas, ambas formas de crecimiento son muy similares similares.. Comprob Comprobar ar esta esta afirmac afirmaci´ i´ on on y comprobar tambi´ en en que en el crecimi crec imiento ento log´ıstico ıst ico y(t) tiende asint´oticamente oticamente a la recta y = r/c. r/c.
−
EJERCICIOS
13
(3) Si un cultivo de bacterias crece con un coeficiente de variaci´on on proporcional a la cantidad existente y se sabe adem´as as que la poblaci´on on se duplica en 1 hora ¿Cu´anto anto habr´ a aumentado en 2 horas?. (4) Verifique que las siguientes ecuaciones son homog´eneas eneas de grado cero y resuelva: x+t (a) tx = x + 2t 2t exp( x/t) x/t) (b) txx = 2x2 t2 (c) x = , x(1) = 0 t
−
−
(5) Demuestr Demuestree que la sustituci´ sustituci´ on on y = at + bx + c cambia x = f ( f (at + bx + c) en una ecuaci´on on con variables separables y aplique este m´etodo etodo para resolver las ecuaciones siguientes: (a) x = (x + t)2
(b) x = sen2 (t
− x + 1)
(6) (6) (a) (a) Si Si ae = bd demuestre que pueden elegirse constantes h, k de modo que las sustituciones t = s h, x = y k reducen la ecuaci´on: on:
−
−
dx at + bx + c = F dt dt + ex + f
a una ecuaci´on on homog´ hom og´enea. ene a. (b) Resuelva las ecuaciones: x+t+4 a) x = x+t 6
b) x =
−
x+t+4 t x 6
− −
(7) Resuelv Resuelvaa las ecuaciones ecuaciones siguientes: siguientes: (x + y3 )dy = 0 − x3)dx + (x (c) (3x (3x2 − y2 ) dy − 2xy dx = 0
(b) (b) cos cos x cos2 y dx
(a) (y (y
− 2sen x sen y cos y dy = 0
(d) x dy = (x5 + x3 y 2 + y) dx
(e) 3y 3y dx + x dy = 0
→R son continuas,
(8) Si a, b : [x1 , x2 ] (a) probar:
(i) y(x) = ke
−
x x1
a(t)dt
(k
∈ R) son todas las soluciones de
y + a(x)y = 0 en en [x [x1 , x2 ] (ii) y(x) =
−e
−
x x1
a(t)dt
x
b(t)e
t x1
a(s)ds
dt es una soluci´ on on de
x1
y + a(x)y + b(x) = 0 en [x [x1 , x2 ] (b) describir todas las soluciones de: y + a(x)y + b(x) = 0 en [x [x1 , x2 ] (c) Comprobar que y1 R existe una unica u ´ nica soluci´ on on de la ecuaci´ ecuacion o´n y + a(x)y + b(x) = 0 en [x [x1 , x2 ] tal que y(x1 ) = y1 , y que cualquier soluci´on on de la ecuaci´ ecuaci´ on on homo ho mog´ g´enea ene a y + a(x)y = 0 que se anule en un punto x [x1 , x2 ] es id´enticamente enticame nte nula.
∀ ∈
∈
(9) Hallar la ecuaci´ on de una curva tal que la pendiente de la recta tangente en un punto on cualquiera es la mitad de la pendiente de la recta que une el punto con el origen.
´ 1. INTRODUCCION
14
(10) Hallar la ecuaci´ on de las curvas tales que la normal en un punto cualquiera pasa por el on origen. (11) Demostrar que la curva para la cual la pendiente de la tangente en cualquier punto es proporcional a la abscisa del punto de contacto es una par´abola. abola. (12) (a) Hal Hallar lar las las soluc solucion iones es de: de: i. y + y = sen x ii. y + y = 3 cos( cos(22x) (b) Halle las soluciones de y + y = sen x + 3 cos(2 cos(2x x) cuya gr´afica afica pase por el origen (Piense, y no haga cuentas de m´as). as). (13) (13 ) Dada la ecuaci´ ecua ci´on on no homog´ hom og´enea enea y + a(x)y = b(x) donde a, b : R per´ıo do p > 0 y b 0. (a) Pruebe que una soluci´on on Φ de esta ecuaci´on on verifica:
→R son continuas con
≡
Φ(x Φ(x + p) p) = Φ(x Φ(x)
∀x ∈ R ⇔
Φ(0) = Φ( p) p)
(b) Encuentre las soluciones de per p er´´ıodo 2π para las ecuaciones: ecuaciones: y + 3y 3y = cos x,
y + cos(x cos(x)y = sen(2x sen(2x)
(14) Suponga que el ritmo al que se enfr´ enfr´ıa un cuerpo caliente es proporcional a la diferencia de temperatura entre el e´l y el ambiente que lo rodea (ley de enfriamiento de Newton). Un cuerpo se calienta calienta 110 ◦ C y se expone al aire libre a una temperatura de 10 ◦C. Al cabo de una hora su temperatura es de 60 ◦ C. ¿Cu´ anto tiempo adicional debe transcurrir anto ◦ para que se enfr´ enfr´ıe a 30 C? (15) Si la resistencia del aire que act´ ua sobre un cuerpo de masa m en ca´ ua ca´ıda libre ejerce una fuerza retardadora sobre el mismo proporcional a la velocidad (= kv) kv ) , la ecuaci´ on on diferencial del movimiento es: d2 y dy dv = g c , o bien = g cv dt2 dt dt donde c = k/m. k/m. Supongamos v = 0 en el instante t = 0, y c > 0. Encontrar limt→∞ v (t) (llamada velocidad terminal). Si la fuerza retardadora es proporcional al cuadrado de la velocidad, la ecuaci´on on se convierte en: dv = g cv 2 dt Si v (0) = 0, encuentre la velocidad terminal en este caso.
−
−
−
−
(16) La ecuaci´ on on y + P ( P (x)y = Q(x)yn , que se conoce como la ecuaci´on on de Bernoulli, es lineal cuando n = 0, 1. Demuestre que se puede reducir a una ecuaci´on on lineal para cualquier 1 n − valor de n = 1 por el cambio de variable z = y , y aplique apliqu e este est e m´etodo etod o para resolver las ecuaciones siguientes:
(a) xy + y = x4 y 3 (b) xy2 y + y 3 = x cos x (c) xy 3y = x4
−
CAP´ıTULO ıT ULO 2
Existencia y unicidad de soluci´ on on En este cap´ cap´ıtulo daremos resultados de existencia y unicidad local lo cal de soluci´ on on para sistemas de 1er. orden de la forma X = F ( F (t, X ). Necesitamos Necesitamos ciertas propiedades propiedades del campo F , F , a saber ´n Definicion o
∈
∈
2.1. Sea F ( F (t, X ) definida definida para t I y X Ω donde I es un intervalo de la recta n y Ω es un abierto de R . Decimos que F es Lipschitz en la variable X en I Ω si F es continua en las variables t y X en I Ω y existe una constante L tal que para todo t I , X, Y Ω,
×
×
∈
∈
F (t, X ) − F ( F (t, Y ) Y ) ≤ LX − Y . F (
×
Decimos que F es localmente localmente Lipschitz en la variable X en I Ω si para todo intervalo cerrado y acotado J contenido en I y todo conjunto cerrado y acotado Ω contenido en Ω se tiene que F es Lipschitz en J Ω .
×
´n Observacion o
2.1. La condici´ on on de Lipschitz local dice que hay una constante como la L de arriba para cada subconjunto J Ω como los descriptos, pero la constante puede ser distinta para distintas elecciones de los conjuntos. Adem´ as as es esencial esencia l aqu´ı que los conjuntos J y Ω sean acotados y est´en en contenidos, junto con sus bordes, b ordes, en I y Ω respectivamente .
×
× →
2.1. Sean I y Ω inter interv valo aloss de la recta. recta. Si f : I Ω R es continua y existe f x (t, x) y es continua en I Ω, se sigue que f es localmente Lipschitz en la variable x en I Ω. Ejemplo
×
×
En efecto, sea J un intervalo cerrado y acotado contenido en I y sea Ω un intervalo cerrado y aco acotad tadoo conten contenido ido en el interv intervalo alo Ω. Sea L = max f x (t, x) : t J , x Ω . Si t J , x, y Ω se tiene f ( f (t, x) f ( f (t, y) = f x (t, θ) (x y) L x y,
∈
{|
|
−
| |
|
∈
∈ }
∈
− |≤ | − |
ya que θ es un punto en el intervalo de extremos x e y y por lo tanto pertenece al intervalo Ω . ´n Observacion o
2.2. El ejemplo 2.1 se generaliza a
Rn
pero no haremos los detalles aqu´ aqu´ı.
2.2. Sea F ( F (t, X ) = A(t) X + X + b(t) con A(t) Rn×n y b(t) Rn . Si los coeficientes aij (t) , bi (t) de la matriz A(t) y el vector b(t) son funciones continuas de la variable t en un intervalo I se sigue que F es localmente Lipschitz en la variable X en I Rn . Si el intervalo I es cerrado y acotado, F es Lipschitz en la variable X en I Rn .
∈
Ejemplo
×
15
∈ ×
´ 2. EXISTENCIA Y UNICIDAD DE SOLUCI ON
16
En efecto, s´olo olo tenemos que ver esta ultima u ´ ltima afirmaci´ on. o n. Sea Sea K una constante mayor que aij (t) para t I y para todo i, j = 1, 1 , . . . , n. n. Entonces,
|
|
∈
n
A(t)X − A(t)Y 2 = A(t)(X )(X − Y ) Y )2 = n
≤ C n
|
aij (t) 2 (x j
|
i,j=1 i,j =1
n
i=1
aij (t)(x )(x j
j=1 j =1
2
− y j )
− y j )2 ≤ C nK 2n X − Y 2.
Por lo tanto A(t)X A(t)Y L X Y donde L2 = C n K 2 n y se sigue que F ( F (t, X ) es Lipschitz con constante L ya que F ( F (t, X ) F ( F (t, Y ) Y ) = A(t)X A(t)Y . Y .
−
≤ − −
−
Estamos en condiciones de enunciar el teorema de existencia de soluci´ on on para un sistema de 1er. orden. 2.1. Sean I un intervalo de la recta y Ω un abierto de Rn . Sea F Sea F ((t, X ) un campo localmente Lipschitz en la variable X en I Ω. Sean τ Sean τ I y ξ Ω. Si τ es interior a I , existen λ > 0 y una funci´ on continuamente diferenciable X : [τ λ, τ + λ] I Ω tales que Teorema
×
∈ ∈ − ⊂ → para todo t ∈ [τ − λ, τ + λ],
X (t) = F ( F (t, X (t)), )), X (τ ) τ ) = ξ.
Si τ es el extremo izquierdo de I , existen λ > 0 y una funci´ on continuamente diferenciable X : [τ, τ + λ] I Ω tales que
⊂ →
X (t) = F ( F (t, X (t)), )),
para todo t
X (τ ) τ ) = ξ.
∈ [τ, τ + λ],
Se tiene el resultado an´ alogo si τ es el extremo derecho del intervalo I . ´n Observacion o
2.3. Este Este teorem teoremaa no lo demost demostrar raremo emoss con esta genera generalid lidad. ad. Daremo Daremoss la demostraci´ on on de una versi´on on muy simplificada para el caso de una ecuaci´on. on. A saber, 2.2. Sea I Sea I un intervalo de la recta. Sea f Sea f ((t, x) una funci´ on Lipschitz en la variable x en I R. Sean τ Sean τ I , ξ R. Si τ es interior a I , existen λ > 0 y una funci´ on continuamente diferenciable x : [τ λ, τ + λ] I R tales que Teorema
×
∈ −
(2.1)
∈
⊂ →
x (t) = f ( f (t, x(t)), )),
para todo t
x(τ ) τ ) = ξ.
∈ [τ − λ, τ + λ],
Se tienen los resultados an´alogos alogos si s i τ es un extremo del intervalo I . ´ n. Demostracion. o
Supongamos que x(t) es una soluci´ soluc i´on on del problema. probl ema. Integrando Integrand o la ecuaci´ ecuacion o´n diferencial a partir de τ y usando la condici´ on on inicial tenemos t
(2.2)
x(t) = ξ +
f ( f (s, x(s)) ds,
para t en [τ
τ
−
− λ, τ + λ].
Rec´ıpro ıpro cam camente, ente, si x(t) es una funci´on on continua en [τ [τ λ, τ + τ + λ] y es soluci´on on de la ecuaci´on on integral (2.2) se sigue que x es continuamente diferenciable y que x = f ( f (t, x). Por Por otro lado, lado, evaluando en t = τ en la ecuaci´ on integral (2.2) vemos que x(τ ) on τ ) = ξ . Por Por lo tanto tanto (2.2) (2.2) es
´N 2. EXISTENCIA Y UNICIDAD DE SOLUCI O
17
equivalente equivalente a (2.1). El m´etodo etodo de construcci´ on on de una soluci´on on ser´a, a, por lo tanto, buscar una soluci´ on on de la ecuaci´ ecu aci´on on integral. integ ral. Y ´esto esto lo l o haremos hare mos por p or un m´etodo etod o iterativo. ite rativo. A saber, sab er, definirem d efiniremos os inductivamente x0 (t) = ξ, t
x1 (t) = ξ + e, inductivamente,
τ
f ( f (s, x0 (s)) ds,
t
(2.3)
xk (t) = ξ +
τ
f ( f (s, xk−1 (s)) ds.
Observemos que estas funciones est´an an definidas en I .
−
Si probamos que la sucesi´on on de funciones continuas xk (t) converge uniformemente en [τ [ τ λ, τ + λ] a una funci´on on x(t) se tendr´a, a, por un lado, que x(t) es continua en [τ [τ λ, τ + λ] y, por el otro, que
−
t
τ
t
f ( f (s, xk (s)) ds
→
f ( f (s, x(s)) ds.
τ
Por lo tanto x ser´ a una soluci´ on continua de (2.2) y en consecuencia una soluci´on on on de (2.1). Veamos entonces que la sucesi´on on de funciones as´ as´ı construida converge converge uniformemente en [τ λ, τ + λ] para alg´ un un λ > 0. Para Para eso, eso, veamos veamos que existe existe λ > 0 tal que la sucesi´on o n es uniformemente de Cauchy en [τ [τ λ, τ + λ]. Esto Esto signifi significa ca que satisface satisface que para todo ε > 0, existe k0 tal que para todo t [τ λ, τ + λ] y k, j k0 ,
−
− ∈ −
≥ |xk (t) − x j (t)| < ε.
Con este fin acotaremos acotaremos primero las diferencias diferencias de dos t´erminos erminos consecutiv consecutivos. os. Para Para esto necesitarem necesitaremos os utilizar utilizar la condici´ condici´ on de Lipschitz en la variable x de f en I R. Sea L la on constante de Lipschtiz, entonces como
×
t
xk+1 (t) = ξ + xk (t) = ξ + restando ambas ecuaciones vemos que
τ t τ
f ( f (s, xk (s)) ds, f ( f (s, xk−1 (s)) ds,
t
xk+1 (t)
≥ τ ,τ , |xk+1(t) − xk (t)| ≤
− xk (t) =
De aqu´ı que, para t
τ
[f ( f (s, xk (s))
f (s, xk−1 (s))] ds. − f (
t
| τ
t
f ( f (s, xk (s))
f (s, xk−1 (s))| ds ≤ L − f (
| τ
xk (s)
− xk−1(s)| ds.
≤ 21L tal que I λ := [τ − λ, τ + λ] ⊂ I . Entonces, max{|xk+1 (t) − xk (t)| , t ∈ [τ, τ + λ]} ≤ L |t − τ | max{|xk (t) − xk−1 (t)| , t ∈ [τ, τ + λ]} ≤ 12 max{|xk (t) − xk−1(t)| , t ∈ [τ, τ + λ]}.
Sea ahora λ
´ 2. EXISTENCIA Y UNICIDAD DE SOLUCI ON
18
− λ, τ ].]. Llamemos ahora mk = max{|xk (t) − xk−1 (t)| , t ∈ I λ }. Tenemos entonces, 1 mk+1 ≤ mk . 2
A una desigualdad an´aloga aloga se llega en el intervalo [τ [τ
Iterando esta desigualdad, obtenemos mk
≤ 2k1−1 m1.
Finalmente, si j = k + m, m−1
|xk (t) − x j (t)| = |xk (t) − xk+m(t)| = |
i=0 m−1
m−1
≤ i=0
Por lo tanto, dado ε > 0 si j, k
(xk+i (t)
mk+i+1
≤ m1
1
i=0
2k+i
m1 = k 2
− xk+i+1(t))| m−1
i=0
1 2i
≤ 2mk−11 .
1 < ε se tiene para t ∈ [τ − λ, τ + λ], ≥ k0 donde 2m k −1 0
|x j (t) − xk (t)| < ε. Claramente, de la demostraci´ on on se ve que si τ es el extremo izquiedo de I y λ que [τ, [τ, τ + λ] I , se tiene una soluci´on on en el intervalo intervalo [τ, τ + λ].
⊂
An´ alogamente, alogamente, si τ es el extremo derecho de I y λ una soluci´ soluci ´on on en el intervalo [τ λ, τ ]. τ ].
−
≤ 1/2L es tal
≤ 1/2L es tal que [τ [τ − λ, τ ] τ ] ⊂ I , se tiene
Antes de discutir cu´al al es el mayor intervalo que contiene a τ donde hay definida una soluci´ on on del problema, veamos un resultado de continuidad de las soluciones respecto del dato inicial que implicar´ a la unicidad de soluci´ on. on. Teorema
2.3. Sean I y f como en el Teorema 2.2. Sean τ
Sean x1 , x2 : [τ, η]
⊂ I → R soluciones de
(2.4)
x = f ( f (t, x)
∈ I y ξ1, ξ2 ∈ R.
en [τ, η],
con xi (τ ) τ ) = ξi , i = 1, 2. Existe una constante C (η ) dependiente de η tal que
|x1(t) − x2(t)| ≤ C (η) |ξ1 − ξ2| para t ∈ [τ, η]. Se tiene el mismo resultado si x1 , x2 : [η, τ ] ⊂ I → R son soluciones de
(2.5)
x = f ( f (t, x)
con xi (τ ) τ ) = ξi , i = 1, 2.
en [η, τ ] τ ],
´N 2. EXISTENCIA Y UNICIDAD DE SOLUCI O ´ n. Demostracion. o
19
Lo demostraremos en el caso que η > τ . τ . La demostraci´ on on del otro caso es
enteramente an´aloga. aloga. Sea L la constante de Lipschitz de f en I t, tenemos
on (2.4) para x1 de τ a × R. Integrando la ecuaci´on t
x1 (t) = ξ1 +
f ( f (s, x1 (s)) ds.
τ
An´ alogamente, integrando (2.4) para x2 obtenemos alogamente, t
x2 (t) = ξ2 +
f ( f (s, x2 (s)) ds.
τ
Restando ambas ecuaciones
t
x1 (t)
− x2(t) = ξ1 − ξ2 +
y por lo tanto,
τ
[f ( f (s, x1 (s))
− f ( f (s, x2 (s))] ds
t
|x1(t) − x2(t)| ≤ |ξ1 − ξ2| + L
(2.6)
|
|
Para obtener de (2.6) una desigualdad para x1 (t) Lema
x1 (s)
τ
− x2(s)| ds.
− x2(t)| usaremos usaremos el
≥ 0 continua en un intervalo que contiene a τ a τ .. Si t g (t) ≤ A + B g (s) ds ,
2.1 (de Gronwall). Sea g
(2.7)
τ
se sigue que g (t)
≤ AeB|t−τ |.
|
Suponiendo probado el Lema de Gronwall, podemos aplicarlo con g (t) = x1 (t) obtenemos
|x1(t) − x2(t)| ≤ |ξ1 − ξ2|eL(t−τ )τ ) ≤ C (η) |ξ1 − ξ2|
si t
− x2(t)| y
∈ [τ, η].
τ ) . donde C (η ) = eL(η−τ )
Probemos ahora el Lema de Gronwall. ´ n del Lema de Gronw Demostracion o Gronwall.
Lo haremos en el caso t
en el otro caso es totalmente an´ aloga. aloga.
≥ τ .τ . La demostraci´on on
t
Sea G(t) =
g (s) ds. ds. Entonces (2.7) nos dice que
τ
O, lo que es lo mismo,
G (t)
≤ A + B G(t).
G (t)
− B G(t) ≤ A.
τ ) . Tenemos Multipliquemos la inecuaci´ on on por e−B (t−τ )
τ ) τ ) e−B (t−τ ) G(t) = e−B (t−τ ) G (t)
− B G(t) ≤ A e−B(t−τ )τ ).
´ 2. EXISTENCIA Y UNICIDAD DE SOLUCI ON
20
Integrando de τ a t y usando que G(τ ) τ ) = 0, tenemos e
−B (t−τ ) τ )
t
G(t)
≤A
τ ) e−B (s−τ ) ds =
τ
− BA
De aqu aq u´ı que G(t)
≤ BA
Como la desigualdad (2.7) dice que g (t) se sigue que g (t)
≤ A + B BA
que es lo l o que q ue quer qu er´´ıamos demostrar. demostr ar.
τ ) eB (t−τ )
τ ) e−B (t−τ )
−1 ,
−1 .
≤ A + B G(t),
τ ) eB (t−τ )
−1
τ ) = AeB (t−τ ) .
Como corolario del Teorema 2.3 se obtiene el resultado de unicidad. A saber, 2.4. Sean I y f como como en el Teor eorema 2.2. Sean Sean τ τ I y xi : J i R para i = 1, 2 tales que
Teorema
τ
∈ J 2 ⊂
→
xi = f ( f (t, xi )
∈ I y ξ ∈ R. Sean τ ∈ J 1 ⊂ I ,
en J i ,
xi (τ ) τ ) = ξ.
∈ J 1 ∩ J 2. ´ n. Sea τ ∈ [t1 , t2 ] ⊂ J 1 ∩ J 2 . Demostracion. o
Entonces, x1 (t) = x2 (t) si t
Probar Probaremo emoss que x1 = x2 en [t1 , t2 ]. Como Como el el intervalo es arbitrario, se sigue que x1 = x2 en J 1 J 2 .
∩
Probamos primero que x1 = x2 en [τ, t2 ]. (Podr´ıamos ıamos tener τ = t1 y no habr´ıa ıa que probar nada m´ as). En efecto, por el Teorema 2.3 tenemos as).
|x1(t) − x2(t)| ≤ C (t2) |ξ − ξ| = 0
en [τ, t2 ].
|x1(t) − x2(t)| ≤ C (t1) |ξ − ξ| = 0
en [t1 , τ ] τ ].
An´ alogamente, alogamente, El teorema est´a demostrado. demostrado.
´n Observacion o
2.4. Estos resultados de continuidad respecto de los datos iniciales y de unicidad son v´ alidos alidos bajo las condiciones condiciones del Teorema Teorema 2.1. Pero Pero no daremos daremos sus demostraciones demostraciones aqu´ aqu´ı aunque son enteramente enteramente similares simi lares a las demostraciones de los Teoremas 2.3 y 2.4. ´n Observacion o
2.5 (Prolongaci´ (Prolongaci´ on on de soluciones). A partir del Teorema 2.4 vemos que si tenemos una soluci´on on x1 del problema (2.8)
x = f ( f (t, x), x(τ ) τ ) = ξ,
´N 2. EXISTENCIA Y UNICIDAD DE SOLUCI O
∈
⊂
en un intervalo τ J 1 I y una soluci´on on x2 del problema (2.8) en τ realidad una soluci´ on on x¯ de (2.8) en J 1 J 2 . En efecto, si definimos
∪
x ¯(t) = se tiene que x ¯ est´ a bien definida en J 1
x1 (t) x2 (t)
si t si t
21
∈ J 2 ⊂ I , tenemos en
∈ J 1, ∈ J 2,
∪ J 2 y es soluci´on on de (2.8) (2. 8) ah´ı. ı.
Podemos por lo tanto definir la soluci´ on maximal de maximal de (2.8). A saber, sea τ intervalo donde hay definida una soluci´on, on, es decir
∈ J ⊂ I el mayor
J = ∪{J : τ ∈ J y J es un intervalo en el que hay definida una soluci´on on de (2.8)}. Entonces hay una soluci´on on definida en J , esta soluci´ on o n es unica u ´ nica y no es posible prolongarla a un intervalo m´ as as grande que J . Esta es la llamada soluci´ on maximal . ´n Observacion o
2.6. Supongamos que f ( f (t, x) es Lipschitz en la variable x en [t1 , t2 ] R para todo intervalo cerrado y acotado [t [t1 , t2 ] contenido en I . Veam eamos os que la soluci soluci´ on o´n maximal est´ a definida en todo I .
×
∈
⊂
En efecto, sea τ [t1 , t2 ] I . Aplicando el Teorema 2.2 con I reemplazado por el intervalo [τ, t2 ], la construcci´ on on nos da una soluci´on on en todo [τ, [τ, t2 ] si t2 τ 1/2L. Si no, obtenemos una soluci´ on on x1 en [τ, τ + 1/ 1/2L]. Consideremos ahora, para τ 1 = τ + 1/ 1/2L el problema
− ≤
x = f ( f (t, x)
en [τ 1 , τ 1 + λ],
x(τ 1 ) = x1 (τ 1 ).
− ≤
Por el Teorema 2.2, si t2 τ 1 1/2L, existe una soluci´on on x2 en [τ 1 , t2 ]. Esta Esta soluci soluci´ on o´n se pega bien con x1 en τ 1 y por lo tanto obtenemos una soluci´on on en [τ, [τ, t2 ].
−
Si por el contrario, t2 τ 1 > 1/2L, existe una soluci´on on x2 en [τ 1 , τ 1 + 1/ 1/2L], y por lo tanto 1 tenemos una soluci´on o n en [τ, [τ, τ + 2 2L ]. Siguiendo Sigui endo as´ as´ı, como existe un k N tal que τ + k 21L > t2 , vemos que en un n´ umero umero finito de pasos debemos tener una soluci´ on on definida en todo [τ, [τ, t2 ].
∈
An´ alogamente alogamente se ve que hay una soluci´on on definida en [t [t1 , τ ]. τ ]. Por lo tanto, hay una soluci´ on o n definida en [t [t1 , t2 ]. De dond donde, e, la soluc soluci´ i´ on on maximal est´ a definida ah´ ah´ı. Como el intervalo intervalo τ [t1 , t2 ] I es arbitrario, se sigue que la soluci´on on maximal est´ a definida en todo I .
∈
⊂
´n Observacion o
2.7. Cuando la soluci´ on on maximal est´ a definida en todo I , decimos que la soluci´ on on es global . Por la Observ Observaci´ aci´ on 2.6, vemos que bajo las condiciones del Teorema 2.2, la on soluci´ on on maximal maximal es global. Este resultado resultado tambi´ tambi´en en es cierto cierto si, en el Teorema eorema 2.2, en lugar de una funci´on on f ( f (t, x) tenemos un campo F ( F (t, X ) Lipschitz en la variable X en J Rn para todo intervalo cerrado y acotado J I . Es decir, decir, el resulta resultado do de existe existenci nciaa glo global bal es v´ alido alido para sistemas de n ecuaciones con n inc´ ognitas ognitas de la forma
×
⊂
X = F ( F (t, X ) si F es Lipschitz en la variable X en J
× Rn para todo intervalo cerrado y acotado J ⊂ I .
´ 2. EXISTENCIA Y UNICIDAD DE SOLUCI ON
22
En particular, las soluciones de un sistema lineal con coeficientes continuos en un intervalo abierto I son globales, tanto en el caso homog´eneo eneo como no homog´eneo. eneo. Enunciaremos este resultado para referencia posterior.
⊂
2.5. Sea I Sean aij (t), bi (t) R un intervalo abierto (posiblemente todo R). Sean n funciones continuas en I para i, j = 1, , n. Sean Sean τ I , ξ R . Existe Existe entonc entonces es una unica ´ soluci´ on X = (x1 , , xn ) de Teorema
···
(2.9)
···
∈
∈
x1 = a11 x1 + a12 x2 +
· · · + a1nxn + b1, x2 = a21 x1 + a22 x2 + · · · + a2n xn + b2 ,
.. .
xn = an1 x1 + an2 x2 +
· · · + annxn + bn,
que satisface X (τ ) τ ) = ξ . Esta soluci´ on est´ a definida en todo el intervalo I . ´n Observacion o
2.8. La condici´ on de Lipschitcianidad global no puede relajarse a Lipschiton cianidad local. En efecto, como vimos en la introducci´ on, on, la soluci´ on on del problema x = x2 x(0) = 1 es x(t) =
1
−∞
. Por lo tanto, el intervalo maximal de existencia para este problema es ( , 1), 1 t mientras que I = R; ya que, en este ejemplo, f ( f (t, x) = x2 es localmente Lipschitz en la variable x en R R (pero no lo es globalmente).
−
×
Ejercicios
(1) Sea F : R
on localmente Lipschitz. Lipschitz. Probar Probar que si x(t) es una soluci´on on de → R una funci´on x (t) = F ( F (x(t))
que verifica x(t0 + T ) T ) = x(t0 ) para alg´ un un t0 de per pe r´ıodo ıo do T . T .
∈ R entonces x(t) es una funci´on on peri´odica odica
(2) Sea F : R funci´on on localmente Lipschitz. Verificar que si x1 (t) y x2 (t) son dos R una funci´ soluciones de x (t) = F ( F (x(t)) tales que x1 (t1 ) = x2 (t2 ) para ciertos t1 , t2 R , entonces Im( Im (x1 ) = I m(x2 ).
→
∈
(3) Sea F : R on localmente Lipschitz tal que F ( on F ( p1 ) = F ( F ( p2 ) = 0 con R una funci´ p1 < p2 . Probar que si x0 ( p1 , p2), entonces la soluci´ on on de
→
∈
x (t) = F ( F (x(t)), )), est´ a definida para todo t (4) Sea F : R
x(0) = x0
∈ R.
→ R una funci´ f unci´on on lo calmente Lipschitz y sea x(t) una soluci´ on on de x (t) = F ( F (x(t))
tal que x(t)
→ x0 cuando t → +∞ . Probar que F ( F (x0 ) = 0.
CAP´ıTULO ıT ULO 3
Sistema Sistemass lineal lineales es de 1er. 1er. orden orden y ecuaci ecuacione oness lineal lineales es de orden orden
n
En este cap´ cap´ıtulo estudiaremos el conjunto de soluciones de sistemas lineales li neales de n ecuaciones con n inc´ ognitas. ognitas. Probaremos, en el caso homog´eneo, eneo, que el conjunto de soluciones es un espacio vectorial lo que nos dice c´omo omo ser´ a la soluci´ on on general general.. Dada Dada la relaci relaci´ on o´n entre ecuaciones de orden n y sistemas de n ecuaciones, deduciremos resultados an´alogos alogos para el caso de ecuaciones. Finalmente, Finalm ente, daremos da remos un m´etodo etod o para encontrar e ncontrar la l a soluci´ soluc i´ on general de sistemas (resp. ecuaciones) on no homog´ ho mog´eneas eneas a partir par tir de d e la soluci´ solucion o´n general del sistema (resp. ecuaci´on) on) homog´eneo eneo asociado. asoc iado. 1. Generalidades Generalid ades y sistemas si stemas homog´ eneos eneos
Empecemos probando que el conjunto de soluciones de un sistema lineal homog´eneo eneo forma un espacio vectorial de dimensi´ on on n.
⊂
3.1. Sea I intervalo lo abierto abierto.. Sean Sean aij (t) funciones funciones continua continuass en I . R un interva × n n Sea A(t) = (aij (t)) R . El conjunto conjunto de las soluciones soluciones del sistema sistema lineal lineal homog´ homog´ eneo eneo de n ecuaciones con n inc´ in c´ogni og nita tass Teorema
∈
X = A(t)X
(3.1) es un espacio vectorial de dimensi´ on n. ´ n. Demostracion. o
Recordemos que todas las soluciones est´an an definidas en todo el intervalo I . Por lo tanto, podemos sumarlas y multiplicarlas multiplicarlas por escalares. Lo que tenemos que ver, para ver que es un espacio vectorial, es que al sumar dos soluciones obtenemos otra soluci´on on y lo mismo sucede sucede si multiplica multiplicamos mos una soluci´ on por un escalar. Esto es consecuencia de la linealidad de la on operaci´on on de derivaci´on on y de la funci´on on A(t)X (en la variable X ). ). En efecto, sean X 1 y X 2 dos n soluciones y X = X 1 + X 2 es decir, X es la funci fu nci´´on on de I en R definida por X (t) = X 1 (t) + X 2 (t) para cada t I . Veamos que X es tambi´ tamb i´en en soluci´ sol uci´on on de (3.1). Tenemos,
∈
X (t) = X 1 (t) + X 2 (t) = A(t)X 1 (t) + A(t)X 2 (t) = A(t)(X )(X 1 (t) + X 2 (t)) = A(t)X (t). Por lo tanto X satisface (3.1). Sea ahora c
∈ R y sea X = c X 1, entonces
X (t) = cX 1 (t) = cA( cA(t)X 1 (t) = A(t)(cX )(cX 1 (t)) = A(t)X (t).
y nuevamente obtenemos que X es soluci´ on on de (3.1). Veamos ahora que hay una base del espacio de soluciones formada por exactamente n soluciones. Para esto, aplicamos el Teorema 2.5 con τ I cualquiera fijo. Sean X i , i = 1, . . . , n, n, las soluciones maximales de (3.1) que verifican X i (τ ) τ ) = ei , donde e1 , . . . , en es la base can´onica onica n de R .
∈
23
{
}
24
3. SIST SISTEM EMAS AS LINE LINEAL ALES ES DE 1ER. 1ER. OR ORDE DEN N Y ECUA ECUACI CION ONES ES LINE LINEAL ALES ES DE OR ORDE DEN N n
Obtene Obt enemo moss as´ı n soluciones. soluciones. Veamos que son linealmen linealmente te independien independientes, tes, ´esto esto es, que la u unica ´ nica manera de obtener la funci´on on 0 al hacer una combinaci´ on lineal de estas soluciones es on tomando todos los coeficientes iguales a 0. Supongamos entonces que tenemos constantes c1 , . . . , cn tales que
· · · + cnX n(t) = 0 para ara tod todoo t ∈ I. Debemos probar que necesariamente c1 = c2 = · · · = cn = 0. En efecto, tomando t = τ en (3.2) obtenemos c1 e1 + · · · + cn en = 0. Como {e1 , . . . , en } son linealmente independientes, se sigue que necesariamente c1 = c2 = · · · = (3.2)
c1 X 1 (t) +
cn = 0 como quer´ quer´ıamos demostrar. demost rar.
{
}
Resta ver que X 1 , . . . , Xn genera generan n el espacio espacio de soluci solucione oness de (3.1). Es decir, decir, que toda soluci´ on puede escribirse como combinaci´ on on on lineal de estas n soluciones. soluciones. En efecto, sea X una soluci´ on o n de (3.1) y sea ξ = X (τ ). τ ). Como Como e1 , . . . , en es una base de Rn , existen constantes c1 , . . . , cn tales que ξ = c1 e1 + + cn en .
{
}
···
··· ···
∈
Construyamos ahora la siguiente funci´ on: on: Y ( Y (t) = c1 X 1 (t) + + cn X n (t) para t I . Entonces, como el conjunto de soluciones de (3.1) es un espacio vectorial, se sigue que Y es tambi´ ta mbi´en en una soluci´ on on de (3.1). Pero, por la elecci´on on de las constantes c1 , , cn , se tiene que Y ( Y (τ ) τ ) = ξ . De modo que tenemos dos soluciones de (3.1) con el mismo dato inicial ξ en t = τ . τ . Por el teorema de unicidad de soluci´on on se sigue que X = Y . Y . Recordando Recordan do qui´en en es Y vemos que X (t) = c1 X 1 (t) +
· · · + cnX n(t)
par para todo todo t
∈ I,
como quer´ quer´ıamos demostrar. demostr ar.
Un resultado importante, del cu´ al esencialmente ya probamos una parte en el teorema anal terior es el siguiente ´n Proposicion o
3.1. Sean {X 1 , . . . , Xn } soluciones de (3.1) y sea τ ∈ I cualquiera. Entonces {X 1, . . . , Xn } son linealmente independientes como funciones de t en I si y s´ olo si los vectores n {X 1(τ ) τ ), . . . , Xn (τ ) τ )} son linealmente independientes en R . ´ n. Demostracion. o
Como en la demostraci´ on del Teorema 3.1 supongamos que los vectores on X 1 (τ ) τ ), . . . , Xn (τ ) τ ) son linealmente independientes en Rn .
{
}
{
}
Veamos que las funciones X 1 , . . . , Xn son linealmente independientes. En efecto, si c1 , . . . , cn son constantes tales que la funci´on on c1 X 1 (t) + + cn X n (t) es la funci´on on 0, veamos que c1 = c2 = = cn = 0. En efecto, evaluando en t = τ vemos que
···
c1 X 1 (τ ) τ ) +
{
···
· · · + cnX n(τ ) τ ) = 0
}
y como X 1 (τ ) τ ), . . . , Xn (τ ) τ ) son linealmente independientes en cn = 0.
{
Rn ,
}
se sigue que c1 = c2 =
··· =
Rec´ Rec´ıprocamente, supongamos que las soluciones X 1 , . . . , Xn son linealmente independientes y veamos que los vectores X 1 (τ ) τ ), . . . , Xn (τ ) τ ) tambi´ tambi´ en en lo son. En efecto, efecto, supongamos supongamos que c1 X 1 (τ ) τ ) + + cn X n (τ ) τ ) = 0.
{
}
···
´ NEOS 1. GENERALIDADES Y SISTEMAS HOMOG E
···
25
∈
Sea Y ( Y (t) = c1 X 1 (t) + + cn X n (t) para t I . Entonces, Entonces, Y es una soluci´ on on de (3.1) con dato inicial Y ( Y (τ ) τ ) = 0. Por el teorema de unicidad de soluci´on on se sigue que Y = 0. Esto es, c1 X 1 (t) +
{
· · · + cnX n(t) = 0
para ara tod todoo t
∈ I.
}
Como las funciones X 1 , . . . , Xn son linealmente independientes obtenemos que necesariamente c1 = c2 = = cn = 0 como quer´ quer´ıamos demostrar. demost rar.
···
Tenemos inmediatamente el siguiente corolario
{
}
∈
3.1. Sean X 1 , . . . , Xn soluciones de (3.1) y sea τ, η I cualesquier cualesquiera. a. Entonces los vectores X 1 (τ ) τ ), . . . , Xn (τ ) τ ) son linealmente independientes si y s´ olo si los vectores X 1 (η ), . . . , Xn (η ) lo son. Corolario
{ }
{
´n Observacion o
}
{
}
3.1. Sea X 1 , . . . , Xn una base de soluciones soluciones de (3.1). Construyam Construyamos os una matriz ubicando a estos vectores como columnas y llam´emosla emosla Q(t). A una matriz de este tipo la llamamos matriz fundamental . Es f´acil acil ver que Q (t) = A(t)Q(t) para para tod todo t
∈ I.
ya que las columnas de A(t)Q(t) son los vectores A(t)X j (t). Como el determinan determinante te de una matriz es no nulo si y s´ olo olo si sus columnas son linealment linealmentee independientes, el Corolario 3.1 dice que
det Q(τ ) τ ) = 0 para un τ
∈ I si y s´olo 0 olo si det Q(t) =
para para todo todo t
∈ I.
Observemos adem´as as que como toda soluci´ on de (3.1) es combinaci´ on on lineal de las columnas on de Q, se sigue que toda soluci´on on es de la forma X (t) = Q(t)C para alg´ un un vector C =
c1 .. .
.
cn
Observemos por otro lado, que dada una matriz U ( U (t) Rn×n cualquiera (es decir, si no pedimos que las columnas sean soluci´on on de un mismo sistema lineal homog´eneo), eneo), no tiene por p or qu´ e ser cierto que el determinan determinante te es distinto distinto de cero en un punto punto si y s´ olo olo si lo es en todos los puntos. Por ejemplo, si t 0 U ( U (t) = , 0 t
∈
se tiene que det U (0) (0) = 0 y det U (1) U (1) = 1,
A partir de los resultados anteriores sobre el conjunto de soluciones de sistemas lineales, y dada la equivalencia de una ecuaci´on on de orden n con un sistema de primer orden de n ecuaciones con n inc´ognitas, ognitas, obtenemos resultados sobre el conjunto de soluciones de una ecuaci´on on lineal de orden n. En efecto, Consideremos Consid eremos la ecuaci´on on (3.3)
x(n) + an−1 (t)x(n−1) + an−2 (t)x(n−2) +
· · · + a1(t)x + a0(t)x = f ( f (t).
26
3. SIST SISTEM EMAS AS LINE LINEAL ALES ES DE 1ER. 1ER. OR ORDE DEN N Y ECUA ECUACI CION ONES ES LINE LINEAL ALES ES DE OR ORDE DEN N n
Como vimos en la introducci´ on, on, si x es soluci´ on on de (3.3) se sigue que X = (x, x , x , . . . , x(n−1) ) es soluci´ on on del sistema
(3.4)
x0 = x1 , x1 = x2 , .. . xn −2 = xn−1 , xn −1 =
−a0(t)x0 − a1(t)x1 − · · · − an−1(t)xn−1 + f ( f (t). Rec´ıpro ıpro cam camente, ente, sea X = (x0 , x1 , · · · , xn−1 ) una soluci´ on de (3.4), entonces la funci´on on on x(t) =
x0 (t) es soluci´ on on de (3.3). En efecto, efecto, la primer ecuaci´ ecuaci´on on del sistema sistema (3.4) dice que x1 = x , la segunda dice que x2 = x1 = x . La tercer tercera a dice que x3 = x2 = x . Y as as´´ı hasta la pen´ p en´ ultima ultima − ( n 1) que dice que xn−1 = xn−2 = x . Finalmente Finalmente,, con esta informaci´ on o n la ultima u ´ ltima ecuaci´ ecuaci´ on on dice que x(n) = xn −1 = a0 (t)x0 a1 (t)x1 an−1 (t)xn−1 + f ( f (t)
−
=
−
−···−
−an−1(t)x(n−1) − an−2(t)x(n−2) − · · · − a1(t)x − a0(t)x + f ( f (t).
Es decir, x es soluci´ on on de (3.3). Se tiene, por lo tanto, el siguiente resultado Teorema
3.2.
∈
−
(1) Sea I un intervalo abierto de la recta y sea τ I . Sean Sean ai (t), i = 1, . . . , n 1 y f ( f (t) funciones continuas en I en I . Para cada n-upla (y0 , y1 , . . . , yn−1 ) existe una unica ´ soluci´ on de (3.3) que satisface x(τ ) τ ) = y0 , x (τ ) τ ) = y1 , x (τ ) τ ) = y2 ,
· · · x(n−1)(τ ) τ ) = yn−1 .
Adem´ as la soluci´ on est´ a definida en todo el intervalo I . (2) Sea I Sea I un intervalo abierto de la recta y sean ai (t), i = 1, . . . , n 1 funciones continuas en I . El conjun conjunto to de soluci solucione oness de (3.3) cuando cuando f = 0 – es decir en el caso de la ecuaci´ on lineal homog´enea enea de orden n – es un espacio vectorial de dimensi´ on n.
−
{
}
(3) Bajo las mismas hip´ otesis que en (2), un conjunto x1 , . . . , xn de soluciones de (3.3) es linealmente independiente – y por ende una base de soluciones – si y s´ olo si
W ( W (x1 , x2 , . . . , xn )(τ )(τ )) := det
x1 (τ ) τ ) x1 (τ ) τ ) x1 (τ ) τ ) .. .
x2 (τ ) τ ) x2 (τ ) τ ) x2 (τ ) τ ) .. .
··· ··· ···
(n−1) x1 (τ ) τ )
(n−1) x2 (τ ) τ )
···
..
.
xn (τ ) τ ) xn (τ ) τ ) xn (τ ) τ ) .. . (n−1)
xn
=0
(τ ) τ )
W ( W (x1 , x2 , . . . , xn ) se llama el Wronskiano de las soluciones x1 , . . . , xn y, dado un con junto de soluciones de la ecuaci´ on lineal (3.3) en el caso homog´ h omog´ eneo eneo f = 0, se tiene que W ( W (x1 , x2 , . . . , xn )(τ )(τ )) = 0 para alg´ un τ I si y s´ olo si W ( W (x1 , x2 , . . . , xn )(t )(t) = 0 para todo t I .
∈
∈
´ NEOS 1. GENERALIDADES Y SISTEMAS HOMOG E ´ n. Demostracion. o
27
Probemos Probemos (1). Veamos primero primero la existencia de soluci´ on on para cada n-upla
de datos iniciales. Sea ξ = (y0 , y1 , . . . , yn−1 ) y sea X = (x0 , x1 , . . . , xn−1 ) la soluci´ on o n de (3.4) con dato inicial X (τ ) τ ) = ξ . Como Como vimo vimoss ante antes, s, x = x0 es soluci´ on o n de (3.3) y adem´as as x = x1 , , x = tanto, x(τ ) por lo lo x2 , , . . . , x(n−1) = xn−1 . Por lo tanto, τ ) = y0 , x (τ ) τ ) = y1 , . . . , x(n−1) (τ ) τ ) = yn−1 . Y por tanto tanto existe soluci´ on, on, como quer´ quer´ıamos demostrar. demost rar. Probemos ahora la unicidad de soluci´ on. on. En efecto, si x y x ˜ son dos soluciones de (3.3) con los mismos datos iniciales en t = τ , τ , veamos que son iguales. ˜ = (˜ ˜ son soluciones Sean X = (x, ( x, x , x , . . . , x(n−1) ) y X x, x, x ˜ , x ˜ , . . . , x˜(n−1) ). Entonces X y X ˜ (τ ). ˜ (t) para de (3.4) y X (τ ) τ ) = X τ ). Por el teorema de unicidad de soluci´ on, on, se sigue que X (t) = X todo t I . En particular, x(t) = x ˜(t) para todo t I , como afirmamos.
∈
∈
Probemos (2). Recordemos Recordemos que en este caso f = 0. Por Por lo tanto tanto se trata de una ecuaci ecuaci´ on o´n homog´enea enea y el sistema sistem a lineal equivalente tambi´en en es homog´eneo. eneo. Es f´ acil ver que el conjunto de soluciones es un espacio vectorial. Para ver que el espacio vectorial de soluciones tiene dimensi´ on on n, basta demostrar que hay 1 2 n n soluciones x , x , . . . , x linealmente independientes tales que cualquier soluci´on on se expresa 1 2 n en la forma x = c1 x + c2 x + + cn x para alguna elecci´ elecci´ on on de constantes c1 , . . . , cn . Veamos enton ent onces ces que qu e ´esto es to es as´ı. ı.
···
Fijemos τ I . Para Para cada cada i = 1, . . . , n, n, sea X i = (xi0 , xi1 , . . . , xni −1 ) la soluci´ on on de (3.4) con n X i (τ ) τ ) = ei donde e1 , . . . , en es la base can´onica onica de R . Sabem Sabemos os que que X 1 , . . . , Xn es una base del conjunto de soluciones de (3.4).
∈
{
}
{
}
Sea ahora x una soluci´ on on de (3.3) y llamemos X = (x, ( x, x , x , . . . , x(n−1) ). Como X es soluci´on on de (3.4), existen constantes c1 , . . . , cn tales que X = c1 X 1 + + cn X n , en particular x = c1 x10 + c2 x20 + con x10 , . . . , xn0 soluciones de (3.3).
···
· · · + cnxn0 .
Veamos que las soluciones x10 , . . . , xn0 son linealmente independientes. En efecto, si tuvi´eeramos x = c1 x10 + c2 x20 + + cn xn0 = 0, tend te ndr´ r´ıamo ıa moss x = c1 (x10 ) + c2 (x20 ) + + cn (xn0 ) = 0,
{
}
··· ··· x = c1 (x10 ) + c2 (x20 ) + · · · + cn (xn0 ) = 0, .. .
x(n−1) = c1 (x10 )(n−1) + c2 (x20 )(n−1) +
· · · + cn(xn0 )(n−1) = 0.
Como el sistema (3.4) dice que (x (xi0 )(k) = xki , se sigue que X := (x, x , x , . . . , x(n−1) ) = c1 X 1 + c2 X 2 +
{
}
y como X 1 , . . . , Xn son linealmente independientes, c1 = c2 =
Con lo cual hemos demostrado el punto (2).
· · · = cn = 0.
· · · + cnX n = 0
28
3. SIST SISTEM EMAS AS LINE LINEAL ALES ES DE 1ER. 1ER. OR ORDE DEN N Y ECUA ECUACI CION ONES ES LINE LINEAL ALES ES DE OR ORDE DEN N n
{
}
Finalmente Finalmente,, probemos el punto (3). Sean x1 , x2 , . . . , xn soluciones soluciones de (3.3). Veamos que son linealmente independientes si y s´olo olo si X 1 , X 2 , . . . , Xn lo son, donde
{
}
(n−1)
X i = (xi , xi , xi , . . . , xi
(3.5)
)
{
es soluci´ on on del sistem sistemaa equiv equivale alent ntee (3.4). (3.4). En efecto, efecto, supongam supongamos os que X 1 , X 2 , . . . , Xn linealmente independientes, y que se tiene para ciertas constantes c1 x1 + c2 x2 +
···
} son
· · · + cnxn = 0
y veamos que c1 = c2 = = cn = 0. En efecto, efecto, como en la demost demostrac raci´ i´ on on del punto (2), si llamamos x = c1 x1 + c2 x2 + + cn xn vemos, derivando sucesivamente, que
···
x = c1 x1 + c2 x2 + .. . (n−1)
x(n−1) = c1 x1
· · · + cnxn = 0,0 , (n−1)
+ c2 x2
+
Por lo tanto, c1 X 1 + c2 X 2 + + cn X n = 0. 0. independien independientes, tes, se sigue que c1 = c2 = = cn = 0.
···
{
···
· · · + cnxn(n−1) = 0.0 . Como {X 1 , X 2 , . . . , Xn } son linealment linealmentee
}
Supongamos ahora que x1 , x2 , . . . , xn son linealmente independientes, y que se tiene para ciertas constantes (3.6)
X = c1 X 1 + c2 X 2 +
y veamos que c1 = c2 = del vector X , que
· · · = cn = 0.
En efecto, efecto, se sigue de (3.6), (3.6), mirand mirandoo la primer primer entra entrada da
c1 x1 + c2 x2 +
{
· · · + cnX n = 0
· · · + cnxn = 0.
}
Como x1 , x2 , . . . , xn son linealmente independientes, se sigue que c1 = c2 =
· · · = cn = 0.
Ahora, recordemos que un conjunto de n soluciones de un sistema lineal homog´eneo eneo de n ecuaciones con n inc´ognitas ognitas es linealmente independiente si y s´olo olo si, para la matriz Q(t) cuyas columnas son las n soluciones del sistema, se tiene det Q(τ ) τ ) = 0
para ara alg´ alg´ un un τ
∈ I.
y que
∈ I si y s´olo 0 para todo t ∈ I. olo si det Q(t) = Como en nuestro caso, las soluciones {X 1 , . . . , Xn } del sistema (3.4) vienen dadas por (3.5), det Q(τ ) τ ) = 0 para alg´ un un τ
la matriz Q(τ ) τ ) resulta
Q(τ ) τ ) =
x1 (τ ) τ ) x1 (τ ) τ ) x1 (τ ) τ ) .. . (n−1)
x1
x2 (τ ) τ ) x2 (τ ) τ ) x2 (τ ) τ ) .. . (n−1)
(τ ) τ ) x2
(τ ) τ )
··· ··· ··· ..
.
···
xn (τ ) τ ) xn (τ ) τ ) xn (τ ) τ ) .. . (n−1)
xn
(τ ) τ )
.
Por lo que det(Q det(Q(τ )) τ )) = W ( W (x1 , x2 , . . . , xn )(τ )(τ )) y se tiene lo enunciado en el punto (3).
´ NEOS 2. SISTEMAS NO HOMOG E
29
2. Sistemas no homog´ eneos eneos
Analizaremos ahora el caso de sistemas lineales no homog´eneos eneos y veremos un m´etodo etodo – el m´etodo etodo de variac var iaci´ i´ on de par´ ametros o de constantes – que nos permite hallar las soluciones de un sistema sistema lineal no homog´ homog´eneo eneo si conocemos conocemos una base del conjunto conjunto de soluciones soluciones del sistema sistema homog´ hom og´eneo ene o aso a socia ciado. do. Teorema
3.3. La soluci´ on general del sistema lineal (2.9) tiene la siguiente forma
(3.7)
X (t) = X p (t) + X H H (t)
donde X p es una soluci´ on on particular del sistema (2.9) y X H on general del sistema H es la soluci´ lineal homog´ h omog´ eneo eneo asociado, as ociado, es decir, de d e (2.9) (2 .9) pero con b con bi = 0. ´ n. Demostracion. o
X
−
Sea X p una soluci´ on particular de (2.9) y sea X otra soluci´ on on. o n. Sea Sea Y = X p . Entonces, si A = (aij ) es la matriz asociada al sistema (2.9) se sigue que Y = X
− X p = A(t)X + X + b(t) −
A(t)X p + b(t) = A(t)(X )(X
− X p) = A(t)Y.
Por lo tanto, Y es una soluci´ on on del sistema homog´eneo eneo asociado. aso ciado. Es decir, una soluci´ on on de Y = A(t)Y.
(3.8)
Sea ahora X = X p + Y donde Y es una soluci´on on de (3.8), entonces X = X p + Y = A(t)X p + b(t) + A(t)Y ( Y (t) = A(t)(X )(X p + Y ) Y ) + b(t) = A(t)X + X + b(t). Es decir, X es soluci´ on on de (2.9).
Veamos ahora el m´etodo eto do de variaci´ on on de constantes. Rn×n continua en un intervalo abierto I . Sea X 1 , 3.4. Sea A(t) , X n n una base del conjunto conjunto de soluciones soluciones de (3.8). Sea Sea bb(t) R continuo en I . Existen Existen funciones funciones c1 (t), , cn (t) continuamente diferenciables en I tales que Teorema
···
∈
∈
X p (t) = c1 (t)X 1 (t) +
{ ···
}
· · · + cn(t)X n(t)
es una soluci´ on particular del sistema X = A(t)X + X + b(t). M´ as precisamente, las funciones ci (t) son primitivas de las soluciones ci (t) del sistema lineal de ecuaciones (para cada t)
Q(t)
c1 (t) c2 (t) .. .
= b(t),
cn (t)
donde Q(t) es la matriz formada por los vectores X j (t) puestos como columnas. ´ n. Demostracion. o
Recordemos que en la observaci´on on 3.1 vimos que si tomamos la matriz Q(t) cuyas columnas son los vectores X j (t) – llamada matriz fundamental – se tiene Q (t) = A(t)Q(t) y la soluci´on on general del sistema (3.8) es X H H (t) = Q(t)C,
30
3. SIST SISTEM EMAS AS LINE LINEAL ALES ES DE 1ER. 1ER. OR ORDE DEN N Y ECUA ECUACI CION ONES ES LINE LINEAL ALES ES DE OR ORDE DEN N n
donde C es un vector constante. (Esto es exactamente decir que X H H (t) = c1 X 1 (t)+ Por lo tanto, lo que se est´a proponiendo es tomar
· · ·+cnX n(t)).
X p (t) = Q(t)C (t), donde ahora reemplazamos el vector constante C por un vector cuyas componentes son funciones continuamente diferenciables en I . Para ver que de esta manera podemos hallar una soluci´on on de X = A(t)X + X + b(t), simplemente derivamos derivamos y vemos qu´ e tiene que satisfacer C (t). Por la regla de deriv derivaci´ aci´ on on del producto (que sigue valiendo en este caso de producto de una matriz por un vector), X p (t) = Q (t)C (t) + Q(t)C (t) = A(t)Q(t)C (t) + Q(t)C (t) = A(t)X p (t) + Q(t)C (t) = A(t)X p (t) + b(t), si Q(t)C (t) = b(t).
(3.9)
∈
Recordemos que det Q(t) = 0 para todo t I por ser una matriz fundamental (ver Observaci´ on 3.1). Podemos por lo tanto invertir la matriz Q(t) y obtenemos on C (t) = Q(t)−1 b(t).
(3.10)
De esta manera obtenemos funciones continuas – las componentes del vector Q(t)−1 b(t) – que deber´ deber´ıan ser las componentes del vector C (t) para que X p sea una soluci´on on del sistema no homo ho mog´ g´eneo en eo X = A(t)X + X + b(t). Lo unico u ´ nico que resta ahora para encontrar las funciones ci (t) (componentes del vector C ) es integrar componente a componente (3.10) para obtener funciones continuamente diferenciables en I . Observemos que al integrar uno tiene constantes de integraci´on on arbitrarias. arbitrarias. Podemos Podemos simplemente tomarlas igual a 0 para obtener una soluci´on on particular. Si las dejamos en la expresi´ on de las funciones ci (t) obtenemos directamente la soluci´ on on general del sistema no homog´eneo eneo ya ¯ que el t´ermino ermino adicional adicio nal Q(t)C que obtenemos es la soluci´ on on general del sistema homog´eneo eneo asociado. Veremos ejemplos de aplicaci´on on de este m´etodo eto do al a l final fin al del cap´ cap´ıtulo siguiente sigui ente donde do nde aprenderemos a hallar la soluci´ on on general de sistemas lineales homog´eneos eneos con coeficientes constantes. En la pr´ actica actica lo que se hace hace es plant plantear ear (3.9) – que es, para para cada cada t, un sistema lineal de ecuaciones con inc´ognitas ognitas c1 (t), , cn (t) –, resolver el sistema y finalmente integrar el resultado para obtener obtener c1 (t), , cn (t).
···
···
Veamos ahora c´ omo se aplica omo ap lica el m´etodo etod o de variaci´ on on de par´ametros ametros para resolver ecuaciones lineales lineal es no n o homog´ ho mog´eneas eneas de orden o rden n. Consideremos Consid eremos la ecuaci´on on (3.11)
x(n) + an−1 (t)x(n−1) + an−2 (t)x(n−2) +
{
}
· · · + a1(t)x + a0(t)x = f ( f (t)
y supongamos que tenemos una base x1 , . . . , xn de soluciones de la ecuaci´on on homog´ hom og´enea enea asoaso ciada.
´ NEOS 2. SISTEMAS NO HOMOG E (n−1)
31
Tomemos las funciones X i (t) = (xi , xi , . . . , xi ). Sabemos que X 1 , . . . , Xn es una base de soluciones del sistema lineal homog´eneo eneo de 1er. orden asociado al siguiente sistema
{
}
x0 = x1 , x1 = x2 , .. .
(3.12)
xn −2 = xn−1 , xn −1 =
−a0(t)x0 − a1(t)x1 − · · · − an−1(t)xn−1 + f ( f (t).
El sistema (3.12) es un sistema no homog´eneo eneo con no homogeneidad
b(t) =
0 0 .. .
.
0 f ( f (t)
El m´etodo eto do de variaci´ varia ci´on on de par´ametros ametros para el sistema (3.12) dice que una soluci´ on particular tiene la forma X p (t) = c1 (t)X 1 (t) + + cn (t)X n (t),
···
donde las funciones c1 , . . . , cn son soluci´ on on del sistema (3.13)
Q(t)
c1 (t) .. .
=
cn (t)
0 .. .
.
f ( f (t)
Aqu´ı la matriz mat riz Q(t) tiene por columnas los vectores X j (t). Por lo tanto el sistema (3.13) es c1 (t)x1 (t) + c2 (t)x2 (t) + (3.14)
· · · + cn (t)xn(t) = 0, c1 (t)x1 (t) + c2 (t)x2 (t) + · · · + cn (t)xn (t) = 0, c1 (t)x1 (t) + c2 (t)x2 (t) + · · · + cn (t)xn (t) = 0, .. .
(n−1)
c1 (t)x1
(n−1)
(t) + c2 (t)x2
(t) +
f (t). · · · + cn (t)xn(n−1)(t) = f (
Como la primer componente del vector X p (t) es una soluci´on on particular de (3.11) se sigue que si las derivadas de las funciones c1 , , cn satisfacen satis facen (3.14), (3.14 ), la funci´on on
···
x p (t) = c1 (t)x1 (t) +
· · · + cn(t)xn(t)
es soluci´ on on de (3.11). Tenemos por lo tanto el siguiente teorema 3.5 (M´etodo etod o de variaci´ on on de par´ametros ametros para una ecuaci´on on lineal linea l no homog´enea enea de orden n). La soluci´ on general de la ecuaci´ on (3.11) tiene la forma Teorema
x(t) = x p (t) + xH (t),
32
3. SIST SISTEM EMAS AS LINE LINEAL ALES ES DE 1ER. 1ER. OR ORDE DEN N Y ECUA ECUACI CION ONES ES LINE LINEAL ALES ES DE OR ORDE DEN N n
donde x p es una soluci´ on particular de (3.11) y xH es la soluci´ on general de la ecuaci´ on homog´enea en ea asocia as ociada da.. Una soluci´ on on particular tiene la forma x p (t) = c1 (t)x1 (t) +
· · · + cn(t)xn(t),
{ · · · , xn} es
donde las derivadas de las funciones c1 , . . . , cn satisfacen el sistema (3.14) y x1 , una base de soluciones de la ecuaci´ on homog´ hom og´ enea enea asociada asoci ada..
CAP´ıTULO ıT ULO 4
Resoluci´ on on de sistemas lineales con coeficientes constantes En este cap´ cap´ıtulo buscaremos soluciones de sistemas lineales con coeficientes constantes. En el caso de sistemas de 2 2 hallaremos la soluci´ on on general. En dimensiones dimensiones may mayores, ores, dejaremos dejaremos indicada la idea de c´ omo son las soluciones. Comenzaremos por omo p or el caso homog´eneo eneo y despu´ es es veremos c´omo omo encontrar la soluci´ on on general del no homog´eneo eneo a partir de la soluci´ on on general del homo ho mog´ g´eneo en eo..
×
En todo el cap´ cap´ıtulo estudiaremos el sistema X = AX,
(4.1) donde A
∈ Rn×n es una matriz constante.
En el caso n = 1 ya sabemos c´ omo es la soluci´ omo solucion o´n general general.. En efecto, efecto, el sistem sistemaa (4.1) (4.1) se reduce a x = ax cuya soluci´ on on general es x(t) = ceat con c
∈ R arbitraria.
Motivados por el caso unidimensional, veamos si para sistemas hay soluciones de la forma X (t) = ξe λt ,
λ
∈ R,
ξ
∈ Rn.
Calculemos: X = ξλe λt ,
AX = (Aξ) Aξ )eλt
de modo que X = AX si y s´olo ol o si Aξ = λξ. λξ . Para tener soluci´ on on no trivial quiero que ξ = 0. De mod modo o que X (t) = ξe λt es soluci´ on on no trivial de (4.1) si y s´olo olo si λ es autovalor de A y ξ es un autovector asociado al autovalor λ.
De modo que si existe una base de Rn formada por autovectores de A, es decir, si existen ξ1 , . . . , ξn linealmente independientes tales que
{
}
Aξi = λi ξi
∈
para ciertos λi R (en otras palabras, la matriz A es diagonalizable), se tiene una base de soluciones del sistema lineal homog´eneo eneo (4.1) formada por p or las funciones X i (t) = ξi eλi t . Observemos que los λi no tienen por qu´ e ser distintos distintos para distintos distintos i. Las Las solu soluci cion ones es λ t X i = ξi e i son linealmente independientes porque lo son en t = 0 ya que X i (0) = ξi . En este caso la soluci´ on on general de (4.1) es n
X (t) =
i=1
con ci
∈ R constantes arbitrarias.
33
ci eλi t ξi
´ DE SISTEMAS 4. RESOLUCION SISTEMAS LINEALES LINEALES CON COEFICIEN COEFICIENTES TES CONSTANTES CONSTANTES
34
Un caso particular en donde estaremos en esta situaci´on on es si todos los autovalores de A son reales y distintos ya que autovectores correspondientes a autovalores distintos siempre son linealmente independientes. Pero, como dijimos antes, podr´ıamos ıamos tener una base de autovectores aunque A no tenga todos los autovalores distintos. Ejemplo
4.1. Sea A=
Hallar las soluciones del sistema X = AX .
1 2 . 2 1
∈
Para hallar los autovalores debemos buscar los λ
− λ
−2 λ−1
1
−2
tales que
R
x1 x2
=0
tenga una soluci´ on on no trivial. Para ´esto esto es necesario y suficiente que p( p(λ) = det
− λ
−2 λ−1
1
−2
= 0.
−
Recordemos que p(λ) se llama el polinomio polin omio caracter´ıstico ısti co de la matriz mat riz A A ( p( p(λ) = det(λI det(λI A)). En este caso, p(λ) = (λ
− 1)2 − 4. Por lo tanto, los autovalores son λ1 = 3 y λ2 = −1.
Busquemos ahora autovectores asociados a estos dos autovalores. Para λ1 = 3, un autovector ser´ a soluci´ on on de (3I (3I A)ξ = 0, es decir, tendr´ a componentes (x1 , x2 ) tales que
−
Es decir, x1 la solu so luci ci´´on on
2 2
−
−2 2
x1 x2
= 0. 0.
− x2 = 0. Por lo tanto, un autovector asociado a λ1 = 3 es ξ1 = (1, (1, 1) y tendremos X 1 (t) = e
3t
1 . 1
−1 debemos resolver el sistema (−I − A)ξ = 0, es decir −2 −2 x1 = 0,0 , −2 −2 x2 que equivale a x1 + x2 = 0 y tiene por soluci´on on a ξ2 = (1, (1, −1), y nos da la soluci´ on on Para el caso de λ2 =
X 2 (t) = e−t
1 . 1
−
De aqu´ aqu´ı que la soluci´ soluci on o´n general es (4.2)
X (t) = c1 e3t
1 + c2 e−t 1
En este ejemplo el sistema que quer´ıamos ıamos resolver es: x1 = x1 + 2x 2x2 , x2 = 2x1 + x2 .
1 . 1
−
´ N DE SIST 4. RESOLUCIO SISTEM EMAS AS LIN LINEALE EALES S CON CON COE COEFICI FICIEN ENTE TES S CONS CONST TANTES NTES
35
De la f´ormula ormula (4.2) encontramos que la soluci´ on general en coordenadas x1 , x2 es on x1 = c1 e3t + c2 e−t , x2 = c1 e3t con c1 y c2 constantes arbitrarias.
− c2e−t,
En much muchas as situac situacion iones es ocurre ocurre que los autov autovalo alores res de la matriz matriz son n´ umeros umeros complejos complejos n (aunque la matriz sea real) y por ende, sus autovectores asociados pertenecen a C . En es esn n ta situaci´ situaci´ on conviene pensar que estamos trabajando en C en lugar de R . Es decir, admitimos on soluciones que tomen valores complejos. La ventaja es que cuando A tiene autovalores comple jos, puede ser diagonalizable en Cn pero no lo ser´a en Rn . Para Para ver que podemos podemos trabajar trabajar en n definamos la exponencial de un n´ umero umero complejo. C ´n Definicion o
4.1. Sea λ = α + iβ
∈ C, definimos
eλ = eα (cos β + β + i sen β ). Es f´acil acil ver que la exponencial exponencial as´ as´ı definida definida sigue satisfaciendo satisfaciendo las propiedades propiedades de la expoλ + λ λ λ nencial real, por ejemplo, e = e e . Adem´ as as se tiene 1
2
1
2
d λt e = λeλt . dt En efecto, d λt d αt e = e (cos βt + i sen βt) βt ) = αeαt (cos βt + i sen βt) βt ) + eαt β ( sen βt + i cos βt) βt ) dt dt = eαt (α cos βt β sen β sen βt) βt ) + i(α sen βt + β cos β cos βt) βt ) = eαt (α + iβ )(cos )(cos βt + i sen βt) βt )
= λeλt .
−
−
Por lo tanto, si hay una base de Cn formada por autovectores de A, tendremos una base del conjunto de soluciones complejas del sistema X = AX construida como en el caso real. Ejemplo
4.2. Hallar las soluciones de x1 = x1
− x3,
x2 = x1 , x3 = x1 En este caso A=
1 1 1
− x2.
0 0 1
−
Para hallar los autovalores planteamos
− −− λ
λI
−A=
−1 0 0
.
1 0 1 1 λ 0 1 1 λ
,
´ DE SISTEMAS 4. RESOLUCION SISTEMAS LINEALES LINEALES CON COEFICIEN COEFICIENTES TES CONSTANTES CONSTANTES
36
con lo cual resulta p(λ) = (λ 1)λ 1)λ2 1 + λ = (λ 1)λ 1)λ2 + (λ 1) = (λ (λ 1)(λ 1)(λ2 + 1) y los autovalores (las raices de p(λ)) son λ1 = 1, λ2 = i y λ3 = i. Por lo tanto, al ser los 3 distintos, hay una base de Cn formada por autovectores de A y podremos encontrar una base de soluciones del sistema diferencial si admitimos que tomen valores complejos.
−
−
−
−
−
−
Empecemos por buscar autovectores asociados a los autovalores hallados. Para λ1 = 1
−−
− −−
0 0 1 1 1 0 1 1 1
x1 x2 x3
= 0.
Es decir, x3 = 0, x1 = x2 . Una soluci´ on on es (1, (1, 1, 0). Por lo tanto una soluci´ on on del sistema es X 1 (t) = e Para λ2 = i i
Es decir
t
1 0 1 1 i 0 1 1 i
1 1 0
.
x1 x2 x3
= 0.
− 1)x 1)x1 + x3 = 0, − x1 + ix2 = 0,
(i
ya que la tercer ecuaci´on on es combinaci´ on lineal de las 2 primeras. Una soluci´ on solucion o´n es x1 = i, x2 = 1, x3 = 1 + i que nos da la soluci´on on compleja del sistema diferencial X 2 (t) = eit
i 1 1+i
= (cos t + isen t)
Finalmente, para el autovalor λ3 =
− − −− i
1 1
i 1 1+i
−sen t + i cos t
=
(cos t
−i, 1
−
0 i 1
1 0 i
−
x1 x2 x3
Es decir ( i 1)x 1)x1 + x3 = 0, x1 ix2 = 0,
−− − −
que tiene por soluci´on on a x1 =
−
−i, x2 = 1, x3 = 1 − i.
Antes de proseguir, observemos que el autovector
− − i 1
1
i
= 0.
cos t + isen t sen t) + i(cos t + sen t)
.
´ N DE SIST 4. RESOLUCIO SISTEM EMAS AS LIN LINEALE EALES S CON CON COE COEFICI FICIEN ENTE TES S CONS CONST TANTES NTES
37
−
asociado al autovalor i es el conjugado del autovector que hallamos asociado al autovalor i. Y el autovalor i es el conjugado del autovalor i. De modo, que la soluci´on on que estamos encontrando
−
−it
X 3 (t) = e
− − i 1
1
,
i
es la conjugada de X 2 (t) (Recordemos que t es real). real). Por lo tanto, tanto, si escr escrib ibim imos os X 2 (t) = X R (t) + iX I I (t) con X R y X I I reales, tendremos que X 3 (t) = X R (t) iX I I (t).
−
Como
1 X R = (X 2 + X 3 ) 2 se sigue que X R y X I I son soluciones reales.
y X I I =
{
1 (X 2 2i
− X 3),
}
Es f´acil acil ver en este ejemplo que X 1 , X R , X I I son linealmente linealmente independien independientes tes y por lo tanto forman forman una base base del espacio espacio de soluci solucione oness reales reales del sistem sistemaa difere diferenci ncial. al. Ensegu Enseguida ida verem veremos os que ´este este es un hecho general y que (cuando A es real) siempre podemos encontrar una base de soluciones reales a partir de una base de soluciones complejas ya que ´estas estas vendr´ an an de a pares conjugados. Para terminar de encontrar la soluci´ on general en este ejemplo observemos que on X R (t) =
− − − − sen t cos t cos t sen t
,
X I I (t) =
cos t sen t cos t + sen t
Por lo tanto la soluci´ on on general real es X (t) = c1 et En coordenadas,
1 1 0
+ c2
x1 = c1 et
sen t cos t cos t sen t
+ c3
.
cos t sen t cos t + sen t
.
− c2sen t + c3 cos t,
x2 = c1 et + c2 cos t + c3 sen t, x3 = (c2 + c3 )cos t
− (c2 − c3)sen t.
Veamos entonces que lo observado en el ejemplo 4.2 es un hecho general y que el mismo procedimiento se puede aplicar a toda matriz real A que posea un autovalor complejo.
En efecto, supongamos que λ = α + iβ con β = 0 es un autovalor de A y sea ξ = w + iv, iv , con w y v vectores reales, un autovector asociado. Es decir, Aξ = λξ. Conjugando componente a componente, y como A es real, tenemos ¯ ξ. ¯ Aξ¯ = λ ¯ es un autovalor de A asociado al autovector ξ¯. Esto nos da dos soluciones complejas Es decir, λ conjugadas, a saber ¯ ¯ 1 (t). X 1 (t) = eλt ξ , X 2 (t) = eλt ξ¯ = X
38
´ DE SISTEMAS 4. RESOLUCION SISTEMAS LINEALES LINEALES CON COEFICIEN COEFICIENTES TES CONSTANTES CONSTANTES
Como en el ejemplo 4.2, podemos realizar combinaciones lineales de X 1 = X R (t) + iX I I (t) y X 2 = X R (t) iX I I (t) para obtener soluciones reales. Tenemos que 1 1 X R (t) = (X 1 (t) + X 2 (t)), )), X I I (t) = (X 1 (t) X 2 (t)), )), 2 2i son soluciones soluciones del sistema. Veamos que X R y X I I son linealmente independientes (recordemos que X 1 y X 2 lo son porque son independientes en t = 0 al ser ξ y ξ¯ autovectores correspondientes a autovalores distintos).
−
−
Veamos que la unica u ´nica combinaci´ on on lineal de X R y X I I que da la funci´on on 0 es aquella que tiene las dos constantes nulas. En efecto, si c1 X R + c2 X I I = 0 se tiene c1 c2 c1 c2 c1 c2 0 = (X 1 + X 2 ) + (X 1 X 2 ) = ( + )X 1 + ( )X 2 . 2 2i 2 2i 2 2i Por lo tanto, 1 1 c1 + c2 = 0, 2 2i 1 1 c1 c2 = 0. 2 2i
−
−
−
Es f´acil acil ver que la ´unica unica soluci´ sol uci´on on posible p osible de este sistema para las constantes c1 , c2 es c1 = c2 = 0. De esta manera sabemos c´omo omo encontrar una base de soluciones reales en el caso de dimensi´on on 2 si los dos autovalores de la matriz A son distintos. La misma idea funciona en m´as as dimensiones. dimensiones. Los autovalores autovalores complejos complejos vienen de a pares conjugados y si en una base de soluciones complejas reemplazamos un par de soluciones con jugadas por la parte real y la parte imaginaria de una de ellas, se sigue teniendo una base de soluciones. Esto es lo que hicimos en el ejemplo. En dimensi´on o n 2 el unico u ´ nico caso que nos queda por analizar es el de un autovalor real de multiplicidad 2. Sea λ este autovalor. autovalor. Si hay una base de autovectores autovectores,, estamos estamos en el caso que sabemos resolver resolver.. Pero Pero en realidad realidad es un caso muy simple, ya que tendremos tendremos A = λI y por lo tanto es un sistema desacoplado x1 = λx1 , x2 = λx2 , cuya soluci´ on on general es x1 (t) = c1 eλt , x2 (t) = c2 eλt con c1 y c2 constantes arbitrarias. Si no hay una base de autovectores, s´olo olo tenemos una soluci´ on o n de la forma X (t) = eλt ξ . ¿C´ omo encontrar otra soluci´ omo on linealmente independiente? on Para tratar de entender cu´al al puede ser la forma de la segunda soluci´ on on del sistema, estudiemos tudiemos primero primero el siguiente siguiente ejemplo: Ejemplo
4.3. Hallar una base de soluciones del sistema λ 0 X = X 1 λ
En este caso, la matriz tiene un s´olo olo autovalor λ con autovector asociado (0, (0, 1) y no hay ning´ un un otro autovector linealmente independiente. Con lo cual el m´etodo etodo desarrollado hasta el momento no nos permite hallar todas las soluciones de la ecuaci´on. on.
´ N DE SIST 4. RESOLUCIO SISTEM EMAS AS LIN LINEALE EALES S CON CON COE COEFICI FICIEN ENTE TES S CONS CONST TANTES NTES
39
Sin embargo, en este ejemplo, el sistema se “desacopla” de la siguiente manera. La ecuaci´ on on para x1 es x1 = λx1 , con lo cual x1 (t) = c1 eλt . Ahora, la ecuaci´on on para x2 resulta x2 = x1 + λx2 = c1 eλt + λx2 , que es una ecuaci´ on on lineal no homog´enea enea de orden 1. Es sencillo entonces hallar la forma general para x2 y la misma es x2 = (c1 t + c2 )eλt . En s´ıntesis, ıntesis, tenemos que la l a soluci´ so luci´on on general de la ecuaci´on on es X =
c1 eλt (c1 t + c2 ) eλt
x1 x2
=
Es decir la base de soluciones es
eλt
λt
= c2 e
0 , eλt 1
0 + c1 eλt 1
0 1 t+ 1 0
0 1 t+ 1 0
.
.
Volviendo al caso general de una matriz con un ´unico unico autovalor y espacio de autovectores de dimensi´ dimensi´ on 1, y basados en el ejemplo 4.3, buscamos una segunda soluci´on on on en la forma X (t) = eλt (ξ1 t + ξ2 ), para ciertos vectores ξ1 , ξ2 . Derivando tenemos X (t) = λeλt (ξ1 t + ξ2 ) + eλt ξ1 . Por otro lado,
AX (t) = eλt Aξ1 t + eλt Aξ2 .
Para que X (t) = AX (t) para todo t se debe tener Aξ1 = λξ1 , Aξ2 = λξ2 + ξ1 . Como λ es autovalor de A podemos tomar como ξ1 un autovector asociado asoci ado y con ´esto se satisface la primer ecuaci´on. on. Elegido ξ1 , la segunda ecuaci´on on pasa a ser un sistema lineal no homog´eneo eneo con matriz del sistema sistem a A λI singular.
−
Veamos que este sistema siempre tiene una soluci´on. o n. En efec efecto to,, sea η un vector tal que linealmente te indepen independie dient ntes. es. Por Por lo tanto, tanto, η, ξ1 es una base de R2 . Enton Entonces ces,, ξ1 y η son linealmen Aη = c1 η + c2 ξ1 . De aqu´ı que la matriz de A en la base η, ξ1 es
{ {
A{η,ξ
1
}
=
c1 0 . c2 λ
} }
Como λ es el unico u ´ nico autovalor de A, el polin´ omio omio caracter caract er´´ıstico ıstic o de A es p(r) = (r lo tanto, c1 = λ. Finalmente, tomamos ξ2 = η/c2 y tenemos 1 1 Aξ2 = Aη = (λη + c2 ξ1 ) = λξ2 + ξ1 , c2 c2 comoo quer´ıamos. com ıam os.
− λ)2. Por
40
´ DE SISTEMAS 4. RESOLUCION SISTEMAS LINEALES LINEALES CON COEFICIEN COEFICIENTES TES CONSTANTES CONSTANTES
Observemos que c2 = 0, y por p or lo l o tanto ta nto podem p odemos os dividir di vidir por ´el, el, ya que si no fuera fu era as´ as´ı se tendr´ıa ıa que η ser´ıa ıa tambi´ tamb i´en en autovec a utovector tor y habr ha br´´ıa un base b ase de autovec a utovector tores, es, lo que hab´ıamos ıam os supuest supu esto o que qu e no pasaba. ´n Observacion o
Ejemplo
4.1. Lo que hemos hecho es hallar una base de Jordan de la matriz A.
4.4. Hallar una base de soluciones para el sistema
x1 = x1
x2 = x1 + 3x 3x2 .
La matriz del sistema es A= que tiene polin´omio omio caracter caracter´istico p( p(λ) = det
− λ
1
−1
1 λ
−3
− x2,
= (λ
− 1 1
1 , 3
− 1)(λ 1)(λ − 3) + 1 = λ2 − 4λ + 4 = (λ (λ − 2)2 .
Por lo tanto A tiene un unico u ´ nico autovalor λ = 2. Busque Busquemos mos los los autove autovecto ctores res.. Debemos Debemos resolver el sistema x 1 1 x1 (2I (2I A) 1 = = 0. x2 1 1 x2 Es decir, x1 + x2 = 0. Por lo tanto un autovector asociado es
− − −
−
1 1
ξ1 =
y tenemos la soluci´on on
1 . 1
X 1 (t) = e2t
−
Como no hay una base de autovectores y el autovalor es doble, buscamos la otra soluci´on on 2 t de la forma X 2 (t) = e (ξ1 t + ξ2 ) donde ξ1 es el autovector encontrado y ξ2 es soluci´ on o n de Aξ2 = 2ξ2 + ξ1 . Es decir, ξ2 es soluci´ on on de (A es decir, x1 + x2 =
− 2I )
− − − − x1 x2
=
1 1
1 1
x1 x2
1 , 1
=
−1. Por ejemplo, podemos tomar
0 , 1
ξ2 =
con lo cual tenemos la soluci´ on on
X 2 (t) = e2t
1 1
−
0 1
t+
−
.
La soluci´ on on general es entonces 2t
X (t) = c1 e
1 + c2 e2t 1
−
1 1
−
t+
0 1
−
.
´ N DE SIST 4. RESOLUCIO SISTEM EMAS AS LIN LINEALE EALES S CON CON COE COEFICI FICIEN ENTE TES S CONS CONST TANTES NTES
41
En coo coordenadas rdenadas ´esto esto es, x1 = e2t (c1 + c2 t), x2 =
Podemos aplicar esta misma idea en m´ as as dimensiones. Ejemplo
x1 = x1 + x2 , x2 = x2
− x3,
x3 = x2 + 3x 3x3 .
A=
−
1 1 0 1 0 1
0 1 3
que tiene polin´omio omi o caract car acter er´´ıstico ıst ico p(λ) igual a λ
det
0 0
1
4.5. Hallar una base de soluciones para el sistema
La matriz del sistema es
−
−e2t (c1 + c2) + c2t .
−1 0 λ−1 1 −1 λ − 3
= (λ
,
(λ − 1) = (λ (λ − 1)(λ 1)(λ2 − 4λ + 4) = (λ (λ − 1)(λ 1)(λ − 2)2 . − 1)2(λ − 3 ) + (λ
Por lo tanto los autovalores son λ1 = 1 y λ2 = 2, este ultimo u ´ ltimo con multiplicidad 2. Busquemos los autovec autovectores tores asociados a λ1 . Debemos resolver (I
− A)
x1 x2 x3
=
−1 0 0 1 −1 −2
0 0 0
x1 x2 x3
= 0,
lo que equivale a x2 = x3 = 0. Por lo tanto un autovector es ξ1 = que da la soluci´on on
1 0 0
X 1 (t) = e Para el autovalor λ2 = 2 debemos resolver (2I (2I
−
− A)
x1 x2 x3
=
1 0 0
,
1 0 0
t
.
−1 0 1 1 −1 −1
x1 x2 x3
= 0,
lo que equivale a x1 x2 = 0, x2 + x3 = 0. Por Por lo tanto, tanto, no hay hay 2 autov autovect ectore oress linealme linealment ntee independientes asociados a este autovalor doble. Un autovector es ξ2 =
− 1 1 1
,
42
´ DE SISTEMAS 4. RESOLUCION SISTEMAS LINEALES LINEALES CON COEFICIEN COEFICIENTES TES CONSTANTES CONSTANTES
que da la soluci´on on 2t
X 2 (t) = e
− 1 1 1
.
Para hallar otra soluci´on on linealmente independiente la buscamos de la forma X 3 (t) = e2t (ξ2 t + ξ3 ),
− − − − − − −
donde ξ2 es el autovector que hallamos y ξ3 es soluci´ sol uci´on on de (A
−
1 0 0
lo que equivale a
1
0 1 1
−1 1
−x1 + x2 = 1, x2 + x3 =
x1 x2 x3
2t
X 3 (t) = e
=
,
1. Una soluci´ on on es 0 1 2
ξ3 = y obtenemos la soluci´on on
1 1 1
2I )ξ3 = ξ2, es decir soluci´on on de
1 1 1
0 1 2
t+
.
De este modo encontramos encontramos 3 soluciones soluciones linealmente linealmente independientes. independientes. La soluci´ on on general es entonces 1 1 1 0 t 2t 2t 1 + c3 e 1 t+ 1 X (t) = c1 e 0 + c2 e . 0 1 1 2 lo que en coordenadas da x1 = c1 et + (c (c2 + c3 t)e2t
−
− −
x2 = (c2 + c3 ) + c3 t e2t x3 =
− (c2 + 2c 2c3 ) + c3 t e2t
×
En el caso de sistemas de 2 2 ya hemos considerado considerado todos los casos posibles. Pero Pero en el caso de sistemas de 3 3 todav´ todav´ıa podr´ podr´ıamos tener un autov a utovalor alor triple. En alguno de estos casos a´ un sabemos c´ omo omo resolver el sistema. Esto es as´ as´ı si hay 3 autovectores linealmente independientes o tambi´ tamb i´en en si s´ olo hay dos autovectores linealmente independientes ya que este caso lo tratamos olo como arriba. arriba. Sin embargo, embargo, este caso es m´ as complicado que el correspondiente en dimensi´on as on 2 o el del ultimo u ´ ltimo ejemplo ya que el espacio de autovectores tiene dimensi´ on on 2 y no est´a claro cu´ al al es el autovector que multiplica a t en la tercer soluci´on. on.
×
Pero bien podr´ıa ıa suceder que s´olo olo podamos encontrar 1 autovector y todos los dem´ as as sean m´ ultiplos ultiplos de ´el. el. No vamos a discutir discutir este caso, pero se ve que a medida medida que crece la dimensi´ dimensi´ on on son m´as as los casos que pueden aparecer. Para analizar el caso general en dimensi´on on n necesit nece sitamo amoss m´as as algebra a´lgebra lineal y m´ as as tiempo.
´ N DE SIST 4. RESOLUCIO SISTEM EMAS AS LIN LINEALE EALES S CON CON COE COEFICI FICIEN ENTE TES S CONS CONST TANTES NTES
43
Veamos ahora c´omo omo funciona el m´etodo etodo de variaci´ on de constantes para resolver sistemas lineales lineal es no homog´eneos. eneos. Ejemplo
4.6. Hallar la soluci´ on on general del siguiente sistema no homog´eneo eneo
x1 = x1 + 2x 2x2 + et , x2 =
t
e −2x1 + x2 + cos2t . cos2t
Seg´ un un el m´etodo etod o de variaci´on on de constantes tenemos que encontrar la soluci´on on general del homog´eneo eneo para despu´es es encontrar una soluci´ solucion o´n particular. particular. Buscamos Buscamos entonces entonces los autov autovalores de la matriz del sistema, es decir, las ra´ ra´ıces del polinomio
−
−2 = (λ − 1)2 + 4.4. p( p(λ) = det 2 λ−1 Los autovalores son, entonces, λ1 = 1 + 2i 2i, λ2 = 1 − 2i. λ
1
Busquemos un autovector asociado al autovalor λ1 . Esto es, resolvamos el sistema
(1 + 2i 2i)I
Esto es, 2ix 2ix1
− x1 x2
−A
=
2i 2
2 2i
x1 x2
= 0.
− 2x2 = 0 que tiene por soluci´on on al vector ξ1 =
y nos da la soluci´on on compleja compleja
1 i
(1+2i)t X 1 = e(1+2i
1 . i
Vimos que una base de soluciones reales est´a formada por la parte real y la parte imaginaria de X 1 . Por lo tanto, vamos a encontrarlas.
1 X 1 (t) = e (cos2t (cos2t + i sen2t sen2t) = et i t
Por lo tanto X R (t) = et
cos2t cos2t , sen2t sen2t
−
X H H (t) = c1 e
−
X I I (t) = et
y la soluci´on on general del sistema homog´eneo eneo asociado es t
cos2t cos2t + i sen2t sen2t . sen2t sen2t + i cos2t cos2t
sen2t sen2t cos2t cos2t
cos2t cos2t sen2t sen2t + c2 et . sen2t sen2t cos2t cos2t
−
Ahora buscamos una soluci´on on particular del sistema no homog´eneo eneo de la forma X p (t) = c1 (t)et
cos2t cos2t sen2t sen2t + c2 (t)et . sen2t sen2t cos2t cos2t
−
Las funciones c1 (t) y c2 (t) las buscamos de modo que la funci´on on C (t) =
c1 (t) , c2 (t)
´ DE SISTEMAS 4. RESOLUCION SISTEMAS LINEALES LINEALES CON COEFICIEN COEFICIENTES TES CONSTANTES CONSTANTES
44
satisfaga – con Q(t) la matriz cuyas columnas son las soluciones X R (t) y X I I (t) –
Q(t)C (t) = e
t
1 . 1/ cos2t cos2t
Es decir, c1 (t) y c2 (t) deben ser soluci´on on de
c1 (t) c2 (t)
et cos2t cos2t et sen2t sen2t t t e sen2t sen2t e cos2t cos2t
−
t
=e
1 . 1/ cos2t cos2t
Es decir, d ecir, soluci´ soluci ´on on del d el sistema sistem a et cos2t cos2t c1 + et sen2t sen2t c2 = et , t
e −etsen2t sen2t c1 + et cos2t cos2t c2 = . cos2t cos2t Multiplicand Multiplicandoo la primer primer ecuaci´ ecuaci´ on o n por cos2t cos2t, la segunda segunda por sen 2t y restando obtenemos t
2
e (cos 2t + sen
2
2t) c1 (t)
t
=e
Por lo tanto, c1 (t) = cos cos 2t
cos2t cos2t
−
sen2t sen2t . cos2t cos2t
sen2t , − sen2t cos2t cos2t
de donde, integrando, obtenemos 1 1 c1 (t) = sen2t sen2t + log(cos log(cos 2t). 2 2 Por otro lado, multiplicando la primer ecuaci´on on por po r sen 2t, la segunda segunda por cos 2t y sumando obtenemos et (sen 2 2t + cos2 2t) c2 (t) = et (sen2t (sen2t + 1). 1). Por lo tanto, c2 (t) = sen sen 2t + 1, 1, de donde, integrando, obtenemos c2 (t) =
− 12 cos2t cos2t + t.
As´ As´ı obtenemos obtene mos una soluci´ solucion o´n particular, a saber
−
1 1 cos2t cos2t X p (t) = sen2t sen2t + log(cos log(cos 2t) et + sen2t sen2t 2 2 1 cos2t cos2t log(cos log(cos 2t) 2 = et . 1 1 sen2t sen2t log(cos log(cos 2t) + t cos2t cos2t 2 2
−
−
−
1 sen2t sen2t cos2t cos2t + t et cos2t cos2t 2
EJERCICIOS
45
Finalmente, la soluci´ on general se obtiene como X (t) = X p (t) + X H on tanto, o, la H (t). Por lo tant soluci´ on on general es 1 1 1 cos2t cos2t sen2t sen2t X (t) = sen2t sen2t + log(cos log(cos 2t) + c1 et + cos2t cos2t + t + c2 et sen2t sen2t cos2t cos2t 2 2 2 1 cos2t cos2t log(cos log(cos 2t) + c1 cos2t cos2t + c2 sen2t sen2t 2 = et . 1 1 sen2t sen2t log(cos log(cos 2t) + t cos2t cos2t c1 sen2t sen2t + c2 cos2t cos2t 2 2
Es decir
−
−
−
−
−
1 x1 (t) = e cos2t cos2t log(cos log(cos 2t) + c1 cos2t cos2t + c2 sen2t sen2t , 2 1 1 x2 (t) = et + sen2t sen2t log(cos log(cos 2t) t cos2t cos2t + c1 sen2t sen2t 2 2 t
−
−
− c2 cos2t cos2t
.
Ejercicios
(1) Hallar la soluci´ on general de los siguientes sistemas on (a)
(c)
x1 = x2 x2 = 2x1 + 3x 3x2
−
x1 x2
= =
−4x1 + 3x 3x2 −2x1 + x2
(b)
(d)
x1 = 8x1 5x2 x2 = 10 10x x1 + 7x 7x2
− −
x1 = x2 = x3 =
−x1 + 3x 3x2 − 3x3 −2x1 + x2 −2x1 + 3x 3x2 − 2x3
En cada caso, hallar el conjunto de datos iniciales tales que la soluci´on on correspondiente tiende a 0 cuando t tiende a + . Idem con t tendiendo tendiendo a . (2) Inicialmente el tanque I contiene 100 litros de agua salada a una concentraci´on on de 1 kg por litro; el tanque II tiene 100 litros li tros de agua pura. El l´ l´ıquido se bombea b ombea del tanque I al tanque II a una velocidad de 1 litro por minuto, y del tanque II al I a una velocidad de 2 litros por minuto. Los tanques tanques se agitan constantemen constantemente. te. ¿Cu´ al es la concentraci´on on en el tanque I despu´ es es de 10 minutos? (3) Hallar la soluci´ on general de los siguientes sistemas on
∞
(a)
(c)
x1 = x1 x2 x2 = x1 + x2
(b)
x1 = 2x1 + x2 x2 = 2x2
(d)
−
−∞
x1 = 2x1 x2 x2 = 4x1 + 2x 2x2
−
x1 = x2 =
−5x1 + 9x 9x2 −4x1 + 7x 7x2
(4) Hallar la soluci´ on general de los siguientes sistemas on (a)
x1 = x2 + 2 x2 = 2x 2 x1 + 3x 3x2 + t
−
(b)
x1 = 2x1 x2 + e2t x2 = 4x1 + 2x 2x2 + 4
−
CAP´ıTULO ıT ULO 5
Ecuaciones lineales de orden
n
con coeficientes constantes
En este cap´ cap´ıtulo encontraremos la l a soluci´ on general de ecuaciones lineales con coeficientes on constantes de orden n. Comenzaremos con el caso homog´eneo. eneo. Teniendo en cuenta la relaci´on on entre ecuaciones de orden n y sistemas de n ecuaciones con n inc´ognitas ognitas vemos que la expresi´on on de la soluci´ on on de la ecuaci´on on (5.1)
x(n) + an−1 (t)x(n−1) + an−2 (t)x(n−2) +
· · · + a1(t)x + a0(t)x = 0.
depender´a de los autovalores de la matriz
(5.2)
A=
0 0 .. . .. . .. .
1 0
0 a0
0 a1
−
··· ···
0 1
0 0
··· 0 1 · · · −an−2 −an−1
0 a2
−
0 0
−
.
El polinomio p olinomio caracter´ caracter´ıstico de esta matriz es el polinomio p olinomio p( p(λ) = λn + an−1 λn−1 + an−2 λn−2 +
· · · + a1λ + a0,
que es el polinomio que se obtiene cuando se reemplaza en la ecuaci´ on on difere diferenci ncial al (5.1) (5.1) la derivaci´on on por la correspondiente potencia de λ. Lla Llamar maremo emoss a este este polinomio polinomio el polinomi polinomio o caracter caracte r´ıstico ısti co de la ecuaci´ ecuac i´ on . Por lo tanto, ta nto, la expresi´on on de las soluciones soluci ones depender´ dep ender´ a de las ra´ ra´ıces ıce s del polin po linomi omio o caract ca racter er´´ıstico ıst ico de la ecuaci´on. on. Por lo que vimos vi mos para par a sistemas, sistem as, si s i todas to das las ra r a´ıces λ1 , . . . , λn son reales y distintas, la soluci´on on general tiene la forma x(t) = c1 eλ t + c2 eλ t + + cn eλn t , 1
2
···
pues ´esta esta es la forma que tiene la l a primer componente comp onente de la soluci´on on general del sistema asociado. Esto nos dice que eλ t , eλ t , . . . , eλn t es una base de soluciones.
{
1
2
}
La misma expresi´ on on tiene la soluci´ on on general compleja si todas to das las l as ra´ıces ıces son s on distintas di stintas aunque algunas alguna s pueden ser s er complejas. compl ejas. Las ra´ıces ıces complejas co mplejas del polinom po linomio io caracter´ cara cter´ıstico ıstic o vienen de d e a pares conjugados. Por lo tanto, si λ = α + iβ con β = 0 es una ra´ız ız compleja, en la base compleja de soluciones aparecen dos soluciones conjugadas a saber
x ¯(t) = e(α−iβ )t .
x(t) = e(α+iβ )t , 47
48
5. ECUACIONES LINEALES DE ORDEN n CON COEFICIENTES CONSTANTES
Podemos reemplazar este par por un par de soluciones reales, a saber, la parte real y la parte imaginaria de x(t) ya que no perdemos la independencia lineal de las soluciones con este reemplazo (la justificaci´ on es la misma que para los sistemas). on Es decir, reemplazamos el par de soluciones complejas e(α+iβ )t , e(α−iβ )t por el par de soluciones reales Re( Re(e(α+iβ )t ) = eαt cos β t, I m(e(α+iβ )t ) = eαt sen βt. En el caso de ecuaciones de segundo orden, o bien se tiene dos ra´ ra´ıces reales distintas, disti ntas, o bien dos ra´ıces ıces complejas conjugadas (y por lo l o tanto distintas), o bien una unica ´ ra´ ra´ız de multiplicidad multipl icidad 2. Este ultimo u ´ ltimo es el unico u ´nico caso que resta analizar. Como sugiere la resoluci´ on de los sistemas, en este caso la soluci´ on on on general es x(t) = (c1 + c2 t)eλt , donde λ es la unica u ´nica ra´ ra´ız del polinomio caracter´ caracter´ıstico de la ecuaci´ on. on. Esto es as´ as´ı porque los λt λt autovalores de la matriz (5.2) son todos simples. Por lo tanto, e , te es una base de soluciones en este caso.
{
}
De modo que sabemos hallar una base de soluciones para cualquier ecuaci´ on on de orden 2. Despu´ es es continuaremos con ecuaciones de orden superior. sup erior. Veamos ahora algunos ejemplos Ejemplo
5.1. Hallar las soluciones de la ecuaci´on on x
− 2x + x = 0.
El polinomio pol inomio caracter caract er´´ıstico ıstic o es p( p(λ) = λ2 cuya unica u ´ni ca ra´ız ız es λ = 1.
− 2λ + 1 = (λ (λ − 1)2 ,
Por lo tanto, la soluci´ on on general es x(t) = (c1 + c2 t)et . Ejemplo
5.2. Hallar las soluciones de la ecuaci´on on x + x = 0. 0.
El polinomio pol inomio caracter caract er´´ıstico ıstic o es p( p(λ) = λ2 + 1, 1, cuyas cu yas ra´ıces ıc es son so n λ1 = i, λ2 = i. Por lo tanto tenemos el par de soluciones complejas conjugadas it it − e , e . Pero Pero podemos reempla reemplazar zarlas las por un par de solucion soluciones es reales, reales, la parte real y la parte it imaginaria imaginaria de e a saber cos t, sen t. Por lo tanto, la soluci´ on general en este caso es on
−
x(t) = c1 cos t + c2 sen t.
El caso general de la ecuaci´ on on de orden n (5.1) escapa, por problemas de tiempo, a nuestras posibili posibilidad dades es en el curso. curso. Sin embarg embargo, o, para para com comple pletitu titud d de estas notas y para para los alumnos alumnos interesados desarrollaremos este caso a continuaci´ on. on. El enfoque utilizado utilizado es independiente independiente de
5. ECUACIONE ECUACIONES S LINEALES LINEALES DE ORDEN n CON COEFICIENTES CONSTANTES
49
´ los sistemas de orden n para los que, como dijimos, hace falta m´ as as conocimientos de Algebra Lineal. Con este fin introduzcamos la noci´on on de operador diferencial. Primero definamos el operador de derivaci´on on D por la relaci´on on Dx = x . El operador D se aplica a una funci´on on y da por resultado otra funci´on. on. Este Este operador operador se puede componer para obtener derivadas de cualquier orden. A saber, Dk x = D(D(D(
· · · (Dx Dx)))) )))) = x(k)
Otro operador que a una funci´on on le hace corresponder otra es la multiplicaci´ on o n por una funci´ on on (po ( podr dr´´ıa ser s er una constante). consta nte). Y todav to dav´´ıa otro o tro operador opera dor es es aDk x := ax(k) donde a es una funci´on on continua (posiblemente constante). Dos operadores pueden sumarse para dar otro operador: Si O1 y O2 son dos operadores (O1 + O2 )x := O1 x + O2 x. Es decir, el resultado es la suma de los resultados. Tambi´ en en pueden componerse para dar como resultado otro operador, a saber sab er O1 O2 x = O1 (O2 x). De este modo, la ecuaci´on on (5.1) puede verse de la siguiente manera: Lx := (Dn + an−1 Dn−1 +
· · · + a1D + a0) x = 0.
Est´ a claro de esta expresi´on on lo que entend´ entend´ıamos cuando dec´ dec´ıamos que el polinomio caracter´ ter´ıstico de la ecuaci´on on se obtiene al reemplazar derivaciones por potencias de λ. Si λ1 , . . . , λk son las raices (reales y complejas) del polinomio p olinomio caracter´ caracter´ıstico con multiplicidades n1 , . . . , nk respectivamente (con lo cual n1 + + nk = n), se sigue que p( p(λ) = (λ
···
− λ1)n (λ − λ2)n · · · (λ − λk )n . 1
2
k
La misma factorizaci´ factorizacion o´n es v´alida alida para el operador L a saber Lx = (D Adem´ as as los operadores (D (D (D
− λ1)n (D − λ2)n · · · (D − λk )n x. 1
− λ j )n
j
2
k
conmutan, conmutan, es decir
− λ j )n (D − λi)n x = (D − λi)n (D − λ j )n x. j
i
i
j
− λ j )n
Por lo tanto si sabemos encontrar una base de soluciones de cada operador (D (D tendremos n soluciones de la ecuaci´ on on (5.1).
− λ1)n x = 0, tendremos · · · (D − λk )n x = (D − λ2)n · · · (D − λk )n
En efecto, si por ejemplo, (D (D Lx = (D
− λ1)n (D − λ2)n 1
2
j
1
k
2
k
(D
− λ1)n x 1
= 0.
Adem´ as as las n j soluciones correspondientes a la raiz λ j son linealmente independientes y se puede ver que las n son linealmente independientes. Veamos entonces c´omo omo es la soluci´on on general de una ecuaci´on on cuando cuando el operador operador diferencial diferencial n n asociado es (D (D λ) , es decir, cuando el polinomio caracter´ caracter´ıstico es p(r) = (r λ) .
−
−
50
5. ECUACIONES LINEALES DE ORDEN n CON COEFICIENTES CONSTANTES
Para ´esto esto observemos que (D Iterando, (D
− λ)(e )(eλt y ) = D(eλt y ) − λeλt y = λeλt y + eλt y − λeλt y = eλt y .
− λ)2(eλty) = (D − λ) (D − λ)(e )(eλt y )
De modo que iterando llegamos a (D
= (D
− λ) eλty
= eλt y .
− λ)n(eλty) = eλty(n)
Como cualquier funci´on on x puede escribirse como x = eλt y (basta tomar y = e−λt x), se ve que (D λ)n x = 0 si y solo o´lo si x = eλt y con y(n) = 0. Es decir, si y s´olo olo si x = eλt y con y = c1 + c2 t + + cn−1 tn−1 . En efecto,
−
···
0 = y(n) = D(y(n−1) ) si y solo o´lo si y(n−1) = kn para una constante kn . Integando obtenemos y(n−2) = kn t + kn−1 y finalmente y= y se tiene lo afirmado.
kn (n
− 1)!
y (n−3) =
,
tn−1 +
kn 2 t + kn−1 t + kn−2 2
kn−1 n−2 t + (n 2)!
−
Por lo tanto la soluci´ on on general de (D (D polinomio de grado a lo sumo n 1.
−
,...
· · · + k1t + k0
− λ)nx = 0 es x(t) = eλt pn(t) donde pn(t) es un
Volviendo al operador general con polinomio caracter´ caracter´ıstico p( p(λ) = (λ
− λ1)n (λ − λ2)n · · · (λ − λk )n , 1
2
k
vemos que cualquier funci´on on de la forma
x(t) = pn (t)eλ t +
(5.3)
1
1
con pnj polinomio de grado a lo sumo n j
· · · + pn (t)eλ t k
k
− 1 es soluci´on on de la ecuaci´on on (5.1).
Por lo tanto hemos encontrado n soluciones de (5.1), (5.4)
{eλ t, teλ t, t2eλ t, . . . , tn −1eλ t , . . . , eλ t, teλ t , . . . , tn −1eλ t}, 1
1
1
1
1
k
k
k
k
que se puede ver que son linealmente independientes. De modo que (5.3) da la soluci´ on on general general de (5.1) y (5.4) es una base de soluciones. Con ´esto encontramos encontram os la soluci´ solucion o´n general compleja. compleja. La soluci´ on general real se obtiene reemon plazando en la base (5.4) los pares de soluciones conjugadas por parte real y parte imaginaria. Por ejemplo, si λ j = α j + iβ j con β j = 0, tomamos en lugar de las 2n 2n j soluciones complejas
¯
¯
¯
¯
eλj t , teλj t , t2 eλj t , . . . , tnj −1 eλj t , eλj t , teλj t , t2 eλj t , . . . , tnj −1 eλj t , las 2n 2n j soluciones reales eαj t cos β j t, eαj t sen β j t,teαj t cos β j t,teαj t sen β j t , . . . , tnj −1 eαj t cos β j t, tnj −1 eαj t sen β j t.
5. ECUACIONE ECUACIONES S LINEALES LINEALES DE ORDEN n CON COEFICIENTES CONSTANTES Ejemplo
51
5.3. Hallar las soluciones de la ecuaci´on on x(5)
− x(4) + 2x 2x − 2x + x − x = 0.
El polinomio caracter´ caracter´ıstico de la ecuaci´ on on es p( p(λ) = λ5
− λ4 + 2λ 2λ3 − 2λ2 + λ − 1 = (λ − 1)(λ 1)(λ4 + 2λ 2λ2 + 1) = (λ (λ − 1)(λ 1)(λ2 + 1)2 , cuyas cu yas ra´ıces ıc es son λ1 = 1, con multiplicidad 1 y λ2 = i, λ3 = −i con multiplicidad 2. Por lo tanto una base de soluciones es {et , cos t, sen t, t cos t, t sen t}.
y la soluci´ on on general es
x(t) = c1 et + (c (c2 + c3 t)cos t + (c (c4 + c5 t)sen t.
Utilicemos Utilicemos ahora el m´ etodo etodo de variaci´ ariacion o´n de par´ametro ametross para para hallar hallar la soluci soluci´ on o´n de una ecuaci´on on lineal linea l no homog´enea. enea. Ejemplo
5.4. Hallar las soluciones de la ecuaci´on on x
− 2x + x = t.
Hallamos primero una base de soluciones de la ecuaci´on on homog´ ho mog´enea enea asoci a sociada ada cuyo c uyo polin´ po lin´ omio omio caract car acter er´´ıstico ıst ico es p( p(λ) = λ2 2λ + 1 = (λ (λ 1)2 . Por lo tanto la unica u ´ni ca ra´ız ız es λ = 1 y una base de soluciones es
−
−
{et, tet}. Buscamos una soluci´ on particular de la ecuaci´ on on on no n o homog´ h omog´enea enea de la l a forma f orma x(t) = c1 (t)et + c2 (t)tet donde las derivadas de las funciones c1 , c2 satisfacen el sistema c1 et + c2 tet = 0, c1 et + c2 (et + tet ) = t. De donde, c2 = te−t y c1 = Integrando obtenemos c2 = y
c1 =
te
−t
−
−c2 t = −t2e−t. dt = 2 −t
t e
−t
−te
+
2 −t
dt = t e
De modo que la soluci´on on general es x(t) = (t2 + 2t 2t + 2)
e−t dt =
− 2
1)e−t −te−t − e−t = −(t + 1)e
te−t dt = t2 e−t + 2(t 2(t + 1)e 1)e−t .
− (t + 1)t 1)t + c1 et + c2 tet = t + 2 + c1 et + c2 tet .
52
5. ECUACIONES LINEALES DE ORDEN n CON COEFICIENTES CONSTANTES
Ejercicios
(1) (a) Encontrar Encontrar (i) y (ii) y (iii) y
una base base de soluciones soluciones reales de las siguien siguientes tes ecuaciones ecuaciones:: 8y + 16y 16y = 0 2y + 10y 10y = 0 y 2y = 0
− − − −
(b) En cada uno de los casos anteriores encontrar una soluci´ on on exacta de la ecuaci´on on x x − no homog´enea enea correspondiente corresp ondiente con t´ermino ermino independiente indep endiente x, e , 1 y e . a (2) Sean (a (a1 , b1 ) y (a2 , b2 ) dos puntos del plano tales que a − no es un n´ umero umero entero. π (a) Probar que existe exactamente una soluci´on on de la ecuaci´ on on diferencial y + y = 0 cuya gr´ g r´afica afica pasa por esos puntos. (b) ¿Se cumple en alg´ un un caso la parte (a) si a1 a2 es un m´ ultiplo ultiplo entero de π ? (c) Generalizar el resultado de (a) para la ecuaci´on on y + k2 y = 0. Discutir Discut ir tambi´en en el caso k = 0. 1
2
−
(3) Hallar todas las soluciones de y
− y − 2y = 0 y de y − y − 2y = e−x que verifiquen:
(a) y(0) = 0, 0, y (0) = 1
(b) y(0) = 1, 1, y (0) = 0
(c) y (0) = 0, 0, y (0) = 0
(d) limx→+∞ y(x) = 0
(e) y (0) = 1
(f ) y (0) = 1
(4) En el interior de la Tierra, la fuerza fuerza de gravedad gravedad es proporcional proporcional a la distancia al centro. centro. Si se perfora un orificio que atraviese la Tierra pasando por el centro, y se deja caer una piedra en el orificio, ¿con qu´ e velocidad llegar´ llega r´ a al centro?. (5) La ecuaci´on on x2 y + pxy + qy = 0 ( p,q constantes) se denomina ecuaci´ on on de Euler. t (a) Demuestre que el cambio de variables x = e transforma la ecuaci´on on en una con coeficientes constantes. (b) Aplique (a) para resolver en R>0 las ecuaciones: i. x2 y + 2xy 2xy 6y = 0 ii. x2 y xy + y = 2x
−
−
(6) Vibraciones en sistemas mec´ anicos: anicos: Una carreta de masa M est´ a sujeta sujeta a una pared por medio de un resorte, que no ejerce fuerza cuando la carreta est´a en la posici´ on on de equilibrio x = 0. Si la carreta se desplaza a una distancia x, el resorte ejerce una fuerza de restauraci´on on igual a κx, κx, donde κ es una constante constante positiva positiva que mide la rigidez del resorte. Por la segunda ley del movimiento de Newton, se tiene que:
−
d2 x (1) M 2 = κx o bien x + a2 x = 0, a = κ/M κ/M dt (a) Si la carreta se lleva a la posici´on on x = x0 y se libera sin velocidad inicial en el instante t = 0, hallar la funci´on on x(t). Verifica erificarr que se trata trata de una funci´ funcion o´n peri´ odica. odica. Calcular Calcul ar su per´ıodo ıodo τ , τ , y su frecuencia f = 1/τ (la cantidad cantidad de ciclos por unidad de tiempo). Verificar que la frecuencia de vibraci´on on aumenta al aumentar
−
EJERCICIOS
53
la rigidez del resorte, o al reducir la masa de la carreta ( como dice el sentido com´ un) un) y que la amplitud de esta oscilaci´on on es x0 . Si se produce produce una amo amorti rtigua guaci´ ci´ o n que se opone al movimiento, y de magnitud on proporcional a la velocidad (= c dx on on (1) que dt ) debida al rozamiento, la ecuaci´ describe el movimiento de la carreta en funci´on on del tiempo se convierte en:
−
M
d2 x dx + c + κx = 0 dt2 dt
o bien: d2 x dx c κ 2 (2) + 2b 2 b + a x = 0 b = a = dt2 dt 2M M (b) Si b > a (la fuerza fuerza de fricci´ fricci´on on debida debida al rozaminto rozaminto es grande grande en comparaci´ comparaci´ on con la rigidez del resorte), encontrar la soluci´ on de (2) que verifique como antes on x(0) = x0 , x (0) (0) = 0. Prob Probar ar que no hay hay ning ningun unaa vibr vibrac aci´ i´ on o n y que la carreta vuelve simplemente a su posici´ on on de equilib equilibrio rio.. Se dice que el mo movim vimien iento to est´ a sobreamortiguado. (c) Si b = a, ver que tampoco hay vibraci´ on y que el comportamiento es similar al del on caso anterior. Se dice que el movimiento es cr´ıticamente ıticamente amortiguado. (d) Si ahora b < a (caso (caso subamo subamortig rtiguad uado), o), probar probar que la soluci´ soluci´ on de (2) con las condiciones iniciales x(0) = x0 , x (0) = 0 es:
√α2 + b2
e−bt cos(αt cos(αt θ) α donde α = a2 b2 , y tan θ = b/α. b/α. Esta funci´on on oscila con una amplitud que se reduce exponencialmente. Su gr´afica afica cruza la posici´ on on de equilibrio x = 0 a intervalos regulares, reg ulares, aunque no es peri´odica. odica. Hacer Hac er un dibujo dibujo.. Probar Probar que el tiempo tiempo requer requerido ido para volve volverr a la posici´ on o n de equilibrio es: 2π T = x(t) = x0
−
√ −
− κ M
c2 4M 2
y su “frecuencia” est´a dada por f = 1/T llamada frecuencia natural del sistema. Notar que esta frecuencia disminuye al disminuir la constante de amortiguaci´on on c. Hasta ahora hemos considerado vibraciones libres, porque s´olo olo act´ uan uan fuerzas internas al sistema. Si una fuerza F ( F (t) act´ ua sobre la carreta, la ecuaci´on ua on ser´a: a: d2 x dx +c + κx = F ( F (t) 2 dt dt (e) Si esta fuerza es peri´odica odica de la forma F ( F (t) = F 0 cos ωt, ωt, con F 0 , ω constantes, hallar x(t). Al valor ω/2 ω/2π se lo llama frecuen frecuencia cia impresa impresa al sistema. sistema. ωc Si tan φ = κ−ω M , probar que la soluci´on on general de (3), con F ( F (t) = F 0 cos ωt puede escribirse: F 0 x(t) = e−bt (C 1 cos(αt cos(αt)) + C 2 sen(αt sen(αt)) )) + cos(ωt cos(ωt φ) (κ ω 2 M )2 + ω 2 c2 (3)
M
2
−
→ ∞
−
El primer t´ermino ermino tiende a cero para t + , luego es “transitorio”, es decir, a medida que pasa el tiempo, la soluci´on on se parece m´as a s y m´ as as al segundo sumando.
54
5. ECUACIONES LINEALES DE ORDEN n CON COEFICIENTES CONSTANTES
Notar que la frecuencia de esta funci´on on es la frecuencia frecuencia impresa al sistema, sistema, y que la F 0 amplitud es el coeficiente . ¿Qu´e pasa cuando la frecuencia fre cuencia (κ ω 2 M ) M )2 + ω 2 c2 impresa impresa se acerca acerca a la frecuencia frecuencia natural del sistema? (Este fen´ omeno se conoce con el nombre de resonancia). (f) Si b < a (caso subamortiguado) hallar la frecuencia impresa ω que provoca amplitud m´ axima. ¿Siempre existe este valor? Este valor de frecuencia impresa (cuando axima. existe) se denomina frecuencia de resonancia. Demostrar que la frecuencia de resonancia es siempre menor que la frecuencia natural.
−
(7) Hallar la soluci´ on general de las siguientes ecuaciones, empleando la soluci´on on on dada: (a) (b) (c) (d)
xy + 2y 2y + xy = 0 xy y 4x3 y = 0 xy y 4x3 y = 0 (1 x2 )y 2xy + 2y 2y = 0
− − − − − −
I = R>0 I = R>0 I = R<0 I = ( , 1), 1), ( 1, 1), 1), (1, (1,
y1 (x) = senx x y1 (x) = exp(x exp(x2 ) y1 (x) = exp(x exp(x2 ) ) y1 (x) = x
−∞ − − ∞ Este ultimo u ´ ltimo es un caso especial especial de la ecuaci´ on on (1 − x2 )y − 2xy + p( p( p+ p + 1)y 1)y = 0 (ecuaci´ on on de Legendre), correspondiente al caso p = 1, en los intevalos en que la ecuaci´on o n es normal.
∀ ∈
(8) Sabiendo que y1 (x) = 1, x R es soluci´ soluci´ on on de la ecuaci´on on homog´enea enea asociada, aso ciada, hallar halla r todas las soluciones de y + xy = 3x. (9) Probar que las funciones φ1 (t) =
t2 t 0 t
≤0 ≥0
φ2 (t) =
0 t t2 t
≤0 ≥0
son linealmente independientes en R pero que W ( W (φ1 , φ2 )(0) = 0. ¿Existe alg´ un un sistema lineal normal de orden 2 definido en alg´un un intervalo ( , ) que admita a φ1 , φ2 como base de soluciones?
−
{
}
CAP´ıTULO ıT ULO 6
Comportamiento asint´ otico de las soluciones otico En los cap´ cap´ıtulos previos, hemos visto que bajo condiciones condiciones muy generales, generales, una ecuaci´ ecuaci´ on on diferencial admite soluci´ on on unica u ´ nica y que dicha soluci´on on es continua con respecto a las condiciones iniciales. Por otro lado, hemos estudiado varias situaciones en donde la soluci´on on puede ser calculada de manera expl´ expl´ıcita. ıcita . Sin embargo, embarg o, en la l a mayor´ mayor´ıa de los l os casos caso s (por (po r ejemplo ejempl o si estamos es tamos estudiando estudia ndo un sistema de ecuaciones no lineal) esa soluci´ on, que existe, no puede ser calculada. De todas formas, on, es mucha la informaci´on on que uno puede llegar a dar sobre la soluci´on on (sobre el comportamiento de la soluci´ on) on) a´ un un sin conocer la f´ormula ormula de la misma y ese es el ob jetivo de este cap´ cap´ıtulo. Por mayor simplicidad, simplicidad, y dado que en la may mayor or´´ıa de las aplicaciones aplicaciones sucede, de ahora en m´ as vamos a suponer que el sistema de ecuaciones diferenciales bajo consideraci´on as on es aut´ onomo, onomo, es decir, consideraremos ecuaciones de la forma X = F ( F (X )
(6.1) donde X : I
⊂ R → Rn y F : Ω ⊂ Rn → Rn es un campo C 1.
En esta situaci´ on, on, es decir si la ecuaci´ on on es aut´ onoma, se tiene la siguiente propiedad que onoma, nos ser´a de d e gran utilidad. utili dad. ´n Proposicion o
6.1. Sean X 1 : I 1 de (6.1). (6.1). Entonces se tiene
⊂ R → Rn y X 2 : I 2 ⊂ R → Rn dos soluciones maximales
{X 1(t) | t ∈ I 1} ∩ {X 2(t) | t ∈ I 2} = ∅
o
{X 1(t) | t ∈ I 1} = {X 2(t) | t ∈ I 2}.
En otras palabras, lo que dice la proposici´ on 6.1 es que dos trayectorias dadas, o bien son on id´enticas, enticas, o bien bi en no n o se cruzan. ´ n. Demostracion. o
Este hecho hecho es una consecuencia consecuencia de la unicidad unicidad de soluciones. soluciones. Supongamos Supongamos que existen existen t1 I 1 y t2 I 2 tales que X 1 (t1 ) = X 2 (t2 ) = X 0 . Entonces, si definimos la funci´ on on ˜ ˜ X (t) = X 2 (t t1 + t2 ), las funciones X 1 y X son soluciones de
∈ −
∈
X = F ( F (X ), X (t1 ) = X 0 ,
de donde, por la unicidad, concluimos que X 1
≡ X ˜ .
˜ ). Ahora la proposici´on on queda demostrada, observando que I m(X 2 ) = I m(X ).
Otra propiedad importante de los sistemas aut´onomos, onomos, es que uno puede suponer siempre que las condic condicion iones es inicial iniciales es est´ estan a´n dadas en el instante t0 = 0. En efe efect cto, o, si X (t) es una 55
´ 6. COMPORTAMIENTO ASINTOTICO DE LAS SOLUCIONES
56
soluci´ on de (6.1), razonando de manera similar a en la prueba de la Proposici´ on on on 6.1, definimos ¯ X (t) = X (t + t0 ). Luego si X verifica
˜ verifica tenemos que X
X = F ( F (X ), X (t0 ) = X 0 , ˜ = F ( ˜ ), X F (X ˜ (0) X (0) = X 0 .
1. Diagramas Diagramas de fases fases
Una forma muy habitual y gr´ afica de entender el comportamiento asint´ afica otico otico o din´ amico de las soluciones de una ecuaci´ on on de la forma (6.1) es a trav´ es es del llamado diagrama de fases. fases. Seg´ un un la Proposici´ on 6.1, si tomamos dos soluciones distintas de la ecuaci´ on on on diferencia diferencial, l, eso n determina trayectorias disjuntas en el “plano de fases” R . Luego, el diagrama de fases consiste precisamente en graficar en Rn una colecci´ colecci´ on de tales trayectorias y las mismas nos dar´an on an una idea genera generall del comportam comportamien iento to de todas las traye trayecto ctoria riass del sistema. sistema. Esto Esto lo hacemo hacemoss fundamentalmente en el caso en que n = 2 ya que es el caso en el que podemos graficar bien. La idea de esbozar el diagrama de fases es el poder predecir el comportamiento asint´otico otico (para el tiempo tendiendo a infinito) de las soluciones dependiendo de donde se encuentran inicialmente. Dibujando “suficientes” trayectorias deber´ deber´ıa ser posible saber sab er si las soluciones tienden a estabilizarse en un cierto punto, si oscilan (soluciones peri´ odicas), si se vuelven infinitamente odicas), grandes, etc... A veces se puede tener una idea de c´omo omo ser´ a el diagrama de fases dibujando en “muchos” puntos del plano la flecha que corresponde al campo vectorial F ( F (x, y ) = f 1 (x, y), f 2 (x, y) . Recordemos que las curvas soluci´on on del sistema X = F ( F (X ) son las trayectorias del campo F , F , es decir, son curvas x = x(t)
y = y(t) tales que el vector velocidad – tangente a la curva – en el punto (x ( x0 , y0 ), con x0 = x(t0 ), y0 = y(t0 ), es el vector f 1 (x0 , y0 ), f 2 (x0 , y0 ) .
Por lo tanto, tanto, si dibujamos muchos muchos de estos vectores vectores podremos tratar de adivinar adivinar c´ omo omo ser´ an an las trayectorias buscando curvas que al pasar por un punto pasen con tangente igual a la flecha que dibujamos en ese punto. Veamos un ejemplo. Ejemplo
6.1. Consideremos el siguiente modelo de dos poblaciones simbi´ oticas, oticas,
− y˙ = (−γ + γ + δx) δx) y
x˙ = ( α + βy) βy ) x
con α, β,γ,δ > 0, que se obtiene al considerar dos poblaciones x e y que decrecen con raz´on on de crecimiento constante en ausencia de la otra especie y cuya raz´on on de crecimiento crece en forma proporcional a la otra poblaci´ on. on.
1. DIAGRAMAS DE FASES
57
Cuando una trayectoria pasa por un punto con y = α/β lo hace en forma perpendicular perpendicular a esa recta porque su vector velocidad (0, (0, ( γ + γ + δx) δx ) y) es vertical. Cuando cruza la recta x = γ/δ lo hace en forma horizontal porque su vector velocidad es (( α + βy) βy ) x, 0). Las flechas que indican el campo F ( F (x, y) = ( α + βy) βy ) x, ( γ + γ + δx) δx ) y en cada c ada punto son aproximadamente aproxima damente as´ as´ı.
−
−
−
−
α / β
γ / δ
Cuando estamos cerca del (0, (0, 0) o del (γ/δ, (γ/δ, α/β α/β ) las flechas son muy cortas. Esto indica que si estamos cerca de estos puntos pasaremos con velocidad muy chica y nos quedaremos cerca por mucho tiempo. El caso l´ımite lo tenemos tenemos en estos puntos puntos en los que no hay flecha. flecha. Lo que sucede es que el campo F se anula anul a ah a h´ı. ¿Qu´e suced s ucedee si s i en e n alg a lg´ un u´n instante t0 estamos en un punto (x0 , y0 ) donde F = (f 1 , f 2 ) se anula? anula? Para contestar contestar esta pregunta pregunta observemos observemos que la funci´ funcion o´n x(t) x0 , y(t) y0 es soluci´ on on de
≡
≡
x˙ = f 1 (x, y),
x(t0 ) = x0
y˙ = f 2 (x, y ),
y(t0 ) = y0
Por unicidad de soluci´ solucion o´n tendremos que ´esta esta es la unica u´nica soluci´ on que pasa por este punto. Pero on esta es una soluci´ on constante a la que por lo tanto llamamos Soluci´ on on estacionaria . Es f´ facil a´cil ver que las ´unicas unicas soluciones estacionarias son las soluciones constantes igual a un valor (x (x0 , y0 ) tal que F ( F (x0 , y0 ) = 0. Esto se debe a que una soluci´ on on constante tendr´a derivada nula en todo tiempo y por lo tanto el campo F , F , que es igual a la derivada, debe anularse. A estos puntos los llamamos puntos de equilibrio o puntos cr´ıticos. ıticos. Tenemos entonces unos puntos especiales en el plano de fases, los ceros del campo F que corresponden a soluciones estacionarias, y vemos que cerca de esos puntos las trayectorias pasan muy muy despacio despacio.. Pero, Pero, ¿qu´ ¿qu´e es lo que hacen? hacen? ¿Se acerca acercan? n? ¿Se alejan? alejan? ¿Se quedan quedan cerca cerca sin acercarse? Puede pasar cualquiera de estas cosas y es algo que uno trata de observar en el diagrama de fases. Una observaci´on on importante es que s´olo olo se llega a un equilibrio en tiempo infinito (+ o infinito) como demostraremos en el Lema 6.1 al final de esta secci´on. on. Otra cosa que observamos en el diagrama de flechas es que sobre los ejes coordenados, las flechas apuntan en la direcci´on on del eje. Esto refleja refleja el hecho hecho de que si inicialmente inicialmente estamos estamos sobre uno de los ejes permaneceremos permaneceremos ah´ ah´ı por todo tiempo. En efecto, efecto, supongamos supongamos que inicialmente inicialmente
58
´ 6. COMPORTAMIENTO ASINTOTICO DE LAS SOLUCIONES
y0 = 0. Sea x(t) la soluci´ on on de la ecuaci´on on x˙ =
−α x
con x(0) = x0 . Entonces X (t) = (x(t), 0) es la soluci´ on del sistema que inicialmente vale (x on (x0 , 0). Y an´alogamente alogamente con datos iniciales sobre el otro eje x = 0. En particular, los semiejes corresponden a trayectorias y como las trayectorias no se cortan se sigue que las trayectorias que se inician en el primer cuadrante no salen de ´el. el. Esto dice que si inicialmente las dos poblaciones son positivas, lo ser´ an para todo tiempo. Que x(t) e y (t) no se an vuelvan vuelvan negativas negativas es importante importante para que el sistema sistema refleje refleje la realidad realidad (son poblaciones). poblaciones). Pero Pero lo que acabamos de observar es que ninguna poblaci´on on se extingue en tiempo finito. Volvamos al ejemplo y tratemos de dibujar trayectorias correspondientes a este campo.
α / β
γ / δ
Da la impresi´on on de que las trayectorias que comienzan cerca de (0, (0 , 0) tienden a (0, (0, 0) cuando el tiempo tiende a + mient mientras ras que cerca cerca del otro punto punto de equilib equilibrio rio,, el (γ/δ, α/β ), ) , hay trayectorias que se acercan pero p ero la mayor´ mayor´ıa parece alejarse. a lejarse.
∞
Conjeturamos que el diagrama de fases para este sistema es
α / β
γ / δ
1. DIAGRAMAS DE FASES
59
Pero, Pero, ¿c´ omo omo estar seguros seguros de que es as´ as´ı? En este curso curso lo que vamos vamos a poder asegurar asegurar es c´omo omo es el diagrama diagrama de fases fases cerca cerca de los punto punto de equilib equilibrio rio.. El diagrama diagrama completo completo s´ olo lo podremos podremos conjetur conjeturar. ar. La raz´ on on es que cerca de la may mayor or´´ıa de los equilibrios equilibrios los sistemas sistemas no lineales tienen diagramas de fase muy parecidos a los de un sistema lineal con coeficientes consta constant ntes. es. En este caso, como conocemos conocemos las f´ ormulas que nos dan las soluciones, podemos ormulas saber exactamente c´omo omo es el diagrama de fases. Veremos ´esto esto en detalle en la pr´ oxima oxima secci´ on. on. Pero antes veamos qu´ e tipo de informaci´ on podemos inferir del diagrama de fases. on Observ Observamos amos que hay dos equilibr equilibrios ios.. Uno es muy muy claro, claro, si inicia inicialme lment ntee no hay ning´ ning´ un un miembro en ninguna de las dos poblaciones, p oblaciones, ´esto seguir´ a asi por siempre. M´ as as a´ un, un, del diagrama de fases vemos que si alguna de las poblaciones es “chica”, ambas poblaciones desaparecer´an an (en realidad no lo har´ an an porque s´olo olo se vuelven nulas cuando el tiempo se vuelve infinito, pero eventualmente ambas poblaciones ser´an an extremadamente chicas.) Hay otro equilibrio equilibrio que es que la poblaci´ on on x sea γ/δ y la poblaci´ on on y sea α/β . Esta sta eess una situaci´ on on muy inestable. inestable. Vemos que s´ olo olo si hay una relaci´ on muy particular entre ambas on poblaciones poblaciones se acercar´ acercar´ an al equilibrio. En cualquier otro caso, por m´as an as cerca que se encuentren del equilibrio se alejar´an an a medida que pase en tiempo. M´as a s a´ un, hay poblaciones tan cercanas un, al equilibrio como se quiera que tender´an an a desaparecer y otras que crecer´an an sin l´ımite. ımi te. Demostremos ahora un lema que dice que s´olo olo se llega a un punto de equilibrio (o punto cr´ cr´ıtico) ıtico ) en tiempo tiemp o infinito. infinit o. Lema
6.1.
→
→ ∈
≡
(1) Si X (t) X 0 cuando t t0 R y X (t) X 0 , se sigue que F ( F (X 0 ) = 0. Es decir decir,, si una trayectoria trayectoria tiende a un punto cr´ cr´ıtico lo har´ a en tiempo infinito (+ o - infinito). (2) Si X (t) X 0 cuando t + se sigue que F ( F (X 0 ) = 0. An´ alogamente con t . Es decir, si una trayectoria trayectoria tiene un l´ımite para para tiempo tendiendo a infinito, ese l´ımite debe ser un punto cr´ cr´ıtico.
→
→ ∞
→ −∞
Demostremos (1) Sabemos que si F ( F (X 0 ) = 0, la unica u ´ni ca soluci´ sol uci´on on de X = F ( F (X ) que satisface X (t0 ) = X 0 es la funci´on on id´entic ent icam ament entee X 0 . Como X (t0 ) = X 0 , se siguir´ sig uir´ıa ıa que X (t) X 0 lo que hab´ hab´ıamos supuesto que no pasaba. Por lo tanto, F ( F (X 0 ) = 0. 0. ´ n. Demostracion. o
≡
Ahora demostremos (2). Tenemos X (t + 1)
−
x1 (t + 1) X (t) = x2 (t + 1)
x1 (σ1 ) x2 (σ2 )
− x1(t) − x2(t)
=
=
f 1 x1 (σ1 ), x2 (σ1 ) f 2 x1 (σ2 ), x2 (σ2 )
donde t < σi < t + 1, i = 1, 2 y el campo F tiene componentes (f (f 1 , f 2 ).
Como X (t) = F ( F (X (t)) F ( F (X 0 ) cuando t + se sigue que X (t + 1) X (t) cuando t + . Pero, X (t + 1) X (t) X 0 X 0 = 0. Por lo tanto, F ( F (X 0 ) = 0.
→ ∞
→
−
→
→ ∞ −
−
→ F ( F (X 0 )
´ 6. COMPORTAMIENTO ASINTOTICO DE LAS SOLUCIONES
60
2. Diagramas Diagramas de fases de sistemas sistemas lineales lineales a coeficiente coeficientess constan constantes tes
Consideremos el sistema X = AX con A
∈ R2×2 y X =
x1 . x2
Queremos dibujar aproximadamente las trayectorias del sistema dependiendo de c´omo omo son los autovalores λ1 y λ2 de A. Supongamos que 0 no es autovalor. En este caso (0, (0 , 0) es el unico u ´ nico equilibrio equilibrio del sistema. sistema. Caso I. λ1 > 0 > λ 2 . Sean ξ1 y ξ2 los autovectores correspondientes a los λ1 y λ2 respectiv respectivamen amente. te. Entonces, Entonces, la soluci´ on on general es X (t) = c1 eλ t ξ1 + c2 eλ t ξ2 . 1
2
Observemos que si X (0) ( 0) = cξ1 , es decir, si inicialmente estoy en la recta de autovectores asociados a λ1 , se sigue que c1 = c y c2 = 0. Por Por lo tant tanto, o, X (t) = ceλ t ξ1 y en todo instante premanezco en esa recta. Adem´ as, as, como λ1 > 0, se tiene que X (t) cuando t + y X (t) 0 cuando t . 1
→
→∞
→ −∞
→ ∞
An´ alogamente, si inicialmente estoy en la recta de autovectores asociados a λ2 tengo X (0) alogamente, (0) = λ t cξ2 . Por lo tanto, c1 = 0 y c2 = c y tengo X (t) = ce ξ2 . De donde en todo instante premanezco premanezco en esa recta. Adem´as, as, como λ2 < 0, se tiene que X (t) cuando t y X (t) 0 cuando t + . 2
→∞
→ ∞
→ −∞
→
Si inicialmente no estoy en ninguna de las dos rectas de autovectores, tanto c1 como c2 son no nulos. Llamemos y1 (t) e y2 (t) a los coeficientes del vector X (t) en la base ξ1 , ξ2 . Es decir, escribamos X (t) = y1 (t)ξ1 + y2 (t)ξ2 .
{
}
Por la forma que tiene la soluci´on on general vemos que y1 (t) = c1 eλ t e y2 (t) = c2 eλ t donde c1 y c2 con las componentes del vector X (0) (0) en la base ξ1 , ξ2 . Es deci decir, r, X (0) (0) = c1 ξ1 + c2 ξ2 . Tenemos 1
{
y1 (t) y2 (t)
→ +∞, → 0,
(t (t
→ +∞) → +∞)
y1 (t) y2 (t)
2
}
→ 0, → +∞,
(t (t
→ −∞), → −∞).
Por lo tanto, X (t) se acerca m´as a s y m´as as a la recta generada por ξ1 cuando t recta generada por ξ2 cuando t .
→ −∞
→ +∞ y a la
Este an´ alisis alisis ya permitir p ermitir´´ıa intentar esbozar el diagrama de fases. Sin embargo ser´ a m´ as as sencillo hacer lo siguiente. Dibujemos primero las curvas (y (y1 (t), y2 (t)) que dicen c´omo omo cambian los coeficientes de la soluci´ on on en la base ξ1 , ξ2 . Una vez que tenemos claro este dibujo, observamos que x1 y X = =Q 1 x2 y2 donde Q = [ξ1 ξ2 ] es la matriz cuyas columnas son los vectores ξ1 y ξ2 .
{
}
Por lo tanto, el diagrama de fases en el plano (x ( x1 , x2 ) se obtiene del dibujo de las curvas (y1 (t), y2 (t)) en el plano (y (y1 , y2 ) mediante mediante una transformac transformaci´ i´ on on lineal (de matriz Q).
2. DIAG DIAGRA RAMA MAS S DE FASES ASES DE SIST SISTEM EMAS AS LINE LINEAL ALES ES A COEF COEFIC ICIE IENT NTES ES CONS CONST TANTE ANTES S
61
Sigamos entonces estas ideas y dibujemos primero las curvas (y (y1 (t), y2 (t)). Observemos que λ t λ t como y1 = c1 e , y2 (t) = c2 e se sigue que 1
2
y1 c1
e = y y2 (t) = c2 Sea α =
λ2 < 0. Entonces y2 = k y1 λ1
1/λ1
t
y1 c1
λ2 /λ1
| |λ /λ
= k y1
2
1
| |α y el gr´afico afico de las curvas (y (y1 (t), y2 (t)) es y
2
y
1
y el diagrama de fases en el plano (x (x1 , x2 ) es c x
ξ
1
2
c
ξ
x
1
2
´ 6. COMPORTAMIENTO ASINTOTICO DE LAS SOLUCIONES
62
Caso II. λ1 > λ 2 > 0 Como en el caso anterior tenemos X (t) = y1 (t)ξ1 + y2 (t)ξ2 donde ξ1 y ξ2 son autovectores correspondientes a los autovalores λ1 y λ2 respectivamente. Como antes se tiene
| |α
y2 = k y1
0<α=
λ2 <1 λ1
Por lo tanto las curvas (y (y1 (t), y2 (t)) son y
2
y
1
y el e l diagrama diagra ma de fases ser´a x 2
c
ξ1
x 1
c
t
ξ2
Observemos que en este caso, dado el signo de λ1 y λ2 se tiene que y1 (t), y2 (t) + y y1 (t), y2 (t) 0 cuando t .
→ ∞
→
→ −∞
→ +∞ cuando
2. DIAG DIAGRA RAMA MAS S DE FASES ASES DE SIST SISTEM EMAS AS LINE LINEAL ALES ES A COEF COEFIC ICIE IENT NTES ES CONS CONST TANTE ANTES S
63
Caso Caso III. λ1 < λ2 < 0 Este caso es exactamente como el anterior. Se tiene
| |α
y2 = k y1
λ2 <1 λ1
0<α=
s´olo olo que en este caso las flechas se invierten y los diagramas de curvas en el plano ( y1 , y2 ) y en el plano de fases es igual al Caso II con las flechas invertidas.
Caso IV. λ1 = λ2 = λ = 0. Supongamos que A = λI . El caso A = λI lo dejamos como ejercicio. La soluci´ solucion o´n general es X (t) = c1 eλt ξ1 + c2 eλt (ξ1 t + ξ2 ) donde ξ1 es un autovector correspondiente al autovalor λ. Por lo tanto, X (t) = y1 (t)ξ1 + y2 (t)ξ2 con
y1 (t) = (c1 + c2 t)eλt
y2 (t) = c2 eλt .
Observemos que y1 (0) = c1 , y2 (0) = c2 . Por lo tanto, si X (0) (0) = cξ1 se sigue que c1 = c y c2 = 0. λt Por lo tanto, X (t) = c1 e ξ1 y se ve que la recta generada por ξ1 (la recta de autovectores asociados al unico u´nico autovalor λ) es invari invarian ante. te. En este caso, caso, la recta genera generada da por ξ2 no es invariante. Tenemos
e
λt
y2 = c2
1 y2 t = log . λ c2
⇒
Por lo tanto, y2 c2 y2 y1 = c1 + log c2 λ c2
1 k1 + log y2 . λ
| |
= y2
|
−∞
|
| |
∞
Adem´ as as y2 no cambia de signo pero como log y2 tiende a cuando y2 tiende a 0 y a + cuando y2 tiende a + se ve que y1 s´ı cambia de signo. Dibujemos las curvas (y (y1 (t), y2 (t)) en el caso en que λ < 0.
| |
∞
64
´ 6. COMPORTAMIENTO ASINTOTICO DE LAS SOLUCIONES
y
2
y
1
Lo que da como diagrama de fases
x 2
c
ξ
1
x 1
c
ξ
2
2. DIAG DIAGRA RAMA MAS S DE FASES ASES DE SIST SISTEM EMAS AS LINE LINEAL ALES ES A COEF COEFIC ICIE IENT NTES ES CONS CONST TANTE ANTES S
Caso V. λ1 = α + iβ , λ2 = α
− iβ con β < 0.
65
En este caso la soluci´on o n real es X (t) = c1 Re e(α+iβ )t ξ1 + c2 I m e(α+iβ )t ξ1 donde ξ1 = v1 + iv2 es un autovector asociado a λ1 . Tenemos por lo tanto X (t) = c1 eαt (cos βt v1
{
− sen βt v2) + c2eαt(sen βt v1 + cos βt v2).
}
Escrito en la base v1 , v2 , X (t) = eαt (c1 cos βt + c2 sen βt) βt )v1 + eαt ( c1 sen βt + c2 cos βt) βt )v2 = y1 (t)v1 + y2 (t)v2 .
−
Escribamos (c (c1 , c2 ) en la forma c1 = r cos θ c2 = rsen θ y recordemos que X (0) (0) = c1 v1 + c2 v2, por lo tanto y1 (0) = c1 , y2 (0) = c2 . Tenemos y1 (t) = eαt r(cos θ cos βt + sen θsen βt) βt ) = eαt r cos(θ cos(θ
− βt) βt ) y2 (t) = eαt r(− cos θ sen βt + sen θ cos βt) βt ) = eαt r sen(θ sen(θ − βt) βt ). Por lo tanto en el plano (y (y1 , y2 ) la curva (y (y1 (t), y2 (t)) se obtiene rotando el punto (c (c1 , c2 ) un αt angulo ´angulo βt > 0 y luego expandiendo (o contrayendo) contrayendo) su m´odulo odulo por p or un factor e . Esto da los siguientes diagramas dependiendo de α.
−
α=0 y
2
y
1
´ 6. COMPORTAMIENTO ASINTOTICO DE LAS SOLUCIONES
66
α<0 y
2
y
1
α>0 y
2
y
1
En el plano de fases tenemos
x1 (t) x2 (t)
=Q
y1 (t) y2 (t)
2. DIAG DIAGRA RAMA MAS S DE FASES ASES DE SIST SISTEM EMAS AS LINE LINEAL ALES ES A COEF COEFIC ICIE IENT NTES ES CONS CONST TANTE ANTES S
67
donde Q = [v1 v2 ] es la matriz cuyas columnas son los vectores v1 y v2 . Como una transformaci transformaci´ on o´n lineal transforma c´ c´ırculos en elipses, el diagrama de fases queda, dependiendo dependiendo de α de la siguiente manera.
α=0 x 2
x 1
α<0 x
2
x
1
68
´ 6. COMPORTAMIENTO ASINTOTICO DE LAS SOLUCIONES
α>0 x
2
x
1
De los diagra diagramas mas de fases fases que hemos esbozado esbozado podemos podemos conclu concluir, ir, en partic particula ular, r, que si todos los autovalores de A tienen parte real negativa, todas las trayectorias tienden a 0 cuando el tiempo tiende a + . Si todos todos los autov autoval alor ores es de A tienen parte real positiva, todas las trayectorias tienden a 0 cuando el tiempo tiende a , es decir, se alejan de 0 cuando el tiempo avanza. avanza. Esto es as´ as´ı tanto tanto si los autov autovalores son reales como si son complejos complejos conjugados. conjugados. La diferencia es el modo en que se acercan o alejan de 0.
∞
−∞
Si un autov autovalor alor es positiv positivoo y el otro otro negativ negativo, o, hay hay exacta exactamen mente te dos trayect trayectori orias as que se acercan a 0 cuando el tiempo tiende a + . Estas Estas trayect trayectori orias as son las que correspon corresponden den a la recta de autovec autovectores tores asociados al autov autovalor negativo. Y hay exactamente exactamente dos trayectorias trayectorias que se acercan a 0 cuando el tiempo tiende a (se alejan del origen cuando el tiempo avanza). Estas trayectorias son las que corresponden a la recta de autovectores asociados al autovalor positivo. positivo. Todas las dem´ as trayectorias se alejan del origen tanto hacia el futuro como hacia el as pasado.
∞ −∞
3. Linearizac Linearizaci´ i´ on on
Ahora que entendemos los diagramas de fases de sistemas lineales con coeficientes constantes estamos en condiciones de esbozar los diagramas de fases de sistemas no lineales cerca de puntos de equilibrio. equilibrio. En esta secci´ seccion o´n supondremos que el campo F es continuamente diferenciable en n. R Recordemos que un punto de equilibrio es un cero de la funci´on on F ( F (X ). ). Sea entonces X 0 un equilibrio, tenemos para X cerca de X 0 , F ( F (X )
∼ DF ( DF (X 0 )(X )(X − X 0 ).
´N 3. LINEARIZA LINEARIZACI CIO
69
− X 0, tendremos tendremos Y = (X − X 0 ) = X = F ( F (X ) ∼ DF ( DF (X 0 )(X )(X − X 0 ) = DF ( DF (X 0 )Y
Por lo tanto, si llamamos Y = X para Y
∼ 0. El sistema
Y = DF ( DF (X 0 )Y es un sistema lineal con coeficientes constantes (con matriz A = DF ( DF (X 0 )). Lo que que tenem tenemos os entonces es que para X (t) cerca de X 0 , X (t) X 0 es parecido a la soluci´on on Y ( Y (t) de este sistema cerca de Y 0 = 0.
−
Enunciemos ahora sin demostraci´on on el resultado que asegura que en efecto esto es as´ as´ı si los autovalores de la matriz DF ( DF (X 0 ) tienen parte real no nula. 6.1 (Estabilidad Lineal). Sea F un campo C 1 en R2 . Sea Sea X 0 un cero de F . F . Si DF ( DF (X 0 ) no tiene autovalores con parte real 0, el diagrama de fases del sistema Teorema
X = F ( F (X ) en un entorno de X 0 es muy parecido al diagrama de fases del sistema Y = DF ( DF (X 0 )Y cerca de Y 0 = 0. Con ´esto esto queremos queremos decir decir que hay una biyec biyecci´ on continuamente diferenciable con inversa continuamente diferenciable entre un entorno de X 0 y un entorno del 0 que manda trayectorias del sistema X = F ( F (X ) en trayectorias del sistema Y = DF ( DF (X 0 )Y preservando su orientaci´ on. En particular, si todos los autovalores de DF ( DF (X 0 ) tienen parte real negativa se sigue que todas las trayectorias que pasan cerca de X 0 tienden a X 0 cuando t tiende a + . Y si todos todos tienen parte real positiva, todas las trayectorias se alejan de X 0 cuando el tiempo crece ( X X (t) X 0 cuando t ).
∞
→ −∞
→
M´ as a´ un, si DF ( DF (X 0 ) tiene un autovalor λ1 > 0 y un autovalor λ2 < 0 hay dos trayectorias del sistema X = F ( F (X ) que tienden a X a X 0 cuando el tiempo tiende a + a + . Estas trayectorias son tangentes en X 0 a la recta de autovectores asociados a λ2 . An´ alogamente, hay dos trayectorias del sistema X = F ( F (X ) que tienden a X 0 cuando el tiempo tiende a . Estas traye trayectoria ctoriass son tangentes en X en X 0 a la recta de autovectores asociados a λ a λ1 . Todas las dem´ as trayectorias que pasan cerca de X 0 se alejan de X 0 tanto hacia el futuro como hacia el pasado.
∞ −∞
6.2. Apliquemos este resultado para estudiar el diagrama de fases del sistema simbi´ otico otico del comienzo del cap´ cap´ıtulo cerca de los equilibrios (0, (0, 0) y (γ/δ, (γ/δ, α/β ). ). Ejemplo
−
En este caso el campo F es ( α + βy) βy )x, ( γ + γ + δx) δx)y . De aqu´ı que DF ( DF (x, y) = Por lo tanto,
−
α + βy δy
DF (0 DF (0,, 0) =
−
−
−
−
α 0
−
βx γ + γ + δx
0 γ
−
que tiene autovalores λ1 = α, λ2 = γ . Por lo tanto todas las trayectorias cercanas al (0, (0, 0) se acercan a ´el el cuando el tiempo ti empo tiende a + como suger´ suger´ıa el diagrama que obtuvimos siguiendo las flechas, a saber
∞
´ 6. COMPORTAMIENTO ASINTOTICO DE LAS SOLUCIONES
70
Por otro lado,
0 DF ( DF (γ/δ,α/β ) = δα/β
βγ/δ 0
√
−√
que tiene como autovalores las raices del polinomio λ2 αγ . Es decir, λ1 = αγ , λ2 = αγ . De nuevo, vemos que el diagrama de fases cerca de este equilibrio es como lo esbozamos al comienzo del cap´ cap´ıtulo ya que ser´a parecido al del Caso I de la secci´ on on anterior. anterior. M´ as a s a´ un, un, podemos saber c´omo omo van a salir (o entrar) las trayectorias al equilibrio (δ/γ, (δ/γ, α/β α/β ) ya que son tangentes a las rectas de autovectores de la matriz 0 βγ/δ δα/β 0
−
Para encontrar estas rectas debemos resolver por un lado x1 αγ βγ/δ =0 δα/β αγ x2 es decir
√
√−
−
√αγ x1 − βγ x2 = 0 δ
lo que da la recta
x2 =
α δ x1 . γ β
Por lo tanto, hay una trayectoria tangente a esta recta que sale del equilibrio (δ/γ,α/β ( δ/γ,α/β ). ).
−√αγ tenemos que resolver −√αγ −√βγ/δ x1 = 0 −δα/β − αγ x2
An´ alogamente, para el autovalor alogamente,
es decir
√αγ x1 + βγ x2 = 0 δ
lo que da la recta x2 =
α δ x1 . γ β
−
Por lo tanto, hay una trayectoria tangente a esta recta que entra al equilibrio (δ/γ, ( δ/γ, α/β α/β ). ). Esto se ve en el diagrama de fases que esbozamos en el ejemplo 6.1 6 .1 al comienzo del cap´ cap´ıtulo. Cerca del equilibrio el diagrama de fases es
´N 3. LINEARIZA LINEARIZACI CIO
71
6.3. La ecuaci´ on on del p´ endulo endulo simple amortiguado es g x + sen x + cx = 0 L donde c > 0, g > 0 es la constante de la gravedad, L > 0 es la longitud del p´endulo endulo y x es el angulo ´angulo que el p´endulo endulo forma con la vertical que apunta hacia aba jo. Esta ecuaci´ on on es equivalente al siguiente sistema en el plano de fases (x, (x, y) donde y representa la velocidad de variaci´on on del angulo, ´angulo, x = y g y = sen x cy L g Este es un sistema de la forma X = F ( F (X ) con F = (y, ( y, sen x cy). cy). Los equilibrios equilibrios del sistema sistema L g (los ceros de F ) F ) son los punto (x (x0 , y0 ) tales que y0 = 0, sen x0 cy0 = 0. Es deci decir, r, y0 = 0, L sen x0 = 0. Para ´angulos angulos x0 entre π y π los equilibrios son ( π, 0), (π, (π, 0), (0, (0, 0). Observemos que f´ısicamente los puntos ( π, 0) y (π, (π, 0) represen representan tan la misma posici´ on on y velocidad velocidad en el espacio. Ejemplo
−
−
− −
−
−
−
−
Analicemos el diagrama de fases cerca de los equilibrios. Para ´esto esto veamos cual es el sistema linearizado X = DF ( DF (X 0 )X. Se tiene F = (y,
− Lg sen x − cy). cy). Por lo tanto DF ( DF (x, y ) =
−
0 g cos x L
1 c
−
y DF (0 DF (0,, 0) =
−
0 g L
1 c
−
−
c2 4g/L . 2 Se tiene que ambas raices tienen parte real negativa. Por lo tanto todas las trayectorias cercanas al (0, (0, 0) tienden a (0, (0, 0) cuando el tiempo tiende a + . Pero la forma en que tienden depende de la magnitud de la amortiguaci´ on dada por la constante c. En efecto, si c2 4g/L, on g/L, los autovalores de la matriz del sistema linearizado son reales y las trayectorias tienden al equilibrio sin oscilar alrededor al rededor de ´el. En cambio, si c2 < 4g/L, g/L, los autovalores son complejos conjugados de parte real negativa y las trayectorias se acercan al equilibrio en forma espiral.
−c ± que tiene por autovalores las raices del polinomio λ2 +cλ+ cλ+g/L = 0, es decir λ1,2 = ∞
≥
72
´ 6. COMPORTAMIENTO ASINTOTICO DE LAS SOLUCIONES y
x
Diagrama de fases del p´ endulo endulo subamortiguado cerca del (0, (0 , 0)
Con respecto al ´angulo angulo x del p´endulo, endulo, ´esto esto dice que cerca de la posici´ posi ci´ on on de equilibrio x = 0, velocidad x = 0, el angulo ´angulo tender´a a cero sin oscilaciones si est´ a muy amortiguado y oscilando indefinidamente alrededor de la posici´on on de equilibrio si est´a subamortiguado. subamortiguado. Esto es exactamente lo que se observa para un resorte, lo que no es casual porque la ecuaci´on on del resorte es la ecuaci´on on linearizada alrededor de x = 0. 0. Analicemos Analic emos qu´e pasa para el equilibri equ ilibrioo (π, (π, 0). El sentido com´ un nos dice que es un equilibrio un muy inestable y que en esa posici´on o n con velocidad no nula por chica que sea nos vamos a alejar de ah´ ah´ı y tambi´ en en que el p´endulo endulo va a caer de cualquier posici´ on on por cercana que sea. Estas son situaciones en las que el vector velocidad apunta alej´andose andose del equili equilibri brio. o. Pero Pero tambi´ en en podemos imaginarnos que si le damos el impulso correcto prodr pro dr´´ıamos llegar arriba (a la posici´ on on x = π o x = π ) con velocidad 0. Por supuesto que el impulso debe ser el correcto y una peque˜ na na variaci´ on on podr p odr´´ıa hacer que no lleguemos o que nos pasemos (que lleguemos con velocidad positiva). p ositiva). Veamos que ´esto esto se ve linearizando alrededor de (π, 0). Es decir, que vemos que casi todas las trayectorias se alejan de este equilibrio y que hay exactamente una que se acerca con x < π. Para ´esto esto veamos que los autovalores autovalores de DF ( DF (π, 0) son de distinto signo. En ese caso, el diagrama de fases cerca del equilibrio (π, (π, 0) es similar al del Caso I de los sistemas lineales con coeficientes constantes que da exactamente la situaci´on on que describimos (s´olo olo que en este caso nos restringimos a la regi´on on x < π que corresponde a x < 0 para el problema linearizado). En efecto,
−
DF ( DF (π, 0) =
0 g L
1 c
−
4g/L c2 + 4g/L 2 una de las cuales cuales es positiva positiva y la otra otra negati negativ va. Y se tiene tiene lo afirmad afirmado. o. Si queremo queremoss ver ver c´ omo omo entra la trayectoria con la que nos acercamos al equilibrio desde x < π, debemos encontrar la c c2 + 4g/L 4g/L recta de autovectores asociados al autovalor negativo ya que la trayectoria es 2 tangente a esta recta en (π, (π, 0). Para ´esto esto debemos debem os resolver re solver
−c ± que tiene por autovalores las raices del polinomio λ2 +cλ−g/L = 0, es decir λ1,2 = −−
(c +
c2 + 4g/L 4g/L))/2 g/L
( c+
−
1 2 c + 4g/L 4g/L))/2
x y
=0
´N 3. LINEARIZA LINEARIZACI CIO
73
c2 + 4g/L 4g/L Es decir, y = x. Esto nos da asint´ oticamente oticamente la velocidad que el p´endulo endulo 2 debe tener al pasar por el ´angulo angulo x si queremos llegar a la posici´on on x = π con velocidad 0.
−
c+
El diagrama de fases cerca del (π, (π, 0) ser´a
o π
2.5
Diagrama de fases del p´ endulo endulo amortiguado cerca del ( π, 0)
¿Qu´ e pasa si tratamos t ratamos de analizar anali zar la ecuaci´ on on del d el p´endulo endulo sin amortiguaci amort iguaci´ on? o´n? Esto es, g x + sen x = 0. L En este caso el sistema equivalente en el plano de fases es
x = y
− Lg sen x
y =
que tiene los mismos equilibrios del caso amortiguado. Sin embargo, en este caso DF ( DF (x, y) =
−
0 g cos x L
1 0
y DF (0 DF (0,, 0) =
−
0 g L
1 0
que tiene por autovalores las raices del polinomio λ2 + g/L = 0. Es decir, decir, λ1,2 = i g/L que tienen parte real 0. Por lo tanto el teorema de estabilidad lineal no nos dice nada en este caso. Sin embargo, hay otra forma de analizar la ecuaci´on on del p´endulo endulo ya que es un caso particular de sistema conservativo.
±
Un sistema es conservativo cuando tiene la forma X = U ( U (X ). ). En el caso particular de inc´ ognita ognita escalar – x = U (x) – esta ecuaci´on on es equivalente a un sistema de 2 2 para el cu´al al podemos dibujar el diagrama de fases. Esto ser´a objeto de la pr´oxima oxima secci´ on. on.
−
−
×
´ 6. COMPORTAMIENTO ASINTOTICO DE LAS SOLUCIONES
74
4. Sistemas Conservati Conservativos vos
En esta secci´on on estudiaremos ecuaciones de la forma x = U (x). La funci´on on U ( U (x) se llama potencial del campo conservativo conservativo y representa una forma de energ´ energ´ıa llamada Energ´ Energ´ıa Potencial. La energ ener g´ıa mec´anica anica del sistema es la suma de las energ´ıas ıas cin´etica etica y potencial. Para sistemas conservativos la energ´ıa ıa mec´anica anica se conserv conservaa (de ah´ ah´ı el nombre). En efecto, si x(t) es una trayectoria, la energ´ energ´ıa sobre esa trayectoria trayectoria es
−
1 E (t) = x˙ 2 (t) + U (x(t)), )), 2 por lo tanto d E (t) = x˙ (t)x ¨(t) + U (x(t))x˙ (t) = x˙ (t) x ¨(t) + U (x(t)) = 0. dt De aqu aq u´ı que la energ´ ene rg´ıa ıa E permanezca constante sobre cada trayectoria.
En el plano de fases se tiene el sistema equivalente
x = y y =
−U (x)
1 2 y + U ( U (x). Por 2 lo visto arriba (conservaci´ on on de la energ´ energ´ıa) se tiene que E (x(t), y(t)) es constante sobre cada trayectoria de este sistema. Esto nos dice que las trayectorias est´an an contenida contenidass en los conjuntos conjuntos de nivel de la funci´on on E (x, y ). La energ energ´ıa se represen representa ta en el plano de fases por la funci´ funci´ on on E (x, y) =
Los equilibrios son los puntos (x (x0 , y0 ) tales que y0 = 0 y U (x0 ) = 0. Es decir, decir, los punt puntos os de la forma (x (x0 , 0) con x0 punto cr´ cr´ıtico de la funci´on on U . Observemos que E (x, t) = (U (x), y) se anula s´olo olo en los equilibrios. Por lo tanto, los conjuntos de nivel que no contienen equilibrios est´an an formados por curvas suaves y cada una de estas curvas es una trayectoria del sistema.
Una forma f´acil acil de ver c´omo omo graficar el diagrama de fases es observar que si una trayectoria pasa en alg´ un momento por un punto de la forma (¯ un x, x, 0), la energ´ energ´ıa sobre esa trayectoria ser´a igual a U (x¯). Entonces, 1 U (¯ U (¯ x) = E (x(t), y(t)) = y2 (t) + U ( U (x(t)) 2
≥ U ( U (x(t))
para to do t.
Es decir, si en alg´ un instante una trayectoria pasa por el punto (¯ un x, x, 0), en todo instante permanece en el pozo potencial U (x) U (x ¯)
≤
y U ( U (x(t)) no puede volver a ser igual a U (x ¯) hasta que no se tenga simult´ aneamente aneamente y(t) = 0. Observemos que de todos modos puede haber trayectorias que nunca corten al eje x. Esto es as´ı si la funci´ func i´on on potencial potencial U es acotada superiormente. superiormente. En efecto, supongamos que est´ a acotada superiormente por la constante M . Si la energ´ energ´ıa sobre la traye t rayectoria ctoria es mayor que M , (y ´esto est o es as´ as´ı si en alg alg´ u´n instante la velocidad y es muy grande), no podremos tener nunca y = 0. Es un decir, no se corta al eje x.
||
Por otro lado observemos que a medida que y crece, U (x) decrece d ecrece y rec re c´ıprocamente, ıpro camente, cuando y decrece, U (x) crece.
||
4. SISTEMAS CONSERVATIVOS
75
Finalmente, observemos que por la paridad de la funci´on on E (x, y) en la variable y se tiene que los conjuntos de nivel de E son sim´etricos etricos respecto respe cto del eje x. Estamos entonces en condiciones de graficar el diagrama de fases (en forma aproximada porque no conocemos exactamente los conjuntos de nivel de la funci´on on E ). ). Vamos a hacerl hacerloo g primero en el caso del p´endulo. endulo. En este caso el potencial es U (x) = (1 cos x). Para esbozar L el diagrama de fases conviene dibujar simult´aneamente aneamente el potencial U y el diagrama ya que, como vimos, hay una gran correspondencia entre los dos gr´aficos. aficos.
−
U(x)
−π
π
x
π
x
y
o
o
o
−π
En el gr´afico afico observamos que hay dos trayectorias que conectan los equilibrios inestables ( π, 0) y (π, (π, 0). En una de esas trayect trayectori orias as se sale de ( π, 0) en tiempo t = y se llega a (π, 0) en tiempo t = + . En la otra, otra, se sale sale de (π, (π, 0) en tiempo t = y se llega a ( π, 0) en tiempo t = + . Todas las otras trayectori trayectorias as corresponden corresponden a curvas curvas cerradas. cerradas. Esto es claro en las que se encuentran dentro de las trayectorias que conectan ( π, 0) y (π, (π, 0). Pero tambi´ tamb i´en en las que est´an an fuera corresponden a curvas cerradas si recordamos que estamos identificando los
−
∞
∞
−
−∞
−
−∞
−
´ 6. COMPORTAMIENTO ASINTOTICO DE LAS SOLUCIONES
76
−
angulos ´angulos π y π y que los conjuntos de nivel de E (x, y ) son en este caso sim´etricos etricos respecto resp ecto del g eje y porque la funci´on on potencial U ( U (x) = (1 cos x) es par. Estas trayectorias, que no cortan L el eje x, corresponden corresp onden a niveles n iveles de energ´ıa ıa mayores que el m´aximo aximo de U ( U (x). Para Para estos niveles niveles de energ e nerg´´ıa, el p´endulo endulo da vueltas vuelta s sin si n parar. p arar. Para niveles de energ´ıa ıa menores que el m´ aximo aximo de U ( U (x), el p´endulo endulo alcanza un angulo a´ngulo m´ aximo aximo x¯ tal que U (¯ U (¯ x) = E = max U ( U (x(t)). )). es el valo valorr m´ aximo aximo que alcanza U sobre la trayectoria trayectoria e igual al a l nivel de energ´ energ´ıa de esa trayectoria. Cuando llega a este valor valor lo hace con velocidad velocidad nula. Entonces Entonces el movimiento movimiento cambia cambia de sentido sentido y lo que observamos, observamos, en definitiva definitiva es el movimiento oscilatorio oscilat orio que asociamos con la idea de un p´endulo. endulo.
−
Adem´ as, as, vemos que el (0, (0, 0) es un equilibrio estable en el sentido de que las trayectorias que comienzan comienzan cerca de ´el el permanecen permanecen cerca, cerca, pero no tienden tienden a (0, (0, 0) como ocurre en el caso del p´endulo endu lo amo amorti rtigua guado. do. Los conjuntos de nivel correspondientes a niveles de energ´ energ´ıa menores que max U ( U (x(t)) constan de una sola componente que es una trayectoria. En cambio para niveles de enrg´ enrg´ıa mayores se tienen dos componentes componentes (dos trayectori trayectorias). as). Observe Observen n que para el nivel nivel de energ´ energ´ıa igual al maximo de U hay 4 trayectorias, las dos que describimos antes que conectan los equilibrios ( π, 0) y (π, (π, 0) y las dos trayectorias estacionarias correspondientes a estos equilibrios.
−
Para afianzar las ideas veamos otros ejemplos. 6.4. Supongamos que nos dan el gr´afico afico del potencial U y esbocemos el diagrama de fases correspondiente. Sea entonces el gr´afico afico de U Ejemplo
U(x)
x
El diagrama diagrama de fases correspondien correspondiente te dibujando dibujando tambi´ tambi´ en en el potencial potencial en el mismo gr´ afico para ayudarnos ser´a
4. SISTEMAS CONSERVATIVOS
77
U(x)
x
y
o
o
o x
Para dibujar el diagrama de fases observamos que, como U (x) + cuando x todos los conjuntos de nivel cortan al eje x. Por lo tanto, lo m´as as f´ acil es dibujarlos a partir de un acil punto en este eje. Adem´as, as, como son sim´etricos etricos respecto de este eje, los dibujamos para y > 0 y los completamos comple tamos en forma sim´etrica etrica despu´es. es.
→ ∞
→ ±∞
Empezamos entonces en un punto (¯ x, x, 0) y dibujamos la parte de la trayectoria que pasa por ese punto contenida en y > 0. Como Como y > 0 se tiene que U ( U (x(t)) < U (¯ U (¯ x). Por lo tanto tanto la trayectoria debe moverse hacia valores de x en los lo s que ´esto esto suceda. s uceda. M´ as as a´ un, un, y crece cuando U ( U (x) decrece y cuando U ( U (x) comienza a crecer y decrece. Eventualmente, Eventualmente, U ( U (x) alcanza nuevamente el valor U (¯ U (¯ x) y en ese punto la trayectoria vuelve a cortar al eje x. Veamos otro ejemplo. Ejemplo
6.5. Supongamos ahora que el potencial tiene el siguiente gr´afico. afico.
U(x)
x
78
´ 6. COMPORTAMIENTO ASINTOTICO DE LAS SOLUCIONES
En este caso, como U est´ a acotada superiormente, hay trayectorias que no cortan al eje x. De todos modos, siguen las subidas y bajadas del potencial U con decrecimiento y crecimiento, respectivamente de y . Otras trayectorias trayectoria s s´ı cortan c ortan el eje x. Las que lo hacen para valores de x ¯ menores que el punto donde alcanza el mayor m´aximo aximo relativo o para valores mayores que el punto donde alcanza el menor m´ aximo relativo, son trayectorias contenidas en x x aximo ¯yx x ¯ respectivamente ya que s´olo olo hacia esos lados decrece U . U .
≤
≥
El diagrama diagrama de fases correspondien correspondiente te dibujando dibujando tambi´ tambi´ en en el potencial potencial en el mismo gr´ afico ser´ a U(x)
x
y
o
o
o x
El conjunto de nivel correspondiente al m´aximo aximo de U est´a formado por el equilibrio (x (x0 , 0) con U ( U (x0 ) el m´ aximo aximo de U , U , una trayectoria que entra a (x (x0 , 0) desde valores de x menores que x0 y una desde valores mayores, una que sale hacia valores menores que x0 y una hacia valores mayores. Correspondiente al otro m´aximo aximo relativo se tiene un equilibrio, una trayectoria que sale y vuelve a entrar al equilibrio con valores de x menores y dos trayectorias con valores de x mayores, una que entra y una que sale.
EJERCICIOS
79
Las u unicas ´ nicas trayectorias cerradas se encuentran alrededor del equilibrio correspondiente al m´ınimo ıni mo relati rel ativo vo de U que resulta ser el ´unico unico equilibrio estable. Las flechas que indican el sentido de recorrido de las soluciones se obtienen de la observaci´on general para sistemas conservativos que mientras m ientras se est´e en el semiplano superior x crece y en el inferior x decrece por ser y = x˙ .
Ejercicios
(1) Realice un gr´afico afico aproximado de las l as l´ıneas ıneas de flujo de los siguientes si guientes campos vectoriales:
−
(a) F ( F (x, y) = (y, x)
(c) F ( F (x, y) = (x, x2 )
−
(b) F ( F (x, y ) = (x, y)
(2) Para el siguiente sistema de dos poblaciones que compiten por un mismo alimento, realice un gr´afico afico aproximado del diagrama de fases a partir del dibujo del campo vectorial. x˙ = (2 x y )x
− − y˙ = (3 − x − 2y)y
(3) Esboce el diagrama de fases de los sistemas lineales de los siguientes ejercicios del Cap Ca p´ıtul ıt uloo 4: Ejer Ejerci cici cioo 1, (b) (b) y (c); (c);
Ejer Ejerci cici cioo 3, (b) (b) y (d) (d)
(4) Sea A R2×2 una matriz cuyos autovalores son λ y µ. Esbozar Esbozar el diagrama diagrama de fases correspondiente al sistema X = AX si
∈
(a) λ > µ > 0 (c) λ > 0 > µ (e) 0 = λ = µ
(b) 0 > λ > µ (d) λ = α + iβ, µ = α (f) (f ) λ = 0, µ > 0
∈R
− iβ con β = 0
(5) Para los siguientes sistemas hallar los puntos de equilibrio y esbozar el diagrama de fases cerca de cada uno de ellos. x˙ = xey x˙ = ex−y 1 (a) (b) y˙ = sen x + y 1 y˙ = xy 1
−
−
−
(6) Considerar los siguientes sistemas de depredador–presa con crecimiento log log´´ıstico.
− − − − Hallar los puntos de equilibrio (x (x0 , y0 ) con x0 ≥ 0, y0 ≥ 0 y esbozar el diagrama de (a)
− − − −
x˙ = (2 ( 2 x y )x y˙ = ( 1 + x y )y
(b)
x˙ = (1 x y)x y˙ = ( 2 + x y)y
fases en la cercan´ cercan´ıa de los mismos. Tratar de ver c´omo omo es el diagrama global. Observ Observar que el comportamien comportamiento to puede ser muy distinto dependiendo de los par´ametros ametro s de la ecuaci´on. on. (7) Considerar el sistema de ecuaciones diferenciales que modela el flujo en R2 asociado a un campo vectorial gradiente, es decir un sistema de la forma X = V ( V (X ) con 2 2 V C (R ). Hallar los puntos puntos de equilibrio equilibrio del sistema e investigar investigar su estabilidad estabilidad si los extremos locales de la funci´on on V son no degenerados (es decir, si los autovalores del Hessia Hessiano no de V en los extremo extremoss locales locales de V son no nulos). nulos). ¿Puede ¿Puede decir decir si hay hay alguna relaci´ on on entre entre las trayector trayectorias ias del sistema y las l´ıneas de nivel nivel de la funci´ funci´ on on
∈
−
´ 6. COMPORTAMIENTO ASINTOTICO DE LAS SOLUCIONES
80
V ( V (X ) = V ( V (x1 , x2 )? Utilice Utilice esta esta informa informaci´ ci´ on para esbozar el diagrama de fases del on 2 sistema X = V ( V (X ) si V ( V (x1 , x2 ) = x1 + 4x 4x22 .
−
(8) Si la fuerza de atracci´on on entre dos masas que se encuentran a distancia r es de la forma k F = r2 (a) Hallar la energ´ energ´ıa potencial potencial y esbozar el diagrama de fases si el potencial potencial tiende a cero cuando r = + . (b) Si a una distancia r0 la energ´ ene rg´ıa ıa cin´ ci n´etic et icaa es T 0 < U ( U (r0 ). ¿Cu´ ¿Cu´ al al es la distancia m´ axima axima de separaci´on on posible de estas masas? (c) ¿Qu´e sucede si la energ´ıa ıa total es i. posi positi tiv va, ii. ii. nega negati tiv va, iii. iii. nula? la?
−
∞
−
(9) Supongamos que la fuerza de atracci´ on on entre los ´atomos atomos de una mol´ecula ecula diat´ omica omica es 1 a de la forma F ( F (x) = + 3 donde x es la distancia entre los mismos. 2 x x (a) Hallar la energ´ energ´ıa potencial potencial y esbozar el diagrama de fases si el potencial potencial tiende a 0 cuando x tiende a infinito. (b) Utilizando el diagrama de fases, observar que (i) la distancia entre los ´atomos atomos permanece p ermanece constante si y s´ olo olo si en alg´ un un momento se encuentran a distancia a y velocidad 0. (ii) Si la energ´ıa ıa total to tal E 0 es negativa, la distancia entre los atomos a´tomos crece y decrece en forma oscilatoria entre dos valores m´aximo aximo y m´ınimo dependientes s´ olo olo de E 0 . (iii) Si la energ´ıa ıa total es no negativa, la distancia entre los atomos a´tomos tiende a infinito cuando el tiempo tiende a infinito aunque pueden acercarse inicialmente. (iv) (iv ) ¿Cu´al al es la energ ene rg´´ıa m´ınima ıni ma posib po sible le del sistem sis tema, a, E min min ? (v) En todos los casos, ¿cu´ al al es la distanci di stanciaa m´ınima ınima entre los lo s atomos a´to mos si la energ ener g´ıa de la mol´ecula ecul a es E 0 E min min ?
−
≥
Agradecimientos Quisiera agradecer profundamente a Juli´ an an Fern´andez andez Bonder por sus comentarios y sugerencias respecto del material y la presentaci´ on on de estas notas. Tambi´ en en mi agradecimiento a Gabriel Acosta, Gabriela Armentano, Javier Etcheverry, Pablo Groisman y Ariel Lombardi por su inestimable ayuda con los gr´aficos. aficos.
81
Bibli ib liog ogra raff´ıa
[1] Amann, H., “Ordinary differential equations”, Walter Walter de Gruyter, Berlin, 1990. [2] Ayres, F. “Ecuaciones Diferenciales”, Diferenciales”, Colecci´ on on Schaum. 1969. [3] Birkhoff, G. and Rota, G.C. “Ordinary Differential Differential equations”, Ginn & Company, 1962. [4] Coddington, E.A., “Introducci´ “Introducci´ on a las ecuaciones diferenciales ordinarias”, Compa˜ on n´ n´ıa ıa Editori Edi torial al Continent Cont inental, al, SA, 1ra. ed. en espa˜ nol, nol, 7ma. ed. en e n ingl´ ing l´es, es, 1968. [5] Coddington, E.A. & Levinson, N. “Theory of ordinary differential differential equations”, Mc-Graw Hill, Hill, 1955. [6] de Figueiredo, Figueiredo, D. & Neves, Neves, A., “Equac˜ oes Diferenciais Aplicadas”, Colec˜ oes ao ao Matem´ atica atica Universit´ aria, aria, Instituto de Matem´ atica Pura e Aplicada, Rio de Janeiro, 1997. atica [7] Hirsch, M. W. & Smale, S., “Differential “Differential equations, dynamical systems and linear algebra”, Academic Press, Press, New York, 1974. [8] Hurewicz, W. “Sobre ecuaciones diferenciales diferenciales ordinarias”, Ediciones RIALP, RIALP, Madrid, 1966. [9] INCE, E.L. “Integraci´ “Integraci´ on de ecuaciones diferenciales ordinarias”, Editorial Dossat, S.A., Madrid, 1939. on [10] Perko, Perko, L., “ Differential equations and dynamical systems”, Springer-Verlag, Springer-Verlag, New York, 1991. [11] Pontryagin, Pontryagin, L.S. “Ordinary Differential Equations”, Addison-Wesley Addison-Wesley,, 1962.
83