Universidad Técnica Federico Santa Marı́a

Departamento de Matemática Versión preliminar

Matemática III (MAT023)

Ejercicios certamen 3
Temas: Ecuaciones de orden superior, transformada de Laplace

Soluciones
1. [Ecuación no homogénea con coeficientes variables 1] Muestre que ϕ (x) = eαx
d2 y dy
es una solución de la ecuación x dx 2 − (x + 2) dx + 2y = 0 para un valor de α adecuado y

determinar la solución general de

d2 y dy
x 2
− (x + 2) + 2y = x3
dx dx
para x > 0.

Solución. Si ϕ (x) = eαx es solución de la EDO, se debe cumplir

d2 (eαx ) d (eαx )
x − (x + 2) + 2 (eαx ) = 0 para toda x
dx2 dx
luego
α2 xeαx − (x + 2) αeαx + 2eαx = 0 para todo x
como eαx > 0 se debe cumplir

α2 x − (x + 2) α + 2 = 0

factorizando y agrupando términos

α (α − 1) x + 2 (1 − α) = 0

esto se cumple para α = 1. Tenemos que ϕ (x) = ex es solución de la homogénea.

Usando la fórmula de Abel podemos obtener una segunda solución
R −(x+2)
e− dx
Z
x
x
ϕ2 (x) = e dx
e2x
2
ex+ln x
Z
x
= e dx
e2x
Z
= e x
x2 e−x dx

= − x2 + 2x + 2

la ecuación homogénea tiene solución

yh (x) = αex + β x2 + 2x + 2

1
MAT023 (Apuntes de clases)
Universidad Técnica Federico Santa Marı́a
Departamento de Matemática Versión preliminar

Ahora buscamos una solución particular de la EDO no homogénea usando variación de

parámetros
yp (x) = C1 (x) ex + C2 (x) x2 + 2x + 2

Por el método tenemos

0 (x2 + 2x + 2)
x2 2x + 2
C10 (x) = = e−x x2 + 2x + 2


x 2
e x + 2x + 2
ex 2x + 2

e−x (x2 + 2x + 2) dx = −e−x (x2 + 4x + 6). De manera similar

R
integrando C1 (x) =
ex 0
e x x2
C20 (x) = = −1
ex x2 + 2x + 2
ex 2x + 2
integrando C2 (x) = −x. la solución particular propuesta por el método es
yp (x) = −e−x x2 + 4x + 6 ex + (−x) x2 + 2x + 2

= −x3 − 3x2 − 6x − 6
(notar que −3x2 − 6x − 6 esta en el kernel del operador diferencial), la solución general de la
EDO
d2 y dy
x 2 − (x + 2) + 2y = x3 para x > 0
dx dx
es yG (x) = αex +β (x2 + 2x + 2)−x3 −3x2 −6x−6 ( o bien yG (x) = αex +µ (x2 + 2x + 2)−x3 ).

2. [Ecuación no homogénea con coeficientes variables 2] Busque una solución polino-

mial de la ecuación
1 − x2 y 00 − 2xy 0 + 2y = 0

y resolver el P.V.I.
1 − x2 y 00 − 2xy 0 + 2y = x

y (2) = 0
y 0 (2) = 1

Solución. Buscamos una solución de la forma y (x) = xn , reemplazando en la EDO tenemos

1 − x2 n (n − 1) xn−2 − 2xnxn−1 + 2xn = 0

agrupoando términos
n (n − 1) xn−2 − n (n − 1) xn − 2nxn + 2xn = 0
⇐⇒
n−2 n
n (n − 1) x − x (n + 2) (n − 1) = 0

2
MAT023 (Apuntes de clases)
Universidad Técnica Federico Santa Marı́a
Departamento de Matemática Versión preliminar

tomando n = 1 se cumple la igualdad. La función y (x) = x es una solución.

Para resolver el PVI
1 − x2 y 00 − 2xy 0 + 2y = x

y (2) = 0
y 0 (2) = 1
primero determinamos la solución general de la EDO no homogénea. Usando la fórmula de
Abel podemos obtener una segunda solución de la homogénea de la forma
Z − R −2x2 dx
e 1−x
y2 (x) = x dx
x2
 
1 x−1
= x ln +1
2 x+1
se sigue que la solución de la homogenea en x > 1 es dada por
   
1 x−1
yh (x) = αx + β x ln +1
2 x+1
vamos a determinar una solución particular usando variación de parámetros
   
1 x−1
yp (x) = C1 (x) x + C2 (x) x ln +1
2 x+1
donde
1
x ln x−1


0 2 x+1

+1

x d 1 x−1
1−x2 dx 2
x ln x+1 + 1 1 x−1
C10 (x) = 1
= x + x2 ln
x ln x−1


x 2 x+1

+1
2 x+1
d 1 x−1
1 dx 2
x ln x+1 + 1
integrando
Z  
1 2 x−1
C1 (x) = x + x ln dx
2 x+1
1 3 1 1 1 1
= x ln (x − 1) − ln (x + 1) − ln (x − 1) − x3 ln (x + 1) + x2
6 6 6 6 3
mientras que
x 0
x
0
1 1−x2
C2 (x) = 1 x−1

= −x2
x 2
x ln x+1
+ 1

d
1 dx 12 x ln x−1
x+1
+1
3
integrando C2 (x) = − x3 , la solución particular que propone el método es
    
1 3 1 1 1 3 1 2 1 x−1
yp (x) = x ln (x − 1) − ln (x + 1) − ln (x − 1) − x ln (x + 1) + x x + x ln
6 6 6 6 3 2 x+1
1
= − x ln x2 − 1


6
3
MAT023 (Apuntes de clases)
Universidad Técnica Federico Santa Marı́a
Departamento de Matemática Versión preliminar

La solución general para x > 1 es

   
1 x−1 1
+ 1 − x ln x2 − 1


yG (x) = αx + β x ln
2 x+1 6
ahora buscamos las constantes para cumplir con las condiciones iniciales
   
1 2−1 1
0 = α2 + β (2) ln + 1 − (2) ln (4 − 1)
2 2+1 6
calculamos la derivada
     
0 d 1 x−1 1 2


yG (x) = αx + β x ln + 1 − x ln x − 1
dx 2 x+1 6
     
d 1 x−1 d 1 2


= α+β x ln +1 + − x ln x − 1
dx 2 x+1 dx 6
        
1 x−1 1 2 1 2

1 2x
= α+β ln + + − ln x − 1 + − x 2
2 x+1 2 x2 − 1 6 6 (x − 1)
        
1 2−1 1 2 1 1 8
1=α+β ln + + − ln (4 − 1) + −
2 2+1 2 4−1 6 6 (4 − 1)
resolvemos el sistema
   
1 2−1 1
0 = α2 + β (2) ln + 1 − (2) ln (4 − 1)
2 2+1 6
        
1 2−1 1 2 1 1 8
1 = α+β ln + + − ln (4 − 1) + −
2 2+1 2 4−1 6 6 (4 − 1)
la solución encontrada es α = 13
3
− 25
6
ln 3, β = − 26
3
y reemplazamos para obtener la solución
del PVI:
      
13 25 26 1 x−1 1
+ 1 − x ln x2 − 1


yG (x) = − ln 3 x + − x ln
3 6 3 2 x+1 6

3. [Ecuaciones de coeficientes constantes 1]

a) Hallar una ecuación diferencial lineal con solución general

c1 e3t cos (2t) + c2 e3t sin (2t) + t3
Solución. Las funciones e3t cos (2t) y e3t sin (2t) están en el núcleo del operador (D − 3)2 +
22 = D2 − 6D + 13 y
d2 3
d 3

2
t − 6 t + 13t3 = 13t3 − 18t2 + 6t
dt dt
se sigue que c1 e3t cos (2t) + c2 e3t sin (2t) + t3 es la solución general de
y 00 − 6y 0 + 13y = 13t3 − 18t2 + 6t

4
MAT023 (Apuntes de clases)
Universidad Técnica Federico Santa Marı́a
Departamento de Matemática Versión preliminar

b) Considere la ecuación y 00 + q (x) y = 0 donde q es continua en todo R y dos soluciones

y1 , y2 las cuales satisfacen y1 (0) = 1, y10 (0) = 0 y y2 (0) = 3, y20 (0) = 1. Si W (x) es el
Wronskiano de y1 e y2 en x demuestre que W (x) = 1 para todo x donde las soluciones
estén definidas.
Solución. El teorema de Abel para el Wronskiano nos dice que, el Wronskiano de la
ecuación y 00 + Rq (x) y = 0 es constante (pues la función que acompaña a y 0 es nula luego
W (x) = Ke− 0dx ). Luego
y1 (x) y2 (x)
K=
y10 (x) y20 (x)
si reemplazamos x = 0 podemos obtener el valor de la constante

1 3
K= = 1.
0 1

4. [Ecuaciones de coeficientes constantes 2] Considere la EDO homogénea de coeficientes

constantes

d2 y dy
α2 + 1 + β +y =0
dx2 dx
dividir el plano αβ en regiones según la naturaleza de las soluciones de la EDO.

Solución. Las soluciones de una EDO de coeficientes constantes dependen de las raices de
la ecuación caracterı́stica, en este caso

α2 + 1 m2 + βm + 1 = 0

las cuales dependen del discriminante de la ecuación que corresponde a

∆ = β 2 − 4 α2 + 1

se presentan tres casos:

β 2 − 4 (α2 + 1) > 0, las raı́ces son reales y distintas

β 2 − 4 (α2 + 1) = 0, raı́z real repetida
β 2 − 4 (α2 + 1) < 0 , raı́ces complejas conjugadas

Note que la curva que separa las regiones es

β 2 − 4 α2 + 1


= 0
⇐⇒
β 2 − 4α2 = 4
⇐⇒
β2
− α2 = 1
4
que corresponde a una hipérbola en el plano αβ, Identificar las regiones de los puntos
anteriores.

5
MAT023 (Apuntes de clases)
Universidad Técnica Federico Santa Marı́a
Departamento de Matemática Versión preliminar

5. [Ecuaciones de coeficientes constantes 3] Determine los valores de α ∈ R para los

cuales, las soluciones y = ϕ (x) de la ecuación

d2 y 4 + α2
 
dy
+α + y=0
dx2 dx 4

cumplen lı́mx→+∞ e−4x ϕ (x) = 0.

Solución. La ecuación caracterı́stica de la edo es

4 + α2
 
2
m + αm + =0
4

la cual tiene raı́ces

q
4+α2


−α ± α2 − 4 4
m =
√ 2
−α ± −4
=
2
α
= − +i
2
α α α
luego las soluciones tienen la forma ϕ (x) = c1 e− 2 x cos (x) + c2 e− 2 x sin (x) = e− 2 x (c1 cos (x) + c2 sin (x)),
note que
α
lı́m e−4x ϕ (x) = lı́m e−(4+ 2 )x (c1 cos (x) + c2 sin (x))
x→+∞ x→+∞
α
este lı́mite es cero si 4 + 2
> 0 esto es α > −8.

6. [Ecuaciones de coeficientes constantes 4] Resolver la ecuacion diferencial

xy 00 + 2xy 0 + xy = e−x

para x > 0.

Solución. La ecuación es
e−x
y 00 + 2y 0 + y =
x
cuya solución general tiene la forma yG (x) = yh (x) + yp (x).

Vamos a resolver la EDO homogenea asociada

y 00 + 2y 0 + y = 0

que se puede escribir en la forma (usando operadores)

(D + 1)2 y = 0

esta tiene solución yh (x) = c1 e−x + c2 xe−x

6
MAT023 (Apuntes de clases)
Universidad Técnica Federico Santa Marı́a
Departamento de Matemática Versión preliminar

Para la solución particular usamos variación de parámetros

yp (x) = C1 (x) e−x + C2 (x) xe−x

donde
0 xe−x
e−x d
(xe−x )
C10 (x) = x dx
= −1
e−x xe−x
−e−x d
dx
(xe−x )
integrando C1 (x) = −x. Por otro lado

e−x 0
e−x
−e−x 1
C20 (x) = x
=
e−x xe −x x
−e−x d
dx
(xe−x )

integrando C2 (x) = ln x, se sigue

yp (x) = −xe−x + x ln xe−x

como −xe−x esta en el kernel del operador, la solución general se puede escribir en la forma

yG (x) = αe−x + βxe−x + x ln xe−x con α, β ∈ R.

√
7. [Cambio de variable 1] Use el cambio de variable t = x para resolver la EDO
√

4xy 00 + 2 − 8 x y 0 − 5y = 1

para x > 0.
Solución. Segun la indicación, definamos

Q (t) = y t2

se sigue que

Q0 (t) = y 0 t2 2t

Q00 (t) = y 00 t2 4t2 + y 0 t2 2

√
reemplazando el cambio x = t2 en la ecuación 4xy 00 (x) + (2 − 8 x) y 0 (x) − 5y (x) = 1 se tiene

4t2 y 00 t2 + (2 − 8t) y 0 t2 − 5y t2 = 1

la que podemos agrupar en la forma

 2 00 2

4t y t + 2y 0 t2 − 4 y 0 t2 2t − 5 5y t2 = 1

 
 


que corresponde a
Q00 (t) − 4Q0 (t) − 5Q (t) = 1

7
MAT023 (Apuntes de clases)
Universidad Técnica Federico Santa Marı́a
Departamento de Matemática Versión preliminar

resolvemos esta EDO de coeficientes constantes no homogénea para obtener

1
Q (t) = αe5t + βe−t −
5
√
y finalmente, como y (t2 ) = Q (t) se sigue y (x) = Q ( x), esto es
√ √ 1
y (x) = αe5 x
+ βe− x
− con α, β ∈ R
5

8. [Cambio de variable 2]

a) Encuentre la solución general de

d2 y
dy
x 2
− 6x2 + 1 + 8x3 y = 0 para x > 0,
dx dx
usando el cambio t = x2 .
Solución. Definamos √ 
Q (t) = y t
entonces
√  1
Q0 (t) = y 0 t t−1/2
2
√  1 √  −1
Q00 (t) = y 00 t t−1 + y 0 t t−3/2
4 4
√
reemplazando el cambio x = t en la ecuación nos queda
√ 00 √  √  √ 
ty t − (6t + 1) y 0 t + 8t3/2 y t =0

multiplicando por t−3/2

√ 
√  √ 
t−1 y 00 t − 6t−1/2 + t−3/2 y 0 t + 8y t =0

agrupamos términos
√ t−3/2 0 √  1 −1/2 0 √  h √ i
 −1     
t 00
4 y t − y t − 12 t y t +8 y t =0
4 4 2

se sigue
4Q00 (t) − 12Q0 (t) + 8Q (t) = 0
esta ecuación tiene solución Q (t) = αe2t + βet . Como y (x) = Q (x2 ) se sigue que
2 2
y (x) = αe2x + βex
d y 2 dy
2 3
es solución de la EDO homogénea x dx 2 − (6x + 1) dx + 8x y = 0 para x > 0.

8
MAT023 (Apuntes de clases)
Universidad Técnica Federico Santa Marı́a
Departamento de Matemática Versión preliminar

b) Encuentre la solución general de la ecuación diferencial:

d2 y
dy
x 2 − 6x2 + 1 + 8x3 y = 8x5 − 12x3
dx dx
Solución. Usar la parte anterior con variación de parámetros para obtener la solución
2 2
particular yp (x) = x2 , con esto la solución general es dada por yG (x) = αe2x + βex + x2
con α, β ∈ R.

9. [Cambio de variable 3] Resolver la ecuación diferencial

d2 y dy
x2 2
+x + y = tan (ln x)
dx dx
para 1 < x < eπ/2 .
Solución. Es una EDO de Euler-Cauchy, usar el cambio x = et y Q (t) = y (et ) para
transformar la EDO en
Q00 (t) + Q (t) = tan t
resolver esta EDO de coeficientes constantes y encontrar la solución particular con variación
de parámetros, se obtiene
 
1 1 − sin t
Q (t) = α cos t + β sin t + ln cos t
2 1 + sin t
y volvemos a la variable original
s
1 − sin (ln x)
y (x) = α cos (ln x) + β sin (ln x) + cos (ln x) ln
1 + sin (ln x)

con α, β ∈ R.

10. [Laplace 1]

a) Determine la transformada de Laplace inversa de:

π
1) 2
− arctan s
Solución. Supongamos
π
L [f (t)] (s) = − arctan s
2
entonces
d 1
L [f (t)] (s) = − 2
ds s +1
d 1
pero ds L [f (t)] (s) = L [−tf (t)] (s) y − s2 +1 = L [− sin t] (s), se sigue
L [−tf (t)] (s) = L [− sin t] (s)
y finalmente
sin t
f (t) =
t
9
MAT023 (Apuntes de clases)
Universidad Técnica Federico Santa Marı́a
Departamento de Matemática Versión preliminar

 
s2 +1
2) ln s2 −4
Solución.
  2 
d s +1 d
ln s2 + 1 − ln s2 − 4




ln 2 =
ds s −4 ds
2s 2s
= 2 − 2
s +1 s −4
s 1 1
= 2 2 − −
s + 1 s − 2 s +  2
2t −2t
= L 2 cos t − e − e
 2 
Si L [f (t)] (s) = ln ss2 −4
+1 d
entonces ds L [f (t)] (s) = L [−tf (t)] (s), se sigue

−tf (t) = 2 cos t − e2t − e−2t

de donde
2 cos t − e2t − e−2t
f (t) =
−t
2t −2t
e + e − 2 cos t
=
t
1
3) s2 +αs+β
donde α, β ∈ R (varios casos distintos)
Supongamos que s2 + αs + β = 0 tiene raı́ces reales y distintas r1 y r2 entonces
1 1
=
s2 + αs + β (s − r1 ) (s − r2 )
  
1 1 1
= −
r1 − r2 s − r1 s − r2
en este caso
 
−1 1 1 1
L 2
(t) = er1 t − er2 t
s + αs + β r1 − r2 r1 − r2

Si s2 + αs + β = 0 tiene raı́ces reales iguales r entonces

1 1
=
s2 + αs + β (s − r)2

en tal caso  
−1 1
L (t) = tert
s2 + αs + β

10
MAT023 (Apuntes de clases)
Universidad Técnica Federico Santa Marı́a
Departamento de Matemática Versión preliminar

Si s2 + αs + β = 0 tiene raı́ces complejas, completamos cuadrados

1 1
= 2
s2 + αs + β
q
α α2 α2



s2 + 2s 2
+ 4
+ β− 4

1
= q 2
α 2 α2


s+ 2
+ β− 4
q
2
1 β − α4
= q q 2
α2
β− 4
α 2


s+ 2 + β− α2
4

entonces
  r !
1 1 −α α2
L−1 2
(t) = q e2
t
sin β− t
s + αs + β β− α2 4
4

b) Determine la transformada de Laplace de:

sin t
1) t
(ver anterior)
2) f (t) = |t−2|
t−2
Solución.

1 si t > 2
f (t) =
−1 si t < 2
= −1 + 2H (t − 2)
se sigue
e−2s
 
|t − 2| −1
L (s) = +2
t−2 s s
Rt
3) t 0 et−u sin udu
Solución. Note que la función es t (et ∗ sin t) luego
d 
L t et ∗ sin t (s) = − L et ∗ sin t (s)

 
ds  
d 1 1
= −
ds s − 1 s2 + 1
3s2 − 2s + 1
=
(s − 1)2 (s2 + 1)2

11. [Laplace 2] Use la transformada de Laplace para resolver la ecuación diferencial

d2 y dy
t 2 + (1 − 2t) − 2y = 0
dt dt
con las condiciones iniciales y (0) = 1, y 0 (0) = 2.

11
MAT023 (Apuntes de clases)
Universidad Técnica Federico Santa Marı́a
Departamento de Matemática Versión preliminar

Solución. Aplicamos directamente la transformada

 2     
d dy dy d dy
− L 2
(s) + L (s) + 2 L (s) − 2L [y] (s) = 0
ds dt dt ds dt
ahora usamos las propiedades de las derivadas
d 2 d
s F (s) − sy (0) − y 0 (0) + (sF (s) − y (0)) + 2 (sF (s) − y (0)) − 2F (s) = 0


−
ds ds
donde F (s) = L [y] (s)
d 2
d
− s F (s) − s − 2 + (sF (s) − 1) + 2 (sF (s) − 1) − 2F (s) = 0
ds ds
⇐⇒
2 0 0
−2sF (s) − s F (s) + 1 + sF (s) − 1 + 2 (F (s) + sF (s)) − 2F (s) = 0

agrupamos términos
F 0 (s) −s2 + 2s + F (s) (−s) = 0

esto es
F 0 (s) −1
=
F (s) s−2
integrando ln F (s) = ln (s − 2)−1 se sigue
K
F (s) =
s−2
de donde y (t) = Ke2t , como y (0) = 1, se sigue K = 1 y por tanto y (t) = e2t es la solución
buscada.

12. [Laplace 3] Resolver el problema de valor inicial usando transformada de Laplace

Z t
d2 y
= f (t) + (t − u) y (u) du
dt2 0
y (0) = 0
y 0 (0) = 0

donde 
0 si t ≤ π
f (t) =
sin t si t ≥ π
Solución. Podemos escribir la ecuación en la forma
d2 y
= H (t − π) sin t + t ∗ y (t)
dt2
aplicamos transformada
1
s2 L [y (t)] (s) = e−πs L [sin (t + π)] (s) + L [y (t)] (s)
s2
12
MAT023 (Apuntes de clases)
Universidad Técnica Federico Santa Marı́a
Departamento de Matemática Versión preliminar

como sin (t + π) = − sin t se sigue

 
1 1
s − 2 L [y (t)] (s) = −e−πs 2
2
s s +1

despejando
−πs −s2
L [y (t)] (s) = e
(s4 − 1) (s2 + 1)
descomponemos en fracciones parciales

−s2 1 1 1 1
= − + −
(s4 − 1) (s2 + 1) 8 (s + 1) 8 (s − 1) 4 (s2 + 1) 2 (s2 + 1)2

y tomamos Laplace inversa

−s2
 
1 −t 1 t 1
= L e − e + t cos t (s)
(s4 − 1) (s2 + 1) 8 8 4

se sigue  
−πs 1 −t 1 t 1
L [y (t)] (s) = e L e − e + t cos t (s)
8 8 4
se sigue  
1 −(t−π) 1 t−π 1
y (t) = H (t − π) e − e + (t − π) cos (t − π)
8 8 4

13
MAT023 (Apuntes de clases)