Ecuaciones Diferenciales Ordinarias: Guía Completa
Ecuaciones Diferenciales Ordinarias: Guía Completa
Índice
1. Ecuaciones diferenciales de primer orden y problemas de valor
inicial. 2
1.1. Introducción a las ecuaciones diferenciales . . . . . . . . . . . . . 2
1.1.1. De…nición y clasi…cación. . . . . . . . . . . . . . . . . . . 2
1.2. Problemas de valor inicial. . . . . . . . . . . . . . . . . . . . . . . 4
1.3. Ecuaciones separables y de coe…cientes homogéneos. . . . . . . . 11
3. Transformadas de Laplace 58
3.1. De…nición, propiedades de linealidad y formulas básicas de trans-
formada de Laplace. . . . . . . . . . . . . . . . . . . . . . . . . . 58
3.2. De…nición, propiedades de linealidad y formulas básicas de trans-
formada inversa de Laplace. . . . . . . . . . . . . . . . . . . . . . 63
3.3. Teoremas y propiedades. . . . . . . . . . . . . . . . . . . . . . . . 65
1
3.3.1.
Primer teorema de traslación. . . . . . . . . . . . . . . . . 65
3.3.2.
Transformada de Laplace de la función escalón. . . . . . . 66
3.3.3.
Segundo teorema de traslación. . . . . . . . . . . . . . . . 67
3.3.4.
Teorema de convolución. (Investigación por parte del alum-
no) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
3.4. Solución de ecuaciones diferenciales usando transformada de Laplace. 68
d2 u d2 u @2u @2u
5. dx2 = dt2 2 du
dt , @x2 = @t2 2 @u
@t , uxx = utt 2ut
du dv @u @v
6. dy = dx , @y = @x
7. (x + y) dx + (x y) dy = 0 (Forma diferencial)
(x + y) dx + (x y) dy = 0
dy
(x + y) + (x y) = 0
dx
dy
(x y) +y = x
dx
2
son ejemplos de ecuaciones diferenciales.
Las ecuaciones diferenciales se clasi…can según el tipo, el orden, el grado y
según la linealidad.
y 00 = x cos y
5
c) ln x + y 00 = y 000 + 7 (y 0 ) , es de tercer orden.
iii) Clasi…cación por grado
El grado de una ED se re…ere al grado o exponente de la mayor derivada.
dn y dn 1 y dy
an (x) + an 1 (x) + + a1 (x) + a0 (x) y = f (x)
dxn dxn 1 dx
3
De acuerdo a la forma de la ecuación anterior, podemos decir que una
EDO es lineal si se cumplen las siguientes condiciones:
a) La variable dependiente y todas sus derivadas son de primer
grado, es decir, la potencia de cada variable dependiente o sus
derivadas es 1:
b) Los coe…cientes a0 , a1 , . . . , an dependen solo de la variable inde-
pendiente y también el miembro del lado derecho.
Ejemplo:
a) y 00 2y 0 + y = 0, es lineal
b) (1 y) y 0 + 2y = ex , no es lineal.
2
d y dy
c) 1 x2 dx 2 +x dx xy = 3x2 ex , EDO, segundo orden, n = 2, a2 (x) = 1 x2 ,
a1 (x) = x, a0 (x) = x, f (x) = 3x2 ex
2
d y dy
a2 (x) dx2 + a1 (x) dx + a0 (x) y = f (x), sí es lineal.
y (x0 ) = y0
y 0 (x0 ) = y1
..
.
yn 1
(x0 ) = yn 1
4
Comprobar si
y = cos x
es solución de la ED.
Solución. Derivando
y0 = sin x
y 00 = cos x
entonces
y 00 + y = cos x + cos x
= 0
Ejemplo. Sea
y 00 + y = 0
Comprobar si
y = sin x
es solución de la ED.
Solución. Derivando
y0 = cos x
y 00 = sin x
y 00 + y = sin x + sin x
= 0
sí es una solución.
La solución general es
y = c1 cos x + c2 sin x
5
Solución. Usando derivación implícita
d h 2
i d
(x 5) + y 2 = (1)
dx dx
d h 2 di
(x 5) y2+ = 0
dx dx
dy
2 (x 5) + 2y = 0
dx
dy
(x 5) + y = 0
dx
dy x 5
=
dx y
2 2
dy (x 5)
=
dx y2
2 2
dy (x 5)
+1 = +1
dx y2
2 2
dy (x 5) y2
+1 = +
dx y2 y2
2 2
dy (x 5) + y 2
+1 =
dx y2
2
dy 1
+1 =
dx y2
2
Otra forma es usar el teorema de la función implícita, con F (x; y) = (x 5) +
y2 1 = 0
@F
dy @x
= @F
dx @y
2 (x 5)
=
2y
x 5
=
y
Ejemplo. Considerese el PVI siguiente, compruebe que y = ex + C es la
solución general y resuelva el PVI
dy
= ex
dx
y (0) = 4
6
Solución.
y = ex + C
dy
= ex
dx
y (0) = e0 + C = 4
1+C = 4
C = 5
y = x2 + C
dy
= 2x
dx
dy
x = 2y
dx
x (2x) = 2 x2 + C
2x2 6= 2x2 + 2C
y = Cx2
dy
= 2Cx
dx
Ahora sustituyendo
dy
x = 2y
dx
x (2Cx) = 2 Cx2
2Cx2 = 2Cx2
y = Cx2
7
Ahora, consideremos la CI
y (0) = 8
2
y (0) = C (0) = 8
0 =
6 8
y (4) = 8
2
y (4) = C (4)
= 8
16C = 8
1
C =
2
La solución particular al PVI es
1 2
y= x
2
Ejemplo. Sea
y 00 + 5y 0 + 6y = 0
y (0) = 1
0
y (0) = 6
y = c1 e 2x + c2 e 3x
y0 = 2c1 e 2x 3c2 e 3x
y 00 = 4c1 e 2x + 9c2 e 3x
Sustituyendo en la ecuación
y 00 + 5y 0 + 6y = 0
8
4c1 e 2x + 9c2 e 3x
+5 2c1 e 2x 3c2 e 3x
2x 3x
+6 c1 e + c2 e
= 4c1 e 2x + 9c2 e 3x
10c1 e 2x 15c2 e 3x
+6c1 e 2x + 6c2 e 3x
= 4c1 e 2x 10c1 e 2x + 6c1 e 2x
= 0
y = c1 e 2x + c2 e 3x
y (0) = c1 + c2 = 1
y0 = 2c1 e 2x 3c2 e 3x
c1 + c2 = 1
2c1 3c2 = 6
y 00 2y 0 + y = 0
y (0) = 3
y 0 (0) = 5
cuya solución es
y = c1 ex + c2 xex
Comprobar que la solución dada, realmente es la solución y encontrar los valores
de c1 y c2 .
Solución.
y = c1 ex + c2 xex
y0 = c1 ex + c2 (xex + ex )
= c1 ex + c2 (x + 1) ex
y 00 = c1 ex + c2 (x + 1) ex + c2 ex
= c1 ex + c2 (x + 2) ex
9
Sustituyendo en la ED
y 00 2y 0 + y = 0
c1 ex + c2 (x + 2) ex
2 [c1 ex + c2 (x + 1) ex ]
+c1 ex + c2 xex
= c1 ex + c2 xex + 2c2 ex
2c1 ex 2c2 xex 2c2 ex
+c1 ex + c2 xex
= 2c1 ex 2c1 ex
+2c2 xex 2c2 xex
+2c2 ex 2c2 ex
= 0
Ahora, tenemos que
y = c1 ex + c2 xex
y (0) = c1 = 3
y 0 = c1 ex + c2 (x + 1) ex
y 0 (0) = c1 + c2 = 5
3 + c2 = 5
c2 = 8
Por lo tanto la solución al PVI es
y = 3ex 8xex
Ejemplo. Compruebe que la expresión dada es una solución implícita de la
ED
(1 xy) y 0 = y2
y = exy
Solución. Usando el teorema de la función implícita
@F
y0 = @x
@F
@y
10
@F
y0 = @x
@F
@y
yexy
=
1 xexy
yexy
=
1 xexy
(1 xexy ) y 0 = yexy
(1 xy) y 0 = y2
dy
De…nición. Si la ecuación diferencial dx = f (x; y) puede ser escrita como
dy
= g (x) h (y)
dx
decimos que la ecuación es de variables separables, donde g (x) es una función
que depende solo de x y h (y) es una función que depende solo de y.
En caso de estar escrita en la forma diferencial
M (x; y) dx + N (x; y) dy = 0
entonces es necesario que tanto M (x; y) = g1 (x) h1 (y) y N (x; y) = g2 (x) h2 (y).
El método de solución consta de escribir la ecuación diferencial de tal forma
que queden las variables separadas para poder integrar, sin olvidar la constante
de integración. La solución puede aparecer de forma explícita o implícita.
Ejemplo. Resolver la ED
dy
ex+y =x
dx
Solución. Reescribiendo la ED
dy
ex ey = x
dx
dy x y
= xe e
dx
ey dy = xe x dx
Z Z
ey dy = xe x dx
ey = e x (x + 1) + C
y = ln e x (x + 1) + C
ey + e x
(x + 1) = C
11
Comprobación. Usando el teorema de la función implícita (TFI)
@F
dy @x
= @F
dx @y
donde
F (x; y) = ey + e x
(x + 1) C=0
@F x x
= e e (x + 1)
@x
= e x xe x
e x
x
= xe
@F
= ey
@y
@F
dy @x
= @F
dx @y
x
dy xe
=
dx ey
dy x y
= xe e
dx
dy (x+y)
= xe
dx
dy
ex+y =x
dx
Ejemplo. Resolver el PVI
xyy 0 = 1 + y 2 ; y (1) = 3
12
Solución. Reescribiendo
dy
xy = 1 + y2
dx
dy 1 1 + y2
=
dx x y
ydy dx
=
1 + y2 x
Z Z
ydy dx
=
1 + y2 x
Z Z
1 2ydy dx
=
2 1 + y2 x
1 1
ln 1 + y 2 = ln x + ln C
2 2
ln 1 + y 2 = 2 ln x + ln C
2
ln 1 + y = ln x2 + ln C
ln 1 + y 2 ln x2 = ln C
2
1+y
ln = ln C
x2
1 + y2
= C
x2
x 2 1 + y2 = C
Comprobación.
@F
dy @x
= @F
dx @y
2(1+y 2 )
x3
= 2y
x2
2 2
2x 1+y
=
2x3 y
dy 1 + y2
=
dx xy
dy
xy = 1 + y2
dx
y (1) = 3
1 + 32
= C
12
C = 10
13
Por lo tanto, la solución particular es
1 + y2
= 10
x2
Ejemplo: Resolver la ED
sin x cos2 ydx cos x sin ydy = 0
Solución.
sin x cos2 ydx cos x sin ydy = 0
sin x cos2 ydx = cos x sin ydy
sin x sin y
dx= dy
cos x cos2 y
Z Z
sin x sin y
dx= dy
cos x cos2 y
Z Z
sin x 1 sin y
dx= dy
cos x cos y cos y
Z
ln (cos x) = sec y tan ydy
ln (cos x) = sec y + C1
ln (cos x) + sec y = C
Comprobación. Usando el TFI
@F
dy @x
= @F
dx @y
sin x
dy cos x
=
dx sec y tan y
sin x
dy cos x
= 1 sin y
dx cos y cos y
sin x
dy cos x
= sin y
dx cos2 y
2
dy sin x cos y
=
dx cos x sin y
cos x sin ydy = sin x cos2 ydx
sin x cos2 ydx cos x sin ydy = 0
Ecuaciones diferenciales de coe…cientes homogéneos
De…nición. (Función homogénea) Sea z = f (x; y) una función de dos
variables, diremos que f es homogénea de grado n, si para 2 R, tenemos que
n
f ( x; y) = f (x; y)
14
Ejemplos. Determine si las funciones siguientes son homogéneas o no, de
serlo, determine el grado.
1. f (x; y) = x3 + 2x2 y + 3xy 2 + y 3
p
2. f (x; y) = x2 + y 2
x
3. f (x; y) = ln y
p
4. f (x; y) = x2 + xy
Solución.
1.
3 2 2 3
f ( x; y) = ( x) + 2 ( x) ( y) + 3 ( x) ( y) + ( y)
3 3
= x + 2 2 x2 y + 3 x 2 y 2 + 3 y 3
3 3
= x + 3 2x2 y + 3 3xy 2 + 3 y 3
3
= x3 + 2x2 y + 3xy 2 + y 3
3
= f (x; y)
Sí es homogénea de grado 3.
2.
q
2 2
f ( x; y) = ( x) + ( y)
q
2 2
= x + 2 y2
q
2
= (x2 + y 2 )
p p
2
= x2 + y 2
p
= x2 + y 2
1
= f (x; y)
Sí es homogénea de grado 1
3.
x
f ( x; y) = ln , 6= 0
y
x
= ln
y
0 x
= ln
y
0
= f (x; y)
Sí es homogénea de grado 0.
15
4.
2
p
f ( x; y) = ( x) + ( x) ( y)
q
2 2 2
= x + xy
p p
2 2 2
= x + xy
2 2 p
= x + xy
p
= x2 + xy
n
6 = f (x; y)
No es homogénea.
p
4. f (x; y) = x + y
Solución.
p
f ( x; y) = x+ y
p
= (x + y)
p p
= x+y
p
= f (x; y)
1
= 2 f (x; y)
Sí es homogénea de grado 12 .
De…nición. Sea
M (x; y) dx + N (x; y) dy = 0
diremos que la ecuación diferencial es una ecuación de coe…cientes homogé-
neos si tanto M y N son funciones homogéneas del mismo grado.
El método de solución para una ecuación diferencial de coe…cientes homogé-
neos, es hacer el cambio de variable x = vy o y = ux y por lo tanto
dx = vdy + ydv
dy = udx + xdu
x2 + y 2 dx xydy = 0
16
dy
x2 + y 2 xy = 0
dx
dy x2 + y 2
=
dx xy
x2 +y 2
dy x2
= xy
dx x2
2
dy 1 + xy
= y
dx x
2
du 1+u
u+x =
dx u
du 1 + u2
x = u
dx u
du 1 + u2 u2
x =
dx u
du 1
x =
dx u
dx
udu =
x
dy du
=u+x
dx dx
Solución. Sea y = ux, entonces dy = udx + xdu, al sustituir en la ecuación
diferencial, tenemos que
x2 + y 2 dx xydy = 0
h i
2
x2 + (ux) dx x (ux) (udx + xdu) = 0
h i
2
x2 + (ux) x2 u2 dx x3 udu = 0
x2 + u2 x2 u2 x2 dx x3 udu = 0
2 3
x dx x udu = 0
dx xudu = 0
dx = xudu
dx
= udu
x
17
Z Z
dx
= udu
x
1 2
ln x = u +C
2
1 y 2
ln x = +C
2 x
y 2
2 ln x = +C
x
y 2
ln x2 = +C
x
y 2
x2 = e( x ) +C
y 2
x2 = e( x ) eC
y2
x2 = Ce x2
y2
x2 e x2 =C
Comprobación. Usando el TFI
@F
dy @x
= @F
dx @y
2 y2 y2
x2 2 xy 3 e x2 + 2xe x2
= y2
2y
x2 x2 e x2
y2
y2
2 x +x e x2
= y2
2ye x2
y2
x +x
=
y
x2 +y 2
x
=
y
x + y2
2
=
xy
dy x2 + y 2
=
dx xy
xydy = x + y 2 dx
2
x2 + y 2 dx xydy = 0
Ejemplo: Resuelva la ED
2xydx x2 y 2 dy = 0
18
Solución. Sea x = vy, dx = vdy + ydv, sustituyendo en la ED
h i
2
2 (vy) y [vdy + ydv] (vy) y 2 dy = 0
2vy 3 dv v2 y2 y2 2v 2 y 2 dy = 0
3 2 2 2
2vy dv y v y dy = 0
3 2 2 2
2vy dv + y + v y dy = 0
3 2 2
2vy dv + y 1+v dy = 0
2vy 3 dv = y 2 1 + v 2 dy
2v 1
dv = dy
1 + v2 y
Z Z
2vdv dy
=
1 + v2 y
ln 1 + v 2 = ln y + ln C
2
ln 1 + v + ln y = ln C
2
ln y 1 + v = ln C
" #
2
x
ln y 1 + = ln C
y
" #
2
x
y 1+ = C
y
x2
y 1+ = C
y2
x2 + y 2
y = C
y2
x2 + y 2
= C
y
19
x2 +y 2
Comprobación. Usando el TFI con F (x; y) = y C = 0, tenemos que
@F
dy @x
= @F
dx @y
2x
y
= 2y 2 (x2 +y 2 )
y2
2x
y
= 2y 2 x2 y 2
y2
2x
y
= x2 y2
y2
2xy 2
=
y (x2 y 2 )
2xy
=
x2 y 2
dy 2xy
= 2
dx x y2
x2 y 2 dy = 2xydx
2xydx x2 y 2 dy = 0
20
Teorema. Una ecuación diferencial es exacta si
@M @N
=
@y @x
@2f @ @f
=
@x@y @x @y
@
= x2 + 4xy
@x
= 2x + 4y
@2f @ @f
=
@y@x @y @x
@
= 2xy + 2y 2
@y
= 2x + 4y
2x2 y + xy dx + xy 2 xy dy = 0
Solución.
My = 2x2 + x
Nx = y2 y
Como
My 6= Nx
entonces la ecuación diferencial no es exacta.
Ejemplo. Sea
8x 2y 2 + 2 dx 4xydy = 0
Determine si la ecuación diferencial es exacta.
Solución.
My = 4y
Nx = 4y
21
Solución.
My = x cos y
Nx = y sin x
Como My 6= Nx , la ecuación diferencial no es exacta.
Ejemplo. Sea ey dx + xey dy = 0, determine si la ecuación diferencial es
exacta o no.
Solución.
ey dx + xey dy = 0
dx + xdy = 0
My = ey
Nx = ey
Como My = Nx , la ecuación diferencial sí es exacta.
Ejemplo. Demuestre que la ecuación 2xydx + 1 + x2 dy = 0 es exacta y
determine la
solución general.
Solución.
My = 2x = Nx
Sí es exacta. Ahora como la ecuación diferencial es exacta, existe una función
f , tal que la
ecuación es la diferencial df = @f @f
@x dx + @y dy , es decir que
@f
= M (x; y)
@x
@f
= 2xy
@x
df = 2xydx
Z x
f (x; y) = 2xydx
= x2 y + h (y)
Ahora, tenemos que
@f
= N (x; y)
@y
= 1 + x2
pero
@f @
= x2 y + h (y)
@y @y
= x2 + h0 (y)
22
Entonces
x2 + h0 (y) = 1 + x2
h0 (y) = 1
h (y) = y
Por lo tanto
f (x; y) = x2 y + y = C
La solución general es
x2 y + y = C
Comprobación. Una forma es usar el TFI
@F
dy @x
= @F
dx @y
2xydx + x2 + 1 dy = 0
My = 6x 4y
Nx = 6x 4y
Entonces
23
de aquí que
f (x; y) = 3x2 y 2xy 2 + c1
Por lo tanto la solución general es
3x2 y 2xy 2 = C
My = 6x 4y
Nx = 6x 4y
Entonces
Z x Z y
f (x; y) = 6xy 2y 2 dx + 3x2 4xy dy = C
= 3x2 y 2xy 2 = C
3x2 y 2xy 2 = C
Ejemplo. Sea
8x 2y 2 + 2 dx 4xydy = 0
Determine si la ecuación diferencial es exacta y su solución.
Solución.
My = 4y
Nx = 4y
2xy + 3x2 dx + 3y 2 + x2 dy = 0
Solución.
My = 2x
Nx = 2x
24
como My = Nx , entonces la ecuación sí es exacta.
Para solucionarla, sea f (x; y) una función tal que
@f
= M
@x
@f
= N
@y
es decir
@f
= 2xy + 3x2
@x
f (x; y) = x2 y + x3 + h (y)
pero
@f
= x2 + h0 (y)
@y
= 3y 2 + x2
entonces
h0 (y) = 3y 2
h (y) = y 3
Por lo que
f (x; y) = x2 y + x3 + y 3
Por lo tanto la solución de la edo es
x2 y + x3 + y 3 = C
Z x Z y
2xy + 3x2 dx + 3y 2 + x2 dy = 0
x2 y + x3 + y 3 = C
Factores de integración.
De…nición. Sea
M (x; y) dx + N (x; y) dy = 0
una función (x; y) se dice que es un factor integrante de la ED si al multi-
plicar la ecuación por la función , esta se hace exacta.
Ejemplo. (x; y) = ey es un factor integrante de
dx + xdy = 0
@M
= 0
@y
@N
= 1
@x
25
ey dx + xey dy = 0
@M
= ey
@y
@N
= ey
@x
Consideremos el cómo obtener un factor integrante, la ecuación diferencial
no exacta
M dx + N dy = 0
sería exacta si al multiplicar por (x; y), obtenemos que
M dx + N dy = 0
@ @
( M) = ( N)
@y @x
@M @ @N @
+M = +N
@y @y @x @x
@M @N @ @
= N M
@y @x @x @y
Para simpli…car, supongamos que := (x), entonces
@M @N @
= N
@y @x @x
@M @N d
= N
@y @x dx
@M @N
d @y @x
= dx
N
( @M
@y
@N
@x )
Ahora supongamos que existe g (x) = N , entonces
Z Z @M @N
d @y @x
= dx
N
Z
ln = g (x) dx
R
g(x)dx
(x) = e
en donde
@M @N
@y @x
g (x) =
N
Analogamente, tenemos que
R
h(y)dy
(y) = e
26
en donde
@N @M
@x @y
h (y) :=
M
Ejemplo. Determine el factor integrante de
Solución.
My = 3x2
Nx = 0
My Nx 3x2 0
=
N y
2
3x
= 6= g (x)
y
Pero
Nx My 0 3x2
=
M 3x2 y
3x2
=
3x2 y
1
= = h (y)
y
Por lo tanto el factor integrante está dado por
R
h(y)dy
(y) = e
R dy
= e y
ln y
= e
1
= eln y
= y 1
1
=
y
Veri…quemos que esto es así,
1 1
3x2 ydx + ydy = (0)
y y
3x2 dx + dy = 0
Como
My = 0 = N x
27
la ecuación sí es exacta. Su solución es
x3 + y = C
y= x3 + C
Comprobación.
Solución.
My = 0
Nx = y sin x
My = 0 = N x
28
Z x Z y
x sec x tan xdx ydy = C
y2
x sec x ln (sec x + tan x) = C
2
Ahora usando la C.I. x = 0 y y = 2
22
0 sec 0 ln (sec 0 + tan 0) = C
2
C = 2
2
y
x sec x ln (sec x + tan x) = 2
2
Z Z
x sec x tan xdx = x sec x sec xdx
u = x
du = dx
dv = sec x tan xdx
v = sec x
3xy + y 2 + x2 + xy y 0 = 0
Solución.
3xy + y 2 + x2 + xy y 0 = 0
dy
3xy + y 2 + x2 + xy = 0
dx
3xy + y 2 dx + x2 + xy dy = 0
My = 3x + 2y
Nx = 2x + y
29
Entonces, existe un factor integrante
R
g(x)dx
(x) = e
R dx
= e x
= eln x
= x
Multiplicando por el factor integrante
x 3xy + y 2 dx + x2 + xy dy = 0
2 2 3 2
3x y + xy dx + x + x y dy = 0
Ahora veri…quemos que la ecuación diferencial ya es exacta
My = 3x2 + 2xy
Nx = 3x2 + 2xy
La solución es
1
x3 y + x2 y 2 = C1
2
2x3 y + x2 y 2 = C
Comprobación.
6x2 y + 2xy 2 dx + 2x3 + 2x2 y dy = 0
2 2 3 2
3x y + xy dx + x + x y dy = 0
2 2
3xy + y dx + x + xy dy = 0
dy
1+x = 0
dx
dy 1
+ = 0
dx x
dy 1
=
dx x
30
donde P (x) = 0 y Q (x) = x1 .
El método de solución consta de encontrar un factor integrante := (x),
de tal forma que
la ecuación lineal se vuelva exacta. Consideremos que
d dy d
( y) = +y
dx dx dx
y comparemosla con
dy
+P y = Q
dx
entonces
d
= P
dx
separando las variables
d
= P dx
integrando
Z
ln = P dx
R
P dx
(x) = e
Integrando Z
R R
P dx P dx
e y= e Qdx + C
R
Z R
P dx P dx
y=e e Qdx + C
Ejemplo. Resolver
dy + 3x2 y x2 dx = 0
31
la ecuación diferencial es lineal, con P (x) = 3x2 y Q (x) = x2 . Ahora susti-
tuimos en la fórmula
R
Z R
P dx
y=e e P dx Qdx + C
R
Z R
3x2 dx 3x2 dx 2
y = e e x dx + C
Z
x3 3
= e ex x2 dx + C
Z
x3 1 3
= e ex 3x2 dx + C
3
3 1 x3
= e x e +C
3
3 1
= Ce x +
3
3 1
y = Ce x +
3
Ejemplo. Determine la solución de la ecuación
1
y0 =
x sin y + 2 sin 2y
Solución. Reescribiendo la ecuación
dy 1
=
dx x sin y + 2 sin 2y
dx x sin y + 2 sin 2y
=
dy 1
dx
= x sin y + 2 sin 2y
dy
dx
+ ( sin y) x = 2 sin 2y
dy
La ecuación diferencial es lineal, pero en y, entonces P (y) = sin y y Q (y) =
2 sin 2y. Por lo tanto
R
Z R
x = e P dy e P dy Qdy + C
R
Z R
x = e ( sin y)dy e ( sin y)dy 2 sin 2ydy + C
R
Z R
x = e sin ydy e sin ydy 2 sin 2ydy + C
Z
x = e cos y ecos y 2 sin 2ydy + C
cos y
x = e [ 4ecos y (cos y 1) + C]
32
Z
ecos y 2 sin 2ydy
Recordemos que
Sea
u = cos y
du = sin ydy
Entonces Z Z
cos y
4 e cos y ( sin ydy) = 4 eu udu
= 4ecos y (cos y 1)
Solución. Reescribiendo
2
xy 0 = (x + y)
dy 2
x = (x + y)
dx
dy
x = x2 + 2xy + y 2
dx
dy 1 2
= x + 2y + x y
dx
dy
x = x2 + 2xy + y 2
dx
dy x2 + 2xy + y 2
=
dx x
dx x
=
dy x2 + 2xy + y 2
33
La ecuación no es lineal.
dy 1 2
= x + 2y + x y
dx
dy 1 2
2y = x+ y
dx x
dy x2 + y 2
2y =
dx x
Ecuación de Bernoulli
De…nición. Una ecuación diferencial de la forma
dy
+ P (x) y = Q (x) y n , n 6= 0; 1
dx
se le llama la ecuación de Bernoulli.
Su solución se obtiene al proponer un cambio de variable u = y 1 n
, entonces
du n dy
= (1 n) y
dx dx
n
Multiplicando la ecuación de Bernoulli por (1 n) y , obtenemos que
n dy n n
(1 n) y + (1 n) y P (x) y = (1 n) y Q (x) y n
dx
du
+ (1 n) P (x) u = (1 n) Q (x)
dx
Esta ecuación es lineal en u
du
+ P1 (x) u = Q1 (x)
dx
con P1 (x) = (1 n) P (x) y Q1 (x) = (1 n) Q (x). Su solución está dada
por Z
R R
P1 dx P1 dx
u=e e Q1 dx + C
34
vemos que es una ecuación de Bernoulli, en donde P (x) = x1 , Q (x) = x y
n = 2. Por lo tanto, al sustituir en la fórmula, tenemos que
R dx
Z R dx
y 1 2 = e(2 1) x (1 2) e(1 2) x xdx + C
R dx
Z R dx
y 1 = e x e x xdx + C
Z
y 1 = eln x e ln x xdx + C
Z
1
y 1 = x eln x xdx + C
Z
1
y = x x 1 xdx + C
Z
1
y = x dx + C
1
y = x ( x + C)
1
y =
x (C x)
Ejemplo. Resuelve la ecuación
3
y0 y = y5
x
Solución. Usando la fórmula
R
Z R
y1 n
= e(n 1) P dx
(1 n) e(1 n) P dx
Qdx + C
Z
4 12 ln x
y = e 4 e12 ln x dx + C
Z
4 12
y = x 4 x12 dx + C
4 12 1 13
y = x 4 x +C
13
4 4 12
y = x + Cx
13
4 13
x12 y 4
= x +C
13
4 13
x12 y 4
x + = C
13
x12 4
4
+ x13 = C
y 13
13x12 + 4x13 y 4
= C
13y 4
35
Comprobación. Tarea
Tarea 5.
y dy y
22. x e x 1 dx = e x (y x); y (1) = 0
Solución. Haciendo u = xy , y = ux
dy du
=u+x
dx dx
y sustituyendo
y dy y
x ex 1 = e x (y x)
dx
du
x (eu 1) u + x = eu (ux x)
dx
du
x (eu 1) u + x = xeu (u 1)
dx
du
(eu 1) u + x = eu (u 1)
dx
du eu (u 1)
u+x =
dx eu 1
du eu (u 1)
x = u
dx eu 1
du eu (u 1) u (eu 1)
x =
dx eu 1
du ueu eu ueu + u
x =
dx eu 1
u
du u e
x =
dx eu 1
dx eu 1
= du
x u eu
Z Z
dx eu 1
= du
x u eu
Z u
e 1
ln x = du
eu u
ln x = ln (eu u) + ln C
ln x + ln (eu u) = ln C
ln [x (eu u)] = ln C
x (eu u) = C
36
y y
x ex = C
x
y
xe x y = C
y
:K V !V
( ; ~v ) 7! ~v
cumple que
(vi) es cerrada ~v 2 V .
(vii) ( ~v ) = ( ) ~v :
(viii) ( + ) ~v = ~v + ~v :
(ix) (~u + ~v ) = ~u + ~v :
37
(x) 1 ~v = ~v
K es un campo. (K = R)
Llamamos a V un espacio vectorial real.
Ejemplos. Rn es un espacio vectorial.
Mn m (R) es un espacio vectorial.
Pn [x] es un espacio vectorial
C [a; b] es un espacio vectorial
C 1 [a; b] es un espacio vectorial
Subespacio vectorial
De…nición. Sea S V , diremos que S es un subespacio de V (S V ) si
S 6= ? y (S; +; ) es un espacio vectorial.
Teorema. S V ssi
(i) ~u + ~v 2 S, 8~u; ~v 2 S:
(ii) ~u 2 S, 8 2 R y ~u 2 S.
Combinación lineal
De…nición. Sea ~u 2 V , diremos que ~u es una combinación lineal del conjunto
S = f~v1 ; ~v2 ; : : : ; ~vn g si se puede expresar de la forma
n
X
~u = i~
vi
i=1
38
Base de un espacio vectorial
De…nición. Una base B es un subconjunto del espacio vectorial tal que
gen (B) = V
y B es linealmente independiente.
Dimensión de un espacio vectorial
De…nición. Sea V un espacio vectorial
dim V := #B
en donde B es una base para V .
De…nición. Sea V = C n (R), recordemos que C n (R) es un espacio vectorial
C n (R) := ff jf : R ! Rg
entonces ci = 0 8i.
Ejemplo. Sea f1 (x) ; f2 (x) 2 C n (R), con f1 = 3x + 4 y f2 = 4x + 3, decir
si el conjunto S = ff1 ; f2 g es o no linealmente independiente.
Solución. Sean c1 ,c2 2 R
c1 f1 + c2 f2 = 0
es decir
c1 (3x + 4) + c2 (4x + 3) = 0
(3c1 + 4c2 ) x + (4c1 + 3c2 ) = 0
3c1 + 4c2 = 0
4c1 + 3c2 = 0
3 4
A =
4 3
3 4
det A = det
4 3
= (3) (3) (4) (4)
= 9 16
= 7
39
la solución es c1 = 0; c2 = 0: Entonces S es linealmente independiente.
Ejemplo. Sean f1 = 3 + 2x 5x2 , f2 = 3 5x y f3 = 5x 2x2 , decir si el
conjunto S = ff1 ; f2 ; f3 g es o no linealmente independiente.
Solución.
c1 f1 + c2 f2 + c3 f3 = 0
2
c1 3 + 2x 5x + c2 (3 5x) + c3 5x 2x2 = 0
(3c1 + 3c2 + 0c3 ) + (2c1 5c2 + 5c3 ) x + ( 5c1 + 0c2 2c3 ) x2 = 0
2.1.2. Wronskiano.
De…nición. Sea S = ff1 ; f2 ; f3 ; : : : ; fn g, de…nimos el determinante
f1 f2 fn
f10 f20 fn0
W (f1 ; f2 ; : : : ; fn ) := .. .. .. ..
. . . .
(n 1) (n 1) (n 1)
f1 f2 fn
al cual le llamaremos el Wronskiano de las funciones fi .
Teorema. S = ff1 ; f2 ; f3 ; : : : ; fn g es linealmente independiente en un in-
tervalo (a; b), si y sólo si W (f1 ; : : : ; fn ) 6= 0, para
x 2 (a; b).
Ejemplo. Demuestre que el conjunto S = ff1 ; f2 ; f3 g con f1 = 3 + 2x 5x2 ,
f2 = 3 5x y f3 = 5x 2x2 , es linealmente
independiente en todo R.
Solución. Calculemos el Wronskiano de las funciones
3 + 2x 5x2 3 5x 5x 2x2
W (f1 ; f2 ; f3 ) = 2 10x 5 5 4x
10 0 4
3 5x 5x 2x2 3 + 2x 5x2 3 5x
= 10 4
5 5 4x 2 10x 5
= 10 (3 5x) (5 4x) ( 5) 5x 2x2 4 3 + 2x 5x2 ( 5) (2 10x) (3 5x)
2 2 2
= 10 15 12x 25x + 20x + 25x 10x 4 15 10x + 25x 6 + 10x + 30x 50x2
= 10 15 12x + 10x2 4 21 + 30x 25x2
= 150 + 120x 100x2 + 84 120x + 100x2
= 66 6= 0
40
Como W (f1 ; f2 ; f3 ) 6= 0, concluimos que S es linealmente independiente
para todo R.
Ejemplo. Demuestre que f1 (x) = ex y f2 (x) = e x son linealmente inde-
pendientes.
Solución. Calculemos el Wronskiano de las funciones f1 y f2
ex e x
W ex ; e x
=
ex e x
= (ex ) e x
(ex ) e x
= 1 1
= 2= 6 0
41
Para solucionar este tipo de ecuaciones, asociamos a la ecuación diferencial
una ecuación algebraica, denominada la ecuación característica y depen-
diendo de las soluciones de dicha ecuación es como será la solución de la
ecuación diferencial.
x3 x2 4x + 4 = 0
(x 1) (x + 2) (x 2) = 0
x1 = 1, x2 = 2, x3 = 2
1 1 1 4 4
1 2 2
1 2 2 6
1 1 1 4 4
1 0 4
2 1 0 4 0
2 4
2 1 2 0
2
1 0
Tarea de investigación. ¿Cómo resolver o encontrar los ceros de
un polinomio o encontrar las raíces de una ecuación?
Teorema del factor
Teorema del residuo
Ceros posibles racionales
Método de Ru¢ ni o Horner
División sintética
Teorema fundamental del álgebra.
Regla de los signos de Descártes.
Criterio de acotación de las raíces.
42
Encontrar su ecuación característica.
Solución.
9 2 + 18 + 10 = 0
Ejemplo. Sea
y (5) + 4y (4) + 5y 000 6y 0 4y = 0
Encontrar su ecuación característica.
Solución.
5
+4 4+5 3 6 4=0
1x 2x
yh = c1 e + c2 e
1x 1x
yh = c1 e + c2 xe
x x
yh = c1 e cos x + c2 e sin x
en donde
1;2 = i 2C
Ejemplo. Encuentre la solución para la ecuación
y 00 2y 0 3y = 0
de donde
( + 1) ( 3) = 0
1 = 1, 2 =3
Por lo tanto la solución es
x
yh = c1 e + c2 e3x
y 00 + 2y 0 + y = 0
43
Solución. La ecuación característica asociada a la EDO es
2
+2 +1 = 0
2
( + 1) = 0
9y 00 + 18y 0 + 10y = 0
x 1 x 1
yh = c1 e cos x + c2 e sin x
3 3
Ejemplo. Resolver la ecuación diferencial
Solución.
44
Considerando la ecuación auxiliar
5
+4 4+5 3 6 4 = 0
2
( 1) ( + 2) ( + 1) +2 +2 = 0
, Solution is: 2; 1; 1; 1 + i; 1 i
2x x
yh = c1 e + c2 e + c3 ex + c4 e x
cos x + c5 e x
sin x
y 00 + 2y 0 + 4y = 5x4 + 3x2 x
x
p x
p
yh = c1 e cos 3x + c2 e sin 3x
Entonces proponemos
Derivándola
yp0 = 4Ax3 + 3Bx2 + 2Cx + D
45
y
yp00 = 12Ax2 + 6Bx + 2C
Sustituyendo en la ecuación
y 00 + 2y 0 + 4y =
12Ax2 + 6Bx + 2C +
2 4Ax3 + 3Bx2 + 2Cx + D
+4 Ax4 + Bx3 + Cx2 + Dx + E
= 12Ax2 + 6Bx + 2C
8Ax3 + 6Bx2 + 4Cx + 2D
+4Ax4 + 4Bx3 + 4Cx2 + 4Dx + 4E
4A + 0B + 0C + 0D + 0E = 5
8A + 4B + 0C + 0D + 0E = 0
12A + 6B + 4C + 0D + 0E = 3
0A + 6B + 4C + 4D + 0E = 1
0A + 0B + 2C + 2D + 4E = 0
cuya solución es A = 45 ; B = 25 ; C = 34 ; D = 11
4 ;E =
7
4.
Por lo que
5 5 3 3 2 11 7
yp = x4 x + x + x
4 2 4 4 4
Por lo tanto la solución general es
y = yh + yp
p p
= c1 e x cos 3x + c2 e x sin 3x +
5 5 3 3 2 11 7
+ x4 x + x + x
4 2 4 4 4
Ejemplo. Encontrar la solución de la ecuación
y 00 y=8
yp = A
46
de donde
yp0 = 0
yp00 = 0
sustituyendo en la ED
0 A = 8
A = 8
A = 8
de aquí que
yp = 8
Por lo tanto la solución general es
x
y = c1 e + c2 ex 8
Ejemplo. Resolver
y 00 + 4y = 2e x
Solución. Como
yh = c1 cos 2x + c2 sin 2x
Proponemos que
x
yp = Ae
de donde
yp0 = Ae x
yp00 = Ae x
al sustituir
x
Ae + 4Ae x = 2e x
5Ae x = 2e x
5A = 2
2
A =
5
Por lo que
2 x
yp = e
5
y la solución general es
2 x
y = c1 cos 2x + c2 sin 2x + e
5
Ejemplo. Resolver
y 00 y = 8e x
47
Solución. Tenemos que
x
yh = c1 e + c2 ex
x
Como yp = Ae no es una buena propuesta, proponemos
x
yp = Axe
yp0 = Axe x
+ Ae x
x
= A (1 x) e
yp00 = A (1 x) e x
Ae x
x
= A (x 2) e
sustituyendo en la ED
x
A (x 2) e Axe x = 8e x
2Ae x = 8e x
2A = 8
A = 4
de aquí que
x
yp = 4xe
Por lo tanto, la solución general es
x
y = c1 e + c2 ex 4xe x
y 000 2y 00 y 0 + 2y = 8ex + 6e x
yp = Aex + Be x
48
yp0 = x Aex Be x
+ Aex + Be x
= (x + 1) Aex + (1 x) Be x
yp00 = (x + 2) Aex + (x 2) Be x
yp000 = (x + 3) Aex (x 3) Be x
Sustituyendo en la ecuación
y 000 2y 00 y 0 + 2y =
(x + 3) Aex (x 3) Be x
2 (x + 2) Aex 2 (x 2) Be x
(x + 1) Aex + (x 1) Be x
+2x Aex + Be x
= 2Aex + 6Be x
= 8ex + 6e x
2A = 8
6B = 6
A = 4
B = 1
Por lo que
yp = x 4ex + e x
x
= (x + c1 ) e + (4x + c2 ) ex + c3 e2x
y 00 + p (x) y 0 + q (x) y = 0
Supongamos que y = y1 (x) es una solución ya conocida. Ahora supongamos
que y2 (x) = u (x) y1 (x) y derivando
y20 = uy10 + u0 y1
49
y volviendo a derivar, tenemos que
y 00 + p (x) y 0 + q (x) y = 0
uy100 + 2u0 y10 + u00 y1 + p (uy10 + u0 y1 ) + quy1 = 0
uy100 + 2u0 y10 + u00 y1 + puy10 + pu0 y1 + quy1 = 0
(uy100 + puy10 + quy1 ) + 2u0 y10 + u00 y1 + pu0 y1 = 0
u (y100 + py10 + qy1 ) + 2u0 y10 + u00 y1 + pu0 y1 = 0
v 0 y1 + (2y10 + py1 ) v = 0
v 0 y1 = (2y10 + py1 ) v
v0 2y10 + py1
=
v y1
Z
2y10 + py1
ln v = dx
y1
R 0 +py
2y1 1
y1 dx
v = e
R y0
2 y1 +p dx
1
v = e
R
2 ln y1
v = e e pdx
2 R
v = eln yR1 e pdx
e pdx
v =
y12
R
pdx
0 e
u =
y12
entonces Z R !
pdx
e
u= dx
y12
50
como y2 = y1 u, entonces
Z R
pdx
e
y2 (x) = y1 (x) dx
y12
x2 y 00 2xy 0 4y =
0, x > 0
1
y1 (x) =
x
= x 1
2 1 1
x2 2x 4 = 0
x3 x2 x
2 2 4
+ = 0
x x x
4 4
= 0
x x
0 = 0
x2 y 00 2xy 0 4y = 0
00 2 0 4
y y y = 0
x x2
51
Z R
1 e ( 2
x )dx
y2 (x) = dx
x 1 2
x
Z R dx
1 e2 x
= 1 dx
x x2
Z 2 ln x
1 e
= 1 dx
x x2
Z ln x2
1 e
= 1 dx
x x2
Z !
x2
1 1
= 1 dx
x x2
Z
1
= x4 dx
x
1 x5
=
x 5
x4
=
5
c1
y= + c2 x4
x
Tarea: Comprobar que es solución y que son linealmente independientes y1
y y2 .
y y2
W (y1 ; y2 ) = 10
y1 y20
1 x4 1 x4
W ; = x 5
x 5 1 4x3
x2 5
3
1 4x 1 x4
=
x 5 x2 5
4x2 x2
= +
5 5
= x2 6= 0, porque x > 0
52
Supongamos que la forma de la solución es
yp = uy1 + vy2
yp0 = uy10 + u0 y1 + vy20 + v 0 y2
yp00 = u y100 + u0 y10 + u0 y10 + u00 y1
+v y200 + v 0 y20 + v 0 y20 + v 00 y2
= u y100 + 2u0 y10 + u00 y1
+v y200 + 2v 0 y20 + v 00 y2
d
p (u0 y1 + v 0 y2 ) + (u0 y10 + v 0 y20 ) + (u0 y1 + v 0 y2 ) = f (x)
dx
Sea
u 0 y1 + v 0 y2 = 0
u0 y10 + v 0 y20 = f (x)
y1 y2
W (y1 ; y2 ) =
y10 y20
53
0 y2
f (x) y20
u0 =
W (y1 ; y2 )
y2 f (x)
=
W (y1 ; y2 )
Z
y2 f (x)
u = dx
W (y1 ; y2 )
y1 0
y10 f (x)
v0 =
W (y1 ; y2 )
y1 f (x)
v0 =
W (y1 ; y2 )
Z
y1 f (x)
v = dx
W (y1 ; y2 )
entonces
x x
y1 = e y y2 = xe
y
x
f (x) = 4e
usando la fórmula de variación de parámetros, tenemos que
Z Z
y2 f (x) y1 f (x)
yp = y 1 dx + y2 dx
W (y1 ; y2 ) W (y1 ; y2 )
y
x x e x xe x
W e ; xe =
e x (x 1) e x
2x
= (x 1) e + xe 2x
2x 2x
= xe +e + xe 2x
= e 2x
54
Z x x Z x
x xe (4e ) x e (4e x )
yp = e dx + xe dx
e 2x e 2x
Z Z
= 4e x xdx + 4xe x dx
x x2 x
= 4e + 4xe (x)
2
= 2x2 e x
+ 4x2 e x
= 2x2 e x
Por lo tanto, la solución general es
x x
y = c1 e + c2 xe + 2x2 e x
y 00 + 4y = cos2 2x
Solución. La ecuación característica es
2
+4 = 0
= 2i
yh = c1 cos 2x + c2 sin 2x
de aquí que y1 = cos 2x y y2 = sin 2x, f (x) = cos2 2x, encontremos W (y1 ; y2 )
cos 2x sin 2x
W (y1 ; y2 ) =
2 sin 2x 2 cos 2x
= 2 cos2 2x + 2 sin2 2x
= 2 cos2 2x + sin2 2x
= 2
Z Z
y2 f (x) y1 f (x)
yp = y1 dx + y2 dx
W (y1 ; y2 ) W (y1 ; y2 )
Z Z
sin 2x cos2 2x cos 2x cos2 2x
= cos 2x dx + sin 2x dx
2 2
Z
cos 2x 1 1 sin 2x
= cos 2x cos 6x + cos3 2xdx
4 16 48 2
cos2 2x cos 2x cos 6x 3 sin2 2x 1
= + + sin 2x sin 6x
64 192 16 48
11 5 3
= cos 8x cos 4x
128 384 32
La solución general es
11 5 3
y = c1 cos 2x + c2 sin 2x + cos 8x cos 4x
128 384 32
55
Generalizando para una edo de tercer orden
0 y 2 y3
0 y20 y30
f (x) y200 y300
u0 =
W (y1 ; y2 ; y3 )
Z
(y2 y30 y3 y20 ) f (x)
u = dx
W (y1 ; y2 ; y3 )
y1 0 y3
y10 0 y30
y1 f (x) y300
00
v0 =
W (y1 ; y2 ; y3 )
Z
(y3 y10 y1 y30 ) f (x)
v = dx
W (y1 ; y2 ; y3 )
y 1 y2 0
y10 y20 0
y100 y200 f (x)
w0 =
W (y1 ; y2 ; y3 )
Z
(y1 y20 y2 y10 ) f (x)
w = dx
W (y1 ; y2 ; y3 )
y1 = e x , y2 = xe x , y3 = x2 e x
y10 = e x , y20 = (x 1) e x , y30 = x2 2x e x
y100 = e x
, y200 = (x 2) e x
, y300 = x2 4x + 2 e x
x
f (x) = e
e x xe x x2 e x
x x 2 x
W e ; xe ;x e = e x (x 1) e x 2
x 2x e x
e x (x 2) e x x2
4x + 2 e x
3x
= 2e
56
Z x
xe x2 2x e x
x2 e x
[ (x 1) e x
] e x
u = dx
2e 3x
Z
x x2 2x e 2x 2
x [ (x 1)] e 3x
= dx
2e 3x
Z 3 2 3
x + 2x + x x2
= dx
2
Z
x2
= dx
2
Z
1
= x2 dx
2
1 3
= x
6
Z
(y3 y10 y1 y30 ) f (x)
v = dx
W (y1 ; y2 ; y3 )
Z
x2 e x ( e x ) (e x ) x2 2x e x
e x
= dx
2e 3x
Z
x2 e 2x x2 2x e 2x e x
= dx
2e 3x
Z
2x2 + 2x e 3x
= dx
2e 3x
Z
= x2 + 2 dx
1 3
= x + 2x
3
Z
(y1 y20 y2 y10 ) f (x)
w = dx
W (y1 ; y2 ; y3 )
Z
f(e x ) [ (x 1) e x ] (xe x
)( e x
)g e x
= 3x
dx
2e
Z 3x
( x + 1 + x) e
= dx
2e 3x
Z
dx
=
2
1
= x
2
x
yp = e u + xe x v + x2 e x w
1 3 1 3 1
yp = e x x + xe x x + 2x + x2 e x
x
6 3 2
1 4 x 2 3 x
= x e + x e + 2x2 e x
3 3
57
Por lo tanto, la solución general es
y = yh + yp
x x 1 4 2
= c1 e + c2 xe + c3 x2 e x
x e x
+ x3 e x
+ 2x2 e x
3 3
x x 1 4 2
= c1 e + c2 xe + c3 x2 e x
x e x
+ x3 e x
3 3
3. Transformadas de Laplace
3.1. De…nición, propiedades de linealidad y formulas bási-
cas de transformada de Laplace.
De…nición. Sea f (t) una función de orden exponencial y seccionalmente
continua
Z 1
L ff (t)g : = f (t) e st dt
0
= F (s)
1. f (t) = 1
2. f (t) = t
3. f (t) = t2
4. f (t) = t3
5. f (t) = eat
Solución.
58
1.
Z 1
st
L f1g = 1e dt
0
Z 1
st
= e dt
0
1
1 st
= e
s 0
1
= , Re (s) > 0
s
2. Z 1
st
L ftg = te dt
0
usando integración por partes
st
t e
1 st
1 se
1 st
0 s2 e
Z 1 1
st t st 1 st
te dt = e e
0 s s2 0
1
=
s2
3. Z 1
2
L t = t2 e st
dt
0
Integrando por partes, tenemos que
t2 e st
1 st
2t se
1 st
2 s2 e
1 st
0 s3 e
Z 1
L t2 = t2 e st
dt
0
1
t2 2t 2 st
= e
s s2 s3 0
2
=
s3
59
4. Z 1
L t3 = t3 e st
dt
0
Integrando por partes, tenemos que
t3 e st
1
3t2 se
st
1 st
6t s e
2
1 st
6 s3 e
1 st
0 s4 e
Z 1
3
L t = t3 e st
dt
0
1
t3 3t2 6t 6 st
= e
s s2 s3 s4 0
6
=
s4
Generalizando
Z 1
n
L ft g = tn e st
dt
0
n!
= , n 2 N0
sn+1
5.
Z 1
at
L e = eat e st
dt
Z0 1
(s a)t
= e dt
0
1
= ,s>a
s a
Teorema de linealidad
Sean f (t) y g (t) dos funciones, , 2 R
60
Ejemplo: Calcular la transformada de Laplace de la función f (t) = 3e 2t +
2t
4e .
Solución: Usando el teorema de linealidad y L feat g = s 1 a , tenemos que
2t
L ff (t)g = L 3e + 4e2t
= 3L e 2t + 4L e2t
3 4
= +
s+2 s 2
7s + 2
=
s2 4
Ejemplo: Calcular la transformada de Laplace para f (t) = 4t5 2e3t + 5t.
Solución.
1
L ei!t =
s i!
1 s + i!
=
s i! s + i!
s + i!
=
s2 + ! 2
s !
= 2 2
+i 2
s +! s + !2
61
Por lo tanto
s !
L fcos !tg + iL fsin !tg = +i 2
s2 + ! 2 s + !2
!
L fsin !tg =
s2 + ! 2
s
L fcos !tg =
s2 + ! 2
3t
Ejemplo: Calcular L 3e + 2 cos 3t .
Solución.
3t
L 3e + 2 cos 3t = 3L e 3t + 2L fcos 3tg
3 2s
= +
s + 3 s2 + 9
5s2 + 6s + 27
=
s3 + 3s2 + 9s + 27
5s2 + 6s + 27
F (s) =
s3 + 3s2 + 9s + 27
1 5s2 + 6s + 27
L fF (s)g = L 1
s + 3s2 + 9s + 27
3
5s2 + 6s + 27
= L 1
(s + 3) (s2 + 9)
2s 3
= L 1 +
s2 + 9 s + 3
2s 3
= L 1 +L 1
s2 + 9 s+3
3t
= 2 cos 3t + 3e
62
Z 1
st
L ff (t)g = f (t) e dt
0
Z 5 Z 10 Z 1
st st st
= 2e dt + (2 t) e dt + (cos t) e dt
0 5 10
10 10 1
2 st 5 2 st t 1 st e st ( s cos t + sin t)
= e + e + e +
s 0 s 5 s s 5 s2 + 1 10
2 5s 2 2 10s 2 5s 10 10s 1 10s 5 5s
= e + e + e e e + e
s s s s s s s
10s
1 e ( s cos 10 + sin 10)
+ e 5s
s s2 + 1
2 6 5s 13 10s e 10s ( s cos 10 + sin 10)
= + e e
s s s s2 + 1
Solución.
1 5 6s 3 1 5 1 6s
L + 2 = L +L
s 6 s2 + 9 2s + 8s + 10 s 6 s2 +9
1 3
+L
2s2 + 8s + 10
1 1 1 s
= 5L + 6L
s 6 s2 + 32
3 1 1
+ L
2 s2
+ 4s + 4 + 1
( )
3 1
= 5e6t + 2 cos 3t + L 1 2
2 (s + 2) + 1
3
= 5e6t + 2 cos 3t + e 2t
sin t
2
n o
1 5
Ejemplo: Determinar L (s+2)4
.
63
Solución.
( ) ( )
1 5 1 1
L 4 = 5L 4
(s + 2) (s + 2)
5t3 e 2t
=
6
7s 1
Ejemplo: Determine la transformada inversa de F (s) = s3 +6s2 +11s+6 :
Solución.
1 7s 1 1 7s 1
L = L
s3 + 6s2 + 11s + 6 (s + 3) (s + 2) (s + 1)
11 15 4
= L 1 +
s+3 s+2 s+1
11 15
= L 1 +L 1
s+3 s+2
4
+L 1
s+1
1 1
= 11L 1 + 15L 1
s+3 s+2
1
4L 1
s+1
3t 2t t
= 11e + 15e 4e
2
Ejemplo: Determine la transformada inversa de F (s) = s3s+s+9s+2
2 5s+3 .
Solución.
( )
1 s2 + 9s + 2 1 s2 + 9s + 2
L = L 2
s3 + s2 5s + 3 (s + 3) (s 1)
( )
1 1 1 1
= L +2 +3 2
s+3 s 1 (s 1)
( )
1 1 1 1 1 1
= L + 2L + 3L 2
s+3 s 1 (s 1)
3t
= e + 2et + 3tet
s2 + 9s + 2 A B C
2 = + +
(s + 3) (s 1) s + 3 s 1 (s 1)2
2
A (s 1) + B (s + 3) (s 1) + C (s + 3) = s2 + 9s + 2
2
A s 2s + 1 + B s2 + 2s 3 + C (s + 3) = s2 + 9s + 2
(A + B) s2 + ( 2A + 2B + C) s + (A 3B + 3C) = s2 + 9s + 2
64
A+B = 1
2A + 2B + C = 9
A 3B + 3C = 2
n dn
L ftn f (t)g = ( 1) F (s)
dsn
Ejemplo: Determine L ft sin 3tg.
Solución. Como L fsin 3tg = s23+9
1 d 3
L ft sin 3tg = ( 1)
ds s2 + 9
3 (2s)
= 1 2
(s2 + 9)
6s
= 2
(s2 + 9)
en general
n o
L f (n) (t) = sn F (s) sn 1
f (0) sn 2 0
f (0) sn 3 00
f (0) f (n 1)
(0)
65
Solución.
L ff (t)g = L e6t t2
= L t2 s !s 6
2
=
s3 s !s 6
2
= 3
(s 6)
2t
Ejemplo. Si f (t) = e sin 3t, calcular L ff (t)g.
Solución.
L ff (t)g = L e 2t sin 3t
= L fsin 3tgjs !s+2
3
=
s2 +9 s !s+2
3
= 2
(s + 2) + 9
3
= 2
s + 4s + 13
f (t) = U (t 3) + U (t 2)
L ff (t)g = L fU (t 3) + U (t 2)g
= L fU (t 3)g + L fU (t 2)g
1 3s 1 2s
= e + e
s s
1
= e 3s + e 2s
s
Observación. La función escalonada
66
1. Si a < b, se tiene que
8
< 0, t < a
U (t a) U (t b) = 1, a t < b
:
0, t b
2.
0; 0 t < a
f (t) U (t a) =
f (t) , t a
3. 8
< 0, t < a
f (t) [U (t a) U (t b)] = f (t) , a t < b
:
0, t b
Ejemplo. Sea 8
>
> 0, 0 < t < 1
<
1, 1 t < 2
f (t) =
>
> 1, 2 t < 3
:
0, t 3
Solución.
a)
L ff (t)g = L fU (t 1) 2U (t 2) + U (t 3)g
= L fU (t 1)g 2L fU (t 2)g + L fU (t 3)g
1
= e s 2e 2s
+e 3s
s
L fU (t a) g (t a)g = e as G (s)
1
L e as G (s) = U (t a) g (t a)
0, 0 t < 1
f (t) = 2
(t 1) , t 1
67
2
Solución. Como a = 1 y g (t 1) = (t 1) , entonces g (t) = t2 , expresando
la función f (t) con el escalón unitario
2
f (t) = U (t 1) (t 1)
y 00 2y 0 + 5y = 8e t
y (0) = 2
y 0 (0) = 12
L fy 00 2y 0 + 5yg = L 8e t
L fy 00 g 2L fy 0 g + 5L fyg = 8L e t
8
s2 Y sy (0) y 0 (0) 2 [sY y (0)] + 5Y =
s+1
Aplicando las condiciones iniciales
8
s2 Y 2s 12 2 (sY 2) + 5Y =
s+1
8
s2 Y 2s 12 2sY + 4 + 5Y =
s+1
8 2s + 8
s2 2s + 5 Y = +
s+1 1
2s (s + 5)
s2 2s + 5 Y =
s+1
2s2 + 10s
Y =
(s + 1) (s2 2s + 5)
68
1
y (t) = L fY g
1 2s2 + 10s
= L
(s + 1) (s2 2s + 5)
1 3s + 5 1
= L
s2 2s + 5 s + 1
( )
1 3s 5 1
= L 2 + 2
(s 1) + 22 (s 1) + 22 s 1
( ) ( )
1 s 1 1 2
= 3L 2 +L 2
(s 1) + 22 (s 1) + 22
1 1
L
s 1
( ) ( )
1 s 1 1 2
= 3L 2 +L 2
(s 1) + 22 (s 1) + 22
1 1
L
s 1
= 3et cos 2t + et sin 2t et
s4 4 Y (s) = s3 2s
3
s 2s
Y (s) =
s4 4
s3 2s
Y (s) = 2
(s 2) (s2 + 2)
s s2 2
Y (s) = 2
(s 2) (s2 + 2)
s
Y (s) =
s2 + 2
1 s
y (t) = L
s2 +2
p
= cos 2t
69
Ejercicio 7.8 Tarea 23.
< 0, 0 t <
y 00 + y = 1, t < 2 ; y (0) = 0, y 0 (0) = 0.
:
0; 2 t
Solución. Expresando la función f (t) en términos del escalón unitario, ten-
emos que 8
>
> 0, t < 0
<
0, 0 t <
f (t) =
>
> 1, t<2
:
0; 2 t
8
< 0, t < a
[U (t a) U (t b)] = 1, a t < b
:
0, t b
f (t) = U (t ) U (t 2 )
y 00 + y = U (t ) U (t 2 )
Aplicando la transformada de Laplace
L fy 00 g + L fyg = L fU (t )g L fU (t 2 )g
s 2 s
e e
s2 Y (s) sy (0) y 0 (0) + Y (s) =
s s
2 e s e 2 s
s Y (s) + Y (s) =
s
s 2 s
e e
s2 + 1 Y (s) =
s
s
e e 2 s
Y (s) =
s (s2 + 1)
s
1 e e 2 s
y (t) = L
s (s2 + 1)
( )
s 2 s s 2 s
1 e e s e e
= L
s s2 + 1
s 2 s
1 e s
1 e s
= L L 1 L e s +L 1 e 2 s
s s s2 + 1 s2 + 1
= U (t ) U (t 2 ) U (t ) cos (t ) + U (t 2 ) cos (t 2 )
= [cos (t 2 ) 1] U (t 2 ) + [1 cos (t )] U (t )
= [cos (t) 1] U (t 2 ) + [1 cos (t )] U (t )
1 as
L fe G (s)g = U (t a) g (t a)
70
4. Sistemas de ecuaciones diferenciales lineales
4.1. Teoría de sistemas de ecuaciones diferenciales.
De…nición. Un sistema de n ecuaciones diferenciales de la forma
dy1
= f1 (x; y1 ; y2 ; : : : ; yn )
dx
dy2
= f2 (x; y1 ; y2 ; : : : ; yn )
dx
..
.
dyn
= fn (x; y1 ; y2 ; : : : ; yn )
dx
es llamado sistema de ecuaciones diferenciales ordinarias de primer orden, o
simplemente sistema de primer orden.
De…nición. Un sistema de primer orden de la forma
dy1
= f1 (x; y1 ; y2 ; : : : ; yn )
dx
dy2
= f2 (x; y1 ; y2 ; : : : ; yn )
dx
..
.
dyn
= fn (x; y1 ; y2 ; : : : ; yn )
dx
junto con n condiciones iniciales
y1 (x0 ) = y10
y2 (x0 ) = y20
..
.
yn (x0 ) = yn0
f x; y; y 0 ; : : : ; y (n) = 0
y (n) = g x; y; y 0 ; : : : ; y (n 1)
71
introducimos n nuevas variables y1 ; y2 ; : : : ; yn de…nidas como sigue
y1 = y
y2 = y0
..
.
yn = y (n 1)
72
4.2. Método con transformada de Laplace para un sis-
tema.
73