Guía Completa de Ecuaciones Diferenciales
Guía Completa de Ecuaciones Diferenciales
2023
Índice general
3 Transformada de Laplace 83
3.1 Transformada de Laplace . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
3.1.1 Transformada de funciones elementales . . . . . . . . . . . . . . . . . . . . . . 85
3.1.2 Propiedades de la transformada de Laplace . . . . . . . . . . . . . . . . . . . . 86
3.2 Transformada inversa y convolución . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
3.2.1 Transformada inversa de Laplace . . . . . . . . . . . . . . . . . . . . . . . . . . 91
3.2.2 Convolución . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 96
3.3 Solución de ecuaciones diferenciales . . . . . . . . . . . . . . . . . . . . . . . . . . . . 101
3.3.1 Ecuaciones integrales e integrodiferenciales . . . . . . . . . . . . . . . . . . . . 104
3.4 Transformada de funciones continuas a pedazos . . . . . . . . . . . . . . . . . . . . . 106
3.4.1 Transformada de una función periódica . . . . . . . . . . . . . . . . . . . . . . 113
3
3.5 Transformada de funciones impulso . . . . . . . . . . . . . . . . . . . . . . . . . . . . 115
1.1. Introducción
Una ecuación diferencial es una ecuación en la que interviene una variable dependiente (función
incógnita) y sus derivadas con respecto a una o más variables independientes. En general, las ecua-
ciones diferenciales se clasifican según su tipo, orden o linealidad.
Si la función incógnita depende de una sola variable independiente, decimos que la ecuación es
una ecuación diferencial ordinaria, por otro lado, si ésta depende de dos o más variables, la ecua-
ción se denomina ecuación diferencial a derivadas parciales. Por ejemplo, la ecuación
d2 x dx
2
+ 4 + 5x = t + 4
dt dt
es una ecuación diferencial ordinaria, mientras que
∂u ∂2 u
=4 2
∂t ∂x
es una ecuación diferencial a derivadas parciales.
El orden de una ecuación diferencial es el orden de la derivada más alta que aparece en la ecuación
diferencial. Por ejemplo
dy
x − y cos x = 0
dx
es una ecuación diferencial ordinaria de primer orden, mientras que
d2 y dy
2
− 4 + 5y = e x cos x
dx dx
es una ecuación diferencial ordinaria de segundo orden.
Una ecuación diferencial ordinaria de orden n, en general puede escribirse en la forma
y(n) = f x, y, y0 , . . . , y(n−1) (1.1)
1
2 CAPÍTULO 1. ECUACIONES DIFERENCIALES DE PRIMER ORDEN
Una ecuación diferencial ordinaria de orden n es lineal, si ésta, puede escribirse en la forma
a n ( x ) y ( n ) + · · · + a1 ( x ) y 0 + a0 ( x ) y = g ( x )
xy0 − y = e x
es una ecuación diferencial no lineal, ello debido a la presencia del producto yy0 .
Una solución de la ecuación diferencial 1.1 sobre el intervalo I es una función y = y( x ) tal que
existen y0 , y00 , . . . , y(n) y se satisface
y(n) ( x ) = f x, y( x ), y0 ( x ), . . . , y(n−1) ( x ) , para todo x ∈ I
1
y0 − y = 1 − x
x
en el intervalo (0, ∞). En efecto, reemplazando la función en la ecuación diferencial se tiene
1 0 1
y10 − y1 = x ln x − x2 − x ln x − x2
x x
= ln x + 1 − 2x − ln x + x
= 1−x
y0 ( x ) = f ( x ) g(y) (1.2)
1 0
y (x) = f (x)
g(y)
1 0
Z Z
y ( x )dx = f ( x )dx
g(y)
1.2. ECUACIONES SEPARABLES 3
1 + y2 = C1 x
ln |y| = x2 + ln x + C
donde C1 es una constante no nula. Por otro lado, la solución nula y( x ) = 0 también satisface la
ecuación diferencial. Por tanto, el conjunto de soluciones de la ecuación diferencial es
2
y( x ) = C1 xe x , y( x ) = 0
Definición 1.1 Un problema a valor inicial consiste de una ecuación diferencial junto con una
condición inicial
y0 ( x ) = f ( x, y), x∈I
y ( x0 ) = y0
y0 − 2xy = x, y (0) = 1
y0 − 2xy = x
de donde
y0
=x
1 + 2y
1.2. ECUACIONES SEPARABLES 5
Ejercicios propuestos
En los problemas 1 a 8, resuelva la ecuación diferencial
1. y0 + 2xy2 = 0
2. y0 = 1 + x + y2 + xy2
4. y0 = x2 − x2 y2 , x > 0
p
5. y0 = x3 (1 − y), y(0) = 3
6. y0 = 2 y + 1 cos x, y(π ) = 0
p
7. x2 + yy0 = 0, y (0) = 2
10. Si se intenta determinar la forma de un cable flexible y no extensible suspendido entre dos
puntos A y B a la misma altura, se pueden analizar las fuerzas que actúan sobre el cable y
obtener la ecuación diferencial
" 2 #1/2
d2 y dy
2
= k 1+
dx dx
dy
donde k es una constante positiva. Utilice la sustitución z( x ) = para reducir la ecuación
dx
diferencial de segundo orden a una ecuación diferencial de primer orden. Resuelva y exprese
la solución en términos de funciones exponenciales.
y0 + P( x )y = Q( x ) (1.3)
donde P( x ) y Q( x ) son funciones reales. Para la resolución de una ecuación diferencial lineal, con-
sideramos el método del factor integrante, el cual, consiste en multiplicar la ecuación diferencial
1.3 por el factor
R
P( x )dx
µ( x ) = e
1.3. ECUACIONES DIFERENCIALES LINEALES 7
y de esta forma, escribir el lado izquierdo de la ecuación 1.3 como una derivada. En efecto, multi-
plicando la ecuación por el factor integrante, obtenemos
R R
P( x )dx
y0 + P( x )y = e P( x )dx
e Q( x )
R 0 R
e P(x)dx y = e P( x )dx
Q( x )
2
y0 − y = x2 cos x, x>0
x
Solución Se trata de una ecuación diferencial lineal, y además, el factor integrante de la ecuación
es
(− 2x )dx −2
R
µ( x ) = e = e−2 ln x = eln x = x −2
x −2 y = sen x + C
De donde
y0 + 2xy = x3
Para calcular la integral indefinida del lado derecho, hacemos el cambio de variable
u = x2 =⇒ du = 2xdx
Luego,
1
Z Z
2
x3 e x dx = ueu du
2
Integrando por partes la integral del lado derecho, se tiene
Z
ueu du = ueu − eu + C = (u − 1)eu + C
Luego,
2 1 1 2
ex y = ((u − 1)eu + C ) = ( x2 − 1)e x + C1
2 2
Por tanto, la solución general de la ecuación diferencial es
1 2 2
y( x ) = ( x − 1) + C1 e−x
2
Ejemplo 1.7 Resuelva la ecuación diferencial
dy
( x 2 + 1) = 1 − 4xy
dx
Solución Reescribiendo la ecuación, tenemos
dy 4x 1
+ 2 y= 2 .
dx x + 1 x +1
La ecuación es lineal, y su factor integrante es
2
µ( x ) = e (4x/(1+ x ))dx = e2 ln(1+ x ) = (1 + x2 )2
R 2
1
(1 + x 2 ) y = x + x 3 + C
3
Despejando y, obtenemos la solución general de la ecuación lineal
C + x + 13 x3
y( x ) =
1 + x2
Ejemplo 1.8 Resuelva la siguiente ecuación diferencial
dx 1 − 2x
− = 4t + et
dt t
Solución Reescribiendo la ecuación, tenemos
dx 2 1
+ x = + 4t + et
dt t t
que es una ecuación lineal cuyo factor integrante es µ(t) = t2 . Multiplicando la ecuación por este
factor, obtenemos
2 dx 2 2 1 t
t + x =t + 4t + e
dt t t
(t2 x )0 = t + 4t3 + t2 et
Integrando con respecto a t, se tiene
1 2
Z
t2 x = t + t4 + t2 et dt
2
Como
Z Z
t2 et dt = t2 et − 2tet dt = t2 et − (2tet − 2et ) = (t2 − 2t + 2)et
Tenemos
1 2
t2 x = t + t4 + (t2 − 2t + 2)et + C
2
Despejando x (t), obtenemos la solución general de la ecuación diferencial
1 2 2 C
x (t) = + t + 1 − + 2 et + 2
2
2 t t t
Ejercicios propuestos
En los problemas 1 a 5, resuelva la la ecuación diferencial lineal
1. ( xy0 − 1) ln x = 2y
2. y0 sen x + y cos x = x sen x
10 CAPÍTULO 1. ECUACIONES DIFERENCIALES DE PRIMER ORDEN
3. ty0 = −3y + t3 − t2
4. (1 + x2 )y0 − 2xy = (1 + x2 )2
1 − 2x
0
5. y + y=1
x2
6. Resuelva el problema a valor inicial
1
y0 + y = 2 cos t, y (0) = −1
2
7. Resuelva el problema a valor inicial
y0 − y = 1 + 3 sen t, y (0) = y0
donde y0 es una constante.
8. Considere el siguiente método para resolver ecuaciones diferenciales de primer orden
y0 + P( x )y = Q( x ) (1.5)
a) Si Q( x ) = 0, muestre que la solución es
R
y( x ) = Ce− P( x )dx
y + xy0 = 0
es una ecuación diferencial exacta, pues la función φ( x, y) = xy satisface:
∂φ ∂φ
= y, =x
∂x ∂y
M ( x, y) + N ( x, y)y0 = 0 (1.8)
es exacta si, y solamente si
∂M ∂N
( x, y) = ( x, y), ∀( x, y) ∈ R
∂y ∂x
En otras palabras, existe una función diferenciable φ = φ( x, y) tal que
φx ( x, y) = M ( x, y) y φy ( x, y) = N ( x, y).
Solución. Sean
Como
My = −2x = Nx
por el teorema anterior, la ecuación diferencial es exacta, luego, existe una función φ( x, y) tal que:
φ( x, y) = x3 − x2 y + 2x + g(y)
donde g(y) es una función diferenciable a determinarse. Reemplazando esta expresión en la se-
gunda ecuación, se tiene
∂ 3
( x − x2 y + 2x + g(y)) = 6y2 − x2 + 3
∂y
− x2 + g0 (y) = 6y2 − x2 + 3
de donde
g0 (y) = 6y2 + 3
g(y) = 2y3 + 3y
Por tanto:
φ( x, y) = x3 − x2 y + 2x + 2y3 + 3y
x3 − x2 y + 2x + 2y3 + 3y = C
dy
3x2 y + 8xy2 + ( x3 + 8x2 y + 12y2 ) =0
dx
Solución Denotemos por
Como
∂M ∂N
= 3x2 + 16xy y = 3x2 + 16xy,
∂y ∂x
∂φ ∂φ
= M( x, y) = 3x2 y + 8xy2 y = N ( x, y) = x3 + 8x2 y + 12y2
∂x ∂y
φ( x, y) = x3 y + 4x2 y2 + 4y3 + g( x )
De donde
g0 ( x ) = 0 =⇒ g( x ) = C1
Y, por tanto,
φ( x, y) = x3 y + 4x2 y2 + 4y3 + C1
x3 y + 4x2 y2 + 4y3 = C
Solución Sean
M( x, y) = 2y2 + ye xy y N ( x, y) = 4xy + xe xy + 2y
My = 4y + e xy + xye xy y Nx = 4y + e xy + xye xy
son iguales, la ecuación diferencial es exacta, luego, existe una función φ = φ( x, y) tal que
φx = 2y2 + ye xy y φy = 4xy + xe xy + 2y
φ( x, y) = 2xy2 + e xy + g(y)
14 CAPÍTULO 1. ECUACIONES DIFERENCIALES DE PRIMER ORDEN
2xy2 + e xy + y2 = C
M( x, y) + N ( x, y)y0 = 0
es una ecuación diferencial no exacta. Consideramos una función µ = µ( x, y) tal que
µM( x, y) + µN ( x, y)y0 = 0
sea exacta, es decir:
(µM)y = (µN ) x
Derivando parcialmente cada una de las expresiones, obtenemos
µy M + µMy = µ x N + µNx
de donde:
µy M − µ x N = ( Nx − My )µ (1.9)
Si conseguimos encontrar una función µ que satisfaga la ecuación diferencial parcial 1.9, entonces
la ecuación diferencial
µM( x, y) + µN ( x, y)y0 = 0
es exacta. La ecuación diferencial parcial 1.9 no siempre es posible resolverla, pero en algunos
casos, el problema de determinar el factor integrante es sencillo. Por ejemplo, si
Nx − My
M
es una expresión que solo depende de y, entonces considerando µ = µ(y), la ecuación diferencial
parcial se transforma en
dµ Nx − My
= µ
dy M
1.4. ECUACIONES DIFERENCIALES EXACTAS Y FACTORES INTEGRANTES 15
Nx − My
N
es una expresión que solo depende de x, entonces considerando µ = µ( x ), la ecuación diferencial
parcial se simplifica a
dµ Nx − My
=− µ
dx N
la cual nuevamente, es una ecuación diferencial lineal.
Solución Sean
M ( x, y) = y y N ( x, y) = 2xy − e−2y
Como
My = 1 6= 2y = Nx
la ecuación diferencial no es exacta. Por tanto, consideramos una función µ = µ( x, y) tal que
Derivando, obtenemos
dµ
Si asumimos µ = µ(y), entonces µ x = 0 y µy = . Luego,
dy
dµ
y + µ = 2µy
dy
de donde
dµ 2y − 1
= µ
dy y
1 dµ 1
= 2−
µ dy y
16 CAPÍTULO 1. ECUACIONES DIFERENCIALES DE PRIMER ORDEN
y2 dy
− 2ye x + (y − e x ) =0
2 dx
Solución Sean
y2
M ( x, y) = − 2ye x y N ( x, y) = y − e x
2
Entonces
My = y − 2e x y Nx = −e x
Por tanto, la ecuación diferencial es inexacta. Suponga que existe µ = µ( x, y) tal que
2
y dy
µ( x, y) − 2ye + µ( x, y)(y − e x )
x
=0 (1.11)
2 dx
dµ
µ(y − 2e x ) = (y − e x ) − µe x
dx
dµ
µ(y − e x ) = (y − e x )
dx
de donde
dµ
=µ =⇒ µ( x ) = e x
dx
Reemplazando esta expresión en la ecuación 1.11, tenemos que
1 2 x dy
y e − 2ye2x + (ye x − e2x ) =0
2 dx
es una ecuación diferencial exacta. Luego, existe una función φ( x, y) tal que
1 2 x
φx = y e − 2ye2x y φy = ye x − e2x
2
Integrando la segunda ecuación con respecto a y, tenemos
1 2 x
φ( x, y) = y e − ye2x + g( x )
2
18 CAPÍTULO 1. ECUACIONES DIFERENCIALES DE PRIMER ORDEN
x2 y0 + xy = −y−3/2
( xy + y−3/2 ) + x2 y0 = 0
x a yb ( xy + y−3/2 ) + x a yb x2 y0 = 0
( x a+1 yb+1 + x a yb−3/2 ) + x a+2 yb y0 = 0
Para que la ecuación resultante sea exacta, debemos tener
Como esta ecuación diferencial es exacta, existe una función φ = φ( x, y) tal que
2 5/2 5/2
φ( x, y) = x y + g( x )
5
para alguna función g( x ). Reemplazando en la primera ecuación
2 5/2 5/2
x y + g( x ) = x3/2 y5/2 + x1/2
5 x
x3/2 y5/2 + g0 ( x ) = x3/2 y5/2 + x1/2
g0 ( x ) = x1/2
Integrando, tenemos
2 3/2
g( x ) = x
3
Luego,
3( xy)5/2 + 5x3/2 = C1
Ejercicios propuestos
1. Verifique que las siguientes ecuaciones diferenciales son exactas y luego determine sus solu-
ciones.
y
a) + 6x + (ln x − 2)y0 = 0
x
b) (9x2 + y − 1) − (4y − x )y0 = 0
x y dy
c) 2 2 3/2
+ 2 2 3/2
=0
(x + y ) ( x + y ) dx
d) (2x − y) + (2y − x )y0 = 0; y (1) = 3
20 CAPÍTULO 1. ECUACIONES DIFERENCIALES DE PRIMER ORDEN
3. Encuentre el valor de b para el cual, la siguiente ecuación diferencial es exacta, y luego re-
suelva usando el valor encontrado de b
( xy2 + bx2 y) + ( x + y) x2 y0 = 0
4. Encuentre todas las funciones f ( x ) para los cuales, la siguiente la ecuación diferencial
y2 sen x + y f ( x )y0 = 0
sea exacta. Resuelva la ecuación diferencial para las funciones encontradas.
5. Considere la ecuación
y2 + 2xy − x2 y0 = 0.
a) Demuestre que esta ecuación no es exacta
b) Demuestre que multiplicando ambos lados de la ecuación por y−2 se obtiene una ecua-
ción exacta.
c) Use la solución de la ecuación resultante para resolver la ecuación original.
d) ¿Se perdieron soluciones en el proceso?
6. Encuentre un factor integrante y resuelva la ecuación diferencial dada
a) (3x2 y + 2xy + y3 ) + ( x2 + y2 )y0 = 0
b) 1 + ( x/y − sen y)y0 = 0
7. La ecuación diferencial f ( x )y0 + x2 + y = 0 tiene como factor integrante a la función
µ( x ) = x. Encuentre todas las funciones posibles f ( x ).
8. La ecuación diferencial e x sec y − tan y + y0 = 0 tiene un factor integrante de la forma
µ( x, y) = e−ax cos y para alguna constante a. Encuentre el valor de a, luego resuelva la ecua-
ción diferencial.
9. Verifique que la ecuación diferencial
( x2 + 2xy − y2 ) + (y2 + 2xy − x2 )y0 = 0
no es exacta. Multiplique la ecuación diferencial por µ( x, y) = ( x + y)−2 y verifique que la
nueva ecuación diferencial es exacta y resuélvalo.
10. (a) Demuestre que si ( Nx − My )/( xM − yN ) depende solo del producto xy, esto es,
Nx − My
= H ( xy),
xM − yN
entonces la ecuación M ( x, y) + N ( x, y)y0 = 0 tiene un factor integrante de la forma
µ( xy). Dar la fórmula general para µ( xy).
a) Use el resultado del inciso (a) para encontrar una solución implı́cita a
(3y + 2xy2 ) + ( x + 2x2 y)y0 = 0,
satisfaciendo la condición inicial y(1) = 1.
1.5. SUSTITUCIONES DIVERSAS 21
dv 3 1 dv v2 − 1
x + v = v − v −1 =⇒ x = .
dx 2 2 dx 2v
Separando variables, tenemos
2v dv 1
=
v2 − 1 dx x
22 CAPÍTULO 1. ECUACIONES DIFERENCIALES DE PRIMER ORDEN
ln |v2 − 1| = ln x + C
v2 = 1 + C1 x
y2 = x2 (1 + C1 x )
Finalmente, la solución general de la ecuación diferencial esta conformada por las dos familias de
funciones
p p
y( x ) = x 1 + C1 x, y( x ) = − x 1 + C1 x
xy0 = x cos(y/x ) + y
dy y y
= cos +
dx x x
Haciendo el cambio de variable u = y/x, tenemos
dy du
y = ux =⇒ = x +u
dx dx
Reemplazando lo obtenido en la ecuación diferencial, tenemos
du
x + u = cos(u) + u
dx
De donde obtenemos
du
x = cos u.
dx
Separando variables e integrando cada miembro de la ecuacion, tenemos
du 1
sec(u) = =⇒ ln | sec u + tan u| = ln | x | + C
dx x
y
Como u = , eliminando logarı́tmos de la ecuación anterior, concluimos que las soluciones de la
x
ecuación diferencial satisfacen
y y
sec + tan = Cx
x x
1.5. SUSTITUCIONES DIVERSAS 23
dy
y−α + P ( x ) y 1− α = Q ( x ) (1.12)
dx
Haciendo el cambio de variable v = y1−α , obtenemos
dv dy dy 1 dv
= (1 − α ) y − α =⇒ y−α =
dx dx dx 1 − α dx
Reemplazando en la ecuación 1.12, finalmente obtenemos
1 dv
+ P( x )v = Q( x )
1 − α dx
la cual es una ecuación diferencial lineal.
dy 4
− y = −6ty2 , t>0
dt t
Solución Claramente la ecuación es de Bernoulli. Dividiendo la ecuación por y2 , tenemos
dy 4 −1
y −2 − y = −6t
dt t
Haciendo el cambio de variable v = y−1 , se tiene
dv dy dy dv
= − y −2 =⇒ y −2 =−
dt dt dt dt
Reemplazando
dv 4
− − v = −6t
dt t
o bien,
dv 4
+ v = 6t
dt t
El factor integrante de la ecuación es
4
R
e t dt = e4 ln t = t4
24 CAPÍTULO 1. ECUACIONES DIFERENCIALES DE PRIMER ORDEN
y 0 + P ( x ) y + Q ( x ) y2 = R ( x )
v = y − y1 =⇒ y = y1 + v
1 4
y 0 + y − y2 = − 2 , x>0
x x
2
Sabiendo que y1 ( x ) = es una solución particular de la misma.
x
Solución Se trata de una ecuación de Ricatti. Por tanto, haciendo el cambio de variable
2 2
v = y− =⇒ y= +v
x x
y reemplazando en la ecuación diferencial, obtenemos
2
2 0 1 2 2 4
− 2 +v + +v − +v = − 2
x x x x x
2 2 1 4 4 4
− 2 + v 0 + 2 + v − 2 − v − v2 = − 2
x x x x x x
0 3 2
v − v=v
x
2 Si en lugar de hacer el cambio de variable y = y1 + v hacemos
1
y = y1 + ,
v
la ecuación resultante es una ecuación diferencial lineal de primer orden.
26 CAPÍTULO 1. ECUACIONES DIFERENCIALES DE PRIMER ORDEN
w = v −1 =⇒ w 0 = − v −2 v 0 =⇒ v 0 = − v2 w 0
Luego,
3 3
− v2 w 0 − v = v2 =⇒ w 0 + w = −1
x x
El factor integrante de la ecuación lineal es
3
R
e x dx = e3 ln x = x3
Multiplicando la ecuación por este factor, tenemos
0 3
x w + w = − x3
3
x
( x3 w)0 = − x3
Integrando, obtenemos
1 C1 − x4
x3 w = − x4 + C =⇒ w=
4 4x3
Reemplazando w = v−1 , tenemos
C1 − x4 4x3
v −1 = =⇒ v=
4x3 C1 − x4
2
Finalmente, como v = y − , tenemos:
x
2 4x3 2 4x3
y− = =⇒ y( x ) = +
x C1 − x4 x C1 − x4
2
Observe que y1 ( x ) = no pertenece a familia de soluciones obtenida, por tanto, el conjunto de
x
soluciones de la ecuación diferencial es:
2 4x3 2
y( x ) = + , y( x ) = .
x C1 − x4 x
1 2 1 4
y0 = y − y−
2x x x
Solución Observamos que y1 ( x ) = −2 es una solución particular de la ecuación diferencial. Hace-
mos el cambio de variable
1
y( x ) = −2
z( x )
1.5. SUSTITUCIONES DIVERSAS 27
Ejercicios propuestos
En los problemas 1 a 9, resuelva las siguientes ecuaciones diferenciales mediante un cambio de
variables adecuado
dy
1. y + ( x − 2y) =0
dx
28 CAPÍTULO 1. ECUACIONES DIFERENCIALES DE PRIMER ORDEN
dy
2. xy2 = y3 − x 3 , y (1) = 2
dx
p dy
3. (t − ty) =y
dt
dy
4. = y( xy3 − 1)
dx
dy
5. xy = y2 + x
dx
dy 1
6. x2 − 2xy = 3y4 , y (1) =
dx 2
dy p
7. = 2 + y − 2x + 3
dx
dy
8. = tan2 ( x + y)
dx
dy t+y+1
9. =
dt t−y+3
10. Resuelva las siguientes ecuaciones de Ricatti conociendo una de sus soluciones particulares
dy
a) = y2 − y sen x + cos x, y1 ( x ) = sen x
dx
dy
b) + y2 − (1 + 2e x )y + e2x = 0, y1 ( x ) = e x
dx
y0 = f ( x, y), y ( x0 ) = y0
es importante saber de antemano si esta tiene solución y si esa solución es la única. Por ejemplo, el
problema valor inicial
√
y0 = 2 y, y (0) = 0
∂f
Teorema 1.3 (Existencia y unicidad de soluciones) Sean f ( x, y) y ( x, y) funciones continuas so-
∂y
bre un rectángulo R centrado en ( x0 , y0 ). Entonces, existe h > 0 tal que el problema a valor inicial
y0 = f ( x, y), y ( x0 ) = y0
tiene una única solución definida en el intervalo ( x0 − h, x0 + h).
Ejemplo 1.22 Muestre que el problema a valor inicial tiene solución única definida en un intervalo
alrededor de su valor inicial dado
x
y 0 = x 2 − y2 + 8 , y (3) = −1
y
x
Solución Sea f ( x, y) = x2 − y2 + 8 , entonces
y
∂f x
= −2y − 8 2
∂y y
Como ambas funciones son continuas en un rectángulo centrado en (3, −1), por el teorema ante-
rior, concluimos que el problema a valor inicial tiene solución única definida en un intervalo de
(3 − h, 3 + h) para algun h > 0.
Ejemplo 1.23 Muestre que el problema a valor inicial, tiene infinitas soluciones
xy0 = y, y (0) = 0
Determine una región del plano xy donde el problema a valor inicial tiene solución única.
Solución La ecuación diferencial puede escribirse en la forma
1 0 1
y =
y x
Integrando, tenemos
ln |y| = ln | x | + C
de donde y( x ) = Cx. Claramente, esta familia de funciones, son soluciones del problema a valor
inicial. Esto se debe a lo siguiente. Reescribiendo la ecuación diferencial en la forma estándar
x
y0 =
y
tenemos que
x ∂f x
f ( x, y) = y ( x, y) = − 2
y ∂y y
no son continuas en ninguna región que contenga el punto (0, 0). Por otro lado, las funciones
f ( x, y) y f y ( x, y) son continuas, ya sea en el semiplano y > 0 o bien en el semiplano y < 0, y por
tanto, si la condición inicial ( x0 , y0 ) satisface y0 6= 0, el problema a valor inicial
xy0 = y. y ( x0 ) = y0
tiene solución única definida en algún intervalo centrado en x0 .
30 CAPÍTULO 1. ECUACIONES DIFERENCIALES DE PRIMER ORDEN
f ( x, y, C ) = 0 (1.13)
fx
f x + f y y0 = 0 =⇒ y0 ( x ) = − ( x, y)
fy
Por tanto, si y = y( x ) es una trayectoria ortogonal a la familia de curvas 1.13, entonces su derivada
debe satisfacer la ecuación diferencial
1 fy
y0 ( x ) = − fx
( x, y) = ( x, y)
− fx
fy
x2 + y2 + Cx = 0
C + 2x
2x + 2yy0 + C = 0 =⇒ y0 = − , y 6= 0
2y
Por tanto, las trayectorias ortogonales a la familia de curvas, son soluciones de la ecuación
diferencial
2y
y0 =
C + 2x
2y 2xy
y0 = = 2
x2 +y 2 x − y2
− + 2x
x
La ecuación diferencial obtenida es una ecuación diferencial homogénea, por tanto, haciendo el
y
cambio de variable v = , obtenemos
x
dv 2v dv v + v3
x +v = =⇒ x =
dx 1 − v2 dx 1 − v2
1.7. APLICACIONES DE LAS ECUACIONES DIFERENCIALES DE PRIMER ORDEN 31
1 − v2 dx
Z Z
3
dv =
v + v x
Z
1 2v dx
Z
− dv =
v 1 + v2 x
ln |v| − ln(1 + v2 ) = ln | x | + C
Despejando v, obtenemos
v
= C1 x
1 + v2
y
Finalmente, reemplazando v = , obtenemos
x
y
x = C2 x
y2
1+ 2
x
que es equivalente a
x2 + y2 = C3 y
y
4
−2 2 4x
−2
−4
Figura 1.1: Trayectorias ortogonales
32 CAPÍTULO 1. ECUACIONES DIFERENCIALES DE PRIMER ORDEN
Ejemplo 1.25 Un objeto de masa 5 kg se libera desde el reposo a 1000 m sobre el suelo y se le per-
mite caer bajo la influencia de la gravedad. Suponiendo que la fuerza debida a la resistencia del
aire es proporcional a la velocidad del objeto, con constante de proporcionalidad b = 50 N-s/m,
determine la ecuación de movimiento del objeto ¿Cuándo tocará el objeto al suelo?
Solución Por la segunda ley de Newton, la velocidad v = v(t) del objeto satisface (considere la
dirección del movimiento como positiva)
dv
m = mg − bv, v (0) = 0
dt
donde m = 5, g = 9.81 y b = 50. Reescribiendo la ecuación diferencial
dv b
+ v=g
dt m
vemos que es una ecuación diferencial lineal, cuyo factor integrante es
R
(b/m)dt
µ(t) = e = e(b/m)t
Multiplicando la ecuación diferencial por este factor,obtenemos
(b/m)t dv b
e + v = ge(b/m)t
dt m
d
e(b/m)t v = ge(b/m)t
dt
1.7. APLICACIONES DE LAS ECUACIONES DIFERENCIALES DE PRIMER ORDEN 33
mg m2 g
x (t) = t + 2 e−(b/m)t + C
b b
Considerando como punto de referencia el punto donde se libera el objeto, se tiene x (0) = 0.
Reemplazando esta condición inicial en la ecuación anterior
mg m2 g m2 g
x (0) = · 0 + 2 e−(b/m)·0 + C = 0 =⇒ C=−
b b b2
Luego
mg m2 g m2 g
x (t) = t + 2 e−(b/m)t − 2
b b b
Reemplazando los datos, se tiene
5 · 9, 81 52 · 9, 81 −(50/5)t
−10t
x (t) = t+ e − 1 = 0, 981t + 0, 0981 e − 1
50 502
El objeto tocará el suelo cuando x (t) = 1000, por tanto, debemos encontrar el valor t tal que
0, 981t + 0, 0981 e−10t − 1 = 1000
B Q 5 (1, t)
P 5 (x, y)
x
0 (1, 0)
dy
Haciendo w = , la ecuación puede escribirse como
dx
dw βp
( x − 1) =− 1 + w2
dx α
O bien
1 dw β 1
√ =
1 + w2 dx α1−x
Integrando, se tiene
β
senh−1 (w) = − ln(1 − x ) + C
α
dy
Por otro lado, cuando x = 0, se tiene w(0) = dx (0) = 0. Luego
β
senh−1 (0) = − ln 1 + C =⇒ C = 0
α
De donde, despejando w, obtenemos
β
w = senh − ln(1 − x )
α
De donde, la ecuación de la curva de persecución y = y( x ) satisface la ecuación diferencial
dy 1
= (1 − x )− β/α − (1 − x ) β/α (1.16)
dx 2
Ejemplo 1.26 En el caso del problema anterior, determine la posición donde el barco A alcanza al
barco B, si α > β.
Solución Integrando la ecuación 1.16, tenemos
1 1 1+ β/α 1 1− β/α
y( x ) = (1 − x ) − (1 − x ) +C
2 1 + β/α 1 − β/α
Como y = 0 cuando x = 0, se tiene
1 1 1 αβ
0= − +C =⇒ C=
2 1 + β/α 1 − β/α α2 − β2
Luego, la posición del barco es dado por
1 1 1+ β/α 1 1− β/α αβ
y( x ) = (1 − x ) − (1 − x ) + 2
2 1 + β/α 1 − β/α α − β2
Finalmente, el barco A alcanza al barco B, cuando x = 1, de donde
αβ
y (1) =
α2 − β2
6 L/min
0.1 kg/L x(t)
1000 L
x(0) 5 0 kg 6 L/min
El objetivo es encontrar una expresión para la cantidad de sal x (t) que hay en el tanque en el ins-
tante t.
dx
En efecto, la razón de cambio de la cantidad de sal en el tanque, en el instante t, dt , debe ser igual
a la razón a la que la sal entra menos la razón a la que sale, es decir
dx
= re − rs
dt
Donde la razón de entrada es dada por el producto de la concentración de sal por la tasa de flujo
re = (0.1 kg/L) · (6L/min) = 0.6 kg/min
Y de forma análoga, la razón de salida es
x (t) 3
rs = kg/L · (6 L/min) = x (t) kg/min
1000 500
Por tanto, la cantidad de sal presente en el tanque en el instante t es dada por
dx 3
= 0.6 − x ( t ), x (0) = 0
dt 500
la cual, es una ecuación lineal de primer orden
dx 3
+ x = 0.6
dt 500
cuyo factor integrante es e3t/500 . Multiplicando la ecuación por este factor y luego integrando,
tenemos
3t/500 dx 3
e + x = 0.6e3t/500
dt 500
d 3t/500
e x = 0.6e3t/500 =⇒ e3t/500 x = 100e3t/500 + C
dt
Despejando x (t), se tiene
Ejemplo 1.27 Un tanque con una capacidad de 500 L contiene originalmente 200 L con 50 kg de
sal en solución. Se hace entrar 12 kg de sal por litro, a razón de 3 L por minuto, y se deja que la
mezcla salga del tanque a razón de 2 L/min. Encuentre la cantidad de sal presente en el tanque en
cualquier instante antes del momento en que la solución comienza a derramarse.
Solución Sea x (t) la cantidad de sal presente en el tanque en el instante t. Entonces
dx
= re − rs
dt
donde
1 3
re = ·3 =
2 2
x (t) 2
rs = ·2 = x (t)
200 + (3 − 2)t 200 + t
Luego, x (t), satisface la ecuación diferencial
dx 2 3
+ x= , x (0) = 50
dt 200 + t 2
La ecuación diferencial es lineal, con factor integrante
R
2/(200+t)dt
e = e2 ln(200+t) = (200 + t)2
Multiplicando la ecuación diferencial por este factor, y luego integrando, obtenemos
2 dx 2 3
(200 + t) + x = (200 + t)2
dt 200 + t 2
d 3 1
(200 + t)2 x (t) = (200 + t)2 =⇒ (200 + t)2 x (t) = (200 + t)3 + C
dt 2 2
de donde
1 C
x (t) = (200 + t) +
2 (200 + t)2
Por otro lado, cuando t = 0, x = 50. Luego
C
50 = 100 + =⇒ C = −50 · 2002
2002
Finalmente, la cantidad de sal presente en el tanque antes de que comience a derramarse es
2
1 200
x (t) = (200 + t) − 50
2 200 + t
38 CAPÍTULO 1. ECUACIONES DIFERENCIALES DE PRIMER ORDEN
Problemas propuestos
1. Encuentre las trayectorias ortogonales a la familia de curvas dada
a) y = cx2
b) y2 − x2 = c
c) y = ecx
2. Hallar una familia de trayectorias oblicuas que corte a la familia de circunferencias x2 + y2 =
c2 formando un ángulo de 45o .
3. Un objeto de 5 kg recibe una velocidad inicial hacia abajo de 50 m/s y luego se le permite caer
bajo la influencia de la gravedad. Suponga que la fuerza en newtons debida a la resistencia
del aire es de −10v, donde v es la velocidad del objeto en m/s. Determine la ecuación de
movimiento del objeto. Si el objeto está inicialmente a 500 m sobre el suelo, determine el
momento en que el objeto golpeará el suelo.
4. Un objeto de masa 8 kg recibe una velocidad inicial hacia arriba de 20 m/s y luego se le
permite caer bajo la influencia de la gravedad. Suponga que la fuerza en newtons debida a la
resistencia del aire es −16v, donde v es la velocidad del objeto en m/s. Determine la ecuación
de movimiento del objeto. Si el objeto está en un principio a 100 m sobre el suelo, determine
el momento en que el objeto golpeará el suelo.
5. Un cohete con masa inicial m0 kg se lanza de manera vertical desde el suelo. El cohete arroja
gas a la razón constante de α kg/s y a una velocidad constante de β m/s con respecto del
cohete. Suponga que la magnitud de la fuerza gravitacional es proporcional a la masa, con
constante de proporcionalidad g. Como la masa no es constante, la segunda ley de Newton
conduce a la ecuación
dv
(m0 − αt) − αβ = − g(m0 − αt)
dt
donde v = dx/dt es la velocidad del cohete, x es la altura sobre el suelo y m0 − αt es la
masa del cohete t segundos después del lanzamiento. Si la velocidad inicial es cero, resuelva
la ecuación anterior para determinar la velocidad del cohete y su altura sobre el suelo para
0 ≤ t < m0 /α.
6. Una paracaidista cuya masa es de 75 kg se arroja de un helicóptero que vuela a 2000 m sobre
el suelo y este cae hacia abajo, bajo la influencia de la gravedad. Suponga que la fuerza debida
a la resistencia del aire es proporcional a la velocidad de la paracaidista, con la constante de
proporcionalidad b1 = 30 N-s/m cuando el paracaı́das está cerrado y b2 = 90 N-s/m cuando
se abre. Si el paracaı́das no se abre hasta que la velocidad de la paracaidista es de 20 m/s, ¿
después de cuántos segundos llegará ella al suelo?
7. Un conejo parte del origen y corre por el eje y positivo con velocidad a. Al mismo tiempo, un
perro que corre con velocidad b sale del punto (c, 0) y persigue al conejo. ¿ Qué trayectoria
sigue el perro?
a) Suponga que a < b. ¿ Qúe distancia recorre el conejo antes de que el perro lo alcance?
b) Suponga que a = b. ¿ A qué distancia del conejo consigue acercarse el perro?
8. Un gran tanque está parcialmente lleno con 400 L de fluido en los que se disolvieron 5 kg de
sal. La salmuera tiene 0.25 kg de sal por litro que entra al tanque a razón de 20 L/min. La
solución bien mezclada sale del tanque a razón de 15 L/min. Determine la cantidad de kg de
sal que hay en el tanque después de 30 minutos.
1.7. APLICACIONES DE LAS ECUACIONES DIFERENCIALES DE PRIMER ORDEN 39
9. Un tanque grande contiene inicialmente 800 litros de salmuera en la que se han disuelto 8
kilogramos de sal. A partir del instante t = 0 comienzan a entrar 3.5 L/min de salmuera
conteniendo 250 gramos de sal por litro. La mezcla se conserva uniforme y abandona el
tanque a razón de 8 L/min.
a) Hallar la cantidad de sal en el tanque al cabo de una hora
b) Hallar la cantidad de sal en el tanque cuando éste contiene solamente 200 litros de sal-
muera.
10. Una solución salina entra a una razón constante de 6 litros/minuto en un tanque de gran
tamaño que en un principio contenı́a 50 litros de solución salina en que se habı́an disuelto
0.5 kg de sal. La solución dentro del tanque se mantiene bien revuelta y sale del tanque
con la misma razón. Si la concentración de sal en la solución que entra al tanque es de 0.05
kg/litro, determine la masa de sal en el tanque después de t minutos. ¿ Cuándo llegará la
concentración de sal en el tanque a 0.03 kg/litro?
40 CAPÍTULO 1. ECUACIONES DIFERENCIALES DE PRIMER ORDEN
Capı́tulo 2
En este capı́tulo vamos a estudiar ecuaciones diferenciales lineales de orden superior a coeficientes
constantes.
a n ( x ) y ( n ) + a n −1 ( x ) y ( n −1) + · · · + a 1 ( x ) y 0 + a 0 ( x ) y = b ( x )
x2 y00 + 6xy0 − 4y = 0
es una ecuación diferencial lineal homogénea de segundo orden. Mientras que la ecuación
diferencial
a n ( x ) y ( n ) + a n −1 ( x ) y ( n −1) + · · · + a 1 ( x ) y 0 + a 0 ( x ) y = b ( x )
y ( x 0 ) = y 0 , y 0 ( x 0 ) = y 1 , · · · , y ( n −1) ( x 0 ) = y n −1
41
42 CAPÍTULO 2. ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR
e x , e2x , e3x
e x e2x e3x
x 2x 3x
W (e , e , e ) = e x 2e2x 3e3x
e x 4e2x 9e3x
1 1 1
= e6x 1 2 3
1 4 9
= 2e3x
Como W (e x , e2x , e3x ) 6= 0 para todo x, concluimos que las funciones e x , e2x y e3x son linealmente
independientes.
Por otra lado, las funciones
1, e− x , e1− x
1 e− x e 1− x
−x 1− x
W (1, e ,e )= 0 − e − x − e 1− x
0 e− x e 1− x
=0
a n ( x ) y ( n ) + a n −1 ( x ) y ( n −1) + · · · + a 1 ( x ) y 0 + a 0 ( x ) y = 0
Teorema 2.5 (Existencia de un conjunto fundamental) Toda ecuación diferencial lineal homogénea
posee un conjunto fundamental de soluciones.
44 CAPÍTULO 2. ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR
a n ( x ) y ( n ) + a n −1 ( x ) y ( n −1) + · · · + a 1 ( x ) y 0 + a 0 ( x ) y = 0
y ( x ) = c1 y1 ( x ) + · · · + c n y n ( x )
xy000 − y00 = 0
en el intervalo (0, ∞). En efecto, inialmente verificamos que las funciones son soluciones de la
ecuación diferencial dada
x (1)000 − (1)00 = 0 − 0 = 0
x ( x )000 − ( x 00 ) = 0 − 0 = 0
x ( x3 )000 − ( x3 )00 = 6x − 6x = 0
1 x x3
3
W (1, x, x ) = 0 1 3x2 = 6x > 0, si x > 0
0 0 6x
y ( x ) = c1 + c2 x + c3 x 3 .
a n ( x ) y ( n ) + a n −1 ( x ) y ( n −1) + · · · + a 1 ( x ) y 0 + a 0 ( x ) y = b ( x )
y ( x ) = c1 y1 ( x ) + · · · + c n y n ( x ) + y p ( x )
y00 + y = 2e x
2.1. ECUACIONES DIFERENCIALES LINEALES 45
y00 + y = 0
y( x ) = c1 cos x + c2 sen x + e x
a n ( x ) y ( n ) + a n − 1 ( x ) y ( n − 1 ) + · · · + a 1 ( x ) y 0 + a 0 ( x ) y = bi ( x ) , i = 1, 2, . . . , k
y00 + 2y = 4 − 2x
y00 + 2y = 4
y00 + 2y = −2x
y00 + 2y = 4 − 2x
Ejercicios propuestos
1. La familia de funciones y( x ) = c1 e x + c2 e− x es la solución general de la ecuación diferencial
y00 − y = 0
sobre el intervalo (−∞, ∞). Encuentre un miembro de la familia que satisfaga las condiciones
iniciales y(0) = 0 y y0 (0) = 1.
3. Si se especifican los valores de una solución a una ecuación diferencial en dos puntos dis-
tintos, estas condiciones se llaman condiciones en la frontera. (En contraste, las condiciones
iniciales especifican valores de una función y su derivada en el mismo punto). La finalidad
de este ejercicio es mostrar que los problemas con valores en la frontera no tienen un teore-
ma de existencia y unicidad similar al de los problemas con valores iniciales. Dado que toda
solución de
y00 + y = 0 (2.1)
tiene la forma y(t) = c1 cos t + c2 sen t, donde c1 y c2 son constantes arbitrarias, muestre que
a) Existe una única solución de 2.1 que satisface las condiciones en la frontera y(0) = 2 y
y(π/2) = 0.
b) No existe una solución de 2.1 que satisfaga y(0) = 2 y y(π ) = 0.
c) Existe una infinidad de soluciones de 2.1 que satisfacen y(0) = 2 y y(π ) = −2
a) f ( x ) = cos x, g( x ) = cos 2x
b) f ( x ) = e− x , g( x ) = xe− x
c) f ( x ) = e x , g( x ) = cosh x, h( x ) = senh x
10. Suponga que y1 = e x y y2 = e− x son dos soluciones de una ecuación diferencial homogénea.
Explique por qué y3 = cosh x y y4 = senh x son también soluciones de la ecuación.
y00 + P( x )y0 + Q( x )y = 0
Si y1 ( x ) es una solución particular dada, entonces una segunda solución y2 ( x ) linealmente inde-
pendiente, puede ser determinada planteando
y2 ( x ) = u ( x ) y1 ( x )
donde u( x ) satisface una ecuación diferencial de primer orden, de esta forma podemos determinar
un conjunto fundamenta de soluciones de la ecuación de segundo orden. Este método se denomi-
na, método de reducción de orden. En efecto, calculando las primeras dos derivadas de y2 , se tiene:
y20 = u0 y1 + uy10
y200 = u00 y1 + 2u0 y10 + uy100
y00 − 4y0 + 4y = 0
y2 ( x ) = u( x )e2x
u00 = 0 =⇒ u 0 = c1 =⇒ u ( x ) = c1 x + c2
y2 ( x ) = (c1 x + c2 )e2x
y2 ( x ) = xe2x
y( x ) = c3 e2x + c4 xe2x
1
y00 + y0 = 0,
x
sabiendo que y1 = ln x es una solución particular.
y2 = u( x ) ln x
2.2. REDUCCIÓN DE ORDEN 49
u
y20 ( x ) = u0 ln x +
x
00 00 u0 u0 u
y2 ( x ) = u ln x + + − 2
x x x
2u0 u
= u00 ln x + − 2
x x
Reemplazando en la ecuación diferencial, se tiene
u0 u 1 0 u
u00 ln x + 2 − 2+ u ln x + =0
x x x x
2 ln x
(ln x )u00 + + u0 = 0
x x
o bien,
00 2 1
u + + u0 = 0
x ln x x
Esta, es una ecuación diferencial lineal en la variable u0 , cuyo factor integrante es:
Z
2 1
+ dx
e x ln x x = e2 ln(ln x ) + ln x = x ln2 x
c1 c1
u0 = 2
=⇒ u=− + c2
x ln x ln x
Considerando c1 = −1 y c2 = 0, tenemos
1
y2 ( x ) = ln x = 1
ln x
y1 ( x ) = ln x, y2 ( x ) = 1
forman un conjunto fundamental de soluciones de la ecuación diferencial, por tanto, una solución
general es dada por
y( x ) = c3 ln x + c4 .
50 CAPÍTULO 2. ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR
Ejercicios propuestos
En los problemas 1 a 8, use el método de reducción de orden para determinar una segunda solución
linealmente independiente y luego halle una solución general para la misma
2. y00 − y = 0, y1 = cosh x
3. x2 y00 + 3xy0 + y = 0, y 1 = x −1
4. x2 y00 + 2xy0 − 6y = 0, y1 = x 2
5. (1 − x2 )y00 + 2xy0 = 0, y1 = 1
6. x2 y00 + ( x2 − x )y0 + (1 − x )y = 0, y1 = x
7. xy00 − ( x + 3)y0 + 3y = 0, y1 = e x
9. Una ecuación diferencial de segundo orden y00 = f ( x, y0 ) que no contiene la variable depen-
diente y, puede ser resuelta reduciendo a la ecuación diferencial de primer orden
du
= f ( x, u)
dx
a) x2 y00 + (y0 )2 = 0
b) e− x y00 = (y0 )2
c) y0 y00 − x = 0 y(1) = 2, y0 (1) = 1
10. Una ecuación diferencial de segundo orden y00 = f (y, y0 ) que no contiene la variable inde-
pendiente x, puede ser resuelta reduciendo a la ecuación diferencial de primer orden
du
u = f (y, u)
dy
a) yy00 + (y0 )2 + 1 = 0
b) y00 + 2y(y0 )3 = 0
c) y0 y00 = 2, y(0) = 1, y0 (0) = 2
52 CAPÍTULO 2. ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR
a n y ( n ) + a n −1 y ( n −1) + · · · + a 1 y 0 + a 0 y = 0 (2.2)
donde an , an−1 , . . . , a0 son constantes con an 6= 0. Planteamos como solución de la ecuación dife-
rencial 2.2 una función de la forma
y( x ) = eλx
donde λ es una constante a ser determinada. Reemplazando esta función en la ecuación diferencial,
obtenemos
(n) ( n −1) 0
an eλx + an−1 eλx + · · · + a1 eλx + a0 eλx = 0
( an λn + an−1 λn−1 + · · · + a1 λ + a0 )eλx = 0
de donde concluimos que λ debe satisfacer la ecuación polinomial denominada ecuación carac-
terı́stica asociada a la ecuación diferencial 2.2
a n λ n + a n −1 λ n −1 + · · · + a 1 λ + a 0 = 0 (2.3)
y ( x ) = e λ0 x
forman un conjunto fundamental de soluciones de la ecuación diferencial 2.2 y por tanto, una
solución general de la misma es
y ( x ) = c 1 e λ1 x + c 2 e λ2 x + · · · + c n e λ n x
y00 − y0 − 2y = 0
y( x ) = c1 e2x + c2 e− x .
λ3 − 4λ2 + 3λ = 0
λ(λ2 − 4λ + 3) = 0
λ(λ − 1)(λ − 3) = 0
cuyas raı́ces son
λ1 = 0, λ2 = 1, λ3 = 3
Por tanto, la solución general de la ecuación diferencial homogénea es
y( x ) = c1 e0x + c2 e x + c3 e3x
= c1 + c2 e x + c3 e3x
λ2 + 2λ − 8 = 0
(λ + 4)(λ − 2) = 0 =⇒ λ1 = −4 y λ2 = 2
Luego una solución general es dada por
y( x ) = c1 e−4x + c2 e2x
Para hallar la solución que satisface las condiciones iniciales, reemplazamos las condiciones inicia-
les en la solución general y su derivada
y00 + 6y0 + 9y = 0
y( x ) = c1 e−3x + c2 xe−3x
y000 − 3y0 − 2y = 0
y( x ) = c1 e2x + c2 e− x + c3 xe− x
e(a+ib)x y e(a−ib)x
y1 ( x ) = e ax cos(bx ) y y2 ( x ) = e ax sen(bx )
Ejercicios propuestos
En los ejercicios 1 a 5, halle una solución general para las ecuaciones diferenciales
1. y00 − 5y0 + 6y = 0
2. 4y00 − 4y0 + y = 0
3. y000 + 3y00 − 4y0 − 12y = 0
4. y00 − 2y0 + 26y = 0
5. y000 + 2y00 + 5y0 − 26y = 0
En los ejercicios 6 a 8, halle la solución al problema a valor inicial
6. y00 − 4y0 + 3y = 0; y(0) = 1, y0 (0) = 1/3
7. y00 + 2y0 + 17y = 0; y(0) = 1, y 0 (0) = −1
8. y000 − 4y00 + 7y0 − 6y = 0; y(0) = 1, y0 (0) = 0, y00 (0) = 0
9. Una clase importante de ecuaciones diferenciales con coeficientes variables son aquellas que
pueden expresarse en la forma
d2 y dy
ax2 + bx + cy = h( x ) (2.4)
dx2 dx
donde a, b y c son constantes y y( x ) es la función incógnita. Tales ecuaciones reciben el nom-
bre de ecuaciones de Cauchy-Euler o equidimensionales. El objetivo de este problema es
mostrar que la sustitución x = et transforma la ecuación 2.4 en una ecuación con coeficientes
constantes en la nueva variable independiente t.
a) Use la regla de la cadena para mostrar que para x = et ,
dy dy
=x
dt dx
b) Derive la ecuación en la parte (a) con respecto a t para obtener
d2 y dy 2
2d y
= + x
dt2 dt dx2
c) Sustituya las fórmulas de las partes (a) y (b) en la ecuación 2.4 para obtener la ecuación
a coeficientes constantes
d2 y dy
a 2
+ (b − a) + cy = h(et )
dt dt
1. Use la técnica del problema anterior para obtener una solución general de las siguientes
ecuaciones de Cauchy-Euler para x > 0.
a) x2 y00 − 3xy0 + 6y = 0
b) x2 y00 + 9xy0 + 17y = 0
c) x2 y00 + 3xy0 + 5y = 0
2.4. ECUACIONES DIFERENCIALES LINEALES NO HOMOGÉNEAS 57
a n y ( n ) + a n −1 y ( n −1) + · · · + a 1 y 0 + a 0 y = b ( x )
donde b( x ) es una función de tipo polinomial, exponencial, combinación lineal de senos y cosenos
o una combinación de los anteriores mediantes sumas y productos.
y p ( x ) = A k x k + A k −1 x k −1 + · · · + A 1 x + A 0
donde Ak , Ak−1 , . . . , A0 son constantes para los cuales y p ( x ) es una solución particular de la ecua-
ción diferencial no homogénea dada.
y00 − 4y0 + 5y = 5x − 4
Solución Planteamos como solución particular a la función
y p ( x ) = Ax + B
cuya primera y segunda derivada son respectivamente
y0p ( x ) = A
y00p ( x ) = 0
y00 − y = x2 + 1
Solución Primero resolvamos la ecuación diferencial
y00 − y = 0
cuya ecuación carecterı́stica es
λ2 − 1 = 0 =⇒ (λ + 1)(λ − 1) = 0
y cuyas raı́ces son: λ1 = −1 y λ2 = 1. Por tanto, la solución complementarı́a es
y c ( x ) = c1 e − x + c2 e x
Por otro lado,como la parte no homogénea de la ecuación diferencial es un polinomio de segundo
grado, planteamos como solución particular a
y p ( x ) = Ax2 + Bx + C
y0p ( x ) = 2Ax + B
y00p ( x ) = 2A
y00p − y p = x2 + 1
2A − ( Ax2 + Bx + C ) = x2 + 1
− Ax2 − Bx + (2A − C ) = x2 + 1
Igualando coeficientes, se tiene
−A = 1
−B = 0
2A − C = 1
cuyas soluciones son: A = −1, B = 0 y C = −3. Luego, la función
y p ( x ) = − x2 − 3
y p ( x ) = x m ( A k x k + A k −1 x k −1 + · · · + A 1 x + A 0 )
y00 − 2y0 = 4x − 6
Solución Primero hallemos una solución general para la ecuación diferencial homogénea asociada
y00 − 2y0 = 0
cuya ecuación caracterı́stica es
λ2 − 2λ = 0
λ ( λ − 2) = 0
yc ( x ) = c1 + c2 e2x
y p ( x ) = x ( Ax + B) = Ax2 + Bx
y0p ( x ) = 2Ax + B
y00p ( x ) = 2A
y00p − 2y0p = 4x − 6
2A − 2(2Ax + B) = 4x − 6
−4Ax + (2A − 2B) = 4x − 6
Igualando coeficientes, obtenemos el sistema lineal
−4A = 4
2A − 2B = −6
y p ( x ) = − x2 + 2x
y( x ) = yc ( x ) + y p ( x )
= c1 + c2 e2x − x2 + 2x
60 CAPÍTULO 2. ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR
y p ( x ) = Aeλx
donde A es una constante para el cual, y p ( x ) es una solución particular de la ecuación diferencial
no homogénea.
y00 − 4y = 3e− x
y p ( x ) = Ae− x
y0p ( x ) = − Ae− x
y00p ( x ) = Ae− x
y00p − 4y p = 3e− x
Ae− x − 4 Ae− x = 3e− x
−3Ae−x = 3e−x
y p ( x ) = −e− x
λ2 − 2λ − 3 = 0 =⇒ (λ − 3)(λ + 1) = 0
yc ( x ) = c1 e3x + c2 e− x
Por otro lado, planteamos como solución particular de la ecuación no homogénea a la función
y p ( x ) = Ae−2x
2.4. ECUACIONES DIFERENCIALES LINEALES NO HOMOGÉNEAS 61
y0p ( x ) = −2Ae−2x
y00p ( x ) = 4Ae−2x
Por tanto, y p ( x ) = e−2x es una solución particular de la ecuación no homogénea. Luego, la solución
general es dada por
y( x ) = c1 e3x + c2 e− x + e−2x
c1 + c2 = −2
3c1 − c2 = 2
y( x ) = −2e− x + e−2x
y p = Ax m eλx
y00 − 2y0 + y = 8e x
y00 − 2y0 + y = 0
λ2 − 2λ + 1 = 0
( λ − 1)2 = 0 =⇒ λ1 = 1 (se repite dos veces)
62 CAPÍTULO 2. ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR
yc (t) = c1 e x + c2 te x
y p ( x ) = x2 ( Ae x ) = Ax2 e x
y00p − 2y0p + y p = 8e x
A( x2 + 4x + 2)e x − 2A( x2 + 2x )e x + Ax2 e x = 8e x
2Ae x = 8e x =⇒ A=4
y p ( x ) = 4x2 e x
y( x ) = yc ( x ) + y p ( x )
= c1 e x + c2 xe x + 4x2 e x
y p = A cos(kx ) + B sen(kx )
y00 + y0 − 2y = 10 sen x
y00 + y0 − 2y = 0
λ2 + λ − 2 = 0
(λ + 2)(λ − 1) = 0 =⇒ λ1 = −2 y λ2 = 1
2.4. ECUACIONES DIFERENCIALES LINEALES NO HOMOGÉNEAS 63
yc ( x ) = c1 e−2x + c2 e x
y p ( x ) = A cos x + B sen x
−3A + B = 0
− A − 3B = 10
y p ( x ) = − cos x − 3 sen x
Observación 2.11 Si ki es una raı́z compleja de la ecuación caracterı́stica con multiplicidad m, en-
tonces debe plantear como solución particular a
y p = x m ( A cos(kx ) + B sen(kx ))
y00 + 4y = 12 cos 2x
Solución Primero hallemos, una solución complementaria para la ecuación diferencial. Para ello,
considere la ecuación caracterı́stica asociada a la ecuación
λ2 + 4 = 0 =⇒ (λ + 2i )(λ − 2i ) = 0
cuyas soluciones son λ1 = −2i y λ2 = 2i. Luego, la solución complementaria es dada por
yc ( x ) = c1 cos 2x + c2 sen 2x
64 CAPÍTULO 2. ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR
Para hallar una solución particular para la ecuación no homogénea, observamos que la parte no
homogénea es b( x ) = 12 cos 2x y que λ = 2i es una raı́z de multiplicidad uno, por ello, planteamos
como solución particular a
y p ( x ) = x ( A cos 2x + B sen 2x )
y00p + 4y p = 12 cos 2x
(4B − 4Ax ) cos 2x − (4A + 4Bx ) sen 2x + 4x ( A cos 2x + B sen 2x ) = 12 cos 2x
4B cos 2x − 4A sen 2x = 12 cos 2x
4B = 12
−4A = 0
y00 − y = 2x sen x
y00p − y p = 2x sen x
(2C − B − Ax ) cos x − (2A + D + Cx ) sen x − ( Ax + B) cos x − (Cx + D ) sen x = 2x sen x
(2C − 2B − 2Ax ) cos x − (2A + 2D + 2Cx ) sen x = 2x sen x
(2C − 2B) cos x − 2Ax cos x − (2A + 2D ) sen x − 2Cx sen x = 2x sen x
Solución La parte no homogénea es producto de un polinomio de primer grado, por ello deberia-
mos considerar como solución particular una función de la forma ( Ax + B)e− x , pero, observamos
que e− x es una solución de la ecuación homogénea, por tanto, planteamos como solución particular
a
y p ( x ) = x ( Ax + B)e− x = ( Ax2 + Bx )e− x
−8A = 4
2A − 4B = 1
4A − 2B = 5
2A + 4B = 0
1
y p ( x ) = e− x cos x − e− x sen x
2
es una solución particular de la ecuación no homogénea.
2.4. ECUACIONES DIFERENCIALES LINEALES NO HOMOGÉNEAS 67
Ejercicios propuestos
En los ejercicios 1 a 6, resuelva la ecuación diferencial por el método de coeficientes indeterminados
1. y00 + y0 − 6y = 6x − 2
2. y00 − 2y0 + y = 8e x
3. y00 + 4y = −2 sen 2x
4. y000 − 2y00 = 2 + e− x
5. y00 + y = 2x cos x
6. y000 − 3y00 + 3y0 − y = 3x − 12e x
En los ejercicios 7 a 10, resuelva el problema a valor inicial usando el método de coeficientes inde-
terminados
7. y00 + y0 = −6x, y(0) = 0, y0 (0) = −1.
8. y00 − y = −4e− x , y(0) = −1, y0 (0) = 1
9. y00 + 2y0 + y = 2x sen x, y(0) = 1, y0 (0) = 1
10. 2y00 − y0 = xe x − 1, y(0) = −1, y0 (0) = 0
68 CAPÍTULO 2. ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR
y00 + 4y = 2 tan x
y00 + 4y = 0
λ2 + 4 = 0 =⇒ λ1 = 2i, λ2 = −2i
yc ( x ) = c1 cos 2x + c2 sen 2x
y p ( x ) = c1 ( x ) cos 2x + c2 ( x ) sen 2x
c10
cos 2x sen 2x 0
=
−2 sen 2x 2 cos 2x c20 2 tan x
0 sen 2x
2 tan x 2 cos 2x −2 tan x sen 2x
c10 = = = − tan x sen 2x
cos 2x sen 2x 2
−2 sen 2x 2 cos 2x
cos 2x 0
−2 sen 2x 2 tan x 2 tan x cos 2x
c20 = = = tan x cos 2x
cos 2x sen 2x 2
−2 sen 2x 2 cos 2x
70 CAPÍTULO 2. ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR
y00 + 2y0 − 3y = 0
es
r2 + 2r − 3 = 0
(r + 3)(r − 1) = 0
y cuyas raı́ces son λ1 = −3 y λ2 = 1. Por tanto, la solución complementaria de la ecuación 3w
yc ( x ) = c1 e−3x + c2 e x
2.4. ECUACIONES DIFERENCIALES LINEALES NO HOMOGÉNEAS 71
y p = c1 ( x )e−3x + c2 ( x )e x
luego, la función
1 4x 1 4x −3x
y p ( x ) = − xe + e e + x2 e x
2 8
1 1 x
= x2 − x + e
2 8
es una solución particular. Finalmente, la solución general de la ecuación diferencial no homogénea
es
−3x x 2 1 1 x
y ( x ) = c1 e + c2 e + x − x + e
2 8
1
= c1 e−3x + c3 e x + x2 − x e x
2
x2 y00 − 2xy0 + 2y = x3 ln x
72 CAPÍTULO 2. ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR
y c ( x ) = c1 x + c2 x 2 .
Para hallar una solución particular, planteamos
y p ( x ) = c1 ( x ) x + c2 ( x ) x 2
de donde
0 x2
x ln x 2x − x3 ln x
c10 = = = − x ln x
x x2 x2
1 2x
x 0
1 x ln x x2 ln x
c20 = = = ln x
x x2 x2
1 2x
Integrando, obtenemos
x2 1
Z
1 2 1 1
Z
c1 = (− x ln x )dx = − x ln x − − dx = − x2 ln x + x2
2 2 x 2 4
1
Z Z
c2 = ln xdx = x ln x − x dx = x ln x − x
x
Por tanto, una solución particular de la ecuación diferencial es dada por
1 2 1 2
y p ( x ) = − x ln x + x x + ( x ln x − x ) x2
2 4
1 3
= x3 ln x − x3
2 4
Finalmente, la solución general de la ecuación diferencial es
1 3
y( x ) = c1 x + c2 x2 + x3 ln x − x3
2 4
Ejemplo 2.27 Encuentre la solución general de la ecuación diferencial
π π
y000 − y00 + y0 − y = sec x, − <x<
2 2
2.4. ECUACIONES DIFERENCIALES LINEALES NO HOMOGÉNEAS 73
y000 − y00 + y0 − y = 0
Considerando la ecuación la ecuación caracterı́stica y resolviendo esta, se tiene
λ3 − λ2 + λ − 1 = 0
λ2 ( λ − 1) + ( λ − 1) = 0
(λ2 + 1)(λ − 1) = 0 =⇒ λ1 = i, λ2 = −i, λ3 = 1
De donde, al solución complementaria es
yc ( x ) = c1 cos x + c2 sen x + c3 e x
Para hallar una solución particular, planteamos
y p ( x ) = c1 ( x ) cos x + c2 ( x ) sen x + c3 ( x )e x
Ejercicios propuestos
En los problemas 1 a 6, determine una solución general de la ecuación diferencial usando el método
de variación de parámetros
1. x 00 − 2x 0 + x = t−1 et
2. x 00 + x = csc t
3. y00 + y = tan2 x
4. y00 + 2y0 + y = e− x ln x
5. y00 + y = sec3 x
6. y00 + y = 3 sec x + 1 − x2
7. Use el método de variación de parámetros para mostrar que
Z t
y(t) = c1 cos t + c2 sen t + f ( x ) sen(t − s)ds
0
es una solución general de la ecuación diferencial y00 + y = f (t), donde f (t) es una función
continua en (−∞, ∞).
8. Dado que y1 ( x ) = 1 + x y y2 ( x ) = e x son soluciones de la ecuación homogénea asociada a
xy00 − (1 + x )y0 + y = x2 e2x , x > 0,
determine una solución particular
9. La ecuación de Bessel de orden un medio,
2 00 0 1 2
x y + xy + x − y=0
4
tiene dos soluciones linealmente independientes,
y1 ( x ) = x −1/2 cos x y y2 ( x ) = x −1/2 sen x.
Determine una solución general de la ecuación no homogénea
1
x2 y00 + xy0 + x2 − y = x5/2 , x>0
4
0 = mg − kL =⇒ mg = kL (2.7)
donde k > 0 es la constante del resorte. Ahora bien, si u(t) denota el desplazamiento de la masa
hacia abajo desde su posición de equilibrio en el instante t, entonces por la segunda ley de Newton
donde w es el peso de la masa dada por w = mg. La fuerza del resorte Fr (t) es una fuerza que se
opone al movimiento de la masa, y es dada por
Fr (t) = −k ( L + u(t))
donde γ > 0 se denomina constante de amortiguamiento. Y finalmente, Fe (t) es una fuerza externa
aplicada al sistema resorte-masa. Considerando, lo anterior, la ecuación 2.8 puede escribirse como:
Tomando en cuenta la ecuación 2.7, la ecuación del movimiento de la masa se puede expresar como
l
l+L+u
m
m
u
Ejemplo 2.28 (Movimiento libre no amortiguado) Una masa que pesa 2 lb alarga 6 pulg. un re-
sorte. Si se tira hacia abajo la masa otras 3 pulg. más y luego se suelta, y si no hay resistencia del
aire, determine la posición de la masa en cualquier instante t.
76 CAPÍTULO 2. ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR
mu00 + ku = 0
w 2 lb 1
w = mg =⇒ m= = = lb-s2 /pie
g 32 pie/s2 16
w 2 lb
w = kL =⇒ k= = = 4 lb/pie
L 0.5 pies
1 00
u + 4u = 0 =⇒ u00 + 64u = 0
16
1
u (0) = pies, u0 (0) = 0 pies/s
4
Considerando estas condiciones iniciales, concluimos que la ecuación del movimiento de la masa
es
1
u(t) = cos(8t)
4
u(t)
0.2
0.1
2 4 6 8 10 12 14 16 18 t
−0.1
−0.2
Ejemplo 2.29 (Movimiento libre amortiguado) Una masa que pesa 16 lb estira un resorte 3 pulg.
La masa está sujeta a un amortiguador viscoso con constante de amortiguamiento de 2 lb-s/pie.
Si la masa se pone en movimiento desde su posición de equilibrio con una velocidad hacia abajo
de 3 pulg./s, encuentre su posición en cualquier instante t. Determine cuándo la masa vuelve por
primera vez a su posición de equilibrio.
mu00 + γu0 + ku = 0
w 16 lb 1
w = mg =⇒ m= = 2
= lb-s2 /pie
g 32 pies/s 2
w 16 lb
w = kL =⇒ k= = = 64 lb/pie
L 0.25 pies
Además, la constante de amortiguamiento es dada por γ = 2 lb-s/pie. Por tanto, la ecuación del
movimiento se reduce a
1 00
u + 2u0 + 64u = 0
2
O bien
1 √
u(t) = √ e−2t sen(2 31t)
8 31
de donde
√ π
2 31t = π =⇒ t = √ ≈ 0.282 s
2 31
78 CAPÍTULO 2. ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR
u(t) ·10−2
1.5
0.5
−0.5
Ejemplo 2.30 (Movimiento forzado no amoortiguado) Una masa que pesa 8 lb alarga 6 pulg. un
resorte. Sobre el sistema actúa una fuerza externa de 8 sen 8t lb. Si se tira hacia abajo de la masa
3 pulg. y luego se suelta, determine la posición de la masa en cualquier instante. Determine las
primeras cuatro veces en las que la velocidad de la masa es cero.
Solución La ecuación del movimiento de la masa es
mu00 + ku = Fe (t)
donde la masa m y la constante del resorte son dados por
w 8 lb 1
w = mg =⇒ m= = = lb-s2 /pie
g 32 pies/s2 4
w 8 lb
w = kL =⇒ k= = = 16 lb/pie
L 1/2 pies
Y además, la fuerza externa es Fe = 8 sen 8t lb. Luego
1 00
u + 16u = 8 sen 8t
4
o bien
u00 + 64u = 32 sen 8t
cuya solución general es
u(t) = c1 cos 8t + c2 sen 8t − 2t cos 8t
Por otro lado, las condiciones iniciales del problema son
1
u (0) = pies, u0 (0) = 0 pies/s
4
2.5. APLICACIONES DE LAS ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR 79
u(t)
40
20
5 10 15 20 25t
−20
−40
Ejemplo 2.31 (Movimiento forzado amortiguado) Un resorte se alarga 6 pulg. por una masa que
pesa 8 lb. La masa está sujeta a un mecanismo de amortiguamiento de 0.25 lb-s/pie y sobre ella
actúa una fuerza externa de 4 cos 2t lb. Determine la respuesta de estado estable de este sistema.
Solución La ecuación diferencial que describe el movimiento de la masa es
mu00 + γu0 + ku = Fe (t)
donde la masa y la constante del resorte son dados por
w 8 lb 1
w = mg =⇒ m= = 2
= lb-s2 /pie
g 32 pies/s 4
w 8 lb
w = kL =⇒ k= = = 16 lb/pie
L 0.5 pies
80 CAPÍTULO 2. ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR
1 00
u + 0.25u0 + 16u = 4 cos 2t =⇒ u00 + u0 + 64u = 16 cos 2t
4
u ( t ) = c 1 e λ1 t + c 2 e λ2 t + u p ( t )
donde λ1 y λ2 son las raı́ces complejas de la ecuación caracterı́stica asociada a la ecuación ho-
mogénea cuya parte real es negativa, debido a la constante de amortiguamiento, por tanto, con-
cluimos que la solución complementaria, denominada respuesta de estado transitorio satisface
y por tanto
Por ello, la solución particular u p (t) se denomina respuesta de estado estable del sistema. Para
hallar la solución particular, planteamos
60A + 2B = 16
−2A + 60B = 0
240 8
u p (t) = cos 2t + sen 2t
901 901
8
= (30 cos 2t + sen 2t)
901
2.5. APLICACIONES DE LAS ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR 81
Ejercicios propuestos
1. Una masa que pesa 3 lb. estira 3 pulg. un resorte. Si la masa se empuja hacia arriba, contra-
yendo el resorte una distancia de 1 pulg. y luego se pone en movimiento con una velocidad
hacia abajo de 2 pies/s, y no hay resistencia del aire, encuentre la posición de la masa en cual-
quier instante t. Determine la frecuencia, el periodo, la amplitud y la fase del movimiento.
2. Una masa que pesa 10 lb estira un resorte 3 pulg. La masa está sujeta a un amortiguador
viscoso con constante de amortiguamiento de 2 lb-s/pie. Si la masa se pone en movimiento
desde su posición de equilibrio con una velocidad hacia abajo de 3 pulg./s, encuentre su
posición en cualquier instante t. Determine cuándo la masa vuelve por primera vez a su
posición de equilibrio.
Transformada de Laplace
1 −(s−2)t A
= lı́m − e
A→∞ s−2 0
1 −(s−2) A 1
= lı́m − e +
A→∞ s−2 s−2
1
= , si s > 2
s−2
Ejemplo 3.2 Determine la transformada de Laplace de la función
1 − t, si 0 ≤ t < 1
n
f (t) = 0, si t ≥ 1
83
84 CAPÍTULO 3. TRANSFORMADA DE LAPLACE
Solución Observe que ésta, es una función continua a pedazos en el intervalo [0, ∞). Aplicando la
definición de la transformada de Laplace, tenemos
Z ∞
L[ f ](s) = f (t)e−st dt
0
Z 1 Z ∞
−st
e−st dt
= (1 − t)e dt + 0 ·
0 1
− 1
st e−st
Z 1
e
= (1 − t ) − (−1) ds, Integrando por partes
−s 0 0 −s
1
1 e−st
= + 2
s s 0
−
1 e s 1
= + 2 − 2, s 6= 0
s s s
Por otro lado, para s = 0 tenemos
Z ∞
L[ f ](0) = f (t)e−0t dt
0
Z 1
= (1 − t)dt
0
1
1
= − (1 − t )2
2 0
1
=
2
Por tanto, la transformada de Laplace de la función f (t) es la función continua (verifique)
−s
1− 1 +e ,
s 6= 0
L[ f ](s) = s s2 s2
1
, s=0
2
2
Ejemplo 3.3 Muestre que la función f (t) = et , t ≥ 0 no tiene transformada de Laplace.
Solución En efecto, usando la definición
Z ∞
2
L[ f ](s) = et e−st dt
Z0 ∞
2 − st
= et dt
0
Z ∞ 2
s2 s
= e− 4 e(t− 2 ) dt
0
Z ∞
2 2
Como e(t−s/2) dt = ∞ para todo s ∈ R, concluimos que la función f (t) = et no tiene transfor-
0
mada de Laplace.
A continuación, se establecen condiciones suficientes para la existencia de la transformada de La-
place de una función.
3.1. TRANSFORMADA DE LAPLACE 85
En efecto,
Z ∞ Z ∞
f (t)e−st dt ≤ | f (t)|e−st dt
M M
Z ∞
≤ Ke at e−st dt
M
Z ∞
=K e−(s−a)t dt
M
Como la integral impropia del lado derecho converge para todo s > a, concluimos que la transfor-
mada de Laplace L[ f ](s) existe para todo s > a.
k
L[k](s) = , s>0
s
Función exponencial
Si k es una constante, entonces
1
L[ekt ](s) = , s>k
s−k
Función potencia
Si n es un entero positivo, entonces
n!
L[tn ](s) = , s>0
s n +1
86 CAPÍTULO 3. TRANSFORMADA DE LAPLACE
Como
Z A
k k cos( Ak) + s · sen( Ak) −sA
sen(kt)e−st dt = − e
0 s2 +k 2 s2 + k 2
y además
k cos( Ak) + s · sen( Ak) −sA |k| + |s| −sA
2 2
e ≤ 2 e ≤ e−sA
s +k s + k2
Concluimos que
Z A
k
lı́m sen(kt)e−st dt = , si s > 0
A→∞ 0 s2 + k2
Por tanto,
k
L[sen(kt)](s) = , s>0
s2 + k2
Propiedad de linealidad
Sean f 1 y f 2 funciones cuyas transformadas de Laplace existen para s > s1 y s > s2 respectivamente
y c es una constante. Entonces,
L[ f 1 + f 2 ](s) = L[ f 1 ](s) + L[ f 2 ](s), s > máx{s1 , s2 }
L[c · f ](s) = c · L[ f ](s), s > s1
Ejemplo 3.5 Determine la transformada de Laplace de la función
f (t) = 2 + 5 cos(2t)
Solución Tenemos
L[ f ](s) = L[2 + 5 cos(2t)](s)
= L[2](s) + 5L[cos(2t)](s)
2 s
= +5 2
s s +4
7s2 + 8
= , s>0
s ( s2 + 4)
Ejemplo 3.6 Si k > 0, demuestre que
s
L[cosh(kt)](s) = , s>k
s2 − k2
Solución En efecto, como
ekt + e−kt
cosh(kt) =
2
Usando la propiedad de linealidad, tenemos
1
L[cosh(kt)](s) = L[ekt + e−kt ](s)
2
1 1 1
= +
2 s−k s+k
s
= 2 , s>k
s − k2
Por tanto,
d
L[t f (t)](s) = − L[ f ](s)
ds
Ejemplo 3.9 Determine la transformada de Laplace de la función
f (t) = t sen(3t)
Solución Tenemos
L[ f ](s) = L[t sen(3t)](s)
d
= − L[sen(3t)](s)
ds
d 3
=−
ds s2 + 9
6s
= 2 , s>0
( s + 9)2
Ejemplo 3.10 Determine la transformada de Laplace de la función usando la propiedad de la de-
rivada
f ( t ) = t2 e − t
Solución. Tenemos
L[ f ](s) = L[t2 e−t ](s)
d2
= 2 L[e−t ](s)
ds
d2
1
= 2
ds s + 1
d 1
= −
ds ( s + 1)2
2
= , s > −1
( s + 1)3
Ejemplo 3.11 Determine la transformada de Laplace de la función
Ejercicios propuestos
Ejer. 1-5 Usando la definición de la transformada de Laplace, encuentre la transformada de la
función
1. f (t) = (t + 1)e−at , t ≥ 0
2. f (t) = senh( at), t ≥ 0
t, 0≤t<1
n
3. f (t) = 1, t≥1
0, 0≤t<2
n
4. f (t) = 1, t>2
0, 0≤t<2
(
5. f (t) = 1, 2<t<4
0, t>4
a) f (t) = e−2t t3
b) f (t) = tet + et sen 3t + e2t cos 5t
c) f (t) = ebt g(t), donde
0, 0≤t<a
n
g(t) = 1, t>a
3.2. TRANSFORMADA INVERSA Y CONVOLUCIÓN 91
a) f (t) = t cos 2t
b) f (t) = t2 senh t
c) f (t) = te−3t cos 3t
1
L[ f ( at)](s) = F (s/a), s > as0
a
donde a > 0. Utilice el resultado anterior, para encontrar la transformada de Laplace de las
funciones
a) f (t) = ekt
b) f (t) = sen kt
c) f (t) = 1 − cos kt
sen kt
a) f (t) =
t
1 − cos kt
b) f (t) =
t
En particular, tenemos:
−1 k
L (t) =k
s
−1 1
L (t) = ekt
s−k
n!
L −1 n +1 ( t ) = tn
s
−1 s
L (t) = cos(kt)
s2 + k 2
k
L −1 2 (t) = sen(kt)
s + k2
−1 k
L (t) = senh(kt)
s2 − k 2
s
L −1 2 (t) = cosh(kt)
s − k2
1
F (s) =
s−4
Solución Como
1
L[e4t ](s) =
s−4
concluimos que
−1 1
L (t) = e4t
s−4
3
F (s) =
s2 − 4s + 13
Solución Como
3 3
=
s2 − 4s + 13 ( s − 2)2 + 9
3.2. TRANSFORMADA INVERSA Y CONVOLUCIÓN 93
Tenemos
−1 −1 3
L [ F (s)](t) = L (t)
(s − 2)2 + 32
3
= e2t L−1 2 (t)
s + 32
= e2t sen(3t)
Ejemplo 3.15 Determine la transformada inversa de Laplace de la función
2 12 3
F (s) = − 3 +
s s (s − 1)2 + 16
Solución En efecto
2 12 3
L−1 [ F ](t) = L−1 − 3 + (t)
s s (s − 1)2 + 16
−1 2 −1 2! 3 −1 4
=L ( t ) − 6L (t) + L (t)
s s3 4 ( s − 1 ) 2 + 42
3 4
= 2 − 6t2 + et L−1 2 (t)
4 s + 42
3
= 2 − 6t2 + et sen(4t)
4
Ejemplo 3.16 Determine la transformada inversa de la función
2s
F (s) =
(s + 2)(s − 1)
Solución Para hallar la transformada inversa, descomponemos la función racional en fracciones
parciales. En efecto
2s A B
= +
(s + 2)(s − 1) s+2 s−1
Usando el método del residuo, tenemos:
2s 2s 4
A = lı́m (s + 2) = lı́m =
s→−2 (s + 2)(s − 1) s→−2 s − 1 3
2s 2s 2
B = lı́m (s − 1) = lı́m =
s →1 (s + 2)(s − 1) s →1 s + 2 3
Luego,
4 1 2 1
F (s) = +
3s+2 3s−1
Tomando la transformada inversa, obtenemos
−1 4 −1 1 2 −1 1
L [ F (s)](t) = L (t) + L (t)
3 s+2 3 s−1
4 2
= e−2t + et
3 3
94 CAPÍTULO 3. TRANSFORMADA DE LAPLACE
De donde
3 −3 5
F (s) = + +
s − 1 s − 2 ( s − 2)2
Tomando la transformada inversa, tenemos
−1 −1 1 −1 1 −1 1
L [ F (s)](t) = 3L ( t ) − 3L ( t ) + 5L (t)
s−1 s−2 ( s − 2)2
1
= 3et − 3e2t + 5e2t L−1 2 (t)
s
= 3et − 3e2t + 5e2t t
Ejemplo 3.18 Determine la transformada inversa de la función
s−3
F (s) =
( s2 + 4)(s + 1)
Solución Descomponemos la expresión racional en fracciones parciales
s−3 A Bs + C
= + 2
(s2 + 4)(s + 1) s+1 s +4
2
A(s + 4) + ( Bs + C )(s + 1)
=
(s + 1)(s2 + 4)
( A + B)s2 + ( B + C )s + (4A + C )
=
(s + 1)(s2 + 4)
3.2. TRANSFORMADA INVERSA Y CONVOLUCIÓN 95
A+B = 0
B+C = 1
4A + C = −3
4
5A = −4 =⇒ A=−
5
Reemplazando en la primera y segunda ecuación, obtenemos
4 1
B= , C=
5 5
Luego, la transformada se puede escribir como
4 1
4 1 5s + 5
F (s) = − + 2
5s+1 s +4
4 1 4 s 1 2
=− + +
5 s + 1 5 s2 + 4 10 s2 + 4
Tomando la transformada de Laplace inversa, se obtiene
4 4 1
L−1 [ F (s)](t) = − e−t + cos 2t + sen 2t
5 5 10
Ejemplo 3.19 Determine la transformada inversa de la función
1
F (s) =
s ( s2 + 1)2
1 A Bs + C Ds + E
= + 2 + 2
s ( s2 + 1) 2 s s +1 ( s + 1)2
A(s2 + 1)2 + ( Bs + C )s(s2 + 1) + ( Ds + E)s
=
s ( s2 + 1)2
( A + B)s4 + Cs3 + (2A + B + D )s2 + (C + E)s + A
=
s ( s2 + 1)2
A+B =0
C =0
2A + B + D =0
C+E =0
A =1
96 CAPÍTULO 3. TRANSFORMADA DE LAPLACE
B = −1, D = −1 y E = 0
Luego,
1 1 s s
= − 2 − 2
s ( s2 + 1) 2 s s + 1 ( s + 1)2
Como
d d 1 2s
L[t sen t](s) = − L[sen t](s) = − =
ds ds s2 + 1 ( s2 + 1)2
tenemos que
−1 s 1
L 2 2
(t) = t sen t
( s + 1) 2
3.2.2. Convolución
Sean f (t) y g(t) funciones continuas a pedazos y definidas para t ≥ 0. La convolución de las
funciones f y g denotado por f ∗ g, es definida mediante
Z t Z t
( f ∗ g)(t) = f (t − τ ) g(τ )dτ = f (τ ) g(t − τ )dτ
0 0
1. f ∗ g = g ∗ f
2. f ∗ ( g1 + g2 ) = f ∗ g1 + f ∗ g2
3. ( f ∗ g) ∗ h = ( f ∗ g) ∗ h
4. f ∗ 0 = 0 ∗ f = 0
Solución Como
f (t) = et ∗ sen t,
98 CAPÍTULO 3. TRANSFORMADA DE LAPLACE
F (s) = L[ f (t)](s)
= L[et ∗ sen t](s)
= L[et ](s) · L[sen t](s)
1 1
= · 2
s−1 s +1
1
=
(s − 1)(s2 + 1)
t
Z
1
L f (τ )dτ (s) = F (s)
0 s
t
Z
L f (τ )dτ (s) = L[1 ∗ f (t)](s)
0
= L[1](s) · L[ f ](s)
1
= F (s)
s
El corolario 3.1, puede ser bastante útil para encontrar la transformada de Laplace inversa de al-
gunos productos
1
F (s) =
(s − 1)(s + 2)
Solución Como
1 1
F (s) = ·
s−1 s+2
3.2. TRANSFORMADA INVERSA Y CONVOLUCIÓN 99
se tiene,
−1 1 1 −1 1 −1 1
L · (t) = L (t) ∗ L (t)
s−1 s+2 s−1 s+2
= et ∗ e−2t
Z t
= eτ e−2(t−τ ) dτ
0
Z t
= e−2t+3τ dτ
0
1 −2t+3τ t
= e
3 0
1 t 1 −2t
= e − e
3 3
1
F (s) =
(s + 1)(s − 1)2
Solución Como
1 1
F (s) =
s + 1 ( s − 1)2
tenemos
−1 −1 1 −1 1
L [ F ](t) = L (t) ∗ L (t)
s+1 ( s − 1)2
= e−t ∗ tet
= tet ∗ e−t
Z t
= τeτ e−(t−τ ) dτ
0
Z t
= τe2τ −t dτ
0
t
e2τ −t
Z t
1
=τ − e2τ −t dτ
2 0 2 0
t
1 1 2τ −t
= tet − e
2 4 0
1 t 1 t 1 −t
= te − e + e
2 4 4
1
F (s) =
s ( s2 + 1)
100 CAPÍTULO 3. TRANSFORMADA DE LAPLACE
Solución En efecto
1 1
L−1 [ F ](t) = L−1 [ ](t) ∗ L−1 [ 2 ](t)
s s +1
= 1 ∗ sen t
= sen t ∗ 1
Z t
= sen τ · 1dτ
0
= (− cos τ )|0t
= − cos t + 1
Ejercicios propuestos
Ejer. 1-5 Determine la transformada inversa de Laplace de las siguientes funciones
1 1 1
1. F (s) = − +
s2 s s−2
s+1
2. F (s) =
s2− 4s
1
3. F (s) =
s2 + 4s + 5
1
4. F (s) =
s2 ( s2+ 1)
1
5. F (s) =
s(s2 + 2s + 2)
6. Demuestre las propiedades conmutativa, distributiva y asociativa de la convolución
a) f ∗ g = g ∗ f
b) f ∗ ( g1 + g2 ) = f ∗ g1 + f ∗ g2
c) f ∗ ( g ∗ h) = ( f ∗ g) ∗ h
7. Verifique
a) 1 ∗ 1 = t
1
b) 1 ∗ cos( at) = a sen( at)
t2 4 at
c) t2 ∗ sen( at) = sen2
a − a3 2
1
a) F (s) =
s3 ( s − 1)
1
b) F (s) =
s ( s + a2 )2
2
1
c) F (s) =
( s + 1) ( s2 + 4)
2
Teorema 3.3 Suponga que f es continua y f 0 es continua por pedazos para t ≥ 0 y además que f
es una función de orden exponencial cuando t → ∞. Entonces
L[ f 0 ](s) = sL[ f ](s) − f (0)
Solución En efecto
Z ∞
L[ f 0 (t)](s) = f 0 (t)e−st dt
0
Z A
= lı́m f 0 (t)e−st dt
A→∞ 0
Tomando el lı́mite cuando A → ∞ y usando el hecho que la función f (t) es de orden exponencial,
tenemos
L[ f 0 (t)](s) = sL[ f (t)](s) − f (0)
De forma análoga, si f y f 0 son funciones continuas y f 00 es continua a pedazos para t ≥ 0 y además
f y f 0 son de orden exponencial cuando t → ∞. Entonces
y0 − 4y = 1 + 2e−t ; y (0) = −1
donde Y (s) = L[y](s) y F (s) = L[ f ](s). Reemplazando los valores iniciales, tenemos
Despejando Y (s):
s+3 1
Y (s) = + F (s)
(s + 1)(s + 2) (s + 1)(s + 2)
Descomponiendo en fracciones parciales las expresiones racionales del lado derecho, tenemos
2 1 1 1
Y (s) = − + − F (s)
s+1 s+2 s+1 s+2
s2 + 1
F (s) =
(s − 2)(s − 1)2
4 2 5
F (s) = − − +
s − 1 ( s − 1)2 s−2
1 s
sF (s) − f (0) = + F (s) 2
s2 s +1
donde F (s) = L[ f ](s). Despejando F (s), tenemos
1 1
F (s) = 5
+ 3
s s
Reescribiendo la función de forma adecuada, tenemos
1 4! 1 2!
F (s) = +
4! s5 2! s3
Finalmente, tomando la transformada inversa obtenemos la solución a la ecuación integro-diferencial
1 4 1 2
f (t) = t + t
24 2
Ejemplo 3.32 Resuelva la ecuación integral
Z t
−2t
f (t) = cos(t) + e f (τ )e2τ dτ
0
Ejercicios propuestos
Ejer. 1-5 Usando la transformada de Laplace, resuelva las siguientes ecuaciones diferenciales
1. x 0 − 4x = sen(2t); x (0) = 0
2. x 00 + 8x 0 + 16x = 1 + 2t; x (0) = 0, x 0 (0) = −1
3. x 00 + 9x = e−t ; x (0) = 1, x 0 (0) = −1
4. y000 + 2y00 − y0 − 2y = sen 3t, y(0) = 0, y0 (0) = 0 y y00 (0) = 1
5. x 0 + y0 + 2x + y = 0, y0 + 5x + 3y = 1; x (0) = −1, y(0) = 4
6. Si u = u(t) es una función tal que u(0) = 0, resuelva la siguiente ecuación diferencial en
términos de la función u(t)
y00 + 4y = u(t) − u0 (t); y(0) = 0, y0 (0) = 1
f (t)
f (t)
t
a t
Las función de Heaviside se puede usar para escribir una función continua a pedazos de forma
compacta. Por ejemplo, si f (t), t ≥ 0 es una función continua por pedazos definido por
f 1 ( t ), 0≤t<a
(
f (t) = f 2 ( t ), a≤t<b
f 3 ( t ), b≤t
entonces, la función puede escribirse en la forma
f (t) = f 1 (t)[ H (t) − H (t − a)] + f 2 (t)[ H (t − a) − H (t − b)] + f 3 (t) H (t − b)
o de forma equivalente
f (t) = f 1 (t) + ( f 2 (t) − f 1 (t)) H (t − a) + ( f 3 (t) − f 2 (t)) H (t − b)
f (t)
1
1 2 3 t
−1
1
L[ H (t − a)](s) = e−as , s>0
s
Demostración Usando la definición de la transformada de Laplace, tenemos
Z ∞
L[ H (t − a)](s) = H (t − a)e−st dt
Z0 a Z ∞
= 0 · e−st dt + 1 · e−st dt
0 a
Z ∞
= e−st dt
a
t 0<t<1
(
f (t) = 2−t 1<t<2
0 2<t
y
y
t a t
f (t)
1
1 2 t
−1
1 2 −s 1 −2s 1
L[ f ](s) = − e + e = (1 − e − s )2
s s s s
Una de las ventajas de usar la transformada de Laplace, es que nos permite resolver de forma
sencilla ecuaciones diferenciales cuyas funciones entrada son funciones continuas a pedazos.
y 0 + y = f ( t ); y (0) = 0
f (t)
1 2 t
y 0 + y = ( t − 1) H ( t − 1) − ( t − 2) H ( t − 2)
1 1
Y (s) = e−s − 2 e−2s
s2 ( s + 1) s ( s + 1)
Como
1 1 1 1
=− + 2+
s2 ( s + 1) s s s+1
3.4. TRANSFORMADA DE FUNCIONES CONTINUAS A PEDAZOS 111
y ( t ) = − 1 + t + e − t t → t −1 H ( t − 1 ) − − 1 + t + e − t t → t −2 H ( t − 2 )
= t − 2 + e 1− t H ( t − 1 ) − t − 3 + e 2− t H ( t − 2 )
O bien,
0, si 0 < t < 1
y(t) = t − 2 + e 1− t , si 1 < t < 2
1 + e 1− t − e 2− t , si 2 < t
y(t)
0.5
f (t)
π 2π t
−π
-
112 CAPÍTULO 3. TRANSFORMADA DE LAPLACE
f (t) = t · [ H (t) − H (t − π )] + (t − 2π ) · [ H (t − π ) − H (t − 2π )] + 0 · H (t − 2π )
= tH (t) − 2πH (t − π ) − (t − 2π ) H (t − 2π )
1
s2 Y (s) − sy(0) − y0 (0) + Y (s) = − e−πs L[2π ](s) − e−2πs L[t](s)
s2
1 2π −πs 1
( s 2 + 1 )Y ( s ) = 2 − e − 2 e−2πs
s s s
1 1 1
Y (s) = − 2π 2 e−πs − 2 2 e−2πs
s2 ( s2+ 1) s ( s + 1) s ( s + 1)
Descomponiendo en fracciones parciales cada una de las expresiones del lado derecho, se tiene
1 1 1 s −πs 1 1
Y (s) = 2 − 2 − 2π − 2 e − 2
− 2 e−2πs
s s +1 s s +1 s s +1
y(t)
π 2π 3π 4π 5π 6π t
−π
3.4. TRANSFORMADA DE FUNCIONES CONTINUAS A PEDAZOS 113
Ahora bien, si n es un número entero positivo, entonces haciendo el cambio de variable τ = t − nT,
tenemos
Z ( n +1) T Z T Z T
−st −s(τ +nT ) −snT
f (t)e dt = f (τ + nT )e dτ = e f (τ )e−sτ dτ
nT 0 0
Luego,
Z ∞
L[ f ](s) = f (t)e−st dt
0
Z T Z T Z T
= f (t)e−st dt + e−sT
f (t)e−st dt + e−2sT f (t)e−st dt + · · ·
0 0 0
2 Z T
= 1 + e−sT + e−sT + · · · f (t)e−st dt
0
Z T
1
= f (t)e−st dt, s>0
1 − e−sT 0
1, 0≤t<1
n
f (t) = −1, 1≤t<2 ; f ( t + 2) = f ( t )
f (t)
1 2 3 4 5 t
−1
114 CAPÍTULO 3. TRANSFORMADA DE LAPLACE
1
Z 2
L[ f ](s) = f (t)e−st dt
1 − e−2s 0
Z 1 Z 2
1 −st −st
= 1 · e dt + (−1) · e dt
1 − e−2s 0 1
!
1 1 −st 1 1 −st 2
= − e + e
1 − e−2s s 0 s 1
Ejercicios propuestos
Ejer. 1-3 Escriba cada función en términos de la función de Heaviside, luego encuentre su trans-
formada de Laplace
2, 0≤t<3
n
1. f (t) = −2, t≥3
t, 0≤t<2
n
2. f (t) = 0, t≥2
0, 0 ≤ t < 3π/2
3. f (t) = sen t, t ≥ 3π/2
e−2s − 3e−4s
4. F (s) =
s+2
e−3s (s − 5)
5. F (s) =
s2 + 3s + 2
Ejer. 6-8 Resuelva las siguientes ecuaciones diferenciales
sen t , 0 ≤ t < π
n
6. y00 + y = 0 , t≥π , y(0) = 0, y0 (0) = 0
f (t)
2 8 10 16 18 t
10. El tanque mezclador de la figura siguiente contiene inicialmente 500 litros de una solución
salina, con una concentración de sal de 0.02 kg/litro. Durante los primeros 10 minutos de
operación, se abre la válvula A, añadiendo 12 litros/min de una solución con una concentra-
ción de sal de 0.04 kg/litro. Después de 10 minutos, se cambia a la válvula B, la cual agrega
una concentración de 0.06 kg/litro a 12 litro/min. La válvula C elimina 12 litro/min, mante-
niendo constante el volumen. Determine la concentración de sal en el tanque como función
del tiempo.
donde f (t) es una función continua en un intervalo abierto que contiene t = a. La función delta de
Dirac δ(t − a), es representada geométricamente de la siguiente forma
En particular, si consideramos la función f (t) = e−st , obtenemos el siguiente resultado
f (t)
a t
de donde
s−2
Y (s) = e−s
s(s2 + 4s)
s − 2 −s
= 2 e
s ( s + 4)
como
s−2 31 1 1 3 1
= − −
s2 ( s + 4) 8 s 2 s2 8 s + 4
Por tanto,
31 1 1 3 1
Y (s) = − 2− e−s
8s 2s 8s+4
Ejercicios propuestos
Ejer. 1-6 Resuelva los siguientes problemas a valor inicial
1. y0 + y = δ(t − 1); y (0) = 2
2. y00 + 2y0 + y = δ(t − 1); y(0) = 0, y0 (0) = 0
3. y00 + 2y0 + 2y = δ(t − π ); y(0) = 1, y0 (0) = 1
8. Una viga uniforme de longitud L tiene una carga concentrada w0 en x = 21 L (Vea figura 3.5).
Use la transformada de Laplace para resolver la ecuación diferencial
d4 y
1
EI 4 = w0 δ x − L , 0 < x < L,
dx 2
sujeto a las condiciones iniciales y(0) = 0, y0 (0) = 0, y00 (0) = 0 y y000 (0) = 0.
118 CAPÍTULO 3. TRANSFORMADA DE LAPLACE
w0
x
L
y
Figura 3.5
9. Una masa unida a un resorte se libera desde el reposo, un metro por debajo de la posición
de equilibrio para el sistema masa-resorte y comienza a vibrar. Después de π/2 segundos,
la masa es golpeada por un martillo que ejerce un impulso sobre la masa. El sistema queda
descrito mediante el problema simbólico con valores iniciales
d2 x π
+ 9x = − 3δ t − ; x (0) = 1, x 0 (0) = 0,
dt2 2
donde x (t) denota el desplazamiento con respecto a la posición de equilibrio en el instante t.
¿Qué le ocurre a la masa después de ser golpeada?
10. a) Por el método de variación de parámetros muestre que la solución del problema a valor
inicial
y00 + 2y0 + 2y = f (t); y(0) = 0, y0 (0) = 0
es
Z t
y(t) = e−(t−τ ) f (τ ) sen(t − τ )dτ
0
c) Use la transformada de Laplace para resolver el problema a valor inicial con f (t) =
δ(t − π ) y confirme que la solución coincide con el resultado en el inciso (b).
Capı́tulo 4
En este capı́tulo, vamos a estudiar métodos matriciales para la resolución de sistemas de ecuacio-
nes diferenciales lineales de primer orden.
x0 ( t ) = A ( t )x ( t ) + b ( t ), t∈I (4.1)
donde I es un intervalo abierto, x(t), b(t) ∈ Rn y A(t) ∈ Mn×n . Si b(t) = 0, el sistema es llamado
homogéneo, por otro lado, si b(t) 6= 0, el sistema se denomina no homogéneo.
Teorema 4.1 Si x1 (t) y x2 (t) son soluciones del sistema de ecuaciones diferenciales homogéneo
x0 ( t ) = A ( t )x ( t ). (4.2)
Teorema 4.2 Si x1 (t), x2 (t), . . . , xn (t) son soluciones linealmente independientes del sistema ho-
mogéneo 4.2, entonces la solución general del sistema lineal homogéneo 4.2 es
Teorema 4.3 Si x1 (t), x2 (t), · · · , xn (t) es un conjunto fundamental de soluciones del sistema ho-
mogéneo 4.2 y x p (t) es una solución particular del sistema no homogéneo 4.1, entonces
x ( t ) = c 1 x1 ( t ) + c 2 x2 ( t ) + · · · + c n x n ( t ) + x p ( t )
119
120 CAPÍTULO 4. SISTEMAS DE ECUACIONES DIFERENCIALES
Suponga que x(t) = veλt es una solución del sistema 4.3, donde v ∈ Rn y λ ∈ R entonces,
sustituyendo en la ecuación 4.3, tenemos
(veλt )0 = A(veλt )
λveλt = Aveλt =⇒ ( Av − λv)eλt = 0
de donde
Av = λv
x(t) = veλt
es una solución particular del sistema 4.3. Analicemos los siguientes casos
es un conjunto fundamental de soluciones y por tanto, la solución general del sistema 4.3 es dada
por
(5 − λ)(1 − λ) − 3(−1) = 0
λ2 − 6λ + 8 = 0
(λ − 4)(λ − 2) = 0
4.2. SISTEMAS LINEALES HOMOGÉNEOS CON COEFICIENTES CONSTANTES 121
λ1 = 4, λ2 = 2
5 − 4 −1 x 0
3 1−4 y = 0
1 −1 x 0
3 −3 y = 0
que es equivalente a
x−y = 0
1
v1 = 1
5 − 2 −1 x 0
3 1−2 y = 0
3 −1 x 0
3 −1 y = 0
que es equivalente a
3x − y = 0
1
v2 = 3
1 1
x1 ( t ) = 1 e4t , x2 ( t ) = 3 e2t
1 1 1
!
x0 ( t ) = 2 1 −1 x( t )
−8 −5 −3
122 CAPÍTULO 4. SISTEMAS DE ECUACIONES DIFERENCIALES
Solución Los autovalores de la matriz asociada son dadas por la ecuación caracterı́stica
1−λ 1 1
!
det 2 1−λ −1 =0
−8 −5 −3 − λ
(1 − λ)2 (−3 − λ) − 10 + 8 + 8(1 − λ) − 2(−3 − λ) − 5(1 − λ) = 0
(1 − 2λ + λ2 )(−3 − λ) − 2 + 8 − 8λ + 6 + 2λ − 5 + 5λ = 0
−3 − λ + 6λ + 2λ2 − 3λ2 − λ3 + 7 − λ = 0
−λ3 − λ2 + 4λ + 4 = 0
− λ2 ( λ + 1) + 4( λ + 1) = 0
(λ + 1)(4 − λ2 ) = 0
(λ + 1)(2 − λ)(2 + λ) = 0
De donde, los autovalores de la matriz asociada al sistema son
λ1 = −1, λ2 = 2, λ3 = −2
Para hallar un autovector asociado a λ1 = −1, resolvemos el sistema
1 − (−1) 1 1 x 0
! ! !
2 1 − (−1) −1 y = 0
−8 −5 −3 − (−1) z 0
2 1 1 x 0
! ! !
2 2 −1 y = 0
−8 −5 −2 z 0
Que es equivalente a
2 1 1 x 0
! ! !
0 1 −2 y = 0
0 −1 2 z 0
2 0 3 x 0
! ! !
0 1 −2 y = 0
0 0 0 z 0
O bien,
2x + 3z = 0
y − 2z = 0
Tomando z = 2, obtenemos y = 4 y x = −3. Por tanto, un autovector asociado al autovalor λ1 es
−3
!
v1 = 4
2
Para hallar un autovector asociado a λ2 = 2, resolvemos
1−2 1 1 x 0
! ! !
2 1−2 −1 y = 0
−8 −5 −3 − 2 z 0
−1 1 1 x 0
! ! !
2 −1 −1 y = 0
−8 −5 −5 z 0
4.2. SISTEMAS LINEALES HOMOGÉNEOS CON COEFICIENTES CONSTANTES 123
que es equivalente a
−1 1 1 x 0
! ! !
0 1 1 y = 0
0 −13 −13 z 0
−1 0 0 x 0
! ! !
0 1 1 y = 0
0 0 0 z 0
O bien,
−x =0
y+z = 0
0
!
v2 = −1
1
1 − (−2) 1 1 x 0
! ! !
2 1 − (−2) −1 y = 0
−8 −5 −3 − (−2) z 0
3 1 1 x 0
! ! !
2 3 −1 y = 0
−8 −5 −1 z 0
Que es equivalente a
1 −2 2 x 0
! ! !
2 3 −1 y = 0
−8 −5 −1 z 0
1 −2 2 x 0
! ! !
0 7 −5 y = 0
0 −21 15 z 0
x 0
! ! !
1 0 4/7
0 7 −5 y = 0
0 0 0 z 0
o bien,
4
x + z=0
7
7y − 5z = 0
−4
!
v3 = 5
7
124 CAPÍTULO 4. SISTEMAS DE ECUACIONES DIFERENCIALES
−4e−2t
−3e−t 0 !
x(t) = c1 4e−t + c2 −e2t + c3 5e−2t
− 2t
2e t e 7e−2t
Solución Primero debemos hallar la solución general del sistema de ecuaciones. En efecto, los
autovalores de la matriz asociada al sistema de ecuaciones son raı́ces de la ecuación caracterı́stica
−2 − λ 1
det −5 4−λ =0
(−2 − λ)(4 − λ) + 5 = 0
−8 + 2λ − 4λ + λ2 + 5 = 0
λ2 − 2λ − 3 = 0
(λ − 3)(λ + 1) = 0
Por tanto, los autovalores de la matriz asociada al sistema son
λ1 = 3 y λ2 = −1
que es equivalente a
−5x + y = 0
Haciendo y = 5, obtenemos x = 1. . Por tanto, un autovector asociado a λ1 es
1
v1 = 5
que es equivalente a
−x + y = 0
1
v2 = 1
Para hallar la solución al problema a valor inicial, reemplazamos las condiciones iniciales en la
expresión anterior anterior
1 1 1
3 = c1 5 + c2 1
que es equivalente a
1 1 c1 1
5 1 c2 = 3
1 1 c1 1
0 −4 c2 = −2
de donde c2 = 1/2 y c1 = 1/2. Por tanto, la solución del problema a valor inicial es
−t
e + e−t
3t 3t
1 e 1 e 1
x( t ) = + =
2 5e3t 2 e−t 2 5e3t + e−t
Ejercicios propuestos
Ejer. 1-3 Resuelva el sistema lineal dado
1 2
1. x0 = 4 3 x
1 3
2. x0 = 12 1 x
2 10
3. x0 = −1 −5 x
1 −1 4
!
4. x0 = 3 2 −1 x
2 1 −1
−1 1 0
!
5. x0 = 1 2 1 x
0 3 −1
126 CAPÍTULO 4. SISTEMAS DE ECUACIONES DIFERENCIALES
1/2 0 3
6. x0= 1 −1/2 x, x ( 0 ) = 5
6 −3 −10
7. x0 = 2 1 x, x ( 0 ) = −6
0 0 −1 7
! !
0
8. x = 2 0 0 x, x(0) = 5
−1 2 4 5
9. Suponga que
x0 (t) = Ax(t) + b(t) (4.4)
es un sistema de ecuaciones diferenciales lineales, cuya matriz asociada A es diagonalizable,
luego, existe una matriz invertible P tal que
A = PDP−1
donde D es una matriz diagonal, cuyas entradas de la diagonal son los autovalores de la
matriz A. Reemplazando en la ecuación 4.4, tenemos
x0 (t) = PDP−1 x(t) + b(t)
Mediante el cambio de variable y(t) = P−1 x(t), la ecuación anterior puede escribirse como
y0 (t) = Dy(t) + b(t)
donde b(t) = P−1 b(t). Como la matriz D es diagonal, el sistema de ecuaciones obtenido,
puede escribirse en la forma
y10 (t) = λ1 y1 (t) + b1 (t)
y20 (t) = λ2 y2 (t) + b2 (t)
..
.
y0n (t) = λn yn (t) + bn (t),
que es un sistema de ecuaciones diferenciales lineales no acoplados, por tanto, cada una de
las ecuaciones del sistema se pueden resolver separadamente. Una vez resuelto, la solución
del sistema de ecuciones diferenciales es dada por
x(t) = Py(t)
Usando el método descrito anteriormente, resuelva el sistema de ecuaciones
2 −1
e t
0
x ( t ) = 3 −2 x( t ) +
t
10. Como se muestra en la figura siguiente, dos grandes tanques de mezcla conectados A y
B inicialmente contienen 100 litros de salmuera. El lı́quido se bombea dentro y sale de los
tanques como se indica en la figura 4.1; la mezcla bombeada que entra y sale de un tanque se
supone que está bien agitada.
a) Construya un modelo matemático en forma de sistema lineal de ecuaciones diferencia-
les de primer orden para el número de kilogramos x1 (t) y x2 (t) de sal en los tanques A
y B, respectivamente, en el tiempo t. Escriba el sistema en forma matricial.
b) Utilice el método de autovalores de esta sección para resolver el sistema lineal del inciso
a) sujeto a x1 (0) = 20, x2 (0) = 5.
4.2. SISTEMAS LINEALES HOMOGÉNEOS CON COEFICIENTES CONSTANTES 127
es una solución del sistema de ecuaciones 4.5. Usando la fórmula de Euler1 , la solución anterior
puede escribirse como
x(t) = (v1 eαt cos( βt) − v2 eαt sen( βt)) − i (v1 eαt sen( βt) + v2 eαt cos( βt))
Como la combinación lineal de dos soluciones es también una solución de un sistema homogéneo,
concluimos que la parte real e imaginaria de la solución compleja
x1 (t) = v1 eαt cos( βt) − v2 eαt sen( βt), x2 (t) = v1 eαt sen( βt) + v2 eαt cos( βt)
3 −2
x0 ( t ) = 4 −1 x( t )
(3 − λ)(−1 − λ) + 8 = 0
λ2 − 2λ + 5 = 0
λ2 − 2λ = −5
λ2 − 2λ + 12 = −5 + 12
( λ − 1)2 = −4
De donde
λ − 1 = 2i o λ − 1 = −2i
λ = 1 + 2i o λ = 1 − 2i
Es suficiente considerar λ = 1 + 2i. Para hallar un autovector asociado a este autovalor, resolvemos
el sistema
3 − (1 + 2i ) −2 x 0
4 −1 − (1 + 2i ) y = 0
2 − 2i −2 x 0
4 −2 − 2i y = 0
2 − 2i −2 x 0
0 0 y = 0
(2 − 2i ) x − 2y = 0
Luego, la función
1
x( t ) = 1−i e(1+2i)t
4.2. SISTEMAS LINEALES HOMOGÉNEOS CON COEFICIENTES CONSTANTES 129
es una solución compleja del sistema de ecuaciones. Ahora, determinemos la parte real e imagina-
ria de la solución anterior
1 t t
x( t ) = 1 − i (e cos 2t + ie sen 2t)
et cos 2t + iet sen 2t
=
et (cos 2t + sen 2t) + iet (sen 2t − cos 2t)
son soluciones linealmente independientes del sistema. Finalmente, la solución general del sistema
de ecuaciones diferenciales es
cos 2t sen 2t
x(t) = c1 et cos 2t + sen 2t + c2 et sen 2t − cos 2t
1 −1 2
!
x0 ( t ) = −1 1 0 x( t )
−1 0 1
Solución Para hallar los autovalores de la matriz asociada al sistema, resolvemos la ecuación ca-
racterı́stica
1 − λ −1 2
!
det −1 1 − λ 0 =0
−1 0 1−λ
−[−2(1 − λ)] + (1 − λ)[(1 − λ)2 − 1] = 0
(1 − λ)(λ2 − 2λ + 2) = 0
de donde
λ1 = 1, λ2 = 1 + i, λ3 = 1 − i
1 − 1 −1 2 x 0
! ! !
−1 1 − 1 0 y = 0
−1 0 1−1 z 0
0 −1 2 x 0
! ! !
−1 0 0 y = 0
−1 0 0 z 0
que es equivalente a
x =0
−y + 2z = 0
130 CAPÍTULO 4. SISTEMAS DE ECUACIONES DIFERENCIALES
0
!
v1 = 2
1
De donde, la función
0 0
! !
x1 ( t ) = 2 t
e = 2et
1 et
1 − (1 + i ) −1 2 x 0
! ! !
−1 1 − (1 + i ) 0 y = 0
−1 0 1 − (1 + i ) z 0
− i −1 2 x 0
! ! !
−1 − i 0 y = 0
−1 0 −i z 0
que es equivalente a
−1 0 −i x 0
! ! !
−1 − i 0 y = 0
− i −1 2 z 0
−1 0 −i x 0
! ! !
0 −i i y = 0
0 −1 1 z 0
−1 0 −i x 0
! ! !
0 −1 1 y = 0
0 0 0 z 0
−x − iz = 0
−y + z = 0
−i
!
v2 = 1
1
y además,
−i
!
x2 (t) = 1 e (1+ i ) t
1
4.2. SISTEMAS LINEALES HOMOGÉNEOS CON COEFICIENTES CONSTANTES 131
es una solución compleja del sistema de ecuaciones dado. Escribiendo la solución obtenida en la
forma x2 (t) = u1 (t) + iu2 (t), tenemos
−i
!
x2 ( t ) = 1 (et cos t + iet sen t)
1
e sen t − iet cos t
t
(6 − λ)(4 − λ) + 5 = 0
λ2 − 10λ + 29 = 0 =⇒ ( λ − 5)2 = −4
de donde los autovalores son
λ1 = 5 + 2i, λ2 = 5 − 2i
Para calcular un autovector asociado al autovalor λ1 , resolvemos
6 − (5 + 2i ) −1 x 0
5 4 − (5 + 2i ) y = 0
1 − 2i −1 x 0
5 −1 − 2i y = 0
1 − 2i −1 x 0
0 0 y = 0
que es equivalente a
(1 − 2i ) x − y = 0
Haciendo x = 1, obtenemos y = 1 − 2i. De donde
1
v1 = 1 − 2i
132 CAPÍTULO 4. SISTEMAS DE ECUACIONES DIFERENCIALES
es un autovalor asociado a λ1 y
1
x1 (t) = 1 − 2i e(5+2i)t
es una solución compleja del sistema. Escribiendo en la forma x1 (t) = u1 (t) + iu2 (t), tenemos
1 5t 5t
x1 ( t ) = 1 − 2i (e cos 2t + ie sen 2t)
e5t cos 2t + ie5t sen 2t
=
e5t (cos 2t + 2 sen 2t) + ie5t (sen 2t − 2 cos 2t)
e5t cos 2t e5t sen 2t
= +i
e5t (cos 2t + 2 sen 2t) e5t (sen 2t − 2 cos 2t)
Para hallar la solución al problema a valor inicial, reemplazamos las condiciones iniciales
−2 1 0
8 = c 1 1 + c 2 −2
que es equivalente a
1 0 c1 −2
1 −2 c2 = 8
Resolviendo el sistema lineal, obtenemos c1 = −2 y c2 = −5. Por tanto, la solución del problema a
valor inicial es
e5t cos 2t e5t sen 2t
x( t ) = −2 −5
e5t (cos 2t + 2 sen 2t) e5t (sen 2t − 2 cos 2t)
−2 cos 2t − 5 sen 2t
= e5t 8 cos 2t − 9 sen 2t
Ejercicios propuestos
Resuelva los siguientes sistemas de ecuaciones diferenciales:
1 2
1. x0 = −5 −1 x
6 −1
2. x0 = 5 2 x
1 1
3. x0 = −2 −1 x
1 0 0
!
4. x0 = 2 1 −2 x
3 2 1
4.2. SISTEMAS LINEALES HOMOGÉNEOS CON COEFICIENTES CONSTANTES 133
2 5 1
!
5. x0 = −5 −6 4 x
0 0 2
1 −5 1
6. x0 = 1 −3 x, x(0) = 1
4 −5 0
7. x0 = 5 −4 x, x(0) = −1
1 −12 −14 4
! !
8. x0 = 1 2 −3 x, x(0) = 6
1 1 −2 −7
9. El movimiento del sistema masa-resorte que se ilustra en la figura siguiente se describe me-
diante el sistema de segundo orden
m1 x100 = −k1 x1 + k2 ( x2 − x1 ),
m2 x200 = −k2 ( x2 − x1 ) − k3 x2
10. El sistema de los tanques de mezclado que se muestra en la figura siguiente, es un sistema ce-
rrado. Los tanques A, B y C inicialmente contienen el número de litros de salmuera indicados
en la figura.
134 CAPÍTULO 4. SISTEMAS DE ECUACIONES DIFERENCIALES
x0 (t) = Ax(t)
tiene un autovalor λ con multiplicidad k (que se repite k veces). Si para este autovalor podemos
encontrar k autovectores v1 , v2 , . . . , vk . Entonces las funciones
( A − λI )2 v = 0 pero ( A − λI )v 6= 0
x(t) = (v + t( A − λI )v)eλt
es una solución del sistema de ecuaciones. Si mismo asi, no se tiene k soluciones linealmente inde-
pendientes que corresponden al autovalor λ, buscamos vectores v tales que
( A − λI )3 v = 0 pero ( A − λI )2 v 6= 0
4.2. SISTEMAS LINEALES HOMOGÉNEOS CON COEFICIENTES CONSTANTES 135
que es equivalente a
1 −2 1 x 0
! ! !
−2 1 1 y = 0
1 1 −2 z 0
1 −2 1 x 0
! ! !
0 −3 3 y = 0
0 3 −3 z 0
1 0 −1 x 0
! ! !
0 −1 1 y = 0
0 0 0 z 0
o bien,
x −z = 0
−y + z = 0
136 CAPÍTULO 4. SISTEMAS DE ECUACIONES DIFERENCIALES
1
!
v1 = 1
1
es un autovector asociado a λ1 = 2 y
1 e2t
!
x1 ( t ) = 1 e2t = e2t
1 e2t
0 − (−1) 1 1 x 0
! ! !
1 0 − (−1) 1 y = 0
1 1 0 − (−1) z 0
1 1 1 x 0
! ! !
1 1 1 y = 0
1 1 1 z 0
que es equivalente a
x+y+z = 0
1 1
! !
v2 = −1 , v3 = 0
0 −1
Luego,
1 1
! !
−t
x2 (t) = −1 e y x3 (t) = 0 e−t
0 −1
3 −4
x0 ( t ) = 1 −1 x( t )
4.2. SISTEMAS LINEALES HOMOGÉNEOS CON COEFICIENTES CONSTANTES 137
(3 − λ)(−1 − λ) + 4 = 0
−3 − 3λ + λ + λ2 + 4 = 0
λ2 − 2λ + 1 = 0
( λ − 1)2 = 0
3−1 −4 x 0
1 −1 − 1 y = 0
2 −4 x 0
1 −2 y = 0
que es equivalente a
x − 2y = 0
2
v1 = 1
es una solución del sistema de ecuaciones diferenciales. Requerimos calcular una segunda solución
linealmente independiente. Por tanto, calculamos un autovector generalizado v2 , que es solución
del sistema
( A − λI )2 v = 0 con ( A − λI )v 6= 0
2 −4 2 x 0 2 −4 x 0
1 −2 y = 0 con 1 −2 y 6= 0
0 0 x 0
0 0 y = 0 con x − 2y 6= 0
1
v2 = 0
138 CAPÍTULO 4. SISTEMAS DE ECUACIONES DIFERENCIALES
El sistema es equivalente a
−1 1 1 x 0
! ! !
0 1 1 y = 0
0 −1 −1 z 0
−1 0 0 x 0
! ! !
0 1 1 y = 0
0 0 0 z 0
4.2. SISTEMAS LINEALES HOMOGÉNEOS CON COEFICIENTES CONSTANTES 139
que es equivalente a
−x =0
y+z = 0
Tomando z = −1, obtenemos y = 1 y x = 0. Luego,
0
!
v1 = 1
−1
es un autovalor asociado a λ = 2 y
0
! 0 !
x1 (t) = 1 2t
e = e2t
−1 −e2t
es una solución del sistema de ecuaciones diferenciales. Para hallar otras dos soluciones lineal-
mente independientes, determinamos los autovectores generalizados. Consideramos, por tanto,
las soluciones del sistema
( A − λI )2 v = 0 con ( A − λI )v 6= 0
donde A es la matriz asociada al sistema y λ = 2. En efecto,
−1 1 1 2 x 0 −1 1 1 x 0
! ! ! ! ! !
2 −1 −1 y = 0 con 2 −1 −1 y 6= 0
−3 2 2 z 0 −3 2 2 z 0
0 0 0 x 0
! ! !
−1 1 1 y = 0 con − x 6= 0 o y + z 6= 0
1 −1 −1 z 0
Que es equivalente a
x−y−z = 0 con − x 6= 0 o y + z 6= 0
En este caso, podemos considerar el vector
1
!
v2 = 0
1
Luego, una segunda solución del sistema es dada por
Para hallar una tercera solución linealmente independiente, consideramos una solución del siste-
ma
( A − λI )3 v = 0 con ( A − λI )2 v 6= 0
En efecto,
−1 1 1 3 x 0 −1 1 1 2 x 0
! ! ! ! ! !
2 −1 −1 y = 0 con 2 −1 −1 y 6= 0
−3 2 2 z 0 −3 2 2 z 0
0 0 0 x 0
! ! !
0 0 0 y = 0 con x − y − z 6= 0
0 0 0 z 0
x
!
y ∈ R3 con x − y − z 6= 0
z
1
!
v3 = 0
0
Luego
1 2 2
x3 (t) = v3 + t( A − λI )v3 + t ( A − λI ) v3 eλt
2
1 −1 1 1 1 0 0 0 1
" ! ! ! ! !#
1
= 0 +t 2 −1 −1 0 + t2 −1 1 1 0 e2t
0 −3 2 2 0 2 1 −1 −1 0
1 −1 0
" ! ! !#
1
= 0 +t 2 + t2 −1 e2t
0 −3 2 1
de donde
1−t
2t − 1 t2 2t
x3 ( t ) = 2 e
1
−3t + t2
2
Finalmente, la solución general del sistema de ecuaciones es
0
!
1
! 2 − 2t !
x(t) = c1 e2t 1 + c2 e2t t + c3 e2t 4t − t2
−1 1−t −6t + t2
Ejercicios propuestos
Ejer. 1-7 Resuelva el sistema de ecuaciones diferenciales dado
3 −1
1. x0 = 9 −3 x
4.3. SISTEMAS DE ECUACIONES LINEALES NO HOMOGÉNEOS 141
−6 5
2. x0 = −5 4 x
3 −1 −1
!
3. x0 = 1 1 −1 x
1 −1 1
1 1 1
!
4. x0 = 2 1 −1 x
0 −1 1
4 1 0
!
5. x0 = 0 4 1 x
0 0 4
1 0 0
!
6. x0 = 2 2 −1 x
0 1 0
5 −3 −2
!
7. x0 = 8 −5 −4 x
−4 3 3
1 −4 3
8. x0 = 4 −7 x, x(0) = 2
2 4 −1
9. x0 = −1 6 x, x(0) = 6
0 0 1 1
! !
10. x0 = 0 1 0 x, x(0) = 2
1 0 0 5
que es equivalente a
x−y = 0
tomando y = 1, obtenemos x = 1 y
1
v1 = 1
que es equivalente a
x+y = 0
4.3. SISTEMAS DE ECUACIONES LINEALES NO HOMOGÉNEOS 143
1
v2 = −1
e−3t
x2 (t) =
−e−3t
e−t e−3t
xc ( t ) = c 1 + c2
e−t −e−3t
0
Aa + 1 =b
Ab = −a
b = − A −1 a
0
Aa + = − A −1 a
1
0
( A + A −1 )a = − 1
Tenemos
0
( A + A −1 )a = −1
−2 1 −2/3 −1/3 0
1 −2 + −1/3 −2/3 a = −1
−8/3 2/3 0
2/3 −8/3 a = −1
despejando a, tenemos
−1
−8/3 2/3 0
a= 2/3 −8/3 −1
1 −8/3 −2/3 0
= 60 −2/3 −8/3 −1
9
3 2/3
=
20 8/3
1 1
=
10 4
Reemplazando el resultado obtenido en la ecuación b = − A−1 a, tenemos
1 −2 −1 1 1 1 6 1 2
b=− = =
3 −1 −2 10 4 30 9 10 3
Reemplazando los valores obtenidos, obtenemos la solución particular x p (t)
1 1 1 2 1 cos t + 2 sen t
x p (t) = cos t + sen t =
10 4 10 3 10 4 cos t + 3 sen t
Finalmente, la solución general del sistema homogéneo es
x ( t ) = xc ( t ) + x p ( t )
−t −3t
e e 1 cos t + 2 sen t
= c1 + c 2 +
e−t −e−3t 10 4 cos t + 3 sen t
Ejemplo 4.11 Halle la solución general del sistema de ecuaciones diferenciales
2 −1
e t
0
x ( t ) = 3 −2 x( t ) +
t
Solución Primero, hallamos la solución general del sistema homogéo asociado
2 −1
x0 ( t ) = 3 − 2 x ( t )
Los autovalores de sus matriz asociada son obtenidas a partir de la su ecuación caracterı́stica
2−λ −1
det 3 −2 − λ = 0
(2 − λ)(−2 − λ) + 3 = 0
−4 − 2λ + 2λ + λ2 + 3 = 0
λ2 − 1 = 0
(λ + 1)(λ − 1) = 0
4.3. SISTEMAS DE ECUACIONES LINEALES NO HOMOGÉNEOS 145
λ1 = −1 y λ2 = 1
que es equivalente a
3x − y = 0
1
v1 = 3
es un autovector asociado a λ1 y
e−t
x1 (t) =
3e−t
que es equivalente a
x−y = 0
1
v2 = 1
es un autovector asociado a λ2 y
et
x2 (t) =
et
es la solución general del sistema de ecuaciones homogéno que lo denominamos solución comple-
mentaria. Como la parte no homogénea del sistema se puede escribir en la forma
e t 0
1 t
b( t ) = = 1 t + 0 e
t
x p (t) = a + bt + (c + dt)et
para alguna constante k a ser determinada. Reemplazando esta expresión en la tercera ecuación,
tenemos
1 1
Ac + 0 = c + k 1
k−1
Ac − Ic = k
k−1
( A − I )c = k
1 −1 c1 k−1
3 −3 c = k
2
1
c1 − c2 =
2
en este último caso, podemos considerar c1 = 1/2 y c2 = 0, con esto, tenemos que
1/2 3/2
c= 0 y d= 3/2
0
Ab = −1 luego,
0
b = A −1 − 1
1 −2 1 0
=
−1 −3 2 −1
1
= 2
1
Aa = 2 luego,
1
a = A −1 2
1 −2 1 1
=
−1 −3 2 2
0
= −1
0
1 1 1 3 t
x p (t) = −1 + 2 t + 0 + 3 t e
2
1 2t + (1 + 3t) et
=
2 −2 + 4t + 3tet
x ( t ) = xc ( t ) + x p ( t )
−t t
e e 1 2t + (1 + 3t) et
= c1 + c2 +
3e−t et 2 −2 + 4t + 3tet
Ejercicios propuestos
Ejer. 1-6 Usando el método de coeficientes indeterminados, resuelva el sistema lineal no homogéneo
2 3 −7
1. x0 = −1 −2 x + 5
148 CAPÍTULO 4. SISTEMAS DE ECUACIONES DIFERENCIALES
3 −4 1
2. x0 = 1 −1 x+ 1 et
1 3−2t2
3. x0 = 3 1 x+
t+5
−1 5 sen t
4. x0 = −1 1 x + −2 cos t
1 1 1 1
! !
5. x0 = 0 2 3 x + −1 e4t
0 0 5 2
1 3 2 sen t
! !
6. x0 = −1 2 1 x+ 0
4 −1 −1 0
Ejer. 7-9 Resuelva el problema a valor inicial
−1 −2 3 −4
7. x0 = 3 4 x + 3 e −t , x ( 0 ) = 5
1 −1 t 3
8. x0 = 1 3 x + t+1 e −t , x ( 0 ) = 2
−1 −1 −2 1 0
! ! !
9. x0 = 1 1 1 x+ 0 e−t , x(0) = 0
2 1 3 0 0
10. Considere los grandes tanques de mezclado que se muestran en la figura 4.4. Suponga que
ambos tanques A y B inicialmente contienen 100 litros de salmuera. El lı́quido se bombea
dentro y se saca de los tanques como se indica en la figura; se supone que la mezcla bombeada
que entra y sale está bien agitada.
A B
Además, suponga que x1 (t), x2 (t), . . . , xn (t) es un conjunto fundamental de soluciones del sistema
homogéneo asociado a 4.6. Entonces la solución general del sistema homogéneo (solución comple-
mentaria)
x0 (t) = Ax(t)
x c ( t ) = c 1 x1 ( t ) + c 2 x2 ( t ) + · · · + c n x n ( t )
c1
| | ··· |
!
c2
= x1 (t) x2 (t) · · · xn (t) ..
| | ··· |
.
cn
= Φ ( t )c
Φ0 (t) = AΦ(t)
Definición 4.2 La matriz Φ(t) cuyas columnas forman un conjunto fundamental de soluciones del
sistema homogéneo es llamada matriz fundamental de soluciones
Para hallar una solución particular del sistema no homogéneo 4.6, planteamos
x p ( t ) = Φ ( t )c( t )
Φ ( t )c0 ( t ) = b ( t )
Como las columnas de la matriz fundamental son linealmente independientes, esta matriz es in-
vertible, luego
−1
c0 (t) = [Φ(t)] b( t )
x ( t ) = xc ( t ) + x p ( t )
Z
= Φ ( t )c + Φ ( t ) [Φ(t)]−1 b(t)dt
El método descrito arriba para hallar la solución particular es denominado método de variación
de parámetros.
t −3
4 −2
0
x ( t ) = 8 −4 x( t ) + , t>0
− t −2
Solución Primero, hallemos la matriz fundamental asociado al sistema homogéneo
4 −2
x0 ( t ) = 8 − 4 x ( t )
(4 − λ)(−4 − λ) + 16 = 0
−16 − 4λ + 4λ + λ2 + 16 = 0
λ2 = 0
Por tanto, λ = 0 es el único autovalor, y además tiene multiplicidad 2. Para hallar los autovectores
asociados, resolvemos
4−0 −2 x 0
8 −4 − 0 y = 0
4 −2 x 0
8 −4 y = 0
que es equivalente a
2x − y = 0
es un autovector asociado a λ = 0 y
1 1
x1 (t) = 2 e0t = 2
4.3. SISTEMAS DE ECUACIONES LINEALES NO HOMOGÉNEOS 151
es una solución del sistema de ecuaciones. Para hallar la segunda solución, determinamos un au-
tovector generalizado v2 dado por el sistema condicionado
( A − λI )2 v2 = 0 con ( A − λI )v2 6= 0
donde A es la matriz asociada al sistema y λ el autovalor de A. Reemplazando, tenemos
4 −2 2 x 0 4 −2 x 0
8 −4 y = 0 con 8 −4 y 6= 0
0 0 x 0
0 0 y = 0 con 2x − y 6= 0
x
y ∈ R con 2x − y 6= 0
2
x p ( t ) = Φ ( t )c( t )
Luego,
t −3
Z
−4t 1/2 + 2t
c( t ) = dt
4 −2 − t −2
−4t−2 − 1/2t−2 − 2t−1
Z
= dt
4t−3 + 2t−2
9/2t−1 − 2 ln t
=
−2t−2 − 2t−1
Luego, la solución particular es
9/2t−1 − 2 ln t
1 1/4 + t
x p (t) = 2 2t −2t−2 − 2t−1
9/2t−1 − 2 ln t − 1/2t−2 − 1/2t−1 − 2t−1 − 2
=
9t−1 − 4 ln t − 4t−1 − 4
−2 + 2t−1 − 1/2t−2 − 2 ln t
=
−4 + 5t−1 − 4 ln t
Finalmente, la solución general del sistema no homogéneo es
x ( t ) = xc ( t ) + x p ( t )
−2 + 2t−1 − 1/2t−2 − 2 ln t
1 1/4 + t c1
= 2 c2 +
2t −4 + 5t−1 − 4 ln t
−2 + 2t−1 − 1/2t−2 − 2 ln t
c1 + 1/4c2 + tc2
= +
2c1 + 2c2 t −4 + 5t−1 − 4 ln t
Ejemplo 4.13 Resuelva el sistema de ecuaciones diferenciales
1 −1
− t 2
0
x (t) = 1 3 x( t ) + 2t
Solución Primero hallemos un conjunto fundamental de soluciones del sistema homogéneo
1 −1
x0 ( t ) = 1 3 x( t )
En efecto, los autovalores de la matriz asociada son soluciones de
1 − λ −1
det 1 3−λ =0
(1 − λ)(3 − λ) + 1 = 0
3 − λ − 3λ + λ2 + 1 = 0
λ2 − 4λ + 4 = 0
( λ − 2)2 = 0
De donde λ = 2 es el único autovalor de la matriz, pero con multiplicidad 2. Para hallar los
autovectores asociados, resolvemos
1 − 2 −1 x 0
1 3−2 y = 0
−1 −1 x 0
1 1 y = 0
4.3. SISTEMAS DE ECUACIONES LINEALES NO HOMOGÉNEOS 153
que es equivalente a
x+y = 0
1
v1 = −1
es un autovector asociado a λ = 2, por tanto,
e2t
x1 ( t ) =
−e2t
es una solución del sistema homogéneo. Para hallar una segunda solución, determinamos un au-
tovector generalizado v2 que es solución del sistema condicionado
( A − λI )2 v2 = 0 con ( A − λI )v2 6= 0
Donde A es la matriz asocidada al sistema y λ es el autovalor repetido. Luego
−1 −1 2 x 0 −1 −1 x 0
1 1 y = 0 con 1 1 y 6= 0
0 0 x 0
0 0 y = 0 con x + y 6= 0
x
y ∈ R con x + y 6= 0
2
1
v2 = 0
Por otro lado, como solución particular del sistema no homogéneo, planteamos
e2t (1 − t)e2t
x p (t) = c( t )
−e2t te2t
Como
1 te2t −(1 − t)e2t t t−1
−1
Φ (t) = 4t = e−2t
e e2t e2t 1 1
tenemos
2
t t−1
Z
− t
−2t
c( t ) = e 1 1 dt
2t
−t3 + 2t2 − 2t
Z
= e−2t dt
−t2 + 2t
(−t3R+ 2t2 − 2t)e−2t dt
R
=
(−t2 + 2t)e−2t dt
Integrando por partes cada uno de los componentes del vector c(t), obtenemos
1 4t3 − 2t2 + 6t + 3
c(t) = e−2t
8 4t2 − 4t − 2
x ( t ) = xc ( t ) + x p ( t )
2t 1
e (1 − t)e2t c1 6t2 + 4t + 1
= c2 + 8
−e2t te2t −2t2 − 8t − 3
4.4. SISTEMAS DE ECUACIONES LINEALES Y LA TRANSFORMADA DE LAPLACE 155
Ejercicios propuestos
Ejer.1-7 Utilice variación de parámetros para resolver el sistema dado
3 −3 4
1. x0 = 2 −2 x + −1
2 −1 0
2. x0 = 3 −2 x+ 4t
1 −2 tan t
3. x0 = 1 −1 x+ 1
0 1 0
4. x0 = −1 0 x+
sec t tan t
0 −1 sec t
5. x0 = 1 0 x + 0
e−2t
2 −2 1
6. x0 = 8 −6 x+ 3 t
3 −1 −1 0
! !
7. x0 = 1 1 −1 x+ t
1 −1 1 2et
Ejer. 8-10 Utilice variación de parámetros para resolver el problema a valor inicial
1 −1 1/t 2
8. x0 = 1 −1 x + 1/t , x(1) = −1
3 −4 1 1
9. x0 = 1 −1 x+ 1 et , x(0) = 1
2 0 1 1 1
! ! !
10. x0 = 0 2 0 x+ 0 e2t , x(0) = 1
0 1 3 1 1
donde x(t), b(t) ∈ Rn y A ∈ Mn,n . Aplicando la transformada de Laplace a cada una de las
ecuaciones lineales, obtenemos el sistema lineal
Si s > max{|λ| : det(λI − A) = 0}, entonces la matriz sI − A es invertible y por tanto, podemos
despejar X(s):
Aplicando la transformada inversa a cada uno de los componentes de la función X(s), obtenemos
la solución del sistema de ecuaciones diferenciales 4.7
0 −1 t 1
x0 ( t ) = 1 0 x( t ) + 1 ; x(0) = −1
1
0 −1
sX(s) − x(0) = X ( s ) + s2
1 0 1
s
1
1 0 0 −1 1
s X( s ) − X( s ) = −1 + s2
0 1 1 0 1
s
!
s2 +1
s 1
X( s ) = s2
−1 s 1− s
s
cos t + sen t
x( t ) = t − cos t + sen t
Finalmente, la solución del problema a valor inicial se obtiene tomando la transformada inversa
de X(s)
1 1 t 3 −t 5 3t
− 3 − 4 e + 8 e + 24 e
x( t ) = 4 3 −t 5 3t
3 − 4e + 12 e
Ejercicios propuestos
Resuelva los siguientes problemas a valor inicial mediante la transformada de Laplace
1 −3 0
1. x0 = −2 2 x, x(0) = 5
1 −1 1
2. x0 = 5 −3 x, x(0) = 2
3 −2
t
2
3. x0 = x+ , x(0) =
2 −2 3et 1
21 1 1
4. x0 = −1 e
4 1
t, x (x+
0 ) = 1
4 5 t 1
5. x0 = −2 −2 x + 4e cos t , x(0) = 1
0
158 CAPÍTULO 4. SISTEMAS DE ECUACIONES DIFERENCIALES
3 −2 1 − H (t − π ) 0
6. x0= 2 −2 x+ 0 , x ( 0 ) = 0
2 −2 0 1
0
7. x = 4 −2 x + δ(t − π ) , x(0) = 0
1 2 −3 1
! !
0
8. x = 1 1 2 x, x(0) = 0
1 −1 4 0
2 0 1 1 1
! ! !
9. x0 = 0 2 0 x + 0 e2t , x(0) = 1
0 0 3 1 1
1 0 0 0 0
! ! !
10. x0 = 2 1 −2 x + 0 , x(0) = 1
3 2 1 et cos 2t 1
Capı́tulo 5
159
160 CAPÍTULO 5. ECUACIONES DIFERENCIALES PARCIALES
Ejemplo 5.2 Aplique el método de separación de variables para hallar las soluciones producto de
la ecuación diferencial parcial
ut = u xx + u
Solución Reemplazando u( x, t) = X ( x ) T (t) en la ecuación diferencial, tenemos
T0 X 00 + X
XT 0 = X 00 T + XT =⇒ = =λ
T X
donde λ es una constante. De lo anterior, obtenemos
T 0 = λT, X 00 + (1 − λ) X = 0
Para resolver las ecuaciones ordinarias, analicemos los siguientes casos:
Si λ = 1, se tiene
T0 = T =⇒ T ( t ) = c1 e t
X 00 = 0 =⇒ X ( x ) = c2 + c3 x
Por tanto
u( x, t) = et (C1 + C2 x )
son soluciones de la ecuación diferencial.
Si λ > 1, podemos considerar λ = 1 + α2 , donde α es una constante positiva, entonces
2 )t
T 0 = (1 + α2 ) T =⇒ T ( t ) = c 1 e (1+ α
X 00 − α2 X = 0 =⇒ X ( x ) = c2 cosh(αx ) + c3 senh(αx )
Por tanto
2
u( x, t) = e(1+α )t (C1 cosh(αt) + C2 senh(αt))
son soluciones de la ecuación diferencial.
Finalmente, si λ = 1 − α2 , donde α es una constante positiva, entonces
2 )t
T 0 = (1 − α2 ) T =⇒ T ( t ) = c 1 e (1− α
X 00 + α2 X = 0 =⇒ X ( x ) = c2 cos(αx ) + c3 sen(αx )
Por tanto
2
u( x, t) = e(1−α )t (C1 cos(αx ) + C2 sen(αx ))
tambien, son soluciones de la ecuación diferencial parcial.
5.1. SEPARACION DE VARIABLES 161
Ejercicios propuestos
En los ejercicios 1 a 5, utilice separación de variables para encontrar soluciones producto para la
ecuación diferencial parcial dada
1. xu x = yuy
2. yu xy + u = 0
∂u ∂2 u
3. = k 2, k>0
∂t ∂x
∂2 u 2
2∂ u
4. = a
∂t2 ∂x2
∂2 u ∂2 u
5. + 2 =0
∂x2 ∂y
6. Utilice el método de separación de variables para resolver el problema con valores iniciales
ut = uy + u, u(0, y) = 2e−y − e2y
7. Demuestre que la ecuación diferencial parcial dada tiene la solución producto indicada
∂2 u 1 ∂u 1 ∂2 u
+ + =0
∂r2 r ∂r r2 ∂θ 2
u(r, θ ) = (c1 cos αθ + c2 sen αθ )(c3 r α + c4 r −α )
8. Demuestre que la ecuación diferencial parcial dada tiene la solución producto indicada
2
∂u ∂ u 1 ∂u
=k +
∂t ∂r2 r ∂r
2
u(r, t) = e−kα t (c1 J0 (αr ) + c2 Y0 (αr ))
1 L nπ
Z
an = f ( x ) cos x dx, n = 0, 1, . . .
L −L L
1 L nπ
Z
bn = f ( x ) sen x dx, n = 1, 2, . . .
L −L L
f ( x ) = x, −π < x < π
Solución Los coeficientes de Fourier son dados por:
1
Z π
an = x cos(nx )dx
π −π
= 0, n = 0, 1, . . .
1 π
Z
bn = x sen(nx )dx
π −π
2
Z π
= x sen(nx )dx
π 0
π
2 x 1 π
Z
= − cos(nx ) + cos(nx )dx
π n 0 n 0
π
2 π 1
= − cos(nπ ) + 2 sen(nx )
π n n 0
2
= − (−1)n , n = 1, 2, . . .
n
Por tanto, la serie de Fourier de f ( x ) es
∞
2
f (x) ∼ ∑ n
(−1)n+1 sen(nx )
n =1
2 1
= 2 sen x − sen(2x ) + sen(3x ) − sen(4x ) + · · ·
3 2
Ejemplo 5.4 Calcule la serie de Fourier de la función
0, −π < x < 0
n
f ( x ) = x, 0<x<π
5.2. SERIES DE FOURIER 163
1 π
Z
an = f ( x ) cos(nx )dx
π −π
1 π
Z
= x cos(nx )dx
π 0
π
1 x 1 π
Z
= sen(nx ) − sen(nx )dx
π n 0 n 0
π
1 1
= 0 + 2 cos(nx )
π n 0
(−1)n − 1
= , n = 1, 2, . . .
πn2
1 π
Z
a0 = f ( x )dx
π −π
1
Z π
= xdx
π 0
π
=
2
1 π
Z
bn = f ( x ) sen(nx )dx
π −π
1
Z π
= x sen(nx )dx
π 0
π
1 x 1 π
Z
= − cos(nx ) + cos(nx )dx
π n 0 n 0
1
= − (−1)n , n = 1, 2, . . .
n
∞
(−1)n − 1
π 1
f (x) ∼ + ∑ 2
cos(nx ) + (−1)n+1 sen(nx )
4 n =1 πn n
2 1 − (−1)2k
n = 2k =⇒ bn = b2k = =0
π 2k
2 1 − (−1)2k−1 4 1
n = 2k − 1 =⇒ bn = b2k−1 = =
π 2k − 1 π 2k − 1
Luego, la serie de Fourier de f ( x ) es
∞ ∞
2 1 − (−1)n 4 1
f (x) ∼ ∑π n
sen ( nx ) = ∑ π 2k − 1
sen(2k − 1) x
n =1 k =1
El siguiente teorema da las condiciones suficientes para la convergencia de una serie de Fourier1
Teorema 5.1 (Convergencia de las series de Fourier) Si f y f 0 son continuas por partes en el inter-
f (x+ ) + f (x− )
valo [− L, L], entonces para cualquier x ∈ (− L, L), la serie de Fourier converge a .
2
f (− L+ ) + f ( L− )
Para x = ± L, la serie converge a .
2
En particular, si f es continua en x, entonces f ( x + ) = f ( x − ) = f ( x ), luego
∞
a0 nπ nπ
+ ∑ an cos x + bn sen x = f (x)
2 n =1
L L
2n − 1
como sen π = (−1)n+1 , la ecuación anterior puede escribirse como
2
1 1 1 π
1− + − +··· =
3 5 7 4
donde
1 L nπ
Z
an = f P ( x ) cos x dx
L −L L
2 L nπ
Z
= f ( x ) cos x dx, n = 0, 1, . . .
L 0 L
Como f P ( x ) = f ( x ) en el intervalo (0, L), entonces
∞
a0 nπ
f (x) ∼ + ∑ an cos x , x ∈ (0, L) (5.2)
2 n =1
L
La serie trigonométrica 5.2, se denomina serie de cosenos de Fourier de f en el intervalo (0, L).
166 CAPÍTULO 5. ECUACIONES DIFERENCIALES PARCIALES
Por otro lado, tambien podemos extender la función f sobre el intervalo (− L, L) de la siguiente
forma (extensión impar)
− f (− x ), −L < x < 0
f I (x) =
f ( x ), 0<x<L
donde
1 L nπ
Z
bn = f I ( x ) sen x dx
L −L L
2 L nπ
Z
= f ( x ) sen x dx, n = 1, 2, . . .
L 0 L
La serie trigonométrica 5.3, se denomina serie de senos de Fourier de f en el intervalo (0, L).
f ( x ) = π − x, 0<x<π
donde
π
2 2 ( π − x )2
Z π
a0 = (π − x )dx = − =π
π 0 π 2 0
2 π
Z
an = (π − x ) cos(nx )dx
π 0
π
2 π−x
1 π
Z
= sen(nx ) + sen(nx )dx
π n 0 n 0
π
2 cos(nx )
= −
π n2 0
2 1 − (−1)n
= , n = 1, 2, . . .
π n2
5.2. SERIES DE FOURIER 167
Luego
∞
π 2 1 − (−1)n
f (x) = +∑ cos(nx )
2 n =1 π n2
∞
π 4 1
= +∑ cos(2n − 1) x
2 n=1 π (2n − 1)2
π 4 1 1 1
= + cos x + cos 3x + cos 5x + cos 7x + · · ·
2 π 9 25 49
f (x)
1 2 x
x, 0<x<1
n
f (x) = 1, 1≤x<2
∞
a0 nπ
f (x) = + ∑ an cos x
2 n =1
2
donde
Z 2 Z 1 Z 2
2 3
a0 = f ( x )dx = xdx + (1)dx =
2 0 0 1 2
Z 2
2 πn
an = f ( x ) cos
x dx
2 0 2
Z 1 πn Z 2 πn
= x cos x dx + cos x dx
0 2 1 2
2 nπ 4 nπ 2 nπ
= sen + 2 2 cos −1 − sen
nπ 2 n π 2 nπ 2
4 nπ
= 2 2 cos −1
n π 2
168 CAPÍTULO 5. ECUACIONES DIFERENCIALES PARCIALES
f ( x ) = x (1 − x ), 0<x<1
donde
Z 1
2
bn = x (1 − x ) sen(nπx )dx
1 0
1 1
x (1 − x )
1
Z
=2 − cos(nπx ) + (1 − 2x ) cos(nπx )dx
nπ nπ 0
0
1 Z 1
1 − 2x
2
2
= sen ( nπx ) + sen ( nπx ) dx
nπ nπ 0 nπ 0
1
!
4 cos(nπx )
= 2 2 −
n π nπ 0
4 1 − (−1)n
=
π3 n3
Luego
∞
4 1 − (−1)n
f (x) = ∑ π3 n3
sen(nπx )
n =1
∞
8 1
= ∑
π (2n − 1)3
3
sen(2n − 1)πx
n =1
8 1 1
= 3 sen(πx ) + sen(3πx ) + sen(5πx ) + · · ·
π 27 125
1, 0<x<1
n
f ( x ) = 2 − x, 1≤x<2
5.2. SERIES DE FOURIER 169
f (x)
1 2 x
donde
Z 2 nπ
2
bn = f ( x ) sen x dx
2 0 2
Z 1 nπ Z 2 nπ
= sen x dx + (2 − x ) sen x dx
0 2 1 2
2 nπ 2 nπ 2 2 nπ
=− cos −1 + cos + sen
nπ 2 nπ 2 nπ 2
2
2 2 nπ
= + sen
nπ nπ 2
Luego
∞ 2 !
2 2 nπ nπ
f (x) = ∑ nπ
+
nπ
sen
2
sen
2
x
n =1
Ejercicios propuestos
En los problemas 1 a 3, calcule la serie de Fourier de la función dada
1. f ( x ) = | x |, −π < x < π
x + l, −l < x < 0
2. f ( x ) = x − l, 0<x<l
1, −2 < x < 0
n
3. f ( x ) = x, 0<x<2
En los problemas 4 a 6, calcule la serie de senos de Fourier de la función dada
4. f ( x ) = −1, 0<x<1
5. f ( x ) = cos x, 0<x<π
170 CAPÍTULO 5. ECUACIONES DIFERENCIALES PARCIALES
x, 0 < x < π/2
6. f ( x ) = π − x, π/2 < x < π
ut = ku xx (5.4)
u(0, t) = 0, u( L, t) = 0, t > 0
u( x, 0) = f ( x ), 0 < x < L.
T0 X 00
XT 0 = kX 00 T =⇒ = = −λ
kT X
de donde obtenemos dos ecuaciones diferenciales ordinarias
X 00 + λX = 0
T 0 + λkT = 0
X (0) = 0, X ( L) = 0
5.3. ECUACIÓN DEL CALOR 171
2t
T 0 + k (nπ/L))2 T = 0 =⇒ T (t) = c1 e−k(nπ/L)
Por tanto, las funciones
2
nπ
un ( x, t) = e−k(nπ/L) t sen x , n = 1, 2, . . .
L
forman un conjunto fundamental de soluciones que satisfacen la ecuación diferencial y las condi-
ciones de frontera. Para hallar una solución que satisfaga la condición inicial, suponemos que la
solución tiene la forma
∞ nπ
u( x, t) = ∑ bn e−k(nπ/L) t sen
2
x
n =1
L
de donde bn son los coeficientes de Fourier del desarrollo en senos de la función f ( x ), es decir
Z L nπ
2
bn = f ( x ) sen x dx
L 0 L
172 CAPÍTULO 5. ECUACIONES DIFERENCIALES PARCIALES
2
Z π
bn = x2 sen(nx )dx
π 0
2 (2 − π 2 n2 )(−1)n − 2
bn = , n = 1, 2, . . .
π n3
Luego, la solución del problema con valores iniciales y en la frontera es
∞
2 (2 − π 2 n2 )(−1)n − 2 −n2 t
u( x, t) = ∑π n3
e sen(nx )
n =1
X 0 (0) = 0, X 0 (1) = 0
Si λ = 0, se tiene
X 00 = 0 =⇒ X ( x ) = c1 x + c2
X ( x ) = c2
son soluciones del problema 5.6. Por otro lado, si λ = 0, reemplazando y resolviendo la ecuación
diferencial T 0 + 2λT = 0, obtenemos
T0 = 0 =⇒ T ( t ) = c3
X 00 + α2 X = 0 =⇒ X ( x ) = c1 cos αx + c2 sen αx
Considerando las condiciones en la frontera X 0 (0) = X 0 (1) = 0, las soluciones son no triviales si
α = nπ, y en este caso
X ( x ) = c1 cos(nπx )
X 00 − α2 X = 0 =⇒ X ( x ) = c1 cosh(αx ) + c2 senh(αx )
Considerando las condiciones en la frontera X 0 (0) = X 0 (1) = 0, concluimos que para este caso,
solo se tiene como solución la solución nula.
Finalmente, el conjunto fundamental de soluciones de la ecuación diferencial con condiciones en
la frontera está formada por las funciones
2
u0 ( x, t) = 1, un ( x, t) = cos(nπx )e−2(nπ ) t , n = 1, 2, . . .
Ejercicios propuestos
1. Resuelva el problemas con valores iniciales y en la frontera
ut = ku xx , 0 < x < L, t > 0
u(0, t) = 0, u( L, t) = 0, t > 0
u( x, 0) = x ( L − x ), 0 < x < L
2. Determine una solución formal del problema con valores iniciales y en la frontera
ut = 4u xx , 0 < x < π, t > 0
u x (0, t) = 0, u(π, t) = 0, t > 0
u( x, 0) = x, 0 < x < π
3. Determine una solución formal del problema de flujo de calor descrito por el problema con
valores iniciales y en la frontera
ut = 3u xx , 0 < x < π, t > 0
u x (0, t) = 0, u x (π, t) = 0, t > 0
u( x, 0) = x, 0 < x < π
5. Determine una solución formal del problema con valores iniciales y en la frontera
ut = u xx , 0 < x < π, t > 0
u(0, t) = 0, u(π, t) + u x (π, t) = 0, t > 0
u( x, 0) = f ( x ), 0 < x < π
5.4. ECUACIÓN DE LA ONDA 175
u(0, t) = u( L, t) = 0, t > 0,
u( x, 0) = f ( x ), ut ( x, 0) = 0 0 < x < L,
T 00 X 00
XT 00 = a2 X 00 T =⇒ = = −λ
a2 T X
de donde obtenemos dos ecuaciones diferenciales ordinarias
X 00 + λX = 0, 0<x<L
00 2
T + λa T = 0, t>0
X (0) = 0, X ( L) = 0
T 0 (0) = 0
X 00 + λX = 0, X (0) = X ( L ) = 0 (5.8)
En la sección anterior, se demostró que este problema tiene soluciones no triviales si λ = (nπ/L)2 .
Luego
nπ
X ( x ) = c sen x , n = 1, 2, . . .
L
son soluciones del problema con valores en la frontera 5.8. Ahora, procedemos a resolver el pro-
blema
T 00 + (nπ/L)2 a2 T = 0, T 0 (0) = 0 (5.9)
Para encontrar la solución del problema que satisface la condición inicial u( x, 0) = f ( x ), supone-
mos que la solución tiene la forma
∞ nπa nπ
u( x, t) = ∑ bn cos L
t sen
L
x (5.10)
n =1
∞ nπ
u( x, 0) = ∑ bn sen
L
x = f (x)
n =1
de donde bn son los coeficientes de Fourier del desarrollo en senos de la función f ( x ), es decir
Z L nπ
2
bn = f ( x ) sen x dx
L 0 L
Ejemplo 5.13 Determine la solución del problema de la cuerda vibrante descrita por el problema
con valores iniciales y en la frontera
∞
u( x, t) = ∑ bn cos(nπt) sen(nπx)
n =1
∞
4 1 − (−1)n
u( x, t) = ∑ π3 n3
cos(nπt) sen(nπx )
n =1
∞
8 1
=
π3 ∑ (2n − 1)3
cos ((2n − 1)πt) sen ((2n − 1)πx )
n =1
5.4. ECUACIÓN DE LA ONDA 177
u(0, t) = 0, u(π, t) = 0,
u( x, 0) = 0, ut ( x, 0) = sen x,
X (0) = 0, X (π ) = 0, T (0) = 0
Como en el ejemplo anterior, las soluciones del problema a valor con valores en la frontera
X 00 + λX = 0, X (0) = X ( π ) = 0
X ( x ) = c sen(nx ), n = 1, 2, . . .
T 00 + n2 a2 T = 0, T (0) = 0
forman un conjunto fundamental de soluciones de la ecuación diferencial que satisfacen las condi-
ciones en la frontera
u(0, t) = 0, u(π, t) = 0, u( x, 0) = 0, 0 < x < π, 0 < t
Para encontrar una solución del problema que satisfaga la condición inicial ut ( x, 0) = sen x, plan-
teamos
∞
u( x, t) = ∑ cn sen(nat) sen(nx)
n =1
178 CAPÍTULO 5. ECUACIONES DIFERENCIALES PARCIALES
En este caso, no es necesario encontrar la serie de senos de Fourier, pues por comparación
1
c1 =
a
cn = 0, n = 2, 3, . . .
Por tanto, la solución al problema con condiciones iniciales y en la frontera es
1
u( x, t) = sen( at) cos( x )
a
Ejercicios propuestos
1. Resuelva la ecuación de la onda
utt = a2 u xx , 0 < x < 1, t>0
sujeta a las condiciones
u(0, t) = 0, u(1, t) = 0, t > 0
u( x, 0) = 0, ut ( x, 0) = x (1 − x ), 0<x<1
4. Una cuerda vibrante queda descrita mediante el problema con valores iniciales y en la fron-
tera 5.12.
utt = α2 u xx , 0 < x < L, t > 0
u(0, t) = 0, u( L, t) = 0, t > 0 (5.12)
u( x, 0) = f ( x ), ut ( x, 0) = g( x ), 0 < x < L
u xx + uyy = 0 (5.13)
en el rectángulo 0 < x < a, 0 < y < b, y que satisfaga las condiciones de frontera (problema de
Dirichlet)
u( x, 0) = 0, u( x, b) = 0, 0 < x < a,
u(0, y) = 0, u( a, y) = f (y), 0 < y < b
Y 00 + λY = 0, Y (0) = Y ( b ) = 0 (5.14)
X 00 − (nπ/b)2 X = 0, X (0) = 0
u( x, 0) = u( x, b) = 0, 0 < x < a
u(0, y) = 0, 0 < y < b
Para encontrar una solución que satisfaga la condición u( a, y) = f (y), suponemos que la solución
tiene la forma
∞ ∞ nπ nπ
u( x, y) = ∑ cn un (x, y) = ∑ cn senh b
x sen
b
y
n =1 n =1
donde
y, 0<y<1
n
f (y) = 2 − y, 1≤y<2
Solución Por lo expuesto arriba, la solución a la ecuación de Laplace con condiciones en la frontera
es dada por
∞ nπ nπ
u( x, t) = ∑ cn senh 2
x sen
2
y
n =1
donde cn , satisface
nπ 2 2 nπ
Z
cn senh 3 = f (y) sen y dy
2 2 0 2
Despejando cn , se tiene
Z 2 nπ
1
cn = f (y) sen y dy
senh(3nπ/2) 0 2
Calculemos la integral del lado derecho
Z 2 nπ Z 1 nπ Z 2 nπ
f (y) sen y dy = y sen y dy + (2 − y) sen y dy
0 2 0 2 1 2
4 nπ 2 nπ 4 nπ 2 nπ
= 2 2 sen − cos + 2 2 sen + cos
π n 2 nπ 2 π n 2 πn 2
8 nπ
= 2 2 sen .
π n 2
de donde
8 sen(nπ/2)
cn =
π 2 n2 senh(3nπ/2)
Finalmente, la solución del problema de Dirichlet es
∞
8 sen(nπ/2) nπ nπ
u( x, y) = ∑ 2 2 senh x sen y
n =1
π n senh(3nπ/2) 2 2
Ejemplo 5.16 Determine una solución formal al problema con valores en la frontera
u xx + uyy = 0, 0 < x < π, 0 < y < π
u x (0, y) = u x (π, y) = 0, 0 ≤ y ≤ π
u( x, 0) = cos x − 2 cos(4x ), 0 ≤ x ≤ π
u( x, π ) = 0, 0 ≤ x ≤ π
Solución Suponga que la solución tiene la forma u( x, y) = X ( x )Y (y). Reemplazando en la ecua-
ción diferencial, se tiene
X 00 Y 00
X 00 Y + XY 00 = 0 =⇒ =− = −λ
X Y
182 CAPÍTULO 5. ECUACIONES DIFERENCIALES PARCIALES
X 0 (0) = 0, X 0 (π ) = 0, Y (π ) = 0
X 00 + λX = 0, X 0 (0) = X 0 ( π ) = 0
Este problema tiene como soluciones no triviales cuando λ = 0 y λ = n2 , los cuales son respecti-
vamente
X0 = c, Xn = c cos(nx )
Y 00 − λY = 0, Y (π ) = 0
Si λ = 0, se tiene
Y 00 = 0 Y (y) = c1 y + c2
Y0 (y) = c1 y − c1 π = c1 (y − π )
Por tanto
Yn (y) = −c2 tanh(nπ ) cosh(ny) + c2 senh(ny)
= cn senh(n(y − π ))
u0 ( x, y) = y − π, un ( x, y) = senh(n(y − π )) cos(nx ), n = 1, 2, . . .
forman un conjunto fundamental de solución de la ecuación de Laplace y que satisfacen las condi-
ciones en la frontera
u x (0, y) = u x (π, y) = 0
u( x, π ) = 0
5.5. ECUACIÓN DE LAPLACE 183
Para encontrar una solución que tambien satisfaga la condición u( x, 0) = cos x − 2 cos(4x ), consi-
deramos
∞
u( x, y) = a0 (y − π ) + ∑ an senh(n(y − π )) cos(nx)
n =1
Ejercicios propuestos
1. a) Encuentre la solución u( x, y) de la ecuación de Laplace en el rectángulo 0 < x < a,
0 < y < b, que satisface las condiciones de frontera
u( x, 0) = 0, u( x, b) = g( x ), 0<x<a
u(0, y) = 0, u( a, y) = 0, 0 < y < b
b) Encuentre la solución si
x, 0 < x < a/2
g( x ) = a − x, a/2 ≤ y < a.
uy ( x, 0) = 0, uy ( x, 1) = 0, 0<x<1
u x (0, y) = 0, u x (1, y) = sen(2πy), 0 < y < 1
5. Determine una solución del problema mixto con valores en la frontera para un rectángulo
u xx + uyy = 0, 0 < x < π, 0 < y < 1
u( x, 0) = cos x − cos(3x ), u( x, 1) = cos(2x ), 0<x<π
u x (0, y) = 0, u x (π, y) = 0, 0 < y < 1
Bibliografı́a
[1] B OYCE , W.; D IPRIMA , R.; M EADE , D.Elementary Differential Equations and Boundary Value Pro-
blems, 2017, John Wiley & Sons. Inc.
[2] N AGLE , K.; S AFF , E.; S NIDER , A. Ecuaciones diferenciales y problemas con valores en la frontera.
México, 2005. Cuarta edición, Pearson Educación.
[3] Z ILL , D. Ecuaciones diferenciales con valores en la frontera. México, 2018. Novena edición. Cengage
Learning Editores, S.A. de C.V.
185