Ecuaciones Diferenciales Ordinarias: Claudio Fern Andez y Rolando Rebolledo Con La Colaboraci On de Felipe Montoya
Ecuaciones Diferenciales Ordinarias: Claudio Fern Andez y Rolando Rebolledo Con La Colaboraci On de Felipe Montoya
Ecuaciones Diferenciales Ordinarias: Claudio Fern Andez y Rolando Rebolledo Con La Colaboraci On de Felipe Montoya
Ordinarias
3
4 PREFACIO
Prefacio 3
Captulo 1. Introduccion 7
1. El cambio incesante 7
2. La Mecanica segun Newton 8
3. Las Ecuaciones Diferenciales en otras ciencias 13
Captulo 2. Conceptos basicos 15
1. Terminologa basica 15
2. Reduccion de ecuaciones al primer orden 18
3. Un algoritmo para construir soluciones de la ecuacion
fundamental 21
4. Ejercicios propuestos 22
Captulo 3. Ecuaciones de Primer Orden 25
1. Introduccion 25
2. Ecuaciones con variables separables 25
3. Diferenciales exactos 32
4. Factores Integrantes 36
5. Curvas integrales: caso general 38
6. Reduccion de ecuaciones al caso lineal de primer orden 40
7. Sistemas Dinamicos 42
8. Ejercicios propuestos 44
Captulo 4. Ecuaciones Diferenciales Lineales 49
1. La ecuacion lineal homogenea con coeficientes constantes 49
2. Sistemas no homogeneos con coeficientes constantes 58
3. Sistemas homogeneos con coeficientes variables 67
4. Sistemas no homogeneos con coeficientes variables 74
5. Ejercicios propuestos 79
Captulo 5. La transformacion integral de Laplace 91
1. Transformacion de Laplace 91
2. Aplicacion a la resolucion de ecuaciones diferenciales 94
3. Truncamientos y desplazamientos 101
4. Producto de convolucion 104
5
6 CONTENIDOS
Introduccion
1. El cambio incesante
Si usted deseara explicar a un ser inteligente, venido de un mundo
plano, plano el como su mundo, lo que es una esfera, seguramente se
vera en serios problemas de comunicacion. El problema se presentara
en realidad con cada nocion que requiriese de una dimension superior
a 2. En efecto, como representar en un plano una figura que proviene
de un espacio de tres dimensiones? Los pintores enfrentan a menudo
este problema y lo resuelven con la perspectiva. Pero, esta ultima hace
uso de la educacion de nuestros sentidos realizada en un ambiente tridi-
mensional. Nuestro individuo plano no dispondra de ese bagaje. Una
solucion alternativa consistira en la siguiente: se le podra definir el
concepto de esfera diciendole que corresponde a una familia infinita de
circulos concentricos que contiene como punto de partida un punto y
luego el radio aumenta hasta alcanzar un maximo, reduciendose luego
hasta un punto final. Implcitamente estamos definiendo as una esfera
por un mecanismo evolutivo. En efecto, para nosotros, seres que vivi-
mos en un universo con una dimension mayor, la definicion anterior
esta inspirada de la idea de cortar una esfera por un numero infinito
de planos paralelos que podran ser entendidos como el universo plano
del principio, en movimiento de traslacion paralela.
Pero, mejor aun, sabemos que el universo en el cual nosotros nos
movemos tiene mas de tres dimensiones. Podemos decir, para fijar
ideas, que necesitamos tres coordenadas espaciales y una temporal.
Ahora bien, como percibimos nosotros los objetos propios de este es-
pacio de cuatro dimensiones? Copiando la solucion entregada al in-
dividuo plano, podramos decir que percibimos esos objetos mediante
la traza que deja su evolucion en el espacio de tres dimensiones en el
curso de los tiempos. Es por ejemplo la observacion de los albumes de
fotos: observaciones de una persona en distintos instantes de su vida
se considera como representacion de su ser.
En nuestro mundo todo esta sujeto al cambio incesante (evolucion)
y es esa idea fundamental que ha inspirado muchas teoras fsicas y
matematicas. En un principio, ya desde la epoca de Aristoteles, la
7
8 1. INTRODUCCION
luego
t(t + 1)
(2.5) x(t) = x(0) + v(0)t + a , para todo t 1.
2
En la ecuacion (2.5) se puede notar que el termino t(t + 1)/2 corre-
sponde a t(t + t)/2. De modo que si hacemos tender t a cero, en vez
de mantenerlo fijo e igual a 1, entonces x(t) queda representado por
t2
x(t) = x(0) + v(0)t + a .
2
Los dos ejemplos anteriores muestran casos particulares de Ecua-
ciones con Diferencias Finitas. Estas ultimas nos permiten representar
el movimiento de un cuerpo cuando el tiempo es discreto.
Lo que hasta aqu hemos hecho corresponde a la Cinematica de un
punto material. Es decir, el estudio del movimiento de un cuerpo sin
considerar la causa del mismo. El estudio de las causas del movimiento
corresponde a la Dinamica. La Segunda Ley de Newton afirma que
la causa del movimiento es la fuerza externa: ella es proporcional a la
aceleracion alcanzada por el movil. La constante de proporcionalidad
es la masa del cuerpo:
(2.6) F (t) = ma(t).
En los casos anteriores, la fuerza externa es nula, en el primero de
ellos, y es constante e igual a ma en el segundo.
Veamos otro ejemplo de una ecuacion en diferencias finitas que
puede ser resuelta facilmente:
10 1. INTRODUCCION
(2.7) Z(0) = 1,
Observar que
(2.11) Zn (0) = 1
(2.12) Zn (t) = aZn (t 1)t,
cuya solucion es
a nt
(2.13) Zn (t) = (1 + ) .
n
Notese que si n entonces Zn (t) exp(at) que satisface la
ecuacion diferencial
d
(2.14) Z(t) = aZ(t)
dt
con la condicion inicial
(2.15) Z(0) = 1
Ejemplo 2.6.
Un ladrillo de masa m esta sujeto por un resorte que a su vez tiene la
segunda extremidad empotrada en un muro vertical. El ladrillo reposa
sobre una superficie plana, que genera una fuerza de friccion lineal. El
resorte ejerce una fuerza proporcional al desplazamiento con respecto a
la posicion de equilibrio. El origen del sistema de coordenadas se fija en
la posicion de equilibrio del sistema. De modo que si el desplazamiento
es x(t), entonces, la fuerza ejercida por el resorte es F (t) = kx(t),
donde k es constante (constante de elasticidad). Plantear las ecuaciones
del movimiento.
Por la Segunda Ley de Newton:
(2.20) mx (t) = x (t) kx(t)
Y suponemos que las condiciones iniciales son:
(2.21) x(0) = x0 , v(0) = v0
Introduciendo las constantes
k
= = ,
m m
la ecuacion se escribe
(2.22) x + x + 2 x = 0
Veremos mas adelante la forma de resolver esta ecuacion en general.
Por de pronto senalamos que, en el caso de ausencia de roce, se tiene
v0
(2.23) x(t) = sen t + x0 cos t.
3. LAS ECUACIONES DIFERENCIALES EN OTRAS CIENCIAS 13
Conceptos basicos
1. Terminologa basica
Es el momento de precisar algunas de las nociones que han sido
introducidas, progresivamente, a traves de los ejemplos del captulo
anterior. Las definiciones principales las resumimos a continuacion.
Definicion 1.1. A lo largo de todo este texto, consideraremos que
la variable independiente es t R, a menos que lo contrario se explicite.
En la mayora de los casos, esta variable representara el tiempo.
Una ecuacion diferencial ordinaria es una que involucra a t
y una funcion desconocida x(t), as como a algunas de sus derivadas.
Las derivadas se designan por los smbolos usuales x , x(k) . El supremo
de los k tales que x(k) aparece en la ecuacion, se llama orden de la
ecuacion. As la forma general de la ecuacion de orden k es
(1.1) F (t, x(t), x (t), . . . , x(k) (t)) = 0.
Diremos que la ecuacion diferencial de orden k se expresa en forma
normal si podemos despejar x(k) de la expresion (1.1), vale decir:
(1.2) x(k) (t) = f (t, x(t), x (t), . . . , x(k1) (t)).
Se puede observar que F en (1.1) es una funcion de k + 2 variables,
en tanto que f depende de k+1 variables. Supondremos, en general, que
x(t) es una funcion con valores en Rd , de modo que sus derivadas en
cada punto t son tambien vectores de Rd . En consecuencia, para que
(1.1) tenga sentido, es necesario que F este definida en un dominio
D(F ) contenido en Rd(k+1)+1 y con valores en R. Por su parte, f debe
estar definida en un dominio D(f ) Rdk+1 con valores en Rd .
El paso de la forma (1.1) a (1.2) se puede hacer gracias al Teorema
de la Funcion Implcita.
Teorema 1.1. Sea F una funcion definida y de clase C 1 sobre un
abierto U Rn , con valores reales. Sea a U , tal que F (a) = 0. Se
supone ademas que
F
(a) = 0.
xn
15
16 2. CONCEPTOS BASICOS
Observar que
y2 = y2 y3
y4 = y12 + y42 ,
de modo que se obtiene finalmente
(2.12) y = f (y),
donde y es la funcion vectorial de componentes y1 , . . . , y4 y f la funcion
definida por
y2
y2 y3
f (y) =
y4
.
2 2
y1 + y4
4. Ejercicios propuestos
(1) Verificar que la funcion x(t) = t es una solucion del problema
con valores iniciales:
(4.1) x = (x t)2/3 + 1
(4.2) x(1) = 1.
Haciendo el reemplazo y(t) = x(t) t, demuestre que y
verifica:
(4.3) y = y 2/3
(4.4) y(1) = 0.
Compruebe enseguida que
( )3
t
(4.5) y(t) = , (t R),
3
es solucion del problema con valores iniciales para y y deter-
mine luego una segunda solucion para el problema con valo-
res iniciales en x. Que ocurre en este caso con el algoritmo
de resolucion de ecuaciones diferenciales tratado en la ultima
seccion?
(2) Elimine las constantes A y B de cada una de las siguientes
expresiones, obteniendo las ecuaciones diferenciales del menor
orden posible que satisface y(x) o x(t):
(i) y = Ax + B
(ii) y = Ae3x + Be2x
(iii) Ax2 + By 2 = 4
(iv) x = A cos(2t) + B sen(2t)
(v) x = A cosh(3t) + B senh(3t)
(vi) x = (A + Bt)e10t
4. EJERCICIOS PROPUESTOS 23
1. Introduccion
En este captulo nos concentraremos en el estudio de la ecuacion de
primer orden
(1.1) x = f (t, x),
donde t vara sobre R y la funcion desconocida x(t) toma tambien sus
valores en dicho espacio en un primer caso: as entonces, D(f ) R2 .
Posteriormente, permitiremos a x(t) tomar sus valores sobre el espacio
de los numeros complejos C. En ese caso D(f ) es un subconjunto de
R C.
En la proxima seccion estudiaremos el caso en que f (t, x) se de-
scompone en el producto de una funcion del tiempo t y de una de la
variable espacial x.
Vale decir,
1 1
= t t0 ,
(t) (t0 )
2. ECUACIONES CON VARIABLES SEPARABLES 27
En efecto,
t por la nota anterior, basta observar que si llamamos
(t) = 0 a(s)ds su derivada en todo punto t existe por el Teorema
Fundamental del Calculo y se tiene (t) = a(t), de donde:
d d d
(2.28) z(t) = z = a(t)z(t), (t R).
dt d dt
Usando el ejercicio anterior tenemos, entonces, la proposicion que
sigue.
Proposicion 2.2. Dada una funcion a : R C continua, la forma
general de las soluciones de la ecuacion lineal homogenea de primer
orden
(2.29) z = a(t)z
2. ECUACIONES CON VARIABLES SEPARABLES 31
es
( t )
(2.30) (t) = C exp a(s)ds , (t R),
t0
3. Diferenciales exactos
Consideremos una funcion de dos variables U (t, x) que sea diferen-
ciable. Hemos visto en los cursos de Calculo Diferencial, que el dife-
rencial de U se calcula como
U U
(3.1) dU = (t, x)dt + (t, x)dx.
t x
Que significa (3.1)? Significa que en una vecindad del punto (t, x)
podemos aproximar los incrementos infinitesimales de la funcion U por
una funcion lineal de variables locales dt, dx. Supongamos ahora que
en vez de plantear (3.1) en cada punto (t, x), los escogemos segun el
grafico de una funcion x = (t) de modo que sobre ese grafico U (t, x) =
U (t, (t)) sea constante. En ese caso tendremos dU = 0, vale decir
U U
(3.2) (t, x)dt + (t, x)dx = 0,
t x
3. DIFERENCIALES EXACTOS 33
4. Factores Integrantes
Ejemplo 4.1. Resolver la ecuacion diferencial
x2 cos xt
(4.1) x = .
2 sen xt + xt cos xt
Se puede observar por una aplicacion del Teorema 3.1, que la ecuacion
planteada no es exacta. En efecto, la derivada del numerador con re-
specto a x no es igual a la derivada del denominador con respecto a t
Que hacer entonces? No escapara al lector sagaz que si multiplicamos
el numerador y el denominador por x la fraccion no cambia, pero si
ahora denotamos
M (t, x) = x3 cos xt
N (t, x) = 2x sen xt + x2 t cos xt,
entonces,
M N
= 3x2 cos xt x3 t sen xt = .
x t
Y ahora podemos aplicar el Teorema 3.1, introduciendo
U (t, x) = 2x sen xtdx + x2 t cos xtdx + g(t) = x2 sen xt + g(t),
U
= x3 cos xt + g (t) = x3 cos xt,
t
de donde g(t) = c = const. Las soluciones x(t) de la ecuacion propuesta
quedan expresadas, entonces, en forma implcita:
(4.2) x2 (t) sen x(t)t = Const.
4.1. Factores Integrantes. El ejemplo anterior nos sirve de ilus-
tracion para un metodo conocido con el nombre de busqueda de fac-
tores integrantes. El proposito es resolver una ecuacion diferencial de
la forma
P (t, x)
(4.3) x = ,
Q(t, x)
cuando
P Q
(4.4) (t, x) = (t, x),
x t
para algun punto (t, x) en el dominio de definicion comun a P y Q.
4. FACTORES INTEGRANTES 37
donde hemos supuesto por simplicidad que la masa del movil es uni-
taria.
Suponemos que F deriva de un potencial:
(5.4) F (x) = V (x),
donde representa el operador gradiente.
Segun lo visto en la definicion, la integral U debe tener el mismo
dominio de definicion que la funcion que determina la ecuacion dife-
rencial, que en este caso es f (t, x, y) = F (x), (y corresponde a la co-
ordenada que sera llenada por la primera derivada posteriormente).
Introducimos entonces la funcion
1
(5.5) U (t, x, y) = y 2 + V (x).
2
Si x = (t) es una solucion de la ecuacion de Newton, reemplazando
en U tenemos:
d d 1 2
U (t, (t), (t)) = ( (t) + V ((t)))
dt dt 2
= 0
De modo que U es una integral de la ecuacion de Newton: corre-
sponde a la Energa total del sistema. En su expresion, el termino
1 2
2
y corresponde a la energa cinetica del movil; V (x) es su energa
potencial. Se dice que F = U es una fuerza conservativa.
Ejemplo 5.2. Resolver la ecuacion del oscilador armonico sin roce,
vista en el primer captulo: corresponde a la ecuacion del movimiento
de un objeto sobre una superficie plana, adherido a un resorte cuya
segunda extremidad esta empotrada en una pared vertical. Vale decir,
k
(5.6) x + 2 x = 0, ( 2 = ),
m
donde m es la masa del movil y k la constante de elasticidad del resorte
al cual esta adherido.
Aplicando el ejemplo anterior (pero esta vez con la masa m = 1),
tenemos
kx2
V (x) =
2
y la energa cinetica es
1 2
my ,
2
donde y corresponde a la coordenada de la velocidad.
40 3. ECUACIONES DE PRIMER ORDEN
Entonces,
1 kx2
(5.7) U (t, x, y) = my 2 + ,
2 2
es una integral de la ecuacion propuesta: la energa total U debe ser
constante para x = (t), y = (t) donde es solucion de (5.6). Esto
permite despejar en funcion de :
1 k2
m +
2
= const.,
2 2
o sea, debe resolver la ecuacion
(5.8) x = C 2 2 x2 ,
donde C es constante.
Esta ultima ecuacion se resuelve por separacion de variables:
dx
= dt,
C 2 2 x2
y obtenemos
x
Arccos = (t + ),
C
donde es una constante de integracion.
Luego, la solucion de la ecuacion del oscilador armonico se ex-
presa:
C
(5.9) (t) = cos(t + )
(5.10) = c1 cos t + c2 sen t, (t R).
Vemos, as, que la busqueda de las integrales de una ecuacion dife-
rencial ayuda a reducir el orden del problema, pues, sabiendo que U
debe ser constante, dicha relacion provee una ecuacion diferencial de
un orden menor que la inicial.
7. Sistemas Dinamicos
Para terminar este captulo reflexionemos una vez mas sobre el ob-
jeto de este curso. A la luz de las primeras ecuaciones que hemos
aprendido a resolver, podemos reinterpretar el papel que esta Teora
juega en el desarrollo de la Fsica Clasica y de la Mecanica en particu-
lar.
En la Mecanica de Newton un sistema fsico queda representado
por dos datos: su posicion q(t) y su momento p(t) (que es igual a su
masa por la velocidad, de modo que tambien se puede usar q (t) en
vez de p(t)). El par x(t) = (q(t), p(t)) recibe el nombre de estado del
sistema en el tiempo t. Llamemos el espacio de estados, tambien
conocido como espacio de fase. En general, es un subconjunto de
Rd Rd .
La Fsica Clasica asegura, entonces, que el sistema queda comple-
tamente determinado si conocemos
(i) Un estado inicial x(t0 ) = (q(t0 ), p(t0 ));
(ii) su evolucion en un intervalo de tiempo infinitesimal [t, t+dt].
Para escribir rigurosamente la frase (ii) disponemos ahora de los utiles
matematicos apropiados: las ecuaciones diferenciales. En particular,
en el caso de la Mecanica, las ecuaciones de Newton.
7. SISTEMAS DINAMICOS 43
8. Ejercicios propuestos
(1) Resolver las ecuaciones siguientes:
x2
(8.1) x = ,
tx t2
(8.2) x = x x2 .
(2) Resolver la ecuacion diferencial
x
(8.3) x = .
1 + xet
(3) (a) Resolver la ecuacion
(8.4) (2t + 3x2 )dt + 6txdx = 0.
(b) Sean (t, x) y (t, x) dos factores integrantes de la ecuacion:
(8.5) M (t, x)dt + N (t, x)dx = 0.
Suponiendo que (t, x) = (t, x)/(t, x) no es identicamente
constante, demuestre que (t, x) = C = const. define,
implcitamente, la solucion general de (8.5).
(4) La poblacion de una ciudad en 1993 era de 200.000 habitantes
y la tasa de crecimiento estimada es de 4%, medida con re-
specto a la poblacion existente en el instante. Predecir el
tamano de la poblacion al ano 2013.
(5) Cuando un proyectil es disparado en el agua, encuentra una
resistencia proporcional al cuadrado de su velocidad. Si U es
la velocidad inicial y si ella disminuye a U/2, despues de viajar
una distancia S, encontrar la velocidad que el proyectil tiene
al alcanzar una distancia horizontal d y el tiempo que le toma
el alcanzarla, despreciando los efectos de la gravitacion.
(6) La Segunda Ley de Newton, que expresa que la fuerza F es
igual a la masa m por la aceleracion a de un movil, suele
enunciarse tambien en la forma
d
(8.6) F (t) = (m(t)v(t)), (t R),
dt
donde v es la velocidad del cuerpo. Esta formulacion se presta
al analisis de la dinamica cuando la masa no es constante.
Un cohete se lanza verticalmente, propulsado mediante la
eyeccion de material con velocidad constante u (hacia abajo)
relativa al cohete. La masa del cohete en el instante t T ,
(donde T < 1/k), esta dada por m0 (1 kt), donde m0 es la
masa inicial y k es una constante mayor que g/u.
46 3. ECUACIONES DE PRIMER ORDEN
dy y
(ii) + = x3 , y(1) = 2
dx x
dy
(iii) = ey sen x, y(0) = 0
dx
dy
(iv) + y cot x = cos x, y( ) = 1
dx 2
dx
(v) 2t3 = x2 + 3tx2 , x(1) = 1
dt
dx
(vi) t x = t3 cost, x() = 0.
dt
(12) Resuelva las siguientes ecuaciones:
dx
(i) = 8xt
dt
2
dx
(ii) 2 sen2 t = 1
dt
dx
(iii) t + x =0
dt
dx
(iv) t3 = x .
dt
(13) La ecuacion diferencial que modela la deflexion w de un disco
circular cuando la carga esta distribuida simetricamente con
respecto al eje, esta dada por:
[ { }]
1 d d 1 d dw p
r (r ) = ,
r dr dr r dr dr D
donde r es la distancia al eje y D es una constante. Obtenga
la solucion general de esta ecuacion cuando la presion p es
constante.
(14) Si se tiene que:
p + t + r dp
dr
=0
p t = 2
con por determinar, y p(r1 ) = p1 , p(r2 ) = p2 ,
obtenga p(r) y t(r).
(15) Encuentre las soluciones generales de:
t(1 t2 ) dx
dt
+ (2t2 1)x = t3
dx x3 3t2 x
= 3
dt t + 3tx2
dx
= x tant 2 sent
dt
dx
t + x = x2 log t
dt
(x2 3t2 )t dx
dt
= (t2 3x2 )x
dx
t = x + t2 x3
dt
48 3. ECUACIONES DE PRIMER ORDEN
dx t2 + tx x2
=2
dt t2 + x2
(2t 3x + 1)dt + (6x 4t + 3)dx = 0
x(1 + t4 x3 )dt + t d x = 0
(x2 + tx + 2t2 ) dx dt
= x(x + 3t)
[(1 + t ) + x]dt + [(1 + x2 )1/2 + t]dx = 0
2 1/2
(1 dx dt
)tan t = 2x
t dt 2x + t2 x2 = 0
dx
5t4 (1 x) + t2 x3 + (t5 + t3 x2 ) dx dt
=0
d2 x
x dt2 + ( dt ) = 0
dx 3
txdx = t2 dx + x2 dt
t(x tcost) + t = 0
t2 x + tx = 1
(t2 x3 + x)dt = (t3 x2 t)dx
(6t + 4x + 3)dt + (3t + 2x + 2)dx = 0
xx + (x )2 2xx = 0
(x2 3xt 2t2 )dt = (t2 xt)dx
8. EJERCICIOS PROPUESTOS 49
1 + n(x t)
x =
n(x t)
sen x dt + cos x dx et (x sen(xt)dt + t sen(xt)dx) = 0.
CAPTULO 4
Vk = Cjmi , si k = mi1 + j, j mi , i = 2, . . . , q.
De este modo, dado cualquier vector constante C, el se
expresa como una combinacion lineal
C = 1 V1 + . . . + d Vd ,
de donde
(1.13) eAt C = 1 eAt V1 + . . . + d eAt Vd
y las funciones j (t) = exp(At)Vj que se calculan segun (1.12),
constituyen una base de soluciones de la ecuacion diferencial
propuesta.
En otros terminos, si se quiere expresar el caso anterior
mediante matrices de cambio de base, introducimos la matriz
P cuyas columnas son los vectores Vj obtenidos reordenando
los vectores Cji . Ademas, llamamos J a la matriz de Jordan,
que tiene bloques diagonales de la forma Ji = i I + Ni , i =
1, . . . , q y bloques nulos en el complemento. Entonces,
(1.14) exp(At) = P exp(Jt)P 1 ,
donde exp(Jt) es una matriz que tiene bloques diagonales de
la forma exp(Ji t).
Ejercicio 1.1. Encontrar la solucion de cada uno de los sistemas
siguientes:
(1)
{
x = 2x y
y = 2y.
(2)
{
x = 2x y
y = x + 2y.
(3)
{
x = y
y = x.
1. LA ECUACION LINEAL HOMOGENEA CON COEFICIENTES CONSTANTES
57
(4)
x = 2x
y = x 2y
z = y 2z.
(donde
) c1 y c2 son constantes.
p 2
b) 2 < q. Las races son complejas conjugadas:
p p
1 = + i , 2 = i .
2 2
Una base de soluciones esta constituda por las funciones:
p
1 (t) = exp([ + i ]t),
2
p
2 (t) = exp([ i ]t).
2
La solucion general se escribe
pt
(t) = e 2 (c1 eit
+ c2 eit
),
donde c1 y c2 son constantes. De manera equivalente, se
puede usar la representacion de las exponenciales comple-
jas en terminos de senos y cosenos, cambiando las cons-
tantes por otras d1 y d2 de modo de expresar la solucion
general en la forma:
pt
(t) = e 2 [d1 cos(t ) + d2 sen(t )].
Caso 2: = 0. En este caso el polinomio caracterstico posee una raz
doble:
p
= .
2
Una base de soluciones esta dada por:
pt pt
1 (t) = e 2 , 2 (t) = te 2 .
La solucion general se escribe:
pt
(t) = e 2 (c1 + c2 t).
Ejercicio 1.2. Resolver las ecuaciones siguientes:
60 4. ECUACIONES DIFERENCIALES LINEALES
(1) x + 7x + x = 0.
(2) x + 6x + 9x = 0.
(3) x 9x + 24x 16x = 0.
donde [C1 , . . . , Cd ] denota la matriz que tiene por columnas a los vec-
tores Ci . Si hemos escogido los Ci iguales a los de la base canonica de
Cd se tendra [C1 , . . . , Cd ] = I y (t) = exp(At). Buscamos entonces
una solucion p (t) de la ecuacion no homogenea que se escriba de la
forma
t
(2.6) U (t) = 1 (s)b(s)ds
t0
t
(2.7) p (t) = (t) 1 (s)b(s)ds,
t0
u = (A I)1 w.
A. Entonces, si escogemos t0 = 0,
t
p (t) = e At
eAs es vds v es vector propio de B,
0
t
= eAt e(IA)s vds
0
= eAt (I A)1 [e(IA)t I]v
= (I A)1 [et eAt ]v,
pues A y (I A) conmutan.
El calculo anterior nos muestra que p tendra una forma exponen-
cial, determinada por las exponenciales de matrices. En general, si b(t)
es (exp(Bt))v, sus componentes seran de la forma p(t) exp(t) donde
p(t) es un polinomio y es raz del polinomio caracterstico de B.
Asimismo, p (t) tendra una forma similar en dicho caso y en vez de
perdernos en largos calculos, podemos ensayar una solucion p de la
ecuacion no homogenea reemplazando en ella una funcion con com-
ponentes de la forma q(t) exp(t), donde q(t) es un polinomio, cuyos
coeficientes se determinaran por simple identificacion. Ilustremos este
metodo con el siguiente ejemplo.
Ejemplo 2.3. Resolver la ecuacion
(2.27) x 3x + 2x = t2 e5t .
En este caso el polinomio caracterstico 2 3 + 2 tiene races 1 y
2. La solucion general de la ecuacion homogenea es
g (t) = c1 et + c2 e2t .
Ensayemos ahora p (t) = (a + bt + ct2 )e5t . Reemplazando en (2.27),
se observa que es posible simplificar por exp(5t) y se debe tener una
igualdad entre dos polinomios, para que p sea solucion:
12(a + bt + ct2 ) + 7(b + 2ct) + 2c = t2 ,
de modo que los coeficientes deben satisfacer a = 1/72, b = 1/72,
c = 1/12.
Luego, la solucion general de la ecuacion (2.27) es
1 1 1
(2.28) (t) = c1 et + c2 e2t + [ t + t2 ]e5t .
72 72 12
3. Sistemas homogeneos con coeficientes variables
En esta seccion estudiaremos la ecuacion homogenea con coefi-
cientes variables:
(3.1) X = A(t)X,
70 4. ECUACIONES DIFERENCIALES LINEALES
xn (t)
y A(t) es la matriz de coeficientes del sistema:
a11 (t) a12 (t) a1n (t)
A(t) = ..
.
..
.
..
. .
an1 (t) an2 (t) ann (t)
Supondremos en todo lo que sigue que las funciones aij (t) son con-
tinuas, de modo que por la teora que desarrollaremos en el Captulo
7, dado cualquier x0 Cn , el problema con valores iniciales
{
X = A(t)X
(3.2)
X(t0 ) = x0 ,
tiene una unica solucion.
Al igual que en el caso de coeficientes constantes, el conjunto de
todas las soluciones de la ecuacion (3.1) es un espacio vectorial. Mas
aun:
Teorema 3.1. El conjunto de todas las soluciones de la ecuacion
(3.1) es un espacio vectorial de dimension n.
Prueba. Para j = 1, 2, . . . , n sea j (t) la solucion del problema con
valores iniciales
{
X = A(t)X,
(3.3)
X(t0 ) = ej ,
donde
1 0 0
0 1 0
e1 = .
..
, e2 = . , . . . , en = .
.. ..
.
0 0 1
Si 1 , 2 , . . . , n son escalares que verifican
1 1 + 2 2 + + n n = 0,
entonces, evaluando en t0 y usando que {e1 , e2 , . . . , en } es un conjunto
linealmente independiente, obtenemos que 1 = 2 = = n = 0.
Es decir, {1 , 2 , . . . , n } es un conjunto linealmente independiente.
3. SISTEMAS HOMOGENEOS CON COEFICIENTES VARIABLES 71
n
(3.7) ij = aik kj ,
k=1
para cada i = 1, 2, . . . , n.
Usamos ahora una conocida formula que permite calcular la derivada
de un determinante y que puede ser demostrada facilmente usando in-
duccion.
Esta es,
(3.8)
11 12 1n 11 12 1n
d 21 22 2n 21 22 2n
(det (t)) = det
...
+ det
.. +
dt .
n1 n2 nn n1 n2 nn
(3.9)
11 12 1n
21 22 2n
+ det
...
n1 n2 nn
Notemos que la matriz (t) cuyas columnas son 1 (t), 2 (t), . . . , n (t)
es una matriz de soluciones del sistema (3.13).
Definicion 3.3. El Wronskiano de n soluciones x1 (t), x2 (t), . . . , xn (t)
de la ecuacion lineal homogenea (3.11) es:
W (t) = W (x1 (t), x2 (t), . . . , xn (t))
= det (t)
x1 (t) x2 (t) xn (t)
x1 (t) x
(t) xn (t)
2
= det .. .. .. .
. . .
(n1) (n1) (n1)
x1 (t) x2 (t) xn (t)
Gracias al lema anterior, los resultados de esta seccion pueden ser
aplicados al caso que estudiamos, resultando el Teorema que sigue,
cuya demostracion es inmediata.
Teorema 3.3. Sean x1 (t), x2 (t), . . . , xn (t) n solucion de la ecuacion
lineal homogenea (3.11). Entonces, su Wronskiano satisface
t
a1 (s)ds
W (t) = W (t0 )e t0
.
Ademas, estas soluciones son linealmente independientes si y solo si
W (t) = 0, para todo t y esto ultimo equivale a que W (t0 ) = 0, para
algun t0 .
de donde,
t
1
(4.5) X(t) = (t) (t0 )X(t0 ) + (t)1 (s)b(s)ds.
t0
Una vez mas hacemos notar que el calculo de la matriz (t)1 (s)
permite encontrar facilmente una solucion particular p (t) de la ecuacion
no homogenea. Esta solucion satisface ademas la condicion inicial
p (t0 ) = 0.
Definicion 4.1. La funcion de Green para el problema X = A(t)X
es la matriz G(t, s) = (t)1 (s), donde (t) es una matriz fundamen-
tal.
Es posible demostrar que la funcion de Green no depende de la
eleccion particular de la matriz fundamental y satisface
{
d
(4.6) G(t, s) = A(t)G(t, s)
dt
G(s, s) = I.
Observemos tambien, que la solucion general de la ecuacion no
homogenea, dada por la identidad (4.5) es de hecho una version del
metodo de variacion de parametros, puesto que dicha solucion puede
ser escrita como:
t
1
(t) = (t)( (t0 )X(t0 ) + (s)b(s)ds)
t0
= (t)C(t),
donde t
1
C(t) = (t0 )X(t0 ) + 1 (s)b(s)ds
t0
es un vector variable. Ciertamente, cuando C es un vector de constan-
tes, (t) = (t)C es la solucion general de la ecuacion homogenea.
4. SISTEMAS NO HOMOGENEOS CON COEFICIENTES VARIABLES 79
f (t)
Ademas, por el Teorema 4.1, una solucion particular del sistema
(4.9) esta dada por
t
(4.10) p (t) = G(t, s)b(s)ds,
t0
vn
Entonces, la primera componente de G(t, s)b(s) es:
(4.12) x1 (t)v1 (s) + x2 (t)v2 (s) + + vn (t)vn (s).
80 4. ECUACIONES DIFERENCIALES LINEALES
5. Ejercicios propuestos
(1) Encuentre la solucion general de la ecuacion X = AX, donde
la matriz (
de coeficientes
) A esta dada por:
3 2
(a) A = .
1 1
1 5 0
(b) A = 1 3 0 .
0 0 1
1 1 1
(c) A = 2 1 1 .
3 2 4
1 1 2
(d) A = 1 1 1 .
2 1 3
3 2 1
(e) A = 3 2 1 .
3 2 1
(2) En los casos 1b, 1c, 1d , 1e, resuelva el problema de valores
iniciales
X = AX
1
(5.1)
X(0) = 1 .
1
(3) Para cada matriz del ejercicio 1 encuentre exp(At).
82 4. ECUACIONES DIFERENCIALES LINEALES
x = Ax + bu
tiene un polinomio caracterstico
p() = n + an1 n1 + + a1 + a0 ,
entonces, la matriz cambio de base T definida por
a1 a2 a3 an1 1
a2 a3 1 0
[ ] a3
T = b .. Ab .. A2 b .. .. An1 b
. . . .
...
a 1
n1
1 0 0
transforma al sistema en uno de la forma:
0 1 0 0 ... 0 0
0 0 1 0 ... 0 0
b bbU,
x = Ax+ con b=
A
.. .. ..
, bb =
..
. . . .
0 0 0 0 ... 1 0
a0 a1 a2 a3 . . . an1 1
(Indicacion: usar un metodo inductivo).
(b) Dado el sistema:
x1 = 9x1 + 6x2 + 3x3
x2 = 17x1 11x2 5x3
x3 = x1 x3
sin resolverlo, hallar una ecuacion diferencial lineal con
coeficientes constantes que satisfaga x1 .
(19) Se tiene el sistema real lineal:
x = Ax + bu
y = cT x
Donde
8 10 4 1
A = 8 10 5 , b = 1 , c =
7 9 6 1
0
0 ,
1
x1
x = x2
x3
u es una funcion conocida, e y se calcula a partir de x.
(a) Determinar los valores propios del sistema si u 0
5. EJERCICIOS PROPUESTOS 87
k1
(b) Si se hace u = k T x, con k = k2 , determine k tal que
k3
los nuevos
valores propios
del sistema sean 2, 3, 4.
5
(c) Si x(0) = 5 , determine x(t) para los casos a) y
5
b), grafique cada componente en el tiempo y compare lo
obtenido en cada caso.
(20) En el diseno de un controlador, se quiere minimizar la ex-
presion:
J= (xT Qx + ru2 )dt.
0
Siendo Q simetrica positiva definida y r positivo.
(a) Demuestre que si P es solucion de la ecuacion
P A + AT P P b r1 bT P + Q = 0
entonces, el valor de u que minimiza la expresion J es:
u = r1 bT P x.
(b) Dado el sistema x = Ax, con
[ ] [ ]
3/2 1/2 1
A= x(0) = ,
1 0 1
determine y grafique cada componente de x(t).
(c) Dado, el sistema [ ]
1
x = Ax + bu, con A y x(0) de b) y b = , si Q =
0
[ ]
1 0
y r = 4, determine u que minimice J definido en
0 2
a), y calcule y grafique los componentes de x(t). Compare
con (b).
(21) Resolver t3 x + t2 x + tx = 1.
(22) Encontrar una funcion f (x2 + y 2 ) que satisfaga la ecuacion:
2f 2f
+ =0
x2 y 2
con las condiciones iniciales f (1) = 0, f (1) = 1.
(23) Resolver:
(a) x + 7x + 10x = 0
(b) x + 4x + 4x = 0
(c) x + x + x = 0
88 4. ECUACIONES DIFERENCIALES LINEALES
en funcion de T y n.
(27) Dada la ecuacion diferencial parcial:
2f 2f 2 f
2
f f
xy( 2
2
) (x2
y ) +y x =0
x y xy x y
(a) Demostrar que si se hace el cambio
x = r cos
y = r sen ,
entonces,
f
x
= cos f
r
1
r
sen f
(r = 0)
f
y
= sen f
r
+ 1
r
cos f
t
x1 (t)x2 (s) x1 (s)x2 (t)
a(s)ds
x3 (t) = e ds
t0 [x1 (s)x2 (s) x1 (s)x2 (s)]2
para algun t0 .
(c) Hallar la solucion general de la ecuacion:
(t2 t + 1)x (t + t2 )x + (t2 + t + 3)x + (4 + t t2 )x = 0
sabiendo que una solucion es de la forma tk cos(t).
(30) Dadas las ecuaciones:
i) t2 x + tx x = 0
ii) t2 x + tx + x = 0
(a) Resuelva i) sabiendo que una solucion es un polinomio en
t.
(b) Repita a) para la ecuacion ii) Que sucede?
(c) Resuelva ahora i) suponiendo que la solucion de de la
forma xk , con k por determinar.
(d) Utilice la idea dada en c) para resolver ii), y exprese la
solucion como una combinacion lineal de dos soluciones
reales.
(e) Resuelva ahora ambas ecuaciones mediante el cambio de
variable t = eu . Compare con los resultados anteriores.
(31) Exprese, en la forma mas explcita que pueda, la solucion ge-
neral del sistema:
90 4. ECUACIONES DIFERENCIALES LINEALES
1/t 1 t
x = Ax con A = 0 3t2 + 1 0
0 5t 2
(32) Dado el sistema homogeneo de ecuaciones:
x = Ax
(a) Demostrar que si se hace un cambio de base x = P x, los
valores propios del sistema no cambian.
(b) Demostrar que si los valores propios son 1 , 2 , . . . , n ,
(no necesariamente distintos) y A = [aij ], i, j =
1, . . . , n, entonces,
n
n
k a = 0.
k=1 =1
1. Transformacion de Laplace
Sea f una funcion definida sobre [0, [ y con valores en Cd .
Definicion 1.1. Supongamos que la funcion Tf es localmente inte-
grable, vale decir, para todo T > 0, la integral 0 f (t)dt existe. Dado
un complejo = + i, diremos que la transformada de Laplace
de f existe en el punto si la integral impropia
T
t
(1.1) f (t)e dt = lim f (t)et dt,
0 T 0
Como
|(T )eT | Const. e()T
y ademas L() existe si > , de lo anterior se deduce que L ()
existe en el mismo dominio y se tiene, haciendo tender T a ,
(2.6) L () = L() (0).
Hemos probado as el caso k = 1. Supongamos ahora que la
propiedad se cumple para 1 k (r 1) y probemos que es tambien
cierta para k = r. Basta llamar (t) = (r1) y aplicar el primer caso
demostrado: tiene transformada de Laplace y su valor queda dado
por la formula
L () = L() (0)
98 5. LA TRANSFORMACION INTEGRAL DE LAPLACE
vale decir,
M t
|g(t)| (e 1)
y, en consecuencia, g(t) = O(et ).
Enseguida, g es nula en 0 y es derivable, de derivada f . Por lo
tanto, como Lg () = Lg(), se obtiene (2.14).
T
T k
= e et f (t)dt
k=0 0
T t
0
e f (t)dt
= , ( > 0).
1 eT
a)
1
;
2 + 2 + 2
b)
+1
;
2 + 2 + 2
c)
3 + 2
;
2 + 2 + 2
d)
n .
n=1
3. Truncamientos y desplazamientos
En esta seccion estudiaremos la transformacion de Laplace de al-
gunas funciones especiales y, asimismo, el efecto que producen trun-
camientos y desplazamientos sobre la transformada de una funcion
cualquiera. Comencemos por introducir las funciones de Heaviside.
Definicion 3.1. La funcion de Heaviside con salto en 0 se
define como
{
1 si t 0,
(3.1) H(t) =
0 si t < 0.
La funcion de Heaviside con salto en c > 0 es Hc (t) = H(t
c), (t R).
Observacion 3.1. Calculo de transformadas de Laplace de H y
Hc .
Un calculo directo nos da
LH() = et dt
0
1
(3.2) = , ( > 0).
104 5. LA TRANSFORMACION INTEGRAL DE LAPLACE
= f (s)ec es ds
0
c
= e Lf ().
4. Producto de convolucion
Definicion 4.1. Dadas dos funciones f y g integrables sobre [0, [,
definimos su producto de convolucion f g por la expresion:
t
(4.1) f g(t) = f (s)g(t s)ds, (t 0).
0
5. Ejercicios propuestos
(1) Dado el sistema:
x = Ax , x(0) = x0
(a) Demostrar que su solucion se puede expresar de la forma:
x = L1 [(I A)1 ] x0
(b) Demostrar que L[eAt ] = (I A)1
(c) Usar lo anterior para
1 0 0 1
A = 0 0 1 , x0 = 1
0 1 0 1
(2) Resolver la ecuacion:
N
x 2x + x = ent
n=1
con x(0) = 0, x (0) = 1
(3) Dada la ecuacion
x(n) +an1 x(n1) + +a1 x +a0 x = bm u(m) +bm1 u(m1) + +b1 u +b0 u
(a) Demostrar que si u(t) = beat , una solucion particular es
de la forma xp (t) = Aeat . Determinar A y la condicion
para que esto se cumpla.
108 5. LA TRANSFORMACION INTEGRAL DE LAPLACE
(1)n (t/2)2n
es x(t) =
n=0
(n!)2
[Indicacion: usar serie binomial cuando sea pertinente.]
5. EJERCICIOS PROPUESTOS 109
{
t2
f (t0 ) t1 t0 t2
f (t)(t t0 )dt =
t1 0 en caso contrario,
para toda f continua en t0 , entonces se tendra L[] = 1.
Hacemos notar que (t) no es una funcion; es mas bien
una distribucion. Su construccion rigurosa utiliza la Teora
de Distribuciones, al interior de la cual se puede justificar la
propiedad anterior. Esto permite que dicha propiedad sea uti-
lizada en la resolucion formal de algunas ecuaciones diferen-
ciales que aparecen muy corrientemente en Ingeniera.
(10) (a) Usando el ejercicio anterior, resolver formalmente la ecuacion
diferencial:
x(iv) = K(t L), (K constante)
con las condiciones de borde:
x(0) = 0, x (0) = 0, x (2L) = 0, x (2L) = 0
110 5. LA TRANSFORMACION INTEGRAL DE LAPLACE
1
(0.7) a2(m+1) = (1)m+1 , (m 0).
(m + 1)!
Finalmente, la solucion buscada es
t2m
(0.8) (t) = (1)m .
m=0
m!
1. Funciones analticas
Definicion 1.1. Una funcion f definida en un dominio D(f ) de
R y con valores en Cd , es analtica en un punto t0 D(f ) si existe una
sucesion de coeficientes (an )n en Cd y un intervalo abierto alrededor de
t0 , tal que para cada t en dicho intervalo se cumpla la igualdad
(1.1) f (t) = an (t t0 )n .
n=0
2. Puntos ordinarios
En esta seccion comenzamos por estudiar las ecuaciones de la forma:
(2.1) X = A(t)X, (t ], [),
donde A(t) es una funcion analtica sobre todo el intervalo ], [ con
valores matrices de d d. Decimos que en cada punto del intervalo de
analiticidad es un punto ordinario.
Teorema 2.1. Bajo las hipotesis anteriores, dado un punto ordi-
nario t0 ], [, toda solucion de la ecuacion (2.1) se puede expresar
en la forma de una serie de potencias:
(2.2) (t) = ck (t t0 )k ,
k=0
1
k
(2.4) ck+1 = Ar ckr , (0 k d).
k + 1 r=0
1 k
|ck+1 | M |ckr |,
k + 1 r=0
2. PUNTOS ORDINARIOS 115
( 1)( 3) ( + 1 2k)
(2.20) c2k+1 = (1)k 2k c1 , (k 1).
(2k + 1)!
118 6. SERIES DE POTENCIAS Y ECUACIONES DIFERENCIALES
s = s(t) = log(t ),
y
( )
0 1
(3.16) B(t) = ,
b2 (t) 1 b1 (t)
es una funcion analtica con valores matriciales, de modo que la ecuacion
de Euler general es tambien un caso particular de ecuaciones del tipo
(3.8) y (3.9).
Uno de los resultados fundamentales de esta seccion es el siguiente
teorema debido a Frobenius y Fuchs.
Teorema 3.1. Supongamos que B(t) es analtica para t <
y que los valores propios 1 y 2 de B() no son iguales ni difieren por
un entero. Entonces la ecuacion
1
(3.17) X = B(t)X, t <
t
tiene dos soluciones de la forma
(3.18) 1 (t) = (t ) 1
ck (t )k ,
k=0
(3.19) 2 (t) = (t )2 dk (t )k
k=0
y son linealmente independientes sobre t < . Los coeficientes ck
y dk pueden ser evaluados por sustitucion.
Prueba. Comenzamos por hacer la sustitucion X = (t )1 Y en la
ecuacion original, que se transforma en
(3.20) (t )Y = A0 Y + Ak Y (t )k ,
k=1
donde
A0 = B0 1 I, Ak = Bk , para k 1, B(t) = Bk (t )k .
k=0
1
k1
1
(3.23) ck = (A0 kI) Akm ck .
k m=0
Pero esta desigualdad tiene la misma forma que aquella (2.5) usada en
la demostracion del Teorema 2.1. Repitiendo el mismo argumento
de-
sarrollado en ese teorema se concluye entonces que Y1 (t) = k=0 k
c (t
) converge para |t | < ya que L puede ser arbi-
k
4. La ecuacion de Bessel
La Trigonometra puede ser vista como el estudio de las soluciones
de la ecuacion del oscilador armonico
x + 2 x = 0.
De manera similar, la teora de las Funciones de Bessel se reduce
al estudio de las soluciones de la ecuacion
(4.1) t2 x + tx + (t2 2 )x = 0,
donde es una constante. Esta ecuacion se conoce como la Ecuacion
de Bessel de parametro . Tiene un punto singular regular en t = 0
y la teora desarrollada hasta aqu se aplica en el intervalo ]0, [.
Supongamos por simplicidad que es real y 0. La ecuacion
indicial es en este caso
(4.2) 2 2 = 0,
vale decir, sus races son 1 = , 2 = . Luego 1 2 = 2 y
tendremos dos casos notables: el primero, si 2 es un numero entero;
el segundo, si lo anterior no ocurre.
Caso 2 Z :
En este caso hay dos soluciones de la ecuacion de Bessel, linealmente
independientes, que denotamos
(4.3) (t) = ck tk+ ,
k=0
126 6. SERIES DE POTENCIAS Y ECUACIONES DIFERENCIALES
(4.4) (t) = dk tk .
k=0
= (2 + 1)c1 + {[( + k)2 2 ]ck + ck2 }t+k = 0.
k=2
Dado que = 1/2, se tiene c1 = 0. Los otros coeficientes satisfacen
la ecuacion recursiva
ck2
(4.5) ck = , (k 2).
( + k)2 2
Luego, ck = 0 si k es impar. Si k es par, de la forma k = 2m, la
relacion anterior conduce a
(1)m c0
c2m = 2m , (m 1).
2 m!( + 1) . . . ( + m)
Para simplificar la escritura de los coeficientes recurramos a la
funcion Gama, (), definida para > 0 como sigue:
(4.6) () = t1 et dt.
0
(4.7) ( + 1) = (),
que extiende la propiedad de los factoriales. En efecto, de la propiedad
precedente se tiene que (n + 1) = n! sobre los enteros. Asimismo,
la ecuacion (4.7) hace posible tabular la funcion gama para > 1 sin
tener que evaluar la integral de (4.6). En el mismo espritu, uno define
la funcion gama como (1/)( + 1) para negativo y diferente de
todo entero. Se sabe, por ejemplo, que (1/2)
= . De esto se
deduce que (3/2) = /2 y (1/2) = 2 .
Con la funcion gama podemos entones escribir los productos de la
forma:
( + 1) . . . ( + k),
4. LA ECUACION DE BESSEL 127
como
( + k + 1)
.
()
De este modo, si se escoge el coeficiente c0 de segun
1
c0 = ,
2 ( + 1)
la correspondiente solucion es denotada por el smbolo J y es cono-
cida con el nombre de Funcion de Bessel de primera especie de ndice
. Se tiene as:
( ) ( )2k
t (1)k t
(4.8) J (t) = .
2 k=0
k!( + k + 1) 2
Y si se considera la otra raz de la ecuacion indicial, , se llega a
la funcion de Bessel de ndice :
( ) ( )2k
t (1)k t
(4.9) J (t) = .
2 k=0
k!( + k + 1) 2
La solucion general de la ecuacion de Bessel planteada tiene en-
tonces la forma
( )
J (t)
J+1 (t) = 2 J1 (t),
t
para todo = 0, = 1. La existencia de identidades entre las funciones
de Bessel permite tabularlas a partir de los valores de J (t) y J (t)
para 0 1.
Cuando 2 no es entero ni un multiplo entero impar de 1/2, en-
tonces mismo debe ser un entero. El corolario 3.2 se aplica entonces
en toda generalidad y existe una solucion de la forma (4.10) con
c = 0, linealmente independiente de J . Reemplazando dicha forma de
solucion en la ecuacion de Bessel para calcular los coeficientes, obten-
emos:
2cJ (t) + [( 1) ]d1 t
2 2 1
+ {[(k )2 2 ]dk + dk2 }tk = 0.
k=2
= (1 2)d1 t + {[(k )2 2 ]dk + dk2 }tk .
k=2
5. Ejercicios propuestos
(1) Utilizando el metodo de series de potencias, encuentre la solucion
general de las siguientes ecuaciones lineales:
(a) (1 + t2 )x 6tx + 12x = 0
(b) 2(1 + t2 )x 2tx 6x = 0
(c) 6(3 + t2 )x + tx + x = 0
(d) (1 + 21 t2 )x + 81 tx + 16
1
x=0
2
(e) (2 + 7t )x + 35tx + 28x = 0
(f) (3 + 5t2 )x + 20tx + 10x = 0
(g) (2 + 3t2 )x 15tx + 27x = 0
Determine, en cada caso, el radio de convergencia.
(2) En la ecuacion 4x + (4k t2 + 2)x = 0, con k constante:
(a) Demuestre que los coeficientes de la serie estan relaciona-
dos por una recurrencia de tres terminos.
(b) Demuestre que si se hace el cambio de variables x = y
et /4 , la recurrencia tiene ahora dos terminos.
2
(d) t3 x + (sent)x = 0
(e) x + tntx + tx = 0
(f) costx + x + 2x = 0
(4) Dada la ecuacion de Bessel con p = 0:
tx + x + tx = 0
(a) Recordando que la Transformada de Laplace tiene la propiedad
de que:
d
L(t x(t)) = L(x(t)),
ds
demuestre que una solucion de la ecuacion satisface:
1
L(x(t)) =
2
s +1
(b) Usando el metodo de serie de Frobenius, demuestre que
una solucion es:
(1)n (t/2)2n
x(t) =
n=0
(n!)2
(1)n 2n n1 (p 2k 1)
x2 = t + k=0
t2n+1
n=1
(2n + 1)!
(b) Demuestre que si p es un entero, una de las dos soluciones
es un polinomio. Determine estos polinomios para p =
0, 1, 2, 3, 4, 5.
(c) Cuando el termino que contiene a la potencia mas alta
de t es de la forma 2n tn , se conoce como polinomio de
Hermite, Hn (t) y se sabe que estan dados por la formula:
2 dn t2
Hn (t) = (1)n et e
dtn
5. EJERCICIOS PROPUESTOS 133
(1 + t)p = H(p, b, b, t)
n(1 + t) = t H(1, 1, 2, t)
1 1 3
Arcsen t = t H( , , , t2 )
2 2 2
t
et = lim H(a, b, a, )
b b
(e) Demuestre que la solucion general de la ecuacion de Tcheby-
shev:
(1 t2 )x tx + p2 x = 0
se puede expresar, mediante el reemplazo
1
u = (1 t)
2
134 6. SERIES DE POTENCIAS Y ECUACIONES DIFERENCIALES
como
1 1t 1t 1 1 3 1t
x = c1 H(p, p, , ) + c2 ( )H(p + , p + , , )
2 2 2 2 2 2 2
(f) Demuestre que
d ab
H(a, b, c, t) = H(a + 1, b + 1, c + 1, t)
dt c
(g) Demuestre que en la ecuacion de Legendre
(1 t2 )x 2tx + n(n + 1)x = 0
la solucion polinomica, para un entero positivo n dado, es
1
Pn (t) = H(n, n + 1, 1, (1 t))
2
y que esta se puede expresar como
[n/2]
(1)k (2n 2k)!
Pn (t) = tn2k
k=0
2n k!(n k)!(n 2k)!
donde [n/2] es el mayor entero menor o igual a n/2.
(12) En el problema anterior se mencionaron las soluciones de la
ecuacion de Legendre conocidas como Polinomios de Legendre,
Pn (t). Estos polinomios pueden expresarse en forma mas com-
pacta como:
1 dn 2
Pn (t) = (t 1)n
2n n! dtn
(a) Demuestre la propiedad de ortogonalidad, que dice que:
{
1
0 m = n
Pm (t)Pn (t)dt = 2
1 2n+1
m=n
(b) Demuestre que otra forma de escribir la ecuacion de le-
gendre es:
d
((1 t2 )x ) + n(n + 1)x = 0
dt
(c) Utilice el resultado de b) para demostrar que:
1
Pm (t)Pn (t)dt = 0 si m = n
1
x2 x1 1 u1
1+a x3 x2 1 u2
b = (AT A)1 AT .. con A= .. .. ..
c . . . .
xn xN 1 1 uN 1
(e) Obtenga un modelo dinamico lineal de primer orden, del
tipo
x = ax + b + cu
con t = 1, para los siguienes datos:
0.80 1
0.90 0.6
0.99 1.8
0.91 1.5
M =
1.29 2.5
1.11 2
1.41 3.1
0.75 0.5
(f) Usando el modelo discreto, haga una simulacion del mod-
elo a partir de los datos de u, y compare con los datos de
x.
CAPTULO 7
es una funcion (t) que satisface la ecuacion 1.1 y que ademas verifica,
(i) t0 I = Dom();
(ii) (t0 ) = x0 .
El concepto de solucion incluye implcitamente la existencia de un
intervalo en el cual esta verifica la ecuacion. En 1.3 hemos definido los
conceptos de extension de soluciones y el de solucion maximal. Cuando
el intervalo de definicion de una solucion es toda la recta real, ella es
obviamente maximal. En ese caso decimos que la solucion es global.
Aqu, demostraremos solo criterios de existencia local de soluciones,
cuyo dominio I es posiblemente pequeno.
Ejemplo 1.1. La funcion (t) = 1t , (t D() =]0, [) es una
solucion (de hecho es la unica solucion global) del problema con valores
iniciales
{
x (t) = x2 (t),
(1.3)
x(1) = 1.
En este caso, f (t, x) = x2 es una funcion infinitamente diferenciable
en todo R2 , pero la solucion solo tiene sentido para t ]0, [.
El siguiente es un resultado sobre existencia de soluciones locales.
Se conoce como el Teorema de Cauchy-Peano y su demostracion puede
encontrarse por ejemplo en el libro de F. Simmons y esta basada en
una aproximacion de la solucion mediante curvas poligonales. En este
teorema, D Rn+1 es una region dada por :
D = I R
= {(t, x) : t I, x = (x1 , .., xn ), ai xi bi }
donde I = (a, b) es un intervalo en R.
Teorema 1.1. Sea (t0 , x0 ) D y sea f : D Rn una funcion
continua. Entonces, existe > 0 tal que el problema con valor inicial:
{
x (t) = f (t, x),
(1.4)
x(t0 ) = x0 ,
tiene una solucion (t) con dominio en el intervalo ]x0 , x0 + [.
Hacemos notar que la solucion (t) es local y de hecho el intervalo
]x0 , x0 + [ podra ser muy pequeno. Por otra parte, el teorema solo
garantiza la existencia de una solucion. El ejemplo que sigue muestra
que esta, en general, no es unica, aun cuando satisface una condicion
inicial dada.
1. EXISTENCIA Y UNICIDAD DE SOLUCIONES LOCALES 141
{
0, si 0 t
(1.5) (t) =
( 32 (x ))3/2 , si t 1,
{
x (t) = x1/3
(1.6)
x(0) = 0,
t
(1.7) x(t) = x0 + f (s, x(s))ds.
t0
(1.8) M = {v : [, ] Rn : v es continua}.
Dado que una funcion continua definida en un intervalo cerrado y
acotado es automaticamente acotada, podemos usar la expresion:
d(u, v) = 0 si y solo si u = v,
d(u, v) = d(v, u)
t
|T u(t) T v(t)| = | (f (s, u(s)) f (s, v(s))ds|
0
t
|f (s, u(s)) f (s, v(s))|ds
0
t
c|u(s) v(s)|ds
0
t
cd(u, v) ds
0
cd(u, v),
para t 0. Para t 0, esta desigualdad se demuestra de manera
analoga, de modo que obtenemos:
d(T u, T v) cd(u, v)
Ahora, escogemos suficientemente pequeno de modo que: =
c < 1. Entonces,
para todo t [, ].
En otras palabras, tomamos 0 la funcion con valor constante x0 , 1 =
T 0 , 2 = T 1 y, en general, n+1 = T n .
La desigualdad (1.12) da de inmediato que:
d(n+1 , n ) d(n , n1 )
Ademas, usando recursivamente esta relacion, obtenemos:
d(n+1 , n ) n d(1 , 0 ).
Por otra parte,
t
|1 (t) 0 (t)| = | f (s, x0 )ds|
0
t
|f (s, x0 )|ds
0
L,
(1.13) d(n+1 , n ) n L
Por lo tanto,
Analogamente,
d(n+k , n ) (1 + + 2 + k1 )L.
1. EXISTENCIA Y UNICIDAD DE SOLUCIONES LOCALES 145
Pero,
1 + + 2 + k1 1 + + 2 +
= k
k=0
1
= ,
1
donde hemos usado que la serie geometrica converge puesto que
< 1.
Por lo tanto, para todo n, k numeros naturales hemos demostrado
que la sucesion n (t) verifica,
L
(1.14) d(n+k , n ) n1 .
1
Esto demuestra que (n (t)) es una sucesion de Cauchy y por lo
tanto converge a una funcion (t) M .
Solo resta verificar que la funcion lmite (t) es el unico punto fijo
de T .
Tomando lmite cuando k tiende a infinito en (1.14), obtenemos
primero,
L
(1.15) d(, n ) n1 .
1
Usando esta relacion y la estimacion (1.12) sigue que,
d(, ) = d(T , T )
d(, ),
de donde
0 (1 )d(, ) 0
y como < 1 obtenemos que d(, ) = 0, o sea = .
(c)n L
(1.18) |n (t) (t)| ,
1 c
donde L es el supremo de f (t, x).
Esta es una expresion explcita del error cometido al aproximar la
solucion (t) por la n-esima aproximacion n (t). Para estimarla, basta
conocer la constante de Lipschitz c y el supremo L de la funcion f (t, x).
El radio del intervalo de existencia queda determinado por la relacion
c < 1.
3. Ejercicios propuestos
(1) Considerar el problema con valores iniciales
x = y
y = sen(xy)
(3.1)
x(0) = 0
y(0) = 0.
que puede ser escrito en representacion vectorial en la forma
X = F (X),
donde
( ) ( )
x y
F (X) = F = .
y sen(xy)
a) Pruebe que la funcion F satisface una condicion de Lips-
chitz. Deduzca que existe una solucion unica del problema
en [, ] para algun
> 0.
b) Pruebe que C = 6 es una constante de Lipschitz apropi-
ada y que existe una
unica solucion en el intervalo [1/ 6, 1/ 6].
c) Escoja = 1/10 6 y determine una cota para el error
cometido al usar la aproximacion de Picard con n itera-
ciones.
d) Determine el numero n de iteraciones necesarias para tener
un error 0, 01 al usar el metodo de aproximacion de Pi-
card.
3. EJERCICIOS PROPUESTOS 151
X = AX,
donde, ( )
a b
A=
c d
El polinomio caracterstico de la matriz de coeficientes A, puede ser
escrito
L() = 2 tr(A) + det(A).
El discriminante vale entonces:
= (tr(A))2 4det(A)
Caso 1. Los valores propios de A son reales, distintos y tienen el
mismo signo. Esto ocurre cuando (tr(A))2 4det(A) y det(A) > 0
En este caso, existen vectores propios linealmente independientes
v1 , v2 con Av1 = 1 v1 , Av2 = 2 v2 . Ademas, la solucion general tiene
la forma,
(t) = c1 eAt v1 + c2 eAt v2
= c1 e1 t v1 + c2 e2 v2 ,
con c1 y c2 constantes arbitrarias.
El comportamiento asintotico de la solucion general, para valores
grandes de t, depende directamente del signo de los valores propios.
El origen es nodo estable o atractor o foco: 1 , 2 < 0.
En este caso, (t) 0, cuando t y, por lo tanto, toda
solucion es asintoticamente estable, para t positivo.
2. ESTABILIDAD DE SISTEMAS LINEALES HOMOGENEOS EN DIMENSION157
2
=0
Espirales Espirales
Centros
<0,tr=0
Fo
co
s
s
Nodos co Nodos
Fo >0, tr>0
>0, tr<0
tr (A)
Sillas
3. ESTABILIDAD PARA SISTEMAS LINEALES 163
4. El metodo de Liapunov
Presentaremos en lo que sigue, el metodo de Liapunov para analizar
la estabilidad de un punto de equilibrio. Este tiene como base la idea
intuitiva de que cuando la energa de una partcula tiene un valor
mnimo, entonces ella esta en un estado estable. Consideremos el sis-
tema autonomo
{
x1 = f1 (x1 , x2 )
(4.1) ,
x2 = f2 (x1 , x2 )
donde f = (f1 , f2 ) es una funcion definida en una region D R R,
con valores reales.
Suponemos que el origen (0, 0) es un punto crtico del sistema, que
esta en el interior de D.
Una funcion de Liapunov para este sistema , es una funcion E(x1 , x2 ),
de clase C 1 , definida en D R R que verifica,
(i) E(0, 0) = 0 y E(x1 , x2 ) > 0 , para todo (x1 , x2 ) D con
(x1 , x2 ) = (0, 0).
(ii) la funcion F (x1 , x2 ) = x
E
1
E
f1 + x 2
f2 satisface F (0, 0) = 0
y F (x1 , x2 ) 0 , para todo (x1 , x2 ) D con (x1 , x2 ) = (0, 0).
Una funcion E(x1 , x2 ) que satisface la propiedad (i) se llama definida
positiva o de tipo positivo, en cambio una que verifica (ii) se llama neg-
ativa semipositiva o de tipo seminegativo. En forma similar, se definen
las funciones negativa definida y positiva semidefinida.
Intuitivamente, si E(x1 , x2 ) es positiva definida, entonces el grafico,
en el espacio R3 , de la superficie x3 = E(x1 , x2 ) se parece a un paraboloide
que se abre hacia arriba del plano (x1 , x2 ), con vertice en el origen.
La importancia de una funcion de Liapunov para el sistema de ecua-
ciones, proviene del hecho que si (t) = (x1 (t), x2 (t)) es una solucion del
mismo, entonces , a lo largo de la trayectoria (x1 (t), x2 (t)), la funcion
e(t) = E(x1 (t), x2 (t)) tiene derivada,
d E E
E(x1 (t), x2 (t)) = f1 + f2 ,
dt x1 x2
que es siempre negativa semidefinida, por lo que e(t) es decreciente.
Usando este hecho, se puede demostrar el siguiente resultado.
Teorema 4.1. Si existe una funcion de Liapunov para el sistema
(fsisnoli), entonces el origen (0, 0) es estable.
E E
Si ademas la funcion F (x1 , x2 ) = x 1
f1 + x 2
es negativa definida,
entonces el origen es asintoticamente estable.
5. LINEALIZACION 165
5. Linealizacion
Consideremos nuevamente el sistema nolineal
166 8. ESTUDIO CUALITATIVO DE ECUACIONES DIFERENCIALES
{
x1 = f1 (x1 , x2 )
(5.1)
x2 = f2 (x1 , x2 ),
donde F (x1 , x2 ) = (f1 (x1 , x2 ), f2 (x1 , x2 )) es una funcion diferenciable
que satisface F (0, 0) = (0, 0)
Cuando la funcion F (x1 , x2 ) ( o cada una de sus coordenadas) es
analtica, entonces el sistema pude ser escrito en la forma,
{
x1 = ax1 + bx2 + a2 x21 + a3 x22 + a4 x1 x2 +
(5.2)
x2 = cx1 + dx2 + b2 x21 + b3 x22 + b4 x1 x2 + ,
donde las series de potencias en las variables (x1 , x2 ) no tienen terminos
constantes, puesto que F (0, 0) = (0, 0).
Consideramos ahora el sistema lineal asociado,
{
x1 = ax1 + bx2
(5.3)
x2 = cx1 + dx2 ,
para el cual, de acuerdo a un resultado anterior, tenemos un respuesta
completa al estudio de la estabilidad de la solucion trivial (0, 0), en
terminos de los valores propios de la matriz de coeficientes,
( )
a b
A=
c d
Notemos que la diferencia entre ambos sistemas, vale decir, los
terminos nolineales, es grande, pero, para valores de (x1 , x2 ) cercanos
al origen (0, 0), ella es despreciable. Por esta razon, resulta razonable
pensar que la estabilidad de la solucion trivial para la ecuacion no li-
neal dada, este relacionada con su estabilidad para la ecuacion lineal
asociada. Esta ultima depende solo de los valores propios de la matriz
A.
El teorema que sigue establece esencialmente este hecho. En su
enunciado usaremos la notacion
f1 (x) = ax1 + bx2 + a2 x21 + a3 x22 + a4 x1 x2 + = ax1 + bx2 + g1 (x),
f2 (x) = cx1 + dx2 + b2 x21 + b3 x22 + b4 x1 x2 + = cx1 + dx2 + g2 (x).
Es decir, la funcion G(x) = (g1 (x), g2 (x)), con x = (x1 , x2 ), representa
a los terminos nolineales.
Teorema 5.1. Supongamos que la funcion G(x) ||x||
es continua y se
anula en el origen. Entonces, la solucion trivial (0, 0) de la ecuacion
nolineal es asintoticamente estable si la solucion trivial de la ecuacion
lineal asociada tambien es asintoticamente estable.
5. LINEALIZACION 167
para cada instante t para el cual ||(t)|| < r. As, si el dato inicial es
suficientemente pequeno , entonces la norma de la solucion decrece en
el tiempo. Por lo tanto, la solucion trivial es estable. Por otra parte,
de las desigualdades anteriores tambien se obtiene que,
a
||(t)||2 + ||(t)||2 0,
2
de donde se concluye facilmente que la norma de la solucon es expo-
nencialmente decreciente y, por lo tanto, la solucon trivial es, de hecho,
asintoticamente estable.
6. Ejercicios propuestos
(1) Recordando que un punto de equilibrio es una solucion cons-
tante de una ecuacion, encuentre los puntos de equilibiro de
las siguientes ecuaciones y sistemas:
(a) x + 5x + 6x = 1
(b) x1 = sen(x2 )
x2 = ex1 x2 1
6. EJERCICIOS PROPUESTOS 169
x1 + x2
(c) = x1
x1 + x2
x1 x2
= x2
x1 + x2
(d) (x + 1)ex = x 1
(2) Determine si el punto de equilibrio de los sistemas siguientes
es estable o inestable:
(a) x = 4x 6y
y = x + y
(b) x = 2x + y
y = 6x 3y
(c) x = 21 x
y = 41 x + 12 y
(d) x = 7x + 12y + 12z
y = x + 3y + 3z
z = 7x 13y 13z
(e) x = 7y + 13z
y = 4y + 8z
z = 3y 6z
(f) x = x + y
y = y + z
z = 3x 9y 7z 5w
w = 4x + 8y + 8z + 4w
(g) x = 21 x 12 z 12 w
y = 25 x y + 72 z + 32 w
z = x + z + w
w = 27 x + y 32 z 12 w
170 8. ESTUDIO CUALITATIVO DE ECUACIONES DIFERENCIALES
(e) x = x y + z + 1
y = x + y
z = 3x 4y 3z
(f) x = 3x 10y x2
y = x + y + xy 4
(g) x = 10x 3y + y 3 /x
y = 21x + 6y 2x/y
(5) En el caso de sistemas no lineales, del tipo x = F (x), tambien
se puede hacer un analisis de los puntos de equilibrio. Si x0
es tal que F (x0 ) = 0, se puede hacer la traslacion x = x x0
de modo que el nuevo sistema tendra en x = 0 un punto de
equilibrio.
(a) Demuestre que si se define G(x) = F (x + x0 ), entonces
G(x) = G (0)x + 0(x)
con G matriz Jacobiana de G, y 0(x) tal que:
0(x)
lim = 0.
x0 x
( + )x = x + y, (, C, x V );
(3) x = x, ( C, x V ) y ademas 1x = x, 0x = 0, (x V ).
Un subconjunto W de V es un subespacio vectorial (o simple-
mente subespacio) de V , si es un espacio vectorial con las operaciones
de V convenientemente restringidas. Al respecto, notese que W V es
173
174 9. APENDICE: COMPLEMENTOS DE ALGEBRA LINEAL
2. Algebra de Matrices
Las operaciones sobre matrices se definen a partir de aquellas cor-
respondientes a las transformaciones lineales.
x1
X = ...
xn
2. ALGEBRA DE MATRICES 179
Notese que (2.9) puede obtenerse de una manera mas directa usando
(2.1) y (2.7). En efecto, dado cualquier x V , cuyas coordenadas
expresamos por columnas X o X segun la base BV o BV respectiva-
mente, tenemos que si Y designa la columna de coordenadas de Lx en
W , entonces,
Y = AX = A X
pero, X = P X , luego AP X = A X de donde resulta (2.9).
Apliquemos el ultimo metodo anterior para estudiar lo que ocurre
si cambiamos la base de W , sin cambiar aquella de V : entonces las
coordenadas Y de Lx se cambian en Y segun la relacion:
Y = QY .
En consecuencia, si llamamos A la matriz de la representacion de L
segun (BV , BW ), resulta Y = A X y
AX = Y = QY = QA X,
de donde
(2.10) A = Q1 A.
Finalmente, si cambiamos ambas bases (en V y en W ), obtenemos
que la nueva matriz que representa a L es
(2.11) Q1 AP.
En el caso particular de las transformaciones lineales de V en V
(endomorfismos), tenemos que si se cambia la base del espacio segun
una matriz de transferencia P , dada una representacion matricial A de
L L(V, V ) en la primera base, ella se transforma en P 1 AP . Esto
justifica la definicion siguiente:
Definicion 2.2. Dos matrices A y B de n n son equivalentes si
existe una matriz de transferencia P , tal que
B = P 1 AP.
Segun lo anterior, esto significa que ambas matrices representan la
misma transformacion lineal.
3. Funciones de Matrices
Consideremos nuevamente el espacio Md de las matrices comple-
jas de d d. Sobre este espacio vectorial hemos introducido ademas
una operacion de multiplicacion que lo transforma en un algebra no
conmutativa.
Una funcion de variable real con valores matriciales es sim-
plemente una aplicacion M : R Md .
3. FUNCIONES DE MATRICES 183
4. Reduccion de Matrices
Sea V un espacio vectorial de dimension d y consideremos T
L(V, V ). Decimos que un subespacio W de V es invariante por T si
T x W para todo x W (vale decir, T (W ) W ).
Dados dos subespacios de V , sean E y F , la suma directa de ellos,
que escribimos E F es el espacio generado por E F .
Si V se puede escribir como V = W1 . . . Wq donde cada Wi
es invariante para T , entonces podemos considerar las restricciones
Ti = T |Wi de T sobre cada espacio y se escribe:
T = T1 . . . Tq .
Veremos enseguida el caso mas simple de subespacio invariante.
4.1. Valores y Vectores propios. Supongamos en primer lugar
que T admite un subespacio invariante de dimension 1. Esto quiere
decir que existe un vector x V , tal que
(4.1) T x = x,
para algun escalar C. El vector x recibe entonces el nombre de
vector propio y el escalar , es un valor propio.
Si existen d escalares 1 , . . . , d que sean valores propios y para
los cuales se encuentren vectores propios respectivos v1 , . . . , vd que,
ademas, constituyan una base de V , entonces
V = W1 . . . Wd ,
donde Wi es el espacio unidimensional engendrado por vi , (i = 1, . . . , d).
En este caso la aplicacion T se escribe:
(4.2) T = 1 I . . . d I.
Por ende, si la representamos matricialmente segun la base (v1 , . . . , vd )
su matriz D es de tipo diagonal:
1 0 0 ...
0 2 0 . . .
(4.3) D=
.
...
0 0 . . . d
El procedimiento de diagonalizacion de una matriz consiste en
la busqueda de d valores propios y de respectivos d vectores propios,
linealmente independientes con los cuales se pueda construir una base
de V . As por ejemplo, si T tiene una representacion matricial A segun
188 9. APENDICE: COMPLEMENTOS DE ALGEBRA LINEAL
donde los valores i son sus races sobre C y los numeros mi las respec-
tivas multiplicidades. Designando por W (i ) el nucleo de (T i I)mi ,
el espacio V puede descomponerse entonces en la suma
(4.13) V = W (1 ) . . . W (q ).
La restriccion de T i I a cada espacio W (i ) es nilpotente, ((T
i I)mi = 0), y por ende puede representarse por una matriz nilpotente
de mi mi . En consecuencia, la restriccion de T a tales subespacios
queda representada por un bloque de Jordan de la forma Ji .
Existe entonces una matriz de cambio de base P , tal que la matriz
A verifique
A = P 1 JP,
donde J es una matriz que, escrita en bloques, es de la forma
J1 0 ...
0 J2 . . .
(4.14) J =
.
...
0 0 . . . J q
Una aplicacion importante del teorema anterior es que la exponen-
cial de cualquier matriz A se puede calcular en la forma
(4.15)
exp J1 0 ...
0 exp J2 ...
exp A = P 1 (exp J)P, exp J =
.
...
0 ... exp Jq
Ejercicio 5.2. Calcular
eA en los casos siguientes:
0 1 0
(a) A = 0 0 1
0 0 0
0 1 2 4
0 0 3 5
(b) A =
0
0 0 6
0 0 0 0
6. Ejercicios propuestos
(1) Considerar el espacio P OL4 de los polinomios de tercer grado
en la variable x, con coeficientes complejos.
(a) Encontrar una base de P OL4 y demuestre que tiene di-
mension 4.
(b) Representar mediante una matriz la aplicacion lineal T
que a cada polinomio p(x) asocia el polinomio p(x + 1).
(c) Representar mediante una matriz la aplicacion D que al
polinomio p(x) le asocia su derivada.
(d) Obtener finalmente las representaciones de la suma y el
producto de T y D.
(2) Encontrar, si es posible, las matrices B tales que AB = C en
los dos casos siguientes:
6. EJERCICIOS PROPUESTOS 193
2 3 2 1 2 1
(a) A = 4 1 2 C = 2 1 1
0 1 0 1 2 0
2 1 1 1 2 1
(b) A = 8 1 5 C = 2 1 1
4 3 3 5 0 3
(3) Designamos por b un parametro real y se considera la funcion
matricial
( )
cosh t bsenh t
Ab (t) = senh t .
b
cosh t
(a) Encontrar los valores propios de Ab (t) y sus vectores pro-
pios.
(b) Determinar las matrices invertibles P para las cuales B =
P 1 Ab (t)P es una matriz diagonal. Probar que B no de-
pende del parametro b.
(c) Se deja fijo un valor de b. Pruebe que los distintos valores
Ab (t), al hacer variar t R, forman un grupo para el
producto de matrices.
(d) Sea la matriz
( )
M= ,
donde los coeficientes son numeros complejos que verifican
2 = 0.
Utilizar las preguntas anteriores para calcular M n .
(4) Suponga que la matriz A satisface
(6.1) A2 + A = I.
Encuentre exp(At).
(5) Suponga que la matriz A satisface
(6.2) A2 = 2A.
Encuentre exp(At).
(6) Verifique que ( )
0 1
A= ,
1 0
satisface A2 = I y use este hecho para calcular exp(At).
Resultados de ejercicios seleccionados
Cap. II
(2) i) y = 0
ii) y y 6y = 0
iii) xyy + x(y )2 + yy (1 2x) y 2 = 0
iv) x + 4x = 0
v) x 9x = 0
vi) x 20x + 100x = 0
t2
(6) i) x = 70 3 7
t 12
1 6
t + 607 5
t (2 + t2 )n t 2
+ c1
2
+ c2 t + c3
n1
n t
t
ii) x = 3e + 4(1) e + ck tk
k=0
n6
(8) t = 30 n2 77.54
Cap. III
x
(1) i) e t = Cx
Cet
ii) x = 1+Ce t
Cap. IV
cosh( 2t)
2senh( 2t)
(3) a) eAt = et
1
senh( 2t) cosh( 2t)
2
195
196 RESULTADOS DE EJERCICIOS SELECCIONADOS
et (cost + 2sent) 5et sent 0
b) eAt = et sen t et (cost 2sen t) 0
0 0 et
(1 t) t t
c) eAt = e2t (2t t2 ) (1 t + t2 ) (t + t2 )
2 2 2
2 2
3t + 2t (2t t2 ) (1 + 2t t2 )
2 2
1 2t t2 t 2t + t2
d) eAt = et t 1 t
t2 t2
1 (2t + 2 ) t 1 + 2t + 2
2
1 6t
+ 2e 31 + 13 e6t 16 + 16 e6t
e) eAt = 12 + 12 e6t 23 + 13 e6t 16 + 16 e6t
12 + 12 e6t 31 + 13 e6t 56 + 16 e6t
(12) x(t) = 12 cos( 1t ) + 3 sen( 1t ) + 2t1 sen( 1t )
4
zwn t
(14) a) x(t) = wA2 [1 e1z2 sen(wn 1 z 2 t + )]
n
con cos = z
c) z = 12
(15) a) x(t) = 32 e2t + e3t + 21
b) L = 12
3
c) tinf lexion
= n 2
P1 (t) 1 0.3(T t) + 0.045(T t)2
(17) P2 (t) = e0.7(T t) 1 0.3(T t)
P3 (t) 1
(19) a) 1 = 1, 2 =1, 3 = 2
23/15
b) K = 37/15
1
(21) x(t) = c1 cos n t + c2 sen n t + 2t1
(22) f (x2 + y 2 ) = n(x2 + y 2 )
(28) c) x(t) = t12 ((1 2 )cost 2sent + 1)
(29) c) x(t) = c1 t cost + c2 t sen t + c3 et
Cap. V
N
1 1 1
(2) x(t) = tet + ( 2
et + tet + ent )
n=1
(n + 1) (n + 1) (n + 1)2
5
(7) c) F () = 1e15 (e4 ( 45 2 ) + e 2 + 412 )
(10) x(t) = K6 ((t L)3 H(t L) + 3t2 LH(t) t3 H)t))
(11) c) Si u(t) = H(t), y(t) = ( 21 + 2et 32 e2t )H(t)
RESULTADOS DE EJERCICIOS SELECCIONADOS 197
Cap. IX
(1) d) Si > 0
(+ )n +( )n (+ )n ( )n
[ ]n
2 2
=
(+ )n ( )n (+ )n +( )n
2 2
Si < 0
[ ]n cos(nArctg )
sen(nArctg )
n/2
= ( )
2
sen(nArctg
) cos(nArctg
)
REFERENCIAS AL NIVEL INTRODUCTORIO
Bibliografa
[1] Boyce, E.B. and R.C. DiPrima, Elementary Dierential Equations and Bound-
ary Value Problems, New York: John Wiley, 1965.
[2] Brauer, F. and J.A. Nohel, Ordinary Dierential Equations, A First Course,
New York: W.A. Benjamin, 1967.
[3] Kreider, D.L., R.G. Kuller, and D.R. Ostberg, Elementary Dierential Equa-
tions, Reading, Mass.: Addison-Wesley, 1968.
[4] Marcellan, F., Cassasus, L., Zarzo, A. Ecuaciones Diferenciales, Madrid:
McGraw-Hill, 1990.
[5] Ince, E.L., Ordinary Dierential Equations, Dover, Nueva York, 1956.
[6] Simmons, F., Ecuaciones diferenciales con aplicaciones y Notas Historicas,
McGraw-Hill: Mexico, 1977.
[7] Wilson, H.K., Ordinary Dierential Equations, Addison-Wesley Publishing
Company, 1971.
NIVEL AVANZADO
[8] Birkho, G. and G.C. Rota, Ordinary Dierential Equations, Boston: Ginn,
1962.
[9] Hartman, P., Ordinary Dierential Equations, New York: John Wiley, 1964.
[10] Hurewic, W., Lectures on Ordinary Dierential Equations, Cambridge, Mass.:
M.I.T. Press, 1958.
[11] Lefschetz, S., Dierential Equations: Geometric Theory, New York: Inter-
science, 1957.
[12] Nemytskii, V.V. and V.V. Stepanov, Qualitative Theory of Dierential Equa-
tions, Princeton, N.J.: Princeton Univ. Press, 1960.
[13] Piccinini, L.C., Stampacchia, G., Vidossich, G., Ordinary Dierential Equa-
tions in Rn : Problems and Methods, New York: Springer-Verlag, 1984.
[14] Coddington, E.A. y Levinson, N., Theory of Ordinary Dierential Equations,
New York: McGraw-Hill, 1955.
SOBRE ESTABILIDAD
[15] Cesari, L., Asymptotic Behavior and Stability Problems in Ordinary Dieren-
tial Equations, New York: Academic Press, 1963.
[16] Coppel, W.A., Stability and Asymptotic Behavior of Dierential Equations,
Boston: D.C. Heath, 1963.
[17] Hahn, W., Theory and Application of Lyapunovs Direct Methods, Englewood
Clis, N.J.: Prentice-Hall, 1963.
[18] LaSalle, J. and S. Lefschetz, Stability by Liapunovs Direct Method, with Ap-
plications, New York: Academic Press, 1961.
SOBRE ALGEBRA LINEAL
[19] Franklin, J.N., Matrix Theory, Englewood Clis, N.J.: Prentice-Hall, 1968.
199
200 BIBLIOGRAFIA
[20] Halmos, P.R., Finite Dimensional Vector Spaces, New York: Van Nostrand,
1958.
Indice alfabetico
nilpotente, 184
no singular, 177
nodo, 159
Normas, 181